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

The standard cosmological model ($\Lambda$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$\alpha$ detection, motivating the development of modern theoretical predictions. We use hydrodynamical cosmological simulations run with the AREPO code to investigate the potential detectability of large-scale filaments, excluding contributions from the halos embedded in them. We focus on filaments connecting massive ($M_{200c}\sim(1-3)\times10^{12} M_\odot$) halos at z=3, and compare different simulation resolutions, feedback levels, and mock-image pixel sizes. We find increasing simulation resolution does not substantially improve detectability notwithstanding the intrinsic enhancement of internal filament structure. By contrast, for a MUSE integration of 31 hours, including feedback increases the detectable area by a factor of $\simeq$5.5 on average compared with simulations without feedback, implying that even the non-bound components of the filaments have substantial sensitivity to feedback. Degrading the image resolution from the native MUSE scale of (0.2")$^2$ per pixel to (5.3")$^2$ apertures has the strongest effect, increasing the detectable area by a median factor of $\simeq$200 and is most effective when the size of the pixel roughly matches the width of the filament. Finally, we find the majority of Ly$\alpha$ emission is due to electron impact collisional excitations, as opposed to radiative recombination.

The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status.

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

Population III stars forming in minihalos tend to be relatively inefficient, with each minihalo hosting one or a small number of stars which are more massive than local stars but still challenging to observe directly at high redshift. Here we explore a possible mechanisms for the generation of larger clusters of such stars: a nearby ionizing source which ionizes a late forming halo, delaying its collapse until the halo is sufficiently large that the core can self-shield and suffer runaway collapse. We use simulations with a simple but accurate model for the radiative ionizing flux and confirm the basic predictions of previous work: higher ionizing fluxes can delay the collapse to lower redshifts and higher masses, up to an order of magnitude above the atomic cooling limit. In a limited number of runs we also examine the fragmentation of the cores at even higher resolution, using both simple estimates and sink particles to show that the number of fragments is generally small, at most a handful, and that the mass accretion rate on the fragments is of order $10^{-3}$ M$_\odot$ yr$^{-1}$. This rate is sufficiently high that the descent on the main sequence (and hence the suppression of accretion) is delayed until the stellar masses are of order $100-1000$ M$_\odot$, but not high enough to produce direct collapse black holes of mass $\sim 10^5$ M$_{\odot}$. The resulting clusters are larger than those produced in minihalos, but are still likely to fall short of being easily detectable in JWST blind fields.

The Cosmic Dawn and Reionization epochs remain a fundamental but challenging frontier of astrophysics and cosmology. We advocate a large-scale, multi-tracer approach to develop a comprehensive understanding of the physics that led to the formation and evolution of the first stars and galaxies. We highlight the line intensity mapping technique to trace the multi-phase reionization topology on large scales, and measure reionization history in detail. Besides 21cm, we advocate for Lya tomography mapping during the epoch of Wouthuysen-Field coupling as an additional probe of the cosmic dawn era.

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.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

In addition to being a probe of Cosmic Dawn and Epoch of Reionization astrophysics, the 21cm line at $z>6$ is also a powerful way to constrain cosmology. Its power derives from several unique capabilities. First, the 21cm line is sensitive to energy injections into the intergalactic medium at high redshifts. It also increases the number of measurable modes compared to existing cosmological probes by orders of magnitude. Many of these modes are on smaller scales than are accessible via the CMB, and moreover have the advantage of being firmly in the linear regime (making them easy to model theoretically). Finally, the 21cm line provides access to redshifts prior to the formation of luminous objects. Together, these features of 21cm cosmology at $z>6$ provide multiple pathways toward precise cosmological constraints. These include the "marginalizing out" of astrophysical effects, the utilization of redshift space distortions, the breaking of CMB degeneracies, the identification of signatures of relative velocities between baryons and dark matter, and the discovery of unexpected signs of physics beyond the $\Lambda$CDM paradigm at high redshifts.

The "cosmic dawn" refers to the period of the Universe's history when stars and black holes first formed and began heating and ionizing hydrogen in the intergalactic medium (IGM). Though exceedingly difficult to detect directly, the first stars and black holes can be constrained indirectly through measurements of the cosmic 21-cm background, which traces the ionization state and temperature of intergalactic hydrogen gas. In this white paper, we focus on the science case for such observations, in particular those targeting redshifts z $\gtrsim$ 10 when the IGM is expected to be mostly neutral. 21-cm observations provide a unique window into this epoch and are thus critical to advancing first star and black hole science in the next decade.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

Cosmic dawn and the Epoch of Reionization (EoR) are among the least explored observational eras in cosmology: a time at which the first galaxies and supermassive black holes formed and reionized the cold, neutral Universe of the post-recombination era. With current instruments, only a handful of the brightest galaxies and quasars from that time are detectable as individual objects, due to their extreme distances. Fortunately, a multitude of multi-wavelength intensity mapping measurements, ranging from the redshifted 21 cm background in the radio to the unresolved X-ray background, contain a plethora of synergistic information about this elusive era. The coming decade will likely see direct detections of inhomogenous reionization with CMB and 21 cm observations, and a slew of other probes covering overlapping areas and complementary physical processes will provide crucial additional information and cross-validation. To maximize scientific discovery and return on investment, coordinated survey planning and joint data analysis should be a high priority, closely coupled to computational models and theoretical predictions.

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.

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

We present the first simulations of the high-redshift Ly$\alpha$ intensity field that account for scattering in the intergalactic medium (IGM). Using a 3D Monte Carlo radiative transfer code, we find that Ly$\alpha$ scattering smooths spatial fluctuations in the Ly$\alpha$ intensity on small scales and that the spatial dependence of this smoothing depends strongly on the mean neutral fraction of the IGM. Our simulations find a strong effect of reionization on $k=0.1-1~{\rm Mpc^{-1}}$, with $P_{\rm Ly\alpha}\propto k^{-1.75}$ for $\bar{x}_{\rm HI} = 0.63$ and $P_{\rm Ly\alpha} \propto k^{-2.2}$ for $\bar{x}_{\rm HI} = 0.86$ in contrast to $P_{\rm Ly\alpha}\propto k^{-1.5}$ after reionization. At wavenumbers of $k>1 ~ {\rm Mpc^{-1}}$, we find that the signal is sensitive to the emergent Ly$\alpha$ line profiles from galaxies. We also demonstrate that the cross-correlation between a Ly$\alpha$ intensity map and a future galaxy redshift survey could be detected on large scales by an instrument similar to SPHEREx, and over a wide range of scales by a hypothetical intensity mapping instrument in the vein of CDIM.

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

Following the first two annual intensity mapping workshops at Stanford in March 2016 and Johns Hopkins in June 2017, we report on the recent advances in theory, instrumentation and observation that were presented in these meetings and some of the opportunities and challenges that were identified looking forward. With preliminary detections of CO, [CII], Lya and low-redshift 21cm, and a host of experiments set to go online in the next few years, the field is rapidly progressing on all fronts, with great anticipation for a flood of new exciting results. This current snapshot provides an efficient reference for experts in related fields and a useful resource for nonspecialists. We begin by introducing the concept of line-intensity mapping and then discuss the broad array of science goals that will be enabled, ranging from the history of star formation, reionization and galaxy evolution to measuring baryon acoustic oscillations at high redshift and constraining theories of dark matter, modified gravity and dark energy. After reviewing the first detections reported to date, we survey the experimental landscape, presenting the parameters and capabilities of relevant instruments such as COMAP, mmIMe, AIM-CO, CCAT-p, TIME, CONCERTO, CHIME, HIRAX, HERA, STARFIRE, MeerKAT/SKA and SPHEREx. Finally, we describe recent theoretical advances: different approaches to modeling line luminosity functions, several techniques to separate the desired signal from foregrounds, statistical methods to analyze the data, and frameworks to generate realistic intensity map simulations.

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

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

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

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

We propose high-velocity collisions of protogalaxies as a new pathway to form supermassive stars (SMSs) with masses of ~ 10^5 Msun at high redshift (z > 10). When protogalaxies hosted by dark matter halos with a virial temperature of ~ 10^4 K collide with a relative velocity > 200 km/s, the gas is shock-heated to ~ 10^6 K and subsequently cools isobarically via free-free emission and He^+, He, and H line emission. Since the gas density (> 10^4 cm^{-3}) is high enough to destroy H_2 molecules by collisional dissociation, the shocked gas never cools below ~ 10^4 K. Once a gas cloud of ~ 10^5 Msun reaches this temperature, it becomes gravitationally unstable and forms a SMS which will rapidly collapse into a super massive black hole (SMBH) via general relativistic instability. We perform a simple analytic estimate of the number density of direct-collapse black holes (DCBHs) formed through this scenario (calibrated with cosmological N-body simulations) and find n_{DCBH} ~ 10^{-9} Mpc^{-3} (comoving) by z = 10. This could potentially explain the abundance of bright high-z quasars.

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

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

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

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

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

Models and simulations of the epoch of reionization predict that spectra of the 21-cm transition of atomic hydrogen will show a clear fluctuation peak, at a redshift and scale, respectively, that mark the central stage of reionization and the characteristic size of ionized bubbles. This is based on the assumption that the cosmic gas was heated by stellar remnants - particularly X-ray binaries - to temperatures well above the cosmic microwave background at that time (~ 30 K). Here we show instead that the hard spectra (that is, spectra with more high-energy photons than low-energy photons) of X-ray binaries make such heating ineffective, resulting in a delayed and spatially uniform heating that modifies the 21-cm signature of reionization. Rather than looking for a simple rise and fall of the large-scale fluctuations (peaking at several millikelvin), we must expect a more complex signal also featuring a distinct minimum (at less than a millikelvin) that marks the rise of the cosmic mean gas temperature above the microwave background. Observing this signal, possibly with radio telescopes in operation today, will demonstrate the presence of a cosmic background of hard X-rays at that early time.

We present the first complete calculation of the history of the inhomogeneous 21-cm signal from neutral hydrogen during the era of the first stars. We use hybrid computational methods to capture the large-scale distribution of the first stars, whose radiation couples to the neutral hydrogen emission, and to evaluate the 21-cm signal from z ~ 15-35. In our realistic picture large-scale fluctuations in the 21-cm signal are sourced by the inhomogeneous density field and by the Ly-alpha and X-ray radiative backgrounds. The star formation is suppressed by two spatially varying effects: negative feedback provided by the Lyman-Werner radiative background, and supersonic relative velocities between the gas and dark matter. Our conclusions are quite promising: we find that the fluctuations imprinted by the inhomogeneous Ly-alpha background in the 21-cm signal at z ~ 25 should be detectable with the Square Kilometer Array.

The formation of the first stars is an exciting frontier area in astronomy. Early redshifts z ~ 20 have become observationally promising as a result of a recently recognized effect of a supersonic relative velocity between the dark matter and gas. This effect produces prominent structure on 100 comoving Mpc scales, which makes it much more feasible to detect 21-cm fluctuations from the epoch of first heating. We use semi-numerical hybrid methods to follow for the first time the joint evolution of the X-ray and Lyman-Werner radiative backgrounds, including the effect of the supersonic streaming velocity on the cosmic distribution of stars. We incorporate self-consistently the negative feedback on star formation induced by the Lyman-Werner radiation, which dissociates molecular hydrogen and thus suppresses gas cooling. We find that the feedback delays the X-ray heating transition by a Delta z ~ 2, but leaves a promisingly large fluctuation signal over a broad redshift range. The large-scale power spectrum is predicted to reach a maximal signal-to-noise ratio of S/N ~ 3-4 at z ~ 18 (for a projected first-generation instrument), with S/N > 1 out to z ~ 22-23. We hope to stimulate additional numerical simulations as well as observational efforts focused on the epoch prior to cosmic reionization.

Dark and baryonic matter moved at different velocities in the early Universe, which strongly suppressed star formation in some regions. This was estimated to imprint a large-scale fluctuation signal of about 2 mK in the 21-cm spectral line of atomic hydrogen associated with stars at a redshift of 20, although this estimate ignored the critical contribution of gas heating due to X-rays and major enhancements of the suppression. A large velocity difference reduces the abundance of halos and requires the first stars to form in halos of about a million solar masses, substantially greater than previously expected. Here we report a simulation of the distribution of the first stars at z=20 (cosmic age of ~180 Myr), incorporating all these ingredients within a 400 Mpc box. We find that the 21-cm signature of these stars is an enhanced (10 mK) fluctuation signal on the 100-Mpc scale, characterized by a flat power spectrum with prominent baryon acoustic oscillations. The required sensitivity to see this signal is achievable with an integration time of a thousand hours with an instrument like the Murchison Wide-field Array or the Low Frequency Array but designed to operate in the range of 50-100 MHz.

We explore the possibility that the density profiles of elliptical galaxies and cold dark matter (CDM) halos found in cosmological simulations can be understood in terms of the same physical process, collisionless gravitational collapse. To investigate this, we study a simplified model, the collapse of a perfectly cold Plummer sphere. First, we examine an N-body simulation of this model with particles constrained to purely radial orbits. This results in a final state characterized by a profile slightly steeper than \rho \propto r^{-2} at small radii and behaving as \rho \propto r^{-4} at large radii, which can be understood in terms of simple analytic arguments. Next, we repeat our simulation without the restriction of radial orbits. This results in a shallower inner density profile, like those found in elliptical galaxies and CDM halos. We attribute this change to the radial orbit instability (ROI) and propose a form of the distribution function (DF) motivated by a physical picture of collapse. As evidence of the link between our model and CDM halos, we find that our collapse simulation has a final state with pseudo-phase-space density which scales roughly as \rho/\sigma^3 \propto r^{-1.875}, like that observed in CDM halos from cosmological simulations (Navarro et al. 2010). The velocity anisotropy profile is also qualitatively similar to that found near the centers of these halos. We argue that the discrepancy at large radii (where CDM halos scale as \rho \propto r^{-3}) is due to the presence of the cosmological background or continued infall. This leads us to predict that the outer CDM halo density profile is not "universal," but instead depends on cosmological environment (be it an underdense void or overdense region).

Measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect from instruments such as the South Pole Telescope (SPT) and the Atacama Cosmology Telescope (ACT) will soon put improved constraints on reionization. Popular models assume that UV photons alone are responsible for reionization of the intergalactic medium. We explore the effects of a significant contribution of X-rays to reionization on the kSZ signal. Because X-rays have a large mean free path through the neutral intergalactic medium, they introduce partial ionization in between the sharp-edged bubbles created by UV photons. This smooth ionization component changes the power spectrum of the cosmic microwave background (CMB) temperature anisotropies. We quantify this effect by running semi-numerical simulations of reionization. We test a number of different models of reionization without X-rays that have varying physical parameters, but which are constrained to have similar total optical depths to electron scattering. These are then compared to models with varying levels of contribution to reionization from X-rays. We find that models with more than a 10% contribution from X-rays produce a significantly lower power spectrum of temperature anisotropies than all the UV-only models tested. The expected sensitivity of SPT and ACT may be insufficient to distinguish between our models, however, a non-detection of the kSZ signal from the epoch of reionization could result from the contribution of X-rays. It will be important for future missions with improved sensitivity to consider the impact of X-ray sources on reionization.

Visbal & Loeb (2010) have shown that it is possible to measure the clustering of galaxies by cross correlating the cumulative emission from two different spectral lines which originate at the same redshift. Through this cross correlation, one can study galaxies which are too faint to be individually resolved. This technique, known as intensity mapping, is a promising probe of the global properties of high redshift galaxies. Here, we test the feasibility of such measurements with synthetic data generated from cosmological dark matter simulations. We use a simple prescription for associating galaxies with dark matter halos and create a realization of emitted radiation as a function of angular position and wavelength over a patch of the sky. This is then used to create synthetic data for two different hypothetical instruments, one aboard the Space Infrared Telescope for Cosmology and Astrophysics (SPICA) and another consisting of a pair of ground based radio telescopes designed to measure the CO(1-0) and CO(2-1) emission lines. We find that the line cross power spectrum can be measured accurately from the synthetic data with errors consistent with the analytical prediction of Visbal & Loeb (2010). Removal of astronomical backgrounds and masking bright line emission from foreground contaminating galaxies do not prevent accurate cross power spectrum measurements.

We discuss a method for detecting the emission from high redshift galaxies by cross correlating flux fluctuations from multiple spectral lines. If one can fit and subtract away the continuum emission with a smooth function of frequency, the remaining signal contains fluctuations of flux with frequency and angle from line emitting galaxies. Over a particular small range of observed frequencies, these fluctuations will originate from sources corresponding to a series of different redshifts, one for each emission line. It is possible to statistically isolate the fluctuations at a particular redshift by cross correlating emission originating from the same redshift, but in different emission lines. This technique will allow detection of clustering fluctuations from the faintest galaxies which individually cannot be detected, but which contribute substantially to the total signal due to their large numbers. We describe these fluctuations quantitatively through the line cross power spectrum. As an example of a particular application of this technique, we calculate the signal-to-noise ratio for a measurement of the cross power spectrum of the OI(63 micron) and OIII(52 micron) fine structure lines with the proposed Space Infrared Telescope for Cosmology and Astrophysics. We find that the cross power spectrum can be measured beyond a redshift of z=8. Such observations could constrain the evolution of the metallicity, bias, and duty cycle of faint galaxies at high redshifts and may also be sensitive to the reionization history through its effect on the minimum mass of galaxies. As another example, we consider the cross power spectrum of CO line emission measured with a large ground based telescope like CCAT and 21-cm radiation originating from hydrogen in galaxies after reionization with an interferometer similar in scale to MWA, but optimized for post-reionization redshifts.

21cm emission from residual neutral hydrogen after the epoch of reionization can be used to trace the cosmological power spectrum of density fluctuations. Using a Fisher matrix formulation, we provide a detailed forecast of the constraints on cosmological parameters that are achievable with this probe. We consider two designs: a scaled-up version of the MWA observatory as well as a Fast Fourier Transform Telescope. We find that 21cm observations dedicated to post-reionization redshifts may yield significantly better constraints than next generation Cosmic Microwave Background (CMB) experiments. We find the constraints on $\Omega_\Lambda$, $\Omega_{\rm m}h^2$, and $\Omega_\nu h^2$ to be the strongest, each improved by at least an order of magnitude over the Planck CMB satellite alone for both designs. Our results do not depend as strongly on uncertainties in the astrophysics associated with the ionization of hydrogen as similar 21cm surveys during the epoch of reionization. However, we find that modulation of the 21cm power spectrum from the ionizing background could potentially degrade constraints on the spectral index of the primordial power spectrum and its running by more than an order of magnitude. Our results also depend strongly on the maximum wavenumber of the power spectrum which can be used due to non-linearities.

The UV radiation from a quasar leaves a characteristic pattern in the distribution of ionized hydrogen throughout the surrounding space. This pattern or light echo propagates through the intergalactic medium at the speed of light, and can be observed by its imprint on the Ly-alpha forest spectra of background sources. As the echo persists after the quasar has switched off, it offers the possibility of searching for dead quasars, and constraining their luminosities and lifetimes. We outline a technique to search for and characterize these light echoes. To test the method, we create artificial Ly-alpha forest spectra from cosmological simulations at z=3, apply light echoes and search for them. We show how the simulations can also be used to quantify the significance level of any detection. We find that light echoes from the brightest quasars could be found in observational data. With absorption line spectra of 100 redshift z~3-3.5 quasars or galaxies in a 1 square degree area, we expect that ~10 echoes from quasars with B band luminosities L_B=3x10^45 ergs/s exist that could be found at 95% confidence, assuming a quasar lifetime of ~10^7 yr. Even a null result from such a search would have interesting implications for our understanding of quasar luminosities and lifetimes.

Air fluorescence detectors measure the energy of ultra-high energy cosmic rays by collecting fluorescence light emitted from nitrogen molecules along the extensive air shower cascade. To ensure a reliable energy determination, the light signal needs to be corrected for atmospheric effects, which not only attenuate the signal, but also produce a non-negligible background component due to scattered Cherenkov light and multiple-scattered light. The correction requires regular measurements of the aerosol attenuation length and the aerosol phase function, defined as the probability of light scattered in a given direction. At the Pierre Auger Observatory in Malargue, Argentina, the phase function is measured on an hourly basis using two Aerosol Phase Function (APF) light sources. These sources direct a UV light beam across the field of view of the fluorescence detectors; the phase function can be extracted from the image of the shots in the fluorescence detector cameras. This paper describes the design, current status, standard operation procedure, and performance of the APF system at the Pierre Auger Observatory.

The air fluorescence detectors (FDs) of the Pierre Auger Observatory are vital for the determination of the air shower energy scale. To compensate for variations in atmospheric conditions that affect the energy measurement, the Observatory operates an array of monitoring instruments to record hourly atmospheric conditions across the detector site, an area exceeding 3,000 square km. This paper presents results from four instruments used to characterize the aerosol component of the atmosphere: the Central Laser Facility (CLF), which provides the FDs with calibrated laser shots; the scanning backscatter lidars, which operate at three FD sites; the Aerosol Phase Function monitors (APFs), which measure the aerosol scattering cross section at two FD locations; and the Horizontal Attenuation Monitor (HAM), which measures the wavelength dependence of aerosol attenuation.