Using the largest mosaic of Hubble Space Telescope images around a galaxy cluster, we map the distribution of dark matter throughout a $\sim$$6\times6$ Mpc$^2$ area centred on the cluster MS 0451-03 ($z=0.54$, $M_{200}=1.65\times10^{15} \rm{M}_\odot$). Our joint strong- and weak-lensing analysis shows three possible filaments extending from the cluster, encompassing six group-scale substructures. The dark-matter distribution in the cluster core is elongated, consists of two distinct components, and is characterized by a concentration parameter of $c_{200}=3.79\pm0.36$. By contrast, XMM-Newton observations show the gas distribution to be more spherical, with excess entropy near the core, and a lower concentration of $c_{200}=2.35^{+0.89}_{-0.70}$ (assuming hydrostatic equilibrium). Such a configuration is predicted in simulations of major mergers 2-7Gyr after the first core passage, when the two dark-matter halos approach second turnaround, and before their gas has relaxed. This post-merger scenario finds further support in optical spectroscopy of the cluster's member galaxies, which shows that star formation was abruptly quenched 5 Gyr ago. MS 0451-03 will be an ideal target for future studies of the growth of structure along filaments, star-formation processes after a major merger, and the late-stage evolution of cluster collisions.

The dense environment of a galaxy cluster can radically transform the content of in-falling galaxies. Recent observations have found a significant population of active galactic nuclei (AGN) within "jellyfish galaxies," galaxies with trailing tails of gas and stars that indicate significant ram pressure stripping. The relationship between AGN and ram pressure stripping is not well understood. In this letter, we investigate the connection between AGN activity and ram pressure in a fully cosmological setting for the first time using the RomulusC simulation, one of the highest resolution simulations of a galaxy cluster to date. We find unambiguous morphological evidence for ram pressure stripping. For lower mass galaxies (with stellar masses < 10^9.5 solar masses) both star formation and black hole accretion are suppressed by ram pressure before they reach pericenter, whereas for more massive galaxies accretion onto the black hole is enhanced during pericentric passage. Our analysis also indicates that as long as the galaxy retains gas, AGN with higher Eddington ratios are more likely to be the found in galaxies experiencing higher ram pressure. We conclude that prior to quenching star formation, ram pressure triggers enhanced accretion onto the black hole, which then produces heating and outflows due to AGN feedback. AGN feedback may in turn serve to aid in the quenching of star formation in tandem with ram pressure.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

We present a new galaxy cluster lens modeling approach, hybrid-Lenstool, that is implemented in the publicly available modeling software Lenstool. hybrid-Lenstool combines a parametric approach to model the core of the cluster, and a non-parametric (free-form) approach to model the outskirts. hybrid-Lenstool optimizes both strong- and weak-lensing constraints simultaneously (Joint-Fit), providing a self-consistent reconstruction of the cluster mass distribution on all scales. In order to demonstrate the capabilities of the new algorithm, we tested it on a simulated cluster. hybrid-Lenstool yields more accurate reconstructed mass distributions than the former Sequential-Fit approach where the parametric and the non-parametric models are optimized successively. Indeed, we show with the simulated cluster that the mass density profile reconstructed with a Sequential-Fit deviates form the input by $2-3\sigma$ at all scales while the Joint-Fit gives a profile that is within $1-1.5\sigma$ of the true value. This gain in accuracy is consequential for recovering mass distributions exploiting cluster lensing and therefore for all applications of clusters as cosmological probes. Finally we found that the Joint-Fit approach yields shallower slope of the inner density profile than the Sequential-Fit approach, thus revealing possible biases in previous lensing studies.

Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field.

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

In this white paper, we discuss the prospects for characterizing and identifying dark matter using gravitational waves, covering a wide range of dark matter candidate types and signals. We argue that present and upcoming gravitational wave probes offer unprecedented opportunities for unraveling the nature of dark matter and we identify the most urgent challenges and open problems with the aim of encouraging a strong community effort at the interface between these two exciting fields of research.

There exist hitherto unexplained fluctuations in the Cosmic Infrared Background (CIB) on arcminute scales and larger. These have been shown to cross-correlate with the Cosmic X-ray Background (CXB), leading several authors to attribute the excess to a high-redshift growing black hole population. In order to investigate potential sources that could explain this excess, in this paper, we develop a new framework to compute the power spectrum of undetected sources that do not have constant flux as a function of halo mass. In this formulation, we combine a semi-analytic model for black hole growth and their simulated spectra from hydrodynamical simulations. Revisiting the possible contribution of a high-redshift black hole population, we find that too much black hole growth is required at early epochs for z>6 accretion to explain these fluctuations. Examining a population of accreting black holes at more moderate redshifts, z\sim 2-3, we find that such models produce a poor fit to the observed fluctuations while simultaneously overproducing the local black hole mass density. Additionally, we rule out the hypothesis of a missing Galactic foreground of warm dust that produces coherent fluctuations in the X-ray via reflection of Galactic X-ray binary emission. Although we firmly rule out accreting massive black holes as the source of these missing fluctuations, additional studies will be required to determine their origin.

We study the link between supermassive black hole growth and the stellar mass assembly of their host galaxies in the state-of-the-art Romulus suite of simulations. The cosmological simulations Romulus25 and RomulusC employ innovative recipes for the seeding, accretion, and dynamics of black holes in the field and cluster environments respectively. We find that the black hole accretion rate traces the star formation rate among star-forming galaxies. This result holds for stellar masses between 10^8 and 10^12 solar masses, with a very weak dependence on host halo mass or redshift. The inferred relation between accretion rate and star formation rate does not appear to depend on environment, as no difference is seen in the cluster/proto-cluster volume compared to the field. A model including the star formation rate, the black hole-to-stellar mass ratio, and the cold gas fraction can explain about 70 per cent of all variations in the black hole accretion rate among star forming galaxies. Finally, bearing in mind the limited volume and resolution of these cosmological simulations, we find no evidence for a connection between black hole growth and galaxy mergers, on any timescale and at any redshift. Black holes and their galaxies assemble in tandem in these simulations, regardless of the larger-scale intergalactic environment, suggesting that black hole growth simply follows star formation on galactic scales.

The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps.

The Laser Interferometer Space Antenna (LISA) will open three decades of gravitational wave (GW) spectrum between 0.1 and 100 mHz, the mHz band. This band is expected to be the richest part of the GW spectrum, in types of sources, numbers of sources, signal-to-noise ratios and discovery potential. When LISA opens the low-frequency window of the gravitational wave spectrum, around 2034, the surge of gravitational-wave astronomy will strongly compel a subsequent mission to further explore the frequency bands of the GW spectrum that can only be accessed from space. The 2020s is the time to start developing technology and studying mission concepts for a large-scale mission to be launched in the 2040s. The mission concept would then be proposed to Astro2030. Only space based missions can access the GW spectrum between 10 nHz and 1 Hz because of the Earths seismic noise. This white paper surveys the science in this band and mission concepts that could accomplish that science. The proposed small scale activity is a technology development program that would support a range of concepts and a mission concept study to choose a specific mission concept for Astro2030. In this white paper, we will refer to a generic GW mission beyond LISA as bLISA.

TeamIndus' lunar logistics vision includes multiple lunar missions to meet requirements of science, commercial and efforts towards global exploration. The first mission is slated for launch in 2020. The prime objective is to demonstrate autonomous precision lunar landing, and Surface Exploration Rover to collect data on the vicinity of the landing site. TeamIndus has developed various technologies towards lowering the access barrier to the lunar surface. This paper shall provide an overview of design of lander GNC system. The design of the GNC system has been described after concluding studies on sensor and actuator configurations. Frugal design approach is followed in the selection of GNC hardware. The paper describes the constraints for the orbital maneuvers and the lunar descent strategy. Various aspects of the GNC design of autonomous lunar descent maneuver: timeline of events, guidance, inertial and optical terrain-relative navigation schemes are described. The GNC software description focuses on system architecture, modes of operation, and core elements of the GNC software. The GNC algorithms have been tested using Monte-Carlo simulations and Processor-in-Loop runs. The paper concludes with a summary of key risk-mitigation studies for soft landing.

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.

A space-based interferometer operating in the previously unexplored mHz gravitational band has tremendous discovery potential. If history is any guide, the opening of a new spectral band will lead to the discovery of entirely new sources and phenomena. The mHz band is ideally suited to exploring beyond standard model processes in the early universe, and with the sensitivities that can be reached with current technologies, the discovery space for exotic astrophysical systems is vast.

Observational evidence has been mounting for the existence of intermediate mass black holes (IMBHs, 10^2-10^5 Msun), but observing them at all, much less constraining their masses, is very challenging. In one theorized formation channel, IMBHs are the seeds for supermassive black holes in the early universe. As a result, IMBHs are predicted to exist in the local universe in dwarf galaxies, as well as wandering in more massive galaxy halos. However, these environments are not conducive to the accretion events or dynamical signatures that allow us to detect IMBHs. The Laser Interferometer Space Antenna (LISA) will demystify IMBHs by detecting the mergers of these objects out to extremely high redshifts, while measuring their masses with extremely high precision. These observations of merging IMBHs will allow us to constrain the formation mechanism and subsequent evolution of massive black holes, from the 'dark ages' to the present day, and reveal the role that IMBHs play in hierarchical galaxy evolution.

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.

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.

Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence.

The nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

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.

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

There are two big questions cosmologists would like to answer -- How does the Universe work, and what are its origin and destiny? A long wavelength gravitational wave detector -- with million km interferometer arms, achievable only from space -- gives a unique opportunity to address both of these questions. A sensitive, mHz frequency observatory could use the inspiral and merger of massive black hole binaries as standard sirens, extending our ability to characterize the expansion history of the Universe from the onset of dark energy-domination out to a redshift z ~ 10. A low-frequency detector, furthermore, offers the best chance for discovery of exotic gravitational wave sources, including a primordial stochastic background, that could reveal clues to the origin of our Universe.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe.

We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies.

The origin and properties of the initial black hole seeds that grow to produce the observed population of accreting sources remain to be determined. It is a challenge to uniquely disentangle signatures of seeding from fueling and dynamics that shapes the assembly history of growing black holes. To address this, we use a semi-analytic model developed to track the growth of supermassive black holes adopting multiple prescriptions for accretion. In contrast with earlier treatments, we explore the interplay between seeding models and two accretion modes. We find that signatures of the initial seeding do survive in the following observational probes: the black hole occupation fraction; contribution to the unresolved X-ray background; low-luminosity and high-redshift luminosity functions; and in gravitational wave event signatures. We find that the behaviour of the low-mass end of the M-sigma relation is dominated by uncertainties in the adopted accretion prescriptions and does not offer clear discrimination between seeding models. We make concrete predictions for future surveys, particularly for the Lynx X-ray surveyor and LISA (The Laser Interferometer Space Antenna) mission, which will each provide different and yet strong constraints on the seed population. Black hole coalescences detected by LISA and high-redshift quasar luminosity functions observed by Lynx will offer the sharpest seeding discriminants. Although the signatures of the black hole seeding mechanism that persist remain linked to our understanding of black hole accretion and dynamics, we offer new insights on how these upcoming multi-wavelength observations could be leveraged to effectively disentangle them.

We present a gravitational lensing and X-ray analysis of a massive galaxy cluster and its surroundings. The core of MACS\,J0717.5+3745 ($M(R<1\,{\rm Mpc})\sim$\,$2$$\times$$10^{15}\,\msun$, $z$=$0.54$) is already known to contain four merging components. We show that this is surrounded by at least seven additional substructures with masses ranging from $3.8-6.5\times10^{13}\,\msun$, at projected radii $1.6$ to $4.9$\,Mpc. We compare MACS\,J0717 to mock lensing and X-ray observations of similarly rich clusters in cosmological simulations. The low gas fraction of substructures predicted by simulations turns out to match our observed values of $1$--$4\%$. Comparing our data to three similar simulated halos, we infer a typical growth rate and substructure infall velocity. That suggests MACS\,J0717 could evolve into a system similar to, but more massive than, Abell\,2744 by $z=0.31$, and into a $\sim$\,$10^{16}\,\msun$ supercluster by $z=0$. The radial distribution of infalling substructure suggests that merger events are strongly episodic; however we find that the smooth accretion of surrounding material remains the main source of mass growth even for such massive clusters.

In this paper we report a systematic search for an emission line around 3.5 keV in the spectrum of the Cosmic X-ray Background using a total of $\sim$10 Ms Chandra observations towards the COSMOS Legacy and CDFS survey fields. We find a marginal evidence of a feature at an energy of $\sim$3.51 keV with a significance of 2.5-3 $\sigma$, depending on the choice of the statistical treatment. The line intensity is best fit at $8.8\ \pm\ {2.9}\times10^{-7}$ ph cm$^{-2}$s$^{-1}$ when using a simple $\Delta\chi^2$ or $10.2\ ^{+0.2}_{-0.4} \times10^{-7}$ ph cm$^{-2}$s$^{-1}$ when MCMC is used. Based on our knowledge of $Chandra$, and the reported detection of the line by other instruments, an instrumental origin for the line remains unlikely. We cannot though rule out a statistical fluctuation and in that case our results provide a 3$\sigma$ upper limit at 1.85$\times$10$^{-6}$ ph cm$^{-2}$s$^{-1}$. We discuss the interpretation of this observed line in terms of the iron line background; S {\sc XVI} charge exchange as well as potentially from sterile neutrino decay. We note that our detection is consistent with previous measurements of this line toward the Galactic center, and can be modeled as the result of sterile neutrino decay from the Milky Way for the dark matter distribution modeled as an NFW profile. For this case, we estimate a mass m$_{\nu}\sim$7.01 keV and a mixing angle sin$^2$(2$\theta$)= 0.83--2.75 $\times 10^{-10}$. These derived values are in agreement with independent estimates from galaxy clusters; the Galactic center and M31.

Super-massive black holes weighing up to $\sim 10^9 \, \mathrm{M_{\odot}}$ are in place by $z \sim 7$, when the age of the Universe is $\lesssim 1 \, \mathrm{Gyr}$. This implies a time crunch for their growth, since such high masses cannot be easily reached in standard accretion scenarios. Here, we explore the physical conditions that would lead to optimal growth wherein stable super-Eddington accretion would be permitted. Our analysis suggests that the preponderance of optimal conditions depends on two key parameters: the black hole mass and the host galaxy central gas density. In the high-efficiency region of this parameter space, a continuous stream of gas can accrete onto the black hole from large to small spatial scales, assuming a global isothermal profile for the host galaxy. Using analytical initial mass functions for black hole seeds, we find an enhanced probability of high-efficiency growth for seeds with initial masses $\gtrsim 10^4 \, \mathrm{M_{\odot}}$. Our picture suggests that a large population of high-$z$ lower-mass black holes that formed in the low-efficiency region, with low duty cycles and accretion rates, might remain undetectable as quasars, since we predict their bolometric luminosities to be $\lesssim 10^{41} \, \mathrm{erg \, s^{-1}}$. The presence of these sources might be revealed only via gravitational wave detections of their mergers.

We develop a semi-analytic model to explore different prescriptions of supermassive black hole (SMBH) fuelling. This model utilises a merger-triggered burst mode in concert with two possible implementations of a long-lived steady mode for assembling the mass of the black hole in a galactic nucleus. We improve modelling of the galaxy-halo connection in order to more realistically determine the evolution of a halo's velocity dispersion. We use four model variants to explore a suite of observables: the M-sigma relation, mass functions of both the overall and broad-line quasar population, and luminosity functions as a function of redshift. We find that "downsizing" is a natural consequence of our improved velocity dispersion mappings, and that high-mass SMBHs assemble earlier than low-mass SMBHs. The burst mode of fuelling is sufficient to explain the assembly of SMBHs to z=2, but an additional steady mode is required to both assemble low-mass SMBHs and reproduce the low-redshift luminosity function. We discuss in detail the trade-offs in matching various observables and the interconnected modelling components that govern them. As a result, we demonstrate the utility as well as the limitations of these semi-analytic techniques.

Multi-wavelength surveys covering large sky volumes are necessary to obtain an accurate census of rare objects such as high luminosity and/or high redshift active galactic nuclei (AGN). Stripe 82X is a 31.3 deg$^2$ X-ray survey with $Chandra$ and $XMM$-Newton observations overlapping the legacy Sloan Digital Sky Survey (SDSS) Stripe 82 field, which has a rich investment of multi-wavelength coverage from the ultraviolet to the radio. The wide-area nature of this survey presents new challenges for photometric redshifts for AGN compared to previous work on narrow-deep fields because it probes different populations of objects that need to be identified and represented in the library of templates. Here we present an updated X-ray plus multi-wavelength matched catalog, including $Spitzer$ counterparts, and estimated photometric redshifts for 5961 (96% of a total of 6181) X-ray sources, which have a normalized median absolute deviation, $\sigma_{\rm nmad}$ = 0.06 and an outlier fraction, $\eta$ = 13.7%. The populations found in this survey, and the template libraries used for photometric redshifts, provide important guiding principles for upcoming large-area surveys such as $eROSITA$ and 3$XMM$ (in X-ray) and the Large Synoptic Survey Telescope (LSST; optical).

We present new measurements of the large scale clustering component of the cross-power spectra of the source-subtracted Spitzer-IRAC Cosmic Infrared Background (CIB) and Chandra-ACIS Cosmic X-ray Background (CXB) surface brightness fluctuations. Our investigation uses data from the Chandra Deep Field South (CDFS), Hubble Deep Field North (HDFN), EGS/AEGIS field and UDS/SXDF surveys, comprising 1160 Spitzer hours and $\sim$ 12 Ms of Chandra data collected over a total area of 0.3 deg$^2$. We report the first ($>$5$\sigma$) detection of a cross-power signal on large angular scales >20 arcsec between [0.5-2] keV and the 3.6 and 4.5$\mu$m bands, at $\sim$5$\sigma$ and 6.3$\sigma$ significance, respectively. The correlation with harder X-ray bands is marginally significant. Comparing the new observations with existing models for the contribution of the known unmasked source population at $z<$7, we find an excess of about an order of magnitude at 5$\sigma$ confidence. We discuss possible interpretations for the origin of this excess in terms of the contribution from accreting early black holes, including both direct collapse black holes and primordial black holes, as well as from scattering in the interstellar medium and intra-halo light.

Growing the supermassive black holes (~10^9 Msun) that power the detected luminous, highest redshift quasars (z > 6) from light seeds - the remnants of the first stars - within ~ 1 Gyr of the Big Bang poses a timing challenge for growth models. The formation of massive black hole seeds via direct collapse with initial masses ~ 10^4 - 10^5 Msun alleviates this problem. Physical conditions required to form these massive direct collapse black hole (DCBH) seeds are available in the early universe. These viable DCBH formation sites, satellite halos of star-forming galaxies, merge and acquire a stellar component. These produce a new, transient class of objects at high redshift, Obese Black hole Galaxies (OBGs), where the luminosity produced by accretion onto the black hole outshines the stellar component. Therefore, the OBG stage offers a unique way to discriminate between light and massive initial seeds. We predict the multi-wavelength energy output of OBGs and growing Pop III remnants at a fiducial redshift (z = 9), exploring both standard and slim disk accretion onto the growing central black hole for high and low metallicities of the associated stellar population. With our computed templates, we derive the selection criteria for OBGs, that comprise a pre-selection that eliminates blue sources; followed by color-color cuts ([F_{070W} - F_{220W}] > 0; -0.3 < [F_{200W} - F_{444W}] < 0.3) and when available, a high ratio of X-ray flux to rest-frame optical flux (F_X/F_{444W} >> 1) (Abridged).

We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations.

