Direct-collapse black holes (DCBHs) are currently one of the leading contenders for the origins of the first quasars in the universe, over 300 of which have now been found at $z >$ 6. But the birth of a DCBH in an atomically-cooling halo does not by itself guarantee it will become a quasar by $z \sim$ 7, the halo must also be located in cold accretion flows or later merge with a series of other gas-rich halos capable of fueling the BH's rapid growth. Here, we present near infrared luminosities for DCBHs born in cold accretion flows in which they are destined to grow to 10$^9$ M$_{\odot}$ by $z \sim$ 7. Our observables, which are derived from cosmological simulations with radiation hydrodynamics with Enzo, reveal that DCBHs could be found by the James Webb Space Telescope at $z \lesssim$ 20 and strongly-lensed DCBHs might be found in future wide-field surveys by Euclid and the Wide-Field Infrared Space Telescope at $z \lesssim$ 15.

Supermassive black holes (SMBHs) play a crucial role in the evolution of galaxies and are currently detected up to $z\sim 7.5$. Theories describing black hole (BH) growth are challenged by how rapidly seeds with initial mass $M_{\bullet} \lesssim 10^5 \, \mathrm{M_{\odot}}$, formed at $z \sim 20-30$, grew to $M_{\bullet} \sim 10^9 \, \mathrm{M_{\odot}}$ by $z\sim 7$. Here we study the effects of the value of the Hubble parameter, $H_0$, on models describing the early growth of BHs. First, we note that the predicted mass of a quasar at $z=6$ changes by $> 300$% if the underlying Hubble parameter used in the model varies from $H_0 = 65$ to $H_0 = 74$ km s$^{-1}$Mpc$^{-1}$, a range encompassing current estimates. Employing an MCMC approach based on priors from $z \gtrsim 6.5$ quasars and on $H_0$, we study the interconnection between $H_0$ and the parameters describing BH growth: seed mass $M_i$ and Eddington ratio $f_{\rm Edd}$. Assuming an Eddington ratio of $f_{\rm Edd} = 0.7$, in agreement with previous estimates, we find $H_0 = 73.6^{+1.2}_{-3.3}$ km s$^{-1}$Mpc$^{-1}$. In a second analysis, allowing all the parameters to vary freely, we find $\log(M_{i}/\mathrm{M_{\odot}}) > 4.5$ (at 95% CL), $H_0 = 74^{+1.5}_{-1.4}$ km s$^{-1}$Mpc$^{-1}$ and $f_{\rm Edd}=0.77^{+0.035}_{-0.026}$ at 68% CL. Our results on the typical Eddington ratio are in agreement with previous estimates. Current values of the Hubble parameter strongly favour heavy seed formation scenarios, with $M_i \gtrsim 10^4 \, \mathrm{M_{\odot}}$. In our model, with the priors on BH masses of quasars used, light seed formation scenarios are rejected at $\sim 3\sigma$.

Black holes across a broad range of masses play a key role in the evolution of galaxies. The initial seeds of black holes formed at $z \sim 30$ and grew over cosmic time by gas accretion and mergers. Using observational data for quasars and theoretical models for the hierarchical assembly of dark matter halos, we study the relative importance of gas accretion and mergers for black hole growth, as a function of redshift ($0<z<10$) and black hole mass ($ 10^3 \, \mathrm{M_{\odot}} < M_{\bullet} < 10^{10} \, \mathrm{M_{\odot}}$). We find that (i) growth by accretion is dominant in a large fraction of the parameter space, especially at $M_{\bullet} > 10^8 \, \mathrm{M_{\odot}}$ and $z>6$; and (ii) growth by mergers is dominant at $M_{\bullet} < 10^5 \, \mathrm{M_{\odot}}$ and $z>5.5$, and at $M_{\bullet} > 10^8 \, \mathrm{M_{\odot}}$ and $z<2$. As the growth channel has direct implications for the black hole spin (with gas accretion leading to higher spin values), we test our model against $\sim 20$ robust spin measurements available thus far. As expected, the spin tends to decline toward the merger-dominated regime, thereby supporting our model. The next generation of X-ray and gravitational-wave observatories (e.g. Lynx, AXIS, Athena and LISA) will map out populations of black holes up to very high redshift ($z\sim 20)$, covering the parameter space investigated here in almost its entirety. Their data will be instrumental to providing a clear picture of how black holes grew across cosmic time.

In the last two decades, approximately 200 quasars have been discovered at $z>6$, hosting active super-massive black holes with masses $M_{\bullet} \gtrsim 10^9 M_{\odot}$. While these sources reflect only the tip of the iceberg of the black hole mass distribution, their detection challenges standard growth models. The most massive $z>6$ black hole that was inferred thus far (J0100+2802, $M_{\bullet} \approx 1.2\times 10^{10} M_{\odot}$) was recently claimed to be lensed, with a magnification factor $\mu=450$. Here, we perform a consistency check of this claim, finding that the detection of such source requires a bright-end slope $\beta \geq 3.7$ for the intrinsic quasar luminosity function, $\Phi(L) \propto L^{-\beta}$. Commonly used values of $\beta \sim 2.8$ are rejected at $>3\sigma$. If the claim is confirmed, it is very unlikely that all the remaining $51$ sources in the SDSS sample are not magnified. Furthermore, it suffices that $\gtrsim 25\%$ of the remaining sources are lensed for the intrinsic luminosity function to differ significantly (i.e., $>3\sigma$) from the observed one. The presence of additional extremely magnified sources in the sample would lower the requirement to $\sim 4\%$. Our results urge the community to perform more extended multi-wavelength searches targeting $z>6$ lensed quasars, also among known samples. This effort could vitally contribute to solving the open problem of the growth of the brightest $z\sim 7$ quasars.

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

