Deep Neural Networks (DNNs) are powerful algorithms that have been proven capable of extracting non-Gaussian information from weak lensing (WL) data sets. Understanding which features in the data determine the output of these nested, non-linear algorithms is an important but challenging task. We analyze a DNN that has been found in previous work to accurately recover cosmological parameters in simulated maps of the WL convergence ($\kappa$). We derive constraints on the cosmological parameter pair $(\Omega_m,\sigma_8)$ from a combination of three commonly used WL statistics (power spectrum, lensing peaks, and Minkowski functionals), using ray-traced simulated $\kappa$ maps. We show that the network can improve the inferred parameter constraints relative to this combination by $20\%$ even in the presence of realistic levels of shape noise. We apply a series of well established saliency methods to interpret the DNN and find that the most relevant pixels are those with extreme $\kappa$ values. For noiseless maps, regions with negative $\kappa$ account for $86-69\%$ of the attribution of the DNN output, defined as the square of the saliency in input space. In the presence of shape nose, the attribution concentrates in high convergence regions, with $36-68\%$ of the attribution in regions with $\kappa > 3 \sigma_{\kappa}$.

The heaviest neutron stars and lightest black holes expected to be produced by stellar evolution leave the mass-range $2.2$ M$_{\odot}\lesssim m \lesssim 5$ M$_\odot$ largely unpopulated. Objects found in this so-called lower mass gap likely originate from a distinct astrophysical process. Such an object, with mass $2.6$ M$_\odot$ was recently detected in the binary merger GW190814 through gravitational waves by LIGO/Virgo. Here we show that black holes in the mass gap are naturally assembled through mergers and accretion in AGN disks, and can subsequently participate in additional mergers. We compute the properties of AGN-assisted mergers involving neutron stars and black holes, accounting for accretion. We find that mergers in which one of the objects is in the lower mass gap represent up to $4$% of AGN-assisted mergers detectable by LIGO/Virgo. The lighter object of GW190814, with mass $2.6$ M$_\odot$, could have grown in an AGN disk through accretion. We find that the unexpectedly high total mass of 3.4 M$_\odot$ observed in the neutron star merger GW190425 may also be due to accretion in an AGN disk.

The astrophysical origin of gravitational wave (GW) events is one of the most timely problems in the wake of the LIGO/Virgo discoveries. In active galactic nuclei (AGN), binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Previous studies have suggested that binary black hole (BBH) mergers in AGN disks can contribute significantly to BBH mergers observed by GW interferometers. Here we examine the distribution of the effective spin parameter $\chi_\mathrm{eff}$ of this GW source population. We extend our semi-analytical model of binary formation and evolution in AGN disks by following the evolution of the binary orbital angular momenta and black hole (BH) spins. BH spins change due to gas accretion and BH mergers, while the binary orbital angular momenta evolve due to gas accretion and binary-single interactions. We find that the distribution of $\chi_\mathrm{eff}$ predicted by our AGN model is similar to the distribution observed during LIGO/Virgo O1 and O2. On the other hand, if radial migration of BHs is inefficient, $\chi_\mathrm{eff}$ is skewed toward higher values compared with the observed distribution, because of the paucity of scattering events that would randomize spin directions relative to the orbital plane. We suggest that high binary masses and the positive correlation between binary mass and the standard deviation of $\chi_\mathrm{eff}$ for chirp masses up to $\approx 20$ $\mathrm{M}_\odot$, can be possible signatures for mergers originating in AGN disks. Finally, hierarchical mergers in AGN disks naturally produce properties of the recent GW event GW190412, including a low mass ratio, a high primary BH spin, and a significant spin component in the orbital plane.

The formation of supermassive stars has generally been studied under the assumption of rapid accretion of pristine metal-free gas. Recently it was found, however, that gas enriched to metallicities up to $Z \sim 10^{-3}$ Z$_{\odot}$ can also facilitate supermassive star formation, as long as the total mass infall rate onto the protostar remains sufficiently high. We extend the analysis further by examining how the abundance of supermassive star candidate haloes would be affected if all haloes with super-critical infall rates, regardless of metallicity were included. We investigate this scenario by identifying all atomic cooling haloes in the Renaissance simulations with central mass infall rates exceeding a fixed threshold. We find that among these haloes with central mass infall rates above 0.1 M$_{\odot}$ yr$^{-1}$ approximately two-thirds of these haloes have metallicities of $Z > 10^{-3}$ Z$_{\odot}$. If metal mixing within these haloes is inefficient early in their assembly and pockets of metal-poor gas can remain then the number of haloes hosting supermassive stars can be increased by at least a factor of four. Additionally the centres of these high infall-rate haloes provide ideal environments in which to grow pre-existing black holes. Further research into the (supermassive) star formation dynamics of rapidly collapsing haloes, with inhomogeneous metal distributions, is required to gain more insight into both supermassive star formation in early galaxies as well as early black hole growth.

The bright quasar PG1302-102 has been identified as a candidate supermassive black hole binary from its near-sinusoidal optical variability. While the significance of its optical periodicity has been debated due to the stochastic variability of quasars, its multi-wavelength variability in the ultraviolet (UV) and optical bands is consistent with relativistic Doppler boost caused by the orbital motion in a binary. However, this conclusion was based previously on sparse UV data which were not taken simultaneously with the optical data. Here we report simultaneous follow-up observations of PG1302-102 with the Ultraviolet Optical Telescope on the Neil Gehrels Swift Observatory in six optical + UV bands. The additional nine Swift observations produce light curves roughly consistent with the trend under the Doppler boost hypothesis, which predicts that UV variability should track the optical, but with a ~2.2 times higher amplitude. We perform a statistical analysis to quantitatively test this hypothesis. We find that the data are consistent with the Doppler boost hypothesis when we compare the the amplitudes in optical B-band and UV light curves. However, the ratio of UV to V-band variability is larger than expected and is consistent with the Doppler model, only if either the UV/optical spectral slopes vary, the stochastic variability makes a large contribution in the UV, or the sparse new optical data underestimate the true optical variability. We have evidence for the latter from comparison with the optical light curve from ASAS-SN. Additionally, the simultaneous analysis of all four bands strongly disfavors the Doppler boost model whenever Swift V-band is involved. Additional, simultaneous optical + UV observations tracing out another cycle of the 5.2-year proposed periodicity should lead to a definitive conclusion.

We present a new semi-analytic model of the formation of the first stars. Our method takes dark matter halo merger trees (including 3-dimensional spatial information) from cosmological N-body simulations as input and applies analytic prescriptions to compute both the Population III (Pop III) and metal-enriched star formation histories. We have developed a novel method to accurately compute the major feedback processes affecting Pop III star formation: H$_2$ photodissociation from Lyman-Werner (LW) radiation, suppression of star formation due to inhomogeneous reionization, and metal enrichment via supernovae winds. Our method utilizes a grid-based approach relying on fast Fourier transforms (FFTs) to rapidly track the LW intensity, ionization fraction, and metallicity in 3-dimensions throughout the simulation box. We present simulations for a wide range of astrophysical model parameters from $z\approx 30-6$. Initially long-range LW feedback and local metal enrichment and reionization feedback dominate. However, for $z \lesssim 15$ we find that the star formation rate density (SFRD) of Pop III stars is impacted by the combination of external metal enrichment (metals from one halo polluting other pristine halos) and inhomogeneous reionization. We find that the interplay of these processes is particularly important for the Pop III SFRD at $z \lesssim 10$. Reionization feedback delays star formation long enough for metal bubbles to reach halos that would otherwise form Pop III stars. Including these effects can lead to more than an order of magnitude decrease in the Pop III SFRD at $z=6$ compared to LW feedback alone.

The astrophysical origin of gravitational wave (GW) events discovered by LIGO/VIRGO remains an outstanding puzzle. In active galactic nuclei (AGN), compact-object binaries form, evolve, and interact with a dense star cluster and a gas disk. An important question is whether and how binaries merge in these environments. To address this question, we have performed one-dimensional $N$-body simulations combined with a semi-analytical model which includes the formation, disruption, and evolution of binaries self-consistently. We point out that binaries can form in single-single interactions by the dissipation of kinetic energy in a gaseous medium. This ``gas capture'' binary formation channel contributes up to $97\,\%$ of gas-driven mergers and leads to a high merger rate in AGN disks even without pre-existing binaries. We find the merger rate to be in the range $\sim 0.02-60\,\mathrm{Gpc^{-3}yr^{-1}}$. The results are insensitive to the assumptions on gaseous hardening processes: we find that once they are formed, binaries merge efficiently via binary-single interactions even if these gaseous processes are neglected. We find that the average number of mergers per BH is $0.4$, and the probability for repeated mergers in 30 Myr is $\sim 0.21-0.45$. High BH masses due to repeated mergers, high eccentricities, and a significant Doppler drift of GWs are promising signatures which distinguish this merger channel from others. Furthermore, we find that gas-capture binaries reproduce the distribution of LMXBs in the Galactic center, including an outer cutoff at $\sim1$ pc due to the competition between migration and hardening by gas torques.

Among the potential gravitational wave (GW) sources for the upcoming space-based interferometer LISA are extreme- or intermediate-mass ratio inspirals (EMRI/IMRIs). These events involve the coalescence of supermassive black holes in the mass range $10^5M_{\odot}\lesssim M \lesssim 10^7M_{\odot}$ with companion BHs of much lower masses. A subset of E/IMRIs are expected to occur in the accretion discs of active galactic nuclei (AGN). Previous work has shown that torques exerted by the disc can interfere with the inspiral and cause a phase shift in the GW waveform detectable by LISA. Here we use a suite of two-dimensional hydrodynamical simulations with the moving-mesh code DISCO to present a systematic study of disc torques. We measure torques on an inspiraling BH and compute the corresponding waveform deviations as a function of the binary mass ratio $q\equiv M_2/M_1$, the disc viscosity ($\alpha$), and gas temperature (or equivalently Mach number; $\mathcal{M}$). We find that the absolute value of the gas torques is within an order of magnitude of previously determined planetary migration torques, but their precise value and sign depends non-trivially on the combination of these parameters, the inspiral rate, and the accretion rate onto the satellite BH. The gas imprint is generally detectable by LISA for binaries embedded in AGN discs with surface densities above $\Sigma_0\ge10^{4-6} \rm \, g \, cm^{-2}$, depending on $q$, $\alpha$ and $\mathcal{M}$. We find that the deviations are most pronounced in discs with higher viscosities, and for E/IMRIs detected at frequencies where LISA is most sensitive. Torques in colder discs exhibit a noticeable dependence on the GW-driven inspiral rate as well as strong fluctuations at late stages of the inspiral. Our results further suggest that LISA may be able to place constraints on AGN disc parameters and the physics of disc-satellite interaction.

Numerical studies of gas accretion onto supermassive black hole binaries (SMBHBs) have generally been limited to conditions where the circumbinary disk (CBD) is 10-100 times thicker than expected for disks in active galactic nuclei (AGN). This discrepancy arises from technical limitations, and also from publication bias toward replicating fiducial numerical models. Here we present the first systematic study of how the binary's orbital evolution varies with disk scale height. We report three key results: (1) Binary orbital evolution switches from outspiralling for warm disks (aspect ratio ~0.1), to inspiralling for more realistic cooler, thinner disks at a critical aspect ratio ~0.04, corresponding to orbital Mach number ~25. (2) The net torque on the binary arises from a competition between positive torque from gas orbiting close to the black holes, and negative torque from the inner edge of the CBD, which is denser for thinner disks. This leads to increasingly negative net torques on the binary for increasingly thin disks. (3) The accretion rate is modestly suppressed with increasing Mach number. We discuss how our results may influence modeling of the nano-Hz gravitational wave background, as well as estimates of the LISA merger event rate.

We examine the light curves of two quasars, motivated by recent suggestions that a supermassive black hole binary (SMBHB) can exhibit sharp lensing spikes. We model the variability of each light curve as due to a combination of two relativistic effects: the orbital relativistic Doppler boost and gravitational binary self-lensing. In order to model each system we extend previous Doppler plus self-lensing models to include eccentricity. The first quasar is identified in optical data as a binary candidate with a 20-yr period (Ark 120), and shows a prominent spike. For this source, we rule out the lensing hypothesis and disfavor the Doppler-boost hypothesis due to discrepancies in the measured vs. recovered values of the binary mass and optical spectral slope. The second source, which we nickname Spikey, is the rare case of an active galactic nucleus (AGN) identified in Kepler's high-quality, high-cadence photometric data. For this source, we find a model, consisting of a combination of a Doppler modulation and a narrow symmetric lensing spike, consistent with an eccentric SMBHB with a total mass of approximately 30 million solar masses, rest-frame orbital period T=418 days, eccentricity e=0.5, and seen at an inclination of 8 degrees from edge-on. This interpretation can be tested by monitoring Spikey for periodic behavior and recurring flares in the next few years. In preparation for such monitoring we present the first X-ray observations of this object taken by the Neil Gehrels Swift observatory.

We performed a series of numerical experiments to quantify the sensitivity of the predictions for weak lensing statistics obtained in raytracing DM-only simulations, to two hyper-parameters that influence the accuracy as well as the computational cost of the predictions: the thickness of the lens planes used to build past light-cones and the mass resolution of the underlying DM simulation. The statistics considered are the power spectrum and a series of non-Gaussian observables, including the one-point probability density function, lensing peaks, and Minkowski functionals. Counter-intuitively, we find that using thin lens planes ($< 60~h^{-1}$Mpc on a $240~h^{-1}$Mpc simulation box) suppresses the power spectrum over a broad range of scales beyond what would be acceptable for an LSST-type survey. A mass resolution of $7.2\times 10^{11}~h^{-1}\,M_{\odot}$ per DM particle (or 256$^3$ particles in a ($240~h^{-1}$Mpc)$^3$ box) is sufficient to extract information using the power spectrum and non-Gaussian statistics from weak lensing data at angular scales down to 1 arcmin with LSST-like levels of shape noise.

The rotational kinematic Sunyaev-Zeldovich (rkSZ) signal, imprinted on the cosmic microwave background (CMB) by the gaseous halos (spinning "atmospheres") of foreground galaxies, would be a novel probe of galaxy formation. Although the signal is too weak to detect in individual galaxies, we analyze the feasibility of its statistical detection via stacking CMB data on many galaxies for which the spin orientation can be estimated spectroscopically. We use an "optimistic" model, in which fully ionized atmospheres contain the cosmic baryon fraction and spin at the halo's circular velocity $v_{\rm circ}$, and a more realistic model, based on hydrodynamical simulations, with multi-phase atmospheres spinning at a fraction of $v_{\rm circ}$. We incorporate realistic noise estimates into our analysis. Using low-redshift galaxy properties from the MaNGA spectroscopic survey (with median halo mass of $6.6\times10^{11}\,M_{\odot}$), and CMB data quality from Planck, we find that a $3\sigma$ detection would require a few$\times 10^4$ galaxies, even in the optimistic model. This is too high for current surveys, but upcoming higher-angular resolution CMB experiments will significantly reduce the requirements: stacking CMB data on galaxy spins in a $\sim10$ deg$^2$ can rule out the optimistic models, and $\approx$350 deg$^2$ will suffice for a $3\sigma$ detection with ACT. As a proof-of-concept, we stacked Planck data on the position of $\approx2,000$ MaNGA galaxies, aligned with the galaxies' projected spin, and scaled to their halos' angular size. We rule out average temperature dipoles larger than $\approx1.9\,\mu$K around field spiral galaxies.

Approximately two hundred supermassive black holes (SMBHs) have been discovered within the first $\sim$Gyr after the Big Bang. One pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs). A possible obstacle to this scenario is that the collapsing gas fragments and forms a cluster of main-sequence stars. Here we raise the possibility that stellar collisions may be sufficiently frequent and energetic to inhibit the contraction of the massive protostar, avoiding strong UV radiation driven outflows, and allowing it to continue growing into an SMS. We investigate this scenario with semianalytic models incorporating star formation, gas accretion, dynamical friction from stars and gas, stellar collisions, and gas ejection. We find that when the collapsing gas fragments at a density of $\lesssim 3\times 10^{10}\,\mathrm{cm^{-3}}$, the central protostar contracts due to infrequent stellar mergers, and in turn photoevaporates the remaining collapsing gas, resulting in the formation of a $\lesssim 10^4~{\rm M_\odot}$ object. On the other hand, when the collapsing gas fragments at higher densities (expected for a metal-poor cloud with $Z\lesssim10^{-5}\,{\rm Z_\odot}$ with suppressed ${\rm H_2}$ abundance) the central protostar avoids contraction and keeps growing via frequent stellar mergers, reaching masses as high as $\sim 10^5-10^6\,{\rm M_\odot}$. We conclude that frequent stellar mergers represent a possible pathway to form massive BHs in the early universe.

Binary black hole mergers encode information about their environment and the astrophysical processes that led to their formation. Measuring the redshift dependence of their merger rate will help probe the formation and evolution of galaxies and the evolution of the star formation rate. Here we compute the cosmic evolution of the merger rate for stellar-mass binaries in the disks of Active Galactic Nuclei (AGNs). We focus on recent evolution out to redshift $z=2$, covering the accessible range of current Earth-based gravitational-wave observatories. On this scale, the AGN population density is the main contributor to redshift-dependence. We find that the AGN-assisted merger rate does not meaningfully evolve with redshift, differentiating this channel from field binaries and some other dynamical formation scenarios.

Despite the rapidly growing number of stellar-mass binary black hole mergers discovered through gravitational waves, the origin of these binaries is still not known. In galactic centers, black holes can be brought to each others' proximity by dynamical processes, resulting in mergers. It is also possible that black holes formed in previous mergers encounter new black holes, resulting in so-called hierarchical mergers. Hierarchical events carry signatures such as higher-than usual black hole mass and spin. Here we show that the recently reported gravitational-wave candidate, GW170817A, could be the result of such a hierarchical merger. In particular, its chirp mass $\sim40$ M$_\odot$ and effective spin of $\chi_{\rm eff}\sim0.5$ are the typically expected values from hierarchical mergers within the disks of active galactic nuclei. We find that the reconstructed parameters of GW170817A strongly favor a hierarchical merger origin over having been produced by an isolated binary origin (with an Odds ratio of $>10^3$, after accounting for differences between the expected rates of hierarchical versus isolated mergers)

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$\alpha$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $\Gamma$, and the X-ray-to-optical power slopes, $\alpha_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

We study the thermal evolution of UV-irradiated atomic cooling haloes using high-resolution three-dimensional hydrodynamic simulations. We consider the effect of H^- photodetachment by Ly{\alpha} cooling radiation trapped in the optically-thick cores of three such haloes, a process which has not been included in previous simulations. H^- is a precursor of molecular hydrogen, and therefore, its destruction can diminish the H2 abundance and cooling. Using a simple high-end estimate for the trapped Ly{\alpha} energy density, we find that H^- photodetachment by Ly{\alpha} decreases the critical UV flux for suppressing H2-cooling by up to a factor of \approx 5. With a more conservative estimate of the Ly{\alpha} energy density, we find the critical flux is decreased only by ~15-50 percent. Our results suggest that Ly{\alpha} radiation may have an important effect on the thermal evolution of UV-irradiated haloes, and therefore on the potential for massive black hole formation.

The variability of quasars across multiple wavelengths is a useful probe of physical conditions in active galactic nuclei. In particular, variable accretion rates, instabilities, and reverberation effects in the accretion disk of a supermassive black hole (SMBH) are expected to produce correlated flux variations in UV and optical bands. Recent work has further argued that binary quasars should exhibit strongly correlated UV and optical periodicities. Strong UV-optical correlations have indeed been established in small samples of up to approximately 30 quasars with well-sampled light curves, and have extended the "bluer-when-brighter" trend previously found within the optical bands. Here we further test the nature of quasar variability by examining the observed-frame UV-optical correlations in a large sample of 1,315 bright quasars with overlapping UV and optical light curves for the Galaxy Evolution Explorer (GALEX) and the Catalina Real-time Transient Survey (CRTS), respectively. We find that strong correlations exist in this much larger sample, but we rule out, at approximately 95% confidence, the simple hypothesis that the intrinsic UV and optical variations of all quasars are fully correlated. Our results therefore imply the existence of physical mechanism(s) that can generate uncorrelated optical and UV flux variations.

The existence of $\approx$10^9 Msun supermassive black holes (SMBHs) within the first billion year of the universe has stimulated numerous ideas for the prompt formation and rapid growth of BHs in the early universe. Here we review ways in which the seeds of massive BHs may have first assembled, how they may have subsequently grown as massive as $\approx$10^9 Msun, and how multi-messenger observations could distinguish between different SMBH assembly scenarios. We conclude the following: (1) The ultra-rare $\approx$10^9 Msun SMBHs represent only the tip of the iceberg. Early BHs likely fill a continuum from stellar-mass (approx. 10 Msun) to the super-massive ($\approx$10^9 Msun) regime, reflecting a range of initial masses and growth histories. (2) Stellar-mass BHs were likely left behind by the first generation of stars at redshifts as high as z=30, but their initial growth was typically stunted due to the shallow potential wells of their host galaxies. (3) Conditions in some larger, metal-poor galaxies soon became conducive to the rapid formation and growth of massive `seed' holes, via gas accretion and by mergers in dense stellar clusters. (4) BH masses depend on the environment (such as the number and properties of nearby radiation sources and the local baryonic streaming velocity), and on the metal enrichment and assembly history of the host galaxy. (5) Distinguishing between assembly mechanisms will be difficult, but a combination of observations by LISA (probing massive BH growth via mergers) and by deep multi-wavelength electromagnetic observations (probing growth via gas accretion) is particularly promising.

Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $\alpha = 0.03$ to $\alpha = 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.

The coalescence of a compact object with a $10^{4}-10^{7} {\rm M_\odot}$ supermassive black hole (SMBH) produces mHz gravitational waves (GWs) detectable by the future Laser Interferometer Space Antenna (LISA). If such an inspiral occurs in the accretion disc of an active galactic nucleus (AGN), the gas torques imprint a small deviation in the GW waveform. Here we present two-dimensional hydrodynamical simulations with the moving-mesh code DISCO of a BH inspiraling at the GW rate in a binary system with a mass ratio $q\!=\!M_2/M_1\!=\!10^{-3}$, embedded in an accretion disc. We assume a locally isothermal equation of state for the gas (with Mach number $\mathcal{M}=20$) and implement a standard $\alpha$-prescription for its viscosity (with $\alpha = 0.03$). We find disc torques on the binary that are weaker than in previous semi-analytic toy models, and are in the opposite direction: the gas disc slows down, rather than speeds up the inspiral. We compute the resulting deviations in the GW waveform, which scale linearly with the mass of the disc. The SNR of these deviations accumulates mostly at high frequencies, and becomes detectable in a 5-year LISA observation if the total phase shift exceeds a few radians. We find that this occurs if the disc surface density exceeds $\Sigma_0 \gtrsim 10^{2-3}\rm g\,cm^{-2}$, as may be the case in thin discs with near-Eddington accretion rates. Since the characteristic imprint on the GW signal is strongly dependent on disc parameters, a LISA detection of an intermediate mass ratio inspiral would probe the physics of AGN discs and migration.

Based on the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), Kilbinger et al. (2013) obtained a constraint on the amplitude of matter fluctuations of $\sigma_8({\Omega_\mathrm{m}}/0.27)^{0.6}=0.79\pm0.03$ from the two-point correlation function (2PCF). This is $\approx3\sigma$ lower than the value $0.89\pm0.01$ derived from Planck data on cosmic microwave background (CMB) anisotropies. On the other hand, based on the same CFHTLenS data, but using the power spectrum, and performing a different analysis, Liu et al. (2015) obtained the higher value of $\sigma_8({\Omega_\mathrm{m}}/0.27)^{0.64}=0.87^{+0.05}_{-0.06}$. We here investigate the origin of this difference, by performing a fair side-by-side comparison of the 2PCF and power spectrum analyses on CFHTLenS data. We find that these two statistics indeed deliver different results, even when applied to the same data in an otherwise identical procedure. We identify excess power in the data on small scales ($\ell>5,000$) driving the larger values inferred from the power spectrum. We speculate on the possible origin of this excess small-scale power. More generally, our results highlight the utility of analysing the 2PCF and the power spectrum in tandem, to discover (and to help control) systematic errors.

The origins of the stellar-mass black hole mergers discovered by LIGO/Virgo are still unknown. Here we show that, if migration traps develop in the \add{accretion} disks of Active Galactic Nuclei (AGNs) and promote the mergers of their captive black holes, the majority of black holes within disks will undergo hierarchical mergers---with one of the black holes being the remnant of a previous merger. 40% of AGN-assisted mergers detected by LIGO/Virgo will include a black hole with mass $\gtrsim 50$ M$_\odot$, the mass limit from stellar core collapse. Hierarchical mergers at traps in AGNs will exhibit black hole spins (anti-)aligned with the binary's orbital axis, a distinct property from other hierarchical channels. Our results are suggestive, although not definitive (with Odds ratio of $\sim 1$), that LIGO's heaviest merger so far, GW170729, could have originated from this channel.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

The localization of stellar-mass binary black hole mergers using gravitational waves is critical in understanding the properties of the binaries' host galaxies, observing possible electromagnetic emission from the mergers, or using them as a cosmological distance ladder. The precision of this localization can be substantially increased with prior astrophysical information about the binary system. In particular, constraining the inclination of the binary can reduce the distance uncertainty of the source. Here we present the first realistic set of localizations for binary black hole mergers, including different prior constraints on the binaries' inclinations. We find that prior information on the inclination can reduce the localization volume by a factor of 3. We discuss two astrophysical scenarios of interest: (i) follow-up searches for beamed electromagnetic/neutrino counterparts and (ii) mergers in the accretion disks of active galactic nuclei.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data.

The growing number of stellar-mass binary black hole mergers discovered by Advanced LIGO and Advanced Virgo are starting to constrain the binaries' origin and environment. However, we still lack sufficiently accurate modeling of binary formation channels to obtain strong constraints, or to identify sub-populations. One promising formation mechanism that could result in different black hole properties is binaries merging within the accretion disks of Active Galactic Nuclei (AGN). Here we show that the black holes' orbital alignment with the AGN disks preferentially selects heavier black holes. We carry out Monte Carlo simulations of orbital alignment with AGN disks, and find that AGNs harden the initial black hole mass function. Assuming an initial power law mass distribution $M_{\rm bh}^{-\beta}$, we find that the power law index changes by $\Delta \beta\sim1.3$, resulting in a more top-heavy population of merging black holes. This change is independent of the mass of, and accretion rate onto, the supermassive black hole in the center of the AGN. Our simulations predict an AGN-assisted merger rate of $\sim4$Gpc$^{-3}$yr$^{-1}$. With its hardened mass spectra, the AGN channel could be responsible for $10-50$% of gravitational-wave detections.

The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}

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.

We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist.

In this white paper we explore the capabilities required to identify and study supermassive black holes formed from heavy seeds ($\mathrm{M_{\bullet}} \sim 10^4 - 10^6 \, \mathrm{M_{\odot}}$) in the early Universe. To obtain an unequivocal detection of heavy seeds we need to probe mass scales of $\sim 10^{5-6} \, \mathrm{M_{\odot}}$ at redshift $z \gtrsim 10$. From this theoretical perspective, we review the observational requirements and how they compare with planned/proposed instruments, in the infrared, X-ray and gravitational waves realms. In conclusion, detecting heavy black hole seeds at $z \gtrsim 10$ in the next decade will be challenging but, according to current theoretical models, feasible with upcoming/proposed facilities. Their detection will be fundamental to understand the early history of the Universe, as well as its evolution until now. Shedding light on the dawn of black holes will certainly be one of the key tasks that the astronomical community will focus on in the next decade.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

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.

Intermediate-mass black holes (IMBHs) have masses of about 100 to 100,000 solar masses. They remain elusive. Observing IMBHs in present-day globular clusters (GCs) would validate a formation channel for seed black holes in the early universe and inform event predictions for gravitational wave facilities. Reaching a large number of GCs per galaxy is key, as models predict that only a few percent will have retained their gravitational-wave fostering IMBHs. Related, many galaxies will need to be examined to establish a robust sample of IMBHs in GCs. These needs can be meet by using a next-generation Very Large Array (ngVLA) to search for IMBHs in the GCs of hundreds of galaxies out to a distance of 25 Mpc. These galaxies hold tens of thousands of GCs in total. We describe how to convert an ngVLA signal from a GC to an IMBH mass according to a semi-empirical accretion model. Simulations of gas flows in GCs would help to improve the robustness of the conversion. Also, self-consistent dynamical models of GCs, with stellar and binary evolution in the presence of IMBHs, would help to improve IMBH retention predictions for present-day GCs.

Weak gravitational lensing is one of the most promising cosmological probes of the late universe. Several large ongoing (DES, KiDS, HSC) and planned (LSST, EUCLID, WFIRST) astronomical surveys attempt to collect even deeper and larger scale data on weak lensing. Due to gravitational collapse, the distribution of dark matter is non-Gaussian on small scales. However, observations are typically evaluated through the two-point correlation function of galaxy shear, which does not capture non-Gaussian features of the lensing maps. Previous studies attempted to extract non-Gaussian information from weak lensing observations through several higher-order statistics such as the three-point correlation function, peak counts or Minkowski-functionals. Deep convolutional neural networks (CNN) emerged in the field of computer vision with tremendous success, and they offer a new and very promising framework to extract information from 2 or 3-dimensional astronomical data sets, confirmed by recent studies on weak lensing. We show that a CNN is able to yield significantly stricter constraints of ($\sigma_8, \Omega_m$) cosmological parameters than the power spectrum using convergence maps generated by full N-body simulations and ray-tracing, at angular scales and shape noise levels relevant for future observations. In a scenario mimicking LSST or Euclid, the CNN yields 2.4-2.8 times smaller credible contours than the power spectrum, and 3.5-4.2 times smaller at noise levels corresponding to a deep space survey such as WFIRST. We also show that at shape noise levels achievable in future space surveys the CNN yields 1.4-2.1 times smaller contours than peak counts, a higher-order statistic capable of extracting non-Gaussian information from weak lensing maps.

The presence of massive neutrinos affects structure formation, leaving imprints on large-scale structure observables such as the weak lensing field. The common lensing analyses with two-point statistics are insensitive to the large amount of non-Gaussian information in the density field. We investigate non-Gaussian tools, in particular the Minkowski Functionals (MFs)---morphological descriptors including area, perimeter, and genus---in an attempt to recover the higher-order information. We use convergence maps from the Cosmological Massive Neutrino Simulations (\texttt{MassiveNus}) and assume galaxy noise, density, and redshift distribution for an LSST-like survey. We show that MFs are sensitive to the neutrino mass sum, and the sensitivity is redshift dependent and is non-Gaussian. We find that redshift tomography significantly improves the constraints on neutrino mass for MFs, compared to the improvements for the power spectrum. We attribute this to the stronger redshift dependence of neutrino effects on small scales. We then build an emulator to model the power spectrum and MFs, and study the constraints on $[M_{\nu}$, $\Omega_{m}$, $A_{s}]$ from the power spectrum, MFs, and their combination. We show that MFs significantly outperform the power spectrum in constraining neutrino mass, by more than a factor of four. However, a thorough study of the impact from systematics such as baryon physics and galaxy shape and redshift biases will be important to realize the full potential of MFs.

The discovery of the electromagnetic counterparts to the binary neutron star merger GW170817 has opened the era of GW+EM multi-messenger astronomy. Exploiting this breakthrough requires increasing samples to explore the diversity of kilonova behaviour and provide more stringent constraints on the Hubble constant, and tests of fundamental physics. LSST can play a key role in this field in the 2020s, when the gravitational wave detector network is expected to detect higher rates of merger events involving neutron stars ($\sim$10s per year) out to distances of several hundred Mpc. Here we propose comprehensive target-of-opportunity (ToOs) strategies for follow-up of gravitational-wave sources that will make LSST the premiere machine for discovery and early characterization for neutron star mergers and other gravitational-wave sources.

Massive black-hole (MBH) binaries, which are expected to form following the merger of their parent galaxies, produce gravitational waves which will be detectable by Pulsar Timing Arrays at nanohertz frequencies (year periods). While no confirmed, compact MBH binary systems have been seen in electromagnetic observations, a large number of candidates have recently been identified in optical surveys of AGN variability. Using a combination of cosmological, hydrodynamic simulations; comprehensive, semi-analytic binary merger models; and analytic AGN spectra and variability prescriptions; we calculate the expected electromagnetic detection rates of MBH binaries as periodically variable AGN. In particular, we consider two independent variability models: (i) Doppler boosting due to large orbital velocities, and (ii) hydrodynamic variability in which the fueling of MBH accretion disks is periodically modulated by the companion. Our models predict that numerous MBH binaries should be present and distinguishable in the existing data. In particular, our fiducial models produce an expectation value of $0.2$ (Doppler) and $5$ (hydrodynamic) binaries to be identifiable in CRTS, while $20$ and $100$ are expected after five years of LSST observations. The brightness variations in most systems are too small to be distinguishable, but almost $1\%$ of AGN at redshifts $z \lesssim 0.6$ could be in massive binaries. We analyze the predicted binary parameters of observable systems and their selection biases, and include an extensive discussion of our model parameters and uncertainties.

Large cosmological datasets have been probing the properties of our universe and constraining the parameters of dark matter and dark energy with increasing precision. Deep learning techniques have shown potential to be smarter, and to greatly outperform human-designed statistics.

Weak lensing surveys are reaching sensitivities at which uncertainties in the galaxy redshift distributions n(z) from photo-z errors degrade cosmological constraints. We use ray-tracing simulations and a simple treatment of photo-z errors to assess cosmological parameter biases from uncertainties in n(z) in an LSST-like survey. We use the power spectrum and the abundance of lensing peaks to infer cosmological parameters, and find that the former is somewhat more resilient to photo-z errors. We place conservative lower limits on the survey size at which different types of photo-z errors degrade LCDM (wCDM) parameter constraints by 50%. A residual constant photo-z bias of |dz| < 0.003(1+z), satisfying the current LSST requirement, does not significantly degrade constraints for surveys smaller than ~1300 (~490) square degrees using lensing peaks and ~6500 (~4900) square degrees using the power spectrum. Adopting a recent prediction for LSST's full photo-z probability distribution function (PDF), we find that simply approximating n(z) with the photo-z galaxy distribution directly computed from this PDF would degrade surveys as small as ~60 (~65) square degrees using lensing peaks or the power spectrum. Assuming that the centroid bias in each tomographic redshift bin can be removed from the photo-z galaxy distribution, using lensing peaks or the power spectrum still degrades surveys larger than ~200 (~255) or ~248 (~315) square degrees. These results imply that the expected broad photo-z PDF significantly biases parameters, which needs to be further mitigated using more sophisticated photo-z treatments.

The abundance of molecular hydrogen (H2), the primary coolant in primordial gas, is critical for the thermodynamic evolution and star-formation histories in early protogalaxies. Determining the photodissociation rate of H2 by an incident Lyman-Werner (LW) flux is thus crucial, but prohibitively expensive to calculate on the fly in simulations. The rate is sensitive to the H2 rovibrational distribution, which in turn depends on the gas density, temperature, and incident LW radiation field. We use the publicly available cloudy package to model primordial gas clouds and compare exact photodissociation rate calculations to commonly-used fitting formulae. We find the fit from Wolcott-Green et al. (2011) is most accurate for moderate densities n~10^3 cm^{-3} and temperatures, T~10^3K, and we provide a new fit, which captures the increase in the rate at higher densities and temperatures, owing to the increased excited rovibrational populations in this regime. Our new fit has typical errors of a few percent percent up to n =<10^7 cm^{-3}, T =< 8000K, and H2 column density NH2 =<10^{17} cm^{-2}, and can be easily utilized in simulations. We also show that pumping of the excited rovibrational states of H2 by a strong LW flux further modifies the level populations when the gas density is low, and noticeably decreases self-shielding for J_21 > 10^3 and n < 10^2 cm^{-3}. This may lower the "critical flux" at which primordial gas remains H2-poor in some protogalaxies, enabling massive black hole seed formation.

Supermassive black holes (SMBHs) pair and form bound binaries after their host galaxies merge. In a gas-rich merger, accretion discs are expected to form around the binary and its components. These discs control the binary orbital evolution until the system is compact enough for gravitational waves to drive the SMBHs to coalescence. In this work, we implemented a time-dependent 1D model to follow the long-term evolution of the coupled binary + discs system, from a separation of 10^5 down to 20 Schwarzschild radii. We run different models changing the system parameters, including the binary mass ratio q \leq 0.3 and a factor {\gamma} that controls the inflow across the gap created by the secondary. We find that our implementation yields higher residual masses and longer binary residence times than previous studies. Our main conclusion is the non-steady-state nature of the evolution of the system: the properties the disc had when the binary was still at large separations influence its whole evolution. To recover steady state, the binary residence time would have to be much longer than the inflow time-scale of the disc throughout their entire history, which in general is not satisfied.

We study low-density axisymmetric accretion flows onto black holes (BHs) with two-dimensional hydrodynamical simulations, adopting the $\alpha$-viscosity prescription. When the gas angular momentum is low enough to form a rotationally supported disk within the Bondi radius ($R_{\rm B}$), we find a global steady accretion solution. The solution consists of a rotational equilibrium distribution at $r\sim R_{\rm B}$, where the density follows $\rho \propto (1+R_{\rm B}/r)^{3/2}$, surrounding a geometrically thick and optically thin accretion disk at the centrifugal radius, where thermal energy generated by viscosity is transported via strong convection. Physical properties of the inner solution agree with those expected in convection-dominated accretion flows (CDAF; $\rho \propto r^{-1/2}$). In the inner CDAF solution, the gas inflow rate decreases towards the center due to convection ($\dot{M}\propto r$), and the net accretion rate (including both inflows and outflows) is strongly suppressed by several orders of magnitude from the Bondi accretion rate $\dot{M}_{\rm B}$ The net accretion rate depends on the viscous strength, following $\dot{M}/\dot{M}_{\rm B}\propto (\alpha/0.01)^{0.6}$. This solution holds for low accretion rates of $\dot{M}_{\rm B}/\dot{M}_{\rm Edd}< 10^{-3}$ having minimal radiation cooling, where $\dot{M}_{\rm Edd}$ is the Eddington rate. In a hot plasma at the bottom ($r<10^{-3}~R_{\rm B}$), thermal conduction would dominate the convective energy flux. Since suppression of the accretion by convection ceases, the final BH feeding rate is found to be $\dot{M}/\dot{M}_{\rm B} \sim 10^{-3}-10^{-2}$. This rate is as low as $\dot{M}/\dot{M}_{\rm Edd} \sim 10^{-7}-10^{-6}$ inferred for SgrA$^*$ and the nuclear BHs in M31 and M87, and can explain the low luminosities in these sources, without invoking any feedback mechanism.

Gravitational waves (GWs) in the nano-hertz band are great tools for understanding the cosmological evolution of supermassive black holes (SMBHs) in galactic nuclei. We consider SMBH binaries in high-$z$ ultra-luminous infrared galaxies (ULIRGs) as sources of a stochastic GW background (GWB). ULIRGs are likely associated with gas-rich galaxy mergers containing SMBHs that possibly occur at most once in the life of galaxies, unlike multiple dry mergers at low redshift. Adopting a well-established sample of ULIRGs, we study the properties of the GWB due to coalescing binary SMBHs in these galaxies. Since the ULIRG population peaks at $z>1.5$, the amplitude of the GWB is not affected even if BH mergers are delayed by as long as $\sim $ 10 Gyrs. Despite the rarity of the high-$z$ ULIRGs, we find a tension with the upper limits from Pulsar Timing Array (PTA) experiments. This result suggests that if a fraction $f_{\rm m,gal}$ of ULIRGs are associated with SMBH binaries, then no more than $20 f_{\rm m,gal}(\lambda_{\rm Edd}/0.3)^{5/3}(t_{\rm life}/30~{\rm Myr})~\%$ of the binary SMBHs in ULIRGs can merge within a Hubble time, for plausible values of the Eddington ratio of ULIRGs ($\lambda_{\rm Edd}$) and their lifetime ($t_{\rm life}$).

We propose the formation of massive pristine dark-matter (DM) halos with masses of $\sim 10^8~M_\odot$, due to the dynamical effects of frequent mergers in rare regions of the Universe with high baryonic streaming velocity relative to DM. Since the streaming motion prevents gas collapse into DM halos and delays prior star formation episodes, the gas remains metal-free until the halo virial temperatures $\gtrsim 2\times 10^4~{\rm K}$. The minimum cooling mass of DM halos is boosted by a factor of $\sim 10-30$ because frequent major mergers of halos further inhibit gas collapse. We use Monte Carlo merger trees to simulate the DM assembly history under a streaming velocity of twice the root-mean-square value, and estimate the number density of massive DM halos containing pristine gas as $\simeq 10^{-4}~{\rm cMpc}^{-3}$. When the gas infall begins, efficient Ly$\alpha$ cooling drives cold streams penetrating inside the halo and feeding a central galactic disk. When one stream collides with the disk, strong shock forms a dense and hot gas cloud, where the gas never forms H$_2$ molecules due to effective collisional dissociation. As a result, a massive gas cloud forms by gravitational instability and collapses directly into a massive black hole (BH) with $M_\bullet \sim 10^5~M_\odot$. Almost simultaneously, a galaxy with $M_{\star, \rm tot}\sim 10^6~M_\odot$ composed of Population III stars forms in the nuclear region. If the typical stellar mass is as high as $\sim 100~M_\odot$, the galaxy could be detected with the James Webb Space Telescope even at $z\gtrsim 15$. These massive seed BHs would be fed by continuous gas accretion from the host galaxy, and grow to be bright quasars observed at $z\gtrsim 6$.

Observations of high-redshift quasars indicate that super massive black holes (SMBHs) with masses greater than ${\sim}10^9~M_\odot$ were assembled within the first billion years after the Big Bang. It is unclear how such massive black holes formed so early. One possible explanation is that these SMBHs were seeded by `heavy' direct collapse black holes (DCBHs) with masses of $M_{\rm BH}\approx10^5~M_\odot$, but observations have not yet confirmed or refuted this scenario. In this Letter, we utilize a cosmological N-body simulation to demonstrate that before they grow roughly an order of magnitude in mass, DCBHs will have black hole mass to halo mass ratios much higher than expected for black hole remnants of Population III (Pop III) stars which have grown to the same mass (${\sim}10^6~M_\odot$). We also show that when $T_{\rm vir}\approx 10^4~{\rm K}$ halos (the potential sites of DCBH formation) merge with much larger nearby halos ($M_{\rm h} > 10^{10}~M_\odot$), they almost always orbit their larger host halos with a separation of a few kpc, which is sufficient to be spatially resolved with future X-ray and infrared telescopes. Thus, we propose that a future X-ray mission such as Lynx combined with infrared observations will be able to distinguish high-redshift DCBHs from smaller black hole seeds due to the unusually high black hole mass to stellar mass ratios of the faintest observed quasars, with inferred BH masses below ${\sim} 10^6 M_\odot$.