The most massive black holes observed in the Universe weigh up to $\sim 10^{10} \, \mathrm{M_{\odot}}$, nearly independent of redshift. Reaching these final masses likely required copious accretion and several major mergers. Employing a dynamical approach, that rests on the role played by a new, relevant physical scale - the transition radius - we provide a theoretical calculation of the maximum mass achievable by a black hole seed that forms in an isolated halo, one that scarcely merged. Incorporating effects at the transition radius and their impact on the evolution of accretion in isolated haloes we are able to obtain new limits for permitted growth. We find that large black hole seeds ($M_{\bullet} \gtrsim 10^4 \, \mathrm{M_{\odot}}$) hosted in small isolated halos ($M_h \lesssim 10^9 \, \mathrm{M_{\odot}}$) accreting with relatively small radiative efficiencies ($\epsilon \lesssim 0.1$) grow optimally in these circumstances. Moreover, we show that the standard $M_{\bullet}-\sigma$ relation observed at $z \sim 0$ cannot be established in isolated halos at high-$z$, but requires the occurrence of mergers. Since the average limiting mass of black holes formed at $z \gtrsim 10$ is in the range $10^{4-6} \, \mathrm{M_{\odot}}$, we expect to observe them in local galaxies as intermediate-mass black holes, when hosted in the rare haloes that experienced only minor or no merging events. Such ancient black holes, formed in isolation with subsequent scant growth, could survive, almost unchanged, until present.

We estimate the accretion rates onto the supermassive black holes that power 20 of the highest-redshift quasars, at z>5.8, including the quasar with the highest redshift known to date -- ULAS J1120 at z=7.09. The analysis is based on the observed (rest-frame) optical luminosities and reliable "virial" estimates of the BH masses of the quasars, and utilizes scaling relations derived from thin accretion disk theory. The mass accretion rates through the postulated disks cover a wide range, dM_disk/dt~4-190 Msol/yr, with most of the objects (80\%) having dM_disk/dt~10-65 Msol/yr, confirming the Eddington-limited nature of the accretion flows. By combining our estimates of dM_disk/dt with conservative, lower limits on the bolometric luminosities of the quasars, we investigate which alternative values of \eta best account for all the available data. We find that the vast majority of quasars (~85\%) can be explained with radiative efficiencies in the range \eta~0.03-0.3, with a median value close to the commonly assumed \eta=0.1. Within this range, we obtain conservative estimates of \eta>0.14 for ULAS J1120 and SDSS J0100 (at z=6.3), and of >0.19 for SDSS J1148 (at z=6.41; assuming their BH masses are accurate). The implied accretion timescales are generally in the range t_acc = M_BH / dM_BH/dt ~0.1-1 Gyr, suggesting that most quasars could have had ~1-10 mass e-foldings since BH seed formation. Our analysis therefore demonstrates that the available luminosities and masses for the highest-redshift quasars can be explained self-consistently within the thin, radiatively efficient accretion disk paradigm. Episodes of radiatively inefficient, "super-critical" accretion may have occurred at significantly earlier epochs (i.e., z~10 or above).

Using {\em Chandra} observations in the 2.15 deg$^{2}$ COSMOS legacy field, we present one of the most accurate measurements of the Cosmic X-ray Background (CXB) spectrum to date in the [0.3-7] keV energy band. The CXB has three distinct components: contributions from two Galactic collisional thermal plasmas at kT$\sim$0.27 and 0.07 keV and an extragalactic power-law with photon spectral index $\Gamma$=1.45$\pm{0.02}$. The 1 keV normalization of the extragalactic component is 10.91$\pm{0.16}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. Removing all X-ray detected sources, the remaining unresolved CXB is best-fit by a power-law with normalization 4.18$\pm{0.26}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$ and photon spectral index $\Gamma$=1.57$\pm{0.10}$. Removing faint galaxies down to i$_{AB}\sim$27-28 leaves a hard spectrum with $\Gamma\sim$1.25 and a 1 keV normalization of $\sim$1.37 keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. This means that $\sim$91\% of the observed CXB is resolved into detected X-ray sources and undetected galaxies. Unresolved sources that contribute $\sim 8-9\%$ of the total CXB show a marginal evidence of being harder and possibly more obscured than resolved sources. Another $\sim$1\% of the CXB can be attributed to still undetected star forming galaxies and absorbed AGN. According to these limits, we investigate a scenario where early black holes totally account for non source CXB fraction and constrain some of their properties. In order to not exceed the remaining CXB and the $z\sim$6 accreted mass density, such a population of black holes must grow in Compton-thick envelopes with N$_{H}>$1.6$\times$10$^{25}$ cm$^{-2}$ and form in extremely low metallicity environments $(Z_\odot)\sim10^{-3}$.

The massive substructures found in Abell 2744 by Jauzac et al. (2016) present a challenge to the cold dark matter paradigm due to their number and proximity to the cluster centre. We use one of the biggest N-body simulations, the Millennium XXL, to investigate the substructure in a large sample of massive dark matter haloes. A range of effects which influence the comparison with the observations is considered, extending the preliminary evaluation carried out by Jauzac et al. (2016). There are many tens of haloes in the simulation with a total mass comparable with or larger than that of Abell 2744. However, we find no haloes with a number and distribution of massive substructures (> 5 10^13 Msun) that is close to that inferred from the observations of Abell 2744. The application of extreme value statistics suggests that we would need a simulation of at least ten times the volume of the Millennium XXL to find a single dark matter halo with a similar internal structure to Abell 2744. Explaining the distribution of massive substructures in clusters is a new hurdle for hierarchical models to negotiate, which is not weakened by appeals to baryonic physics or uncertainty over the nature of the dark matter particle.

We propose that one of the sources in the recently detected system CR7 by Sobral et al. (2015) through spectro-photometric measurements at $z = 6.6$ harbors a direct collapse blackhole (DCBH). We argue that the LW radiation field required for direct collapse in source A is provided by sources B and C. By tracing the LW production history and star formation rate over cosmic time for the halo hosting CR7 in a $\Lambda$CDM universe, we demonstrate that a DCBH could have formed at $z\sim 20$. The spectrum of source A is well fit by nebular emission from primordial gas around a BH with MBH $\sim 4.4 \times 10^6 \ \rm M_{\odot}$ accreting at a 40 % of the Eddington rate, which strongly supports our interpretation of the data. Combining these lines of evidence, we argue that CR7 might well be the first DCBH candidate.

Direct collapse black holes (DCBH) have been proposed as a solution to the challenge of assembling supermassive black holes by $z>6$ to explain the bright quasars observed at this epoch. The formation of a DCBH seed with $\rm M_{BH}\sim10^{4-5}\ \rm M_{\odot}$ requires a pristine atomic-cooling halo to be illuminated by an external radiation field that is sufficiently strong to entirely suppress H$_{2}$ cooling in the halo. Many previous studies have attempted to constrain the critical specific intensity that is likely required to suppress H$_{2}$ cooling, denoted as $J_{\rm crit}$. However, these studies have typically assumed that the incident external radiation field can be modeled with a black-body spectrum. Under this assumption, it is possible to derive a {unique} value for $J_{\rm crit}$ that depends only on the temperature of the black-body. In this study we consider a more realistic spectral energy distribution (SED) for the external source of radiation that depends entirely on its star formation history and age. The rate of destruction of the species responsible for suppressing molecular hydrogen cooling depends on the detailed shape of the SED. Therefore the value of $J_{\rm crit}$ is tied to the shape of the incident SED of the neighbouring galaxy. We fit a parametric form to the rates of destruction of H$_2$ and H$^-$ that permit direct collapse. Owing to this, we find that $J_{\rm crit}$ is not a fixed threshold but can lie anywhere in the range $J_{\rm crit} \sim 0.5$--$10^{3}$, depending on the details of the source stellar population.

We present a joint optical/X-ray analysis of the massive galaxy cluster Abell 2744 (z=0.308). Our strong- and weak-lensing analysis within the central region of the cluster, i.e., at R<1Mpc from the brightest cluster galaxy, reveals eight substructures, including the main core. All of these dark-matter halos are detected with a significance of at least 5sigma and feature masses ranging from 0.5 to 1.4x10^{14}Msun within R<150kpc. Merten et al. (2011) and Medezinski et al. (2016) substructures are also detected by us. We measure a slightly higher mass for the main core component than reported previously and attribute the discrepancy to the inclusion of our tightly constrained strong-lensing mass model built on Hubble Frontier Fields data. X-ray data obtained by XMM-Newton reveal four remnant cores, one of them a new detection, and three shocks. Unlike Merten et al. (2011), we find all cores to have both dark and luminous counterparts. A comparison with clusters of similar mass in the MXXL simulations yields no objects with as many massive substructures as observed in Abell 2744, confirming that Abell 2744 is an extreme system. We stress that these properties still do not constitute a challenge to $\Lambda$CDM, as caveats apply to both the simulation and the observations: for instance, the projected mass measurements from gravitational lensing and the limited resolution of the sub-haloes finders. We discuss implications of Abell 2744 for the plausibility of different dark-matter candidates and, finally, measure a new upper limit on the self-interaction cross-section of dark matter of sigma_{DM}<1.28cm2/g(68\% CL), in good agreement with previous results from Harvey et al. (2015).

Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes.

The Frontier Fields are a director's discretionary time campaign with HST and the Spitzer Space Telescope to see deeper into the universe than ever before. The Frontier Fields combine the power of HST and Spitzer with the natural gravitational telescopes of massive high-magnification clusters of galaxies to produce the deepest observations of clusters and their lensed galaxies ever obtained. Six clusters - Abell 2744, MACSJ0416.1-2403, MACSJ0717.5+3745, MACSJ1149.5+2223, Abell S1063, and Abell 370 - were selected based on their lensing strength, sky darkness, Galactic extinction, parallel field suitability, accessibility to ground-based facilities, HST, Spitzer and JWST observability, and pre-existing ancillary data. These clusters have been targeted by the HST ACS/WFC and WFC3/IR with coordinated parallels of adjacent blank fields for over 840 HST orbits. The Spitzer Space Telescope has dedicated > 1000 hours of director's discretionary time to obtain IRAC 3.6 and 4.5 micron imaging to ~26.5, 26.0 ABmag 5-sigma point-source depths in the six cluster and six parallel Frontier Fields. The Frontier Field parallel fields are the second-deepest observations thus far by HST with ~29th ABmag 5-sigma point source depths in seven optical - near-infrared bandpasses. Galaxies behind the Frontier Field cluster lenses experience typical magnification factors of a few, with small regions near the critical curves magnified by factors 10-100. Therefore, the Frontier Field cluster HST images achieve intrinsic depths of ~30-33 magnitudes over very small volumes. Early studies of the Frontier Fields have probed galaxies fainter than any seen before during the epoch of reionization 6 < z < 10, mapped out the cluster dark matter to unprecedented resolution, and followed lensed transient events.

Mass determinations from gravitational lensing shear and the higher order estimator flexion are both subject to the mass sheet degeneracy. Mass sheet degeneracy refers to a transformation that leaves the reduced shear and flexion invariant. In general, this transformation can be approximated by the addition of a constant surface mass density sheet. We propose a new technique to break the mass sheet degeneracy. The method uses mass moments of the shear or flexion fields in combination with convergence information derived from number counts which exploit the magnification bias. The difference between the measured mass moments provides an estimator for the magnitude of the additive constant that is the mass-sheet. For demonstrating this, we derive relations that hold true in general for n-th order moments and show how they can be employed effectively to break the degeneracy. We investigate the detectability of this degeneracy parameter from our method and find that the degeneracy parameter can be feasibly determined from stacked galaxy-galaxy lensing data and cluster lensing data. Furthermore, we compare the signal-to-noise ratios of convergence information from number counts with shear and flexion for SIS and NFW models. We find that the combination of shear and flexion performs best on galaxy and cluster scales and the convergence information can therefore be used to break the mass sheet degeneracy without quality loss in the mass reconstruction. In summary, there is power in the combination of shear, flexion, convergence and their higher order moments. With the anticipated wealth of lensing data from upcoming and future satellite missions - EUCLID and WFIRST - this technique will be feasible.

[abridged] We present a strong-lensing analysis of MACSJ0717.5+3745, based on the full depth of the Hubble Frontier Field (HFF) observations, which brings the number of multiply imaged systems to 61, ten of which are spectroscopically confirmed. The total number of images comprised in these systems rises to 165. Our analysis uses a parametric mass reconstruction technique, as implemented in the Lenstool software, to constrain a mass distribution composed of four large-scale mass components + galaxy-scale perturbers. We find a superposition of cored isothermal mass components to provide a good fit to the observational constraints, resulting in a very shallow mass distribution for the smooth (large-scale) component. Given the implications of such a flat mass profile, we investigate whether a model composed of "peaky" non-cored mass components can also reproduce the observational constraints. We find that such a non-cored mass model reproduces the observational constraints equally well. Although the total mass distributions of both models are consistent, as well as the integrated two dimensional mass profiles, we find that the smooth and the galaxy-scale components are very different. We conclude that, even in the HFF era, the generic degeneracy between smooth and galaxy-scale components is not broken, in particular in such a complex galaxy cluster. Consequently, insights into the mass distribution of MACS J0717 remain limited, underlining the need for additional probes beyond strong lensing. Our findings also have implications for estimates of the lensing magnification: we show that the amplification difference between the two models is larger than the error associated with either model. This uncertainty decreases the area of the image plane where we can reliably study the high-redshift Universe by 50 to 70%.

Capitalizing on the all-sky coverage of {\it WISE}, and the 35\% and 50\% sky coverage from SDSS and Pan-STARRS, respectively, we explore the efficacy of $m_{R}$ (optical) - $m_{3.4 \mu m}$ (mid-infrared), hereafter $R-W1$, as a color diagnostic to identify obscured supermassive black hole accretion in wide-area X-ray surveys. We use the $\sim$16.5 deg$^2$ Stripe 82 X-ray survey data as a test-bed to compare $R-W1$ with $R-K$, an oft-used obscured AGN selection criterion, and examine where different classes of objects lie in this parameter space. Most stars follow a well-defined path in $R-K$ vs. $R-W1$ space. We demonstrate that optically normal galaxies hosting X-ray AGN at redshifts $0.5<z<1$ can be recovered with an $R-W1>4$ color-cut, while they typically are not selected as AGN based on their $W1-W2$ colors. Additionally, different observed X-ray luminosity bins favor different regions in $R-W1$ parameter space: moderate luminosity AGN ($10^{43}$ erg s$^{-1} < L_{\rm 0.5-10 keV} < 10^{44}$ erg s$^{-1}$) tend to have red colors while the highest luminosity AGN ($L_{\rm 0.5-10 keV} > 10^{45}$ erg s$^{-1}$) have bluer colors; higher spectroscopic completeness of the Stripe 82X sample is needed to determine whether this is a selection effect or an intrinsic property. Finally, we parameterize X-ray obscuration of Stripe 82X AGN by calculating their hardness ratios (HRs) and find no clear trends between HR and optical reddening. Our results will help inform best-effort practices in following-up obscured AGN candidates in current and future wide-area, shallow X-ray surveys, including the all-sky {\it eROSITA} mission.

We present a high-precision mass model of the galaxy cluster MACSJ1149.6+2223, based on a strong-gravitational-lensing analysis of Hubble Space Telescope Frontier Fields (HFF) imaging data and spectroscopic follow-up with Gemini/GMOS and VLT/MUSE. Our model includes 12 new multiply imaged galaxies, bringing the total to 22, comprised of 65 individual lensed images. Unlike the first two HFF clusters, Abell 2744 and MACSJ0416.1-2403, MACSJ1149 does not reveal as many multiple images in the HFF data. Using the Lenstool software package and the new sets of multiple images, we model the cluster with several cluster-scale dark-matter halos and additional galaxy-scale halos for the cluster members. Consistent with previous analyses, we find the system to be complex, composed of five cluster-scale halos. Their spatial distribution and lower mass, however, makes MACSJ1149 a less powerful lens. Our best-fit model predicts image positions with an RMS of 0.91". We measure the total projected mass inside a 200~kpc aperture as ($1.840\pm 0.006$)$\times 10^{14}$M$_{\odot}$, thus reaching again 1% precision, following our previous HFF analyses of MACSJ0416.1-2403 and Abell 2744. In light of the discovery of the first resolved quadruply lensed supernova, SN Refsdal, in one of the multiply imaged galaxies identified in MACSJ1149, we use our revised mass model to investigate the time delays and predict the rise of the next image between November 2015 and January 2016.

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

Massive black hole binaries (MBHBs) are a natural byproduct of galaxy mergers. Previous studies have shown that flares from stellar tidal disruption events (TDEs) are modified by the presence of a secondary perturber, causing interruptions in the light curve. We study the dynamics of TDE debris in the presence of a milliparsec-separated MBHB by integrating ballistic particle orbits in the time-varying potential of the binary. We find that gaps in the light curve appear when material misses the accretion radius on its first return to pericentre. Subsequent recurrences can be decomposed into "continuous" and "delayed" components, which exhibit different behaviour. We find that this potential can substantially alter the locations of stream self-intersections. When debris is confined to the plane, we find that close encounters with the secondary BH leave noticeable signatures on the fallback rate and can result in significant accretion onto the secondary BH. Tight, equal-mass MBHBs accrete equally, periodically trading the infalling stream.

We use deep Hubble Space Telescope imaging of the Frontier Fields to accurately measure the galaxy rest-frame ultraviolet luminosity function (UV LF) in the redshift range $z \sim 6-8$. We combine observations in three lensing clusters A2744, MACS0416, MACS0717 and their associated parallels fields to select high-redshift dropout candidates. We use the latest lensing models to estimate the flux magnification and the effective survey volume in combination with completeness simulations performed in the source plane. We report the detection of 227 galaxy candidates at $z=6-7$ and 25 candidates at $z \sim 8$. While the total survey area is about 4 arcmin$^{2}$ in each parallel field, it drops to about 0.6 to 1 arcmin$^{2}$ in the cluster core fields because of the strong lensing. We compute the UV luminosity function at $z \sim 7$ using the combined galaxy sample and perform Monte Carlo simulations to determine the best fit Schechter parameters. We are able to reliably constrain the LF down to an absolute magnitude of $M_{UV}=-15.25$, which corresponds to 0.005$L^{\star}$. More importantly, we find that the faint-end slope remains steep down to this magnitude limit with $\alpha=-2.04_{-0.17}^{+0.13}$. Our results confirm the most recent results in deep blank fields but extend the LF measurements more than two magnitudes deeper. The UV LF at $z \sim 8$ is not very well constrained below $M_{UV}=-18$ due to the small number statistics and incompleteness uncertainties. To assess the contribution of galaxies to cosmic reionization we derive the UV luminosity density at $z\sim7$ by integrating the UV LF down to an observationally constrained limit of $M_{UV} = -15$. We show that our determination of Log($\rho_{UV}$)=$26.2\pm0.13$ (erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$) can be sufficient to maintain the IGM ionized.

We present a high-precision mass model of galaxy cluster Abell 2744, based on a strong-gravitational-lensing analysis of the \emph{Hubble Space Telescope Frontier Fields} (HFF) imaging data, which now include both \emph{Advanced Camera for Surveys} and \emph{Wide-Field Camera 3} observations to the final depth. Taking advantage of the unprecedented depth of the visible and near-infrared data, we identify 34 new multiply imaged galaxies, bringing the total to 61, comprising 181 individual lensed images. In the process, we correct previous erroneous identifications and positions of multiple systems in the northern part of the cluster core. With the \textsc{Lenstool} software and the new sets of multiple images, we model the cluster using two cluster-scale dark matter halos plus galaxy-scale halos for the cluster members. Our best-fit model predicts image positions with an \emph{RMS} error of 0.69$\arcsec$, which constitutes an improvement by almost a factor of two over previous parametric models of this cluster. We measure the total projected mass inside a 200~kpc aperture as ($2.162\pm 0.005$)$\times 10^{14}M_{\odot}$, thus reaching 1\% level precision for the second time, following the recent HFF measurement of MACSJ0416.1-2403. Importantly, the higher quality of the mass model translates into an overall improvement by a factor of 4 of the derived magnification factor. % for the high-redshift lensed background galaxies. Together with our previous HFF gravitational lensing analysis, this work demonstrates that the HFF data enables high-precision mass measurements for massive galaxy clusters and the derivation of robust magnification maps to probe the early Universe.