We present a multi wavelength analysis of 28 of the most luminous low-redshift narrow-line, ultra-hard X-ray selected active galactic nuclei (AGN) drawn from the 70 month Swift/BAT all-sky survey, with bolometric luminosities of log(L_bol/erg/s) > 45.25. The broad goal of our study is to determine whether these objects have any distinctive properties, potentially setting them aside from lower-luminosity obscured AGN in the local Universe. Our analysis relies on the first data release of the BAT AGN Spectroscopic Survey (BASS/DR1) and on dedicated observations with the VLT, Palomar, and Keck observatories. We find that the vast majority of our sources agree with commonly used AGN selection criteria which are based on emission line ratios and on mid-infrared colours. Our AGN are predominantly hosted in massive galaxies (9.8 < log(M_*/M_sun) < 11.7); based on visual inspection of archival optical images, they appear to be mostly ellipticals. Otherwise, they do not have distinctive properties. Their radio luminosities, determined from publicly available survey data, show a large spread of almost 4 orders of magnitude - much broader than what is found for lower X-ray luminosity obscured AGN in BASS. Moreover, our sample shows no preferred combination of black hole masses (M_BH) and/or Eddington ratio (lambda_Edd), covering 7.5 < log(M_BH/M_sun) < 10.3 and 0.01 < lambda_Edd < 1. Based on the distribution of our sources in the lambda_Edd-N_H plane, we conclude that our sample is consistent with a scenario where the amount of obscuring material along the line of sight is determined by radiation pressure exerted by the AGN on the dusty circumnuclear gas.

We study the coherence of the near-infrared and X-ray background fluctuations and the X-ray spectral properties of the sources producing it. We use data from multiple Spitzer and Chandra surveys, including the UDS/SXDF surveys, the Hubble Deep Field North, the EGS/AEGIS field, the Chandra Deep Field South and the COSMOS surveys, comprising $\sim$2275 Spitzer/IRAC hours and $\sim$~16 Ms of Chandra data collected over a total area of $\sim$~1~deg$^2$. We report an overall $\sim$5$\sigma$ detection of a cross-power signal on large angular scales $>$ 20$''$ between the 3.6 and 4.5\mum\ and the X-ray bands, with the IR vs [1-2] keV signal detected at 5.2$\sigma$. The [0.5-1] and [2-4] keV bands are correlated with the infrared wavelengths at a $\sim$1$-$3$\sigma$ significance level. The hardest X-ray band ([4-7] keV) alone is not significantly correlated with any infrared wavelengths due to poor photon and sampling statistics. We study the X-ray SED of the cross-power signal. We find that its shape is consistent with a variety of source populations of accreting compact objects, such as local unabsorbed AGNs or high-z absorbed sources. We cannot exclude that the excess fluctuations are produced by more than one population. Because of poor statistics, the current relatively broad photometric bands employed here do not allow distinguishing the exact nature of these compact objects or if a fraction of the fluctuations have instead a local origin.

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.

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.

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.

A major scientific goal of JWST is to probe the epoch of re-ionization of the Universe at z above 6, and up to 20 and beyond. At these redshifts, galaxies are just beginning to form and the observable objects are early black holes, supernovae, and cosmic infrared background. The JWST has the necessary sensitivity to observe these targets individually, but a public deep and wide science enabling survey in the wavelength range from 2-5 $\mu$m is needed to discover these black holes and supernovae and to cover the area large enough for cosmic infrared background to be reliably studied. This enabling survey will also discover a large number of other transients and enable sciences such as supernova cosmology up to z $\sim$ 5, star formation history at high redshift through supernova explosions, faint stellar objects in the Milky Way, and galaxy evolution up to z approaching 10. The results of this survey will also serve as an invaluable target feeder for the upcoming era of ELT and SKA.