Weak lensing maps contain information beyond two-point statistics on small scales. Much recent work has tried to extract this information through a range of different observables or via nonlinear transformations of the lensing field. Here we train and apply a 2D convolutional neural network to simulated noiseless lensing maps covering 96 different cosmological models over a range of {$\Omega_m,\sigma_8$}. Using the area of the confidence contour in the {$\Omega_m,\sigma_8$} plane as a figure-of-merit, derived from simulated convergence maps smoothed on a scale of 1.0 arcmin, we show that the neural network yields $\approx 5 \times$ tighter constraints than the power spectrum, and $\approx 4 \times$ tighter than the lensing peaks. Such gains illustrate the extent to which weak lensing data encode cosmological information not accessible to the power spectrum or even other, non-Gaussian statistics such as lensing peaks.

We quantify the presence of Ly\alpha\ damping wing absorption from a partially-neutral intergalactic medium (IGM) in the spectrum of the $z=7.08$ QSO, ULASJ1120+0641. Using a Bayesian framework, we simultaneously account for uncertainties in: (i) the intrinsic QSO emission spectrum; and (ii) the distribution of cosmic HI patches during the epoch of reionisation (EoR). For (i) we use a new intrinsic Ly\alpha\ emission line reconstruction method (Greig et al.), sampling a covariance matrix of emission line properties built from a large database of moderate-$z$ QSOs. For (ii), we use the Evolution of 21-cm Structure (EOS; Mesinger et al.) simulations, which span a range of physically-motivated EoR models. We find strong evidence for the presence of damping wing absorption redward of Ly\alpha\ (where there is no contamination from the Ly\alpha\ forest). Our analysis implies that the EoR is not yet complete by $z=7.1$, with the volume-weighted IGM neutral fraction constrained to $\bar{x}_{\rm H\,{\scriptsize I}} = 0.40\substack{+0.21 -0.19}$ at $1\sigma$ ($\bar{x}_{\rm H\,{\scriptsize I}} = 0.40\substack{+0.41 -0.32}$ at $2\sigma$). This result is insensitive to the EoR morphology. Our detection of significant neutral HI in the IGM at $z=7.1$ is consistent with the latest Planck 2016 measurements of the CMB Thompson scattering optical depth (Planck Collaboration XLVII).

Black hole mergers detectable with LIGO can occur in active galactic nucleus (AGN) disks. Here we parameterize the merger rates, the mass spectrum and the spin spectrum of black holes (BH) in AGN disks. The predicted merger rate spans $\sim 10^{-4}-10^{4} \rm{Gpc}^{-1} \rm{yr}^{-1}$, so upper limits from LIGO ($<212\rm{Gpc}^{-1}\rm{yr}^{-1}$) already constrain it. The predicted mass spectrum has the form of a broken power-law consisting of a pre-existing BH powerlaw mass spectrum and a harder powerlaw mass spectrum resulting from mergers. The predicted spin spectrum is multi-peaked with the evolution of retrograde spin BH in the gas disk playing a key role. We outline the large uncertainties in each of these LIGO observables for this channel and we discuss ways in which they can be constrained in the future.

Supermassive black hole binaries (SMBHBs) should be common in galactic nuclei as a result of frequent galaxy mergers. Recently, a large sample of sub-parsec SMBHB candidates was identified as bright periodically variable quasars in optical surveys. If the observed periodicity corresponds to the redshifted binary orbital period, the inferred orbital velocities are relativistic (v/c~0.1). The optical and UV luminosities are expected to arise from gas bound to the individual BHs, and would be modulated by the relativistic Doppler effect. The optical and UV light curves should vary in tandem with relative amplitudes which depend on the respective spectral slopes. We constructed a control sample of 42 quasars with aperiodic variability, to test whether this Doppler colour signature can be distinguished from intrinsic chromatic variability. We found that the Doppler signature can arise by chance in ~20% (~37%) of quasars in the nUV (fUV) band. These probabilities reflect the limited quality of the control sample and represent upper limits on how frequently quasars mimic the Doppler brightness+colour variations. We performed separate tests on the periodic quasar candidates, and found that for the majority, the Doppler boost hypothesis requires an unusually steep UV spectrum or an unexpectedly large BH mass and orbital velocity. We conclude that at most ~1/3rd of these periodic candidates can harbor Doppler-modulated SMBHBs.

We have developed a semi-analytic framework to model the large-scale evolution of the first Population III (Pop III) stars and the transition to metal-enriched star formation. Our model follows dark matter halos from cosmological N-body simulations, utilizing their individual merger histories and three-dimensional positions, and applies physically motivated prescriptions for star formation and feedback from Lyman-Werner (LW) radiation, hydrogen ionizing radiation, and external metal enrichment due to supernovae winds. This method is intended to complement analytic studies, which do not include clustering or individual merger histories, and hydrodynamical cosmological simulations, which include detailed physics, but are computationally expensive and have limited dynamic range. Utilizing this technique, we compute the cumulative Pop III and metal-enriched star formation rate density (SFRD) as a function of redshift at $z \geq 20$. We find that varying the model parameters leads to significant qualitative changes in the global star formation history. The Pop III star formation efficiency and the delay time between Pop III and subsequent metal-enriched star formation are found to have the largest impact. The effect of clustering (i.e. including the three-dimensional positions of individual halos) on various feedback mechanisms is also investigated. The impact of clustering on LW and ionization feedback is found to be relatively mild in our fiducial model, but can be larger if external metal enrichment can promote metal-enriched star formation over large distances.

We present the results of 2D, moving-mesh, viscous hydrodynamical simulations of an accretion disc around a merging supermassive black hole binary (SMBHB). The simulation is pseudo-Newtonian, with the BHs modeled as point masses with a Paczynski-Wiita potential, and includes viscous heating, shock heating, and radiative cooling. We follow the gravitational inspiral of an equal-mass binary with a component mass M_bh=10^6 M_sun from an initial separation of 60r_g (where r_g=GM_bh/c^2 is the gravitational radius) to the merger. We find that a central, low-density cavity forms around the binary, as in previous work, but that the BHs capture gas from the circumbinary disc and accrete efficiently via their own minidiscs, well after their inspiral outpaces the viscous evolution of the disc. The system remains luminous, displaying strong periodicity at twice the binary orbital frequency throughout the entire inspiral process, all the way to the merger. In the soft X-ray band, the thermal emission is dominated by the inner edge of the circumbinary disc with especially clear periodicity in the early inspiral. By comparison, harder X-ray emission is dominated by the minidiscs, and the light curve is initially more noisy but develops a clear periodicity in the late inspiral stage. This variability pattern should help identify the EM counterparts of SMBHBs detected by the space-based gravitational-wave detector LISA.

The advent of time domain astronomy is revolutionizing our understanding of the Universe. Programs such as the Catalina Real-time Transient Survey (CRTS) or the Palomar Transient Factory (PTF) surveyed millions of objects for several years, allowing variability studies on large statistical samples. The inspection of $\approx$250k quasars in CRTS resulted in a catalogue of 111 potentially periodic sources, put forward as supermassive black hole binary (SMBHB) candidates. A similar investigation on PTF data yielded 33 candidates from a sample of $\approx$35k quasars. Working under the SMBHB hypothesis, we compute the implied SMBHB merger rate and we use it to construct the expected gravitational wave background (GWB) at nano-Hz frequencies, probed by pulsar timing arrays (PTAs). After correcting for incompleteness and assuming virial mass estimates, we find that the GWB implied by the CRTS sample exceeds the current most stringent PTA upper limits by almost an order of magnitude. After further correcting for the implicit bias in virial mass measurements, the implied GWB drops significantly but is still in tension with the most stringent PTA upper limits. Similar results hold for the PTF sample. Bayesian model selection shows that the null hypothesis (whereby the candidates are false positives) is preferred over the binary hypothesis at about $2.3\sigma$ and $3.6\sigma$ for the CRTS and PTF samples respectively. Although not decisive, our analysis highlights the potential of PTAs as astrophysical probes of individual SMBHB candidates and indicates that the CRTS and PTF samples are likely contaminated by several false positives.

We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10. We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.

The non-zero mass of neutrinos suppresses the growth of cosmic structure on small scales. Since the level of suppression depends on the sum of the masses of the three active neutrino species, the evolution of large-scale structure is a promising tool to constrain the total mass of neutrinos and possibly shed light on the mass hierarchy. In this work, we investigate these effects via a large suite of N-body simulations that include massive neutrinos using an analytic linear-response approximation: the Cosmological Massive Neutrino Simulations (MassiveNuS). The simulations include the effects of radiation on the background expansion, as well as the clustering of neutrinos in response to the nonlinear dark matter evolution. We allow three cosmological parameters to vary: the neutrino mass sum M_nu in the range of 0-0.6 eV, the total matter density Omega_m, and the primordial power spectrum amplitude A_s. The rms density fluctuation in spheres of 8 comoving Mpc/h (sigma_8) is a derived parameter as a result. Our data products include N-body snapshots, halo catalogues, merger trees, ray- traced galaxy lensing convergence maps for four source redshift planes between z_s=1-2.5, and ray-traced cosmic microwave background lensing convergence maps. We describe the simulation procedures and code validation in this paper. The data are publicly available at http://columbialensing.org.

Multi-frequency gravitational wave (GW) observations are useful probes of the formation processes of coalescing stellar-mass binary black holes (BBHs). We discuss the phase drift in the GW inspiral waveform of the merging BBH caused by its center-of-mass acceleration. The acceleration strongly depends on the location where a BBH forms within a galaxy, allowing observations of the early inspiral phase of LIGO-like BBH mergers by the Laser Interferometer Space Antenna (LISA) to test the formation mechanism. In particular, BBHs formed in dense nuclear star clusters or via compact accretion disks around a nuclear supermassive black hole in active galactic nuclei would suffer strong acceleration, and produce large phase drifts measurable by LISA. The host galaxies of the coalescing BBHs in these scenarios can also be uniquely identified in the LISA error volume, without electromagnetic counterparts. A non-detection of phase drifts would rule out or constrain the contribution of the nuclear formation channels to the stellar-mass BBH population.

Using few-body simulations, we investigate the evolution of supermassive black holes (SMBHs) in galaxies ($M_{\star}=10^{10}-10^{12}{\rm M}_{\odot}$ at $z=0$) at $0<z<4$. Following galaxy merger trees from the Millennium simulation, we model BH mergers with two extreme binary decay scenarios for the `hard binary' stage: a full or an empty loss cone. These two models should bracket the true evolution, and allow us to separately explore the role of dynamical friction and that of multi-body BH interactions on BH mergers. Using the computed merger rates, we infer the stochastic gravitational wave background (GWB). Our dynamical approach is a first attempt to study the dynamical evolution of multiple SMBHs in the host galaxies undergoing mergers with various mass ratios ($10^{-4} < q_{\star} < 1$). Our main result demonstrates that SMBH binaries are able to merge in both scenarios. In the empty loss cone case, we find that BHs merge via multi-body interactions, avoiding the `final parsec' problem, and entering the PTA band with substantial orbital eccentricity. Our full loss cone treatment, albeit more approximate, suggests that the eccentricity becomes even higher when GWs become dominant, leading to rapid coalescences (binary lifetime $\lesssim1 {\rm ~Gyr}$). Despite the lower merger rates in the empty loss cone case, due to their higher mass ratios and lower redshifts, the GWB in the full/empty loss cone models are comparable ($0.70\times10^{-15}$ and $0.53\times10^{-15}$ at a frequency of $1~{\rm yr}^{-1}$, respectively). Finally, we compute the effects of high eccentricities on the GWB spectrum.

In symmetric gravitating systems experiencing rapid mass loss, particle orbits change almost instantaneously, which can lead to the development of a sharply contoured density profile, including singular caustics for collisionless systems. This framework can be used to model a variety of dynamical systems, such as accretion disks following a massive black hole merger and dwarf galaxies following violent early star formation feedback. Particle interactions in the high-density peaks seem a promising source of observable signatures of these mass loss events (i.e. a possible EM counterpart for black hole mergers or strong gamma-ray emission from dark matter annihilation around young galaxies), because the interaction rate depends on the square of the density. We study post-mass-loss density profiles, both analytic and numerical, in idealised cases and present arguments and methods to extend to any general system. An analytic derivation is presented for particles on Keplerian orbits responding to a drop in the central mass. We argue that this case, with initially circular orbits, gives the most sharply contoured profile possible. We find that despite the presence of a set of singular caustics, the total particle interaction rate is reduced compared to the unperturbed system; this is a result of the overall expansion of the system dominating over the steep caustics. Finally we argue that this result holds more generally, and the loss of central mass decreases the particle interaction rate in any physical system.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

The optical and UV emission from sub-parsec massive black hole binaries (MBHBs) in active galactic nuclei (AGN) is believed to vary periodically, on timescales comparable to the binary's orbital time. If driven by accretion rate fluctuations, the variability could be isotropic. If dominated by relativistic Doppler modulation, the variability should instead be anisotropic, resembling a rotating forward-beamed lighthouse. We consider the infrared (IR) reverberation of either type of periodic emission by pc-scale circumbinary dust tori. We predict the phase and amplitude of IR variability as a function of the ratio of dust light crossing time to the source variability period, and of the torus inclination and opening angle. We enumerate several differences between the isotropic and anisotropic cases. Interestingly, for a nearly face-on binary with an inclined dust torus, the Doppler boost can produce IR variability without any observable optical/UV variability. Such orphan-IR variability would have been missed in optical searches for periodic AGN. We apply our models to time-domain WISE IR data from the MBHB candidate PG 1302-102 and find consistency with dust reverberation by both isotropically emitting and Doppler-boosted sources in the shorter wavelength W1-W2 (2.8-5.3 micrometer) bands. We constrain the dust torus to be thin (aspect ratio ~0.1), with an inner radius at 1-5 pc. More generally, our dust echo models will aid in identifying new MBHB candidates, determining their nature, and constraining the physical properties of MBHBs and their dust tori.

We explore the sensitivity of weak lensing observables to the expansion history of the universe and to the growth of cosmic structures, as well as the relative contribution of both effects to constraining cosmological parameters. We utilize ray-tracing dark-matter-only N-body simulations and validate our technique by comparing our results for the convergence power spectrum with analytic results from past studies. We then extend our analysis to non-Gaussian observables which cannot be easily treated analytically. We study the convergence (equilateral) bispectrum and two topological observables, lensing peaks and Minkowski functionals, focusing on their sensitivity to the matter density $\Omega_m$ and the dark energy equation of state $w$. We find that a cancelation between the geometry and growth effects is a common feature for all observables, and exists at the map level. It weakens the overall sensitivity by up to a factor of 3 and 1.5 for $w$ and $\Omega_m$, respectively, with the bispectrum worst affected. However, combining geometry and growth information alleviates the degeneracy between $\Omega_m$ and $w$ from either effect alone. As a result, the magnitude of marginalized errors remain similar to those obtained from growth-only effects, but with the correlation between the two parameters switching sign. These results shed light on the origin of cosmology-sensitivity of non-Gaussian statistics, and should be useful in optimizing combinations of observables.

The gravitational waves (GWs) from a binary black hole (BBH) with masses between 10^4 and 10^7 Msun can be detected with the Laser Interferometer Space Antenna (LISA) once their orbital frequency exceeds 10^-4 - 10^-5 Hz. The binary separation at this stage is approximately a=100 R_g (gravitational radius), and the orbital speed is of order v/c=0.1. We argue that at this stage, the binary will be producing bright electromagnetic (EM) radiation via gas bound to the individual BHs. Both BHs will have their own photospheres in X-ray and possibly also in optical bands. Relativistic Doppler modulations and lensing effects will inevitably imprint periodic variability in the EM light-curve, tracking the phase of the orbital motion, and serving as a template for the GW inspiral waveform. Advanced localization of the source by LISA weeks to months prior to merger will enable a measurement of this EM chirp by wide-field X-ray or optical instruments. A comparison of the phases of the GW and EM chirp signals will help break degeneracies between system parameters, and probe a fractional difference difference Delta v in the propagation speed of photons and gravitons as low as Delta v/c = O(10^-17).

Accurate forward modeling of weak lensing (WL) observables from cosmological parameters is necessary for upcoming galaxy surveys. Because WL probes structures in the non-linear regime, analytical forward modeling is very challenging, if not impossible. Numerical simulations of WL features rely on ray-tracing through the outputs of $N$-body simulations, which requires knowledge of the gravitational potential and accurate solvers for light ray trajectories. A less accurate procedure, based on the Born approximation, only requires knowledge of the density field, and can be implemented more efficiently and at a lower computational cost. In this work, we use simulations to show that deviations of the Born-approximated convergence power spectrum, skewness and kurtosis from their fully ray--traced counterparts are consistent with the smallest non-trivial $O(\Phi^3)$ post-Born corrections (so-called geodesic and lens-lens terms). Our results imply a cancellation among the larger $O(\Phi^4)$ (and higher order) terms, consistent with previous analytic work. We also find that cosmological parameter bias induced by the Born approximated power spectrum is negligible even for an LSST-like survey, once galaxy shape noise is considered. When considering higher order statistics such as the $\kappa$ skewness and kurtosis, however, we find significant bias of up to 2.5$\sigma$. Using the LensTools software suite, we show that the Born approximation saves a factor of 4 in computing time with respect to the full ray-tracing in reconstructing the convergence.

We utilize cosmological hydrodynamic simulations to study the formation of Population III (Pop III) stars in dark matter halos exposed to strong ionizing radiation. We simulate the formation of three halos subjected to a wide range of ionizing fluxes, and find that for high flux, ionization and photoheating can delay gas collapse and star formation up to halo masses significantly larger than the atomic cooling threshold. The threshold halo mass at which gas first collapses and cools increases with ionizing flux for intermediate values, and saturates at a value approximately an order of magnitude above the atomic cooling threshold for extremely high flux (e.g. $\approx 5 \times 10^8 ~ M_\odot$ at $z\approx6$). This behavior can be understood in terms of photoheating, ionization/recombination, and Ly$\alpha$ cooling in the pressure-supported, self-shielded gas core at the center of the growing dark matter halo. We examine the spherically-averaged radial velocity profiles of collapsing gas and find that a gas mass of up to $\approx 10^{6}~ M_\odot$ can reach the central regions within $3~{\rm Myr}$, providing an upper limit on the amount of massive Pop III stars that can form. The ionizing radiation increases this limit by a factor of a few compared to strong Lyman-Werner (LW) radiation alone. We conclude that the bright HeII 1640 \AA\ emission recently observed from the high-redshift galaxy CR7 cannot be explained by Pop III stars alone. However, in some halos, a sufficient number of Pop III stars may form to be detectable with future telescopes such as the James Webb Space Telescope (JWST).

Gaseous circumbinary accretion discs provide a promising mechanism to facilitate the mergers of supermassive black holes (SMBHs) in galactic nuclei. We measure the torques exerted on accreting SMBH binaries, using 2D, isothermal, moving-mesh, viscous hydrodynamical simulations of circumbinary accretion discs. Our computational domain includes the entire inner region of the circumbinary disk with the individual black holes (BHs) included as point masses on the grid and a sink prescription to model accretion onto each BH. The BHs each acquire their own well-resolved accretion discs ("minidiscs"). We explore a range of mass removal rates for the sink prescription removing gas from the central regions of the minidiscs. We find that the torque exerted on the binary is primarily gravitational, and dominated by the gas orbiting close behind and ahead of the individual BHs. The torques from the distorted circumbinary disc farther out and from the direct accretion of angular momentum are subdominant. The torques are sensitive to the sink prescription: slower sinks result in more gas accumulating near the BHs and more negative torques, driving the binary to merger more rapidly. For faster sinks, the torques are less negative and eventually turn positive (for unphysically fast sinks). When the minidiscs are modeled as standard alpha discs, our results are insensitive to the choice of sink radius. Scaling the simulations to a binary orbital period tbin = 1yr and background disc accretion rate Mdot = 0.3MEdd in Eddington units, the binary inspirals on a timescale of 3X10^6 years, irrespective of the SMBH masses. For binaries with total mass <10^7Msun, this is shorter than the inspiral time due to gravitational wave (GW) emission alone, implying that gas discs will have a significant impact on the SMBH binary population and can affect the GW signal for Pulsar Timing Arrays.

The Direct Collapse Black Hole (DCBH) scenario provides a solution for forming the massive black holes powering bright quasars observed in the early Universe. A prerequisite for forming a DCBH is that the formation of (much less massive) Population III stars be avoided - this can be achieved by destroying H$_2$ via Lyman-Werner (LW) radiation (E$_{\rm{LW}}$ = 12.6 eV). We find that two conditions must be met in the proto-galaxy that will host the DCBH. First, prior star formation must be delayed; this can be achieved with a background LW flux of J$_{\rm BG} \gtrsim 100\ J_{21}$. Second, an intense burst of LW radiation from a neighbouring star-bursting proto-galaxy is required, just before the gas cloud undergoes gravitational collapse, to finally suppress star formation completely. We show here for the first time using high-resolution hydrodynamical simulations, including full radiative transfer, that this low-level background, combined with tight synchronisation and irradiation of a secondary proto-galaxy by a primary proto-galaxy, inevitably moves the secondary proto-galaxy onto the isothermal atomic cooling track, without the deleterious effects of either photo-evaporating the gas or polluting it by heavy elements. These, atomically cooled, massive proto-galaxies are expected to ultimately form a DCBH of mass $10^4 - 10^5 M_{\odot}$.

Suppression of H2-cooling in early protogalaxies has important implications for the formation of supermassive black holes seeds, the first generation of stars, and the epoch of reionization. This suppression can occur via photodissociation of H2 (by ultraviolet Lyman-Werner [LW] photons) or by photodetachment of H-, a precursor in H2 formation (by infrared [IR] photons). Previous studies have typically adopted idealised spectra, with a blackbody or a power-law shape, in modeling the chemistry of metal-free protogalaxies, and utilised a single parameter, the critical UV flux, or Jcrit, to determine whether H2-cooling is prevented. This can be misleading, and that independent of the spectral shape, there is a critical curve in the (kLW , kH^- ) plane, where kLW and kH^- are the H2-dissocation rates by LW and IR photons, which determines whether a protogalaxy can cool below ~1000 Kelvin. We use a one-zone model to follow the chemical and thermal evolution of gravitationally collapsing protogalactic gas, to compute this critical curve, and provide an accurate analytical fit for it. We improve on previous works by considering a variety of more realistic Pop III or Pop II-type spectra from population synthesis models and perform fully frequency-dependent calculations of the H2-photodissociation rates for each spectrum. We compute the ratio kLW/kH^- for each spectrum, as well as the minimum stellar mass M_star, for various IMFs and metallicities, required to prevent cooling in a neighboring halo a distance d away. We provide critical M_star/d2 values for suppression of H2-cooling, with analytic fits, which can be used in future studies.

Stellar-mass binary black hole mergers are poised to represent the majority of gravitational-wave (GW) observations by Advanced LIGO and Virgo. Probing their origin will be difficult due to the expected lack of electromagnetic emission and limited localization accuracy. Associations with rare host galaxy types -- such as active galactic nuclei (AGN) -- can nevertheless be identified statistically through spatial correlation. We show that (i) fractional contributions $f_{\rm agn}=50-100\%$ from AGN hosts to the total BBH merger rate can be statistically established with 70-300 detected events (expected in 0.5-2 years of observation with Advanced LIGO-Virgo at design sensitivity and current rate estimates); (ii) fractional contributions as low as $f_{\rm agn}=25\%$ can be tested with 1000 events ($\sim$ 5\,years of observation); (iii) the $\sim5\%$ best localized GWs drive these constraints. The presented method and results are generally applicable to binary formation channels with rare host populations.

We explore the evolution of stellar mass black hole binaries (BHBs) which are formed in the self-gravitating disks of active galactic nuclei (AGN). Hardening due to three-body scattering and gaseous drag are effective mechanisms that reduce the semi-major axis of a BHB to radii where gravitational waves take over, on timescales shorter than the typical lifetime of the AGN disk. Taking observationally-motivated assumptions for the rate of star formation in AGN disks, we find a rate of disk-induced BHB mergers ($\mathcal{R} \sim 3~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, but with large uncertainties) that is comparable with existing estimates of the field rate of BHB mergers, and the approximate BHB merger rate implied by the recent Advanced LIGO detection of GW150914. BHBs formed thorough this channel will frequently be associated with luminous AGN, which are relatively rare within the sky error regions of future gravitational wave detector arrays. This channel could also possess a (potentially transient) electromagnetic counterpart due to super-Eddington accretion onto the stellar mass black hole following the merger.

We introduce an intrinsic Ly\alpha\ emission line profile reconstruction method for high-$z$ quasars (QSOs). This approach utilises a covariance matrix of emission line properties obtained from a large, moderate-$z$ ($2 \leq z \leq 2.5$), high signal to noise (S/N > 15) sample of BOSS QSOs. For each QSO, we complete a Monte Carlo Markov Chain fitting of the continuum and emission line properties and perform a visual quality assessment to construct a large database of robustly fit spectra. With this dataset, we construct a covariance matrix to describe the correlations between the high ionisation emission lines Ly\alpha, C IV, Si IV + O IV] and C III], and find it to be well approximated by an $N$-dimensional Gaussian distribution. This covariance matrix characterises the correlations between the line width, peak height and velocity offset from systemic while also allowing for the existence of broad and narrow line components for Ly\alpha\ and C IV. We illustrate how this covariance matrix allows us to statistically characterise the intrinsic Ly\alpha\ line solely from the observed spectrum redward of 1275\AA. This procedure can be used to reconstruct the intrinsic Ly\alpha\ line emission profile in cases where Ly\alpha\ may otherwise be obscured. Applying this reconstruction method to our sample of QSOs, we recovered the Ly\alpha\ line flux to within 15 per cent of the measured flux at 1205\AA\ (1220\AA) ~85 (90) per cent of the time.

The Laser Interferometer Gravitational-Wave Observatory, LIGO, found direct evidence for double black hole binaries emitting gravitational waves. Galactic nuclei are expected to harbor the densest population of stellar-mass black holes. A significant fraction ($\sim30\%$) of these black holes can reside in binaries. We examine the fate of the black hole binaries in active galactic nuclei, which get trapped in the inner region of the accretion disk around the central supermassive black hole. We show that binary black holes can migrate into and then rapidly merge within the disk well within a Salpeter time. The binaries may also accrete a significant amount of gas from the disk, well above the Eddington rate. This could lead to detectable X-ray or gamma-ray emission, but would require hyper-Eddington accretion with a few percent radiative efficiency, comparable to thin disks. We discuss implications for gravitational wave observations and black hole population studies. We estimate that Advanced LIGO may detect $\sim20$ such, gas-induced binary mergers per year.

Weak lensing convergence peaks are a promising tool to probe nonlinear structure evolution at late times, providing additional cosmological information beyond second-order statistics. Previous theoretical and observational studies have shown that the cosmological constraints on $\Omega_m$ and $\sigma_8$ are improved by a factor of up to ~ 2 when peak counts and second-order statistics are combined, compared to using the latter alone. We study the origin of lensing peaks using observational data from the 154 deg$^2$ Canada-France-Hawaii Telescope Lensing Survey. We found that while high peaks (with height $\kappa$ >3.5 $\sigma_\kappa$, where $\sigma_\kappa$ is the r.m.s. of the convergence $\kappa$) are typically due to one single massive halo of ~$10^{15}M_\odot$, low peaks ($\kappa$ <~ $\sigma_\kappa$) are associated with constellations of 2-8 smaller halos (<~$10^{13}M_\odot$). In addition, halos responsible for forming low peaks are found to be significantly offset from the line-of-sight towards the peak center (impact parameter >~ their virial radii), compared with ~0.25 virial radii for halos linked with high peaks, hinting that low peaks are more immune to baryonic processes whose impact is confined to the inner regions of the dark matter halos. Our findings are in good agreement with results from the simulation work by Yang el al. (2011).

Unprecedentedly precise cosmic microwave background (CMB) data are expected from ongoing and near-future CMB Stage-III and IV surveys, which will yield reconstructed CMB lensing maps with effective resolution approaching several arcminutes. The small-scale CMB lensing fluctuations receive non-negligible contributions from nonlinear structure in the late-time density field. These fluctuations are not fully characterized by traditional two-point statistics, such as the power spectrum. Here, we use $N$-body ray-tracing simulations of CMB lensing maps to examine two higher-order statistics: the lensing convergence one-point probability distribution function (PDF) and peak counts. We show that these statistics contain significant information not captured by the two-point function, and provide specific forecasts for the ongoing Stage-III Advanced Atacama Cosmology Telescope (AdvACT) experiment. Considering only the temperature-based reconstruction estimator, we forecast 9$\sigma$ (PDF) and 6$\sigma$ (peaks) detections of these statistics with AdvACT. Our simulation pipeline fully accounts for the non-Gaussianity of the lensing reconstruction noise, which is significant and cannot be neglected. Combining the power spectrum, PDF, and peak counts for AdvACT will tighten cosmological constraints in the $\Omega_m$-$\sigma_8$ plane by $\approx 30\%$, compared to using the power spectrum alone.

Weak gravitational lensing is becoming a mature technique for constraining cosmological parameters, and future surveys will be able to constrain the dark energy equation of state $w$. When analyzing galaxy surveys, redshift information has proven to be a valuable addition to angular shear correlations. We forecast parameter constraints on the triplet $(\Omega_m,w,\sigma_8)$ for an LSST-like photometric galaxy survey, using tomography of the shear-shear power spectrum, convergence peak counts and higher convergence moments. We find that redshift tomography with the power spectrum reduces the area of the $1\sigma$ confidence interval in $(\Omega_m,w)$ space by a factor of 8 with respect to the case of the single highest redshift bin. We also find that adding non-Gaussian information from the peak counts and higher-order moments of the convergence field and its spatial derivatives further reduces the constrained area in $(\Omega_m,w)$ by a factor of 3 and 4, respectively. When we add cosmic microwave background parameter priors from Planck to our analysis, tomography improves power spectrum constraints by a factor of 3. Adding moments yields an improvement by an additional factor of 2, and adding both moments and peaks improves by almost a factor of 3, over power spectrum tomography alone. We evaluate the effect of uncorrected systematic photometric redshift errors on the parameter constraints. We find that different statistics lead to different bias directions in parameter space, suggesting the possibility of eliminating this bias via self-calibration.