Accretion is thought to primarily contribute to the mass accumulation history of supermassive black holes throughout cosmic time. While this may be true at high redshifts, at lower redshifts and for the most massive black holes mergers themselves might add significantly to the mass budget. We evolve SMBHs from $4 > z > 0$ using merger trees derived from hydrodynamical cosmological simulations of a cluster and void region, scaled to the observed value of the stellar mass fraction to account for overcooling. Mass gains from gas accretion proportional to bulge growth and BH-BH mergers are tracked, as are black holes that remain "orbiting" due to insufficient dynamical friction in a merger remnant, as well as those that are ejected due to gravitational recoil. We find that gas accretion remains the dominant source of mass accumulation in almost all SMBHs; mergers contribute $2.5\pm0.1\%$ for all SMBHs in the cluster and $1.0\pm0.1\%$ in the void since $z = 4$. However, mergers are significant for massive SMBHs. The fraction of mass accumulated from mergers for central BHs generally increases for larger values of the host bulge mass: in the void, the fraction is $2\%$ at $M_{*, bul} = 10^{10} M_{\odot}$, increasing to $4\%$ at $M_{*, bul} \gtrsim 10^{11} M_{\odot}$, and in the cluster it is $4\%$ at $M_{*, bul} = 10^{10} M_{\odot}$ and $23\%$ at $10^{12} M_{\odot}$. We find that $40\%$ of SMBHs and $\approx 8\%$ of the total SMBH mass is found orbiting in the cluster region at $z = 0$. The existence of orbiting and ejected SMBHs requires modification of the Soltan argument. We estimate this correction to the integrated accreted mass density of SMBHs to be in the range $6-21\%$, with a mean value of $11\pm3\%$. We also calculate the total energy output and strain from gravitational waves emitted by merging SMBHs, and obtain a signal potentially detectable by pulsar timing arrays.

We use a joint optical/X-ray analysis to constrain the geometry and history of the ongoing merging event in the massive galaxy cluster MACSJ0416.1-2403 (z=0.397). Our investigation of cluster substructure rests primarily on a combined strong- and weak-lensing mass reconstruction based on the deep, high-resolution images obtained for the Hubble Frontier Fields initiative. To reveal the system's dynamics, we complement this lensing analysis with a study of the intra-cluster gas using shallow Chandra data, and a three-dimensional model of the distribution and motions of cluster galaxies derived from over 100 spectroscopic redshifts. The multi-scale grid model obtained from our combined lensing analysis extends the high-precision strong-lensing mass reconstruction recently performed to cluster-centric distances of almost 1 Mpc. Our analysis detects the two well known mass concentrations in the cluster core. A pronounced offset between collisional and collisionless matter is only observed for the SW cluster component, while excellent alignment is found for the NE cluster. Both the lensing analysis and the distribution of cluster light strongly suggest the presence of a third massive structure, almost 2 arcmin SW of the cluster centre. Since no X-ray emission is detected in this region, we conclude that this structure is non-virialised and speculate that it might be part of a large-scale filament almost aligned with our line of sight. Combining all evidence from the distribution of dark and luminous matter, we propose two alternative scenarios for the trajectories of the components of MACSJ0416.1-2403. Upcoming deep X-ray observations that allow the detection of shock fronts, cold cores, and sloshing gas (all key diagnostics for studies of cluster collisions) will allow us to test, and distinguish between these two scenarios.

Exploiting the power of gravitational lensing, the Hubble Frontier Fields (HFF) program aims at observing six massive galaxy clusters to explore the distant Universe far beyond the depth limits of blank field surveys. Using the complete Hubble Space Telescope observations of the first HFF cluster Abell 2744, we report the detection of 50 galaxy candidates at $z \sim 7$ and eight candidates at $z \sim 8$ in a total survey area of 0.96 arcmin$^{2}$ in the source plane. Three of these galaxies are multiply-imaged by the lensing cluster. Using an updated model of the mass distribution in the cluster we were able to calculate the magnification factor and the effective survey volume for each galaxy in order to compute the ultraviolet galaxy luminosity function at both redshifts 7 and 8. Our new measurements extend the $z \sim 7$ UV LF down to an absolute magnitude of $M_{UV} \sim -15.5$. We find a characteristic magnitude of $M^{\star}_{UV}=-20.63^{+0.69}_{-0.56}$ mag and a faint-end slope $\alpha = -1.88^{+0.17}_{-0.20}$ close to previous determinations in blank fields. We show here for the first time that this slope remains steep down to very faint luminosities of 0.01$L^{\star}$. Although prone to large uncertainties, our results at $z \sim 8$ also seem to confirm a steep faint-end slope below 0.1$L^{\star}$. The HFF program is therefore providing an extremely efficient way to study the faintest galaxy populations at $z > 7$ that would otherwise be inaccessible with current instrumentation. The full sample of six galaxy clusters will provide yet better constraints on the build-up of galaxies at early epochs and their contribution to cosmic reionization.

Cluster gravitational lensing surveys like the Hubble Space Telescope Frontier Fields survey will detect distant galaxies 10-50 times fainter than any yet discovered. Using these surveys to measure the luminosity function of such faint, distant galaxies, however, requires that magnification maps built from the constraints of strongly-lensed images be accurate. For models that assume the cluster and nearby (correlated) structures are the only significant sources of lensing, a potential source of error in these maps comes from the fact that light rays also suffer weak deflections by uncorrelated large-scale structure along the line-of-sight, i.e. cosmic weak lensing (CWL). To demonstrate the magnitude of this effect, we calculate the magnification change which results when the same cluster-lens is placed along different lines of sight. Using a simple density profile for a cluster-lens at z~0.3-0.5 and the power spectrum of the matter density fluctuations responsible for CWL, we show that the typical magnifications of ~5(10) of sources at z=6-10 can differ by ~10-20(20-30)% from one line-of-sight to another. However, these fluctuations rise to greater than order unity near critical curves, indicating that CWL tends to make its greatest contribution to the most magnified images. We conclude that the neglect of CWL in determining the intrinsic luminosities of highly-magnified galaxies may introduce errors significant enough to warrant further effort to include this contribution in cluster-lens modeling. We suggest that methods of modeling CWL in galaxy-strong-lensing systems should be generalized to cluster-lensing systems.

Mass accretion by black holes (BHs) is typically capped at the Eddington rate, when radiation's push balances gravity's pull. However, even exponential growth at the Eddington-limited e-folding time t_E ~ few x 0.01 billion years, is too slow to grow stellar-mass BH seeds into the supermassive luminous quasars that are observed when the universe is 1 billion years old. We propose a dynamical mechanism that can trigger supra-exponential accretion in the early universe, when a BH seed is trapped in a star cluster fed by the ubiquitous dense cold gas flows. The high gas opacity traps the accretion radiation, while the low-mass BH's random motions suppress the formation of a slowly-draining accretion disk. Supra-exponential growth can thus explain the puzzling emergence of supermassive BHs that power luminous quasars so soon after the Big Bang.

Extending over three Hubble Space Telescope (HST) cycles, the Hubble Frontier Fields (HFF) initiative constitutes the largest commitment ever of HST time to the exploration of the distant Universe via gravitational lensing by massive galaxy clusters. We here present models of the mass distribution in the six HFF cluster lenses, derived from a joint strong- and weak-lensing analysis anchored by a total of 88 multiple-image systems identified in existing HST data. The resulting maps of the projected mass distribution and of the gravitational magnification effectively calibrate the HFF clusters as gravitational telescopes. Allowing the computation of search areas in the source plane, these maps are provided to the community to facilitate the exploitation of forthcoming HFF data for quantitative studies of the gravitationally lensed population of background galaxies. Our models of the gravitational magnification afforded by the HFF clusters allow us to quantify the lensing-induced boost in sensitivity over blank-field observations and predict that galaxies at $z>10$ and as faint as m(AB)=32 will be detectable, up to 2 magnitudes fainter than the limit of the Hubble Ultra Deep Field.

We present a high-precision mass model of the galaxy cluster MACSJ0416.1-2403, based on a strong-gravitational-lensing analysis of the recently acquired Hubble Space Telescope Frontier Fields (HFF) imaging data. Taking advantage of the unprecedented depth provided by HST/ACS observations in three passbands, we identify 51 new multiply imaged galaxies, quadrupling the previous census and bringing the grand total to 68, comprising 194 individual lensed images. Having selected a subset of the 57 most securely identified multiply imaged galaxies, we use the Lenstool software package to constrain a lens model comprised of two cluster-scale dark-matter halos and 98 galaxy-scale halos. Our best-fit model predicts image positions with an $RMS$ error of 0.68'', which constitutes an improvement of almost a factor of two over previous, pre-HFF models of this cluster. We find the total projected mass inside a 200~kpc aperture to be $(1.60\pm0.01)\times 10^{14}\ M_\odot$, a measurement that offers a three-fold improvement in precision, reaching the percent level for the first time in any cluster. Finally, we quantify the increase in precision of the derived gravitational magnification of high-redshift galaxies and find an improvement by a factor of $\sim$2.5 in the statistical uncertainty. Our findings impressively confirm that HFF imaging has indeed opened the domain of high-precision mass measurements for massive clusters of galaxies.

The Hubble Frontier Fields (HFF) program combines the capabilities of the Hubble Space Telescope (HST) with the gravitational lensing of massive galaxy clusters to probe the distant Universe to an unprecedented depth. Here, we present the results of the first combined HST and Spitzer observations of the cluster Abell 2744. We combine the full near-infrared data with ancillary optical images to search for gravitationally lensed high-redshift (z > 6) galaxies. We report the detection of 15 I814-dropout candidates at z ~ 6-7 and one Y105-dropout at z ~ 8 in a total survey area of 1.43 arcmin^2 in the source plane. The predictions of our lens model allow us to also identify five multiply-imaged systems lying at redshifts between z ~ 6 and z ~ 8. Thanks to constraints from the mass distribution in the cluster, we were able to estimate the effective survey volume corrected for completeness and magnification effects. This was in turn used to estimate the rest-frame ultraviolet luminosity function (LF) at z ~ 6-8. Our LF results are generally in agreement with the most recent blank field estimates, confirming the feasibility of surveys through lensing clusters. Although based on a shallower observations than what will be achieved in the final dataset including the full ACS observations, the LF presented here extends down to Muv ~ -18.5 at z ~ 7 with one identified object at Muv ~ -15 thanks to the highly-magnified survey areas. This early study forecasts the power of using massive galaxy clusters as cosmic telescopes and its complementarity to blank fields.

We constrain the total accreted mass density in supermassive black holes at z>6, inferred via the upper limit derived from the integrated X-ray emission from a sample of photometrically selected galaxy candidates. Studying galaxies obtained from the deepest Hubble Space Telescope images combined with the Chandra 4 Msec observations of the Chandra Deep Field South, we achieve the most restrictive constraints on total black hole growth in the early Universe. We estimate an accreted mass density <1000Mo Mpc^-3 at z~6, significantly lower than the previous predictions from some existing models of early black hole growth and earlier prior observations. These results place interesting constraints on early black growth and mass assembly by accretion and imply one or more of the following: (1) only a fraction of the luminous galaxies at this epoch contain active black holes; (2) most black hole growth at early epochs happens in dusty and/or less massive - as yet undetected - host galaxies; (3) there is a significant fraction of low-z interlopers in the galaxy sample; (4) early black hole growth is radiatively inefficient, heavily obscured and/or is due to black hole mergers as opposed to accretion or (5) the bulk of the black hole growth occurs at late times. All of these possibilities have important implications for our understanding of high redshift seed formation models.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

We predict the existence and observational signatures of a new class of objects that assembled early, during the first billion years of cosmic time: Obese Black-hole Galaxies (OBGs). OBGs are objects in which the mass of the central black hole initially exceeds that of the stellar component of the host galaxy, and the luminosity from black-hole accretion dominates the starlight. From a cosmological simulation, we demonstrate that there are sites where star formation is initially inhibited and direct-collapse black holes (DCBHs) form due to the photo-dissociating effect of Lyman-Werner radiation on molecular hydrogen. We show that the formation of OBGs is inevitable, because the probability of finding the required extra-galactic environment and the right physical conditions in a halo conducive to DCBH formation is quite high in the early universe. We estimate an OBG number density of 0.009/Mpc^3 at z~8 and 0.03/Mpc^3 at z~6. Extrapolating from our simulation volume, we infer that the most luminous quasars detected at z~6 likely transited through an earlier OBG phase. We find that these primordial galaxies start off with an over-massive BH and acquire their stellar component from subsequent merging as well as in-situ star formation. In doing so, they inevitably go through an OBG phase dominated by the accretion luminosity at the Eddington rate or below, released from the growing BH. The OBG phase is characterised by an ultra-violet (UV) spectrum with slope of beta ~ -2.3 and the absence of a Balmer Break. OBGs should also be spatially unresolved, and are expected to be brighter than the majority of known high-redshift galaxies. OBGs could potentially be revealed via HST follow-up imaging of samples of brighter Lyman-break galaxies provided by wide-area ground-based surveys such as UltraVISTA, and should be easily uncovered and studied with instruments aboard JWST...(abridged)

We present late-time Hubble Space Telescope imaging of the fields of six Swift GRBs lying at 5.0<z<9.5. Our data includes very deep observations of the field of the most distant spectroscopically confirmed burst, GRB 090423, at z=8.2. Using the precise positions afforded by their afterglows we can place stringent limits on the luminosities of their host galaxies. In one case, that of GRB 060522 at z=5.11, there is a marginal excess of flux close to the GRB position which may be a detection of a host at a magnitude J(AB)=28.5. None of the others are significantly detected meaning that all the hosts lie below L\star at their respective redshifts, with star formation rates SFR<4Mo/yr in all cases. Indeed, stacking the five fields with WFC3-IR data we conclude a mean SFR<0.17Mo/yr per galaxy. These results support the proposition that the bulk of star formation, and hence integrated UV luminosity, at high redshifts arises in galaxies below the detection limits of deep-field observations. Making the reasonable assumption that GRB rate is proportional to UV luminosity at early times allows us to compare our limits with expectations based on galaxy luminosity functions derived from the Hubble Ultra-Deep Field (HUDF) and other deep fields. We infer that a luminosity function which is evolving rapidly towards steeper faint-end slope (alpha) and decreasing characteristic luminosity (L\star), as suggested by some other studies, is consistent with our observations, whereas a non-evolving LF shape is ruled out at >90% confidence. Although it is not yet possible to make stronger statements, in the future, with larger samples and a fuller understanding of the conditions required for GRB production, studies like this hold great potential for probing the nature of star formation, the shape of the galaxy luminosity function, and the supply of ionizing photons in the early universe.

In this paper, we compare the observationally derived black hole mass function (BHMF) of luminous broad-line quasars (BLQSOs) at 1<z<4.5 drawn from the Sloan Digital Sky Survey (SDSS) presented in Kelly et al. (2010), with models of merger driven BH growth in the context of standard hierarchical structure formation models. In the models, we explore two distinct black hole seeding prescriptions at the highest redshifts: "light seeds" - remnants of Population III stars and "massive seeds" that form from the direct collapse of pre-galactic disks. The subsequent merger triggered mass build-up of the black hole population is tracked over cosmic time under the assumption of a fixed accretion rate as well as rates drawn from the distribution derived by Merloni & Heinz. Our model snapshots are compared to the SDSS derived BHMFs of BLQSOs. Our key findings are that the duty cycle of SMBHs powering BLQSOs increases with increasing redshift for all models and models with Pop III remnants as black hole seeds are unable to fit the observationally derived BHMFs for BLQSOs, lending strong support for the massive seeding model (abridged).

Clusters of galaxies are the most recently assembled, massive, bound structures in the Universe. As predicted by General Relativity, given their masses, clusters strongly deform space-time in their vicinity. Clusters act as some of the most powerful gravitational lenses in the Universe. Light rays traversing through clusters from distant sources are hence deflected, and the resulting images of these distant objects therefore appear distorted and magnified. Lensing by clusters occurs in two regimes, each with unique observational signatures. The strong lensing regime is characterized by effects readily seen by eye, namely, the production of giant arcs, multiple-images, and arclets. The weak lensing regime is characterized by small deformations in the shapes of background galaxies only detectable statistically. Cluster lenses have been exploited successfully to address several important current questions in cosmology: (i) the study of the lens(es) - understanding cluster mass distributions and issues pertaining to cluster formation and evolution, as well as constraining the nature of dark matter; (ii) the study of the lensed objects - probing the properties of the background lensed galaxy population - which is statistically at higher redshifts and of lower intrinsic luminosity thus enabling the probing of galaxy formation at the earliest times right up to the Dark Ages; and (iii) the study of the geometry of the Universe - as the strength of lensing depends on the ratios of angular diameter distances between the lens, source and observer, lens deflections are sensitive to the value of cosmological parameters and offer a powerful geometric tool to probe Dark Energy. In this review, we present the basics of cluster lensing and provide a current status report of the field.

In a recently published article, we quantified the impact of primordial non-Gaussianity on the probability of giant-arc formation. In that work, we focused on the local form of non-Gaussianity and found that it can have only a modest effect given the most recent constraints from Cosmic Microwave Background (CMB) measurements. Here, we present new calculations using a parameterization of scale-dependent non-Gaussianity in which the primordial bispectrum has the equilateral shape and the effective f_NL parameter depends on scale. We find that non-Gaussianity of this type can yield a larger effect on the giant-arc abundance compared to the local form due to both the scale dependence and the relatively weaker constraints on the equilateral shape from CMB measurements. In contrast to the maximum ~40% effect (within the latest CMB constraints) previously found for the local form, we find that the predicted giant-arc abundance for the scale-dependent equilateral form can differ by a factor of a few with respect to the Gaussian case.

We constrain gas and dark matter (DM) parameters of galaxy groups and clusters, by comparing X-ray scaling relations to theoretical expectations, obtained assuming that the gas is in hydrostatic equilibrium with the DM and follows a polytropic relation. We vary four parameters: the gas polytropic index Gamma, its temperature at large radii T_xi, the DM logarithmic slope at large radii zeta and its concentration c_vir. When comparing the model to the observed mass-temperature (M-T) relation of local clusters, our results are independent of both T_xi and c_vir. We thus obtain constraints on Gamma, by fixing the DM profile, and on zeta, by fixing the gas profile. For an NFW DM profile, we find that 6/5<Gamma<13/10, which is consistent with numerical simulations and observations of individual clusters. Taking 6/5<Gamma<13/10 allows the DM profile to be slightly steeper than the NFW profile at large radii. Upon including local groups, we constrain the mass-dependence of Gamma and the value of T_xi. Interestingly, with Gamma=6/5 and zeta=-3, we reproduce the observed steepening/breaking of the M-T relation at low M, if 10^6 K<T_xi<10^7 K, consistent with simulations and observations of the warm-hot intergalactic medium. When extrapolated to high redshift z, the model with a constant Gamma reproduces the expected self-similar behaviour. We also account for the observed, non-self-similar relations provided by some high-z clusters, as they provide constraints on the evolution of Gamma. Comparing our model to the observed luminosity-temperature relation, we discriminate between different M-c_vir relations: a weak dependence of c_vir on M is currently preferred by data. This simple theoretical model accounts for much of the complexity of recent, improved X-ray scaling relations, provided that we allow for a mild dependence of Gamma on M or for T_xi consistent with intercluster values. [abridged]