The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization.

Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-sigma (velocity dispersion) relation. Local samples of active galactic nuclei (AGN) do not suffer from this selection effect, as in these samples Mbh is estimated from megamasers and/or reverberation mapping-based techniques. With the exception of megamasers, Mbh-estimates in these AGN samples are proportional to a virial coefficient fvir. Direct modelling of the broad line region suggests that fvir~3.5. However, this results in a Mbh-Mstar relation for AGN which lies below and is steeper than the one observed for quiescent black hole samples. A similar though milder trend is seen for the Mbh-sigma relation. Matching the high-mass end of the Mbh-Mstar and Mbh-sigma relations observed in quiescent samples requires fvir~15 and fvir~7, respectively. On the other hand, fvir~3.5 yields Mbh-sigma and Mbh-Mstar relations for AGN which are remarkably consistent with the expected `intrinsic' correlations for quiescent samples (i.e., once account has been made of the angular resolution-related selection effect), providing additional evidence that the sample of local quiescent black holes is biased. We also show that, as is the case for quiescent black holes, the Mbh-Mstar scaling relation of AGN is driven by velocity dispersion, thus providing additional key constraints to black hole-galaxy co-evolution models.

We improve the dynamical black hole (BH) mass estimates in three nearby low-mass early-type galaxies--NGC 205, NGC 5102, and NGC 5206. We use new \hst/STIS spectroscopy to fit the star formation histories of the nuclei in these galaxies, and use these measurements to create local color--mass-to-light ratio (\ml) relations. We then create new mass models from \hst~imaging and combined with adaptive optics kinematics, we use Jeans dynamical models to constrain their BH masses. The masses of the central BHs in NGC 5102 and NGC 5206 are both below one million solar masses and are consistent with our previous estimates, $9.12_{-1.53}^{+1.84}\times10^5$\Msun~and $6.31_{-2.74}^{+1.06}\times10^5$\Msun~(3$\sigma$ errors), respectively. However, for NGC 205, the improved models suggest the presence of a BH for the first time, with a best-fit mass of $6.8_{-6.7}^{+95.6}\times10^3$\Msun~(3$\sigma$ errors). This is the least massive central BH mass in a galaxy detected using any method. We discuss the possible systematic errors of this measurement in detail. Using this BH mass, the existing upper limits of both X-ray, and radio emissions in the nucleus of NGC 205 suggest an accretion rate $\lesssim$$10^{-5}$ of the Eddington rate. We also discuss the color--\mleff~relations in our nuclei and find that the slopes of these vary significantly between nuclei. Nuclei with significant young stellar populations have steeper color--\mleff~relations than some previously published galaxy color--\mleff~relations.

The discovery of the first strongly lensed $(\mu \approx 50)$ quasar at $z>6$ (J0439+1634) represents a breakthrough in our understanding of the early Universe. We derive the theoretical consequences of the new discovery. We predict that the observed population of $z > 6$ quasars should contain many sources with magnifications $\mu \lesssim 10$ and with image separations below the resolution threshold. Additionally, current selection criteria could have missed a substantial population of lensed $z > 6$ quasars, due to the contamination of the drop-out photometric bands by lens galaxies. We argue that this predicted population of lensed $z>6$ quasars would be misclassified and mixed up with low-$z$ galaxies. We quantify the fraction of undetected quasars as a function of the slope of the bright end of the quasar luminosity function, $\beta$. For $\beta \lesssim 3.6$, we predict that the undetected lensed quasars could reach half of the population, whereas for $\beta \gtrsim 4.5$ the vast majority of the $z > 6$ quasar population is lensed and still undetected. This would significantly affect the $z > 6$ quasar luminosity function and inferred black hole mass distributions, with profound implications for the ultraviolet, X-ray, and infrared cosmic backgrounds and the growth of early quasars.

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy.

The $M_{\bullet} - \sigma$ relation establishes a connection between central black holes (BHs) and their host spheroids. Supported by observations at $M_{\bullet} \gtrsim 10^5 \, M_{\odot}$, there is limited data on its validity at lower masses. Employing a semi-analytical model to simulate the combined evolution of BHs and their host galaxies, we predict the observational consequences of assuming a bimodality in the accretion efficiency of BHs, with low-mass BHs ($M_{\bullet} \lesssim 10^5 \, M_{\odot}$) accreting inefficiently. We predict a departure from the $M_{\bullet} - \sigma$ relation at a transitional BH mass $\sim 10^5 \, M_{\odot}$, with lower-mass BHs unable to reach the mass dictated by the relation and becoming disconnected from the evolution of the host galaxy. This prediction is an alternative to previous works suggesting a flattening of the relation at $\sim 10^5-10^6 \, M_{\odot}$. Furthermore, we predict a deficit of BHs shining at bolometric luminosities $\sim 10^{42} \, \mathrm{erg \, s^{-1}}$. Joined with a detection bias, this could partly explain the scarce number of intermediate-mass BHs detected. Conversely, we predict an increase in source density at lower bolometric luminosities, $<10^{42} \, \mathrm{erg \, s^{-1}}$. Because our predictions assume a bimodal population of high-redshift BH seeds, future observations of fainter BHs will be fundamental for constraining the nature of these seeds.

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.

JWST was conceived and built to answer one of the most fundamental questions that humans can address empirically: "How did the Universe make its first stars?". Our First Lights At REionization (FLARE) project transforms the quest for the epoch of reionization from the static to the time domain. It targets the complementary question: "What happened to those first stars?". It will be answered by observations of the most luminous events: supernovae and accretion on to black holes formed by direct collapse from the primordial gas clouds. These transients provide direct constraints on star-formation rates (SFRs) and the truly initial Initial Mass Function (IMF), and they may identify possible stellar seeds of supermassive black holes (SMBHs). Furthermore, our knowledge of the physics of these events at ultra-low metallicity will be much expanded. JWST's unique capabilities will detect these most luminous and earliest cosmic messengers easily in fairly shallow observations. However, these events are very rare at the dawn of cosmic structure formation and so require large area coverage. Time domain astronomy can be advanced to an unprecedented depth by means of a shallow field of JWST reaching 27 mag AB in 2 and 4.4 microns over a field as large as 0.1 square degree visited multiple times each year. Such a survey may set strong constraints or detect massive Pop III SNe at redshifts beyond 10, pinpointing the redshift of the first stars, or at least their death. Based on our current knowledge of superluminous supernovae (SLSNe), such a survey will find one or more SLSNe at redshifts above 6 in five years and possibly several direct collapse black holes. Although JWST is not designed as a wide field survey telescope, we show that such a wide field survey is possible with JWST and is critical in addressing several of its key scientific goals.

We use a semi-analytic model, ${\it Delphi}$, that jointly tracks the dark matter and baryonic assembly of high-redshift ($z \simeq 4-20$) galaxies to gain insight on the number density of Direct Collapse Black Hole (DCBH) hosts in three different cosmologies: the standard Cold Dark Matter (CDM) model and two Warm Dark Matter (WDM) models with particle masses of 3.5 and 1.5 keV. Obtaining the Lyman-Werner (LW) luminosity of each galaxy from ${\it Delphi}$, we use a clustering bias analysis to identify all, pristine halos with a virial temperature $T_{vir}>=10^4$ K that are irradiated by a LW background above a critical value as, DCBH hosts. In good agreement with previous studies, we find the DCBH number density rises from $\sim10^{-6.1}$ to $\sim 10^{-3.5}\, \mathrm{cMpc^{-3}}$ from $z\simeq 17.5$ to $8$ in the CDM model using a critical LW background value of $30 J_{21}$ (where $J_{21}= 10^{-21} \, {\rm erg\, s^{-1}\, Hz^{-1} \, cm^{-2} \, sr^{-1}}$). We find that a combination of delayed structure formation and an accelerated assembly of galaxies results in a later metal-enrichment and an accelerated build-up of the LW background in the 1.5 keV WDM model, resulting in DCBH hosts persisting down to much lower redshifts ($z \simeq 5$) as compared to CDM where DCBH hosts only exist down to $z \simeq 8$. We end by showing how the expected colours in three different bands of the Near Infrared Camera (NIRCam) onboard the forthcoming James Webb Space Telescope (${\it JWST}$) can be used to hunt for potential $z \simeq 5-9$ DCBHs, allowing hints on the WDM particle mass.

The detection of gravitational waves (GWs) generated by merging black holes has recently opened up a new observational window into the Universe. The mass of the black holes in the first and third LIGO detections, ($36-29 \, \mathrm{M_{\odot}}$ and $32-19 \, \mathrm{M_{\odot}}$), suggests low-metallicity stars as their most likely progenitors. Based on high-resolution N-body simulations, coupled with state-of-the-art metal enrichment models, we find that the remnants of Pop III stars are preferentially located within the cores of galaxies. The probability of a GW signal to be generated by Pop III stars reaches $\sim 90\%$ at $\sim 0.5 \, \mathrm{kpc}$ from the galaxy center, compared to a benchmark value of $\sim 5\%$ outside the core. The predicted merger rates inside bulges is $\sim 60 \times \beta_{III} \, \mathrm{Gpc^{-3} \, yr^{-1}}$ ($\beta_{III}$ is the Pop III binarity fraction). To match the $90\%$ credible range of LIGO merger rates, we obtain: $0.03 < \beta_{III} < 0.88$. Future advances in GW observatories and the discovery of possible electromagnetic counterparts could allow the localization of such sources within their host galaxies. The preferential concentration of GW events within the bulge of galaxies would then provide an indirect proof for the existence of Pop III stars.

Direct collapse black holes (DCBHs) are excellent candidates as seeds of supermassive black holes (SMBHs) observed at $z \gsim 6$. The formation of a DCBH requires a strong external radiation field to suppress $\rm H_2$ formation and cooling in a collapsing gas cloud. Such strong field is not easily achieved by first stars or normal star-forming galaxies. Here we investigate a scenario in which the previously-formed DCBH can provide the necessary radiation field for the formation of additional ones. Using one-zone model and the simulated DCBH Spectral Energy Distributions (SEDs) filtered through absorbing gas initially having column density $N_{\rm H}$, we derive the critical field intensity, $J_{\rm LW}^{\rm crit}$, to suppress $\rm H_2$ formation and cooling. For the SED model with $N_{\rm H}=1.3\times10^{25}$ cm$^{-2}$, $8.0\times10^{24}$ cm$^{-2}$ and $5.0\times10^{24}$ cm$^{-2}$, we obtain $J_{\rm LW}^{\rm crit}\approx22$, 35 and 54, all much smaller than the critical field intensity for normal star-forming galaxies ($J_{\rm LW}^{\rm crit}\simgt 1000$). X-ray photons from previously-formed DCBHs build up a high-$z$ X-ray background (XRB) that may boost the $J_{\rm LW}^{\rm crit}$. However, we find that in the three SED models $J_{\rm LW}^{\rm crit}$ only increases to $\approx80$, 170 and 390 respectively even when $\dt{\rho}_\bullet$ reaches the maximum value allowed by the present-day XRB level ($0.22, 0.034, 0.006~M_\odot$yr$^{-1}$Mpc$^{-3}$), still much smaller than the galactic value. Although considering the XRB from first galaxies may further increase $J_{\rm LW}^{\rm crit}$, we conclude that our investigation supports a scenario in which DCBH may be more abundant than predicted by models only including galaxies as external radiation sources.

Since the beginning of the new millennium, more than 100 $z\sim 6$ quasars have been discovered through several surveys and followed-up with multi-wavelength observations. These data provided a large amount of information on the growth of supermassive black holes at the early epochs, the properties of quasar host galaxies and the joint formation and evolution of these massive systems. We review the properties of the highest-$z$ quasars known so far, especially focusing on some of the most recent results obtained in (sub-)millimeter bands. We discuss key observational challenges and open issues in theoretical models and highlight possible new strategies to improve our understanding of the galaxy-black hole formation and evolution in the early Universe.

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 compare model results from a semi-analytic (merger-tree based) framework for high-redshift (z ~ 5-20) galaxy formation against reionization indicators, including the Planck electron scattering optical depth and the ionizing photon emissivity, to shed light on the reionization history and sources in Cold (CDM) and Warm Dark Matter (WDM; particle masses of $m_x = 1.5,3$ and 5 keV) cosmologies. This model includes all the key processes of star formation, supernova feedback, the merger/accretion/ejection driven evolution of gas and stellar mass and the effect of the ultra-violet background (UVB) in photo-evaporating the gas content of low-mass galaxies. We find that the delay in the start of reionization in light (1.5 keV) WDM models can be compensated by a steeper redshift evolution of the ionizing photon escape fraction and a faster mass assembly, resulting in reionization ending at comparable redshifts (z~5.5) in all the DM models considered. We find the bulk of the reionization photons come from galaxies with a halo mass $M_h < 10^9 M_\odot$ and a UV magnitude $ -15 < M_{UV} < -10$ in CDM. The progressive suppression of low-mass halos with decreasing $m_x$ leads to a shift in the reionization population to larger halo masses of $M_h > 10^9 M_\odot$ and $ -17 < M_{UV} < -13$ for 1.5 keV WDM. We find that current observations of the electron scattering optical depth and the Ultra-violet luminosity function are equally compatible with all the (cold and warm) DM models considered in this work. We propose that global indicators including the redshift evolution of the stellar mass density and the stellar mass-halo mass relation, observable with the James Webb Space Telescope, can be used to distinguish between CDM and WDM (1.5 keV) cosmologies.

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.

The peculiar emission properties of the $z \sim 6.6$ Ly$\alpha$ emitter CR7 have been initially interpreted with the presence of either a direct collapse black hole (DCBH) or a substantial mass of Pop III stars. Instead, updated photometric observations by Bowler et al. (2016) seem to suggest that CR7 is a more standard system. Here we confirm that the original DCBH hypothesis is consistent also with the new data. Using radiation-hydrodynamic simulations, we reproduce the new IR photometry with two models involving a Compton-thick DCBH of mass $\approx 7 \times 10^6 \, \mathrm{M_{\odot}}$ accreting (a) metal-free ($Z=0$) gas with column density $N_H = 8 \times 10^{25} \, \mathrm{cm^{-2}}$, or (b) low-metallicity gas ($Z = 5 \times 10^{-3} \, \mathrm{Z_{\odot}}$) with $N_H = 3 \times 10^{24} \, \mathrm{cm^{-2}}$. The best fit model reproduces the photometric data to within $1 \sigma$. Such metals can be produced by weak star-forming activity occurring after the formation of the DCBH. The main contribution to the Spitzer/IRAC $3.6 \, \mathrm{\mu m}$ photometric band in both models is due to HeI/HeII $\lambda 4714, 4687$ emission lines, while the contribution of [OIII] $\lambda 4959, 5007$ emission lines, if present, is sub-dominant. Spectroscopic observations with JWST will be required to ultimately clarify the nature of CR7.

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 first black hole seeds, formed when the Universe was younger than 500 Myr, are recognized to play an important role for the growth of early (z ~ 7) super-massive black holes. While progresses have been made in understanding their formation and growth, their observational signatures remain largely unexplored. As a result, no detection of such sources has been confirmed so far. Supported by numerical simulations, we present a novel photometric method to identify black hole seed candidates in deep multi-wavelength surveys. We predict that these highly-obscured sources are characterized by a steep spectrum in the infrared (1.6-4.5 micron), i.e. by very red colors. The method selects the only 2 objects with a robust X-ray detection found in the CANDELS/GOODS-S survey with a photometric redshift z > 6. Fitting their infrared spectra only with a stellar component would require unrealistic star formation rates (>2000 solar masses per year). To date, the selected objects represent the most promising black hole seed candidates, possibly formed via the direct collapse black hole scenario, with predicted mass >10^5 solar masses. While this result is based on the best photometric observations of high-z sources available to date, additional progress is expected from spectroscopic and deeper X-ray data. Upcoming observatories, like the JWST, will greatly expand the scope of this work.

Massive black holes (MBHs) are nowadays recognized as integral parts of galaxy evolution. Both the approximate proportionality between MBH and galaxy mass, and the expected importance of feedback from active MBHs in regulating star formation in their host galaxies point to a strong interplay between MBHs and galaxies. MBHs must form in the first galaxies and be fed by gas in these galaxies, with continuous or intermittent inflows that, at times, can be larger than the Eddington rate. Feedback from supernovae and from the MBHs themselves modulates the growth of the first MBHs. While current observational data only probe the most massive and luminous MBHs, the tip of the iceberg, we will soon be able to test theoretical models of MBH evolution on more "normal" MBHs: the MBHs that are indeed relevant in building the population that we observe in local galaxies, including our own Milky Way.

Massive Black Hole (MBH) seeds at redshift $z \gtrsim 10$ are now thought to be key ingredients to explain the presence of the super-massive ($10^{9-10} \, \mathrm{M_{\odot}}$) black holes in place $ < 1 \, \mathrm{Gyr}$ after the Big Bang. Once formed, massive seeds grow and emit copious amounts of radiation by accreting the left-over halo gas; their spectrum can then provide crucial information on their evolution. By combining radiation-hydrodynamic and spectral synthesis codes, we simulate the time-evolving spectrum emerging from the host halo of a MBH seed with initial mass $10^5 \, \mathrm{M_{\odot}}$, assuming both standard Eddington-limited accretion, or slim accretion disks, appropriate for super-Eddington flows. The emission occurs predominantly in the observed infrared-submm ($1-1000 \, \mathrm{\mu m}$) and X-ray ($0.1 - 100 \, \mathrm{keV}$) bands. Such signal should be easily detectable by JWST around $\sim 1 \, \mathrm{\mu m}$ up to $z \sim 25$, and by ATHENA (between $0.1$ and $10 \, \mathrm{keV}$, up to $z \sim 15$). Ultra-deep X-ray surveys like the Chandra Deep Field South could have already detected these systems up to $z \sim 15$. Based on this, we provide an upper limit for the $z \gtrsim 6$ MBH mass density of $\rho_{\bullet} \lesssim 2.5 \times 10^{2} \, \mathrm{M_{\odot} \, Mpc^{-3}}$ assuming standard Eddington-limited accretion. If accretion occurs in the slim disk mode the limits are much weaker, $\rho_{\bullet} \lesssim 7.6 \times 10^{3} \, \mathrm{M_{\odot} \, Mpc^{-3}}$ in the most constraining case.

CR7 is the brightest $z=6.6 \, {\rm Ly}\alpha$ emitter (LAE) known to date, and spectroscopic follow-up by Sobral et al. (2015) suggests that CR7 might host Population (Pop) III stars. We examine this interpretation using cosmological hydrodynamical simulations. Several simulated galaxies show the same "Pop III wave" pattern observed in CR7. However, to reproduce the extreme CR7 ${\rm Ly}\alpha$/HeII1640 line luminosities ($L_{\rm \alpha/He II}$) a top-heavy IMF and a massive ($>10^{7}{\rm M}_{\odot}$) PopIII burst with age $<2$ Myr are required. Assuming that the observed properties of ${\rm Ly}\alpha$ and HeII emission are typical for Pop III, we predict that in the COSMOS/UDS/SA22 fields, 14 out of the 30 LAEs at $z=6.6$ with $L_{\alpha} >10^{43.3}{\rm erg}\,{\rm s}^{-1}$ should also host Pop III stars producing an observable $L_{\rm He II}>10^{42.7}{\rm erg}\,{\rm s}^{-1}$. As an alternate explanation, we explore the possibility that CR7 is instead powered by accretion onto a Direct Collapse Black Hole (DCBH). Our model predicts $L_{\alpha}$, $L_{\rm He II}$, and X-ray luminosities that are in agreement with the observations. In any case, the observed properties of CR7 indicate that this galaxy is most likely powered by sources formed from pristine gas. We propose that further X-ray observations can distinguish between the two above scenarios.

The observational evidence that Super-Massive Black Holes ($M_{\bullet} \sim 10^{9-10} \, \mathrm{M_{\odot}}$) are already in place less than $1 \, \mathrm{Gyr}$ after the Big Bang poses stringent time constraints on the growth efficiency of their seeds. Among proposed possibilities, the formation of massive ($\sim 10^{3-6} \, \mathrm{M_{\odot}}$) seeds and/or the occurrence of super-Eddington ($\dot{M}>\dot{M}_{Edd}$) accretion episodes may contribute to the solution of this problem. In this work we analytically and numerically investigate the accretion flow onto high-redshift ($z \sim 10$) black holes to understand the physical requirements favoring rapid and efficient growth. Our model identifies a "feeding-dominated" accretion regime and a "feedback-limited" one, the latter being characterized by intermittent (duty cycles ${\cal D} \lesssim 0.5$) and inefficient growth, with recurring outflow episodes. We find that low-mass seeds ($\lesssim 10^{3-4} \, \mathrm{M_{\odot}}$) evolve in the feedback-limited regime, while more massive seeds ($\gtrsim 10^{5-6} \, \mathrm{M_{\odot}}$) grow very rapidly as they are found in the feeding-dominated regime. In addition to the standard accretion model with a fixed matter-energy conversion factor ($\epsilon = 0.1$), we have also explored slim disk models, appropriate for super-Eddington accretion, where radiation is trapped in the disk and the radiative efficiency is reduced ($\epsilon \lesssim 0.04$), which may ensure a continuous growth with $\dot{M} \gg \dot{M}_{Edd}$ (up to $\sim 300\dot{M}_{Edd}$ in our simulations). Under these conditions, outflows play a negligible role and a black hole can accrete $80\%-100\%$ of the gas mass of the host halo ($\sim 10^7 \, \mathrm{M_{\odot}}$) in $\sim 10 \, \mathrm{Myr}$, while in feedback-limited systems we predict that black holes can accrete only up to $\sim 15\%$ of the available mass.

We present a framework for high-redshift ($z \geq 7$) galaxy formation that traces their dark matter (DM) and baryonic assembly in four cosmologies: Cold Dark Matter (CDM) and Warm Dark Matter (WDM) with particle masses of $m_x =$ 1.5, 3 and 5 ${\rm keV}$. We use the same astrophysical parameters regulating star formation and feedback, chosen to match current observations of the evolving ultra violet luminosity function (UV LF). We find that the assembly of observable (with current and upcoming instruments) galaxies in CDM and $m_x \geq 3 {\rm keV}$ WDM results in similar halo mass to light ratios (M/L), stellar mass densities (SMDs) and UV LFs. However the suppression of small-scale structure leads to a notably delayed and subsequently more rapid stellar assembly in the $1.5 {\rm keV}$ WDM model. Thus galaxy assembly in $m_x \leq 2 {\rm keV}$ WDM cosmologies is characterized by: (i) a dearth of small-mass halos hosting faint galaxies; and (ii) a younger, more UV bright stellar population, for a given stellar mass. The higher M/L ratio (effect ii) partially compensates for the dearth of small-mass halos (effect i), making the resulting UV LFs closer to CDM than expected from simple estimates of halo abundances. We find that the redshift evolution of the SMD is a powerful probe of the nature of DM. Integrating down to a limit of $M_{UV} =-16.5$ for the James Webb Space Telescope (JWST), the SMD evolves as $\log$(SMD)$\propto -0.63 (1+z)$ in $m_x = 1.5 {\rm keV}$ WDM, as compared to $\log$(SMD)$\propto -0.44 (1+z)$ in CDM. Thus high-redshift stellar assembly provides a powerful testbed for WDM models, accessible with the upcoming JWST.

The possible formation of Direct Collapse Black Holes (DCBHs) in the first metal-free atomic cooling halos at high redshifts ($z > 10$) is nowadays object of intense study and several methods to prove their existence are currently under development. The abrupt collapse of a massive ($\sim 10^4 - 10^5 \, \mathrm{M_{\odot}}$) and rotating object is a powerful source of gravitational waves emission. In this work, we employ modern waveforms and the improved knowledge on the DCBHs formation rate to estimate the gravitational signal emitted by these sources at cosmological distances. Their formation rate is very high ($\sim 10^4 \, \mathrm{yr^{-1}}$ up to $z\sim20$), but due to a short duration of the collapse event ($\sim 2-30\, \mathrm{s}$, depending on the DCBH mass) the integrated signal from these sources is characterized by a very low duty-cycle (${\cal D}\sim 10^{-3}$), i.e. a shot-noise signal. Our results show that the estimated signal lies above the foreseen sensitivity of the Ultimate-DECIGO observatory in the frequency range $(0.8-300) \, \mathrm{mHz}$, with a peak amplitude $\Omega_{gw} = 1.1 \times 10^{-54}$ at $\nu_{max} = 0.9 \, \mathrm{mHz}$ and a peak Signal-to-Noise Ratio $\mathrm{SNR}\sim 22$ at $\nu = 20 \, \mathrm{mHz}$. This amplitude is lower than the Galactic confusion noise, generated by binary systems of compact objects in the same frequency band. For this reason, advanced techniques will be required to separate this signal from background and foreground noise components. As a proof-of-concept, we conclude by proposing a simple method, based on the auto-correlation function, to recognize the presence of a ${\cal D} \ll 1$ signal buried into the continuous noise. The aim of this work is to test the existence of a large population of high-z DCBHs, by observing the gravitational waves emitted during their infancy.

Theoretical models predict that a population of Intermediate Mass Black Holes (IMBHs) of mass $M_\bullet \approx 10^{4-5} \, \mathrm{M_{\odot}}$ might form at high ($z > 10$) redshift by different processes. Such objects would represent the seeds out of which $z > 6$ Super-Massive Black Holes (SMBHs) grow. We numerically investigate the radiation-hydrodynamic evolution governing the growth of such seeds via accretion of primordial gas within their parent dark matter halo of virial temperature $T_{vir} \sim 10^4 \, \mathrm{K}$. We find that the accretion onto a Direct Collapse Black Hole (DCBH) of initial mass $M_0=10^5 \, \mathrm{M_{\odot}}$ occurs at an average rate $\dot{M}_{\bullet} \simeq 1.35 \, \dot{M}_{Edd} \simeq 0.1 \, \mathrm{M_{\odot} \, yr^{-1}}$, is intermittent (duty-cycle $ < 50\%$) and lasts $\approx 142 \, \mathrm{Myr}$; the system emits on average at super-Eddington luminosities, progressively becoming more luminous as the density of the inner mass shells, directly feeding the central object, increases. Finally, when $\approx 90\%$ of the gas mass has been accreted (in spite of an average super-Eddington emission) onto the black hole, whose final mass is $\sim 7 \times 10^6 \, \mathrm{M_{\odot}}$, the remaining gas is ejected from the halo due to a powerful radiation burst releasing a peak luminosity $L_{peak}\sim 3\times 10^{45} \, \mathrm{erg \, s^{-1}}$. The IMBH is Compton-thick during most of the evolution, reaching a column density $N_H \sim 10^{25} \, \mathrm{cm^{-2}}$ in the late stages of the simulation. We briefly discuss the observational implications of the model.

We present a theoretical model embedding the essential physics of early galaxy formation (z = 5-12) based on the single premise that any galaxy can form stars with a maximal limiting efficiency that provides enough energy to expel all the remaining gas, quenching further star formation. This simple idea is implemented into a merger-tree based semi-analytical model that utilises two mass and redshift-independent parameters to capture the key physics of supernova feedback in ejecting gas from low-mass halos, and tracks the resulting impact on the subsequent growth of more massive systems via halo mergers and gas accretion. Our model shows that: (i) the smallest halos (halo mass $M_h \leq 10^{10} M_\odot$) build up their gas mass by accretion from the intergalactic medium; (ii) the bulk of the gas powering star formation in larger halos ($M_h \geq 10^{11.5} M_\odot$) is brought in by merging progenitors; (iii) the faint-end UV luminosity function slope evolves according to $\alpha = -1.75 \log \,z -0.52$. In addition, (iv) the stellar mass-to-light ratio is well fit by the functional form $\log\, M_* = -0.38 M_{UV} -0.13\, z + 2.4$, which we use to build the evolving stellar mass function to compare to observations. We end with a census of the cosmic stellar mass density (SMD) across galaxies with UV magnitudes over the range $-23 \leq M_{UV} \leq -11$ spanning redshifts $5 < z < 12$: (v) while currently detected LBGs contain $\approx 50$% (10%) of the total SMD at $z=5$ (8), the JWST will detect up to 25% of the SMD at $z \simeq 9.5$.

The cosmological 21cm signal is a physics-rich probe of the early Universe, encoding information about both the ionization and the thermal history of the intergalactic medium (IGM). The latter is likely governed by X-rays from star-formation processes inside very high redshift (z > 15) galaxies. Due to the strong dependence of the mean free path on the photon energy, the X-ray SED can have a significant impact on the interferometric signal from the cosmic dawn. Recent Chandra observations of nearby, star-forming galaxies show that their SEDs are more complicated than is usually assumed in 21cm studies. In particular, these galaxies have ubiquitous, sub-keV thermal emission from the hot interstellar medium (ISM), which generally dominates the soft X-ray luminosity (with energies < 1 keV, sufficiently low to significantly interact with the IGM). Using illustrative soft and hard SEDs, we show that the IGM temperature fluctuations in the early Universe would be substantially increased if the X-ray spectra of the first galaxies were dominated by the hot ISM, compared with X-ray binaries with harder spectra. The associated large-scale power of the 21cm signal would be higher by roughly a factor of three. More generally, we show that the peak in the redshift evolution of the large-scale (k = 0.2 1/Mpc) 21cm power is a robust probe of the soft-band SED of the first galaxies, and importantly, is not degenerate with their bolometric luminosities. On the other hand, the redshift of the peak constrains the X-ray luminosity and halo masses which host the first galaxies.

Recent observations have shown the presence of extra-solar planets in Galactic open stellar clusters, as in the Praesepe (M44). These systems provide a favorable environment for planetary formation due to the high heavy-element content exhibited by the majority of their population. The large stellar density, and corresponding high close-encounter event rate, may induce strong perturbations of planetary orbits with large semimajor axes. Here we present a set of N-body simulations implementing a novel scheme to treat the tidal effects of external stellar perturbers on planetary orbit eccentricity and inclination. By simulating five nearby open clusters we determine the rate of occurrence of bodies extracted from their parent stellar system by quasi-impulsive tidal interactions. We find that the specific free-floating planet production rate (total number of free-floating planets per unit of time, normalized by the total number of stars) is proportional to the stellar density of the cluster, with a constant of proportionality equal to (23 +/- 5)10^-6 pc^3 Myr^-1. For the Pleiades (M45) we predict that about 26% of stars should have lost their planets. This raises the exciting possibility of directly observing these wandering planets with the James Webb Space Telescope in the NIR band. Assuming a surface temperature of the planet of 500 K, a free-floating planet of Jupiter size inside the Pleiades would have a specific flux @4.4 micron of approximately 400 nJy, which would lead to a very clear detection (S/N of order 100) in only one hour of integration.

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