We have investigated a recently proposed halo-based model, Camelus, for predicting weak-lensing peak counts, and compared its results over a collection of 162 cosmologies with those from N-body simulations. While counts from both models agree for peaks with $\mathcal{S/N}>1$ (where $\mathcal{S/N}$ is the ratio of the peak height to the r.m.s. shape noise), we find $\approx 50\%$ fewer counts for peaks near $\mathcal{S/N}=0$ and significantly higher counts in the negative $\mathcal{S/N}$ tail. Adding shape noise reduces the differences to within $20\%$ for all cosmologies. We also found larger covariances that are more sensitive to cosmological parameters. As a result, credibility regions in the $\{\Omega_m, \sigma_8\}$ are $\approx 30\%$ larger. Even though the credible contours are commensurate, each model draws its predictive power from different types of peaks. Low peaks, especially those with $2<\mathcal{S/N}<3$, convey important cosmological information in N-body data, as shown in \cite{DietrichHartlap, Kratochvil2010}, but \textsc{Camelus} constrains cosmology almost exclusively from high significance peaks $(\mathcal{S/N}>3)$. Our results confirm the importance of using a cosmology-dependent covariance with at least a 14\% improvement in parameter constraints. We identified the covariance estimation as the main driver behind differences in inference, and suggest possible ways to make Camelus even more useful as a highly accurate peak count emulator.

Supermassive black hole binaries (SMBHBs) at sub-parsec separations should be common in galactic nuclei, as a result of frequent galaxy mergers. Hydrodynamical simulations of circumbinary discs predict strong periodic modulation of the mass accretion rate on time-scales comparable to the orbital period of the binary. As a result, SMBHBs may be recognized by the periodic modulation of their brightness. We conducted a statistical search for periodic variability in a sample of 35,383 spectroscopically confirmed quasars in the photometric database of the Palomar Transient Factory (PTF). We analysed Lomb-Scargle periodograms and assessed the significance of our findings by modeling each individual quasar's variability as a damped random walk (DRW). We identified 50 quasars with significant periodicity beyond the DRW model, typically with short periods of a few hundred days. We find 33 of these to remain significant after a re-analysis of their periodograms including additional optical data from the intermediate-PTF and the Catalina Real-Time Transient Survey (CRTS). Assuming that the observed periods correspond to the redshifted orbital periods of SMBHBs, we conclude that our findings are consistent with a population of unequal-mass SMBHBs, with a typical mass ratio as low as q = M2/M1 ~ 0.01.

We perform one-dimensional radiation hydrodynamical simulations to solve accretion flows onto massive black holes (BHs) with a very high rate. Assuming that photon trapping limits the luminosity emerging from the central region to $L\lesssim L_{\rm Edd}$, Inayoshi, Haiman & Ostriker (2016) have shown that an accretion flow settles to a "hyper-Eddington" solution, with a steady and isothermal ($T\simeq 8000$ K) Bondi profile reaching $\gtrsim 5000$ times the Eddington accretion rate $\dot{M}_{\rm Edd}\equiv L_{\rm Edd}/c^2$. Here we address the possibility that gas accreting with finite angular momentum forms a bright nuclear accretion disc, with a luminosity exceeding the Eddington limit ($1\lesssim L/L_{\rm Edd}\lesssim 100$). Combining our simulations with an analytic model, we find that a transition to steady hyper-Eddington accretion still occurs, as long as the luminosity remains below $L/L_{\rm Edd}\lesssim 35~(M_{\rm BH}/10^4~M_\odot)^{3/2}(n_\infty/10^5~{\rm cm^{-3}})(T_\infty/10^4~{\rm K})^{-3/2}(r_{\star}/10^{14}~{\rm cm})^{-1/2}$, where $n_\infty$ and $T_\infty$ are the density and temperature of the ambient gas, and $r_\star$ is the radius of the photosphere, at which radiation emerges. If the luminosity exceeds this value, accretion becomes episodic. Our results can be accurately recovered in a toy model of an optically thick spherical shell, driven by radiation force into a collapsing medium. When the central source is dimmer than the above critical value, the expansion of the shell is halted and reversed by ram pressure of the collapsing medium, and by shell's weight. Our results imply that rapid, unimpeded hyper-Eddington accretion is possible even if the luminosity of the central source far exceeds the Eddington limit, and can be either steady or strongly episodic.

The recent discovery of the gravitational wave source GW150914 has revealed a coalescing binary black hole (BBH) with masses of $\sim 30~M_\odot$. Previous proposals for the origin of such a massive binary include Population III (PopIII) stars. PopIII stars are efficient producers of BBHs and of a gravitational wave background (GWB) in the $10-100$ Hz band, and also of ionizing radiation in the early Universe. We quantify the relation between the amplitude of the GWB ($\Omega_{\rm gw}$) and the electron scattering optical depth ($\tau_{\rm e}$), produced by PopIII stars, assuming that $f_{\rm esc}\approx 10\%$ of their ionizing radiation escapes into the intergalactic medium. We find that PopIII stars would produce a GWB that is detectable by the future O5 LIGO/Virgo if $\tau_{\rm e} \gtrsim 0.07$, consistent with the recent Planck measurement of $\tau_e=0.055 \pm 0.09$. Moreover, the spectral index of the background from PopIII BBHs becomes as small as ${\rm d}\ln \Omega_{\rm gw}/{\rm d}\ln f\lesssim 0.3$ at $f \gtrsim 30$ Hz, which is significantly flatter than the value $\sim 2/3$ generically produced by lower-redshift and less-massive BBHs. A detection of the unique flattening at such low frequencies by the O5 LIGO/Virgo will indicate the existence of a high-chirp mass, high-redshift BBH population, which is consistent with the PopIII origin. A precise characterization of the spectral shape near $30-50$ Hz by the Einstein Telescope could also constrain the PopIII initial mass function and star formation rate.

The largest observed supermassive black holes (SMBHs) have a mass of M_BH ~ 10^{10} M_sun, nearly independent of redshift, from the local (z~0) to the early (z>6) Universe. We suggest that the growth of SMBHs above a few 10^{10} M_sun is prevented by small-scale accretion physics, independent of the properties of their host galaxies or of cosmology. Growing more massive BHs requires a gas supply rate from galactic scales onto a nuclear region as high as >10^3 M_sun/yr. At such a high accretion rate, most of the gas converts to stars at large radii (~10-100 pc), well before reaching the BH. We adopt a simple model (Thompson et al. 2005) for a star-forming accretion disk, and find that the accretion rate in the sub-pc nuclear region is reduced to the smaller value of at most a few M_sun/yr. This prevents SMBHs from growing above ~10^{11} M_sun in the age of the Universe. Furthermore, once a SMBH reaches a sufficiently high mass, this rate falls below the critical value at which the accretion flow becomes advection dominated. Once this transition occurs, BH feeding can be suppressed by strong outflows and jets from hot gas near the BH. We find that the maximum SMBH mass, given by this transition, is between M_{BH,max} ~ (1-6) * 10^{10} M_sun, depending primarily on the efficiency of angular momentum transfer inside the galactic disk, and not on other properties of the host galaxy.

We report the formation of intermediate-mass black holes (IMBHs) in suites of numerical $N$-body simulations of Population III remnant black holes (BHs) embedded in gas-rich protogalaxies at redshifts $z\gtrsim10$. We model the effects of gas drag on the BHs' orbits, and allow BHs to grow via gas accretion, including a mode of hyper-Eddington accretion in which photon trapping and rapid gas inflow suppress any negative radiative feedback. Most initial BH configurations lead to the formation of one (but never more than one) IMBH in the center of the protogalaxy, reaching a mass of $10^{3-5}\mathrm{M}_{\odot}$ through hyper-Eddington growth. Our results suggest a viable pathway to forming the earliest massive BHs in the centers of early galaxies. We also find that the nuclear IMBH typically captures a stellar-mass BH companion, making these systems observable in gravitational waves as extreme mass-ratio inspirals (EMRIs) with \textit{eLISA}.

We study circumbinary accretion discs in the framework of the restricted three-body problem (R3Bp) and via numerically solving the height-integrated equations of viscous hydrodynamics. Varying the mass ratio of the binary, we find a pronounced change in the behaviour of the disc near mass ratio $q \equiv M_s/M_p \sim 0.04$. For mass ratios above $q=0.04$, solutions for the hydrodynamic flow transition from steady, to strongly-fluctuating; a narrow annular gap in the surface density around the secondary's orbit changes to a hollow central cavity; and a spatial symmetry is lost, resulting in a lopsided disc. This phase transition is coincident with the mass ratio above which stable orbits do not exist around the L4 and L5 equilibrium points of the R3B problem. Using the DISCO code, we find that for thin discs, for which a gap or cavity can remain open, the mass ratio of the transition is relatively insensitive to disc viscosity and pressure. The $q=0.04$ transition has relevance for the evolution of massive black hole binary+disc systems at the centers of galactic nuclei, as well as for young stellar binaries and possibly planets around brown dwarfs.

We explore the formation of massive high-redshift Population III (Pop III) galaxies through photoionization feedback. We consider dark matter halos formed from progenitors that have undergone no star formation as a result of early reionization and photoevaporation caused by a nearby galaxy. Once such a halo reaches $\approx 10^9~M_\odot$, corresponding to the Jeans mass of the photoheated intergalactic medium (IGM) at $z\approx 7$, pristine gas is able to collapse into the halo, potentially producing a massive Pop III starburst. We suggest that this scenario may explain the recent observation of strong He II 1640~\AA~line emission in CR7, which is consistent with $\sim 10^7~M_\odot$ of young Pop III stars. Such a large mass of Pop III stars is unlikely without the photoionization feedback scenario, because star formation is expected to inject metals into halos above the atomic cooling threshold ($\sim 10^8~M_\odot$ at $z \approx 7$). We use merger trees to analytically estimate the abundance of observable Pop III galaxies formed through this channel, and find a number density of $\approx 10^{-7}~{\rm Mpc^{-3}}$ at $z=6.6$ (the redshift of CR7). This is approximately a factor of ten lower than the density of Ly$\alpha$ emitters as bright as CR7.

We study very-high rate spherically symmetric accretion flows onto a massive black hole (BH; 10^2 < M_BH < 10^6 Msun) embedded in a dense gas cloud with a low abundance of metals, performing one-dimensional hydrodynamical simulations which include multi-frequency radiation transfer and non-equilibrium primordial chemistry. We find that rapid gas supply from the Bondi radius at a hyper-Eddington rate can occur without being impeded by radiation feedback when (n/10^5 cm^-3) > (M_BH/10^4Msun)^{-1}(T/10^4 K)^{3/2}, where n and T are the density and temperature of ambient gas outside of the Bondi radius. The resulting accretion rate in this regime is steady, and larger than 3000 times the Eddington rate. At lower Bondi rates, the accretion is episodic due to radiative feedback and the average rate is limited below the Eddington rate. For the hyper-Eddington case, the steady solution consists of two parts: a radiation-dominated central core, where photon trapping due to electron scattering is important, and an accreting envelope which follows a Bondi profile with T~8000 K. When the emergent luminosity is limited below the Eddington luminosity because of photon trapping, radiation from the central region does not affect the gas dynamics at larger scales. We apply our result to the rapid formation of massive BHs in protogalaxies with a virial temperature of T_vir> 10^4 K. Once a seed BH forms at the center of the galaxy, it can grow up to a maximum ~10^5 (T_vir/10^4 K) Msun via gas accretion independent of the initial BH mass. Finally, we discuss possible observational signatures of rapidly accreting BHs with/without allowance for dust. We suggest that these systems could explain Lya emitters without X-rays and luminous infrared sources with hot dust emission, respectively.

Constraining cosmology using weak gravitational lensing consists of comparing a measured feature vector of dimension $N_b$ with its simulated counterpart. An accurate estimate of the $N_b\times N_b$ feature covariance matrix $\mathbf{C}$ is essential to obtain accurate parameter confidence intervals. When $\mathbf{C}$ is measured from a set of simulations, an important question is how large this set should be. To answer this question, we construct different ensembles of $N_r$ realizations of the shear field, using a common randomization procedure that recycles the outputs from a smaller number $N_s\leq N_r$ of independent ray-tracing $N$--body simulations. We study parameter confidence intervals as a function of ($N_s,N_r$) in the range $1\leq N_s\leq 200$ and $1\leq N_r\lesssim 10^5$. Previous work has shown that Gaussian noise in the feature vectors (from which the covariance is estimated) lead, at quadratic order, to an $O(1/N_r)$ degradation of the parameter confidence intervals. Using a variety of lensing features measured in our simulations, including shear-shear power spectra and peak counts, we show that cubic and quartic covariance fluctuations lead to additional $O(1/N_r^2)$ error degradation that is not negligible when $N_r$ is only a factor of few larger than $N_b$. We study the large $N_r$ limit, and find that a single, 240Mpc$/h$ sized $512^3$-particle $N$--body simulation ($N_s=1$) can be repeatedly recycled to produce as many as $N_r={\rm few}\times10^4$ shear maps whose power spectra and high-significance peak counts can be treated as statistically independent. As a result, a small number of simulations ($N_s=1$ or $2$) is sufficient to forecast parameter confidence intervals at percent accuracy.

In this introductory chapter, we outline expectations for when and how the hydrogen and helium atoms in the universe turned from neutral to ionized, focusing on the earliest, least well understood stages, and emphasizing the most important open questions. We include a historical summary, and highlight the role of reionization as one of the few milestones in the evolution of the universe since the Big Bang, and its status as a unique probe of the beginning stages of structure formation.

We propose an observational test for gravitationally recoiling supermassive black holes (BHs) in active galactic nuclei, based on a correlation between the velocities of BHs relative to their host galaxies, |\Delta v|, and their obscuring dust column densities, \Sigma_{dust} (both measured along the line of sight). We use toy models for the distribution of recoil velocities, BH trajectories, and the geometry of obscuring dust tori in galactic centres, to simulate 2.5x10^5 random observations of recoiling quasars. BHs with recoil velocities comparable to the escape velocity from the galactic centre remain bound to the nucleus, and do not fully settle back to the centre of the torus due to dynamical friction in a typical quasar lifetime. We find that |\Delta v| and \Sigma_ {dust} for these BHs are positively correlated. For obscured (\Sigma_{dust}>0) and for partially obscured (0<\Sigma_{dust}<~2.3 g/m^2) quasars with |\Delta v|>=45 km/s, the sample correlation coefficient between log10(|\Delta v|) and \Sigma_{dust} is r_{45} = 0.28+/-0.02 and r_{45} = 0.13+/-0.02, respectively. Allowing for random +/-100 km/s errors in |\Delta v| unrelated to the recoil dilutes the correlation for the partially obscured quasars to r_{45} = 0.026+/-0.004 measured between |\Delta v| and \Sigma_{dust}. A random sample of >~3,500 obscured quasars with |\Delta v|>=45 km/s would allow rejection of the no-correlation hypothesis with 3 sigma significance 95% of the time. Finally, we find that the fraction of obscured quasars, F_{obs}(|\Delta v|), decreases with |\Delta v| from F_{obs}(<10 km/s)>~0.8 to F_{obs}(>10^3 km/s)<~0.4. This predicted trend can be compared to the observed fraction of type II quasars, and can further test combinations of recoil, trajectory, and dust torus models.

Lensing peaks have been proposed as a useful statistic, containing cosmological information from non-Gaussianities that is inaccessible from traditional two-point statistics such as the power spectrum or two-point correlation functions. Here we examine constraints on cosmological parameters from weak lensing peak counts, using the publicly available data from the 154 deg$^2$ CFHTLenS survey. We utilize a new suite of ray-tracing N-body simulations on a grid of 91 cosmological models, covering broad ranges of the three parameters $\Omega_m$, $\sigma_8$, and $w$, and replicating the Galaxy sky positions, redshifts, and shape noise in the CFHTLenS observations. We then build an emulator that interpolates the power spectrum and the peak counts to an accuracy of $\leq 5\%$, and compute the likelihood in the three-dimensional parameter space ($\Omega_m$, $\sigma_8$, $w$) from both observables. We find that constraints from peak counts are comparable to those from the power spectrum, and somewhat tighter when different smoothing scales are combined. Neither observable can constrain $w$ without external data. When the power spectrum and peak counts are combined, the area of the error "banana'' in the ($\Omega_m$, $\sigma_8$) plane reduces by a factor of $\approx2$, compared to using the power spectrum alone. For a flat $\Lambda$ cold dark matter model, combining both statistics, we obtain the constraint $\sigma_8(\Omega_m/0.27)^{0.63}=0.85\substack{+0.03 \\ -0.03}$.

Because most large galaxies contain a central black hole, and galaxies often merge, black-hole binaries are expected to be common in galactic nuclei. Although they cannot be imaged, periodicities in the light curves of quasars have been interpreted as evidence for binaries, most recently in PG~1302-102, with a short rest-frame optical period of 4 years. If the orbital period matches this value, then for the range of estimated black hole masses the components would be separated by 0.007-0.017 pc, implying relativistic orbital speeds. There has been much debate over whether black hole orbits could be smaller than 1 pc. Here we show that the amplitude and the sinusoid-like shape of the variability of PG~1302-102 can be fit by relativistic Doppler boosting of emission from a compact, steadily accreting, unequal-mass binary. We predict that brightness variations in the ultraviolet light curve track those in the optical, but with a 2-3 times larger amplitude. This prediction is relatively insensitive to the details of the emission process, and is consistent with archival UV data. Follow-up UV and optical observations in the next few years can test this prediction and confirm the existence of a binary black hole in the relativistic regime.

Recently, Planck measured a value of the cosmic microwave background (CMB) optical depth due to electron scattering of $\tau=0.066 \pm 0.016$. Here we show that this low value leaves essentially no room for an early partial reionisation of the intergalactic medium (IGM) by high-redshift Population III (Pop III) stars, expected to have formed in low-mass minihaloes. We perform semi-analytic calculations of reionisation which include the contribution from Pop II stars in atomic cooling haloes, calibrated with high-redshift galaxy observations, and Pop III stars in minihaloes with feedback due to Lyman-Werner (LW) radiation and metal enrichment. We find that without LW feedback or prompt metal enrichment (and assuming a minihalo escape fraction of 0.5) the Pop III star formation efficiency cannot exceed $\sim{\rm a~few}\times 10^{-4}$, without violating the constraints set by Planck data. This excludes massive Pop III star formation in typical $10^6 M_\odot$ minihaloes. Including LW feedback and metal enrichment alleviates this tension, allowing large Pop III stars to form early on before they are quenched by feedback. We find that the total density of Pop III stars formed across cosmic time is $\lesssim 10^{4-5}~M_\odot~{\rm Mpc^{-3}}$ and does not depend strongly on the feedback prescription adopted. Additionally, we perform a simple estimate of the possible impact on reionisation of X-rays produced by accretion onto black hole remnants of Pop III stars. We find that unless the accretion duty cycle is very low ($\lesssim 0.01$), this could lead to an optical depth inconsistent with Planck.

Graham et al. (2015) have detected a 5.2 year periodic optical variability of the quasar PG 1302-102 at redshift $z=0.3$, which they interpret as the redshifted orbital period $(1+z)t_{\rm bin}$ of a putative supermassive black hole binary (SMBHB). Here we consider the implications of a $3-8$ times shorter orbital period, suggested by hydrodynamical simulations of circumbinary discs (CBDs) with nearly equal--mass SMBHBs ($q\equiv M_2/M_1\gtrsim 0.3$). With the corresponding $2-4$ times tighter binary separation, PG 1302 would be undergoing gravitational wave dominated inspiral, and serve as a proof that the BHs can be fueled and produce bright emission even in this late stage of the merger. The expected fraction of binaries with the shorter $t_{\rm bin}$, among bright quasars, would be reduced by 1-2 orders of magnitude, compared to the 5.2 year period, in better agreement with the rarity of candidates reported by Graham et al. (2015). Finally, shorter periods would imply higher binary speeds, possibly imprinting periodicity on the light curves from relativistic beaming, as well as measurable relativistic effects on the Fe K$\alpha$ line. The CBD model predicts additional periodic variability on timescales of $t_{\rm bin}$ and $\approx 0.5 t_{\rm bin}$, as well as periodic variation of broad line widths and offsets relative to the narrow lines, which are consistent with the observations. Future observations will be able to test these predictions and hence the binary+CBD hypothesis for PG 1302.

Constraining the properties of Population III (Pop III) stars will be very challenging because they reside in small galaxies at high redshift which will be difficult to directly detect. In this paper, we suggest that intensity mapping may be a promising method to study Pop III stars. Intensity mapping is a technique proposed to measure large-scale fluctuations of galaxy line emission in three dimensions without resolving individual sources. This technique is well suited for observing many faint galaxies because it can measure their cumulative emission even if they cannot be directly detected. We focus on intensity mapping of He II recombination lines. These lines are much stronger in Pop III stars than Pop II stars because the harder spectra of Pop III stars are expected to produce many He II ionizing photons. Measuring the He II 1640 \AA{} intensity mapping signal, along with the signals from other lines such as Ly$\alpha$, H$\alpha$, and metal lines, could give constraints on the initial mass function (IMF) and star formation rate density (SFRD) of Pop III stars as a function of redshift. To demonstrate the feasibility of these observations, we estimate the strength of the Pop III He II 1640 \AA{} intensity mapping signal from $z=10-20$. We show that at $z\approx10$, the signal could be measured accurately by two different hypothetical future instruments, one which cross-correlates He II 1640 \AA{} with CO(1-0) line emission from galaxies and the other with 21 cm emission from the intergalactic medium (IGM).

Weak gravitational lensing is a powerful cosmological probe, with non--Gaussian features potentially containing the majority of the information. We examine constraints on the parameter triplet $(\Omega_m,w,\sigma_8)$ from non-Gaussian features of the weak lensing convergence field, including a set of moments (up to $4^{\rm th}$ order) and Minkowski functionals, using publicly available data from the 154deg$^2$ CFHTLenS survey. We utilize a suite of ray--tracing N-body simulations spanning 91 points in $(\Omega_m,w,\sigma_8)$ parameter space, replicating the galaxy sky positions, redshifts and shape noise in the CFHTLenS catalogs. We then build an emulator that interpolates the simulated descriptors as a function of $(\Omega_m,w,\sigma_8)$, and use it to compute the likelihood function and parameter constraints. We employ a principal component analysis to reduce dimensionality and to help stabilize the constraints with respect to the number of bins used to construct each statistic. Using the full set of statistics, we find $\Sigma_8\equiv\sigma_8(\Omega_m/0.27)^{0.55}=0.75\pm0.04$ (68% C.L.), in agreement with previous values. We find that constraints on the $(\Omega_m,\sigma_8)$ doublet from the Minkowski functionals suffer a strong bias. However, high-order moments break the $(\Omega_m,\sigma_8)$ degeneracy and provide a tight constraint on these parameters with no apparent bias. The main contribution comes from quartic moments of derivatives.

Graham et al. (2015) discovered a supermassive black hole binary (SMBHB) candidate and identified the detected 5.2 yr period of the optical variability as the orbital period of the binary. Hydrodynamical simulations predict multiple periodic components for the variability of SMBHBs, thus raising the possibility that the true period of the binary is different from 5.2 yr. We analyse the periodogram of PG1302 and find no compelling evidence for additional peaks. We derive upper limits on any additional periodic modulations in the available data, by modeling the light-curve as the sum of red noise and the known 5.2 yr periodic component, and injecting additional sinusoidal signals. We find that, with the current data, we would be able to detect with high significance (false alarm probability <1%) secondary periodic terms with periods in the range predicted by the simulations, if the amplitude of the variability was at least ~0.07 mag (compared to 0.14 mag for the main peak). A three-year follow-up monitoring campaign with weekly observations can increase the sensitivity for detecting secondary peaks ~50%, and would allow a more robust test of predictions from hydrodynamical simulations.

The first generation of stars produces a background of Lyman-Werner (LW) radiation which can photo-dissociate molecular hydrogen, increasing the mass of dark matter halos required to host star formation. Previous studies have determined the critical mass required for efficient molecular cooling with a constant LW background. However, the true background is expected to increase rapidly at early times. Neglecting this evolution could underestimate star formation in small halos that may have started to cool in the past when the LW intensity was much lower. Background evolution is a large source of uncertainty in pre-reionization predictions of the cosmological 21cm signal, which can be observed with future radio telescopes. To address this, we perform zero-dimentional one-zone calculations that follow the density, chemical abundances, and temperature of gas in the central regions of dark matter halos, including hierarchical growth and an evolving LW background. We begin by studying the physics of halos subjected to a background that increases exponentially with redshift. We find that when the intensity increases more slowly than $J_{\rm LW}(z) \propto 10^{-z/5}$, cooling in the past is a relatively small effect. We then self-consistently compute the cosmological LW background over $z=15-50$ and find that cooling in the past due to an evolving background has a modest impact. Finally, we compare these results to three-dimensional hydrodynamical cosmological simulations with varying LW histories. While only a small number of halos were simulated, the results are consistent with our one-zone calculations.

High-redshift quasar observations imply that supermassive black holes (SMBHs) larger than $\sim 10^9 ~ M_\odot$ formed before $z=6$. That such large SMBHs formed so early in the Universe remains an open theoretical problem. One possibility is that gas in atomic cooling halos exposed to strong Lyman-Werner (LW) radiation forms $10^4-10^6 ~ M_\odot$ supermassive stars which quickly collapse into black holes. We propose a scenario for direct collapse black hole (DCBH) formation based on synchronized pairs of pristine atomic cooling halos. We consider halos at very small separation with one halo being a subhalo of the other. The first halo to surpass the atomic cooling threshold forms stars. Soon after these stars are formed, the other halo reaches the cooling threshold and due to its small distance from the newly formed galaxy, is exposed to the critical LW intensity required to form a DCBH. The main advantage of this scenario is that synchronization can potentially prevent photoevaporation and metal pollution in DCBH-forming halos. Since the halos reach the atomic cooling threshold at nearly the same time, the DCBH-forming halo is only exposed to ionizing radiation for a brief period. Tight synchronization could allow the DCBH to form before stars in the nearby galaxy reach the end of their lives and generate supernovae winds. We use N-body simulations to estimate the abundance of DCBHs formed in this way. The largest source of uncertainty in our estimate is the initial mass function (IMF) of metal free stars formed in atomic cooling halos. We find that even for tight synchronization, the density of DCBHs formed in this scenario could explain the SMBHs implied by $z=6$ quasar observations. Metal pollution and photoevaporation could potentially reduce the abundance of DCBHs below that required to explain the observations in other models that rely on a high LW flux.

Residual errors in shear measurements, after corrections for instrument systematics and atmospheric effects, can impact cosmological parameters derived from weak lensing observations. Here we combine convergence maps from our suite of ray-tracing simulations with random realizations of spurious shear. This allows us to quantify the errors and biases of the triplet $(\Omega_m,w,\sigma_8)$ derived from the power spectrum (PS), as well as from three different sets of non-Gaussian statistics of the lensing convergence field: Minkowski functionals (MF), low--order moments (LM), and peak counts (PK). Our main results are: (i) We find an order of magnitude smaller biases from the PS than in previous work. (ii) The PS and LM yield biases much smaller than the morphological statistics (MF, PK). (iii) For strictly Gaussian spurious shear with integrated amplitude as low as its current estimate of $\sigma^2_{sys}\approx 10^{-7}$, biases from the PS and LM would be unimportant even for a survey with the statistical power of LSST. However, we find that for surveys larger than $\approx 100$ deg$^2$, non-Gaussianity in the noise (not included in our analysis) will likely be important and must be quantified to assess the biases. (iv) The morphological statistics (MF,PK) introduce important biases even for Gaussian noise, which must be corrected in large surveys. The biases are in different directions in $(\Omega_m,w,\sigma_8)$ parameter space, allowing self-calibration by combining multiple statistics. Our results warrant follow-up studies with more extensive lensing simulations and more accurate spurious shear estimates.

We propose a new spectral signature for supermassive black hole binaries (SMBHBs) with circumbinary gas disks: a sharp drop in flux blueward of the Lyman limit. A prominent edge is produced if the gas dominating the emission in the Lyman continuum region of the spectrum is sufficiently cold (T < 20,000 K) to contain significant neutral hydrogen. Circumbinary disks may be in this regime if the binary torques open a central cavity in the disk and clear most of the hot gas from the inner region, and if any residual UV emission from the individual BHs is either dim or intermittent. We model the vertical structure and spectra of circumbinary disks using the radiative transfer code TLUSTY, and identify the range of BH masses and binary separations producing a Lyman edge. We find that compact supermassive binaries with orbital periods of ~0.1 - 10 yr, whose gravitational waves (GWs) are expected to be detectable by pulsar timing arrays (PTAs), could have prominent Lyman edges. Such strong spectral edge features are not typically present in AGN spectra and could serve as corroborating evidence for the presence of a SMBHB.

In order to elucidate the origin of spin in both dark matter and baryons in galaxies, we have performed hydrodynamical simulations from cosmological initial conditions. We study atomic cooling haloes in the redshift range $100 > z > 9$ with masses of order $10^9{\rm M_{\odot}}$ at redshift $z=10$. We assume that the gas has primordial composition and that ${\rm H_2}$-cooling and prior star-formation in the haloes have been suppressed. We present a comprehensive analysis of the gas and dark matter properties of four halos with very low ($\lambda \approx 0.01$), low ($\lambda \approx 0.04$), high ($\lambda \approx 0.06$) and very high ($\lambda \approx 0.1$) spin parameter. Our main conclusion is that the spin orientation and magnitude is initially well described by tidal torque linear theory, but later on is determined by the merging and accretion history of each halo. We provide evidence that the topology of the merging region, i.e. the number of colliding filaments, gives an accurate prediction for the spin of dark matter and gas: halos at the center of knots will have low spin while those in the center of filaments will have high spin. The spin of a halo is given by $\lambda \approx 0.05 \times \left(\frac{7.6}{\rm number\,\,\, of \,\,\, filaments}\right)^{5.1}$.

Supermassive stars (SMSs; >10^5 Msun) formed in the first protogalaxies with virial temperature T_vir>10^4 K are expected to collapse into seeds of supermassive black hole (SMBHs) in the high-redshift universe (z>7). Fragmentation of the primordial gas is, however, a possible obstacle to SMS formation. We discuss the expected properties of a compact, metal-free, marginally unstable nuclear protogalactic disk, and the fate of the clumps formed in the disk by gravitational instability. Interior to a characteristic radius R_f=few*10^{-2} pc, the disk fragments into massive clumps with M_c~30 Msun. The clumps grow via accretion and migrate inward rapidly on a timescale of ~10^4 yr, which is comparable or shorter than the Kelvin-Helmholz time >10^4 yr. Some clumps may evolve to zero-age main sequence stars and halt gas accretion by radiative feedback, but most of the clumps can migrate inward and merge with the central protostar before forming massive stars. Moreover, we found that dust-induced-fragmentation in metal-enriched gas does not modify these conclusions unless Z> 3*10^{-4} Zsun, because clump migration below this metallicity remains as rapid as in the primordial case. Our results suggest that fragmentation of a compact, metal--poor disk can not prevent the formation of a SMS.

We present the results of 2D, moving mesh, viscous hydrodynamical simulations of accretion onto merging supermassive black hole (SMBH) binaries. We include viscous heating, shock heating, and radiative cooling, and simulate the transition from the "pre-decoupling" epoch, where the inspiral timescale is longer than the viscous timescale, to the "post-decoupling" epoch, where the inspiral timescale is shorter than the viscous timescale. We find that there is no abrupt halt to the accretion at decoupling, but rather the accretion shows a slow decay, with significant accretion well after the expected decoupling. Moreover, we find that the luminosity in X-rays is significantly higher prior to the merger, as orbital energy from the SMBH binary is converted to heat via strong shocks inside the cavity, and radiated away. Following the merger, the cavity refills viscously and the accretion rate relaxes to the Shakura-Sunyaev value, while the X-ray luminosity drops as the shocks quickly dissipate.

We report the serendipitous discoveries of companion galaxies to two high-redshift quasars. SDSS J025617.7+001904 is a z=4.79 quasar included in our recent survey of faint quasars in the SDSS Stripe 82 region. The initial MMT slit spectroscopy shows excess Lyman alpha emission extending well beyond the quasar's light profile. Further imaging and spectroscopy with LBT/MODS1 confirms the presence of a bright galaxy (i_AB = 23.6) located 2arcsec (12 kpc projected) from the quasar with strong Lyman alpha emission (EW_0 ~ 100Ang) at the redshift of the quasar, as well as faint continuum. The second quasar, CFHQS J005006.6+344522 (z=6.25), is included in our recent HST SNAP survey of z~6 quasars searching for evidence of gravitational lensing. Deep imaging with ACS and WFC3 confirms an optical dropout ~4.5 mag fainter than the quasar (Y_AB=25) at a separation of 0.9 arcsec. The red i_775-Y_105 color of the galaxy and its proximity to the quasar (5 kpc projected if at the quasar redshift) strongly favor an association with the quasar. Although it is much fainter than the quasar it is remarkably bright when compared to field galaxies at this redshift, while showing no evidence for lensing. Both systems may represent late-stage mergers of two massive galaxies, with the observed light for one dominated by powerful ongoing star formation and for the other by rapid black hole growth. Observations of close companions are rare; if major mergers are primarily responsible for high-redshift quasar fueling then the phase when progenitor galaxies can be observed as bright companions is relatively short.

Quadrupole oscillation modes in stars can resonate with incident gravitational waves (GWs), and grow non-linear at the expense of GW energy. Stars near massive black hole binaries (MBHB) can act as GW-charged batteries, cooling radiatively. Mass-loss from these stars can prompt MBHB accretion at near-Eddington rates. GW opacity is independent of amplitude, so distant resonating stars can eclipse GW sources. Absorption by the Sun of GWs from Galactic white dwarf binaries may be detectable with second-generation space-based GW detectors as a shadow within a complex diffraction pattern.

Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk.

Explaining the existence of supermassive black holes (SMBHs) larger than $\sim 10^9 M_\odot$ at redshifts $z >\sim 6$ remains an open theoretical question. One possibility is that gas collapsing rapidly in pristine atomic cooling halos ($T_{\rm vir} >\sim 10^4 \rm{K}$) produces $10^4-10^6 M_\odot$ black holes. Previous studies have shown that the formation of such a black hole requires a strong UV background to prevent molecular hydrogen cooling and gas fragmentation. Recently it has been proposed that a high UV background may not be required for halos that accrete material extremely rapidly or for halos where gas cooling is delayed due to a high baryon-dark matter streaming velocity. In this work, we point out that building up a halo with $T_{\rm vir} >\sim 10^4 \rm{K}$ before molecular cooling becomes efficient is not sufficient for forming a direct collapse black hole (DCBH). Though molecular hydrogen formation may be delayed, it will eventually form at high densities leading to efficient cooling and fragmentation. The only obvious way that molecular cooling could be avoided in the absence of strong UV radiation, is for gas to reach high enough density to cause collisional dissociation of molecular hydrogen ($\sim 10^4 ~ {\rm cm}^{-3}$) before cooling occurs. However, we argue that the minimum core entropy, set by the entropy of the intergalactic medium (IGM) when it decouples from the CMB, prevents this from occurring for realistic halo masses. This is confirmed by hydrodynamical cosmological simulations without radiative cooling. We explain the maximum density versus halo mass in these simulations with simple entropy arguments. The low densities found suggest that DCBH formation indeed requires a strong UV background.

We present the results of a calculation of the thermal spectrum from a 2D, moving mesh, high-accuracy, viscous hydrodynamical simulation of an accreting supermassive black hole binary. We include viscous heating, shock heating, and radiative cooling, evolving for longer than a viscous time so that we reach a quasi-steady accretion state. In agreement with previous work, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We also find that emission from the shock-heated mini-disks and accretion streams prevents any deficit in high-energy emission that may be expected inside the circumbinary cavity, and instead leads to a characteristic brightening of the spectrum beginning in soft X-rays.

The rapid decline in the number of strong Lyman Alpha (Lya) emitting galaxies at z > 6 provides evidence for neutral hydrogen in the IGM, but is difficult to explain with plausible models for reionization. We demonstrate that the observed reduction in Lya flux from galaxies at z > 6 can be explained by evolution in the escape fraction of ionizing photons, f_esc. We find that the median observed drop in the fraction of galaxies showing strong Lya emission, as well as the observed evolution of the Lya luminosity function both follow from a small increase in f_esc of Delta f_esc ~ 0.1 from f_esc ~ 0.6 at z ~ 6. This high escape fraction may be at odds with current constraints on the ionising photon escape fraction, which favor smaller values of f_esc < 20%. However, models that invoke a redshift evolution of f_ esc that is consistent with these constraints can suppress the z~7 Lya flux to the observed level, if they also include a small evolution in global neutral fraction of Delta x_HI ~ 0.2. Thus, an evolving escape fraction of ionising photons can be a plausible part of the explanation for evolution in the Lya emission of high redshift galaxies. More generally, our analysis also shows that the drop in the Lya fraction is quantitatively consistent with the observed evolution in the Lya luminosity functions of Lya Emitters.