One of the key open questions in cosmology today pertains to understanding when, where and how super massive black holes form, while it is clear that mergers likely play a significant role in the growth cycles of black holes, how supermassive black holes form, and how galaxies grow around them. Here, we present Hubble Space Telescope WFC3/IR grism observations of a clumpy galaxy at z=1.35, with evidence for 10^6 - 10^7 Msun rapidly growing black holes in separate sub-components of the host galaxy. These black holes could have been brought into close proximity as a consequence of a rare multiple galaxy merger or they could have formed in situ. Such holes would eventually merge into a central black hole as the stellar clumps/components presumably coalesce to form a galaxy bulge. If we are witnessing the in-situ formation of multiple black holes, their properties can inform seed formation models and raise the possibility that massive black holes can continue to emerge in star-forming galaxies as late as z=1.35 (4.8 Gyr after the Big Bang).

We study the mass distribution of a sample of 28 galaxy clusters using strong and weak lensing observations. The clusters are selected via their strong lensing properties as part of the Sloan Giant Arcs Survey (SGAS) from the Sloan Digital Sky Survey (SDSS). Mass modelling of the strong lensing information from the giant arcs is combined with weak lensing measurements from deep Subaru/Suprime-cam images to primarily obtain robust constraints on the concentration parameter and the shape of the mass distribution. We find that the concentration c_vir is a steep function of the mass, c_vir \propto M_vir^-0.59\pm0.12, with the value roughly consistent with the lensing-bias-corrected theoretical expectation for high mass (10^15 h^-1 M_sun) clusters. However, the observationally inferred concentration parameters appear to be much higher at lower masses (10^14 h^-1 M_sun), possibly a consequence of the modification to the inner density profiles provided by baryon cooling. The steep mass-concentration relation is also supported from direct stacking analysis of the tangential shear profiles. In addition, we explore the two-dimensional shape of the projected mass distribution by stacking weak lensing shear maps of individual clusters with prior information on the position angle from strong lens modelling, and find significant evidence for a large mean ellipticity with the best-fit value of e = 0.47 \pm 0.06 for the mass distribution of the stacked sample. We find that the luminous cluster member galaxy distribution traces the overall mass distribution very well, although the distribution of fainter cluster galaxies appears to be more extended than the total mass.

The formation of the first massive objects in the infant Universe remains impossible to observe directly and yet it sets the stage for the subsequent evolution of galaxies. While some black holes with masses > billion solar masses? have been detected in luminous quasars less than one billion years after the Big Bang, these individual extreme objects have limited utility in constraining the channels of formation of the earliest black holes. The initial conditions of black hole seed properties are quickly erased during the growth process. From deep, optimally stacked, archival X-ray observations, we measure the amount of black hole growth in z=6-8 galaxies (0.7-1 billion years after the Big Bang). Our results imply that black holes grow in tandem with their hosts throughout cosmic history, starting from the earliest times. We find that most copiously accreting black holes at these epochs are buried in significant amounts of gas and dust that absorb most radiation except for the highest energy X-rays. This suggests that black holes grow significantly more than previously thought during these early bursts, and due to obscuration they do not contribute to the re-ionization of the Universe with their ultraviolet emission.

We explore the use of strong lensing by galaxy clusters to constrain the dark energy equation of state and its possible time variation. The cores of massive clusters often contain several multiply imaged systems of background galaxies at different redshifts. The locations of lensed images can be used to constrain cosmological parameters due to their dependence on the ratio of angular diameter distances. We employ Monte-Carlo simulations of cluster lenses, including the contribution from substructures, to assess the feasibility of this potentially powerful technique. At the present, parametric lens models use well motivated scaling relations between mass and light to incorporate cluster member galaxies, and do not explicitly model line-of-sight structure. Here, we quantify modeling errors due to scatter in the cluster galaxy scaling relations and un-modeled line-of-sight halos. These errors are of the order of a few arcseconds on average for clusters located at typical redshifts (z ~ 0.2 - 0.3). Using Bayesian Markov Chain Monte-Carlo techniques, we show that the inclusion of these modeling errors is critical to deriving unbiased constraints on dark energy. However, when the uncertainties are properly quantified, we show that constraints competitive with other methods may be obtained by combining results from a sample of just 10 simulated clusters with 20 families each. Cosmography with a set of well studied cluster lenses may provide a powerful complementary probe of the dark energy equation of state. Our simulations provide a convenient method of quantifying modeling errors and assessing future strong lensing survey strategies.

We present results of a study of the virial state of high redshift dark matter haloes in an N-body simulation. We find that the majority of collapsed, bound haloes are not virialized at any redshift slice in our study ($z=15-6$) and have excess kinetic energy. At these redshifts, merging is still rampant and the haloes cannot strictly be treated as isolated systems. To assess if this excess kinetic energy arises from the environment, we include the surface pressure term in the virial equation explicitly and relax the assumption that the density at the halo boundary is zero. Upon inclusion of the surface term, we find that the haloes are much closer to virialization, however, they still have some excess kinetic energy. We report trends of the virial ratio including the extra surface term with three key halo properties: spin, environment, and concentration. We find that haloes with closer neighbors depart more from virialization, and that haloes with larger spin parameters do as well. We conclude that except at the lowest masses ($M < 10^6 \Msun$), dark matter haloes at high redshift are not fully virialized. This finding has interesting implications for galaxy formation at these high redshifts, as the excess kinetic energy will impact the subsequent collapse of baryons and the formation of the first disks and/or baryonic structures.

In this paper, we examine if the properties of central black holes in galactic nuclei correlate with their host dark matter halo. We analyze the entire sample of galaxies where black hole mass, velocity dispersion, sigma, and asymptotic circular velocity, Vc, have all been measured. We fit M-sigma and M-Vc to a power law, and find that in both relationships the scatter and slope are similar. This model-independent analysis suggests that although the black hole masses are not uniquely determined by dark matter halo mass, when considered for the current sample as a whole, the M-Vc correlation may be as strong (or as weak) as M-sigma. Although the data are sparse, there appears to be more scatter in the correlation for both sigma and Vc at the low--mass end. This is not unexpected given our current understanding of galaxy and black hole assembly. In fact, there are several compelling reasons that account for this: (i) SMBH formation is likely less efficient in low-mass galaxies with large angular momentum content; (ii) SMBH growth is less efficient in low-mass disk galaxies that have not experienced major mergers; (iii) dynamical effects, such as gravitational recoil, increase scatter preferentially at the low-mass end. Therefore, the recent observational claim of the absence of central SMBHs in bulgeless, low mass galaxies, or deviations from the correlations defined by high-mass black holes in large galaxies today is, in fact, predicated by current models of black hole growth. We show how this arises as a direct consequence of the coupling between dark matter halos and central black holes at the earliest epochs.

We track the growth and evolution of high redshift seed black holes over cosmic time. This population of massive, initial black hole seeds form at these early epochs from the direct collapse of pre-galactic gas discs. Populating dark matter halos with seeds formed in this fashion, we follow their mass assembly history to the present time using a Monte-Carlo merger tree approach. Using this formalism, we predict the black hole mass function at the present time; the integrated mass density of black holes in the Universe; the luminosity function of accreting black holes as a function of redshift and the scatter in observed, local Mbh{\sigma}s relation. Signatures of these massive seed models appear predominantly at the low mass end of the present day black hole mass function. In fact, our prediction of the shape of the Mbh{\sigma}s relation at the low mass end and increased scatter has recently been corroborated by observations. These models predict that low surface brightness, bulge-less galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts also has a direct influence on the black hole occupation fraction in galaxies at z = 0. This effect is more pronounced for low mass galaxies today as we predict the existence of a population of low mass galaxies that do not host nuclear black holes. This is the key discriminant between the models studied here and the Population-III remnant seed model.

Tracking the evolution of high redshift seed black hole masses to late times, we examine the observable signatures today. These massive initial black hole seeds form at extremely high redshifts from the direct collapse of pre-galactic gas discs. Populating dark matter halos with seeds formed in this fashion, we follow the mass assembly history of these black holes to the present time using a Monte-Carlo merger tree approach. Utilizing this formalism, we predict the black hole mass function at high redshifts and at the present time; the integrated mass density of black holes in the Universe; the luminosity function of accreting black holes as a function of redshift and the scatter in observed, local M-sigma relation. Comparing the predictions of the `light' seed model with these massive seeds we find that significant differences appear predominantly at the low mass end of the present day black hole mass function. However, all our models predict that low surface brightness, bulge-less galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts has a direct influence on the black hole occupation fraction in galaxies at z=0. This effect is more pronounced for low mass galaxies. This is the key discriminant between the models studied here and the Population III remnant `light' seed model. We find that there exists a population of low mass galaxies that do not host nuclear black holes. Our prediction of the shape of the M-sigma relation at the low mass end and increased scatter has recently been corroborated by observations.

For over a decade, it has been debated whether the concordance LCDM model is consistent with the observed abundance of giant arcs in clusters. While previous theoretical studies have focused on properties of the lens and source populations, as well as cosmological effects such as dark energy, the impact of initial conditions on the giant-arc abundance is relatively unexplored. Here, we quantify the impact of non-Gaussian initial conditions with the local bispectrum shape on the predicted frequency of giant arcs. Using a path-integral formulation of the excursion set formalism, we extend a semi-analytic model for calculating halo concentrations to the case of primordial non-Gaussianity, which may be useful for applications outside of this work. We find that massive halos tend to collapse earlier in models with positive f_NL, relative to the Gaussian case, leading to enhanced concentration parameters. The converse is true for f_NL < 0. In addition to these effects, which change the lensing cross sections, non-Gaussianity also modifies the abundance of supercritical clusters available for lensing. These combined effects work together to either enhance (f_NL > 0) or suppress (f_NL < 0) the probability of giant-arc formation. Using the best value and 95% confidence levels currently available from the Wilkinson Microwave Anisotropy Probe, we find that the giant-arc optical depth for sources at z_s~2 is enhanced by ~20% and ~45% for f_NL = 32 and 74 respectively. In contrast, we calculate a suppression of ~5% for f_NL = -10. These differences translate to similar relative changes in the predicted all-sky number of giant arcs.

Galaxy formation is significantly modulated by energy output from supermassive black holes at the centers of galaxies which grow in highly efficient luminous quasar phases. The timescale on which black holes transition into and out of such phases is, however, unknown. We present the first measurement of the shutdown timescale for an individual quasar using X-ray observations of the nearby galaxy IC 2497, which hosted a luminous quasar no more than 70,000 years ago that is still seen as a light echo in `Hanny's Voorwerp', but whose present-day radiative output is lower by at least 2 and more likely by over 4 orders of magnitude. This extremely rapid shutdown provides new insights into the physics of accretion in supermassive black holes, and may signal a transition of the accretion disk to a radiatively inefficient state.

Current efforts in observational cosmology are focused on characterizing the mass-energy content of the Universe. We present results from a geometric test based on strong lensing in galaxy clusters. Based on Hubble Space Telescope images and extensive ground-based spectroscopic follow-up of the massive galaxy cluster Abell 1689, we used a parametric model to simultaneously constrain the cluster mass distribution and dark energy equation of state. Combining our cosmological constraints with those from X-ray clusters and the Wilkinson Microwave Anisotropy Probe 5-year data gives {\Omega}m = 0.25 +/- 0.05 and wx = -0.97 +/- 0.07 which are consistent with results from other methods. Inclusion of our method with all other techniques available brings down the current 2{\sigma} contours on the dark energy equation of state parameter wx by about 30%.

In this paper, we examine in detail the key structural properties of high redshift dark matter haloes as a function of their spin parameter. We perform and analyze high resolution cosmological simulations of the formation of structure in a LCDM Universe. We study the mass function, ellipticities, shapes, density profiles, rotation curves and virialization for a large sample of dark matter haloes from z = 15 - 6. We also present detailed convergence tests for individual haloes. We find that high spin haloes have stronger clustering strengths (up to 25%) at all mass and redshift ranges at these early epochs. High redshift spherical haloes are also up to 50% more clustered than aspherical haloes. High spin haloes at these redshifts are also preferentially found in high density environments, and have more neighbors than their low spin counterparts. We report a systematic offset in the peak of the circular velocity curves for high and low spin haloes of the same mass. Therefore, estimating halo masses without knowledge of the spin, using only the circular velocity can yield errors of up to 40%. The strong dependence of key structural properties on spin that we report here likely have important implications for studies of star formation and feedback from these galaxies.

We present an analytic formulation for the equilibrium gas density profile of early-type galaxies that explicitly includes the contribution of stars in the gravitational potential. We build a realistic model for an isolated elliptical galaxy and explore the equilibrium gas configurations as a function of multiple parameters. For an assumed central gas temperature k_B*T_0=0.6 keV, we find that neglecting the gravitational effects of stars, which can contribute substantially in the innermost regions, leads to an underestimate of the enclosed baryonic gas mass by up to ~65% at the effective radius, and by up to ~15% at the NFW scale radius, depending on the stellar baryon fraction. This formula is therefore important for estimating the baryon fraction in an unbiased fashion. These new hydrostatic equilibrium solutions, derived for the isothermal and polytropic cases, can also be used to generate more realistic initial conditions for simulations of elliptical galaxies. Moreover, the new formulation is relevant when interpreting X-ray data. We compare our composite isothermal model to the standard beta-model used to fit X-ray observations of early-type galaxies, to determine the value of the NFW scale radius r_s. Assuming a 10% stellar baryon fraction, we find that the exclusion of stars from the gravitational potential leads to (i) an underestimate of r_s by ~80%, and to (ii) an overestimate of the enclosed dark matter at r_s by a factor of ~2, compared to the equivalent beta-model fit results when stars are not taken into account. For higher stellar mass fractions, a beta-model is unable to accurately reproduce our solution, indicating that when the observed surface brightness profile of an isolated elliptical galaxy is found to be well fitted by a beta-model, the stellar mass fraction cannot be much greater than ~10%.

Scaling relations of observed galaxy cluster properties are useful tools for constraining cosmological parameters as well as cluster formation histories. One of the key cosmological parameters, sigma8, is constrained using observed clusters of galaxies, although current estimates of sigma8 from the scaling relations of dynamically relaxed galaxy clusters are limited by the large scatter in the observed cluster mass-temperature (M-T) relation. With a sample of eight strong lensing clusters at 0.3 < z <0.8, we find that the observed cluster concentration-mass relation can be used to reduce the M-T scatter by a factor of 6. Typically only relaxed clusters are used to estimate sigma8, but combining the cluster concentration-mass relation with the M-T relation enables the inclusion of unrelaxed clusters as well. Thus, the resultant gains in the accuracy of sigma8 measurements from clusters are twofold: the errors on sigma8 are reduced and the cluster sample size is increased. Therefore, the statistics on sigma8 determination from clusters are greatly improved by the inclusion of unrelaxed clusters. Exploring cluster scaling relations further, we find that the correlation between brightest cluster galaxy (BCG) luminosity and cluster mass offers insight into the assembly histories of clusters. We find preliminary evidence for a steeper BCG luminosity - cluster mass relation for strong lensing clusters than the general cluster population, hinting that strong lensing clusters may have had more active merging histories.

Despite observed strong correlations between central supermassive black holes (SMBHs) and star-formation in galactic nuclei, uncertainties exist in our understanding of their coupling. We present observations of the ratio of heavily-obscured to unobscured quasars as a function of cosmic epoch up to z~3, and show that a simple physical model describing mergers of massive, gas-rich galaxies matches these observations. In the context of this model, every obscured and unobscured quasar represent two distinct phases that result from a massive galaxy merger event. Much of the mass growth of the SMBH occurs during the heavily-obscured phase. These observations provide additional evidence for a causal link between gas-rich galaxy mergers, accretion onto the nuclear SMBH and coeval star formation.

In this paper, we explore the establishment and evolution of the empirical correlation between black hole mass and velocity dispersion with redshift. We track the growth and accretion history of massive black holes starting from high redshift using two seeding models:(i) Population III remnants, and (ii) massive seeds from direct gas collapse. Although the seeds do not initially satisfy the $M_{\rm BH} - \sigma$ relation, the correlation is established and maintained at all times if self-regulating accretion episodes are associated with major mergers. The massive end of the $M_{\rm BH} - \sigma$ relation is established early, and lower mass MBHs migrate over time. How MBHs migrate toward the relation, the slope and the scatter of the relation all depend critically on the seeding model as well as the adopted self-regulation prescription. We expect flux limited AGN surveys and LISA to select accreting and merging MBHs respectively that have already migrated onto the $\msigma$ relation. This is a consequence of major mergers being more common at high redshift for the most massive, biased, galaxies that anchor the $\msigma$ relation early. We also predict the existence of a large population of low mass `hidden' MBHs at high redshift which can easily escape detection. Additionally, we find that if MBH seeds are massive, $\sim 10^5 M_{\odot}$, the low-mass end of the $\msigma$ flattens towards this asymptotic value, creating a characteristic `plume'.

By stacking an ensemble of strong lensing clusters, we demonstrate the feasibility of placing constraints on the dark energy equation of state. This is achieved by using multiple images of sources at two or more distinct redshift planes. The sample of smooth clusters in our simulations is based on observations of massive clusters and the distribution of background galaxies is constructed using the Hubble Deep Field. Our source distribution reproduces the observed redshift distribution of multiply imaged sources in Abell 1689. The cosmology recovery depends on the number of image families with known spectroscopic redshifts and the number of stacked clusters. Our simulations suggest that constraints comparable to those derived from other competing established techniques on a constant dark energy equation of state can be obtained using 10 to 40 clusters with 5 or more families of multiple images. We have also studied the observational errors in the image redshifts and positions. We find that spectroscopic redshifts and high resolution {\it Hubble Space Telescope} images are required to eliminate confidence contour relaxation relative to the ideal case in our simulations. This suggests that the dark energy equation of state, and other cosmological parameters, can be constrained with existing {\it Hubble Space Telescope} images of lensing clusters coupled with dedicated ground-based arc spectroscopy.

We derive radial mass profiles of four strong lensing selected clusters which show prominent giant arcs (Abell 1703, SDSS J1446+3032, SDSS J1531+3414, and SDSS J2111-0115), by combining detailed strong lens modeling with weak lensing shear measured from deep Subaru Suprime-cam images. Weak lensing signals are detected at high significance for all four clusters, whose redshifts range from z=0.28 to 0.64. We demonstrate that adding strong lensing information with known arc redshifts significantly improves constraints on the mass density profile, compared to those obtained from weak lensing alone. While the mass profiles are well fitted by the universal form predicted in N-body simulations of the LCDM model, all four clusters appear to be slightly more centrally concentrated (the concentration parameters c_vir \sim 8) than theoretical predictions, even after accounting for the bias toward higher concentrations inherent in lensing selected samples. Our results are consistent with previous studies which similarly detected a concentration excess, and increases the total number of clusters studied with the combined strong and weak lensing technique to ten. Combining our sample with previous work, we find that clusters with larger Einstein radii are more anomalously concentrated. We also present a detailed model of the lensing cluster A1703 with constraints from multiple image families, and find the dark matter inner density profile to be cuspy with the slope consistent with -1, in agreement with expectations.

We analyze high resolution, N-body hydrodynamical simulations of fiducial galaxy clusters to probe tidal stripping of the dark matter subhalos. These simulations include a prescription for star formation allowing us to track the fate of the stellar component as well. We investigate the effect of tidal stripping on cluster galaxies hosted in these dark matter subhalos as a function of cluster-centric radius. To quantify the extent of the dark matter halos of cluster galaxies, we introduce the half mass radius r_half as a diagnostic, and study its evolution with projected cluster-centric distance R as a function of redshift. We find a well defined trend for (r_half,R): the closer the galaxies are to the center of the cluster, the smaller the half mass radius. Interestingly, this trend is inferred in all redshift frames examined in this work ranging from z=0 to z=0.7. At z=0, galaxy halos in the central regions of clusters are found to be highly truncated, with the most compact half mass radius of 10 kpc. We also find that r_half depends on luminosity and we present scaling relations of r_half with galaxy luminosity. The corresponding total mass of the cluster galaxies is also found to increase with projected cluster-centric distance and luminosity, but with more scatter than the (r_half,R) trend. Comparing the distribution of stellar mass to total mass for cluster galaxies, we find that the dark matter component is preferentially stripped, whereas the stellar component remains protected by the halo and is much less affected by tidal forces. We compare these results with galaxy-galaxy lensing probes of r_half and find qualitative agreement. (Abridged)

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

Weak gravitational lensing provides a potentially powerful method for the detection of clusters. In addition to cluster candidates, a large number of objects with possibly no optical or X-ray component have been detected in shear-selected samples. We develop an analytic model to investigate the claim of Weinberg & Kamionkowski (2002) that unvirialised protoclusters account for a significant number of these so-called "dark" lenses. In our model, a protocluster consists of a small virialised region surrounded by in-falling matter. We find that, in order for a protocluster to simultaneously escape X-ray detection and create a detectable weak lensing signal, it must have a small virial mass (~10^{13} \Msun) and large total mass (~ 10^{15} \Msun), with a relatively flat density profile outside of the virial radius. Such objects would be characterized by rising tangential shear profiles well beyond the virial radius. We use a semi-analytic approach based on the excursion set formalism to estimate the abundance of lensing protoclusters with a low probability of X-ray detection. We find that they are extremely rare, accounting for less than 0.4 per cent of the total lenses in a survey with background galaxy density n = 30 arcmin^{-2} and an intrinsic ellipticity dispersion of 0.3. We conclude that lensing protoclusters with undetectable X-Ray luminosities are too rare to account for a significant number of dark lenses.

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

Theories of structure formation in a cold dark matter dominated Universe predict that massive clusters of galaxies assemble from the hierarchical merging of lower mass subhalos. Exploiting strong and weak gravitational lensing signals inferred from panoramic HST imaging data, we present a high resolution reconstruction of the mass distribution in the massive, lensing cluster Cl0024+16. Applying galaxy-galaxy lensing techniques we track the fate of dark matter subhalos as a function of projected cluster-centric radius out to 5 Mpc, well beyond the virial radius. We report the first detection of the statistical lensing signal of dark matter subhalos associated with late-type galaxies in clusters. The mass of a typical dark matter subhalo that hosts an L* galaxy increases with projected cluster-centric radius in line with expectations from the tidal stripping hypothesis. Early-type galaxies appear to be hosted on average in more massive dark matter subhalos compared to late-type galaxies. We interpret our findings as evidence for the active assembly of mass via tidal stripping in galaxy clusters. The mass function of dark matter subhalos as a function of projected cluster-centric radius, is compared with an equivalent mass function derived from clusters in the Millenium Run simulation populated with galaxies using semi-analytic models. The shape of the lensing determined mass functions are in very good agreement with the subhalo mass functions derived from the Millenium Run simulation. However, simulated subhalos appear to be more efficiently stripped than lensing observations suggest. (Abridged)

We make a case for the existence for ultra-massive black holes (UMBHs) in the Universe, but argue that there exists a likely upper limit to black hole masses of the order of $M \sim 10^{10} \msun$. We show that there are three strong lines of argument that predicate the existence of UMBHs: (i) expected as a natural extension of the observed black hole mass bulge luminosity relation, when extrapolated to the bulge luminosities of bright central galaxies in clusters; (ii) new predictions for the mass function of seed black holes at high redshifts predict that growth via accretion or merger-induced accretion inevitably leads to the existence of rare UMBHs at late times; (iii) the local mass function of black holes computed from the observed X-ray luminosity functions of active galactic nuclei predict the existence of a high mass tail in the black hole mass function at $z = 0$. Consistency between the optical and X-ray census of the local black hole mass function requires an upper limit to black hole masses. This consistent picture also predicts that the slope of the $M_{\rm bh}$-$\sigma$ relation will evolve with redshift at the high mass end. Models of self-regulation that explain the co-evolution of the stellar component and nuclear black holes naturally provide such an upper limit. The combination of multi-wavelength constraints predicts the existence of UMBHs and simultaneously provides an upper limit to their masses. The typical hosts for these local UMBHs are likely the bright, central cluster galaxies in the nearby Universe.

Deep 1.1 mm continuum observations of 1E0657-56 (the "Bullet Cluster") taken with the millimeter-wavelength camera AzTEC on the 10-m Atacama Submillimeter Telescope Experiment (ASTE), have revealed an extremely bright (S$_{\rm{1.1mm}}=15.9$ mJy) unresolved source. This source, MMJ065837-5557.0, lies close to a maximum in the density of underlying mass-distribution, towards the larger of the two interacting clusters as traced by the weak-lensing analysis of Clowe et al. 2006. Using optical--IR colours we argue that MMJ065837-5557.0 lies at a redshift of $z = 2.7 \pm 0.2$. A lensing-derived mass-model for the Bullet Cluster shows a critical-line (caustic) of magnification within a few arcsecs of the AzTEC source, sufficient to amplify the intrinsic millimeter-wavelength flux of the AzTEC galaxy by a factor of $\gg 20$. After subtraction of the foreground cluster emission at 1.1mm due to the Sunyaev-Zel'dovich effect, and correcting for the magnification, the rest-frame FIR luminosity of MMJ065837-5557.0 is $\le 10^{12} \rm L_{\odot}$, characteristic of a luminous infrared galaxy (LIRG). We explore various scenarios to explain the colors, morphologies and positional offsets between the potential optical and IR counterparts, and their relationship with MMJ065837-5557.0. Until higher-resolution and more sensitive (sub)millimeter observations are available, the detection of background galaxies close to the caustics of massive lensing clusters offers the only opportunity to study this intrinsically faint millimeter-galaxy population.

Clusters of galaxies offer a robust test bed for probing the nature of dark matter that is insensitive to the assumption of the gravity theories. Both Modified Newtonian Dynamics (MOND) and General Relativity (GR) would require similar amounts of non-baryonic matter in clusters as MOND boosts the gravity only mildly on cluster scales. Gravitational lensing allows us to estimate the enclosed mass in clusters on small (20 - 50 kpc) and large (several 100 kpc) scales independent of the assumptions of equilibrium. Here we show for the first time that a combination of strong and weak gravitational lensing effects can set interesting limits on the phase space density of dark matter in the centers of clusters. The phase space densities derived from lensing observations are inconsistent with neutrino masses ranging from 2 - 7 eV, and hence do not support the 2 eV-range particles required by MOND. To survive, plausible modifications for MOND may be either an additional degree of dynamical freedom in a co-variant incarnation or mass-varying theories of neutrinos.

(abridged) The Observatory for Multi-Epoch Gravitational Lens Astrophysics (OMEGA) is a mission concept for a 1.5-m near-UV through near-IR space observatory that will be dedicated to frequent imaging and spectroscopic monitoring of ~100 multiply-imaged active galactic nuclei over the whole sky. Using wavelength-tailored dichroics with extremely high transmittance, efficient imaging in six channels will be done simultaneously during each visit to each target. The separate spectroscopic mode, engaged through a flip-in mirror, uses an image slicer spectrograph. After a period of many visits to all targets, the resulting multidimensional movies can then be analyzed to a) measure the mass function of dark matter substructure; b) measure precise masses of the accreting black holes as well as the structure of their accretion disks and their environments over several decades of physical scale; and c) measure a combination of the local Hubble expansion and cosmological distances to unprecedented precision.

Galaxy formation requires a process that continually heats gas and quenches star formation in order to reproduce the observed shape of the luminosity function of bright galaxies. To accomplish this, current models invoke heating from supernovae, and energy injection from active galactic nuclei. However, observations of radio-loud active galactic nuclei suggest that their feedback is likely to not be as efficient as required, signaling the need for additional heating processes. We propose the self-annihilation of weakly interacting massive particles that constitute dark matter as a steady source of heating. In this paper, we explore the circumstances under which this process may provide the required energy input. To do so, dark matter annihilations are incorporated into a galaxy formation model within the Millennium cosmological simulation. Energy input from self-annihilation can compensate for all the required gas cooling and reproduce the observed galaxy luminosity function only for what appear to be extreme values of the relevant key parameters. The key parameters are: the slope of the inner density profile of dark matter haloes and the outer spike radius. The inner density profile needs to be steepened to slopes of -1.5 or more and the outer spike radius needs to extend to a few tens of parsecs on galaxy scales and a kpc or so on cluster scales. If neutralinos or any thermal relic WIMP with s-wave annihilation constitute dark matter, their self-annihilation is inevitable and could provide enough power to modulate galaxy formation. Energy from self-annihilating WIMPs could be yet another piece of the feedback puzzle along with supernovae and active galactic nuclei.

The recently revised cosmological constraints from the Five-Year WMAP data ameliorate previous tension between cosmological constraints from the microwave background and from cluster abundances. We demonstrate that the revised estimates of cosmological parameters are in excellent agreement with the mass function of X-ray clusters in the Sloan Digital Sky Survey. Velocity segregation between galaxies and the underlying dark matter could cause virial mass estimates to be biased, causing the mass scale of the mass function to be offset from the true value. Modest velocity segregation ($\sigma_{gxy}/\sigma_{DM}$=1.13$^{+0.06}_{-0.05}$) is sufficient to match the mass function to the Five-Year WMAP results. When the new WMAP results are combined with constraints from supernovae and baryon acoustic oscillations, there is no need for velocity segregation ($\sigma_{gxy}/\sigma_{DM}$=1.05$\pm$0.05). This result agrees with expectations for velocity segregation from state-of-the-art numerical simulations of clusters. Together with the improved agreement between the new WMAP results and recent cosmic shear measurements, this result demonstrates that the amplitude of large-scale structure in the nearby universe matches that predicted from the structure seen in the microwave background. The new constraint we place on velocity segregation in clusters indicates that virial mass estimates for clusters are reasonably accurate. This result suggests that future cluster surveys will be able to probe both cosmological parameters and fundamental cluster physics.

We present a reconstruction of the mass distribution of galaxy cluster Abell 1689 at z = 0.18 using detected strong lensing features from deep HST/ACS observations and extensive ground based spectroscopy. Earlier analyses have reported up to 32 multiply imaged systems in this cluster, of which only 3 were spectroscopically confirmed. In this work, we present a parametric strong lensing mass reconstruction using 24 multiply imaged systems with newly determined spectroscopic redshifts, which is a major step forward in building a robust mass model. In turn, the new spectroscopic data allows a more secure identification of multiply imaged systems. The resultant mass model enables us to reliably predict the redshifts of additional multiply imaged systems for which no spectra are currently available, and to use the location of these systems to further constrain the mass model. In particular, we have detected 5 strong galaxy-galaxy lensing systems just outside the Einstein ring region, further constraining the mass profile. Our strong lensing mass model is consistent with that inferred from our large scale weak lensing analysis derived using CFH12k wide field images. Thanks to a new method for reliably selecting a well defined background lensed galaxy population, we resolve the discrepancy found between the strong and weak lensing mass models reported in earlier work. [ABRIDGED]

We present a parametric strong lensing model of the cluster Abell 2218 based on HST ACS data. We constrain the lens model using 8 previously known multiply imaged systems, of which 7 have spectroscopically confirmed redshifts. In addition, we propose five candidate multiply imaged systems and estimate their redshifts using our mass model. The model parameters are optimized in the source plane by a bayesian Monte Carlo Markov Chain as implemented in the the publicly available software Lenstool. We find rms=0."12 for the scatter of the sources in the source plane, which translates into rms=1."49 between the predicted and measured image positions in the image plane. We find that the projected mass distribution of Abell 2218 is bimodal, which is supported by an analysis of the light distribution. We also find evidence for two structures in velocity space, separated by ~1000 km/s, corresponding to the two large scale dark matter clumps. We find that the lensing constraints can not be well reproduced using only dark matter halos associated with the cluster galaxies, but that the dark matter is required to be smoothly distributed in large scale halos. At 100" the enclosed projected mass is 3.8e14 solar masses. At that radius, the large scale halos contribute ~85% of the mass. We find that the model is not very sensitive to the fainter galaxy sized halos, unless they locally perturb a given multiply imaged system. Therefore, dark galaxy sized substructure can be reliably constrained only if it locally perturbs one of the systems. In an appendix we give a self-contained description of the parametric profile we use, the dual pseudo isothermal elliptical mass distribution (dPIE). This profile is a two component pseudo isothermal mass distribution (PIEMD) with both a core radii and a scale radii. (Abridged)

We investigate the evolution of high redshift seed black hole masses at late times and their observational signatures. The massive black hole seeds studied here form at extremely high redshifts from the direct collapse of pre-galactic gas discs. Populating dark matter halos with seeds formed in this way, we follow the mass assembly of these black holes to the present time using a Monte-Carlo merger tree. Using this machinery we predict the black hole mass function at high redshifts and at the present time; the integrated mass density of black holes and the luminosity function of accreting black holes as a function of redshift. These predictions are made for a set of three seed models with varying black hole formation efficiency. Given the accuracy of current observational constraints, all 3 models can be adequately fit. Discrimination between the models appears predominantly at the low mass end of the present day black hole mass function which is not observationally well constrained. However, all our models predict that low surface brightness, bulgeless galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts has a direct influence on the black hole occupation fraction in galaxies at z=0. This effect is more pronounced for low mass galaxies. This is the key discriminant between the models studied here and the Population III remnant seed model. We find that there exists a population of low mass galaxies that do not host nuclear black holes. Our prediction of the shape of the black hole mass - velocity dispersion relation at the low mass end is in agreement with the recent observational determination from the census of low mass galaxies in the Virgo cluster.

The redshift distribution of galaxy lenses in known gravitational lens systems provides a powerful test that can potentially discriminate amongst cosmological models. However, applications of this elegant test have been curtailed by two factors: our ignorance of how galaxies evolve with redshift, and the absence of methods to deal with the effect of incomplete information in lensing systems. In this paper, we investigate both issues in detail. We explore how to extract the properties of evolving galaxies, assuming that the cosmology is well determined by other techniques. We propose a new nested Monte Carlo method to quantify the effects of incomplete data. We apply the lens-redshift test to an improved sample of seventy lens systems derived from recent observations, primarily from the SDSS, SLACS and the CLASS surveys. We find that the limiting factor in applying the lens-redshift test derives from poor statistics, including incomplete information samples, and biased sampling. Many lenses that uniformly sample the underlying true image separation distribution will be needed to use this test as a complementary method to measure the value of the cosmological constant or the properties of evolving galaxies. Planned future surveys by missions like the SNAP satellite or LSST are likely to usher in a new era for strong lensing studies that utilize this test. With expected catalogues of thousands of new strong lenses, the lens-redshift test could offer a powerful tool to probe cosmology as well as galaxy evolution.

The properties of clusters of galaxies offer key insights into the assembly process of structure in the universe. Numerical simulations of cosmic structure formation in a hierarchical, dark matter dominated universe suggest that galaxy cluster concentrations, which are a measure of a halo's central density, decrease gradually with virial mass. However, cluster observations have yet to confirm this correlation. The slopes of the run of measured concentrations with virial mass are often either steeper or flatter than predicted by simulations. In this work, we present the most complete sample of observed cluster concentrations and masses yet assembled, including new measurements for 10 strong lensing clusters, thereby more than doubling the existing number of strong lensing concentration estimates. We fit a power law to the observed concentrations as a function of virial mass, and find that the slope is consistent with the slopes found in simulations, though our normalization factor is higher. Observed lensing concentrations appear to be systematically larger than X-ray concentrations, a more pronounced effect than found in simulations. We also find that at fixed mass, the bulk of observed cluster concentrations are distributed log-normally, with the exception of a few anomalously high concentration clusters. We examine the physical processes likely responsible for the discrepancy between lensing and X-ray concentrations, and for the anomalously high concentrations in particular. The forthcoming Millennium simulation results will offer the most comprehensive comparison set to our findings of an observed concentration-mass power law relation.

In this paper we derive the mass function of seed black holes that result from the central mass concentrated via disc accretion in collapsed haloes at redshift $z\approx 15$. Using standard arguments including stability, we show that these pre-galactic discs can assemble a significant mass concentration in the inner regions, providing fuel for the formation and initial growth of super-massive black holes. Assuming that these mass concentrations do result in central seed black holes, we determine the mass distribution of these seeds as a function of key halo properties. The seed mass distribution determined here turns out to be asymmetric and skewed to higher masses. Starting with these initial seeds, building up to $10^9$ solar masses by $z = 6$ to power the bright quasars is not a problem in the standard LCDM cosmogony. These seed black holes in gas rich environments are likely to grow into the supermassive black holes at later times via mergers and accretion. Gas accretion onto these seeds at high redshift will produce miniquasars that likely play an important role in the reionization of the Universe. Some of these seed black holes on the other hand could be wandering in galaxy haloes as a consequence of frequent mergers, powering the off-nuclear ultra-luminous X-ray sources detected in nearby galaxies.

We present high-resolution mass reconstructions for five massive cluster-lenses spanning a redshift range from $z = 0.18$--0.57 utilising archival {\it Hubble Space Telescope} ({\it HST}) data and applying galaxy-galaxy lensing techniques. These detailed mass models were obtained from the observations by combining constraints from the strong and weak lensing regimes. We ascribe local weak distortions in the shear maps to perturbations induced by the presence of galaxy haloes around individual bright early-type cluster member galaxies. This technique constrains the mass enclosed within an aperture for these subhaloes. We are sensitive to a specific mass range for these subhaloes, $10^{11}$ -- $10^{12.5} \msun$, which we associate with galaxy-scale subhaloes. Adopting a parametric model for the subhaloes, we also derive their velocity dispersion function and the aperture radius function. The mass spectrum of substructure in the inner regions of the observed clusters is directly compared with that in simulated clusters extracted from the {\it Millennium Simulation}. The massfunction, aperture radii and velocity dispersion function are compared in detail. Overall, we find good agreement between the distribution of substructure properties retrieved using the lensing analysis and those obtained from the simulation (truncated abstract).

We present a new determination of the cluster mass function and velocity dispersion function in a volume $\sim10^7 h^3$Mpc$^{-3}$ using the Fourth Data Release of the Sloan Digital Sky Survey (SDSS). We use the caustic technique to remove foreground and background galaxies. The cluster virial mass function agrees well with recent estimates from both X-ray observations and cluster richnesses. The mass function lies between those predicted by the First-Year and Three-Year WMAP data. We constrain the cosmological parameters $\Omega_m$ and $\sigma_8$ and find good agreement with WMAP and constraints from other techniques. With the CIRS mass function alone, we estimate $\Omega_m=0.24^{+0.14}_{-0.09}$ and $\sigma_8=0.92^{+0.24}_{-0.19}$, or $\sigma_8=0.84\pm$0.03 when holding $\Omega_m=0.3$ fixed. We also use the WMAP parameters as priors and constrain velocity segregation in clusters. Using the First and Third-Year results, we infer velocity segregation of $\sigma_{gxy}/\sigma_{DM}\approx0.94\pm$0.05 or 1.28$\pm$0.06 respectively. We compare the velocity dispersion function of clusters to that of early-type galaxies and conclude that clusters comprise the high-velocity end of the velocity dispersion function of dark matter haloes. The evolution of cluster abundances provides constraints on dark energy models; the mass function presented here offers an important low redshift calibration benchmark.