The presence of quasars at redshifts z > 6 indicates the existence of supermassive black holes (SMBHs) as massive as a few times 10^9 Msun, challenging models for SMBH formation. One pathway is through the direct collapse of gas in T_{vir} > 10^4 K halos; however, this requires the suppression of H_2 cooling to prevent fragmentation. In this paper, we examine a proposed new mechanism for this suppression which relies on cold-mode accretion flows leading to shocks at high densities (n > 10^4 cm^{-3}) and temperatures (T > 10^4 K). In such gas, H_2 is efficiently collisionally dissociated. We use high-resolution numerical simulations to test this idea, demonstrating that such halos typically have lower temperature progenitors, in which cooling is efficient. Those halos do show filamentary flows; however, the gas shocks at or near the virial radius (at low densities), thus preventing the proposed collisional mechanism from operating. We do find that, if we artificially suppress H_2 formation with a high UV background, so as to allow gas in the halo center to enter the high-temperature, high-density "zone of no return", it will remain there even if the UV flux is turned off, collapsing to high density at high temperature. Due to computational limitations, we simulated only three halos. However, we demonstrate, using Monte Carlo calculations of 10^6 halo merger histories, that a few rare halos could assemble rapidly enough to avoid efficient H_2 cooling in all of their progenitor halos, provided that the UV background exceeds J_{21} ~ few at redshifts as high as z ~ 20.

We compare the efficiency of moments and Minkowski functionals (MFs) in constraining the subset of cosmological parameters (Omega_m,w,sigma_8) using simulated weak lensing convergence maps. We study an analytic perturbative expansion of the MFs in terms of the moments of the convergence field and of its spatial derivatives. We show that this perturbation series breaks down on smoothing scales below 5', while it shows a good degree of convergence on larger scales (15'). Most of the cosmological distinguishing power is lost when the maps are smoothed on these larger scales. We also show that, on scales comparable to 1', where the perturbation series does not converge, cosmological constraints obtained from the MFs are approximately 1.5-2 times better than the ones obtained from the first few moments of the convergence distribution --- provided that the latter include spatial information, either from moments of gradients, or by combining multiple smoothing scales. Including either a set of these moments or the MFs can significantly tighten constraints on cosmological parameters, compared to the conventional method of using the power spectrum alone.

The weak lensing power spectrum is a powerful tool to probe cosmological parameters. Additionally, lensing peak counts contain cosmological information beyond the power spectrum. Both of these statistics can be affected by the preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions. If not accounted for, these biases introduce systematic errors for cosmological measurements. Here we quantify these systematic errors, using convergence maps from a suite of ray-tracing N-body simulations. At the cut-off magnitude m of on-going and planned major weak lensing surveys, the logarithmic slope of the cumulative number counts s = dlog[n(>m)]/dlog(m) is in the range 0.1 < s < 0.5. At s = 0.2, expected in the I band for LSST, the inferred values of Omega_m, w and sigma_8 are biased by many sigma (where sigma denotes the marginalized error) and therefore the biases will need to be carefully modeled. We also find that the parameters are biased differently in the (Omega_m, w, sigma_8) parameter space when the power spectrum and when the peak counts are used. In particular, w derived from the power spectrum is less affected than w derived from peak counts, while the opposite is true for the best-constrained combination of [sigma_8 Omega_m^gamma] (with gamma=0.62 from the power spectrum and gamma = 0.48 from peak counts). This suggests that the combination of the power spectrum and peak counts can help mitigate the impact of magnification and size biases.

A near-equal mass binary black hole can clear a central cavity in a circumbinary accretion disk; however, previous works have revealed accretion streams entering this cavity. Here we use 2D hydrodynamical simulations to study the accretion streams and their periodic behavior. In particular, we perform a suite of simulations, covering different binary mass ratios $q=M_2/M_1$ in the range $0.01 \leq q \leq 1$. In each case, we follow the system for several thousand binary orbits, until it relaxes to a stable accretion pattern. We find the following results: (i) while the binary is efficient in maintaining a low-density cavity, the time-averaged mass accretion rate into the cavity, through narrow coherent accretion streams, is suppressed by at most a factor of $\sim 5$, compared to a disk with a single BH with the same mass; (ii) the largest suppression occurs for $q\approx 0.05$; binaries whose mass ratios are either lower or higher both suppress accretion less significantly; (iii) for $q \gtrsim 0.05$, the accretion rate is strongly modulated by the binary, and depending on the precise value of $q$, the power spectrum of the accretion rate shows either one, two, or three distinct periods; and (v) for $q \lesssim 0.05$, the accretion rate becomes steady, with no time-variations. Most binaries produced in galactic mergers are expected to have $q\gtrsim 0.05$. If the luminosity of these binaries tracks their accretion rate, then a periodogram of their light-curve could help in their identification, and to constrain their mass ratio and disk properties.

We present the results of 2D hydrodynamical simulations of circumbinary disk accretion using the finite-volume code DISCO. This code solves the 2D viscous Navier-Stokes equations on a high-resolution moving mesh which shears with the fluid flow, greatly reducing advection errors in comparison with a fixed grid. We perform a series of simulations for binary mass ratios in the range 0.026 < q < 1.0, each lasting longer than a viscous time so that we reach a quasi-steady accretion state. In each case, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We find that for q > 0.1, the binary excites eccentricity in the inner region of the circumbinary disk, creating an overdense lump which gives rise to enhanced periodicity in the accretion rate. The dependence of the periodicity on mass ratio may provide a method for observationally inferring mass ratios from measurements of the accretion rate. We also find that for all mass ratios studied, the magnitude of the accretion onto the secondary is sufficient to drive the binary toward larger mass ratio. This suggests a mechanism for biasing mass ratio distributions toward equal mass.