Our aim is to constrain the properties of dark matter halos inhabiting high density environments, such as is the case in massive galaxy clusters. We use galaxy-galaxy lensing techniques that utilize a maximum likelihood method to constrain the parameters of the lenses. It has been demonstrated that such a technique provides strong constraints on the parameters that characterize a galaxy halo, as well as on the aperture mass of these halos. In this analysis, we only use weak shear data and do not include strong lensing constraints. We present the results of a study of galaxy-galaxy lensing in a homogeneous sample of massive X-ray luminous clusters at z~0.2. These have been observed in three bands with the CFH12k instrument. We find dark matter halos in these clusters to be compact compared to those inferred around isolated field galaxies of equivalent luminosity at this redshift: the half mass radius is found to be smaller than 50 kpc, with a mean total mass of order 0.2 10^{12} M_sun. This is in good agreement with previous galaxy-galaxy lensing results and with numerical simulations, in particular with the tidal stripping scenario. We thus provide a strong confirmation of tidal truncation from a homogeneous sample of galaxy clusters. Moreover, it is the first time that cluster galaxies are probed successfully using galaxy-galaxy lensing techniques from ground based data.

The lack of an observed supernova associated with GRB 060614 appears to require a new paradigm for the formation of (a subset of) long-duration GRBs. This requirement is based on the presumed low redshift of the burst, which was inferred from the spatial coincidence of the afterglow with a z=0.125 galaxy. We explore the possibility that this low-redshift galaxy is a chance superposition along the line of sight to GRB 060614. We examine the galaxy distribution of the field of GRB 060614 and find that the probability of a chance association with a galaxy at least as bright as the putative host is only ~ 0.5 - 1.9%. However, for the current ensemble of ~180 Swift GRBs it is likely that several such coincidences have occurred, and given the "non-standard" nature of GRB 060614 it is not implausible that this is one such occurrence. Thus the conclusion that GRB 060614 requires a revision to the formation paradigm for long-duration GRBs should be approached with caution.

We report on the first results of an imaging survey to detect strong gravitational lensing targeting the richest clusters selected from the photometric data of the Sloan Digital Sky Survey (SDSS) with follow-up deep imaging observations from the Wisconsin Indiana Yale NOAO (WIYN) 3.5m telescope and the University of Hawaii 88-inch telescope (UH88). The clusters are selected from an area of 8000 deg^2 using the Red Cluster Sequence technique and span the redshift range 0.1 < z < 0.6, corresponding to a comoving cosmological volume of ~ 2 Gpc^3. Our imaging survey thus targets a volume more than an order of magnitude larger than any previous search. A total of 240 clusters were imaged of which 141 had sub-arcsecond image quality. Our survey has uncovered16 new lensing clusters with definite giant arcs, an additional 12 systems for which the lensing interpretation is very likely, and 9 possible lenses which contain shorter arclets or candidate arcs which are less certain and will require further observations to confirm their lensing origin. The number of new cluster lenses detected in this survey is likely > 30. Among these new systems are several of the most dramatic examples of strong gravitational lensing ever discovered with multiple bright arcs at large angular separation. These will likely become 'poster-child' gravitational lenses similar to Abell 1689 and CL0024+1654. The new lenses discovered in this survey will enable future sysetmatic studies of the statistics of strong lensing and its implications for cosmology and our structure formation paradigm.

We calculate the expected lensing statistics of the galaxy population in large, low-redshift surveys. Galaxies are modeled using realistic, multiple components: a dark matter halo, a bulge component and disc. We use semi-analytic models of galaxies coupled with dark matter haloes in the Millennium Run to model the lens galaxy population. We predict that a fraction of 1.4+/-0.18*10^-3 of radio sources will be lensed by galaxies within a survey like the 2dF below z<0.2. With a simulated sample of lensed radio sources, the predicted lensing galaxy population consists mainly of ellipticals (~80%) with an average lens velocity dispersion of 164+/-3 km/s, producing typical image separations of ~3 arcsec. The lens galaxy population lies on the fundamental plane but its velocity dispersion distribution is shifted to higher values compared to all early-type galaxies. Taking magnification bias into account, we predict that the ratio of 4:2 image systems is 30+/-5%, consistent with the observed ratio found in the Cosmic Lens All-Sky Survey. We also find that the population of 4-image lens galaxies is distinguishable from the population of lens galaxies in 2-image systems. Our key result is the explicit demonstration that the population of lens galaxies differs markedly from the galaxy population as a whole: lens galaxies have a higher average luminosity and reside in more massive haloes than the overall sample of ellipticals. This bias restricts our ability to infer galaxy evolution parameters from a sample of lensing galaxies. (abridged)

Gravitational lensing has now become a popular tool to measure the mass distribution of structures in the Universe on various scales. Here we focus on the study of galaxy's scale dark matter halos with galaxy-galaxy lensing techniques: observing the shapes of distant background galaxies which have been lensed by foreground galaxies allows us to map the mass distribution of the foreground galaxies. The lensing effect is small compared to the intrinsic ellipticity distribution of galaxies, thus a statistical approach is needed to derive some constraints on an average lens population. An advantage of this method is that it provides a probe of the gravitational potential of the halos of galaxies out to very large radii, where few classical methods are viable, since dynamical and hydrodynamical tracers of the potential cannot be found at this radii. We will begin by reviewing the detections of galaxy-galaxy lensing obtained so far. Next we will present a maximum likelihood analysis of simulated data we performed to evaluate the accuracy and robustness of constraints that can be obtained on galaxy halo properties. Then we will apply this method to study the properties of galaxies which stand in massive cluster lenses at z~0.2. The main result of this work is to find dark matter halos of cluster galaxies to be significantly more compact compared to dark matter halos around field galaxies of equivalent luminosity, in agreement with early galaxy-galaxy lensing studies and with theoretical expectations, in particular with the tidal stripping scenario. We thus provide a strong confirmation of tidal truncation from a homogeneous sample of galaxy clusters. Moreover, it is the first time that cluster galaxies are probed successfully using galaxy-galaxy lensing techniques from ground based data.

In this paper we discuss the evolution of gravitationally unstable pre-galactic discs that result from the collapse of haloes at high redshift $z \approx 10$ or so, which have not yet been enriched by metals. In cases where molecular hydrogen formation is suppressed the discs are maintained at a temperature of a few thousand degrees Kelvin. However, when molecular hydrogen is present cooling can proceed down to a few hundred degrees Kelvin. Analogous to the case of the larger scale proto-galactic discs, we assume that the evolution of these discs is mainly driven by angular momentum redistribution induced by the development of gravitational instabilities in the disc. We also properly take into account the possibility of disc fragmentation. We thus show that this simple model naturally predicts the formation of supermassive black holes in the nuclei of such discs and provides a robust determination of their mass distribution as a function of halo properties. We estimate that roughly 5% of discs resulting from the collapse of haloes with $M\approx 10^7 M_{\odot}$ should host a massive black hole with a mass $M_{\rm BH}\approx 10^5 M_{\odot}$. We confirm our arguments with time-dependent calculations of the evolution of the surface density and of the accretion rate in these primordial discs. This mechanism offers an efficient way to form seed black holes at high redshift. The predicted masses for our black hole seeds enable the comfortable assembly of $10^9 M_{\odot}$ black holes powering the luminous quasars detected by the Sloan Digital Sky Survey at $z = 6$ for a concordance cosmology. (abridged)

We present SMARTS consortium optical/IR light curves of SN 2006aj associated with GRB 060218. We find that this event is broadly similar to two previously observed events SN 1998bw/GRB 980425 and SN 2003lw/GRB 031203. In particular, all of these events are greatly under-luminous in gamma-rays compared to typical long-duration GRBs. We find that the observation by Swift of even one such event implies a large enough true event rate to create difficulties in interpreting these events as typical GRBs observed off-axis. Thus these events appear to be intrinsically different from and much more common than high-luminosity GRBs, which have been observed in large numbers out to a redshift of at least 6.3. The existence of a range of intrinsic energies of GRBs may present challenges to using GRBs as standard candles.

Gravitational lensing of a background source by a foreground galaxy lens occasionally produces four images of the source. The cusp and the fold relations impose conditions on the ratios of magnifications of these four-image lenses. In this theoretical investigation, we explore the sensitivity of these relations to the presence of substructure in the lens. Starting with a smooth lens potential, we add varying amounts of substructure, while keeping the source position fixed, and find that the fold relation is a more robust indicator of substructure than the cusp relation for the images. This robustness is independent of the detailed spatial distribution of the substructure, as well as of the ellipticity of the lensing potential and the presence of external shear.

In this paper we examine how well galaxies and intra-cluster gas trace the gravitational potential of clusters. Utilizing mass profiles derived from gravitational lensing and X-ray observations, coupled with measured galaxy velocities, we solve for the velocity anisotropy parameter using the anisotropic Jeans equation. This is done for five clusters, three at low redshift: A2199, A496 and A576 and two at high redshifts: A2390 and MS1358. With X-ray temperature profiles obtained from Chandra and ASCA/ROSAT data, we estimate the ratio of energy in the galaxies compared to the X-ray gas. We find that none of these clusters is strictly in hydro-static equilibrium. We compare the properties of our sample with clusters that form in high-resolution cosmological N-body simulations that include baryonic physics. Simulations and data show considerable scatter both these profiles. We demonstrate the future feasibility and potential for directly comparing the orbital structure of clusters inferred from multi-wavelength observations with high resolution simulated clusters.

We report optical, near-infrared and X-ray observations of the afterglow of GRB 050814, which was seen to exhibit very red optical colours. By modelling its spectral energy distribution we find that z = 5.3 +/- 0.3. We next present a carefully selected sample of 19 Swift GRBs, intended to estimate in an unbiased way the GRB redshift distribution, including the mean redshift (z_mean) as well as constraints on the fraction of high-redshift bursts. We find that z_mean = 2.7 and that at least 5% of the GRBs originate at z > 5. The redshift distribution of the sample is qualitatively consistent with models where the GRB rate is proportional to the star formation rate in the Universe. The high mean redshift of this GRB sample and the wide redshift range clearly demonstrates the suitability of GRBs as efficient probes of galaxies and the intergalactic medium over a significant fraction of the history of the Universe.

We present mid-infrared observations of AGN in the GOODS fields, performed with the Spitzer Space Telescope. These are the deepest infrared and X-ray fields to date and cover a total area of ~0.1 square degrees. AGN are selected on the basis of their hard (2-8 keV) X-ray emission. The median AGN infrared luminosity is at least 10 times larger than the median for normal galaxies with the same redshift distribution, suggesting that the infrared emission is dominated by the central nucleus. The X-ray to infrared luminosity ratios of GOODS AGN, most of which are at 0.5<z<1.5, are similar to the values obtained for AGN in the local Universe. The observed infrared flux distribution has an integral slope of ~1.5 and there are 1000 sources per square degree brighter than ~50 uJy at 3-6 microns. The counts approximately match the predictions of models based on AGN unification, in which the majority of AGN are obscured. This agreement confirms that the faintest X-ray sources, which are dominated by the host galaxy light in the optical, are obscured AGN. Using these Spitzer data, the AGN contribution to the extragalactic infrared background light is calculated by correlating the X-ray and infrared catalogues. This is likely to be a lower limit given that the most obscured AGN are missed in X-rays. We estimate the contribution of AGN missed in X-rays, using a population synthesis model, to be ~45% of the observed AGN contribution, making the AGN contribution to the infrared background at most ~2-10% in the 3-24 micron range, depending on wavelength, lower than most previous estimates. The AGN contribution to the infrared background remains roughly constant with source flux in the IRAC bands but decreases with decreasing flux in the MIPS 24 um band, where the galaxy population becomes more important.

The exceptionally high luminosities of gamma-ray bursts (GRBs), gradually emerging as extremely useful probes of star formation, make them promising tools for exploration of the high-redshift Universe. Here we present a carefully selected sample of Swift GRBs, intended to estimate in an unbiased way the GRB mean redshift (z_mean), constraints on the fraction of high-redshift bursts and an upper limit on the fraction of heavily obscured afterglows. We find that z_mean = 2.8 and that at least 7% of GRBs originate at z > 5. In addition, consistent with pre-Swift observations, at most 20% of afterglows can be heavily obscured. The redshift distribution of the sample is qualitatively consistent with models where the GRB rate is proportional to the star formation rate in the Universe. We also report optical, near-infrared and X-ray observations of the afterglow of GRB 050814, which was seen to exhibit very red optical colours. By modelling its spectral energy distribution we find that z = 5.3 +/- 0.3. The high mean redshift of GRBs and their wide redshift range clearly demonstrates their suitability as efficient probes of galaxies and the intergalactic medium over a significant fraction of the history of the Universe.

In this letter, we calculate the redshift distribution of gamma-ray bursts assuming that they trace (i) the globally averaged star formation rate or (ii) the average metallicity in the Universe. While at redshifts 5 and below, both the star formation rate and metallicity are observationally determined modulo some uncertainties, at higher redshifts there are few constraints. We extrapolate the star formation rate and metallicity to higher redshifts and explore models that are broadly consistent with bounds on the optical depth from WMAP results. In addition, we also include parametric descriptions of the luminosity function, and the typical spectrum for GRBs. With these essential ingredients included in the modeling, we find that a substantial fraction 75% of GRBs are expected to originate at redshifts below 4, in variance with some previous estimates. Conversely, if we assume as expected for the collapsar model that gamma-ray bursts favour a low metallicity environment and therefore, relate the GRB rate to a simple model of the average metallicity as a function of redshift, we find that a higher fraction of bursts, about 40% originate from redshifts greater than 4. We conclude with the implications of SWIFT GRB detections.