Pulsar timing arrays (PTAs) may eventually be able to detect not only the stochastic gravitational-wave (GW) background of SMBH binaries, but also individual, particularly massive binaries whose signals stick out above the background. In this contribution, we discuss the possibility of identifying and studying such `resolved' binaries through their electromagnetic emission. The host galaxies of such binaries are themselves expected to be also very massive and rare, so that out to redshifts z~2 a unique massive galaxy may be identified as the host. At higher redshifts, the PTA error boxes are larger and may contain as many as several hundred massive-galaxy interlopers. In this case, the true counterpart may be identified, if it is accreting gas efficiently, as an active galactic nucleus (AGN) with a peculiar spectrum and variable emission features. Specifically, the binary's tidal torques expel the gas from the inner part of the accretion disk, making it unusually dim in X-ray and UV bands and in broad optical emission lines. The tails of the broad wings of any FeK-alpha emission line may also be `clipped' and missing. The binary's orbital motion, as well as the gas motions it induces, may trigger quasiperiodic variations. These include coherent flux variability, such as luminous, multi-wavelength flares, as well as Doppler shifts of broad emission lines and `see-saw' oscillations in the FeK-alpha line. Additional features, such as evidence for a recent major merger or dual collimated jets, could also corroborate the counterpart. These properties would make resolved PTA sources stand out among AGN with similar overall luminosities and allow their identification.

Observations of quasars at redshifts z > 6 reveal that 10^9 Msol supermassive black holes (SMBHs) had already formed when the Universe was < 0.9 Gyr old. One hypothesis for the origins of these SMBHs is that they grew from the remnants of the first generation of massive stars, which formed in low-mass (~ 10^5 to 10^6 Msol) dark matter minihaloes at z > 20. This is the regime where baryonic streaming motions--the relative velocities of baryons with respect to dark matter in the early Universe--most strongly inhibit star formation by suppressing gas infall and cooling. We investigate the impact of this effect on the growth of the first SMBHs using a suite of high-fidelity, ellipsoidal-collapse Monte Carlo merger-tree simulations. We find that the suppression of seed BH formation by the streaming motions significantly reduces the number density of the most massive BHs at z > 15, but the residual effect at lower redshifts is essentially negligible. The streaming motions can reduce by a factor of few the number density of the most luminous quasars at z ~ 10-11, where such objects could be detected by the James Webb Space Telescope. We conclude, with minor theoretical caveats, that baryonic streaming motions are unlikely to pose a significant additional obstacle to the formation of the observed high-redshift quasar SMBHs. Nor do they appreciably affect the heating and reionization histories of the Universe or the merger rates of nuclear BHs in the mass and redshift ranges of interest for proposed gravitational-wave detectors.

We propose a novel use of high-redshift galaxies, discovered in deep Hubble Space Telescope (HST) fields around strong lensing clusters. These fields probe small comoving volumes (about 1000 cubic Mpc) at high magnification ({\mu} > 10), and can detect otherwise inaccessible ultra-faint galaxies. Even a few galaxies found in such small volumes require a very high number density of collapsed dark matter (DM) halos. This implies significant primordial power on small scales, allowing these observations to rule out popular alternatives to standard cold dark matter (CDM) models, such as warm dark matter (WDM). In this work, we analytically compute WDM halo mass functions at z = 10, including the effects of both particle free-streaming and residual velocity dispersion. The preliminary number density corresponding to the two galaxies at z about 10 already detected by the Cluster Lensing And Supernova survey with Hubble (CLASH) constrains the WDM particle mass to m_x > 1 (0.9) keV at 84% (99.9%) confidence limit (for a thermal relic relativistic at decoupling). This limit depends only on the WDM halo mass function and, unlike previous constraints on m_x, is independent of any astrophysical modeling. The forthcoming HST Frontier Fields can significantly tighten these constraints.

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 observations of a luminous galaxy at redshift z=6.573 --- the end of the reioinization epoch --- which has been spectroscopically confirmed twice. The first spectroscopic confirmation comes from slitless HST ACS grism spectra from the PEARS survey (Probing Evolution And Reionization Spectroscopically), which show a dramatic continuum break in the spectrum at restframe 1216 A wavelength. The second confirmation is done with Keck + DEIMOS. The continuum is not clearly detected with ground-based spectra, but high wavelength resolution enables the Lyman alpha emission line profile to be determined. We compare the line profile to composite line profiles at redshift z=4.5. The Lyman alpha line profile shows no signature of a damping wing attenuation, confirming that the intergalactic gas is ionized at redshift z=6.57. Spectra of Lyman breaks at yet higher redshifts will be possible using comparably deep observations with IR-sensitive grisms, even at redshifts where Lyman alpha is too attenuated by the neutral IGM to be detectable using traditional spectroscopy from the ground.

Galactic nuclei are expected to be densely populated with stellar and intermediate mass black holes. Exploring this population will have important consequences for the observation prospects of gravitational waves as well as understanding galactic evolution. The gas cloud G2 currently approaching Sgr A* provides an unprecedented opportunity to probe the black hole and neutron star population of the Galactic nucleus. We examine the possibility of a G2-black hole encounter and its detectability with current X-ray satellites, such as Chandra and NuSTAR. We find that multiple encounters are likely to occur close to the pericenter, which may be detectable upon favorable circumstances. This opportunity provides an additional, important science case for leading X-ray observatories to closely follow G2 on its way to the nucleus.

Structures in Warm Dark Matter (WDM) models are exponentially suppressed below a certain scale, characterized by the dark matter particle mass, $m_{\rm x}$. Since structures form hierarchically, the presence of collapsed objects at high-redshifts can set strong lower limits on $m_{\rm x}$. We place robust constraints on $m_{\rm x}$ using recent results from the {\it Swift} database of high-redshift gamma-ray bursts (GRBs). We parameterize the redshift evolution of the ratio between the cosmic GRB rate and star formation rate (SFR) as $\propto (1+z)^\alpha$, thereby allowing astrophysical uncertainties to partially mimic the cosmological suppression of structures in WDM models. Using a maximum likelihood estimator on two different $z>4$ GRB subsamples (including two bursts at $z>8$), we constrain $m_{\rm x} \gtrsim 1.6$-1.8 keV at 95% C.L., when marginalized over a flat prior in $\alpha$. We further estimate that 5 years of a SVOM-like mission would tighten these constraints to $m_{\rm x} \gtrsim 2.3 $ keV. Our results show that GRBs are a powerful probe of high-redshift structures, providing robust and competitive constraints on $m_{\rm x}$.

(abridged) When a sufficiently massive satellite (or secondary) black hole is embedded in a gas disk around a (primary) supermassive black hole, it can open an empty gap in the disk. A gap-opening secondary close to the primary will leave an imprint in the broad component of the FeKa emission line, which varies in a unique and predictable manner. If the gap persists into the innermost disk, the effect consists of a pair of dips in the broad line which ripple blue-ward and red-ward from the line centroid energy respectively, as the gap moves closer to the primary. This ripple effect could be unambiguously detectable and allow an electromagnetic monitoring of massive black hole mergers as they occur. As the mass ratio of the secondary to primary black hole increases to q>0.01, we expect the gap to widen, possibly clearing a central cavity in the inner disk, which shows up in the broad FeKa line component. If the secondary stalls at >100r_{g} in its in-migration, due to low co-rotating gas mass, a detectable ripple effect occurs in the broad line component on the disk viscous timescale as the inner disk drains and the outer disk is dammed. If the secondary maintains an accretion disk within a central cavity, due to dam bursting or leakage, a periodic 'see-saw' oscillation effect is exhibited in the observed line profile. Here we demonstrate the range of ripple effect signatures potentially detectable with Astro-H and IXO/Athena, and oscillation effects potentially detectable with XMM or LOFT for a wide variety of merger and disk conditions. Observations of the ripple effect and periodic oscillations can be used to provide an early warning of gravitational radiation emission from the AGN.

We have performed cosmo-hydro simulations using the RAMSES code to study atomic cooling (ACHs) haloes at z=10 with masses 5E7Msun<~M<~2E9Msun. We assume primordial gas and H2-cooling and prior star-formation have been suppressed. We analysed 19 haloes (gas and DM) at a resolution of ~10 (proper) pc, selected from a total volume of ~2E3 (comoving) Mpc3. This is the largest statistical hydro-sim. study of ACHs at z>10 to date. We examine the morphology, angular momentum (AM), thermodynamic, and turbulence of these haloes, in order to assess the prevalence of disks and supermassive black holes (SMBHs). We find no correlation between either the magnitude or the direction of the AM of the gas and its parent DM halo. Only 3 haloes form rotationally supported cores. Two of the most massive haloes form massive, compact overdense blobs. These blobs have an accretion rate ~0.5 Msun/yr (at a distance of 100 pc), and are possible sites of SMBH formation. Our results suggest that the degree of rotational support and the fate of the gas in a halo is determined by its large-scale environment and merger history. In particular, the two haloes forming blobs are located at knots of the cosmic web, cooled early on, and experienced many mergers. The gas in these haloes is lumpy and highly turbulent, with Mach N. >~ 5. In contrast, the haloes forming rotationally supported cores are relatively more isolated, located midway along filaments, cooled more recently, and underwent fewer mergers. Thus, the gas in these haloes is less lumpy and less turbulent (Mach <~ 4), and could retain most of its AM. The remaining 14 haloes have intermediate properties. If verified in a larger sample of haloes and with additional physics, our results will have implications for observations of the highest-redshift galaxies and quasars with JWST.

The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.

The spectra of several high-redshift (z>6) quasars have shown evidence for a Gunn-Peterson (GP) damping wing, indicating a substantial mean neutral hydrogen fraction (x_HI > 0.03) in the z ~ 6 intergalactic medium (IGM). However, previous analyses assumed that the IGM was uniformly ionized outside of the quasar's HII region. Here we relax this assumption and model patchy reionization scenarios for a range of IGM and quasar parameters. We quantify the impact of these differences on the inferred x_HI, by fitting the spectra of three quasars: SDSS J1148+5251 (z=6.419), J1030+0524 (z=6.308), and J1623+3112 (z=6.247). We find that the best-fit values of x_HI in the patchy models agree well with the uniform case. More importantly, we confirm that the observed spectra favor the presence of a GP damping wing, with peak likelihoods decreasing by factors of > few - 10 when the spectra are modeled without a damping wing. We also find that the Ly alpha absorption spectra, by themselves, cannot distinguish the damping wing in a relatively neutral IGM from a damping wing in a highly ionized IGM, caused either by an isolated neutral patch, or by a damped Ly alpha absorber (DLA). However, neutral patches in a highly ionized universe (x_HI < 0.01), and DLAs with the large required column densities (N_HI > few x 10^{20} cm^{-2}) are both rare. As a result, when we include reasonable prior probabilities for the line of sight (LOS) to intercept either a neutral patch or a DLA at the required distance of ~ 40-60 comoving Mpc away from the quasar, we find strong lower limits on the neutral fraction in the IGM, x_HI > 0.1 (at 95% confidence). This strengthens earlier claims that a substantial global fraction of hydrogen in the z~6 IGM is in neutral form.

Peaks in two-dimensional weak lensing (WL) maps contain significant cosmological information, complementary to the WL power spectrum. This has recently been demonstrated using N-body simulations which neglect baryonic effects. Here we employ ray-tracing N-body simulations in which we manually steepen the density profile of each dark matter halo, mimicking the cooling and concentration of baryons into dark matter potential wells. We find, in agreement with previous works, that this causes a significant increase in the amplitude of the WL power spectrum on small scales (spherical harmonic index l>1,000). We then study the impact of the halo concentration increase on the peak counts, and find the following. (i) Low peaks (with convergence 0.02 < kappa_peak < 0.08), remain nearly unaffected. These peaks are created by a constellation of several halos with low masses (10^12-10^13 M_sun) and large angular offsets from the peak center (> 0.5 R_vir); as a result, they are insensitive to the central halo density profiles. These peaks contain most of the cosmological information, and thus provide an unusually sensitive and unbiased probe. (ii) The number of high peaks (with convergence kappa_peak > 0.08) is increased. However, when the baryon effects are neglected in cosmological parameter estimation, then the high peaks lead to a modest bias, comparable to that from the power spectrum on relatively large-scales (l<2000), and much smaller than the bias from the power spectrum on smaller scales (l>2,000). (iii) In the 3D parameter space (sigma_8, Omega_m, w), the biases from the high peaks and the power spectra are in different directions. This suggests the possibility of "self-calibration": the combination of peak counts and power spectrum can simultaneously constrain baryonic physics and cosmological parameters.

We study the interaction of a supermassive black hole (SMBH) binary and a standard radiatively efficient thin accretion disk. We examine steady-state configurations of the disk and migrating SMBH system, self-consistently accounting for tidal and viscous torques and heating, radiative diffusion limited cooling, gas and radiation pressure, and the decay of the binary's orbit. We obtain a "phase diagram" of the system as a function of binary parameters, showing regimes in which both the disk structure and migration have a different character. Although massive binaries can create a central gap in the disk at large radii, the tidal barrier of the secondary causes a significant pile-up of gas outside of its orbit, which can lead to the closing of the gap. We find that this spillover occurs at an orbital separation as large as ~200 M_7^{-1/2} gravitational radii, where M = 10^7 M_7 Msun is the total binary mass. If the secondary is less massive than ~10^6 Msun, then the gap is closed before gravitational waves (GWs) start dominating the orbital decay. In this regime, the disk is still strongly perturbed, but the piled-up gas continuously overflows as in a porous dam, and crosses inside the secondary's orbit. The corresponding migration rate, which we label Type 1.5, is slower than the usual limiting cases known as Type I and II migration. Compared to an unperturbed disk, the steady-state disk in the overflowing regime is up to several hundred times brighter in the optical bands. Surveys such as PanSTARRS or LSST may discover the periodic variability of this population of binaries. Our results imply that the circumbinary disks around SMBHs can extend to small radii during the last stages of their merger, when they are detectable by LISA, and may produce coincident electromagnetic (EM) emission similar to active galactic nuclei (AGN).

Many astrophysical binaries, from planets to black holes, exert strong torques on their circumbinary accretion disks, and are expected to significantly modify the disk structure. Despite the several decade long history of the subject, the joint evolution of the binary + disk system has not been modeled with self-consistent assumptions for arbitrary mass ratios and accretion rates. Here we solve the coupled binary-disk evolution equations analytically in the strongly perturbed limit, treating the azimuthally-averaged angular momentum exchange between the disk and the binary and the modifications to the density, scale-height, and viscosity self-consistently, including viscous and tidal heating, diffusion limited cooling, radiation pressure, and the orbital decay of the binary. We find a solution with a central cavity and a migration rate similar to those previously obtained for Type-II migration, applicable for large masses and binary separations, and near-equal mass ratios. However, we identify a distinct new regime, applicable at smaller separations and masses, and mass ratio in the range 0.001< q < 0.1. For these systems, gas piles up outside the binary's orbit, but rather than creating a cavity, it continuously overflows as in a porous dam. The disk profile is intermediate between a weakly perturbed disk (producing Type-I migration) and a disk with a gap (with Type-II migration). However, the migration rate of the secondary is typically slower than both Type-I and Type-II rates. We term this new regime "Type-1.5" migration.

Observations of high-redshift quasars at z>6 imply that supermassive black holes (SMBHs) with masses over a billion solar masses were in place less than 1 Gyr after the Big Bang. If these SMBHs assembled from "seed" BHs left behind by the first stars, then they must have accreted gas at close to the Eddington limit during a large fraction (>50%) of the time. A generic problem with this scenario, however, is that the mass density in million-solar-mass SMBHs at z=6 already exceeds the locally observed SMBH mass density by several orders of magnitude; in order to avoid this overproduction, BH seed formation and growth must become significantly less efficient in less massive protogalaxies, while proceeding uninterrupted in the most massive galaxies that formed first. Using Monte-Carlo realizations of the merger and growth history of BHs, we show that X-rays from the earliest accreting BHs can provide such a feedback mechanism. Our calculations paint a self-consistent picture of black-hole-made climate change, in which the first miniquasars---among them the ancestors of the z>6 quasar SMBHs---globally warm the IGM and suppress the formation and growth of subsequent generations of BHs. We present two specific models with global miniquasar feedback that provide excellent agreement with recent estimates of the z=6 SMBH mass function. For each of these models, we estimate the rate of BH mergers at z>6 that could be detected by the proposed gravitational-wave observatory eLISA/NGO.

It is commonly believed that the earliest stages of star-formation in the Universe were self-regulated by global radiation backgrounds - either by the ultraviolet Lyman-Werner (LW) photons emitted by the first stars (directly photodissociating H_2), or by the X-rays produced by accretion onto the black hole (BH) remnants of these stars (heating the gas but catalyzing H_2 formation). Recent studies have suggested that a significant fraction of the first stars may have had low masses (a few M_sun). Such stars do not leave BH remnants and they have softer spectra, with copious infrared (IR) radiation at photon energies around 1eV. Similar to LW and X-ray photons, these photons have a mean-free path comparable to the Hubble distance, building up an early IR background. Here we show that if soft-spectrum stars, with masses of a few M_sun, contributed more than 1% of the UV background (or their mass fraction exceeded 90%), then their IR radiation dominated radiative feedback in the early Universe. The feedback is different from the UV feedback from high-mass stars, and occurs through the photo-detachment of H^- ions, necessary for efficient H_2 formation. Nevertheless, we find that the baryon fraction which must be incorporated into low-mass stars in order to suppress H_2-cooling is only a factor of few higher than for high-mass stars.

Supermassive black holes (SMBHs) are common in local galactic nuclei, and SMBHs as massive as several billion solar masses already exist at redshift z=6. These earliest SMBHs may grow by the combination of radiation-pressure-limited accretion and mergers of stellar-mass seed BHs, left behind by the first generation of metal-free stars, or may be formed by more rapid direct collapse of gas in rare special environments where dense gas can accumulate without first fragmenting into stars. This chapter offers a review of these two competing scenarios, as well as some more exotic alternative ideas. It also briefly discusses how the different models may be distinguished in the future by observations with JWST, (e)LISA and other instruments.

Pulsar timing arrays (PTAs) are expected to detect gravitational waves (GWs) from individual low-redshift (z<1.5) compact supermassive (M>10^9 Msun) black hole (SMBH) binaries with orbital periods of approx. 0.1 - 10 yrs. Identifying the electromagnetic (EM) counterparts of these sources would provide confirmation of putative direct detections of GWs, present a rare opportunity to study the environments of compact SMBH binaries, and could enable the use of these sources as standard sirens for cosmology. Here we consider the feasibility of such an EM identification. We show that because the host galaxies of resolved PTA sources are expected to be exceptionally massive and rare, it should be possible to find unique hosts of resolved sources out to redshift z=0.2. At higher redshifts, the PTA error boxes are larger, and may contain as many as 100 massive-galaxy interlopers. The number of candidates, however, remains tractable for follow-up searches in upcoming wide-field EM surveys. We develop a toy model to characterize the dynamics and the thermal emission from a geometrically thin, gaseous disc accreting onto a PTA-source SMBH binary. Our model predicts that at optical and infrared frequencies, the source should appear similar to a typical luminous active galactic nucleus (AGN). However, owing to the evacuation of the accretion flow by the binary's tidal torques, the source might have an unusually low soft X-ray luminosity and weak UV and broad optical emission lines, as compared to an AGN powered by a single SMBH with the same total mass. For sources near z=1, the decrement in the rest-frame UV should be observable as an extremely red optical color. These properties would make the PTA sources stand out among optically luminous AGN, and could allow their unique identification.

Recent studies have shown that the number counts of convergence peaks N(kappa) in weak lensing (WL) maps, expected from large forthcoming surveys, can be a useful probe of cosmology. We follow up on this finding, and use a suite of WL convergence maps, obtained from ray-tracing N-body simulations, to study (i) the physical origin of WL peaks with different heights, and (ii) whether the peaks contain information beyond the convergence power spectrum P_ell. In agreement with earlier work, we find that high peaks (with amplitudes >~ 3.5 sigma, where sigma is the r.m.s. of the convergence kappa) are typically dominated by a single massive halo. In contrast, medium-height peaks (~0.5-1.5 sigma) cannot be attributed to a single collapsed dark matter halo, and are instead created by the projection of multiple (typically, 4-8) halos along the line of sight, and by random galaxy shape noise. Nevertheless, these peaks dominate the sensitivity to the cosmological parameters w, sigma_8, and Omega_m. We find that the peak height distribution and its dependence on cosmology differ significantly from predictions in a Gaussian random field. We directly compute the marginalized errors on w, sigma_8, and Omega_m from the N(kappa) + P_ell combination, including redshift tomography with source galaxies at z_s=1 and z_s=2. We find that the N(kappa) + P_ell combination has approximately twice the cosmological sensitivity compared to P_ell alone. These results demonstrate that N(kappa) contains non-Gaussian information complementary to the power spectrum.

In this paper, we show that Minkowski Functionals (MFs) of weak gravitational lensing (WL) convergence maps contain significant non-Gaussian, cosmology-dependent information. To do this, we use a large suite of cosmological ray-tracing N-body simulations to create mock WL convergence maps, and study the cosmological information content of MFs derived from these maps. Our suite consists of 80 independent 512^3 N-body runs, covering seven different cosmologies, varying three cosmological parameters Omega_m, w, and sigma_8 one at a time, around a fiducial LambdaCDM model. In each cosmology, we use ray-tracing to create a thousand pseudo-independent 12 deg^2 convergence maps, and use these in a Monte Carlo procedure to estimate the joint confidence contours on the above three parameters. We include redshift tomography at three different source redshifts z_s=1, 1.5, 2, explore five different smoothing scales theta_G=1, 2, 3, 5, 10 arcmin, and explicitly compare and combine the MFs with the WL power spectrum. We find that the MFs capture a substantial amount of information from non-Gaussian features of convergence maps, i.e. beyond the power spectrum. The MFs are particularly well suited to break degeneracies and to constrain the dark energy equation of state parameter w (by a factor of ~ three better than from the power spectrum alone). The non-Gaussian information derives partly from the one-point function of the convergence (through V_0, the "area" MF), and partly through non-linear spatial information (through combining different smoothing scales for V_0, and through V_1 and V_2, the boundary length and genus MFs, respectively). In contrast to the power spectrum, the best constraints from the MFs are obtained only when multiple smoothing scales are combined.

The power spectrum of cosmic infrared background (CIB) anisotropies is sensitive to the connection between star formation and dark matter halos over the entire cosmic star formation history. Here we develop a model that associates star-forming galaxies with dark matter halos and their subhalos. The model is based on a parameterized relation between the dust-processed infrared luminosity and (sub)halo mass. By adjusting 3 free parameters, we attempt to simultaneously fit the 4 frequency bands of the Planck measurement of the CIB anisotropy power spectrum. To fit the data, we find that the star-formation efficiency must peak on a halo mass scale of ~ 5x10^12 solar mass and the infrared luminosity per unit mass must increase rapidly with redshift. By comparing our predictions with a well-calibrated phenomenological model for shot noise, and with a direct observation of source counts, we show that the mean duty cycle of the underlying infrared sources must be near unity, indicating that the CIB is dominated by long-lived quiescent star formation, rather than intermittent short "star bursts". Despite the improved flexibility of our model, the best simultaneous fit to all four Planck channels remains relatively poor. We discuss possible further extensions to alleviate the remaining tension with the data. Our model presents a theoretical framework for a future joint analysis of both background anisotropy and source count measurements.

The ability of primordial gas to cool in proto-galactic haloes exposed to Lyman-Werner (LW) radiation is critically dependent on the self-shielding of H_2. We perform radiative transfer calculations of LW line photons, post-processing outputs from three-dimensional adaptive mesh refinement (AMR) simulations of haloes with T_vir > 10^4 K at redshifts around z=10. We calculate the optically thick photodissociation rate numerically, including the effects of density, temperature, and velocity gradients in the gas, as well as line overlap and shielding of H_2 by HI, over a large number of sight-lines. In low-density regions (n<10^4 cm^-3) the dissociation rates exceed those obtained using most previous approximations by more than an order of magnitude; the correction is smaller at higher densities. We trace the origin of the deviations primarily to inaccuracies of (i) the most common fitting formula (Draine & Bertoldi 1996) for the suppression of the dissociation rate and (ii) estimates for the effective shielding column density from local properties of the gas. The combined effects of gas temperature and velocity gradients are comparatively less important, typically altering the spherically averaged rate only by a factor of less than two. We present a simple modification to the DB96 fitting formula for the optically thick rate which improves agreement with our numerical results to within approx. 15 per cent, and can be adopted in future simulations. We find that estimates for the effective shielding column can be improved by using the local Sobolev length. Our correction to the H_2 self-shielding reduces the critical LW flux to suppress H_2-cooling in T_vir>10^4 K haloes by an order of magnitude; this increases the number of such haloes in which supermassive (approx. M=10^5 M_sun) black holes may have formed.

Independent lines of evidence suggest that the first stars, which ended the cosmic dark ages, came in pairs, rather than singly. This could change the prevailing view that the early Universe had a Swiss-cheese-like appearance.

We study the effects of a thin gaseous accretion disk on the inspiral of a stellar--mass black hole into a supermassive black hole. We construct a phenomenological angular momentum transport equation that reproduces known disk effects. Disk torques modify the gravitational wave phase evolution to detectable levels with LISA for reasonable disk parameters. The Fourier transform of disk-modified waveforms acquires a correction with a different frequency trend than post-Newtonian vacuum terms. Such inspirals could be used to detect accretion disks with LISA and to probe their physical parameters.

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.

Clusters of galaxies have been used extensively to determine cosmological parameters. A major difficulty in making best use of Sunyaev-Zel'dovich (SZ) and X-ray observations of clusters for cosmology is that using X-ray observations it is difficult to measure the temperature distribution and therefore determine the density distribution in individual clusters of galaxies out to the virial radius. Observations with the new generation of SZ instruments are a promising alternative approach. We use clusters of galaxies drawn from high-resolution adaptive mesh refinement (AMR) cosmological simulations to study how well we should be able to constrain the large-scale distribution of the intra-cluster gas (ICG) in individual massive relaxed clusters using AMiBA in its configuration with 13 1.2-m diameter dishes (AMiBA13) along with X-ray observations. We show that non-isothermal beta models provide a good description of the ICG in our simulated relaxed clusters. We use simulated X-ray observations to estimate the quality of constraints on the distribution of gas density, and simulated SZ visibilities (AMiBA13 observations) for constraints on the large-scale temperature distribution of the ICG. We find that AMiBA13 visibilities should constrain the scale radius of the temperature distribution to about 50% accuracy. We conclude that the upgraded AMiBA, AMiBA13, should be a powerful instrument to constrain the large-scale distribution of the ICG.

It has been shown that HD molecules can form efficiently in metal-free gas collapsing into massive protogalactic halos at high redshift. The resulting radiative cooling by HD can lower the gas temperature to that of the cosmic microwave background, T_CMB=2.7(1+z)K, significantly below the temperature of a few 100 K achievable via H_2-cooling alone, and thus reduce the masses of the first generation of stars. Here we consider the suppression of HD-cooling by UV irradiation in the Lyman-Werner (LW) bands. We include photo-dissociation of both H_2 and HD, and explicitly compute the self-shielding and shielding of both molecules by neutral hydrogen as well as the shielding of HD by H_2. We use a simplified dynamical collapse model, and follow the chemical and thermal evolution of the gas, in the presence of a UV background. We find that a LW flux of J_crit = 1e-22 erg/cm^2/sr/s/Hz is able to suppress HD cooling and thus prevent collapsing primordial gas from reaching temperatures below 100 K. The main reason for the lack of HD cooling for J>J_crit is the partial photo-dissociation of H_2, which prevents the gas from reaching sufficiently low temperatures (T<150K) for HD to become the dominant coolant; direct HD photo-dissociation is unimportant except for a narrow range of fluxes and column densities. Since the prevention of HD-cooling requires only partial H_2 photo-dissociation, the critical flux J_crit is modest, and is below the UV background required to reionize the universe at redshift z=10-20. We conclude that HD-cooling can reduce the masses of typical stars only in rare halos forming well before the epoch of reionization.

Supermassive black holes (SMBHs) are common in local galactic nuclei, and SMBHs as massive as several billion solar masses already exist at redshift z=6. These earliest SMBHs may arise by the combination of Eddington-limited growth and mergers of stellar-mass seed BHs left behind by the first generation of metal-free stars, or by the rapid direct collapse of gas in rare special environments where the gas can avoid fragmenting into stars. In this contribution, I review these two competing scenarios. I also briefly mention some more exotic ideas and how the different models may be distinguished in the future by LISA and other instruments.

The coalescence of a supermassive black hole binary (SMBHB) is thought to be accompanied by an electromagnetic (EM) afterglow, produced by the viscous infall of the surrounding circumbinary gas disk after the merger. It has been proposed that once the merger has been detected in gravitational waves (GWs) by LISA, follow-up EM searches for this afterglow can help identify the EM counterpart of the LISA source. Here we study whether the afterglows may be sufficiently bright and numerous to be detectable in EM surveys alone. The viscous afterglow, which lasts for years to decades for SMBHBs in LISA's sensitivity window, is characterized by rapid increases in both the bolometric luminosity and in the spectral hardness of the source. If quasar activity is triggered by the same major galaxy mergers that produce SMBHBs, then the afterglow could be interpreted as a signature of the birth of a quasar. Using an idealized model for the post-merger viscous spreading of the circumbinary disk and the resulting light curve, and using the observed luminosity function of quasars as a proxy for the SMBHB merger rate, we delineate the survey requirements for identifying such birthing quasars. If circumbinary disks have a high disk surface density and viscosity, an all-sky soft X-ray survey with a sensitivity of ~<3x10^-14 erg s^-1 cm^-2 and a time resolution of ~months could identify dozens of birthing quasars with sustained brightening rates of >10%/yr. If >1% of the X-ray emission is reprocessed into optical frequencies, birthing quasars could also be identified in optical transient surveys such as the LSST. Distinguishing a birthing quasar from other variable sources may be facilitated by the monotonic hardening of its spectrum, but will likely remain challenging. This reinforces the notion that joint EM-plus-GW observations offer the best prospects for identifying the EM signatures of SMBHB mergers.

It has been proposed that primordial gas in early dark matter halos, with virial temperatures above 10^4 K, can avoid fragmentation and undergo rapid collapse, possibly resulting in a supermassive black hole (SMBH). This requires the gas to avoid cooling and to remain at temperatures near T=10^4 K. We show that this condition can be satisfied in the presence of a sufficiently strong primordial magnetic field, which heats the collapsing gas via ambipolar diffusion. If the field has a strength above B = 3.6 (comoving) nG, the collapsing gas is kept warm (T=10^4K) until it reaches the critical density n_crit=10^3 cm^{-3} at which the roto-vibrational states of H_2 approach local thermodynamic equilibrium. H_2-cooling then remains inefficient, and the gas temperature stays near 10^4K, even as it continues to collapse to higher densities. The critical magnetic field strength required to permanently suppress H_2-cooling is somewhat higher than upper limit of approx. 2 nG from the cosmic microwave background (CMB). However, it can be realized in the rare (2-3)-sigma regions of the spatially fluctuating B-field; these regions contain a sufficient number of halos to account for the z=6 quasar BHs.

The gravitational waves (GWs) emitted by inspiraling binary black holes, expected to be detected by the Laser Interferometer Space Antenna (LISA), could be used to determine the luminosity distance to these sources with the unprecedented precision of <~ 1%. We study cosmological parameter constraints from such standard sirens, in the presence of gravitational lensing by large-scale structure. Lensing introduces magnification with a probability distribution function (PDF) whose shape is highly skewed and depends on cosmological parameters. We use Monte-Carlo simulations to generate mock samples of standard sirens, including a small intrinsic scatter, as well as the additional, larger scatter from lensing, in their inferred distances. We derive constraints on cosmological parameters, by simultaneously fitting the mean and the distribution of the residuals on the distance vs redshift (d_L - z) Hubble diagram. We find that for standard sirens at redshift z ~ 1, the sensitivity to a single cosmological parameter, such as the matter density Omega_m, or the dark energy equation of state w, is ~ 50%-80% tighter when the skewed lensing PDF is used, compared to the sensitivity derived from a Gaussian PDF with the same variance. When these two parameters are constrained simultaneously, the skewness yields a further enhanced improvement (by ~ 120%), owing to the correlation between the parameters. The sensitivity to the amplitude of the matter power spectrum, sigma_8 from the cosmological dependence of the PDF alone, however, is ~ 20% worse than that from the Gaussian PDF. At higher redshifts, the PDF resembles a Gaussian more closely, and the effects of the skewness become less prominent. These results highlight the importance of obtaining an accurate and reliable PDF of the lensing convergence, in order to realize the full potential of standard sirens as cosmological probes.

Finding electromagnetic (EM) counterparts of future gravitational wave (GW) sources would bring rich scientific benefits. A promising possibility, in the case of the coalescence of a super-massive black hole binary (SMBHB), is that prompt emission from merger-induced disturbances in a supersonic circumbinary disk may be detectable. We follow the post-merger evolution of a thin, zero-viscosity circumbinary gas disk with two-dimensional simulations, using the hydrodynamic code FLASH. We analyze perturbations arising from the 530 km/s recoil of a 10^6 M_sun binary, oriented in the plane of the disk, assuming either an adiabatic or a pseudo-isothermal equation of state for the gas. We find that a single-armed spiral shock wave forms and propagates outward, sweeping up about 20% of the mass of the disk. The morphology and evolution of the perturbations agrees well with those of caustics predicted to occur in a collisionless disk. Assuming that the disk radiates nearly instantaneously to maintain a constant temperature, we estimate the amount of dissipation and corresponding post-merger light-curve. The luminosity rises steadily on the time-scale of months, and reaches few times 10^{43} erg/s, corresponding to about 10% of the Eddington luminosity of the central SMBHB. We also analyze the case in which gravitational wave emission results in a 5% mass loss in the merger remnant. The mass-loss reduces the shock overdensities and the overall luminosity of the disk by 15-20%, without any other major effects on the spiral shock pattern.

UV radiation from early astrophysical sources could have a large impact on subsequent star formation in nearby protogalaxies. Here we study the radiative feedback from the first, short-lived stars using hydrodynamical simulations with transient UV backgrounds (UVBs) and persistent Lyman-Werner backgrounds (LWBs) of varying intensity. We extend our prior work in Mesinger et al. (2006), by studying a more typical region whose proto-galaxies form at lower redshifts, z~13-20, in the epoch likely preceding the bulk of reionization. We confirm our previous results that feedback in the relic HII regions resulting from such transient radiation, is itself transient. Feedback effects dwindle away after ~30% of the Hubble time, and the same critical specific intensity of J_UV~0.1 x 10^{-21} ergs/s/cm^2/Hz/sr separates positive and negative feedback regimes. Additionally, we discover a second episode of eventual positive feedback in halos which have not yet collapsed when their progenitor regions were exposed to the transient UVB. This eventual positive feedback appears in all runs, regardless of the strength of the UVB. However, this feedback regime is very sensitive to the presence of Lyman-Werner radiation, and notable effects disappear under fairly modest background intensities of J_LW>10^{-3} x 10^{-21} ergs/s/cm^2/Hz/sr. We conclude that UV radiative feedback in relic HII regions, although a complicated process, seems unlikely to have a major impact on the progress of cosmological reionization, provided that present estimates of the lifetime and luminosity of a PopIII star are accurate. More likely is that the build-up of the LWB ultimately governs the feedback strength until a persistent UV background can be established. [abridged]

Arguably the best hope of understanding the tail end of the reionization of the intergalactic medium (IGM) at redshift z > 6 is through the detection and characterization of the Gunn-Peterson (GP) damping wing absorption of the IGM in bright quasar spectra. However, the use of quasar spectra to measure the IGM damping wing requires a model of the quasar's intrinsic Lyman-alpha emission line. Here we quantify the uncertainties in the intrinsic line shapes, and how those uncertainties affect the determination of the IGM neutral fraction. We have assembled a catalog of high-resolution HST spectra of the emission lines of unobscured low-redshift quasars, and have characterized the variance in the shapes of their lines. We then add simulated absorption from the high-redshift IGM to these quasar spectra in order to determine the corresponding uncertainties in reionization constraints using current and future samples of z > 6 quasar spectra. We find that, if the redshift of the Lyman-alpha emission line is presumed to coincide with the systemic redshift determined from metal lines, the inferred IGM neutral fraction is systematically biased to low values due to a systematic blueshift of the Lyman-alpha line relative to the metal lines. If a similar blueshift persists in quasars at z > 6, this bias strengthens previous claims of a significant neutral hydrogen fraction at z ~ 6. This technique is capable of making a robust distinction between a highly ionized (x_IGM ~ 10^-3) and a neutral (x_IGM = 1) IGM with even a few bright quasars.

Scaling relations among galaxy cluster observables, which will become available in large future samples of galaxy clusters, could be used to constrain not only cluster structure, but also cosmology. We study the utility of this approach, employing a physically motivated parametric model to describe cluster structure, and applying it to the expected relation between the Sunyaev-Zel'dovich decrement (S_\nu) and the emission-weighted X-ray temperature (T_ew). The slope and normalization of the entropy profile, the concentration of the dark matter potential, the pressure at the virial radius, and the level of non-thermal pressure support, as well as the mass and redshift-dependence of these quantities are described by free parameters. With a suitable choice of fiducial parameter values, the cluster model satisfies several existing observational constraints. We employ a Fisher matrix approach to estimate the joint errors on cosmological and cluster structure parameters from a measurement of S_\nu vs. T_ew in a future survey. We find that different cosmological parameters affect the scaling relation differently: predominantly through the baryon fraction (\Omega_m and \Omega_b), the virial overdensity (w_0 and w_a for low-z clusters) and the angular diameter distance (w_0, w_a for high-z clusters; \Omega_DE and h). We find that the cosmology constraints from the scaling relation are comparable to those expected from the number counts (dN/dz) of the same clusters. The scaling relation approach is relatively insensitive to selection effects and it offers a valuable consistency check; combining the information from the scaling relation and dN/dz is also useful to break parameter degeneracies and help disentangle cluster physics from cosmology.

We propose counting peaks in weak lensing (WL) maps, as a function of their height, to probe models of dark energy and to constrain cosmological parameters. Because peaks can be identified in two-dimensional WL maps directly, they can provide constraints that are free from potential selection effects and biases involved in identifying and determining the masses of galaxy clusters. We have run cosmological N-body simulations to produce WL convergence maps in three models with different constant values of the dark energy equation of state parameter, w=-0.8, -1, and -1.2, with a fixed normalization of the primordial power spectrum (corresponding to present-day normalizations of sigma8=0.742, 0.798, and 0.839, respectively). By comparing the number of WL peaks in 8 convergence bins in the range of -0.1 < kappa < 0.2, in multiple realizations of a single simulated 3x3 degree field, we show that the first (last) pair of models can be distinguished at the 95% (85%) confidence level. A survey with depth and area (20,000 sq. degrees), comparable to those expected from LSST, should have a factor of approx. 50 better parameter sensitivity. We find that relatively low-amplitude peaks (kappa = 0.03), which typically do not correspond to a single collapsed halo along the line of sight, account for most of this sensitivity. We study a range of smoothing scales and source galaxy redshifts (z_s). With a fixed source galaxy density of 15/arcmin^2, the best results are provided by the smallest scale we can reliably simulate, 1 arcminute, and z_s=2 provides substantially better sensitivity than z_s< 1.5.

In the absence of H_2 molecules, the primordial gas in early dark matter halos with virial temperatures just above T_vir >~ 10^4 K cools by collisional excitation of atomic H. Although it cools efficiently, this gas remains relatively hot, at a temperature near T ~ 8000 K, and consequently might be able to avoid fragmentation and collapse directly into a supermassive black hole (SMBH). In order for H_2--formation and cooling to be strongly suppressed, the gas must be irradiated by a sufficiently intense ultraviolet (UV) flux. We performed a suite of three--dimensional hydrodynamical adaptive mesh refinement (AMR) simulations of gas collapse in three different protogalactic halos with T_vir >~ 10^4 K, irradiated by a UV flux with various intensities and spectra. We determined the critical specific intensity, Jcrit, required to suppress H_2 cooling in each of the three halos. For a hard spectrum representative of metal--free stars, we find (in units of 10^{-21} erg s^{-1} Hz^{-1} sr^{-1} cm^{-2}) 10^4<Jcrit<10^5, while for a softer spectrum, which is characteristic of a normal stellar population, and for which H^{-} --dissociation is important, we find 30<Jcrit<300. These values are a factor of 3--10 lower than previous estimates. We attribute the difference to the higher, more accurate H_2 collisional dissociation rate we adopted. The reduction in Jcrit exponentially increases the number of rare halos exposed to super--critical radiation. When H_2 cooling is suppressed, gas collapse starts with a delay, but it ultimately proceeds more rapidly. The infall velocity is near the increased sound speed, and an object as massive as M ~ 10^5 solar mass may form at the center of these halos, compared to the M ~ 10^2 solar mass stars forming when H_2--cooling is efficient.

The supermassive black holes (SMBHs) massive enough to power the bright redshift ~6 quasars observed in the Sloan Digital Sky Survey (SDSS) are thought to have assembled by mergers and/or gas accretion from less massive "seed" BHs. If the seeds are the ~100 Msol remnant BHs from the first generation of stars, they must be in place well before redshift z=6, and must avoid being ejected from their parent proto-galaxies by the large (few hundred km/s) kicks they suffer from gravitational-radiation induced recoil during mergers with other BHs. We simulate the SMBH mass function at redshift z=6 using dark matter (DM) halo merger trees, coupled with a prescription for the halo occupation fraction, accretion histories, and radial recoil trajectories of the growing BHs. Our purpose is (i) to map out plausible scenarios for successful assembly of the z~6 quasar BHs by exploring a wide region of parameter space, and (ii) to predict the rate of low-frequency gravitational wave events detectable by the Laser Interferometer Space Antenna (LISA) for each such scenario. Our main findings are as follows: (1) ~100 Msol seed BHs can grow into the SDSS quasar BHs without super--Eddington accretion, but only if they form in minihalos by z~30 and subsequently accrete ~60% of the time; (2) the scenarios with optimistic assumptions required to explain the SDSS quasar BHs overproduce the mass density in lower--mass (10^5 to 10^7 Msol) SMBHs by a factor of 100-1000, unless seeds stop forming, or accrete at a severely diminished rates or duty cycles (e.g. due to feedback), at z < 20-30. We also present several successful assembly models and their LISA detection rates, including a "maximal" model that gives the highest rate (~30/yr at z~6) without overproducing the total SMBH density.

Supermassive black hole binaries (SMBHBs) in galactic nuclei are thought to be a common by-product of major galaxy mergers. We use simple disk models for the circumbinary gas and for the binary-disk interaction to follow the orbital decay of SMBHBs with a range of total masses (M) and mass ratios (q), through physically distinct regions of the disk, until gravitational waves (GWs) take over their evolution. Prior to the GW-driven phase, the viscous decay is in the stalled "secondary-dominated" regime. SMBHBs spend a non-negligible fraction of 10^7 years at orbital periods t_var between a day and a year. A dedicated optical or X-ray survey could identify coalescing SMBHBs statistically, as a population of periodically variable quasars, whose abundance N_var is proportional to t_var^alpha, in a range of periods t_var around tens of weeks. SMBHBs with M < 10^7 M_sun, with 0.5 < alpha < 1.5, would probe the physics of viscous orbital decay, whereas the detection of a population of higher-mass binaries, with alpha=8/3, would confirm that their decay is driven by GWs. The lowest mass SMBHBs (M < 10^{5-6} M_sun) enter the GW-driven regime at short orbital periods, in the frequency band of the Laser Interferometric Space Antenna (LISA). While viscous processes are strongly sub-dominant in the last few years of coalescence, they could reduce the amplitude of any unresolved background of near-stationary LISA sources. We discuss constraints on the SMBHB population available from existing data, and the sensitivity and sky coverage requirements for a detection in future surveys. SMBHBs may also be identified from velocity shifts in their spectra; we discuss the expected abundance of SMBHBs as a function of their orbital velocity.

Supermassive black hole binaries (SMBHBs) produced in galaxy mergers are thought to complete their coalescence, below separations of r_GW=10^{-3} (M_BH/10^8 M_sun)^{3/4} pc, as their orbit decays due to the emission of gravitational waves (GWs). It may be possible to identify such GW-driven inspirals statistically in an electromagnetic (EM) survey for variable sources. A GW-driven binary spends a characteristic time T_GW at each orbital separation r_orb < r_GW that scales with the corresponding orbital time t_orb as T_GW = (const) t_orb^{8/3}. If the coalescing binary produces variations in the EM emission on this timescale, then it could be identified as a variable source with a characteristic period t_var = t_orb. The incidence rate of sources with similar inferred BH masses, showing near-periodic variability on the time-scale t_var, would then be proportional to t_var^{8/3}. Luminosity variations corresponding to a fraction f_Edd<0.01 of the Eddington luminosity would have been missed in current surveys. However, if the binary inspirals are associated with quasars, we show that a dedicated survey could detect the population of SMBHBs with a range of periods around tens of weeks. The discovery of a population of periodic sources whose abundance obeys N_var = (const) t_var^{8/3} would confirm (i) that the orbital decay is indeed driven by GWs, and (ii) that circumbinary gas is present at small orbital radii and is being perturbed by the BHs. Deviations from the t_var^{8/3} power-law could constrain the structure of the circumbinary gas disk and viscosity-driven orbital decay. We discuss constraints from existing data, and quantify the sensitivity and sky coverage that could yield a detection in future surveys.

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.

Cold dark matter (CDM) hierarchical structure formation models predict the existence of large-scale accretion shocks between the virial and turnaround radii of clusters of galaxies. Kocsis et al. (2005) suggest that the Sunyaev-Zel'dovich (SZ) signal associated with such shocks might be observable with the next generation radio interferometer, ALMA. We study the three--dimensional distribution of accretion shocks around individual clusters of galaxies drawn from adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of LCDM (dark energy dominated CDM) models. In relaxed clusters, we find two distinct sets of shocks. One set ("virial shocks"), with Mach numbers of 2.5-4, is located at radii 0.9-1.3 Rvir, where Rvir is the spherical infall estimate of the virial radius, covering about 40-50% of the total surface area around clusters at these radii. Another set of stronger shocks ("external shocks") is located farther out, at about 3 Rvir, with large Mach numbers (~100), covering about 40-60% of the surface area. We simulate SZ surface brightness maps of relaxed massive galaxy clusters drawn from high resolution AMR runs, and conclude that ALMA should be capable of detecting the virial shocks in massive clusters of galaxies. More simulations are needed to improve estimates of astrophysical noise and to determine optimal observational strategies.

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?

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.

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

We conduct a Markov Chain Monte Carlo study of the Dvali-Gabadadze-Porrati (DGP) self-accelerating braneworld scenario given the cosmic microwave background (CMB) anisotropy, supernovae and Hubble constant data by implementing an effective dark energy prescription for modified gravity into a standard Einstein-Boltzmann code. We find no way to alleviate the tension between distance measures and horizon scale growth in this model. Growth alterations due to perturbations propagating into the bulk appear as excess CMB anisotropy at the lowest multipoles. In a flat cosmology, the maximum likelihood DGP model is nominally a 5.3 sigma poorer fit than Lambda CDM. Curvature can reduce the tension between distance measures but only at the expense of exacerbating the problem with growth leading to a 4.8 sigma result that is dominated by the low multipole CMB temperature spectrum. While changing the initial conditions to reduce large scale power can flatten the temperature spectrum, this also suppresses the large angle polarization spectrum in violation of recent results from WMAP5. The failure of this model highlights the power of combining growth and distance measures in cosmology as a test of gravity on the largest scales.

Supermassive black hole binaries (SMBHBs) with masses in the range 10^4-10^7 M_sun/(1+z), produced in galaxy mergers, are thought to complete their coalescence due to the emission of gravitational waves (GWs). The anticipated detection of the GWs by the LISA will constitute a milestone for fundamental physics and astrophysics. While the GW signatures themselves will provide a treasure trove of information, if the source can be securely identified in electromagnetic (EM) bands, this would open up entirely new scientific opportunities, to probe fundamental physics, astrophysics, and cosmology. We discuss several ideas, involving wide-field telescopes, that may be useful in locating electromagnetic counterparts to SMBHBs detected by LISA. In particular, the binary may produce a variable electromagnetic flux, such as a roughly periodic signal due to the orbital motion prior to coalescence, or a prompt transient signal caused by shocks in the circumbinary disk when the SMBHB recoils and "shakes" the disk. We discuss whether these time-variable EM signatures may be detectable, and how they can help in identifying a unique counterpart within the localization errors provided by LISA. We also discuss a possibility of identifying a population of coalescing SMBHBs statistically, in a deep optical survey for periodically variable sources, before LISA detects the GWs directly. The discovery of such sources would confirm that gas is present in the vicinity and is being perturbed by the SMBHB - serving as a proof of concept for eventually finding actual LISA counterparts.

We report on the detection by Swift of GRB 080913, and subsequent optical/near-infrared follow-up observations by GROND which led to the discovery of its optical/NIR afterglow and the recognition of its high-z nature via the detection of a spectral break between the i' and z' bands. Spectroscopy obtained at the ESO-VLT revealed a continuum extending down to lambda = 9400 A, and zero flux for 7500 A < lambda<9400 A, which we interpret as the onset of a Gunn-Peterson trough at z=6.695+-0.025 (95.5% conf. level), making GRB 080913 the highest redshift GRB to date, and more distant than the highest-redshift QSO. We note that many redshift indicators which are based on promptly available burst or afterglow properties have failed for GRB 080913. We report on our follow-up campaign and compare the properties of GRB 080913 with bursts at lower redshift. In particular, since the afterglow of this burst is fainter than typical for GRBs, we show that 2m-class telescopes can identify most high-redshift GRBs.

The earliest generation of stars and black holes must have established an early 'Lyman-Werner' background (LWB) at high redshift, prior to the epoch of reionization. Because of the long mean free path of photons with energies E<13.6 eV, the LWB was nearly uniform. However, some variation in the LWB is expected due to the discrete nature of the sources, and their highly clustered spatial distribution. In this paper, we compute the probability distribution function (PDF) of the LW flux that irradiates dark matter (DM) halos collapsing at high-redshift (z~10). Our model accounts for (i) the clustering of DM halos, (ii) Poisson fluctuations in the number of corresponding star forming galaxies, and (iii) scatter in the LW luminosity produced by halos of a given mass (calibrated using local observations). We find that > 99% of the DM halos are illuminated by a LW flux within a factor of 2 of the global mean value. However, a small fraction, ~1e-8 to 1e-6, of DM halos with virial temperatures above 1e4 K have a close luminous neighbor within < 10 kpc, and are exposed to a LW flux exceeding the global mean by a factor of > 20, or to J_(21,LW)> 1e3 (in units of 1e-21 erg/s/Hz/sr/cm^2). This large LW flux can photo--dissociate H_2 molecules in the gas collapsing due to atomic cooling in these halos, and prevent its further cooling and fragmentation. Such close halo pairs therefore provide possible sites in which primordial gas clouds collapse directly into massive black holes (M_BH~ 1e4 - 1e6 M_sun), and subsequently grow into supermassive (M_BH > 1e9 M_sun) black holes by z~6.

We propose to use a simple observable, the fractional area of "hot spots" in weak lensing mass maps which are detected with high significance, to determine background cosmological parameters. Because these high-shear regions are directly related to the physical non-linear structures of the universe, they derive cosmological information mainly from the non-linear regime of density fluctuations. We show that in combination with future cosmic microwave background (CMB) anisotropy measurements, this method can place constraints on cosmological parameters that are comparable to those from the redshift distribution of galaxy cluster abundances. In particular, adding the information from this simple shear statistic derived from an LSST-like survey to CMB anisotropy information from a Planck-like experiment improves the constraints on the dark energy parameters (energy density and equation of state) and the normalization of the matter power spectrum by over a factor of ten. The main advantage of the statistic proposed in this paper is that projection effects, normally the main source of uncertainty when determining the presence and the mass of a galaxy cluster, here serve as a source of information.

The high sensitivity of JWST will open a new window on the end of the cosmological dark ages. Small stellar clusters, with a stellar mass of several 10^6 M_sun, and low-mass black holes (BHs), with a mass of several 10^5 M_sun should be directly detectable out to redshift z=10, and individual supernovae (SNe) and gamma ray burst (GRB) afterglows are bright enough to be visible beyond this redshift. Dense primordial gas, in the process of collapsing from large scales to form protogalaxies, may also be possible to image through diffuse recombination line emission, possibly even before stars or BHs are formed. In this article, I discuss the key physical processes that are expected to have determined the sizes of the first star-clusters and black holes, and the prospect of studying these objects by direct detections with JWST and with other instruments. The direct light emitted by the very first stellar clusters and intermediate-mass black holes at z>10 will likely fall below JWST's detection threshold. However, JWST could reveal a decline at the faint-end of the high-redshift luminosity function, and thereby shed light on radiative and other feedback effects that operate at these early epochs. JWST will also have the sensitivity to detect individual SNe from beyond z=10. In a dedicated survey lasting for several weeks, thousands of SNe could be detected at z>6, with a redshift distribution extending to the formation of the very first stars at z>15. Using these SNe as tracers may be the only method to map out the earliest stages of the cosmic star-formation history. Finally, we point out that studying the earliest objects at high redshift will also offer a new window on the primordial power spectrum, on 100 times smaller scales than probed by current large-scale structure data.

We propose to use a simple observable, the fractional area of "hot spots" in weak gravitational lensing mass maps which are detected with high significance, to determine background cosmological parameters. Because these high-convergence regions are directly related to the physical nonlinear structures of the universe, they derive cosmological information mainly from the nonlinear regime of density fluctuations. We show that in combination with future cosmic microwave background anisotropy measurements, this method can place constraints on cosmological parameters that are comparable to those from the redshift distribution of galaxy cluster abundances. The main advantage of the statistic proposed in this paper is that projection effects, normally the main source of uncertainty when determining the presence and the mass of a galaxy cluster, here serve as a source of information.

Primordial gas in protogalactic dark matter (DM) halos with virial temperatures Tvir > 10^4 K begins to cool and condense via atomic hydrogen. Provided this gas is irradiated by a strong ultraviolet (UV) flux and remains free of H2 and other molecules, it has been proposed that the halo with Tvir ~10^4 K may avoid fragmentation, and lead to the rapid formation of a supermassive black hole (SMBH) as massive as M=10^5-10^6 Msun. This ``head--start'' would help explain the presence of SMBHs with inferred masses of several x 10^9 Msun, powering the bright quasars discovered in the Sloan Digital Sky Survey at redshift z>~6. However, high-redshift DM halos with Tvir~10^4K are likely already enriched with at least trace amounts of metals and dust produced by prior star-formation in their progenitors. Here we study the thermal and chemical evolution of low-metallicity gas exposed to extremely strong UV radiation fields. Our results, obtained in one-zone models, suggest that gas fragmentation is inevitable above a critical metallicity, whose value is between Zcr~3x10^{-4} Zsun (in the absence of dust) and as low as Zcr~ 5 x 10^{-6} Zsun (with a dust-to-gas mass ratio of about 0.01 Z/Zsun). We propose that when the metallicity exceeds these critical values, dense clusters of low--mass stars may form at the halo nucleus. Relatively massive stars in such a cluster can then rapidly coalesce into a single more massive object, which may produce an intermediate-mass BH remnant with a mass up to M <~10^2-10^3 Msun.

We utilize the local velocity dispersion function (VDF) of spheroids, together with their inferred age--distributions, to predict the VDF at higher redshifts (0<z<6), under the assumption that (i) most of the stars in each nearby spheroid formed in a single episode, and (ii) the velocity dispersion sigma remained nearly constant afterward. We assume further that a supermassive black hole (BH) forms concurrently with the stars, and within ~1 Gyr of the formation of the potential well of the spheroid, and that the relation between the mass of the BH and host velocity dispersion maintains the form M_BH ~ sigma^{beta} with beta~4, but with the normalization allowed to evolve with redshift as ~(1+z)^{alpha}. We compute the BH mass function associated with the VDF at each redshift, and compare the accumulated total BH mass density with that inferred from the integrated quasar luminosity function (LF; the so--called Soltan argument). This comparison is insensitive to the assumed duty cycle or Eddington ratio of quasar activity, and we find that the match between the two BH mass densities favors a relatively mild redshift evolution, with alpha ~ 0.26, with a positive evolution as strong as alpha>1.3 excluded at the 99% confidence level. A direct match between the characteristic BH mass in the VDF--based and quasar LF--based BH mass functions also yields a mean Eddington ratio of lambda ~ 0.5-1 that is roughly constant within 0<z<3. A strong positive evolution in the M_BH-sigma relation is still allowed by the data if galaxies increase, on average, their velocity dispersions since the moment of formation, due to dissipative processes. If we assume that the mean velocity dispersion of the host galaxies evolves as sigma(z)=sigma(0)*(1+z)^{-gamma}, we find a lower limit of gamma>0.23 for alpha>1.5. abridged

The notion that microparsec-scale black holes can be used to probe gigaparsec-scale physics may seem counterintuitive, at first. Yet, the gravitational observatory LISA will detect cosmologically-distant coalescing pairs of massive black holes, accurately measure their luminosity distance and help identify an electromagnetic counterpart or a host galaxy. A wide variety of new black hole studies and a gravitational version of Hubble's diagram become possible if host galaxies are successfully identified. Furthermore, if dark energy is a manifestation of large-scale modified gravity, deviations from general relativistic expectations could become apparent in a gravitational signal propagated over cosmological scales, especially when compared to the electromagnetic signal from a same source. Finally, since inspirals of white dwarfs into massive black holes at cosmological distances may permit pre-merger localizations, we suggest that careful monitoring of these events and any associated electromagnetic counterpart could lead to high-precision cosmological measurements with LISA.

One of the most dramatic signatures of the reionization era may be the enormous ionized bubbles around luminous quasars (with radii reaching ~40 comoving Mpc), which may survive as "fossil'' ionized regions long after their source shuts off. Here we study how the inhomogeneous intergalactic medium (IGM) evolves inside such fossils. The average recombination rate declines rapidly with time, and the brief quasar episode significantly increases the mean free path inside the fossil bubbles. As a result, even a weak ionizing background generated by galaxies inside the fossil can maintain it in a relatively highly and uniformly ionized state. For example, galaxies that would ionize 20-30% of hydrogen in a random patch of the IGM can maintain 80-90% ionization inside the fossil, for a duration much longer than the average recombination time in the IGM. Quasar fossils at z<10 thus retain their identity for nearly a Hubble time, and will appear "gray,'' distinct from both the average IGM (which has a "swiss-cheese" ionization topology and a lower mean ionized fraction), and from fully-ionized bubbles around active quasars. More distant fossils, at z>10 have a weaker galaxy-generated ionizing background and a higher gas density, so they can attain a swiss-cheese topology similar to the rest of the IGM, but with a smaller contrast between the ionized bubbles and the partially neutral regions separating them. Analogous HeIII-fossils should exist around the epoch of HeII/HeIII reionization at z~3, although rapid recombinations inside the HeIII-fossils will be more common. Our model of inhomogeneous recombination also applies to "double reionization'' models and shows that a non-monotonic reionization history is even more unlikely than previously thought.

We develop an analytic approach to study inhomogeneous reionization on large scales by solving the equations of ionization balance and radiative transfer to first order in perturbations. Given the spatial distribution and spectrum of the ionizing sources, our formalism can be used to predict the large scale power spectra of fluctuations in the abundances of HII, HI and radiation. Our approach avoids common approximations/assumptions in existing analytic methods -- for instance, we do not assume a specific ionization topology from the outset; nor do we make a step-function bubble-like approximation to the HII distribution. Applying our formalism to sources biased according to the Press-Schechter prescription, we find: 1. reionization always proceeds ``inside--out'', with dense regions more highly ionized, at least on large scales; 2. on sufficiently large scales, HII, HI and radiation exhibits a scale independent bias relative to dark matter; 3. the bias is suppressed on scales comparable to or smaller than the mean free path of the ionizing photons; 4. if the ionizing source spectrum is sufficiently soft, the HII bias closely tracks the source bias for most of the reionization process but drops precipitously after percolation; 5. if the ionizing source spectrum is sufficiently hard, the HII bias drops in a more steady fashion throughout the reionization process. The tools developed here will be useful for interpreting future 21 cm, CMB and Lyman-alpha forest observations, both to learn about the reionization astrophysics (such as the hardness of the source spectrum and therefore the nature of the ionizing sources) and to possibly extract interesting cosmological information.

High-redshift quasars (z >~ 6) drive ionization fronts into the intergalactic medium (IGM). If the thickness of the front can be measured, it can provide a novel constraint on the ionizing spectral energy distribution (SED). Here we follow the propagation of an I-front into a uniform IGM, and compute its thickness for a range of possible quasar spectra and ages. We also explore the effects of uniform and non-uniform ionizing backgrounds. We find that even for hard spectra, the fronts are initially thin, with a thickness much smaller than the mean free path of ionizing photons, but the thickness increases as the front approaches equilibrium in 10^8 - 10^9 years, and can eventually significantly exceed simple estimates based on the mean free path. With a high intrinsic hydrogen column density obscuring the source (log(N_H/cm^-2) >~ 19.2) or a hard power-law spectrum combined with some obscuration (e.g. dlog(F_\nu)/dlog(\nu) >~ -1.2 at log(N_H/cm^-2) >~ 18.0), the thickness of the front exceeds ~1 physical Mpc and may be measurable from the morphology of its redshifted 21cm signal. We find that the highly ionized inner part of the front, which may be probed by Lyman line absorption spectra, remains sharp for bright quasars unless a large obscuring column (log(N_H/cm^-2) >~ 19.2) removes most of their ionizing photons up to ~40 eV. For obscured sources with log(N_H/cm^-2) >~ 19.8, embedded in a significantly neutral IGM, the black Lyman-alpha trough (where the neutral fraction is ~10^-3) underestimates the size of the HII region by a factor of >~4.

Supermassive black hole binaries (BHBs) produced in galaxy mergers recoil at the time of their coalescence due to the emission of gravitational waves (GWs). We simulate the response of a thin, 2D disk of collisionless particles, initially on circular orbits around a 10^6 M_sun BHB, to kicks that are either parallel or perpendicular to the initial orbital plane. Typical kick velocities (v_k) can exceed the sound speed in a circumbinary gas disk. While the inner disk is strongly bound to the recoiling binary, the outer disk is only weakly bound or unbound. This leads to differential motions in the disturbed disk that increase with radius and can become supersonic at ~700 Schwarzschild radii for v_k ~500 km/s, implying that shocks form beyond this radius. We indeed find that kicks in the disk plane lead to immediate strong density enhancements (within weeks) in a tightly wound spiral caustic, propagating outward at the speed v_k. Concentric density enhancements are also observed for kicks perpendicular to the disk, but are weaker and develop into caustics only after a long delay (>1 year). Unless both BH spins are low or precisely aligned with the orbital angular momentum, a significant fraction (> several %) of kicks are sufficiently large and well aligned with the orbital plane for strong shocks to be produced. The shocks could result in an afterglow whose characteristic photon energy increases with time, from the UV (~10eV) to the soft X-ray (~100eV) range, between one month and one year after the merger. This could help identify EM counterparts to GW sources discovered by LISA.

Electromagnetic (EM) counterparts to supermassive black hole binary mergers observed by LISA can be localized to within the field of view of astronomical instruments ~10 deg^2 hours to weeks prior to coalescence. The temporal coincidence of any prompt EM counterpart with a gravitationally-timed merger may offer the best chance of identifying a unique host galaxy. We discuss the challenges posed by searches for prompt EM counterparts and propose novel observational strategies to address them. In particular, we discuss the size and shape evolution of the LISA localization error ellipses on the sky, and quantify the requirements for dedicated EM surveys of the area prior to coalescence. A triggered EM counterpart search campaign will require monitoring a several-square degree area. It could aim for variability at the 24-27 mag level in optical bands, for example, which corresponds to 1-10% of the Eddington luminosity of the prime LISA sources of 10^6-10^7 Msun BHs at z=1-2, on time-scales of minutes to hours, the orbital time-scale of the binary in the last 2-4 weeks. A cross-correlation of the period of any variable EM signal with the quasi-periodic gravitational waveform over 10-1000 cycles may aid the detection. Alternatively, EM searches can detect a transient signal accompanying the coalescence. We highlight the measurement of differences in the arrival times of photons and gravitons from the same cosmological source as a valuable independent test of the massive character of gravity, and of possible violations of Lorentz invariance in the gravity sector.

Non-gravitational processes, such as feedback from galaxies and their active nuclei, are believed to have injected excess entropy into the intracluster gas, and therefore to have modified the density profiles in galaxy clusters during their formation. Here we study a simple model for this so-called preheating scenario, and ask (i) whether it can simultaneously explain both global X-ray scaling relations and number counts of galaxy clusters, and (ii) whether the amount of entropy required evolves with redshift. We adopt a baseline entropy profile that fits recent hydrodynamic simulations, modify the hydrostatic equilibrium condition for the gas by including approx. 20% non-thermal pressure support, and add an entropy floor K_0 that is allowed to vary with redshift. We find that the observed luminosity-temperature (L-T) relations of low-redshift (z=0.05) HIFLUGCS clusters and high-redshift (z=0.8) WARPS clusters are best simultaneously reproduced with an evolving entropy floor of K_0(z)=341(1+z)^{-0.83}h^{-1/3} keV cm^2. If we restrict our analysis to the subset of bright (kT > 3 keV) clusters, we find that the evolving entropy floor can mimic a self-similar evolution in the L-T scaling relation. This degeneracy with self-similar evolution is, however, broken when (0.5 < kT < 3 keV) clusters are also included. The approx. 60% entropy increase we find from z=0.8 to z=0.05 is roughly consistent with that expected if the heating is provided by the evolving global quasar population. Using the cosmological parameters from the WMAP 3-year data with sigma_8=0.76, our best-fit model underpredicts the number counts of the X-ray galaxy clusters compared to those derived from the 158 deg^2 ROSAT PSPC survey. Treating sigma_8 as a free parameter, we find a best-fit value of sigma_8=0.80+/- 0.02.

We apply the technique of parameter-splitting to existing cosmological data sets, to check for a generic failure of dark energy models. Given a dark energy parameter, such as the energy density Omega_Lambda or equation of state w, we split it into two meta-parameters with one controlling geometrical distances, and the other controlling the growth of structure. Observational data spanning Type Ia Supernovae, the cosmic microwave background (CMB), galaxy clustering, and weak gravitational lensing statistics are fit without requiring the two meta-parameters to be equal. This technique checks for inconsistency between different data sets, as well as for internal inconsistency within any one data set (e.g., CMB or lensing statistics) that is sensitive to both geometry and growth. We find that the cosmological constant model is consistent with current data. Theories of modified gravity generally predict a relation between growth and geometry that is different from that of general relativity. Parameter-splitting can be viewed as a crude way to parametrize the space of such theories. Our analysis of current data already appears to put sharp limits on these theories: assuming a flat universe, current data constrain the difference Omega_Lambda(geom) - Omega_Lambda(grow) to be -0.0044 +/- 0.0058 (68% C.L.); allowing the equation of state w to vary, the difference w(geom) - w(grow) is constrained to be 0.37 +/- 0.37 (68% C.L.). Interestingly, the region w(grow) > w(geom), which should be generically favored by theories that slow structure formation relative to general relativity, is quite restricted by data already. We find w(grow) < -0.80 at 2 sigma. As an example, the best-fit flat Dvali-Gabadadze-Porrati (DGP) model approximated by our parametrization lies beyond the 3 sigma contour for constraints from all the data sets.

At the epoch of reionization, when the high-redshift inter-galactic medium (IGM) is being enriched with metals, the 63.2 micron fine structure line of OI is pumped by the ~ 1300 AA soft UV background and introduces a spectral distortion in the Cosmic Microwave Background (CMB). Here we use a toy model for the spatial distribution of neutral oxygen, assuming metal bubbles surround dark matter halos, and compute the fluctuations of this distortion, and the angular power spectrum it imprints on the CMB. We discuss the dependence of the power spectrum on the velocity of the winds polluting the IGM with metals, the minimum mass of the halos producing these winds, and on the cosmic epoch when the OI pumping occurs. We find that, although the clustering signal of the CMB distortion is weak \delta y_{rms} ~ 10^{-7} (roughly corresponding to a temperature anisotropy of few nK), it may be reachable in deep integrations with high-sensitivity infrared detectors. Even without a detection, these instruments should be able to useful constraints on the heavy element enrichment history of the IGM.

Light travel time delays distort the apparent shapes of HII regions surrounding bright quasars during early stages of cosmic reionization. Individual HII regions may remain undetectable in forthcoming redshifted 21 cm experiments. However, the systematic deformation along the line of sight may be detectable statistically, either by stacking tomographic 21cm images of quasars identified, for example, by JWST, or as small-scale anisotropy in the three-dimensional 21cm power spectrum. Here we consider the detectability of this effect. The anisotropy is largest when HII regions are large and expand rapidly, and we find that if bright quasars contributed to the early stages of reionization, then they can produce significant anisotropy, on scales comparable to the typical sizes of HII regions of the bright quasars (approx. 30 Mpc and below). The effect therefore cannot be ignored when analyzing future 21cm power spectra on small scales. If 10 percent of the volume of the IGM at redshift z=10 is ionized by quasars with typical ionizing luminosity of S= 5 x 10^{56} photons/second, the distortions can enhance by more than 10 percent the 21cm power spectrum in the radial (redshift) direction, relative to the transverse directions. The level of this anisotropy exceeds that due to redshift-space distortion, and has the opposite sign. We show that on-going experiments such as MWA should be able to detect this effect. A detection would reveal the presence of bright quasars, and shed light on the ionizing yield and age of the ionizing sources, and the distribution and small-scale clumping of neutral intergalactic gas in their vicinity.

We investigate the effect of planetary rotation on the transit spectrum of an extrasolar giant planet. During ingress and egress, absorption features arising from the planet's atmosphere are Doppler shifted by of order the planet's rotational velocity (~1-2 km/s) relative to where they would be if the planet were not rotating. We focus in particular on the case of HD209458b, which ought to be at least as good a target as any other known transiting planet. For HD209458b, this shift should give rise to a small net centroid shift of ~60 cm/s on the stellar absorption lines. Using a detailed model of the transmission spectrum due to a rotating star transited by a rotating planet with an isothermal atmosphere, we simulate the effect of the planet's rotation on the shape of the spectral lines, and in particular on the magnitude of their width and centroid shift. We then use this simulation to determine the expected signal-to-noise ratio for distinguishing a rotating from a non-rotating planet, and assess how this S/N scales with various parameters of HD209458b. We find that with a 6 m telescope, an equatorial rotational velocity of ~2 km/s could be detected with a S/N~5 by accumulating the signal over many transits over the course of several years. With a 30 m telescope, the time required to make such a detection reduces to less than 2 months.

Recent observations by Ferrarese et al. (2006) and Wehner et al. (2006) reveal that a majority of galaxies contain a central massive object (CMO), either a supermassive black hole (SMBH) or a compact stellar nucleus, regardless of the galaxy mass or morphological type, and that there is a tight relation between the masses of CMOs and those of the host galaxies. Several recent studies show that feedback from black holes can successfully explain the $\msigma$ correlation in massive elliptical galaxies that contain SMBHs. However, puzzles remain in spirals or dwarf spheroids that do not appear to have black holes but instead harbor a compact central stellar cluster. Here we use three-dimensional, smoothed particle hydrodynamics simulations of isolated galaxies to study the formation and evolution of CMOs in bulgeless disk galaxies, and simulations of merging galaxies to study the transition of the CMO--host mass relation from late-type bulgeless spirals to early-type ellipticals. Our results suggest that the observed correlations may be established primarily by the depletion of gas in the central region by accretion and star-formation, and may hold for all galaxy types. A systematic search for CMOs in the nuclei of bulgeless disk galaxies would offer a test of this conclusion. (Abridged)

The continuous improvement in localization errors (sky position and distance) in real time as LISA observes the gradual inspiral of a supermassive black hole (SMBH) binary can be of great help in identifying any prompt electromagnetic counterpart associated with the merger. We develop a new method, based on a Fourier decomposition of the time-dependent, LISA-modulated gravitational-wave signal, to study this intricate problem. The method is faster than standard Monte Carlo simulations by orders of magnitude. By surveying the parameter space of potential LISA sources, we find that counterparts to SMBH binary mergers with total mass M~10^5-10^7 M_Sun and redshifts z<~3 can be localized to within the field of view of astronomical instruments (~deg^2) typically hours to weeks prior to coalescence. This will allow targeted searches for variable electromagnetic counterparts as the merger proceeds, as well as monitoring of the most energetic coalescence phase. A rich set of astrophysical and cosmological applications would emerge from the identification of electromagnetic counterparts to these gravitational-wave standard sirens.

Recent observations imply that the observed number counts of Lya Emitters (LAEs) evolved significantly between z=5.7 and z=6.5. It has been suggested that this was due to a rapid evolution in the ionisation state, and hence transmission of the IGM which caused Lya flux from z=6.5 galaxies to be more strongly suppressed. In this paper we consider the joint evolution of the Lya and UV luminosity functions (LFs) and show that the IGM transmission evolved between z=6.5 and z=5.7 by a factor 1.1 <R < 1.8 (95% CL). This result is insensitive to the underlying model of the Lya LF (as well as cosmic variance). Using a model for IGM transmission, we find that the evolution of the mean IGM density through cosmic expansion alone may result in a value for the ratio of transmissions as high as R=1.3. Thus, the existing LFs do not provide evidence for overlap. Furthermore, the constraint R<1.8 suggests that the Universe at z=6.5 was more than half ionised by volume, i.e. x_i,V>0.5.

Energetic feedback processes during the formation of galaxy clusters may have heated and ionized a large fraction of the intergalactic gas in proto-cluster regions. When such a highly ionized hot ``super-bubble'' falls along the sightline to a background quasar, it would be seen as a large void, with little or no absorption, in the Lyman alpha forest. We examine the spectra of 137 quasars in the Sloan Digital Sky Survey, to search for such voids, and find no clear evidence of their existence. The size distribution of voids in the range 5-70 Angstrom (corresponding to physical sizes of approximately 3-35 comoving Mpc/h) is consistent with the standard model for the Lyman alpha forest without additional hot bubbles. We adapt a physical model for HII bubble growth during cosmological reionization (Furlanetto, Zaldarriaga and Hernquist 2004), to describe the expected size-distribution of hot super-bubbles at redshift around z = 3. This model incorporates the conjoining of bubbles around individual neighboring galaxies. Using the non-detection of voids, we find that models in which the volume filling factor of hot bubbles exceeds approximately 20 percent at z=3 can be ruled out, primarily because they overproduce the number of large (40-50 Angstrom) voids. We conclude that any pre-heating mechanism that explains galaxy cluster observations must avoid heating the low-density gas in the proto-cluster regions, either by operating relatively recently (z<3) or by depositing entropy in the high-density regions.

We make use of hydrodynamical simulations of the intergalactic medium (IGM) to create model quasar absorption spectra. We compare these model spectra with the observed Keck spectra of three z>6.2 quasars with full Gunn-Peterson troughs: SDSS J1148+5251 (z=6.42), SDSS J1030+0524 (z=6.28), and SDSS J1623+3112 (z=6.22). We fit the probability density distributions (PDFs) of the observed Ly alpha optical depths with those generated from the simulation, by exploring a range of values for the size of the quasar's surrounding HII region, R_S, the volume-weighted mean neutral hydrogen fraction in the ambient IGM, x_H, and the quasar's ionizing photon emissivity, N_Q. In order to avoid averaging over possibly large sightline-to-sightline fluctuations in IGM properties, we analyze each observed quasar independently. We find the following results for J1148+5251, J1030+0524, and J1623+3112: The best-fit sizes R_S are 40, 41, and 29 (comoving) Mpc, respectively. These constraints are tight, with only ~ 10% uncertainties, comparable to those caused by redshift--determination errors. The best-fit values of N_Q are 2.1, 1.3, and 0.9 x 10^{57} s^{-1}, respectively, with a factor of ~ 2 uncertainty in each case. Finally, the best--fit values of x_H are 0.16, 1.0, and 1.0, respectively. The uncertainty in the case of J1148+5251 is large, and x_H is not well constrained. However, for both J1030+0524 and J1623+3112, we find a significant lower limit of x_H >~ 0.033. Our method is different from previous analyses of the GP absorption spectra of these quasars, and our results strengthen the evidence that the rapid end--stage of reionization is occurring near z ~ 6.

Forthcoming large galaxy cluster surveys will yield tight constraints on cosmological models. It has been shown that in an idealized survey, containing > 10,000 clusters, statistical errors on dark energy and other cosmological parameters will be at the percent level. It has also been shown that through "self-calibration", parameters describing the mass-observable relation and cosmology can be simultaneously determined, though at a loss in accuracy by about an order of magnitude. Here we examine the utility of an alternative approach of self-calibration, in which a parametrized ab-initio physical model is used to compute cluster structure and the resulting mass-observable relations. As an example, we use a modified-entropy ("pre-heating") model of the intracluster medium, with the history and magnitude of entropy injection as unknown input parameters. Using a Fisher matrix approach, we evaluate the expected simultaneous statistical errors on cosmological and cluster model parameters. We study two types of surveys, in which a comparable number of clusters are identified either through their X-ray emission or through their integrated Sunyaev-Zel'dovich (SZ) effect. We find that compared to a phenomenological parametrization of the mass-observable relation, using our physical model yields significantly tighter constraints in both surveys, and offers substantially improved synergy when the two surveys are combined. These results suggest that parametrized physical models of cluster structure will be useful when extracting cosmological constraints from SZ and X-ray cluster surveys. (abridged)

We consider the pumping of the $63.2 \mu$m fine structure line of neutral OI in the high--redshift intergalactic medium (IGM), in analogy with the Wouthuysen--Field effect for the 21cm line of cosmic HI. We show that the soft UV background at $\sim 1300$\AA can affect the population levels, and if a significant fraction of the IGM volume is filled with ``fossil HII regions'' containing neutral OI, then this can produce a non--negligible spectral distortion in the cosmic microwave background (CMB). OI from redshift $z$ is seen in emission at $(1+z)63.2\mu$m, and between $7<z<10$ produces a mean spectral distortion of the CMB with a $y$--parameter of $y=(10^{-9} - 3\times10^{-8}) (Z/10^{-3}{\rm Z_{\odot}}) (I_{UV})$, where $Z$ is the mean metallicity of the IGM and $I_{UV}$ is the UV background at 1300\AA in units of $10^{-20}$ erg/s/Hz/cm$^2$/sr. Because O is in charge exchange equilibrium with H, a measurement of this signature can trace the metallicity at the end of the dark ages, prior to the completion of cosmic reionization and is complementary to cosmological 21cm studies. While future CMB experiments, such as Planck could constrain the metallicity to the $10^{-2} Z_{\odot}$ level, specifically designed experiments could potentially achieve a detection. Fluctuations of the distortion on small angular scale may also be detectable.

We study the potential of a large future weak-lensing survey to constrain dark energy properties by using both the number counts of detected galaxy clusters (sensitive primarily to density fluctuations on small scales) and tomographic shear-shear correlations (restricted to large scales). We use the Fisher matrix formalism, assume a flat universe and parameterize the equation of state of dark energy by w(a)=w_0+w_a(1-a), to forecast the expected statistical errors from either observable, and from their combination. We show that the covariance between these two observables is small, and argue that therefore they can be regarded as independent constraints. We find that when the number counts and the shear-shear correlations (on angular scales l < 1000) are combined, an LSST (Large Synoptic Survey Telescope)-like survey can yield statistical errors on (Omega_DE, w_0, w_a) as tight as (0.003, 0.03, 0.1). These values are a factor of 2-25 better than using either observable alone. The results are also about a factor of two better than those from combining number counts of galaxy clusters and their power spectrum.

We use the spectra of ~ 22,000 early-type galaxies, selected from the Sloan Digital Sky Survey, to infer the ages, metallicities and star formation histories of these galaxies. We find clear evidence of "downsizing", i.e. galaxies with larger velocity dispersion have older stellar populations. In particular, most early-type galaxies with velocity dispersion exceeding 200 km s-1 formed more than 90% of their current stellar mass at redshift z > 2.5. Therefore, star formation was suppressed around this redshift. We also show that chemical enrichment was rapid, lasting 1-2 Gyr and find evidence that [Fe/H] is sub-solar. We study the robustness of these results by comparing three different approaches: using (i) Lick absorption line indices; (ii) fitting a single-burst stellar population model to the whole spectrum (lines+continuum); and (iii) reconstructing the star formation and metallicity histories in multiple age-bins, providing a method to measure mass-weighted ages and metallicities. We find good agreement between the luminosity-weighted ages and metallicities computed with these three methods.

We use the spectra of 22,000 nearby early-type galaxies from the Sloan Digital Sky Survey (SDSS) to determine the age distribution of these galaxies as a function of their velocity dispersion sigma_v in the range 100 km/s < sigma_v < 280 km/s. We then combine the inferred age-distributions with the local abundance of spheroids, including early-type galaxies and late-type bulges, to predict the evolution of the quasar luminosity function (LF) in the redshift range 0<z<6. We make the following simple assumptions: (i) the formation of stars in each galaxy, at the epoch identified with the mean mass-weighted stellar age, is accompanied by the prompt assembly of the nuclear supermassive black hole (SMBH); (ii) the mass of the SMBH obeys the M_bh-sigma_v correlation observed in nearby galaxies; (iii) the SMBH radiates at a fraction f_Edd of the Eddington limit for a fixed duration t_Q, and is identified as a luminous quasar during this epoch, (iv) the intrinsic dispersions in the Eddington ratio and the M_bh-sigma_v relation produce a combined scatter of Delta(log L_Q) around the mean logarithmic quasar luminosity <log L_Q> at fixed sigma_v. These assumptions require that the SMBH remnants of quasars with bolometric luminosity below L_bol=10^{12.5} f_Edd L_sun reside predominantly in bulges of late type galaxies. We find that evolution of the observed quasar LF can be fit over the entire redshift range in this simple model, 0<z<6 with the choices of Delta(log L_Q)=0.6-0.9, t_Q= (6-8)x10^7 yr, and <f_Edd>=0.3-0.5. We find no evidence that any of the model parameters evolves with redshift, supporting the strong connection between the formation of stars and nuclear SMBHs in spheroids.

We model the spectra and surface brightness distributions for the Lyman alpha (Lya) radiation expected from protogalaxies that are caught in the early stages of their assembly. We use the results of a companion paper to characterize the radiation emerging from spherically collapsing gas clouds. We then modify the intrinsic spectra to incorporate the effect of subsequent resonant scattering in the intergalactic medium (IGM). Using these models, we interpret a number of recent observations of extended Lya blobs (LABs) at high redshift. We suggest, based on the angular size, energetics, as well as the relatively shallow surface brightness profiles, and double-peaked spectra, that several of these LABs may be associated with collapsing protogalaxies. We suggest two follow-up observations to diagnose the presence of gas infall. High S/N spectra of LABs should reveal a preferential flattening of the surface brightness profile at the red side of the line. Complementary imaging of the blobs at redshifted Balmer alpha wavelengths should reveal the intrinsic Lya emissivity and allow its separation from radiative transfer effects. We show that Lya scattering by infalling gas can reproduce the observed spectrum of Steidel et al's LAB2 as accurately as a recently proposed outflow model. Finally, we find similar evidence for infall in the spectra of point-like Lyman alpha emitters. The presence of scattering by the infalling gas implies that the intrinsic Lya luminosities, and derived quantities, such as the star-formation rate, in these objects may have been underestimated by about an order of magnitude.

We present Monte Carlo calculations of Lyman alpha (Lya) radiative transfer through collapsing gas clouds, representing proto-galaxies that are caught in the process of their assembly. Such galaxies produce Lya flux over an extended solid angle from a spatially extended Lya emissivity and/or from scattering effects. We study the effect of the gas distribution and kinematics, and of the Lya emissivity profile, on the emergent spectrum and surface brightness distribution. The emergent Lya spectrum is typically double-peaked and asymmetric. In practice, the blue peak is significantly enhanced and the red peak, in most cases, will be undetectable. The resulting effective blueshift, combined with scattering in the intergalactic medium, will render extended Lya emission from collapsing protogalaxies difficult to detect beyond redshift z=4. The surface brightness distribution is typically flat, and a strong wavelength dependence of its slope (with preferential flattening at the red side of the line) would be a robust indication that Lya photons are being generated (rather than just scattered) in a spatially extended region around the galaxy. We also find that for self-ionized clouds whose effective Lya optical depth is less than 10^3, infall and outflow models can produce nearly identical spectra and surface brightness distributions, and are difficult to distinguish from one another. The presence of deuterium with a cosmic abundance may produce a narrow but detectable dip in the spectra of systems with moderate hydrogen column densities, in the range 10^18-10^20 cm^-2. Finally, we present a new analytic solution for the emerging Lya spectrum in the limiting case of a static uniform sphere, extending previous solutions for static plane-parallel slabs.

The reionization history of the intergalactic medium (IGM) at high redshift (z > 6) was likely strongly shaped by several global feedback processes. Because the earliest ionizing sources formed at the locations of the rare density peaks, their spatial distribution was strongly clustered. Here we demonstrate that this clustering significantly boosts the impact of feedback processes operating at high redshift. We build a semi-analytical model to include feedback and clustering simultaneously, and apply this model to the suppression of star-formation in minihalos due to photoionization. The model is built on the excursion-set-based formalism of Furlanetto, Zaldarriaga and Hernquist (2004), which incorporates the clustering of ionizing sources, and which we here extend to include suppression of star formation in minihalos. We find that clustering increases the mean HII bubble size by a factor of several, and it dramatically increases the fraction of minihalos that are suppressed, by a factor of up to 60 relative to a randomly distributed population. This enhanced suppression can significantly reduce the electron scattering optical depth, as required by the three-year data from the Wilkinson Microwave Anisotropy Probe (WMAP). We argue that source clustering is likely to similarly boost the importance of a variety of other feedback mechanisms.

We use three-dimensional hydrodynamic simulations to investigate the effects of a transient photoionizing ultraviolet (UV) flux on the collapse and cooling of pregalactic clouds. These clouds have masses in the range 10^5 -10^7 M_sun, form at high redshifts (z>18), are assumed to lie within the short-lived cosmological HII regions around the first generation of stars. In addition, we study the combined effects of this transient UV flux and a persistent Lyman-Werner (LW) background from distant sources. In the absence of a LW background, we find that a critical specific intensity of J_UV ~ 0.1 x 10^-21 ergs s^-1 cm^-2 Hz^-1 sr^-1 demarcates a transition from net negative to positive feedback for the halo population. A weaker UV flux stimulates subsequent star formation inside the fossil HII regions, by enhancing the H_2 molecule abundance. A stronger UV flux significantly delays star-formation by reducing the gas density, and increasing the cooling time, at the centers of collapsing halos. At a fixed J_UV, the sign of the feedback also depends strongly on the density of the gas at the time of UV illumination. Regardless of the whether the feedback is positive or negative, we find that once the UV flux is turned off, its impact stars to diminish after ~30% of the Hubble time. In the more realistic case when a LW background is present, with J_LW > 0.01 x 10^-21 ergs s^-1 cm^-2 Hz^-1 sr^-1, strong suppression persists down to the lowest redshift (z=18) in our simulations. Finally, we find evidence that heating and photoevaporation by the transient UV flux renders the ~10^6 M_sun halos inside fossil HII regions more vulnerable to subsequent H_2 photo-dissociation by a LW background.

The primordial gas in the earliest dark matter halos, collapsing at redshifts around z=20, with masses M_halo=10^6 M_sun, and virial temperatures T_vir<10^4K, relied on the presence of molecules for cooling. Several theoretical studies have suggested that gas contraction and star-formation in these minihalos was suppressed by radiative, chemical, thermal, and dynamical feedback processes. The recent measurement by the Wilkinson Microwave Anisotropy Probe (WMAP) of the optical depth to electron scattering, tau=0.09+/-0.03, provides the first empirical evidence for this suppression. The new WMAP result is consistent with vanilla models of reionization, in which ionizing sources populate cold dark matter (CDM) halos down to a virial temperature of T_vir=10^4K. On the other hand, we show that in order to avoid overproducing the optical depth, the efficiency for the production of ionizing photons in minihalos must have been about an order of magnitude lower than expected and lower than the efficiency in large halos that can cool via atomic hydrogen (T_vir > 10^4K). This conclusion is insensitive to assumptions about the efficiency of ionizing photon production in the large halos, as long as reionization ends by z=6, as required by the spectra of bright quasars at z<6. Our conclusion is strengthened if the clumping of the ionized gas evolves with redshift, as suggested by semi-analytical predictions and three-dimensional numerical simulations.

Four recent observational results have challenged our understanding of high--redshift galaxies, as they require the presence of far more ultraviolet photons than should be emitted by normal stellar populations. First, there is significant ultraviolet emission from Lyman Break Galaxies (LBGs) at wavelenghts shorter than 912\AA. Second, there is strong Lyman alpha emission from extended ``blobs'' with little or no associated apparent ionizing continuum. Third, there is a population of galaxies with unusually strong Lyman-alpha emission lines. And fourth, there is a strong HeII (1640 \AA) emission line in a composite of LBGs. The proposed explanations for the first three observations are internally inconsistent, and the fourth puzzle has remained hitherto unexplained. Here we show that all four problems are resolved simultaneously if 10-30 percent of the stars in many galaxies at z ~ 3-4 are mainly primordial - unenriched by elements heavier than helium ('metals'). Most models of hierarchical galaxy formation assume efficient intra--galactic metal mixing, and therefore do not predict metal-free star formation at redshifts significantly below z ~5. Our results imply that micro-mixing of metals within galaxies is inefficient on a ~ Gyr time-scale, a conclusion that can be verified with higher resolution simulations, and future observations of the HeII emission line.

We model the number of detectable supernovae (SNe) as a function of redshift at different flux thresholds, making use of the observed properties of local SNe, such as their lightcurves, fiducial spectra, and peak magnitude distributions. We obtain a rate of 0.4-2.3 SNe/arcmin^2/yr at z > 5 at the near infrared (4.5 micron) flux density threshold of 3 nJy (achievable with the James Webb Space Telescope (JWST) in a 10^5 s integration). In a hypothetical one-year survey, it should be possible to detect up to several thousand SNe per unit redshift at z > 6. We discuss the possible application of such a large sample of distant SNe as a probe of the epoch of reionization. By heating the intergalactic medium and raising the cosmological Jeans mass, the process of reionization can suppress star formation in low-mass galaxies. This could have produced a relatively sharp drop in the SNR around the redshift of reionization. We quantify the detectability of this feature in future surveys of distant SNe by varying the redshift and duration of the feature, as well as its impact on the SFR in low-mass halos, which results in different redshifts, widths, and sizes of the drop in the expected SFR. We find that the drop can be detected out to z_re ~ 13, as long as (i) the reionization history contains a relatively rapid feature which is synchronized over different regions to within \Delta z < 1 - 3, (ii) the star-formation efficiency in halos that dominate reionization is ~ 10%, and (iii) reionization significantly suppresses the star formation in low-mass halos. Depending on the details of (i) - (iii), this could be achieved with a survey lasting less than two weeks. Detecting this signature would also help elucidate the feedback mechanism that regulates reionization.

The gravitational waves (GW) emitted during the coalescence of supermassive black holes (SMBHs) in the mass range 10^4-10^7 M_sun will be detectable out to high redshifts with LISA. We calculate the size and orientation of the three-dimensional error ellipse in solid angle and redshift within which the LISA event could be localized using the GW signatures alone. We take into account uncertainties in LISA's measurements of the luminosity distance and direction to the source, in the background cosmology, in weak gravitational lensing magnification due to inhomogeneities along the line of sight, and potential source peculiar velocities. We find that weak lensing errors exceed other sources of uncertainties by nearly an order of magnitude. Under the plausible assumption that BH mergers are accompanied by gas accretion leading to Eddington-limited quasar activity, we then compute the number of quasars that would be found in a typical LISA error volume, as a function of BH mass and redshift. We find that low redshifts offer the best opportunities to identify quasar counterparts to cosmological standard sirens, and that the LISA error volume will typically contain a single near-Eddington quasar at z=1. This will allow a straightforward test of the hypothesis that BH mergers are accompanied by bright quasar activity and, if the hypothesis proves correct, will guarantee the identification of a unique quasar counterpart. This would yield unprecedented tests of the physics of SMBH accretion, and offer an alternative method to precisely constrain cosmological parameters [abridged].

We report on i-band snapshot observations of 157 Sloan Digital Sky Survey (SDSS) quasars at 4<z<5.4 using the Advanced Camera for Surveys on the Hubble Space Telescope (HST) to search for evidence of gravitational lensing of these sources. None of the quasars appear to be strongly lensed and multiply imaged at the angular resolution (~0.1") and sensitivity of HST. The non-detection of strong lensing in these systems constrains the z=4-5 luminosity function to an intrinsic slope of beta>-3.8 (3 sigma), assuming a break in the quasar luminosity function at M_1450^star=-24.5. This constraint is considerably stronger than the limit of beta>-4.63 obtained from the absence of lensing in four z>5.7 quasars. Such constraints are important for our understanding of the true space density of high-redshift quasars and the ionization state of the early universe.

The amount and nature of dark energy (DE) can be tightly constrained by measuring the spatial correlation features and evolution of a sample of ~ 100,000 galaxy clusters over the redshift range 0<z < 1.5. Such an X-ray survey will discover all collapsed structures with mass above 3.5e14 solar masss at redshifts z<2 (i.e. the full range where such objects are expected) in the high Galactic latitude sky. Above this mass threshold the tight correlations between X-ray observables and mass allow direct interpretation of the data. We describe the constraints on Dark Energy that can be inferred from such a survey, using powerful self-calibration techniques to relate the X-ray observables (luminosity and temperature) to the underlying mass.

We revisit the possibility of detecting an extrasolar planet around a background star as it crosses the fold caustic of a foreground binary lens. During such an event, the planet's flux can be magnified by a factor of ~100 or more. The detectability of the planet depends strongly on the orientation of its orbit relative to the caustic. If the source star is inside the inter-caustic region, detecting the caustic-crossing planet is difficult against the magnified flux of its parent star. In the more favorable configuration, when the star is outside the inter-caustic region when the planet crosses the caustic, a close-in Jupiter-like planet around a Sun-like star at a distance of 8 kpc is detectable in 8-minute integrations with a 10m telescope at maximal S/N~15 for phase angle ~10 degrees. In this example, we find further that the presence of methane, at its measured abundance in Jupiter, and/or water, sodium and potassium, at the abundances expected in theoretical atmosphere models of close-in Jupiters, could be inferred from a non-detection of the planet in strong broad absorption bands at 0.6-1.4 microns caused by these compounds, accompanied by a S/N~10 detection in adjacent bands. We conclude that future generations of large telescopes might be able to probe the composition of the atmospheres of distant extrasolar planets.

Large future galaxy cluster surveys, combined with cosmic microwave background observations, can achieve a high sensitivity to the masses of cosmologically important neutrinos. We show that a weak lensing selected sample of ~100,000 clusters could tighten the current upper bound on the sum of masses of neutrino species by an order of magnitude, to a level of 0.03 eV. Since this statistical sensitivity is below the best existing lower limit on the mass of at least one neutrino species, a future detection is likely, provided that systematic errors can be controlled to a similar level.

Cosmological models that include suppression of the power spectrum of density fluctuations on small scales exhibit an exponential reduction of high-redshift, non-linear structures, including a reduction in the rate of gamma ray bursts (GRBs). Here we quantify the constraints that the detection of distant GRBs would place on structure formation models with reduced small-scale power. We compute the number of GRBs that could be detectable by the Swift satellite at high redshifts (z > 6), assuming that the GRBs trace the cosmic star formation history, which itself traces the formation of non-linear structures. We calibrate simple models of the intrinsic luminosity function of the bursts to the number and flux distribution of GRBs observed by the Burst And Transient Source Experiment (BATSE). We find that a discovery of high-z GRBs would imply strong constraints on models with reduced small-scale power. For example, a single GRB at z > 10, or 10 GRBs at z > 5, discovered by Swift during its scheduled two-year mission, would rule out an exponential suppression of the power spectrum on scales below R_c=0.09 Mpc (exemplified by warm dark matter models with a particle mass of m_x=2 keV). Models with a less sharp suppression of small-scale power, such as those with a red tilt or a running scalar index, n_s, are more difficult to constrain, because they are more degenerate with an increase in the power spectrum normalization, sigma_8, and with models in which star-formation is allowed in low-mass minihalos. We find that a tilt of \delta n_s ~ 0.1 is difficult to detect; however, an observed rate of 1 GRB/yr at z > 12 would yield an upper limit on the running of the spectral index, alpha = d(n_s)/d(ln k) > -0.05.

We explore the possibility of resolving an image of a damped Lyman alpha (DLA) system in absorption against an extended, diffuse background X-ray source. Typical columns of neutral hydrogen in DLAs are high enough to block out up to ~30% of the soft X-ray flux at an observed photon energy of 0.5 keV, and we find that ~ 1% of the area of extended X-ray sources at z > 1 have their 0.5 keV flux reduced by at least 20%. We discuss the observability of such absorption and find that < 2 arcsecond resolution, and > 300 photons per angular resolution element are required in the 0.3-8 keV band for its detection, and in order to distinguish it from intrinsic surface brightness fluctuations. For the surface brightness of the currently known high-redshift extended X-ray sources, this requires an integration time of a few Msec on Chandra. The detection will be within the reach of a routine observation with a next generation X-ray telescope such as XEUS or Generation X.

We present grism spectra taken with the Advanced Camera for Surveys to identify 29 red sources with (i'-z') > 0.9 in the Hubble Ultra Deep Field (HUDF). Of these 23 are found to be galaxies at redshifts between z=5.4 and 6.7, identified by the break at 1216 angstrom due to IGM absorption; two are late type dwarf stars with red colors; and four are galaxies with colors and spectral shape similar to dust reddened or old galaxies at redshifts z~1-2. This constitutes the largest uniform, flux-limited sample of spectroscopically confirmed galaxies at such faint fluxes z'<27.5. Many are also among the most distant spectroscopically confirmed galaxies (at redshifts up to z=6.7). We find a significant overdensity of galaxies at redshifts z=5.9 +/- 0.2. Nearly two thirds of the galaxies in our sample (15/23) belong to this peak. Taking into account the selection function and the redshift sensitivity of the survey, we get a conservative overdensity of at least a factor of two along the line-of-sight. The galaxies found in this redshift peak are also localized in the plane of the sky in a non-random manner, occupying about half of the ACS chip. Thus the volume overdensity is a factor of four. The star-formation rate derived from detected sources in this overdense region is sufficient to reionize the local IGM.

We show that the abundance and redshift distribution ($dN/dz$) of galaxy clusters in future high--yield cluster surveys, combined with the spatial power spectrum ($P_c(k)$) of the same clusters, can place significant constraints on the evolution of the dark energy equation of state, $w=w(a)$. We evaluate the expected errors on $w_a=-dw/da$ and other cosmological parameters using a Fisher matrix approach, and simultaneously including cluster structure evolution parameters in our analysis. We study three different types of forthcoming surveys that will identify clusters based on their X-ray emission (such as DUO, the Dark Universe Observatory), their Sunyaev--Zel'dovich (SZ) decrement (such as SPT, the South Pole Telescope), or their weak lensing (WL) shear (such as LSST, the Large Synoptic Survey Telescope). We find that combining the cluster abundance and power spectrum significantly enhances constraints from either method alone. We show that the weak-lensing survey can deliver a constraint as tight as $\Delta w_a \sim 0.1$ on the evolution of the dark energy equation of state, and that the X-ray and SZ surveys each yield $\Delta w_a \sim 0.4$ separately, or $\Delta w_a \sim 0.2$ when these two surveys are combined. For the X-ray and SZ surveys, constraints on dark energy parameters are improved by a factor of two by combining the cluster data with cosmic microwave background (CMB) anisotropy measurements by Planck, but degrade by a factor of two if the survey is required to solve simultaneously for cosmological and cluster structure evolution parameters. The constraint on $w_a$ from the weak lensing survey is improved by $\sim 25$% with the addition of Planck data.

We have discovered an extended Lyman alpha plume associated with a compact source at redshift 5.4 in slitless spectroscopic data from the Grism ACS Program for Extragalactic Science (GRAPES) project. The spatial extent of the emission is about 6 x 1.5 kpc (1 x 0.25 arcsec). Combining our grism data and the broadband images from the Hubble UltraDeep Field (UDF) images, we find a Lyman alpha line flux of 2e-17 erg/cm2/s and surface brightness 7e-17 erg/cm2/s/arcsec2. The UDF images show diffuse continuum emission associated with the Lyman alpha plume (hereafter UDF 5225), with three embedded knots. The morphology of UDF 5225 is highly suggestive of a galaxy in assembly. It is moreover possible that the prominent Lyman alpha emission from this object is due to an active nucleus, and that we are seeing the simultaneous growth through accretion of a galaxy and its central black hole. Followup observations at higher spectral resolution could test this hypothesis.

Gravitational lensing can amplify the apparent brightness of distant sources. Images that are highly magnified are often part of multiply-imaged systems, but we consider the possibility of having large magnifications without additional detectable images. In rare but non-negligible situations, lensing can produce a singly highly magnified image; this phenomenon is mainly associated with massive cluster-scale halos (>~1e13.5 Msun). Alternatively, lensing can produce multiply-imaged systems in which the extra images are either unresolved or too faint to be detectable. This phenomenon is dominated by galaxies and lower-mass halos (<~1e12 Msun), and is very sensitive to the inner density profile of the halos. Although we study the general problem, we customize our calculations to four quasars at redshift z~6 in the Sloan Digital Sky Survey (SDSS), for which Richards et al. (2004) have ruled out the presence of extra images down to an image splitting of 0.3" and a flux ratio of f=0.01. We predict that 9-29% of all z~6 quasars that are magnified by a factor of mu>10 would lack detectable extra images, with 5-10% being true singly-imaged systems. The maximum of 29% is reached only in the unlikely event that all low-mass (<~1e10 Msun) halos have highly concentrated (isothermal) profiles. In more realistic models where dwarf halos have flatter (NFW) inner profiles, the maximum probability is ~10%. We conclude that the probability that all four SDSS quasars are magnified by a factor of 10 is <~1e-4. The only escape from this conclusion is if there are many (>10) multiply-imaged z~6 quasars in the SDSS database that have not yet been identified, which seems unlikely. In other words, lensing cannot explain the brightnesses of the z~6 quasars, and models that invoke lensing to avoid having billion-Msun black holes in the young universe are not viable.

At redshifts beyond z>6, as the mean fraction of neutral hydrogen x_HI in the intergalactic medium (IGM) increases, the line flux of Lyman alpha (Lya) emitters can be significantly suppressed, which can result in a decrease in the observed number of emitters at a given Lya flux. However, cosmological HII regions surrounding the Lya emitting galaxies alleviate these effects. We use simple models of the Lya line suppression that incorporate the presence of HII regions to predict the overall effect of the Lya absorption on the Lya luminosity function. We find, in agreement with other recent studies, that unless ionizing sources are unusually strongly clustered, a fully neutral IGM may be inconsistent with the large abundance of confirmed z=6.5 Lya emitters. However, the presence of local HII regions prohibits placing a tight constraint on the mean neutral fraction. We find x_HI < 0.25; the presence of strong winds and/or the clustering of ionizing sources would further weaken this constraint. We conclude that the evolution of the Lya LF is consistent with reionization occurring near this redshift. Finally, we suggest that a measurement of observed Lya line width as a function of the Lya luminosity, in a future, larger sample of Lya emitters, may serve as a robust diagnostic of the neutral fraction in the IGM.

In cosmological structure formation models, massive non-linear objects in the process of formation, such as galaxy clusters, are surrounded by large-scale shocks at or around the expected virial radius. Direct observational evidence for such virial shocks is currently lacking, but we show here that their presence can be inferred from future, high resolution, high-sensitivity observations of the Sunyaev-Zel'dovich (SZ) effect in galaxy clusters. We study the detectability of virial shocks in mock SZ maps, using simple models of cluster structure (gas density and temperature distributions) and noise (background and foreground galaxy clusters projected along the line of sight, as well as the cosmic microwave background anisotropies). We find that at an angular resolution of 2'' and sensitivity of 10 micro K, expected to be reached at ~ 100 GHz frequencies in a ~ 20 hr integration with the forthcoming ALMA instrument, virial shocks associated with massive M ~ 10^15 M_Sun clusters will stand out from the noise, and can be detected at high significance. More generally, our results imply that the projected SZ surface brightness profile in future, high-resolution experiments will provide sensitive constraints on the density profile of cluster gas.

We discuss currently available observational constraints on the reionization history of the intergalactic medium (IGM), and the extent to which accreting black holes (BHs) can help explain these observations. We show new evidence, based on the combined statistics of Lyman alpha and beta absorption in quasar spectra, that the IGM contains a significant amount of neutral hydrogen, and is experiencing rapid ionization around redshift z=6. However, we argue that quasar BHs, even faint ones that are below the detection thresholds of existing optical surveys, are unlikely to drive the evolution of the neutral fraction around this epoch, because they would over-produce the present-day soft X-ray background. On the other hand, the seeds of the quasar BHs around z=6 likely appeared at much earlier epochs (around z=20), and produced hard ionizing radiation by accretion. These early BHs are promising candidates to account for the high redshift ionization (around z=15) implied by the recent cosmic microwave anisotropy data from WMAP. Using a model for the growth of BHs by accretion and mergers in a hierarchical cosmology, we suggest that the early growth of quasars must include a super-Eddington growth phase, and that, although not yet optically identified, the FIRST radio survey may have already detected several thousand BHs at z>6 with Mbh>10^8 Msun.

Which came first, the stars and gas that make up a galaxy, or the giant black hole at its centre? Observations of a distant galaxy, caught as it forms, could help solve this chicken-and-egg problem.

A bright quasar residing in a dense and largely neutral intergalactic medium (IGM) at high redshifts (z > 6) will be surrounded by a large cosmological Stromgren sphere. The quasar's spectrum will then show a sharp increase in resonant Lyman line absorption at wavelengths approaching and shorter than that corresponding to the Stromgren sphere's boundary along the line of sight. We show here that simultaneously considering the measured absorption in two or more hydrogen Lyman lines can provide the dynamical range required to detect this feature. We model broad and robust features of the Lyman alpha and Lyman beta regions of the spectrum of the z=6.28 quasar SDSS J1030+0524, using a hydrodynamical simulation. From the steep wavelength-dependence of the inferred absorption opacity, we detect the boundary of the Stromgren sphere at a proper distance of 6.0 +/- 0.2 Mpc away from the source redshift. From the spectrum alone, we also find that beyond this distance, cosmic hydrogen turns nearly neutral, with a neutral fraction of x_HI > 0.2, and that the ionizing luminosity of this quasar is in the range (5.2 +/- 2.5) times 10^{56} photons/sec. The method presented here, when applied to future quasars, can probe the complex topology of overlapping ionized regions, and can be used to study the details of the reionization process.

Recent improved determinations of the mass density rho_BH of supermassive black holes (SMBHs) in the local universe have allowed accurate comparisons of rho_BH with the amount of light received from past quasar activity. These comparisons support the notion that local SMBHs are ``dead quasars'' and yield a value epsilon >~ 0.1 for the average radiative efficiency of cosmic SMBH accretion. BH coalescences may represent an important component of the quasar mass assembly and yet not produce any observable electromagnetic signature. Therefore, ignoring gravitational wave (GW) emission during such coalescences, which reduces the amount of mass locked into remnant BHs, results in an overestimate of epsilon. Here, we put constraints on the magnitude of this bias. We calculate the cumulative mass loss to GWs experienced by a representative population of BHs during repeated cosmological mergers, using loss prescriptions based on detailed general relativistic calculations. Despite the possibly large number of mergers in the assembly history of each individual SMBH, we find that near--equal mass mergers are rare, and therefore the cumulative loss is likely to be modest, amounting at most to an increase by 20 percent of the inferred epsilon value. Thus, recent estimates of epsilon >~ 0.1 appear robust. The space interferometer LISA should provide empirical constraints on the dark side of quasar evolution, by measuring the masses and rates of coalescence of massive BHs to cosmological distances.

We quantify and discuss the footprints of neutral hydrogen in the intergalactic medium (IGM) on the spectra of high-redshift (z ~ 6) sources, using mock spectra generated from hydrodynamical simulations of the IGM. We show that it should be possible to extract relevant parameters, including the mean neutral fraction in the IGM, and the radius of the local cosmological Stromgren region, from the flux distribution in the observed spectra of distant sources. We focus on quasars, but a similar analysis is applicable to galaxies and gamma ray burst (GRB) afterglows. We explicitly include uncertainties in the spectral shape of the assumed source template near the Lyman alpha line. Our results suggest that a mean neutral hydrogen fraction, x(HI) of unity can be statistically distinguished from x(HI)<0.01, by combining the spectra of tens of bright (M = -27) quasars. Alternatively, the same distinction can be achieved using the spectra of several hundred sources that are ~100 times fainter. Furthermore, if the radius of the Stromgren sphere can be independently constrained to within ~10 percent, this distinction can be achieved using a single source. The information derived from such spectra will help in settling the current debate as to what extent the universe was reionized at redshifts near z=6.

Two spectacular results emerging from recent studies of nearby dead quasars and distant active quasars are (i) the existence of tight relations between the masses of black holes (BHs) and the properties of their host galaxies (spheroid luminosity or velocity dispersion, galaxy mass), and (ii) a consistency between the local mass density in BHs and that expected by summing up the light received from distant active quasars. These results are partly shaped by successive galactic mergers and BH coalescences, since mergers redistribute the population of BHs in galaxies and BH binary coalescences reduce the mass density in BHs through losses to gravitational waves. Here, we isolate and quantify these effects by following the cosmological merger history of a population of massive BHs representing the quasar population between z = 3 and z =0. Our results suggest that the relation between BH mass and host galaxy properties, and inferences on the global efficiency of BH accretion during active quasar phases, could be influenced by the cumulative effect of repeated mergers.

We present deep unbiased spectroscopy of the Hubble Ultra Deep Field (UDF) carried out using the slitless grism spectroscopy mode of the Advance Camera for Surveys on the Hubble Space Telescope (HST). The GRIsm ACS Program for Extragalactic Science (GRAPES) achieves continuum detection as faint as $z_{AB}=27.2$ using 40 orbits ($9.2 \times 10^4$ seconds) on HST. The data were taken at four orientation angles to correct for the overlap of spectra. GRAPES data provide a unique, uninterrupted, low resolution (R=100) spectral coverage for $5500\AA < \lambda < 10500\AA$, and allow us to detect high redshift galaxies at $4 < z < 7$ whether they have $\lya$ lines or just show the Lyman Break, as well as find low luminosity AGNs in an unbiased fashion. This paper describes in detail the observations and the data reduction, and examines the quality of the extracted spectra. Subsequent papers will deal with the analysis of the data. The extracted and calibrated GRAPES spectra will be available from MAST at STScI.

Krstic has carried out the first quantum mechanical calculations near threshold for the charge transfer (CT) process H^+ + H_2(X ^1Sigma_g^+, nu=0, J=0) --> H(1s) + H_2^+. These results are relevant for models of primordial galaxy and first star formation that require reliable atomic and molecular data for obtaining the early universe hydrogen chemistry. Using the results of Krstic, we calculate the relevant CT rate coefficient for temperatures between 100 and 30,000 K. We also present a simple fit which can be readily implemented into early universe chemical models. Additionally, we explore how the range of previously published data for this reaction translates into uncertainties in the predicted gas temperature and H_2 relative abundance in a collapsing primordial gas cloud. Our new data significantly reduce these cosmological uncertainties that are due to the uncertainties in the previously published CT rate coefficients.

Recent studies have shown that during their coalescence, binary supermassive black holes (SMBHs) experience a gravitational recoil with velocities of 100 km/s < v(kick) < 600 km/s. These velocities exceed the escape velocity v(esc) from typical dark matter (DM) halos at high-redshift (z>6), and therefore put constraints on scenarios in which early SMBHs grow at the centers of DM halos. Here we quantify these constraints for the most distant known SMBHs, with inferred masses in excess of 10^9 M(sun), powering the bright quasars discovered in the Sloan Digital Sky Survey at z>6. We assume that these SMBHs grew via a combination of accretion and mergers between pre-existing seed BHs in individual progenitor halos, and that mergers between progenitors with v(esc) < v(kick) disrupt the BH growth process. Our results suggest that under these assumptions, the z=6 SMBHs had a phase during which gained mass significantly more rapidly than under an Eddington-limited exponential growth rate.

The recently discovered z=10 galaxy (Pello et al 2004) has a strong Lyman alpha emission line that is consistent with being surprisingly symmetric, even given the relatively poor quality of its spectrum. The blue wing of a Lyman alpha line originating at high redshift should be strongly suppressed by resonant hydrogen absorption along the line of sight, an expectation borne out by the observed asymmetric shapes of the existing sample of Lyman alpha emitting sources at lower redshifts (3< z < 6.7). Absorption on the blue side of the line of the Pello et al. source could be reduced if the intergalactic medium (IGM) in the vicinity of the galaxy is highly ionized, but we show that this requires an unrealistically high ionizing emissivity. We suggest instead that the Lyman alpha emitting gas be receding relative to the surrounding gas with a velocity of >= 35km/s, a large velocity that is plausible only if the galaxy is part of a larger system (group of galaxies) with a velocity dispersion >= 35km/s. We find that with this velocity shift, the observed strength and shape of the line is still consistent with the galaxy being surrounded by its own Stromgren sphere embedded in a fully neutral IGM. More generally, we predict that at any given redshift, the bright Lyman alpha emitters with broader lines would exhibit stronger asymmetry than fainter ones. Bright galaxies with symmetric Lyman alpha lines may be signposts for groups and clusters of galaxies, within which they can acquire random velocities comparable to or larger than their linewidths.

The first massive astrophysical black holes likely formed at high redshifts (z>10) at the centers of low mass (~10^6 Msun) dark matter concentrations. These black holes grow by mergers and gas accretion, evolve into the population of bright quasars observed at lower redshifts, and eventually leave the supermassive black hole remnants that are ubiquitous at the centers of galaxies in the nearby universe. The astrophysical processes responsible for the formation of the earliest seed black holes are poorly understood. The purpose of this review is threefold: (1) to describe theoretical expectations for the formation and growth of the earliest black holes within the general paradigm of hierarchical cold dark matter cosmologies, (2) to summarize several relevant recent observations that have implications for the formation of the earliest black holes, and (3) to look into the future and assess the power of forthcoming observations to probe the physics of the first active galactic nuclei.

We use a physically motivated semi-analytic model, based on the mass function of dark matter halos, to predict the number of radio-loud quasars as a function of redshift and luminosity. Simple models in which the central BH mass scales with the velocity dispersion of its host halo as M(bh) sigma(halo)^5 have been previously found to be consistent with a number of observations, including the optical and X-ray quasar luminosity functions. We find that similar models, when augmented with an empirical prescription for radio emission, overpredict the number of faint (10 micro-Jy) radio sources by 1-2 orders of magnitude. This translates into a more stringent constraint on the low-mass end of the quasar black hole mass function than is available from the Hubble and Chandra Deep Fields. We interpret this discrepancy as evidence that black holes with masses below 10^7 Msun are either rare or are not as radio-loud as their more massive counterparts. Models that exclude BHs with masses below 10^7 Msun are in agreement with the deepest existing radio observations, but still produce a significant tail of high-redshift objects. In the 1-10GHz bands, at the sensitivity of 10 micro-Jy, we find surface densities of 100, 10, and 0.3 deg^-2 for sources located at z>6, 10, and 15, respectively. The discovery of these sources with instruments such as the Allen Telescope Array (ATA), Extended Very Large Array (EVLA), and the Square Kilometer Array (SKA) would open a new window for the study of supermassive BHs at high redshift. We also find surface densities of 0.1 deg^-2 at z > 6 for mJy sources that can be used to study 21 cm absorption from the epoch of reionization. We suggest that, although not yet optically identified, the FIRST survey may have already detected 10^3-10^4 such sources.

A population of black holes (BHs) at high redshifts (z>6) that contributes significantly to the ionization of the intergalactic medium (IGM) would be accompanied by the copious production of hard (>10keV) X-ray photons. The resulting hard X-ray background would redshift and be observed as a present-day soft X-ray background (SXB). Under the hypothesis that BHs are the main producers of reionizing photons in the high-redshift universe, we calculate their contribution to the present-day SXB. Our results, when compared to the unresolved component of the SXB in the range 0.5-2 keV, suggest that accreting BHs (be it luminous quasars or their lower-mass ``miniquasar'' counterparts) did not dominate reionization. Distant miniquasars that produce enough X-rays to only partially ionize the IGM to a level of at most 50% are still allowed, but could be severely constrained by improved determinations of the unresolved component of the SXB.

We consider microlensing of an elliptical source crossing a fold caustic. We derive a simple expression for the light curve of a source with uniform surface brightness that is accurate to third order in the ellipticity e (yielding errors of <1% for e<0.3). We consider caustic crossings of a rotating star, an oblate giant extrasolar planet, and an inclined standard thin accretion disk, and find the following results. (1) If most stars have rotation periods similar to chromospherically-active (i.e. spotted) giant stars in the bulge, than ~15% should be sufficiently oblate (e>0.25) to produce detectable (>1%) deviations with current observations. The form of the deviation due to ellipticity is qualitatively similar to that due to limb-darkening. Thus stellar oblateness may have to be taken into account in interpretations of precise limb-darkening measurements with microlensing. (2) Giant planets will generally not produce a detectable oblateness signal, either because they rotate too slowly (if they are close-in and tidally locked), or because they are too faint (if they are farther away from the parent star). (3) There is a near-degeneracy between ellipticity and position angle for equal-area sources with scale-free intensity profiles and elliptical isophotes. This degeneracy results in a factor of ~2 uncertainty in the measurement of the scale-length of a standard thin accretion disk using observations of a single, linear fold crossing quasar microlensing event in passbands probing the outer parts of the disk. Higher-frequency observations can reduce this uncertainty.

GRBs are expected to occur at redshifts far higher than the highest quasar redshifts so far detected. And unlike quasars, GRB afterglows may provide "clean" probes of the epoch of reionization because no complications (such as the presence of a strong Ly-alpha emission line) or "proximity effects" (such as the Stromgren sphere produced by ionizing photons from the quasar) are expected. Thus NIR and optical observations of GRB afterglows may provide unique information about the epoch of reionization. In particular, the flux at wavelengths shortward of Ly-alpha provides a direct measure of the density fluctuations of the IGM at the GRB redshift, while the flux at wavelengths longward of Ly-alpha provides an integrated measure of the number of ionizing photons produced by stars in the host galaxy of the GRB up until the burst occurs. A comparison of the sizes of the Stromgren spheres produced by stars in the host galaxies of GRBs and by quasars then provides an estimate of the relative contributions of star formation and quasars to reionization. We use detailed calculations of the expected shape of the GRB afterglow spectrum in the vicinity of Ly-alpha for GRBs at a variety of redshifts to illustrate these points.

We demonstrate that constraints on cosmological parameters from the distribution of clusters as a function of redshift (dN/dz) are complementary to accurate angular diameter distance (D_A) measurements to clusters, and their combination significantly tightens constraints on the energy density content of the Universe. The number counts can be obtained from X-ray and/or SZ (Sunyaev-Zel'dovich effect) surveys, and the angular diameter distances can be determined from deep observations of the intra-cluster gas using their thermal bremsstrahlung X-ray emission and the SZ effect. We combine constraints from simulated cluster number counts expected from a 12 deg^2 SZ cluster survey and constraints from simulated angular diameter distance measurements based on the X-ray/SZ method assuming a statistical accuracy of 10% in the angular diameter distance determination of 100 clusters with redshifts less than 1.5. We find that Omega_m can be determined within about 25%, Omega_Lambda within 20%, and w within 16%. We show that combined dN/dz + D_A constraints can be used to constrain the different energy densities in the Universe even in the presence of a few percent redshift dependent systematic error in D_A. We also address the question of how best to select clusters of galaxies for accurate diameter distance determinations. We show that the joint dN/dz + D_A constraints on cosmological parameters for a fixed target accuracy in the energy density parameters are optimized by selecting clusters with redshift upper cut--offs in the range 0.5 < z < 1.

Over the last few years, the Sloan Digital Sky Survey (SDSS) has discovered several hundred quasars with redshift between 4.0 and 6.4. Including the effects of magnification bias, one expects a priori that an appreciable fraction of these objects are gravitationally lensed. We have used the Advanced Camera for Surveys on the Hubble Space Telescope to carry out a snapshot imaging survey of high-redshift SDSS quasars to search for gravitationally split lenses. This paper, the first in a series reporting the results of the survey, describes snapshot observations of four quasars at z = 5.74, 5.82, 5.99 and 6.30, respectively. We find that none of these objects has a lensed companion within 5 magnitudes with a separation larger than 0.3 arcseconds; within 2.5 magnitudes, we can rule out companions within 0.1 arcseconds. Based on the non-detection of strong lensing in these four systems, we constrain the z~6 luminosity function to a slope of beta>-4.63 (3 sigma), assuming a break in the quasar luminosity function at M_{1450}^*=-24.0. We discuss the implications of this constraint on the ionizing background due to quasars in the early universe. Given that these quasars are not highly magnified, estimates of the masses of their central engines by the Eddington argument must be taken seriously, possibly challenging models of black hole formation.

A remarkable result of the Wilkinson Microwave Anisotropy Probe (WMAP) observations is that the universe was significantly reionized at large redshifts. The standard explanation is that massive stars formed early and reionized the universe around redshift z=17. Here we explore an alternative possibility, in which the universe was reionized in two steps. An early boost of reionization is provided by a decaying sterile neutrino, whose decay products, relativistic electrons, result in partial ionization of the smooth gas. We demonstrate that a neutrino with a mass of m_nu ~ 200 MeV and a decay time of t ~ 4 * 10^{15} s can account for the electron scattering optical depth tau=0.16 measured by WMAP without violating existing astrophysical limits on the cosmic microwave and gamma ray backgrounds. Reionization is then completed by subsequent star formation at lower redshifts. This scenario alleviates constraints on structure formation models with reduced small-scale power, such as those with a running or tilted scalar index, or warm dark matter models.

The cosmic ultraviolet (UV) ionizing background impacts the formation of dwarf galaxies in the low-redshift universe (z=3) by suppressing gas infall into galactic halos with circular velocities up to v(circ)=75 km/s. Using a one-dimensional, spherically symmetric hydrodynamics code (Thoul & Weinberg 1995), we examine the effect of an ionizing background on low-mass galaxies forming at high redshifts (z>10). We find that the importance of photoionization feedback is greatly reduced, because (1) at high redshift, dwarf-galaxy sized objects can self-shield against the ionizing background, (2) collisional cooling processes at high redshift are more efficient, (3) the amplitude of the ionizing background at high redshift is lower, and (4) the ionizing radiation turns on when the perturbation that will become the dwarf galaxy has already grown to a substantial overdensity. We find that because of these reasons, gas can collect inside halos with circular velocities as low as v(circ)=10 km/s. This result has important implications for the reionization history of the universe.

We study contributions from inhomogeneous (patchy) reionization to arcminute scale ($1000 < \ell < 10,000$) cosmic microwave background (CMB) anisotropies. We show that inhomogeneities in the ionization fraction, rather than in the mean density, dominate both the temperature and the polarization power spectra. Depending on the ionization history and the clustering bias of the ionizing sources, we find that rms temperature fluctuations range from 2 $\mu$K to 8 $\mu$K and the corresponding values for polarization are over two orders of magnitude smaller. Reionization can significantly bias cosmological parameter estimates and degrade gravitational lensing potential reconstruction from temperature maps but not from polarization maps. We demonstrate that a simple modeling of the reionization temperature power spectrum may be sufficient to remove the parameter bias. The high-$\ell$ temperature power spectrum will contain some limited information about the sources of reionization.

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

In light of the recent inference of a high optical depth to Thomson scattering, tau, from the WMAP data we investigate the effects of extended periods of partial ionization and ask if the value of tau inferred by assuming a single sharp transition is an unbiased estimate. We construct and consider several representative ionization models and evaluate their signatures in the CMB. If tau is estimated with a single sharp transition we show that there can be a significant bias in the derived value (and therefore a bias in sigma8 as well). For WMAP noise levels the bias in tau is smaller than the statistical uncertainty, but for Planck or a cosmic variance limited experiment the tau bias could be much larger than the statistical uncertainties. This bias can be reduced in the ionization models we consider by fitting a slightly more complicated ionization history, such as a two-step ionization process. Assuming this two-step process we find the Planck satellite can simultaneously determine the initial redshift of reionization to +-2 and tau to +-0.01 Uncertainty about the ionization history appears to provide a limit of about 0.005 on how well tau can be estimated from CMB polarization data, much better than expected from WMAP but significantly worse than expected from cosmic-variance limits.

The cosmic microwave background (CMB) has provided a precise template for features in the linear power spectrum: the matter-radiation turnover, sound horizon drop, and acoustic oscillations. In a two dimensional power spectrum in redshift and angular space, the features appear as distorted rings, and yield simultaneous, purely geometric, measures of the Hubble parameter H(z) and angular diameter distance D_A(z) via an absolute version of the Alcock-Paczynski test. Employing a simple Fisher matrix tool, we explore how future surveys can exploit these rings of power for dark energy studies. High-z CMB determinations of H and D_A are best complemented at moderate to low redshift (z < 0.5) with a population of objects that are at least as abundant as clusters of galaxies. We find that a sample similar to that of the ongoing SDSS Luminous Red Galaxy (LRG) survey can achieve statistical errors at the ~5% level for D_A(z) and H(z) in several redshift bins. This, in turn, implies errors of sigma(w)=0.03-0.05 for a constant dark energy equation of state in a flat universe. Deep galaxy cluster surveys such as the planned South Pole Telescope (SPT) survey, can extend this test out to z~1 or as far as redshift followup is available. We find that the expected constraints are at the sigma(w)=0.04-0.08 level, comparable to galaxies and complementary in redshift coverage.

In cosmological models favored by current observations, the first astrophysical objects formed in dark matter halos at redshifts starting at z>20, and their properties were determined by primordial H_2 molecular chemistry. These protogalaxies were very abundant, but substantially less massive than typical galaxies in the local Universe. Extreme metal-poor stars, and massive black holes in their nuclei reionized the bulk of the hydrogen in the intergalactic medium. Reionization may have taken place over an extended redshift interval, ending around z=7. Observational probes of the process of reionization are afforded by studying the polarization of the cosmic microwave background anisotropies, as well as by studying the spectra and abundance of distant Lyman alpha-emitting galaxies. Here we review theoretical expectations on how and when the first galaxies formed, and summarize future observational prospects of probing hydrogen reionization.

We discuss a simple model for the growth of supermassive black holes (BHs) at the center of spheroidal stellar systems. In particular, we assess the hypotheses that (1) star formation in spheroids and BH fueling are proportional to one another, and (2) the BH accretion luminosity stays near the Eddington limit during luminous quasar phases. With the aid of this simple model, we are able to interpret many properties of the QSO luminosity function, including the puzzling steep decline of the characteristic luminosity from redshift z=2 to to z=0: indeed the residual star formation in spheroidal systems is today limited to a small number of bulges, characterized by stellar velocity dispersions a factor of 2-3 smaller those of the elliptical galaxies hosting QSOs at z > 2. A simple consequence of our hypotheses is that the redshift evolution of the QSO emissivity and of the star formation history in spheroids should be roughly parallel. We find this result to be broadly consistent with our knowledge of the evolution of both the global star formation rate, and of the evolution of the QSO emissivity, but we identify interesting discrepancies at both low and high redshifts, to which we offer tentative solutions. Finally, our hypotheses allow us to present a robust method to derive the duty cycle of QSO activity, based on the observed QSO luminosity function, and on the present-day relation between the masses of supermassive BHs and those of their spheroidal host stellar systems. The duty cycle is found to be substantially less than unity, with characteristic values in the range (3-6)x10^(-3), and we compute that the average bolometric radiative efficiency is epsilon=0.07. Finally, we find that the growth in mass of individual black holes at high redshift (z>2) can be dominated by mergers, and is therefore not necessarily limited by accretion.

The recent measurement by WMAP of a large electron scattering optical depth tau_e = 0.17 +- 0.04 is consistent with a simple model of reionization in which the intergalactic medium (IGM) is ionized at redshift z ~ 15, and remains highly ionized thereafter. Here, we show that existing measurements of the IGM temperature from the Lyman-alpha forest at z ~ 2 - 4 rule out this ``vanilla'' model. Under reasonable assumptions about the ionizing spectrum, as long as the universe is reionized before z = 10, and remains highly ionized thereafter, the IGM reaches an asymptotic thermal state which is too cold compared to observations. To simultaneously satisfy the CMB and forest constraints, the reionization history must be complex: reionization begins early at z >~ 15, but there must have been significant (order unity) changes in fractions of neutral hydrogen and/or helium at 6 < z < 10, and/or singly ionized helium at 4 < z < 10. We describe a physically motivated reionization model that satisfies all current observations. We also explore the impact of a stochastic reionization history and show that a late epoch of (HeII --> HeIII) reionization induces a significant scatter in the IGM temperature, but the scatter diminishes with time quickly. Finally, we provide an analytic formula for the thermal asymptote, and discuss possible additional heating mechanisms that might evade our constraints.

The recent discovery of a high optical depth tau to Thomson scattering from the WMAP data implies that significant reionization took place at redshifts z~15. This discovery has important implications for the sources of reionization, and allows, for the first time, constraints to be placed on physical reionization scenarios out to redshift z~20. Using a new suite of semi-analytic reionization models, we show that the high value of tau requires a surprisingly high efficiency epsilon of the first generation of UV sources for injecting ionizing photons into the intergalactic medium. We find that no simple reionization model can be consistent with the combination of the WMAP result with data from the z<6.5 universe. Satisfying both constraints requires either of the following: (i) H_2 molecules form efficiently at z~20, survive feedback processes, and allow UV sources in halos with virial temperatures below Tvir=10^4 K to contribute substantially to reionization, or (ii) the efficiency epsilon in halos with Tvir>10^4K decreased by a factor of ~ 30 between (z~20) and (z~6). We discuss the relevant physical issues to produce either scenario, and argue that both options are viable, and allowed by current data. In detailed models of the reionization history, we find that the evolution of the ionized fractions in the two scenarios have distinctive features that Planck can distinguish at 3 sigma significance. At the high WMAP value for tau, Planck will also be able to provide tight statistical constraints on reionization model parameters, and elucidate much of the physics at the end of the Dark Ages. The sources responsible for the high optical depth discovered by WMAP should be directly detectable out to z~15 by the James Webb Space Telescope.

We present the discovery of three new quasars at z>6 in 1300 deg^2 of SDSS imaging data, J114816.64+525150.3 (z=6.43), J104845.05+463718.3 (z=6.23) and J163033.90+401209.6 (z=6.05). The first two objects have weak Ly alpha emission lines; their redshifts are determined from the positions of the Lyman break. They are only accurate to 0.05 and could be affected by the presence of broad absorption line systems. The last object has a Ly alpha strength more typical of lower redshift quasars. Based on a sample of six quasars at z>5.7 that cover 2870 deg^2 presented in this paper and in Paper I, we estimate the comoving density of luminous quasars at z 6 and M_{1450} < -26.8 to be (8 +/- 3)x10^{-10} Mpc^{-3} (for H_0 = 50 km/s/Mpc, Omega = 1). HST imaging of two z>5.7 quasars and high-resolution ground-based images (seeing 0.4'') of three additional z>5.7 quasars show that none of them is gravitationally lensed. The luminosity distribution of the high-redshfit quasar sample suggests the bright end slope of the quasar luminosity function at z 6 is shallower than Psi L^{-3.5} (2-sigma), consistent with the absence of strongly lensed objects.

The recent discovery of a Gunn--Peterson (GP) trough in the spectrum of the redshift 6.28 SDSS quasar has raised the tantalizing possibility that we have detected the reionization of the universe. However, a neutral fraction (of hydrogen) as small as 0.1% is sufficient to cause the GP trough, hence its detection alone cannot rule out reionization at a much earlier epoch. The Cosmic Microwave Background (CMB) polarization anisotropy offers an alternative way to explore the dark age of the universe. We show that for most models constrained by the current CMB data and by the discovery of a GP trough (showing that reionization occurred at z > 6.3), MAP can detect the reionization signature in the polarization power spectrum. The expected 1-sigma error on the measurement of the electron optical depth is around 0.03 with a weak dependence on the value of that optical depth. Such a constraint on the optical depth will allow MAP to achieve a 1-sigma error on the amplitude of the primordial power spectrum of 6%. MAP with two years (Planck with one year) of observation can distinguish a model with 50% (6%) partial ionization between redshifts of 6.3 and 20 from a model in which hydrogen was completely neutral at redshifts greater than 6.3. Planck will be able to distinguish between different reionization histories even when they imply the same optical depth to electron scattering for the CMB photons.

Galaxy cluster surveys provide a powerful means of studying the density and nature of the dark energy. The redshift distribution of detected clusters in a deep, large solid angle SZE or X-ray survey is highly sensitive to the dark energy equation of state. Accurate constraints at the 5% level on the dark energy equation of state require that systematic biases in the mass estimators must be controlled at better than the ~10% level. Observed regularity in the cluster population and the availability of multiple, independent mass estimators suggests these precise measurements are possible. Using hydrodynamical simulations that include preheating, we show that the level of preheating required to explain local galaxy cluster structure has a dramatic effect on X-ray cluster surveys, but only a mild effect on SZE surveys. This suggests that SZE surveys may be optimal for cosmology while X-ray surveys are well suited for studies of the thermal history of the intracluster medium.

Recent discoveries by the Sloan Digital Sky Survey (SDSS) of four bright redshift z=6 quasars could constrain the mechanism by which the supermassive black holes powering these sources are assembled. Here we compute the probability that the fluxes of the quasars are strongly amplified by gravitational lensing, and therefore the likelihood that the black hole masses are overestimated when they are inferred assuming Eddington luminosities. The poorly-constrained shape of the intrinsic quasar luminosity function (LF) at redshift z=6 results in a large range of possible lensing probabilities. If the LF is either steep, or extends to faint magnitudes, the probability for amplification by a factor of ten ore more (and with only one image detectable by SDSS) can reach essentially 100 percent. We show that future observations, in particular, either of the current four quasars at the high angular resolution provided by the Hubble Space Telescope, or of an increased sample of about twenty redshift six quasars at the current angular resolution, should either discover several gravitational lenses, or else provide interesting new constraints on the shape of the z=6 quasar LF.

For a source of Ly alpha radiation embedded in a neutral intergalactic medium (IGM) prior to the reionization epoch, the emission line is strongly suppressed by the intervening IGM. The damping wing of the so-called Gunn-Peterson trough can extend to the red side of the emission line, and erase a significant fraction of the total line flux. However, the transmitted fraction increases with the size of the local cosmological HII region surrounding the source, and therefore with the ionizing luminosity and age of the source. Motivated by the recent discovery of a Ly alpha emitting galaxy at a redshift z=6.56 (Hu et al. 2002), possibly prior to the reionization of the IGM, we revisit the effects of a neutral IGM on the Ly alpha emission line. We show that even for faint sources with little ionizing continuum, the emission line can remain observable. In particular, the line detected by Hu et al. is consistent with a source embedded in a neutral IGM. We provide characterizations of the asymmetry and total transmitted flux of the Ly alpha line as functions of the ionizing emissivity of its source. A statistical sample of Ly alpha emitters extending beyond the reionization redshift can be a useful probe of reionization.

The recent discovery of bright quasars around redshift z=6 suggests that black holes (BHs) with masses in excess of 10^9 Msun have already assembled at a very early stage in the evolution of the universe. An alternative interpretation is that these quasars are powered by less massive BHs, but their fluxes are strongly magnified through gravitational lensing by intervening galaxies. Here we analyze the flux distribution of the Ly alpha emission of the quasar with the highest known redshift, SDSS 1030+0524, at z=6.28. We show that this object could not have been magnified by lensing by more than a factor of five. The constraint arises from the large observed size, 30 (comoving) Mpc, of the ionized region around this quasar, and relies crucially only on the assumption that the quasar is embedded in a largely neutral IGM. Based on the line/continuum ratio of SDSS 1030+0524, we argue further that this quasar also cannot be beamed by a significant factor. We conclude that the minimum mass for its resident BH is 4 x 10^8 Msun (for magnification by a factor of five); if the mass is this low, then the quasars had to switch on prior to redshift z=9. From the size of the ionized region, we are also able to place an absolute lower bound on the age of this quasar at t > 2 x 10^7 years.

The first galaxies formed at high redshifts, and were likely substantially less massive than typical galaxies in the local universe. We argue that (1) the reionization of a clumpy intergalactic medium by redshift z=6, (2) its enrichment by metals by z=3 without disturbing the Lyman alpha forest, and (3) the presence of supermassive black holes powering the recently discovered bright quasars at z=6, strongly suggest that a population of low-mass galaxies exists beyond redshifts z>6. Although the first stars could have been born in dark matter halos with virial temperatures as low as Tvir=200K, collapsing as early as z=25, the first galaxies likely appeared in significant numbers only in halos with Tvir>10,000K that collapsed later (z=15). The gas in these more massive halos initially contracts isothermally to high densities by atomic Lyman alpha cooling. H2 molecules can then form efficiently via non-equilibrium gas-phase chemistry, allowing the gas to cool further to T=100K, and fragment on stellar mass scales. These halos can harbor the first generation of ``mini-galaxies'' that reionized the universe. The continuum and line emission from these sources, as well as their Lyman alpha cooling radiation, can be detected in the future by NGST and other instruments.

The first generation of protogalaxies likely formed out of primordial gas via H2-cooling in cosmological minihalos with virial temperatures of a few 1000K. However, their abundance is likely to have been severely limited by feedback processes which suppressed H2 formation. The formation of the protogalaxies responsible for reionization and metal-enrichment of the intergalactic medium, then had to await the collapse of larger halos. Here we investigate the radiative cooling and collapse of gas in halos with virial temperatures Tvir > 10^4K. In these halos, efficient atomic line radiation allows rapid cooling of the gas to 8000 K; subsequently the gas can contract nearly isothermally at this temperature. Without an additional coolant, the gas would likely settle into a locally gravitationally stable disk; only disks with unusually low spin would be unstable. However, we find that the initial atomic line cooling leaves a large, out-of-equilibrium residual free electron fraction. This allows the molecular fraction to build up to a universal value of about x(H2) = 10^-3, almost independently of initial density and temperature. We show that this is a non--equilibrium freezeout value that can be understood in terms of timescale arguments. Furthermore, unlike in less massive halos, H2 formation is largely impervious to feedback from external UV fields, due to the high initial densities achieved by atomic cooling. The H2 molecules cool the gas further to about 100K, and allow the gas to fragment on scales of a few 100 Msun. We investigate the importance of various feedback effects such as H2-photodissociation from internal UV fields and radiation pressure due to Ly-alpha photon trapping, which are likely to regulate the efficiency of star formation.

We present a robust method to derive the duty cycle of QSO activity based on the empirical QSO luminosity function and on the present-day linear relation between the masses of supermassive black holes and those of their spheroidal host stellar systems. It is found that the duty cycle is substantially less than unity, with characteristic values in the range $3-6\times 10^{-3}$. Finally, we tested the expectation that the QSO luminosity evolution and the star formation history should be roughly parallel, as a consequence of the above--mentioned relation between BH and galaxy masses.

Deep mid-infrared (MIR) observations could reveal a population of faint, high-redshift (z>3) dusty starburst galaxies that are the progenitors of present-day spheroids or bulges, and are beyond the reach of current instruments. We utilize a semi-analytic galaxy formation scheme to find an extreme model for the MIR galaxy counts, designed to maximize the number of detectable sources down to a flux level of a few nJy. The model incorporates the formation of heavily dust-enshrouded stellar populations at high redshift, and is consistent with existing observations, including faint counts at 1.6um in the NICMOS Hubble Deep Field, and the upper limit on the extragalactic MIR background from TeV gamma rays. Our models predict upto 0.5 galaxies/sq.arcsec at the threshold of 100 nJy at 6um, with a comparable or larger surface density at longer MIR wavelengths. We conclude that a significant new population of high-redshift galaxies could be detected by the Space Infrared Telescope Facility (SIRTF) and Next Generation Space Telescope (NGST). Such a population would constitute background noise for the Terrestrial Planet Finder (TPF), and could necessitate repeat observations: every {\it TPF} resolution element have a 10 percent chance of being contaminated by a background galaxy.

We present the results from a survey of i-dropout objects selected from ~1550 deg^2 of multicolor imaging data from the Sloan Digital Sky Survey, to search for luminous quasars at z>5.8. Objects with i*-z*>2.2 and z*<20.2 are selected, and follow-up J band photometry is used to separate L and T type cool dwarfs from high-redshift quasars. We describe the discovery of three new quasars, at z=5.82, 5.99 and 6.28, respectively. Their spectra show strong and broad Ly alpha+NV emission lines, and very strong Ly alpha absorption, with a mean continuum decrement D_A > 0.90. The ARC 3.5m spectrum of the z=6.28 quasar shows that over a range of 300 A immediately blueward of the Ly alpha emission, the average transmitted flux is only 0.003 +/-0.020 times that of the continuum level, consistent with zero flux, and suggesting a tentative detection of the complete Gunn-Peterson trough. The existence of strong metal lines suggests early chemical enrichment in the quasar enviornment. The three new objects, together with the previously published z=5.8 quasar form a complete color-selected flux-limited sample at z>5.8. We estimate that at $z=6$, the comoving density of luminous quasars at M_1450 < -26.89 (h=0.5, Omega=1)is 1.1x10^-9 Mpc^-3. This is a factor of ~2 lower than that at z~5, and is consistent with an extrapolation of the observed quasar evolution at low-z. We discuss the contribution of quasars to the ionizing background at z~6. The luminous quasars discussed in the paper have central black hole masses of several times 10^9 M_sun by the Eddington argument. Their observed space density provides a sensitive test of models of quasar and galaxy formation at high redshift. (Abridged)

We show that the counts of galaxy clusters in future deep cluster surveys can place strong constraints on the matter density, Omega_m, the vacuum energy density, Omega_L, and the normalization of the matter power spectrum, sigma_8. Degeneracies between these parameters are different from those in studies of either high--redshift type Ia Supernovae (SNe), or cosmic microwave background (CMB) anisotropies. Using a mass threshold for cluster detection expected to be typical for upcoming SZE surveys, we find that constraints on Omega_m and sigma_8 at the level of roughly 5% or better can be expected, assuming redshift information is known at least to z=0.5 and in the absence of significant systematic errors. Without information past this redshift, Omega_L is constrained to 25%. With complete redshift information, deep (M_{lim}= 10^{14}h^{-1}{M_sun}), relatively small solid angle (roughly 12 {deg}^2) surveys can further constrain Omega_L to an accuracy of 15%, while large solid angle surveys with ground-based large-format bolometer arrays could measure Omega_L to a precision of 4% or better.

An extensive sample of galaxy clusters will be available in the coming years, detected through their Sunyaev-Zeldovich effect (SZE). We use a semi-analytic model to study the scientific yield of combining SZE data with X-ray and optical follow-up observations. If clusters at a given redshift z_o can be identified with virialized, spherical halos, they populate a well-defined "fundamental plane" (FP) in the parameter space of the three observables: virial temperature (T), total Sunyaev-Zeldovich flux decrement (DS), and angular size (theta). The location and orientation of the FP, as well as its redshift-evolution, are sensitive to both the internal evolution of clusters, and to the underlying cosmological parameters. We show that if clusters are not standard candles (e.g. due to feedback, or energy injection), then this can be inferred from the FP. Likewise, we study the dependence of the FP on the cosmological parameters h, sigma_8, and Omega_0, and quantify constraints on these parameters. We also show that in the absence of any non-gravitational effects, the scatter in the the (DS-T) plane is significantly smaller than in either the (theta-T) or the (theta-DS) planes. As a result, the (DS-T) relation can be an exceptionally sensitive probe of both cluster physics and cosmological parameters, and a comparison of the amount of scatter in these three scaling relations will serve as a test of the origin (cosmological vs. stochastic) of the scatter.

We study the constraints that high-redshift structure formation in the universe places on warm dark matter (WDM) dominated cosmological models. We modify the extended Press-Schechter formalism to derive the halo mass function in WDM models. We show that our predictions agree with recent numerical simulations at low redshift over the halo masses of interest. Applying our model to galaxy formation at high redshift, we find that the loss of power on small scales, together with the delayed collapse of low-mass objects, results in strong limits on the root-mean-square velocity dispersion v_rms of the WDM particles at z=0. For fermions decoupling while relativistic, these limits are equivalent to constraints on the mass m_X of the particles. The presence of a 4 billion solar mass black hole at z=5.8, believed to power the quasar SDSS 1044-1215, implies m_X > 0.5 keV (or v_rms < 0.10 km/s), assuming that the quasar is unlensed and radiating at or below the Eddington limit. Reionization by z=5.8 also implies a limit on m_X. If high-redshift galaxies produce ionizing photons with an efficiency similar to their z=3 counterparts, we find m_X > 1.2 keV (or v_rms < 0.03 km/s). However, given the uncertainties in current measurements from the proximity effect of the ionizing background at redshift 3, values of m_X as low as 0.75 keV (v_rms = 0.06 km/s) are not ruled out. The limit weakens further if, instead, the ionizing-photon production efficiency is greater at high z, but this limit will tighten considerably if reionization is shown in the future to have occurred at higher redshifts. WDM models with m_X < 1 keV (v_rms > 0.04 km/s) produce a low-luminosity cutoff in the high-redshift galaxy luminosity function which is directly detectable with the Next Generation Space Telescope (abridged).

The hardness of the ionizing continuum from the first sources of UV radiation plays a crucial role in the reionization of the intergalactic medium (IGM). While usual stellar populations have soft spectra, mini-quasars or metal-free stars with high effective temperatures may emit hard photons, capable of doubly ionizing helium and increasing the IGM temperature. Absorption within the source and in the intervening IGM will render the ionizing continuum of high-redshift sources inaccessible to direct observation. Here we show that HeII recombination lines from the first luminous objects are potentially detectable by the Next Generation Space Telescope. Together with measurements of the hydrogen Balmer alpha emission line, this detection can be used to infer the ratio of HeII to HI ionizing photons. A measurement of this ratio would shed light on the nature and emission mechanism of the first luminous sources, with important astrophysical consequences for the reheating and reionization of the IGM.

The ubiquity of supermassive black holes (SMBHs) at the centers of nearby luminous galaxies can arise from the multiple mergers experienced by dark matter halos in hierarchical structure formation models, even if only a small fraction of these galaxies harbor SMBHs at high redshifts. We illustrate this possibility using cosmological Monte Carlo simulations of the merger history of dark matter halos and their associated SMBHs. In our most extreme models, in order to populate nearly every bright galaxy with a SMBH at z=0, only a few percent of the halos with virial temperatures above 10^4 K are required to harbor a SMBH at high redshift. This possibility must be included in studies of the luminosity function and the clustering properties of quasars. We predict the number of SMBH merger events that are detectable by the gravitational wave experiment LISA, as a function of redshift, out to z=5. Although the event rates can be significantly reduced in scenarios with rare SMBHs, a minimum of \~10 detectable merger events per year is predicted if SMBH binaries coalesce efficiently. The observed distribution of events with redshift could yield valuable information on the SMBH formation process. If SMBH binaries do not coalesce, we find that at least several SMBH slingshot ejections probably occurred from z=5 to the present in each galaxy more massive than ~10^11 Msun at z=0. Although our results are sensitive to the minimum cooling mass assumed for the formation of SMBHs, we expect the qualitative predictions of our models to be robust.

We reassess constraints on the cosmological baryon density from observations of the mean decrement and power spectrum of the Lyman-alpha forest, taking into account uncertainties in all free parameters in the simplest gravitational instability model. The uncertainty is dominated by that of the photoionizing background, but incomplete knowledge of the thermal state of the intergalactic medium also contributes significantly to the error-budget. While current estimates of the baryon fraction from the forest do prefer values that are somewhat higher than the big bang nucleosynthesis value of \Omega_b h^2 = 0.02 \pm 0.001, the discrepancy is at best about 3 \sigma. For instance, assuming the highest estimate of the ionizing background, as indicated by recent measurements of a large escape fraction from Lyman-break galaxies by Steidel, Pettini & Adelberger, we find \Omega_b h^2 = 0.045 \pm 0.008. A recent measurement of the ionizing background from the proximity effect by Scott et al., on the other hand, implies \Omega_b h^2 = 0.03 \pm 0.01. We provide an expression from which future likelihoods for \Omega_b h^2 can be derived as measurements of the ionizing background improve -- consistency among constraints from the forest, nucleosynthesis and the microwave background will provide a powerful test of the gravitational instability model. We also develop a formalism which treats lower bounds on the baryon density in a statistical manner, which is appropriate if only a lower bound on the ionizing background is known. Finally, we discuss the implications of the escape fraction measurement for the age, structure and stellar content of Lyman-break galaxies.

The early stage in the formation of a galaxy inevitably involves a spatially extended distribution of infalling, cold gas. If a central luminous quasar turned on during this phase, it would result in significant extended Lyman alpha emission, possibly accompanied by other lines. For halos condensing at redshifts between 3 < z < 8 and having virial temperatures between 2 x 10^5 K and 2 x 10^6 K, this emission results in a ``fuzz'' of characteristic angular diameter of a few arcseconds, and surface brightness between 10^-18 and 10^-16 erg/s/cm^2/asec^2. The fuzz around bright, high redshift quasars could be detected in deep narrow band imaging with current telescopes, providing a direct constraint on galaxy formation models. The absence of detectable fuzz might suggest that most that most of the protogalaxy's gas settles to a self-gravitating disk before a quasar turns on. However, continued gas infall from large radii, or an on-going merger spreading cold gas over a large solid angle, during the luminous quasar phase could also result in extended Ly$\alpha$ emission, and can be constrained by deep narrow band imaging.

Warm Dark Matter (WDM) models have recently been resurrected to resolve apparent conflicts of Cold Dark Matter (DM) models with observations. Endowing the DM particles with non-negligible velocities causes free-streaming, which suppresses the primordial power spectrum on small scales. The choice of a root-mean-square velocity dispersion v(rms) = 0.05 km/s at redshift z=0 (corresponding to a particle mass of 1 keV if the WDM particles are fermions decoupling while relativistic) helps alleviate most, but probably not all, of the small-scale problems faced by CDM. An important side-effect of the particle velocities is the severe decrease in the number of collapsed halos at high redshift. This is caused both by the loss of small-scale power, and by the delay in the collapse of the smallest individual halos (with masses near the effective Jeans mass of the DM). The presence of early halos is required in order (1) to host either early quasars or galaxies that can reionize the universe by redshift z=5.8, and (2) to allow the growth of the supermassive black hole believed to power the recently discovered quasar SDSS 1044-1215 at this redshift. We quantify these constraints using a modified Press-Schechter formalism, and find v(rms) < 0.04 km/s (or m_X > 1 keV). If future observations uncover massive black holes at z > 10, or reveal that reionization occurred at z > 10, this could conclusively rule out WDM models as the solution to the small-scale crisis of the CDM paradigm.

Over the coming decade, the observational samples available for studies of cluster abundance evolution will increase from tens to hundreds, or possibly to thousands, of clusters. Here we assess the power of future surveys to determine cosmological parameters. We quantify the statistical differences among cosmologies, including the effects of the cosmic equation of state parameter w, in mock cluster catalogs simulating a 12 sq. deg Sunyaev-Zeldovich Effect survey and a deep 10^4 sq. deg X-ray survey. The constraints from clusters are complementary to those from studies of high-redshift Supernovae (SNe), CMB anisotropies, or counts of high-redshift galaxies. Our results indicate that a statistical uncertainty of a few percent on both Omega_m and w can be reached when cluster surveys are used in combination with any of these other datasets.

The lifetime t_Q of the luminous phase of quasars is constrained by current observations to be between 10^6 and 10^8 years, but is otherwise unkown. We model the quasar luminosity function in detail in the optical and X-ray bands using the Press-Schechter formalism, and show that the expected clustering of quasars depends strongly on their assumed lifetime. We quantify this dependence, and find that existing measurements of the correlation length of quasars are consistent with the range 10^6 < t_Q < 10^8 years. We then show that future measurements of the power spectrum of quasars out to z=3, from the 2dF or Sloan Digital Sky Survey, can significantly improve this constraint, and in principle allow a precise determination of t_Q. We estimate the systematic errors introduced by uncertainties in the modeling of the quasar-halo relationship, as well as by the possible existence of obscured quasars.

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 study the expected redshift evolution of galaxy cluster abundance between 0 < z < 3 in different cosmologies, including the effects of the cosmic equation of state parameter w=p/rho. Using the halo mass function obtained in recent large scale numerical simulations, we model the expected cluster yields in a 12 deg^2 Sunyaev-Zeldovich Effect (SZE) survey and a deep 10^4 deg^2 X-ray survey over a wide range of cosmological parameters. We quantify the statistical differences among cosmologies using both the total number and redshift distribution of clusters. Provided that the local cluster abundance is known to a few percent accuracy, we find only mild degeneracies between w and either Omega_m or h. As a result, both surveys will provide improved constraints on Omega_m and w. The Omega_m-w degeneracy from both surveys is complementary to those found either in studies of CMB anisotropies or of high-redshift Supernovae (SNe). As a result, combining these surveys together with either CMB or SNe studies can reduce the statistical uncertainty on both w and Omega_m to levels below what could be obtained by combining only the latter two data sets. Our results indicate a formal statistical uncertainty of about 3% (68% confidence) on both Omega_m and w when the SZE survey is combined with either the CMB or SN data; the large number of clusters in the X-ray survey further suppresses the degeneracy between w and both Omega_m and h. Systematics and internal evolution of cluster structure at the present pose uncertainties above these levels. We briefly discuss and quantify the relevant systematic errors. By focusing on clusters with measured temperatures in the X-ray survey, we reduce our sensitivity to systematics such as non-standard evolution of internal cluster structure.

The recent discoveries of luminous quasars at high redshifts imply that black holes more massive than a few billion solar masses have been assembled already when the universe was less than a billion years old. We show that the existence of these black holes is not surprising in popular hierarchical models of structure formation. For example, the black hole needed to power the quasar SDSS 1044-0125 at redshift z=5.8 can arise naturally from the growth of stellar mass seeds forming at z>10, when typical values are assumed for the radiative accretion efficiency and the bolometric accretion luminosity in Eddington units. Nevertheless, SDSS 1044-0125 yields a non-trivial constraint on a combination of these parameters. Extrapolating our model to future surveys, we derive the highest plausible redshift for quasars which are not lensed or beamed, as a function of their apparent magnitude. We find that at a limiting magnitude of K(AB)=20, quasar surveys can yield strong constraints on the growth of supermassive black holes out to z=10.

The large-scale structure of high-redshift galaxies produces correlated anisotropy in the far-infrared background (FIRB). In regions of the sky where the thermal emission from Galactic dust is well below average, these high-redshift correlations may be the most significant source of angular fluctuation power over a wide range of angular scales, from about 7' to about 3 degrees, and frequencies, from about 400 to about 1000 GHz. The strength of this signal should allow detailed studies of the statistics of the FIRB fluctuations, including the shape of the angular power spectrum at a given frequency and the degree of coherence between FIRB maps at different frequencies. The FIRB correlations depend upon and hence constrain the redshift-dependent spectral energy distributions, number counts, and clustering bias of the galaxies and active nuclei that contribute to the background. We quantify the accuracy to which Planck and a newly proposed balloon-borne mission EDGE could constrain models of the high-redshift universe through the measurement of FIRB fluctuations. We conclude that the average bias of high-redshift galaxies could be measured to an accuracy of less than approximately 1% or, for example, separated into 4 redshift bins with about 10% accuracy.

Ionizing sources embedded in the neutral intergalactic medium (IGM) before cosmological reionization generate discrete HII regions. We show that a sufficiently bright quasar (for example, one tenth as luminous as that recently discovered by Fan et al.) can ionize a large volume, allowing the transmission of a substantial fraction of the flux of its intrinsic Lyman alpha emission line on both sides of the Lyman alpha wavelength. The observed line profile is richly informative. We show that a sufficiently accurate, high spectral resolution ($R\approx 10^4$) measurement of the line profile of a bright quasar is feasible using the Next Generation Space Telescope (NGST) as well as large ground-based telescopes. Such a measurement has two potentially important applications. First, the red side of the Lyman alpha emission line provides a way to measure the quasar lifetime. Second, the blue side provides a direct measure of the density fluctuations of the intergalactic medium at the quasar redshift. The estimate of the absorption of the red side is limited by the accuracy to which the intrinsic profile of the Lyman alpha emission line is known. The blue side, however, does not sensitively depend on the intrinsic profile, because the former is much narrower than the latter.

We present observations of SDSSp J104433.04--012502.2, a luminous quasar at z=5.80 discovered from Sloan Digital Sky Survey (SDSS) multicolor imaging data. This object was selected as an i'-band dropout object, with i*=21.8 +/- 0.2, z*=19.2 +/- 0.1. It has an absolute magnitude M1450 = -27.2 (H_0 =50 km/s/Mpc, q0 = 0.5). The spectrum shows a strong and broad Ly alpha emission line, strong Ly alpha forest absorption lines with a mean continuum decrement D_A = 0.91, and a Lyman Limit System at z=5.72. The spectrum also shows strong OI and SiIV emission lines similar to those of quasars at z<= 5, suggesting that these metals were produced at redshift beyond six. The lack of a Gunn-Peterson trough in the spectrum indicates that the universe is already highly ionized at z ~ 5.8. Using a high-resolution spectrum in the Ly alpha forest region, we place a conservative upper limit of the optical depth due to the Gunn-Peterson effect of tau < 0.5 in regions of minimum absorption. The Ly alpha forest absorption in this object is much stronger than that in quasars at z<= 5. The object is unresolved in a deep image with excellent seeing, implying that it is unlensed. The black hole mass of this quasar is ~3 x 10^9 M_solar if we assume that it is radiating at the Eddington luminosity and no lensing amplification, implying that it resides in a very massive dark matter halo. The discovery of one quasar at M_1450 < -27 in a survey area of 600 deg^2 is consistent with an extrapolation of the observed luminosity function at lower redshift. The abundance and evolution of such quasars can provide sensitive tests of models of quasar and galaxy formation.

The properties of nearby galaxy clusters limit the range of cosmological parameters consistent with our universe. We describe the limits which arise from studies of the intracluster medium (ICM) mass fraction fICM and consideration of the possible sources of systematic error: Omega_M<0.44h_{50}^{-1/2} at 95% confidence. We emphasize that independent of Type Ia supernovae (SNe Ia) observations, this cluster study, taken together with published cosmic microwave background (CMB) anisotropy studies, indicates a non-zero quintessence or dark energy component Omega_Q>0. We then discuss future galaxy cluster surveys which will probe the abundance of galaxy clusters to intermediate and high redshift. We investigate the sensitivity of these surveys to the cosmological density parameter Omega_M and the equation of state parameter w of any quintessence component. In particular, we show that cluster survey constraints from a proposed large solid angle X-ray survey are comparable in precision and complementary in nature to constraints expected from future CMB anisotropy and SNe Ia studies.

The baryons inside high-redshift halos with virial temperatures above T=10^4 K cool radiatively as they condense inside dark matter potential wells. We show that the release of the gravitational binding energy, over the halo assembly time-scale, results in a significant and detectable Ly alpha flux. At the limiting line flux of approximately 10^-19 erg/s/cm^2/asec^2 of the Next Generation Space Telescope, several sufficiently massive halos, with velocity dispersions above sigma=120 km/s, would be visible per 4'x4' field. The halos would have characteristic angular sizes of approximately 10 asec, would be detectable in a broad-band survey out to redshifts z=6-8, and would provide a direct probe of galaxies in the process of forming. They may be accompanied by HeII Ly alpha emission at the 10 percent level, but remain undetectable at other wavelengths. Our predictions are in good agreement with the recent finding of two Ly alpha ``blobs'' at z=3.1 by Steidel et al. (1999).

We review the relevance of H2 molecules for structure formation in cosmology. Molecules are important at high redshifts, when the first collapsed structures appear with typical temperatures of a few hundred Kelvin. In these chemically pristine clouds, radiative cooling is dominated H2 molecules. As a result, H2 ``astro-chemistry'' is likely to determine the epoch when the first astrophysical objects appear. We summarize results of recent three-dimensional simulations. A discussion of the effects of feedback, and implications for the reionization of the universe is also given.

Empirical studies of the first generation of stars and quasars will likely become feasible within the next decade in several different wavelength bands. Microwave anisotropy experiments, such as MAP or Planck, will set constraints on the ionization history of the intergalactic medium due to these sources. In the infrared, the Next Generation Space Telescope (NGST) will be able to directly detect sub-galactic objects at redshifts z>10. In the optical, data from the Hubble Deep Field already places a constraint on the abundance of high-redshift quasars. However, the epoch of the first quasars might be first probed in X-ray bands, by instruments such as the Chandra X-ray Observatory (CXO) and the X-ray Multi-mirror Mission (XMM). In a 500 Ksec integration, CXO reaches a sensitivity of 2x10-16 erg/s cm^2. Based on simple hierarchical CDM models, we find that at this flux threshold approx. 100 quasars might be detectable from redshifts z>5, and approx. 1 quasar at z=10, in each 17x17 arcmin field. Measurement of the power spectrum of the unresolved soft X-ray background will further constrain models of faint, high-redshift quasars.

Recent observations of Type Ia SNe at redshifts 0 < z <1 reveal a progressive dimming which has been interpreted as evidence for a cosmological constant of Omega_Lambda ~ 0.7. An alternative explanation of the SN results is an open universe with Omega_Lambda=0 and the presence of > 0.1 micron dust grains with a mass density of Omega_dust ~ (few) * 10^{-5} in the intergalactic (IG) medium. The same dust that dims the SNe absorbs the cosmic UV/optical background radiation around ~ 1 micron, and re-emits it at far infrared (FIR) wavelengths. Here we compare the FIR emission from IG dust with observations of the cosmic microwave (CMB) and cosmic far infrared backgrounds (FIRB) by the DIRBE/FIRAS instruments. We find that the emission would not lead to measurable distortion to the CMB, but would represent a substantial fraction (> 75%) of the measured value of the FIRB in the 300-1000 micron range. This contribution would be marginally consistent with the present unresolved fraction of the observed FIRB in an open universe. However, we find that IG dust probably could not reconcile the standard Omega=1 CDM model with the SN observations, even if the necessary quantity of dust existed. Future observations able to reliably resolve the FIRB to a flux limit of ~0.5 mJy, along with a more precise measure of the coarse-grained FIRB, will provide a definitive test of the IG dust hypothesis in all cosmologies.

We consider a class of models for the redshift evolution (between $0\lsim z \lsim 4$) of the observed optical and X-ray quasar luminosity functions (LFs), with the following assumptions: (i) the mass-function of dark matter halos follows the Press-Schechter theory, (ii) the black hole (BH) mass scales linearly with the halo mass, (iii) quasars have a constant universal lifetime, and (iv) a thin accretion disk provides the optical luminosity of quasars, while the X-ray/optical flux ratio is calibrated from a sample of observed quasars. The mass accretion rate $\dot{M}$ onto quasar BHs is a free parameter of the models, that we constrain using the observed LFs. The accretion rate $\dot M$ inferred from either the optical or X-ray data under these assumptions generally decreases as a function of cosmic time from $z \simeq 4$ to $z \simeq 0$. We find that a comparable accretion rate is inferred from the X-ray and optical LF only if the X-ray/optical flux ratio decreases with BH mass. Near $z\simeq 0$, $\dot M$ drops to substantially sub-Eddington values at which advection-dominated accretion flows (ADAFs) exist. Such a decline of $\dot M$, possibly followed by a transition to radiatively inefficient ADAFs, could explain both the absence of bright quasars in the local universe and the faintness of accreting BHs at the centers of nearby galaxies. We argue that a decline of the accretion rate of the quasar population is indeed expected in cosmological structure formation models.

We compute the expected angular power spectrum of the cosmic Far Infrared Background (FIRB). We find that the signal due to source correlations dominates the shot--noise for $\ell \la 1000$ and results in anisotropies with rms amplitudes $(\sqrt{\ell(\ell+1)C_\ell/2\pi})$ between 5% and 10% of the mean for $l \ga 150$. The angular power spectrum depends on several unknown quantities, such as the UV flux density evolution, optical properties of the dust, biasing of the sources of the FIRB, and cosmological parameters. However, when we require our models to reproduce the observed DC level of the FIRB, we find that the anisotropy is at least a few percent in all cases. This anisotropy is detectable with proposed instruments, and its measurement will provide strong constraints on models of galaxy evolution and large-scale structure at redshifts up to at least $z \sim5$.

It is currently unknown whether the Universe was reionized by quasars or stars at z>5. We point out that quasars can be best distinguished from stellar systems by their X-ray emission. Based on a simple hierarchical CDM model, we predict the number counts and X-ray fluxes of quasars at high redshifts. The model is consistent with available data on the luminosity function of high-redshift quasars in the optical and soft X-ray bands. The cumulative contribution of faint, undetected quasars in our model is consistent with the unresolved fraction of the X-ray background. We find that the Chandra X-ray Observatory might detect approximately 100 quasars from redshifts z>5 per its 17'x17' field of view at the flux threshold of 2x10^-16 erg/s/cm^2. The redshifts of these faint point-sources could be identified by follow-up infrared observations from the ground or with the Next Generation Space Telescope.

The bulk of the hydrogen in the universe transformed from neutral to ionized somewhere in the redshift interval 5 < z < 40, most likely due to ionizing photons produced by an early generation of stars or mini-quasars. The resulting free electrons, interacting with the CMB photons via Thomson scattering, are a mixed blessing, providing both a probe of the epoch of the first stars and a contaminant to the pristine primary anisotropy. Here we review our current knowledge of reionization with emphasis on inhomogeneities and describe the possible connections to CMB anisotropy.

In hierarchical models of structure formation, an early cosmic UV background (UVB) is produced by the small (T_vir < 10^4 K) halos that collapse before reionization. The UVB at energies below 13.6eV suppresses the formation of stars or black holes inside small halos, by photo-dissociating their only cooling agent, molecular H2. We self-consistently compute the buildup of the early UVB in Press-Schechter models, coupled with H2 photo-dissociation both in the intergalactic medium (IGM), and inside virialized halos. We find that the intergalactic H2 has a negligible effect on the UVB, both because its initial optical depth is small (tau<0.1), and because it is photo-dissociated at an early stage. If the UV sources in the first collapsed halos are stars, then their UV flux suppresses further star-formation inside small halos. This results in a pause in the buildup of the UVB, and reionization is delayed until larger halos (T_vir> 10^4 K) collapse. If the small halos host mini-quasars with hard spectra extending to approximately 1 keV, then their X-rays balance the effects of the UVB, the negative feedback does not occur, and reionization can be caused by the small halos.

We present models for dusty high-redshift Lyman alpha emitting galaxies by combining the Press-Schechter formalism with a treatment of inhomogeneous dust distribution inside galaxies. These models reproduce the surface density of emitters inferred from recent observations, and also agree with previous non-detections. Although a detailed determination of the individual model parameters is precluded by uncertainties, we find that (i) the dust content of primordial galaxies builds up in no more than 5x10^8 years, (ii) the galactic HII regions are inhomogeneous with a cloud covering factor of order unity, and (iii) the overall star formation efficiency is at least 5 percent. Future observations should be able to detect Lyman alpha galaxies upto redshifts of about z=8. If the universe is reionized at z(reion)<8, the corresponding decline in the number of Lyman alpha emitters at z>z(reion) could prove to be a useful probe of the reionization epoch.

Empirical studies of the first generation of stars and quasars in the Universe will likely become feasible over the next decade. The Next Generation Space Telescope will provide direct imaging and photometry of sub-galactic objects at redshifts above z=10, while microwave anisotropy experiments, such as MAP or Planck, will set constraints on the ionization history of the intergalactic medium due to these sources. We describe the expected signals that will be detectable with these future instruments.

We review observable signatures of the first generation of stars and low-luminosity quasars, including the metal enrichment, radiation background, and dust opacity/emission that they produce. We compute the formation history of collapsed baryonic halos based on an extension of the Press-Schechter formalism, incorporating the effects of pressure and H2-dissociation, and calibrate this history to observational data at redshifts z<5. To be consistent with the lack of faint point-sources in the Hubble Deep Field, we introduce a lower limit of 75 km/s for the circular velocities of halos harboring central black holes. Based on our models, we predict that future satellite instruments such as MAP, Planck, and NGST will be able to detect the signatures small star clusters and low-luminosity quasars at redshifts exceeding z=10.

Models are presented for the Lyman alpha emission of dusty high-redshift galaxies by combining the Press-Schechter formalism with a treatment of the inhomogeneous dust distribution inside galaxies. It is found that the amount of Lyman alpha radiation escaping from the galaxies strongly depends on the time over which the dust is produced through stellar activity, and on the ambient inhomogeneity of the HII regions that surround the ionizing OB stars. Good agreement is found with recent observations, as well as previous non-detections. Our models indicate that the dust content builds up in no more than approximately 5x10^8 yr, the galactic HII regions are inhomogeneous with a cloud covering factor of order unity, and the overall star formation efficiency is at least about 5%. It is predicted that future observations can detect these Lyman alpha galaxies upto redshifts of about 8.

The redshift at which the universe was reionized is currently unknown. We examine the optimal strategy for extracting this redshift, z_{reion}, from the spectra of early sources. For a source located at a redshift z_s beyond but close to reionization, 1<[(1+z_s)/(1+z_{reion})]<(32/27), the Gunn-Peterson trough splits into disjoint Lyman-alpha, beta, and possibly higher Lyman series troughs, with some transmitted flux in between these troughs. We show that although the transmitted flux is suppressed considerably by the dense Lyman-alpha forest after reionization, it is still detectable for sufficiently bright sources and can be used to infer the reionization redshift. The Next Generation Space Telescope will reach the spectroscopic sensitivity required for the detection of such sources.

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.

We study the observational signatures of a potential population of low-luminosity quasars at high redshifts in a LambdaCDM cosmology. We derive the evolution of the quasar luminosity function at fainter luminosities and higher redshifts than currently detected, based on three assumptions: (1) the formation of dark-matter halos follows the Press-Schechter theory, (2) the ratio of central black hole mass to halo mass is the same for all halos, and (3) the light-curve of quasars, in Eddington units, is universal. We show that a universal light-curve provides an excellent fit to the observed quasar luminosity function at redshifts 2.6 < z < 4.5. By extrapolating the evolution of this luminosity function to higher redshifts (4.5 < z < 20), we find that the associated early population of low-luminosity quasars reionizes the universe around redshift z=12. The reprocessing of the UV light of these quasars by dust from early type II supernovae, distorts the microwave background spectrum by a Compton y-parameter, y = 1e-5, comparable to the upper limit set by COBE. The Next Generation Space Telescope could detect tens of quasars from redshifts z > 10 per square arcminute, with its proposed 1nJy sensitivity at 1-3.5 microns. Absorption spectra of several such quasars would reveal the reionization history of the universe.

The high ionization level and universal metallicity (1% solar) of the intergalactic gas at redshifts z<5 implies that nonlinear structure had started to form in the universe at earlier times than we currently probe. In Cold Dark Matter (CDM) cosmologies, the first generation of baryonic objects emerges at redshifts z=10-50. Here we examine the observable consequences of the possibility that an early generation of stars reionized the universe and resulted in the observed metallicity of the Lyman-alpha forest. Forthcoming microwave anisotropy experiments will be sensitive to the damping of anisotropies caused by scattering off free electrons from the reionization epoch. For a large range of CDM models with a Scalo stellar mass function, we find that reionization occurs at a redshift z>10 and damps the amplitude of anisotropies on angular scales <10 degrees by a detectable amount, of order 10-25%. However, reionization is substantially delayed if the initial stellar mass function transformed most of the baryons into low mass stars. In this case, the mass fraction of pre-galactic stars could be constrained from the statistics of microlensing events in galactic halos or along lines of sight to quasars. Deep infrared imaging with future space telescopes (such as SIRTF or the Next Generation Space Telescope) will be able to detect bright star clusters at z>5. The cumulative Bremsstrahlung emission from these star clusters yields a measurable distortion to the spectrum of the microwave background.

We calculate the observable signatures of the first generation of stars at high redshift (5<z<100). To determine the cosmic star-formation history, we use an extension of the Press-Schechter formalism for Cold Dark Matter (CDM) cosmologies that incorporates gas pressure. We calibrate the fraction of gas converted into stars to be 6% so as to reproduce the 1% solar C/H ratio observed in the intergalactic medium (IGM) at z=3. With this star-formation efficiency, we find that NGST would be able to image more than 10^4 star clusters from high redshifts (z>10) within its 4'x4' field of view. If stars occupy a region comparable to the virial radius of the cluster, then 1% of these clusters could be resolved. We calculate the expected number-flux relation and angular size distribution for these early star clusters. We also describe the reionization of the IGM due to the first generation of stars, and the consequent damping of the CMB anisotropies on small angular scales. This damping could be detected below 10 degree angular scales by MAP and PLANCK.

We quantify the consequences of intergalactic dust produced by the first Type II supernovae in the universe. The fraction of gas converted into stars is calibrated based on the observed C/H ratio in the intergalactic medium at z=3, assuming a Scalo mass function for the stars. The associated dust absorbs starlight energy and emits it at longer wavelengths. For a uniform mix of metals and dust with the intergalactic gas, we find that the dust distorts the microwave background spectrum by a y-parameter in the range (0.02-2)x10^{-5}(M_d/0.3M_sun), where M_d is the average mass of dust produced per supernova. The opacity of intergalactic dust to infrared sources at redshifts z>10 is significant, tau=(0.1-1)x(M_d/0.3M_sun), and could be detected with the Next Generation Space Telescope. Although dust suppresses the Ly-alpha emission from early sources, the redshifts of star clusters at z=10-35 can be easily inferred from the Lyman-limit break in their infrared spectrum between 1-3.5 micron.

We use a spherical hydrodynamics code to show that in cold dark matter cosmologies, the first stars form at z=50 through the direct collapse of gas in low-mass systems (approx 10^4 solar masses). Photons from the first stars easily photodissociate the molecular hydrogen throughout the universe and so molecular cooling does not affect the subsequent fragmentation of gas clouds into stars. We examine observable signatures of the pre-galactic population of stars. These include the detected metallicity and photo-ionization of the intergalactic medium, and the soon to be detected damping of microwave anisotropies on small angular scales (< 10 degrees). The Next Generation Space Telescope will be able to directly image the pre-galactic star clusters, while the DIMES experiment could detect their Bremsstrahlung emission.

We investigate the ability of primordial gas clouds to retain molecular hydrogen (H_2) during the initial phase of the reionization epoch. We find that before the Stromgren spheres of the individual ionizing sources overlap, the UV background below the ionization threshold is able to penetrate large clouds and suppress their H_2 abundance. The consequent lack of H_2 cooling could prevent the collapse and fragmentation of clouds with virial temperatures T_vir < 10^4 K (or masses 10^8 Msun [(1+z_vir)/10]^{-3/2}). This negative feedback on structure-formation arises from the very first ionizing sources, and precedes the feedback due to the photoionization heating.

We investigate the formation of molecular hydrogen (H_2) in a primordial H+He gas cloud irradiated by a power-law UV flux. We find that at high densities (>1 cm^{-3}) and low temperatures (<10^4 K), the background radiation enhances the formation of H_2 and results in molecular cooling dominating over photoionization heating. This process could accelerate the collapse and fragmentation of dense objects at high redshifts.

We investigate the early formation of bound objects with masses comparable to the cosmological Jeans mass (10^5 solar masses). We follow the growth of isolated spherically symmetric density peaks starting from the linear perturbative regime. The initial parameters correspond to density peaks of various widths and heights in a Cold Dark Matter cosmology. We use a one-dimensional spherical Lagrangian hydrodynamics code to follow the dynamical, thermal, and non-equilibrium chemical evolution of the gas. The system includes a collisionless dark matter component and a baryonic component composed of the nine species H, H^-, H^+, He, He^+, He^{++}, H_2, H_2^+, and e^-. All relevant chemical reactions between these species and their cooling mechanisms are included in the calculations. We find that radiative cooling by H_2 affects the collapse dynamics of the gas only after it has already virialized and become part of the bound object. Therefore, radiative cooling is unlikely to have triggered the initial collapse of perturbations at redshifts z>10. Nevertheless, objects with baryonic masses well below the linear-theory Jeans mass (<10^3 solar masses) collapse due to shell crossing by the dark matter. Such objects could be the progenitors of a primordial population of high-mass stars in the intergalactic medium.