Angular momentum loss to circumbinary gas provides a possible mechanism for overcoming the `last parsec' problem and allowing the most massive black hole binaries formed from galactic mergers to coalesce. Here, we show that if gas disks also catalyze the merger of the somewhat lower mass binaries detectable with the Laser Interferometer Space Antenna (LISA), then there may be a purely gravitational wave signature of the role of gas in the form of a small but finite eccentricity just prior to merger. Numerical simulations suggest that eccentricity, excited by the interaction between the binary and surrounding gas disk, is only partially damped during the final phase of gravitational radiation-driven inspiral. We estimate that the typical eccentricity at one week prior to coalescence is of the order of 0.01. Higher terminal eccentricities, which can approach 0.1, are possible if the binary has an extreme mass ratio. The detection of even a small eccentricity prior to merger by LISA provides a possible discriminant between gas-driven inspirals and those effected by stellar processes.

It is widely recognized that cold dark matter models predict abundant dark matter substructure in halos of all sizes. Galaxy-galaxy lensing provides a unique opportunity to directly measure the presence and the mass of such substructures in clusters of galaxies. Here we present the mass function of substructures obtained from lensing in massive HST cluster lenses, and compare it with that obtained from high-resolution cosmological N-body simulations. We find excellent agreement in the slope and amplitude in the mass range 10^{11} - 10^{12.5} solar masses probed by the observations, highlighting a significant success of the CDM theory. At lower mass scales, the simulations predict a large abundance of small substructures below our detection threshold with galaxy-galaxy lensing, with many of them being plausibly associated with faint cluster galaxies. Our results suggest that the CDM substructure abundance on the scale of clusters is in good agreement with present observational data.

We present high-resolution mass reconstructions for five massive cluster-lenses spanning a redshift range from z = 0.18 - 0.57 utilizing archival Hubble Space Telescope data and applying galaxy-galaxy lensing techniques. These detailed mass models were obtained by combining constraints from the observed strong and weak lensing regimes. Quantifying the local weak distortions in the shear maps in terms of perturbations induced by the presence of galaxy halos around individual bright early-type cluster member galaxies, we estimate the fraction of mass in the central regions of these clusters that can be associated with small scale mass clumps. This technique enables us to directly map the substructure in the mass range 10^{11} - 10^{12.5} solar masses which we associate with galaxy-scale sub-halos. The determination of the mass spectrum of substructure in the inner regions of these clusters is presented. Constraints are thereby obtained on the masses, mass-to-light ratios and truncation radii for these sub-halos. We find that the fraction of total cluster mass associated with individual sub-halos within the inner 0.5h^{-1} - 0.8h^{-1} Mpc of these clusters ranges from 10 - 20%. Our results have important implications for the survival and evolution of substructure in high density cluster cores and are consistent with the theoretical picture of tidal stripping of galaxy-scale halos in high-density cluster environments as expected in hierarchical Cold Dark Matter dominated structure formation scenarios.

Long duration gamma-ray bursts (GRBs) accompany the deaths of some massive stars and hence, since massive stars are short lived, are a tracer of star formation activity. Given that GRBs are bright enough to be seen to very high redshifts, and detected even in dusty environments, they should therefore provide a powerful probe of the global star formation history of the universe. The potential of this approach can be investigated via submm photometry of GRB host galaxies. Submm luminosity also correlates with star formation rate, so the distribution of host galaxy submm fluxes should allow us to test the two methods for consistency. Here, we report new JCMT/SCUBA 850 micron measurements for 15 GRB hosts. Combining these data with results from previous studies we construct a sample of 21 hosts with <1.4 mJy errors. We show that the distribution of apparent 850 micron flux densities of this sample is reasonably consistent with model predictions, but there is tentative evidence of a dearth of submm bright (>4 mJy) galaxies. Furthermore, the optical/infrared properties of the submm brightest GRB hosts are not typical of the galaxy population selected in submm surveys, although the sample size is still small. Possible selection effects and physical mechanisms which may explain these discrepancies are discussed.

We present a maximum-likelihood analysis of galaxy-galaxy lensing effects in galaxy clusters and in the field. The aim is to determine the accuracy and robustness of constraints that can be obtained on galaxy halo properties in both environments - the high density cluster and the low density field. This paper is theoretically motivated, therefore, we work exclusively with simulated data (nevertheless defined to match observations) to study the accuracy with which input parameters for mass distributions for galaxies can be extracted. We model galaxies in the cluster and the field using a wide range of mass profiles: the truncated pseudo isothermal elliptical mass distribution, the Navarro, Frenk and White profile, and a Power Law model with a core radius. We find that independent of the choice of profile the mean mass of galaxies (of the order of 10^{12}Mo) can be estimated to within 15% from ground-based data and with an error of less than 10% with space observations. Additionally robust constraints can be obtained on the mean slope of the mass profile. The two standard parameters that characterise galaxy halo models, the central velocity dispersion and the truncation radius can also be retrieved reliably from the maximum-likelihood analysis. Furthermore, going beyond the usual formulation, we propose a re-parameterisation of the mass models that allows us to put yet stronger constraints on the aperture mass of a galaxy halo (with less than 10% error). The gain in signal to noise using space observations, expected for instance with the proposed SNAP satellite compared to ground based data in terms of accuracy of retrieving input parameters is highly significant.

There is overwhelming evidence for the presence of supermassive black holes (SMBHs) in the centers of most nearby galaxies. The mass estimates for these remnant black holes from the stellar kinematics of local galaxies and the quasar phenomenon at high redshifts point to the presence of assembled SMBHs. The accretion history of SMBHs can be reconstructed using observations at high and low redshifts as model constraints. Observations of galaxies and quasars in the submillimeter, infrared, optical, and X-ray wavebands are used as constraints, along with data from the demography of local black holes. Theoretical modeling of the growth of black hole mass with cosmic time has been pursued thus far in two distinct directions: a phenomenological approach that utilizes observations in various wavebands, and a semi-analytic approach that starts with a theoretical framework and a set of assumptions with a view to matching observations. Both techniques have been pursued in the context of the standard paradigm for structure formation in a Cold Dark Matter dominated universe. Here, we examine the key issues and uncertainties in the theoretical understanding of the growth of SMBHs.

(abridged) We present a comprehensive lensing analysis of the rich cluster Cl0024+1654 based on panoramic sparse-sampled imaging conducted with HST. We demonstrate an ability to detect reliably weak lensing signals to a cluster radius of ~5 Mpc where the mean shear is around 1%. The projected mass distribution reveals a secondary concentration representing 30% of the overall cluster mass, which is also visible in the distribution of cluster member galaxies. We normalize the mass profile determined from the shear by assuming that background galaxies selected with I=23-26 have a redshift distribution statistically similar to that in the HDFs. Combining strong and weak constraints, we are able to probe the mass profile of the cluster on scales of 0.1 to 5 Mpc thus providing a valuable test of the universal form proposed by NFW on large scales. A generalized power law fit indicates a asymptotic density distribution with 3D slope larger than 2.4. An isothermal mass profile is therefore strongly rejected, whereas a NFW profile with M_200= 6.1+/-1 10^14 M_sun provides a good fit to the lensing data. We isolate cluster members according to their optical-near infrared colors; the red cluster light closely traces the dark matter with a mean mass-to-light ratio of M/L_K~ 40 M/L_sun. Similar profiles for mass and light on 1-5 Mpc scales are expected if cluster assembly is largely governed by infalling groups.

We describe a new wide field Hubble Space Telescope survey of the galaxy cluster Cl0024+16 (z~0.4) consisting of a sparse-sampled mosaic of 39 Wide Field and Planetary Camera 2 images which extends to a cluster radius of 5 Mpc. [Abridged] We examine both the morphology-radius (T-R) and morphology-density (T-Sigma) relations and demonstrate sensitivities adequate for measures from the core to a radius of 5 Mpc, spanning over 3 decades in local projected density. The fraction of early-type galaxies declines steeply from the cluster center to 1 Mpc radius and more gradually thereafter, asymptoting towards the field value at the periphery. We discuss our results in the context of three distinct cluster zones, defined according to different physical processes that may be effective in transforming galaxy morphology in each. By treating infalling galaxies as isolated test particles, we deduce that the most likely processes responsible for the mild gradient in the morphological mix outside the virial radius are harassment and starvation. Although more data are needed to pin down the exact mechanisms, starvation seems more promising in that it would naturally explain the stellar and dynamical homogeneity of cluster E/S0s. However, we find significant scatter in the local density at any given radius outside 0.5 Mpc, and that the same T-Sigma relation holds in subregions of the cluster, independent of location. In this hitherto unprobed region, where the potential of the cluster is weak, galaxies apparently retain their identities as members of infalling sub-groups whose characteristic morphological properties remain intact. Only upon arrival in the central regions is the substructure erased, as indicated by the tight correlation between cluster radius and Sigma.

We study the observational signatures of the lensing signal produced by dark matter halos with embedded misaligned disks. This issue is of particular interest at the present time since most of the observed multiple lens systems have magnification ratios and image geometries that are not well-fit by standard mass models. The presence of substructure exterior to the lens has been invoked by several authors in the context of Cold Dark Matter (CDM) models in order to explain the anamolous magnification ratios. We emphasize that the anomalous magnification ratios may be an artifact of the simple one-component mass models currently in use; the inclusion of a misaligned disk may be able to mimic the effect of substructure. These slight spatial offsets between the dark matter halo and the disk, which are likely to occur during or as a consequence of interactions or mergers, lead to complex image configurations and non-standard magnification ratios. We investigate two illustrative lens models: a spiral disk embedded within a dark matter halo, and a compact disk-like component within an elliptical galaxy. The expected fraction of galaxies with a misaligned disk is estimated to be of the order of 10%. In such cases we find that the resultant lensing geometries are unusual, with high image multiplicities. The caustic structures - both radial and tangential - are drastically modified and the magnification ratios differ compared to expectations from standard lens models. The additional parameters required to specify the relative alignment of the mass components in the primary lens introduce yet another source of uncertainty in the mass modeling of gravitational lens systems.

In this letter, we present the results of our study of galaxy-galaxy lensing in massive cluster-lenses spanning $z = 0.17$ to 0.58, utilizing high-quality archival {\it Hubble Space Telescope} ({\it HST}) data. Local anisotropies in the shear maps are assumed to arise from dark matter substructure within these clusters. Associating the substructure with bright early-type cluster galaxies, we quantify the properties of typical $L^*$ cluster members in a statistical fashion. The fraction of total mass associated with individual galaxies within the inner regions of these clusters ranges from 10--20% implying that the bulk of the dark matter in massive lensing clusters is smoothly distributed. Looking at the properties of the cluster galaxies, we find strong evidence ($>3$-$\sigma$ significance) that a fiducial early-type $L^\ast$ galaxy in these clusters has a mass distribution that is tidally truncated compared to equivalent luminosity galaxies in the field. In fact, we exclude field galaxy scale dark halos for these cluster early-types at $>10$-$\sigma$ significance. We compare the tidal radii obtained from this lensing analysis with the central density of the cluster potentials and find a correlation which is in excellent agreement with theoretical expectations of tidal truncation: $\log [r_t*] \propto (-0.6\pm 0.2) \log [\rho_0]$.

The detailed distribution of mass within clusters of galaxies can be used to probe the nature of dark matter. We show that constraints on the extent of the mass distribution around galaxies in the rich cluster Abell 2218 obtained from combining strong and weak lensing observations are consistent with the predictions which assume that the dominant mass component (dark matter) in these halos is collisionless. A strongly interacting (fluid-like) dark matter is ruled-out at a confidence level of more than 5-sigma.

We present the results of a search for submillimetre-luminous host galaxies of optically dark gamma-ray bursts (GRBs) using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). We made photometry measurements of the 850-micron flux at the location of four `dark bursts', which are those with no detected optical afterglow despite rapid deep searches, and which may therefore be within galaxies containing substantial amounts of dust. We were unable to detect any individual source significantly. Our results are consistent with predictions for the host galaxy population as a whole, rather than for a subset of dusty hosts. This indicates that optically dark GRBs are not especially associated with very submillimetre-luminous galaxies and so cannot be used as reliable indicators of dust-enshrouded massive star-formation activity. Further observations are required to establish the relationship between the wider GRB host galaxy population and SCUBA galaxies.

We report the identification of the optical afterglow of GRB 980613 in R- and I-band images obtained between 16 and 48 hours after the gamma-ray burst. Early near-infrared (NIR) H and K' observations are also reported. The afterglow was optically faint (R ~ 23) at discovery but did not exhibit an unusually rapid decay (power-law decay slope alpha < 1.8 at 2 sigma). The optical/NIR spectral index (beta_RH < 1.1) was consistent with the optical-to-X-ray spectral index (beta_RX ~ 0.6), indicating a maximal reddening of the afterglow of ~0.45 mag in R. Hence the dimness of the optical afterglow was mainly due to the fairly flat spectral shape rather than internal reddening in the host galaxy. We also present late-time HST/STIS images of the field in which GRB 980613 occurred, obtained 799 days after the burst. These images show that GRB 980613 was located close to a very compact, blue V = 26.1 object inside a complex region consisting of star-forming knots and/or interacting galaxy fragments. Therefore, GRB 980613 constitutes a strong case for the association of cosmological gamma-ray bursts with star-forming regions.

Recent studies suggest that the initial mass function (IMF) of the first stars was likely to be extremely top-heavy, unlike what is observed at present. We propose a scenario to generate fragmentation to lower masses once the first massive stars have formed and derive constraints on the primordial IMF. We estimate the mass fraction of pair-unstable supernovae, shown to be the dominant sources of the first heavy elements. These metals enrich the gas up to about $10^{-5}$ solar metallicity, when a transition to efficient cooling-driven fragmentation occurs producing 1 solar mass clumps. We argue that the remaining fraction of the first stars ends up in 100 solar mass VMBHs (Very Massive Black Holes). We obtain constraints on the fraction of first stars that contribute to the initial metal enrichment and the transition redshift for primordial IMF away from a top-heavy one, by making various assumptions about the fate of these VMBHs at late times. We conclude with a discussion of several hitherto unexplored implications of a high-mass dominated star formation mode in the early Universe.

As part of the HST/STIS GRB host survey program we present the detection of three faint gamma-ray burst (GRB) host galaxies based on an accurate localisation using ground-based data of the optical afterglows (OAs). A common property of these three hosts is their extreme faintness. The location at which GRBs occur with respect to their host galaxies and surrounding environments are robust indicators of the nature of GRB progenitors. The bursts studied here are among the four most extreme outliers, in terms of relative distance from the host center, in the recent comprehensive study of Bloom et al. 2002. We obtain a revised and much higher probability that the galaxies identified as hosts indeed are related to the GRBs (P(n_chance)=0.78, following Bloom et al. 2002), thereby strengthening the conclusion that GRBs are preferentially located in star-forming regions in their hosts. Apart from being faint, the three hosts consist of multiple structures, indicative of merging and active star-formation. For one of the hosts, GRB 980329, we estimate a photometric redshift of z~3.5.

We study the evolution of disk accretion during the merger of supermassive black hole binaries in galactic nuclei. In hierarchical galaxy formation models, the most common binaries are likely to arise from minor galactic mergers, and have unequal mass black holes. Once such a binary becomes embedded in an accretion disk at a separation of the order of 0.1 pc, the merger proceeds in two distinct phases. During the first phase, the loss of orbital angular momentum to the gaseous disk shrinks the binary on a timescale of ~10 Myr. The accretion rate onto the primary black hole is not increased, and can be substantially reduced, during this disk-driven migration. At smaller separations, gravitational radiation becomes the dominant angular momentum loss process, and any gas trapped inside the orbit of the secondary is driven inwards by the inspiralling black hole. The implied accretion rate just prior to coalescence exceeds the Eddington limit, so the final merger is likely to occur within a common envelope formed from the disrupted inner disk, and be accompanied by high velocity (~10,000 km/s) outflows.

We describe the first results from a comprehensive study of the distant cluster Cl0024+1654 (z=0.39) based upon a pattern of 38 mosaiced HST-WFPC2 images extending to radii ~5Mpc. These are being analysed in conjunction with extensive spectroscopy conducted with the CFHT, WHT, and Keck Telescopes. The overall goal is to understand the morphological transformations and associated short-term star formation histories of representative numbers of infalling field galaxies in the context of the cluster potential as defined by weak lensing studies. Our HST database contains over 2000 galaxy morphologies to I=22.5. Spectroscopic data and HST morphologies are currently available for about 215 members over an unprecedented range of environments. We confirm the existence of a well-defined morphology-density relation over a large dynamic range within a single system at a significant look-back time. Tentative trends in the E/S0 fraction as a function of radius are discussed. A weak lensing signal in the background galaxies has been detected at the cluster periphery and its inversion demonstrates only marginal substructure. A statistically-significant galaxy-galaxy lensing signal has also been seen for the cluster members. Further work will relate radial dependencies in the dark matter and halo masses in the context of spectroscopic and morphological diagnostics of truncated star formation.

In this paper, we study gravitational lensing by groups of galaxies. Since groups are abundant and therefore have a large covering fraction on the sky, lensing by groups is likely to be very important observationally. Besides, it has recently become clear that many lens models for strong lensing by individual galaxies require external shear to reproduce the observed image geometries; in many cases a nearby group is detected that could provide this shear. In this work, we study the expected lensing behavior of galaxy groups in both the weak and strong lensing regime. We examine the shear and magnification produced by a group and its dependence on the detailed mass distribution within the group. We find that the peak value of the weak lensing shear signal is of the order of 3 per cent and varies by a factor of about 2 for different mass distributions. These variations are large enough to be detectable in the Sloan Digital Sky Survey (SDSS). In the strong lensing regime we find that the image geometries and typical magnifications are sensitive to the group properties and that groups can easily provide enough external shear to produce quadruple images. We investigate the statistical properties of lensing galaxies that are near or part of a group and find that statistical lens properties, like the distribution of time delays, are affected measurably by the presence of the group which can therefore introduce an additional systematic error in the measurement of the Hubble constant from such systems. We conclude that both the detection of weak lensing by groups and accurate observations of strong galaxy lens systems near groups could provide important information on the total mass and matter distribution within galaxy groups.

In this paper, we introduce the octopole moment of the light distribution in galaxies as a probe of the weak lensing shear field. While traditional ellipticity estimates of the local shear derived from the quadrupole moment are limited by the width of the intrinsic ellipticity distribution of background galaxies, the dispersion in the intrinsic octopole distribution is expected to be much smaller, implying that the signal from this higher order moment is ultimately limited by measurement noise, and not by intrinsic scatter. We present the computation of the octopole moment and show that current observations are at the regime where the octopole estimates will soon be able to contribute to the overall accuracy of the estimates of local shear fields. Therefore, the prospects for this estimator from future datasets like the Advanced Camera for Survey and the Next Generation Space Telescope are very promising.

We present deep HST/STIS and ground-based photometry of the host galaxy of the gamma-ray burst GRB 980703 taken 17, 551, 710, and 716 days after the burst. We find that the host is a blue, slightly over-luminous galaxy with V_gal = 23.00 +/- 0.10, (V-R)_gal = 0.43 +/- 0.13, and a centre that is approximately 0.2 mag bluer than the outer regions of the galaxy. The galaxy has a star-formation rate of 8-13 M_sun/yr, assuming no extinction in the host. We find that the galaxy is best fit by a Sersic R^(1/n) profile with n ~= 1.0 and a half-light radius of 0.13 arcsec (= 0.72/h_100 proper kpc). This corresponds to an exponential disk with a scale radius of 0.22 arcsec (= 1.21/h_100 proper kpc). Subtracting a fit with elliptical isophotes leaves large residuals, which suggests that the host galaxy has a somewhat irregular morphology, but we are unable to connect the location of GRB 980703 with any special features in the host. The host galaxy appears to be a typical example of a compact star forming galaxy similar to those found in the Hubble Deep Field North. The R-band light curve of the optical afterglow associated with this gamma-ray burst is consistent with a single power-law decay having a slope of alpha = -1.37 +/- 0.14. Due to the bright underlying host galaxy the late time properties of the light-curve are very poorly constrained. The decay of the optical light curve is consistent with a contribution from an underlying Type Ic supernova like SN1998bw, or a dust echo, but such contributions cannot be securely established.

Alignments in the angular momentum vectors of galaxies can induce large scale correlations in their projected orientations. Such alignments arise from the tidal torques exerted on neighboring proto-galaxies by the smoothly varying shear field. Weak gravitational lensing can also induce ellipticity correlations since the images of neighboring galaxies will be distorted coherently by the intervening mass distribution. Comparing these two sources of shape correlations, it is found that for current weak lensing surveys with a median redshift of 1, the intrinsic signal is a contaminant on the order of 1-10% of the measured signal. However, for shallower surveys, the intrinsic correlations dominate over the lensing signal. The distortions induced by lensing are curl-free, whereas those resulting from intrinsic alignments are not. This difference can be used to disentangle these two sources of ellipticity correlations. When the distortions are dominated by lensing, as occurs at high redshifts, the decomposition provides a valuable tool for understanding properties of the noise and systematic errors.

The distortion field defined by the ellipticities of galaxy shapes projected on the sky can be uniquely decomposed into a gradient and a curl component. If the observed ellipticities are induced by weak gravitational lensing, then the distortion field is curl free. Here we show that, in contrast, the distortion field resulting from intrinsic spin alignments is not curl free. This provides a powerful discriminant between lensing and intrinsic contributions to observed ellipticity correlations. We also show how these contributions can be disentangled statistically from the ellipticity correlations or computed locally from circular integrals of the ellipticity field. This allows for an unambiguous detection of intrinsic galaxy alignments in the data. When the distortions are dominated by lensing, as occurs at high redshifts, the decomposition provides a valuable tool for understanding properties of the noise and systematic errors. These techniques can be applied equally well to the polarization of the microwave background, where it can be used to separate curl-free scalar perturbations from those produced by gravity waves or defects.

Large scale correlations in the orientations of galaxies can result from alignments in their angular momentum vectors. These alignments arise from the tidal torques exerted on neighboring proto-galaxies by the smoothly varying shear field. We compute the predicted amplitude of such ellipticity correlations using the Zel'dovich approximation for a realistic distribution of galaxy shapes. Weak gravitational lensing can also induce ellipticity correlations since the images of neighboring galaxies will be distorted coherently. On comparing these two effects that induce shape correlations, we find that for current weak lensing surveys with a median redshift of z_m = 1, the intrinsic signal is of order 1 - 10 percent of the measured signal. However, for shallower surveys with z_m < 0.3, the intrinsic correlations dominate over the lensing signal. The distortions induced by lensing are curl-free, whereas those resulting from intrinsic alignments are not. This difference can be used to disentangle these two sources of ellipticity correlations.

We present HST/STIS observations of ESO 184-G82, the host galaxy of the gamma-ray burst GRB 980425 associated with the peculiar Type Ic supernova SN1998bw. ESO 184-G82 is found to be an actively star forming SBc sub-luminous galaxy. We detect an object consistent with being a point source within the astrometric uncertainty of 0.018 arcseconds of the position of the supernova. The object is located inside a star-forming region and is at least one magnitude brighter than expected for the supernova based on a simple radioactive decay model. This implies either a significant flattening of the light curve or a contribution from an underlying star cluster.

We use simple energetic arguments to estimate the contribution of massive X-ray binaries and supernova remnants to the cosmic X-ray background (XRB) at energies in excess of 2 keV. Recent surveys have shown that AGN probably account for most of the hard XRB (E > 2 keV), but there have been many suggestions that star-forming galaxies could emerge at fainter fluxes and perhaps account for a significant fraction of the soft and hard X-ray energy density. Assuming that the formation rate of massive X-ray binaries (MXRBs) traces the global star-formation rate, we find that their integrated contribution to the hard XRB can be estimated and is shown to be small (at less than the 1% level). Similarly, the integrated flux of SN is also shown to be insignificant, or at most comparable to MXRBs. AGN therefore remain the most viable candidates for producing the hard XRB, unless additional processes can be shown to dominate the global hard X-ray emission in distant starburst galaxies.

We propose a new technique to directly measure the shapes of dark matter halos of galaxies using weak gravitational lensing. Extending the standard galaxy-galaxy lensing method, we show that the shape parameters of the mass distribution of foreground galaxies can be measured from the two-dimensional shear field derived from background galaxies. This enables the comparison of the ellipticity of the mass distribution with that of the light in galaxies, as well as an estimate of the degree of alignment between the stellar component and dark matter. We choose the specific case of an elliptical, isothermal profile and estimate the feasibility and significance of the detection of this signal. The prospects for applying this technique are excellent with large on-going surveys like the Sloan Digital Sky Survey. The expected signal is smaller, but comparable in significance to that of the mass in standard galaxy-galaxy lensing analyses. Since shapes of halos depend on the degree of dissipation and the transfer of angular momentum during galaxy assembly, constraints obtained from our analysis will provide an important input to models of galaxy formation.

Popular models for the origin of gamma-ray bursts (GRBs) include short-lived massive stars as the progenitors of the fireballs. Hence the redshift distribution of GRBs should track the cosmic star formation rate of massive stars accurately. A significant proportion of high-mass star formation activity appears to occur in regions that are obscured from view in the optical waveband by interstellar dust. The amount of dust-enshrouded star formation activity taking place has been estimated by observing the thermal radiation from the dust that has been heated by young stars in the far-infrared and submillimetre wavebands. Here we discuss an alternative probe - the redshift distribution of GRBs. GRBs are detectable at the highest redshifts, and because gamma rays are not absorbed by dust, the redshift distribution of GRBs should therefore be unaffected by dust extinction. At present the redshifts of GRBs can only be determined from the associated optical transient emission; however, useful information about the prevalence of dust-obscured star formation can also be obtained from the ratio of GRBs with and without an associated optical transient. Eight GRBs currently have spectroscopic redshifts. Once about a hundred redshifts are known, the population of GRBs will provide an important test of different models of the star formation history of the Universe.

A synopsis of the recent results from our study of galaxy-galaxy lensing in clusters is presented. We have applied our analysis techniques to the sample of HST cluster-lenses that span a redshift range from redshift 0.18 to 0.58. We find that there is evidence that the total mass of a typical early-type cluster galaxy increases with redshift. For the lowest redshift bin, a sensible comparison can be made with field galaxies and it is found that cluster galaxies are less massive and less extended than equivalent luminosity field galaxies. This agrees broadly with the theoretical picture of tidal stripping of dark matter halos in high density cluster environments. We also find the following trends -- at a fixed luminosity both the mass-to-light ratio and the fiducial truncation radius increase with redshift from redshift 0.18 to 0.58.

A large fraction of known galaxy-lens systems require a component of external shear to explain the observed image geometries. In most cases, this shear can be attributed to a nearby group of galaxies. We discuss how the dark-matter mass distribution of groups of galaxies can influence the external shear for strongly lensed sources and Calculate the expected weak lensing signal from groups for various mass profiles.

In the presence of a magnetic field, rotational energy can be extracted from black holes via the Blandford-Znajek mechanism. We use self-similar advection dominated accretion (ADAF) models to estimate the efficiency of this mechanism for black holes accreting from geometrically thick disks, in the light of recent magnetohydrodynamic disk simulations, and show that the power from electromagnetic energy extraction exceeds the accretion luminosity for ADAFs at sufficiently low accretion rates. We consider the detectability of isolated stellar mass black holes accreting from the ISM, and show that for any rapidly rotating holes the efficiency of energy extraction could reach 0.01. The estimated total luminosity would be consistent with the tentative identification of some EGRET sources as accreting isolated black holes, if that energy is radiated primarily as gamma rays. We discuss the importance of emission from the Blandford-Znajek mechanism for the spectra of other advection dominated accretion flows, especially those in low luminosity galactic nuclei.

Misaligned accretion disks surrounding rotating compact objects experience a torque due to the Lense-Thirring effect, which leads to precession of the inner disk. It has been suggested that this effect could be responsible for some low frequency Quasi-Periodic Oscillations observed in the X-ray lightcurves of neutron star and galactic black hole systems. We investigate this possibility via time-dependent calculations of the response of the inner disk to impulsive perturbations for both Newtonian point mass and Paczynski-Wiita potentials, and compare the results to the predictions of the linearized twisted accretion disk equations. For most of a wide range of disk models that we have considered, the combination of differential precession and viscosity causes the warps to decay extremely rapidly. Moreover, at least for relatively slowly rotating objects, linear calculations in a Newtonian point mass potential provide a good measure of the damping rate, provided only that the timescale for precession is much shorter than the viscous time in the inner disk. The typically rapid decay rates suggest that coherent precession of a fluid disk would not be observable, though it remains possible that the damping rate of warp in the disk could be low enough to permit weakly coherent signals from Lense-Thirring precession.

Warped accretion discs in Active Galactic Nuclei (AGN) exert a torque on the black hole that tends to align the rotation axis with the angular momentum of the outer disc. We compute the magnitude of this torque by solving numerically for the steady state shape of the warped disc, and verify that the analytic solution of Scheuer and Feiler (1996) provides an excellent approximation. We generalise these results for discs with strong warps and arbitrary surface density profiles, and calculate the timescale on which the black hole becomes aligned with the angular momentum in the outer disc. For massive holes and accretion rates of the order of the Eddington limit the alignment timescale is always short (less than a Myr), so that jets accelerated from the inner disc region provide a prompt tracer of the angular momentum of gas at large radii in the disc. Longer timescales are predicted for low luminosity systems, depending on the degree of anisotropy in the disc's hydrodynamic response to shear and warp, and for the final decay of modest warps at large radii in the disc that are potentially observable via VLBI. We discuss the implications of this for the inferred accretion history of those Active Galactic Nuclei whose jet directions appear to be stable over long timescales. The large energy deposition rate at modest disc radii during rapid realignment episodes should make such objects transiently bright at optical and infrared wavelengths.

In this letter we examine the role of the first supernovae in proto-galaxies, their role in feedback and the consequences for disk assembly. Extending the picture proposed by Dekel & Silk (1986), we argue that energetic supernovae winds can expel baryons from all proto-galaxies with varying degrees of efficiency. The fraction of baryons retained and hence available to assemble into the baryonic disk is therefore, a function of the central velocity dispersion of the halo. Such a coupling of the baryonic component to the dark halo leads to the following interesting consequence, a prediction for a weak scaling of the zero-point of the Tully-Fisher relation or alternatively, the mass-to-light ratio with the central velocity dispersion of the halo. On applying to the case of the Milky Way halo, this feedback mechanism implies: (i) that the Milky Way halo lost approximately 10% of its original gas content; (ii) a range in the inferred redshift of formation $z_f$, and the local baryon fraction $f_b$ for the Milky Way that depends on the initial spin parameter of the halo. We find that for a low spin halo - $z_f < 1$, $f_b \sim 2%$; for a median spin halo - $z_f \sim 1 - 2.5$, $f_b \sim 5%$; and for a high spin halo - $z_f \sim 4 - 8$, $f_b \sim 20%$. The observationally determined ages for the oldest disk stars in the Milky Way seem to rule out a low value for the spin parameter. Given the shape of the spin distribution of halos obtained in N-body simulations, while a high value of the spin parameter is not very probable, it is interesting to note that if this is indeed the case for the Milky Way halo, then feedback processes can cause the local baryon fraction to differ significantly from the universal value.

We show that the luminosity function of the actively star-forming Lyman break galaxies and the B-band quasar luminosity function at $z = 3$ can be fit reasonably well with the mass function of collapsed galaxy scale dark matter haloes predicted by viable variants of hierarchical cold dark matter dominated cosmological models for lifetimes $t_Q$ of the optically bright phase of QSOs in the range $10^{6}$ to $10^{8}$ yr. There is a strong correlation between $t_Q$ and the required degree of non-linearity in the relation between black hole and host halo mass. Such a non-linear relation is motivated by suggesting that the mass of supermassive black holes may be limited by the back-reaction of the emitted energy on the accretion flow in a self-gravitating disc. This would imply a relation of black hole to halo mass of the form $M_{\rm bh} \propto v_{\rm halo}^5 \propto M_{\rm halo}^{5/3}$ and a typical duration of the optically bright QSO phase of the order of the Salpeter time, $\sim 10^{7}$ yr. The high integrated local mass density of black holes inferred from recent kinematic determinations of black hole masses in nearby galaxies seem to indicate that the overall efficiency of supermassive black holes for producing blue light is lower than was previously assumed. We discuss three possible accretion modes with low optical emission efficiency: (i) accretion well above the Eddington rate, (ii) accretion obscured by dust, and (iii) accretion below the critical rate leading to an advection dominated accretion flow lasting for a Hubble time. We further argue that accretion with low optical efficiency might be closely related to the origin of the hard X-ray background.

The inner parts of an accretion disk around a spinning black hole are forced to align with the spin of the hole by the Bardeen-Petterson effect. Assuming that any jet produced by such a system is aligned with the angular momentum of either the hole or the inner disk, this can, in principle provide a mechanism for producing steady jets in AGN whose direction is independent of the angular momentum of the accreted material. However, the torque which aligns the inner disk with the hole, also, by Newton's third law, tends to align the spin of the hole with the outer accretion disk. In this letter, we calculate this alignment timescale for a black hole powering an AGN, and show that it is relatively short. This timescale is typically much less than the derived ages for jets in radio loud AGN, and implies that the jet directions are not in general controlled by the spin of the black hole. We speculate that the jet directions are most likely controlled either by the angular momentum of the accreted material or by the gravitational potential of the host galaxy.

Sunyaev-Zeldovich decrements in the temperature of the Cosmic Microwave Background (CMB) are produced by Thomson scattering of the CMB photons by the confined hot electrons in the Intra-Cluster Medium encountered along the line-of-sight. In this paper, we propose an alternate physical process that can produce detectable thermal and kinematic Sunyaev-Zeldovich decrements around quasars or quasar pairs in regions where no assembled, virialized clusters are detected either in the optical or X-ray wave-band. Invoking quasar outflows, we argue that both thermal and kinematic decrements can be produced by the baryons swept up by the mechanical luminosity of quasars. In contrast to the case of decrements produced by the hot electrons confined within the cluster potential, the magnitude of these quasar-outflow induced effects depend primarily on the luminosity of the QSO, and the physical scale of the outflow - which in turn depends both on the mean local density and the density gradient in the vicinity of the quasar. Quasar outflows produce (i) a frequency independent kinematic decrement ${\Delta T/T} \sim {10^{-4}}$; (ii) a frequency dependent fluctuation in the intensity ${\Delta {J_{\nu}}/{J_{\nu}}} \sim $ a 20% variation relative to ${\Delta T/T}$, say at 15 GHz and 30 GHz; (iii) a temperature increment ${\Delta T/T} \sim {10^{-4}}$ on the opposite side of the QSO; (iv) a linear polarization signal over 20'' at $\sim {10^{-7}}$ level; and therefore (v) patches of signal ${\Delta T/T} \sim {10^{-4}}$ on degree scales.

We present a new wide field image of the distant cluster AC114 ($z=0.31$) obtained with Wide Field and Planetary Camera II onboard the Hubble Space Telescope. This image considerably extends our knowledge of the lensing properties of the cluster beyond those derived by Smail \et (1995a) from a single WF/PC-1 pointing. In conjunction with published ground-based spectroscopy, we utilise several newly-discovered multiple images to construct an improved mass model for the central regions of the cluster. Using this model, we apply the methodology introduced by Natarajan & Kneib (1997) to interpret local perturbations to the cluster shear field on small scales resulting from mass associated with individual cluster galaxies. We use the lensing signal to place new constraints on the average mass-to-light ratio and spatial extents of the dark matter halos associated with morphologically-classified early-type cluster members. We find that the total mass of a fiducial $L^*$ cluster spheroidal galaxy is largely contained within $\sim$ 15 kpc radius halo ($\sim$~8--10 $R_e$) with a mass-to-light ratio ${M/L_V} \sim {15^{+10}_{-4}} $ (90% c.l., $h=0.5$) in solar units within this radius. Comparisons with similar estimates for fieldgalaxies indicates that the cluster galaxies in AC114 may possess less extensive and less massive halos. We discuss the consequences of our result in the context of models for the dynamical evolution of cluster galaxies and the observational prospects for extending this analysis.

We propose a scenario for the formation of a population of baryon-rich, dark matter-deficient dwarf galaxies at high redshift that form from the mass swept out in the Intergalactic Medium (IGM) by energetic outflows from luminous quasars. We predict the intrinsic properties of these galaxies, and examine the prospects for their observational detection in the optical, X-ray and radio wavebands. Detectable thermal Sunyaev-Zeldovich decrements (cold spots) on arc-minute scales in the cosmic microwave background radiation maps are expected during the shock-heated expanding phase from these hot bubbles. We conclude that the optimal detection strategy for these dwarfs is via narrow-band Lyman-$\alpha$ imaging of regions around high redshift quasars. An energetically scaled-down version of the same model is speculated upon as a possible mechanism for forming pre-galactic globular clusters.

We demonstrate that the luminosity function of the recently detected population of star-forming galaxies and the QSO luminosity function at z=3 can be matched with the mass function of dark matter haloes predicted by hierarchical cosmogonies for lifetimes of optically bright QSOs in the range 10^6 to 10^8 yr. We suggest that the mass of supermassive black holes may be limited by the back-reaction of the emitted energy on the accretion flow in a self-gravitating disc. This would imply a relation of black hole to halo mass of the form M_bh ~(v_halo)^5 and a typical duration of the optically bright QSO phase of a few times 10^7 yr. The high integrated mass density of black holes inferred from recent black hole mass estimates in nearby galaxies may indicate that the overall efficiency of supermassive black holes for producing blue light is smaller than previously assumed. We discuss three possible accretion modes with low optical emission efficiency: (i) accretion at far above the Eddington rate, (ii) accretion obscured by dust, and (iii) accretion below the critical rate leading to an advection dominated accretion flow lasting for a Hubble time. We further argue that accretion with low optical efficiency might be closely related to the origin of the hard X-ray background and that the ionizing background might be progressively dominated by stars rather than QSOs at higher redshift.

A gamma-ray burst (GRB) releases an amount of energy similar to that of a supernova explosion, which combined with its rapid variability suggests an origin related to neutron stars or black holes. Since these compact stellar remnants form from the most massive stars not long after their birth, gamma-ray bursts should trace the star formation rate in the Universe; we show that the GRB flux distribution is consistent with this. Because of the strong evolution of the star formation rate with redshift, it follows that the dimmest known bursts have z ~ 6, much above the value usually quoted and beyond the most distant quasars. This explains the absence of bright galaxies in well-studied gamma-ray burst error boxes. The increased distances imply a peak luminosity of 8.3e51 erg/s and a rate density of 0.025 per million years per galaxy. These values are 20 times higher and 150 times lower, respectively, than are implied by fits with non-evolving GRB rates. This means that GRBs are either caused by a much rarer phenomenon than mergers of binary neutron stars, or their gamma-ray emission is often invisible to us due to beaming. Precise burst locations from optical transients will discriminate between the various models for GRBs from stellar deaths, because the distance between progenitor birth place and burst varies greatly among them. The dimmest GRBs are then the most distant known objects, and may probe the Universe at an age when the first stars were forming.

We present the results of the application of the methodology introduced by Natarajan & Kneib (1997) to interpret local perturbations to the cluster shear field resulting from mass associated with individual cluster galaxies. The lensing signal is used to place new constraints on the average mass-to-light ratio and spatial extents of the dark matter halos associated with morphologically-classified early-type cluster members in the Abell cluster AC114. The total mass of a fiducial $L^\ast$ cluster spheroidal galaxy is found to be largely contained within $\sim$ 15 kpc radius halo ($\sim$~8--10 $R_e$) with a mass-to-light ratio ${M/L_V} \sim {15^{+10}_{-4}} $ (90 % c.l.) in solar units within this radius. Comparisons with similar estimates for field galaxies suggests that the cluster galaxies in AC114 may possess less extensive and less massive halos. Additionally, there is some indication that, at a fixed luminosity, S0 galaxies are less extended than ellipticals, suggesting a difference in the efficiency of tidal stripping for different galaxy types. These results enable us to probe the variation of the mass-to-light ratio with scale via gravitational lensing methods. Therefore, the prospects for constraining the mass density of the Universe and understanding galaxy evolution using these gravitational lensing methods are very promising.

We use the relationships between galactic HI mass and B-band luminosity determined by Rao & Briggs to recalculate the mass density of neutral gas at the present epoch based on more recent measures of the galaxy luminosity function than were available to those authors. We find $\Omega_{gas}(z=0) \simeq 5 \times 10^{-4}$ in good agreement with the original Rao & Briggs value, suggesting that this quantity is now reasonably secure. We then show that, if the scaling between H I mass and B-band luminosity has remained approximately constant since $z = 1$, the evolution of the luminosity function found by the Canada-France redshift survey translates to an increase of $\Omega_{gas}$ by a factor of $\approx 3$ at $z = 0.5 - 1$ . A similar value is obtained quite independently from consideration of the luminosity function of Mg II absorbers at $z = 0.65$. By combining these new estimates with data from damped \lya systems at higher redshift, it is possible to assemble a rough sketch of the evolution of $\Omega_{gas}$ over the last 90% of the age of the universe. The consumption of H I gas with time is in broad agreement with models of chemical evolution which include the effects of dust, although more extensive samples of damped \lya systems at low and intermediate redshift are required for a quantitative assessment of the dust bias.

Weak shear maps of the outer regions of clusters have been successfully used to map the distribution of mass at large radii from the cluster center. The typical smoothing lengths employed thus far preclude the systematic study of the effects of galactic-scale substructure on the measured weak lensing signal. In this paper, we present two methods to infer the possible existence and extent of dark halos around bright cluster galaxies by quantifying the `local' weak lensing induced by them. The proposed methods are: direct radial averaging of the shear field in the vicinity of bright cluster members and a maximum-likelihood method to extract fiducial parameters characterizing galaxy halos. We demonstrate using simulations that these observed local weak-shear effects on galaxy scales within the cluster can be used to statistically constrain the mean mass-to-light ratio, and fiducial parameters like the halo size, velocity dispersion and hence mass of cluster galaxies. We compare the two methods and investigate their relative drawbacks and merits in the context of feasibility of application to HST cluster data, whereby we find that the prospects are promising for detection on stacking a minimum of 20 WFPC2 deep cluster fields.

We present a useful analytic approximation to the solution of the Lane-Emden equation for infinite polytropic index - the isothermal sphere. The optimized expression obtained for the density profile can be accurate to within 0.04% within 5 core-radii and to 0.1% within 10 core-radii.

We present the results of an attempt to adapt the distribution function formalism to characterize large-scale structures like clusters of galaxies that form in a cosmological N-body simulation. While galaxy clusters are systems that are not strictly in equilibrium, we show that their evolution can nevertheless be studied using a physically motivated extension of the language of equilibrium stellar dynamics. Restricting our analysis to the virialized region, a prescription to limit the accessible phase-space is presented, which permits the construction of both the isotropic and the anisotropic distribution functions $f(E)$ and $f(E,L)$. The method is applied to models extracted from a catalogue of simulated clusters. Clusters evolved in open and flat background cosmologies are followed during the course of their evolution, and are found to transit through a sequence of what we define as `quasi-equilibrium' states. An interesting feature is that the computed $f(E)$ is well fit by an exponential form. We conclude that the dynamical evolution of a cluster, undergoing relaxation punctuated by interactions and violent mergers with consequent energy-exchange, can be studied both in a qualitative and quantitative fashion by following the time evolution of $f(E)$.

We propose a new approach to study the dynamical implications of mass models of clusters for the velocity structure of galaxies in the core. Strong and weak lensing data are used to construct the total mass profile of the cluster, which is used in conjunction with the optical galaxy data to solve in detail for the nature of galaxy orbits and the velocity anisotropy in the central regions. We also examine other observationally and physically motivated mass models, specifically those obtained from X-ray observations and N-body simulations. The aim of this analysis is to understand qualitatively the structure of the core and test some of the key assumptions of the standard picture of cluster formation regarding relaxation, virialization and equilibrium. This technique is applied to the cluster Abell 2218, where we find evidence for an anisotropic core, which we interpret to indicate the existence of a dynamically disturbed central region. We find that the requirement of physically meaningful solutions for the velocity anisotropy places stringent bounds on the slope of cluster density profiles in the inner regions.

Weak shear maps of the outer regions of clusters have been successfully used to map the distribution of mass at large radii. The effects of substructure in clusters on such reconstructions of the total mass have not been systematically studied. We propose a new method to study the effect of perturbers (bright cluster galaxies or sub-groups within the cluster) on the weak shear field. We present some analytic results below.

We report the results of an attempt to study the velocity structure of the core of the cluster A2218 in order to probe its dynamical state. The aim is to understand if the core is indeed virialized? relaxed? or in hydrodynamic equilibrium? We use cluster lensing data in conjunction with the optical galaxy data to solve in detail for the nature of galaxy orbits and the velocity anisotropy in the core. We present our preliminary results below.