Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

We investigate the [CII] line intensity mapping (IM) signal from galaxies in the Epoch of Reionization (EoR) to assess its detectability, the possibility to constrain the $L_{\rm CII}-{\rm SFR}$ relation, and to recover the [CII] luminosity function (LF) from future experiments. By empirically assuming that ${\rm log}L_{\rm CII}={\rm log}A+\gamma {\rm SFR}\pm\sigma_L$, we derive the [CII] LF from the observed UV LF, and the [CII] IM power spectrum. We study the shot-noise and the full power spectrum separately. Although, in general, the shot-noise component has a much higher signal-to-noise ratio than the clustering one, it cannot be used to put independent constraints on log$A$ and $\gamma$. Full power spectrum measurements are crucial to break such degeneracy, and reconstruct the [CII] LF. In our fiducial survey S1 (inspired by CCAT-p/1000 hr) at $z \sim 6$, the shot-noise (clustering) signal is detectable for 2 (1) of the 5 considered $L_{\rm CII}-{\rm SFR}$ relations. The shot-noise is generally dominated by galaxies with $L_{\rm CII}\gtrsim 10^{8-9}~L_\odot$ ($M_{\rm UV}\sim-20$ to $-22$), already at reach of ALMA pointed observations. However, given the small field of view of such telescope, an IM experiment would provide unique information on the bright-end of the LF. The detection depth of an IM experiment crucially depends on the (poorly constrained) $L_{\rm CII}-{\rm SFR}$ relation in the EoR. If the $L_{\rm CII}-{\rm SFR}$ relation varies in a wide log$A - \gamma$ range, but still consistent with ALMA [CII] LF upper limits, even the signal from galaxies with $L_{\rm CII}$ as faint as $\sim 10^7~L_\odot$ could be detectable. Finally, we consider the contamination by continuum foregrounds (CIB, dust, CMB) and CO interloping lines, and derived the requirements on the residual contamination level to reliably extract the [CII] signal.

The neutral hydrogen (HI) and its 21 cm line are promising probes to the reionization process of the intergalactic medium (IGM). To use this probe effectively, it is imperative to have a good understanding on how the neutral hydrogen traces the underlying matter distribution. Here we study this problem using semi-numerical modeling by combining the HI in the IGM and the HI from halos during the epoch of reionization (EoR), and investigate the evolution and the scale-dependence of the neutral fraction bias as well as the 21 cm line bias. We find that the neutral fraction bias on large scales is negative during reionization, and its absolute value on large scales increases during the early stage of reionization and then decreases during the late stage. During the late stage of reionization, there is a transition scale at which the HI bias transits from negative on large scales to positive on small scales, and this scale increases as the reionization proceeds to the end.

Due to ionosphere absorption and the interference by natural and artificial radio emissions, ground observation of the sky at the decameter or longer is very difficult. This unexplored part of electromagnetic spectrum has the potential of great discoveries, notably in the study of cosmic dark ages and dawn, but also in heliophysics and space weather, planets, cosmic ray and neutrinos, pulsar and interstellar medium, extragalactic radio sources, and even SETI. At a forum organized by the International Space Science Institute-Beijing (ISSI-BJ), we discussed the prospect of opening up this window for astronomical observations by using a constellation of small or micro-satellites. We discussed the past experiments and the current ones such as the low frequency payload on Chang'e-4 mission lander, relay satellite and the Longjiang satellite, and also the future DSL mission, which is a linear array on lunar orbit which can make synthesized map of the whole sky as well as measure the global spectrum. We also discuss the synergy with other experiments, including ground global experiments such as EDGES, SARAS, SCI-HI and High-z, PRIZM/Albatros, ground imaging facillities such as LOFAR and MWA, and space experiments such as SUNRISE, DARE/DAPPER and PRATUSH. We also discussed some technical aspects of the DSL concept.

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.

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

We make an analytical estimate of the maximum 21 cm absorption signal from the cosmic dawn, taking into account the inhomogeneity of gas distribution in the intergalactic medium (IGM) due to non-linear structure formation. The gas located near halos is over-dense but adiabatically heated, while the gas far from halos is under-dense and hence cooler. The cumulative effect of adiabatic heating and cooling from this gas inhomogeneity results in a reduction in the maximum global 21 cm absorption depth by about 40% as compared with homogeneous IGM model, assuming saturated coupling between the spin temperature of neutral hydrogen (HI) and the adiabatic gas kinetic temperature.

We consider the reionization process in a cosmological model in which dark matter interacts with dark energy. Using a semi-analytical reionization model, we compute the evolution of the ionized fraction in terms of its spatial average and linear perturbations. We show that certain types of interactions between dark matter and dark energy can significantly affect the reionization history. We calculate the 21 cm signals in the interaction models, and compare the results with the predictions of the ${\rm \Lambda CDM}$ model.

The recently announced result by EDGES points an unexpected excess in the 21 cm global brightness temperature from cosmic dawn at $z\sim 17$, potentially indicating new phenomena beyond the $\Lambda$CDM model. A generic cosmological model which allows conformal and disformal couplings between dark matter and dark energy is employed to investigate the impact on the 21cm absorption signal and understand the EDGES anomaly. After exploring a wide range of parameter space for couplings, we find that the coupling effects can lead to a moderate change in the Hubble parameter while a negligible change in the spin temperature in the early Universe. Consequently, the decrease of the Hubble parameter from the mixed conformal and disformal couplings can reproduce the 21cm absorption approximately in consistent with the EDGES result at $z=17.5$. However, there is still tension in corresponding parameter space between EDGES and other cosmological observations for this model.

The Chinese satellite Wukong, also known as the DArk Matter Particle Explorer (DAMPE), has released its observation data of the cosmic ray (CR) electrons and positrons. The data shows an excess in the energy spectrum up to TeV energy, and possibly a peak-like fine structure at $\sim 1.4 \TeV$. We investigate the scenario that the source of the excess comes from dark matter annihilation or decay. We find that the annihilation or decay of diffuse dark matter particles in the Galactic halo can give excellent ($W^+W^-$ channel) or at least good (double $\tau^+\tau^-$ channel) fits to the broad excess. However, the annihilation cross-section is $10^{-23}\cm^3s^{-1}$, larger than required for getting the correct relic abundance. We then study whether the narrow peak at $\sim 1.4\TeV$ could be explained by a nearby subhalo, which thanks to the smaller distance, could supply $e^+e^-$ within a narrow energy range. We find that in order to produce a peak width less than the energy bin width (0.2 TeV), the source must be located within $r\lsim 0.53~\kpc$. Our global fit models do not produce the peak-like feature, instead at 1.4 TeV the spectrum show either a slope or a cliff-like feature. However, if less than optimal fit is allowed, the peak-like feature could be generated. Furthermore, an excellent fit with peak could be obtained with model B if the background is rescaled. If the dark matter decay and annihilation rates are determined using the broad excess, the required subhalo mass $\sim10^{5}~M_\odot$ for decay model, or $\sim10^{4.5}\Msun$ for annihilation model and a shallower density profile slope $\alpha=1.2$, or $\sim10^{2.5}\Msun$ for the steep density profile $\alpha=1.7$. However, the probability for the existence of a such nearby subhalo as massive as given above is very low.

Ultra-faint galaxies are hosted by small dark matter halos with shallow gravitational potential wells, hence their star formation activity is more sensitive to feedback effects. The shape of the faint-end of the high-$z$ galaxy luminosity function (LF) contains important information on star formation and its interaction with the reionization process during the Epoch of Reionization (EoR). High-$z$ galaxies with $M_{\rm UV}\gtrsim-17$ have only recently become accessible thanks to the Frontier Fields (FFs) survey combining deep {\it HST} imaging and the gravitational lensing effect. In this paper we investigate the faint-end of the LF at redshift $>$5 using the data of FFs clusters Abell 2744 (A2744), MACSJ0416.1-2403 (M0416), MACSJ0717.5+3745 (M0717) and MACSJ1149.5+2223 (M1149). We analyze both an empirical and a physically-motivated LF model to obtain constraints on a possible turn-over of LF at faint magnitudes. In the empirical model the LF drops fast when the absolute UV magnitude $M_{\rm UV}$ is much larger than a turn-over absolute UV magnitude $M_{\rm UV}^{\rm T}$. We obtain $M_{\rm UV}^{\rm T}\gtrsim-14.6 $ (15.2) at 1 (2) $\sigma$ confidence level (C.L.) for $z\sim6$. In the physically-motivated analytical model, star formation in halos with circular velocity below $v_c^*$ is fully quenched if these halos are located in ionized regions. Using updated lensing models and new additional FFs data, we re-analyze previous constraints on $v_c^*$ and $f_{\rm esc}$ presented by Castellano et al. 2016a (C16a) using a smaller dataset. We obtain new constraints on $v_c^*\lesssim 59$ km s$^{-1}$ and $f_{\rm esc}\lesssim 56\%$ (both at 2$\sigma$ C.L.) and conclude that there is no turn-over detected so far from the analyzed FFs data. Forthcoming {\it JWST} observations will be key to tight these constraints further.

Radio astronomical observation below 30 MHz is hampered by the refraction and absorption of the ionosphere, and the radio frequency interference (RFI), so far high angular resolution sky intensity map is not available. An interferometer array on lunar orbit provides a perfect observatory in this frequency band: it is out of the ionosphere and the Moon helps to block the RFIs from the Earth. The satellites can make observations on the far side of the Moon and then send back the data on the near side part of the orbit. However, for such array the traditional imaging algorithm is not applicable: the field of view is very wide (almost whole sky), and for baselines distributed on a plane, there is a mirror symmetry between the two sides of the plane. A further complication is that for each baseline, the Moon blocks part of the sky, but as the satellites orbit the Moon, both the direction of the baseline and the blocked sky change, so even imaging algorithms which can deal with non-coplanar baseline may not work in this case. Here we present an imaging algorithm based on solving the linear mapping equations relating the sky intensity to the visibilities. We show that the mirror symmetry can be broken by the three dimensional baseline distribution generated naturally by the precession of the orbital plane of the satellites. The algorithm is applicable and good maps could be reconstructed, even though for each baseline the sky blocking by the Moon is different. We also investigate how the map-making is affected by inhomogeneous baseline distributions.

The large-scale structure of the ionization field during the epoch of reionization (EoR) can be modeled by the excursion set theory. While the growth of ionized regions during the early stage are described by the "bubble model", the shrinking process of neutral regions after the percolation of the ionized region calls for an "island model". An excursion set based analytical model and a semi-numerical code (islandFAST) have been developed. The ionizing background and the bubbles inside the islands are also included in the treatment. With two kinds of absorbers of ionizing photons, i.e. the large-scale under-dense neutral islands and the small-scale over-dense clumps, the ionizing background are self-consistently evolved in the model.

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.

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.

We present the algorithm and main results of our semi-numerical simulation, islandFAST, which is developed from the 21cmFAST and designed for the late stage of reionization. The islandFAST predicts the evolution and size distribution of the large scale under-dense neutral regions (neutral islands), and we find that the late Epoch of Reionization (EoR) proceeds very fast, showing a characteristic scale of the neutral islands at each redshift. Using islandFAST, we compare the impact of two types of absorption systems, i.e. the large scale under-dense neutral islands versus small scale over-dense absorbers, in regulating the reionization process. The neutral islands dominate the morphology of the ionization field, while the small scale absorbers dominate the mean free path of ionizing photons, and also delay and prolong the reionization process. With our semi-numerical simulation, the evolution of the ionizing background can be derived self-consistently given a model for the small absorbers. The hydrogen ionization rate of the ionizing background is reduced by an order of magnitude in the presence of dense absorbers.

Recent measurements of the Luminosity Function (LF) of galaxies in the Epoch of Reionization (EoR, $z\gsim6$) indicate a very steep increase of the number density of low-mass galaxies populating the LF faint-end. However, as star formation in low-mass halos can be easily depressed or even quenched by ionizing radiation, a turnover is expected at some faint UV magnitudes. Using a physically-motivated analytical model, we quantify reionization feedback effects on the LF faint-end shape. We find that if reionization feedback is neglected, the power-law Schechter parameterization characterizing the LF faint-end remains valid up to absolute UV magnitude $\sim -9$. If instead radiative feedback is strong enough that quenches star formation in halos with circular velocity smaller than 50 km s$^{-1}$, the LF starts to drop at absolute UV magnitude $\sim -15$, i.e. slightly below the detection limits of current (unlensed) surveys at $z\sim5$. The LFs may rise again at higher absolute UV magnitude, where, as a result of interplay between reionization process and galaxy formation, most of the galaxy light is from relic stars formed before the EoR. We suggest that the galaxy number counts data, particularly in lensed fields, can put strong constraints on reionization feedback. In models with stronger reionization feedback, stars in galaxies with absolute UV magnitude higher than $\sim -13$ and smaller than $\sim - 8$ are typically older. Hence, the stellar age - UV magnitude relation can be used as an alternative feedback probe.

Intensity mapping (IM) is sensitive to the cumulative line emission of galaxies. As such it represents a promising technique for statistical studies of galaxies fainter than the limiting magnitude of traditional galaxy surveys. The strong hydrogen Ly-alpha line is the primary target for such an experiment, as its intensity is linked to star formation activity and the physical state of the interstellar (ISM) and intergalactic (IGM) medium. However, to extract the meaningful information one has to solve the confusion problems caused by interloping lines from foreground galaxies. We discuss here the challenges for a Ly-alpha IM experiment targeting z > 4 sources. We find that the Ly-alpha power spectrum can be in principle easily (marginally) obtained with a 40 cm space telescope in a few days of observing time up to z < 8 (z = 10) assuming that the interloping lines (e.g. H-alpha, [O II], [O III] lines) can be efficiently removed. We show that interlopers can be removed by using an ancillary photometric galaxy survey with limiting AB mag 26 in the NIR bands (Y, J, H, or K). This would enable detection of the Ly-alpha signal from 5 < z < 9 faint sources. However, if a [C II] IM experiment is feasible, by cross-correlating the Ly-alpha with the [C II] signal the required depth of the galaxy survey can be decreased to AB mag 24. This would bring the detection at reach of future facilities working in close synergy.

We exploit a sample of ultra-faint high-redshift galaxies (demagnified HST $H_{160}$ magnitude $>30$) in the Frontier Fields clusters A2744 and M0416 to constrain a theoretical model for the UV luminosity function (LF) in the presence of photoionization feedback. The objects have been selected on the basis of accurate photometric redshifts computed from multi-band photometry including 7 HST bands and deep $K_s$ and IRAC observations. Magnification is computed on an object-by-object basis from all available lensing models of the two clusters. We take into account source detection completeness as a function of luminosity and size, magnification effects and systematics in the lens modeling of the clusters under investigation. We find that our sample of high-$z$ galaxies constrain the cut-off halo circular velocity below which star-formation is suppressed by photo-ionization feedback to $v_c^{\rm cut} < 50$ km s$^{-1}$. This circular velocity corresponds to a halo mass of $\approx5.6\times10^9~M_\odot$ and $\approx2.3\times10^9~M_\odot$ at $z=5$ and 10 respectively: higher mass halos can thus sustain continuous star formation activity without being quenched by external ionizing flux. More stringent constraints are prevented by the uncertainty in the modeling of the cluster lens, as embodied by systematic differences among the lens models available.

Emission from high-$z$ galaxies must unquestionably contribute to the near-infrared background (NIRB). However, this contribution has so far proven difficult to isolate even after subtracting the resolved galaxies to deep levels. Remaining NIRB fluctuations are dominated by unresolved low-$z$ galaxies on small angular scales, and by an unidentified component with unclear origin on large scales ($\approx 1000''$). In this paper, by analyzing mock maps generated from semi-numerical simulations and empirically determined $L_{\rm UV} - M_{\rm h}$ relations, we find that fluctuations associated with galaxies at $5 < z < 10$ amount to several percent of the unresolved NIRB flux fluctuations. We investigate the properties of this component for different survey areas and limiting magnitudes. In all cases, we show that this signal can be efficiently, and most easily at small angular scales, isolated by cross-correlating the source-subtracted NIRB with Lyman Break Galaxies (LBGs) detected in the same field by {\tt HST} surveys. This result provides a fresh insight into the properties of reionization sources.

The power spectrum of Near InfraRed Background (NIRB) fluctuations measured at 3.6 $\mu$m by {\tt Spitzer} shows a clustering excess over the known galaxies signal that has been interpreted in terms of early ($z\simgt 13$), accreting (direct collapse) black holes (DCBH) or low-$z$ intrahalo light (IHL). In addition, these fluctuations correlate with the cosmic X-ray background (CXB) measured at (0.5-2) keV, supporting the black hole explanation. This scenario has been questioned by the recent detection of a correlation between the two {\tt CIBER} 1.1/1.6 $\mu$m bands with the 3.6 $\mu$m {\tt Spitzer} one. This correlation is hardly explained by early DCBHs that, due to intergalactic absorption, cannot contribute to the shortest wavelength bands. Here we show that the new correlation is caused instead by a Diffuse Galactic Light (DGL) component arising from Galactic stellar light scattered by dust. The black hole interpretation of the excess remains perfectly valid and, actually, the inclusion of DGL allows less demanding (by up to about 30\%) requirements on the DCBH abundance/mass.

After two ALMA observing cycles, only a handful of [CII] $158\,\mu m$ emission line searches in z>6 galaxies have reported a positive detection, questioning the applicability of the local [CII]-SFR relation to high-z systems. To investigate this issue we use the Vallini et al. 2013 (V13) model, based on high-resolution, radiative transfer cosmological simulations to predict the [CII] emission from the interstellar medium of a z~7 (halo mass $M_h=1.17\times10^{11}M_{\odot}$) galaxy. We improve the V13 model by including (a) a physically-motivated metallicity (Z) distribution of the gas, (b) the contribution of Photo-Dissociation Regions (PDRs), (c) the effects of Cosmic Microwave Background on the [CII] line luminosity. We study the relative contribution of diffuse neutral gas to the total [CII] emission ($F _{diff}/F_{tot}$) for different SFR and Z values. We find that the [CII] emission arises predominantly from PDRs: regardless of the galaxy properties, $F _{diff}/F_{tot}\leq 10$% since, at these early epochs, the CMB temperature approaches the spin temperature of the [CII] transition in the cold neutral medium ($T_{CMB}\sim T_s^{CNM}\sim 20$ K). Our model predicts a high-z [CII]-SFR relation consistent with observations of local dwarf galaxies ($0.02<Z/Z_{\odot}<0.5$). The [CII] deficit suggested by actual data ($L_{CII}<2.0\times 10^7 L_{\odot}$ in BDF3299 at z~7.1) if confirmed by deeper ALMA observations, can be ascribed to negative stellar feedback disrupting molecular clouds around star formation sites. The deviation from the local [CII]-SFR would then imply a modified Kennicutt-Schmidt relation in z>6 galaxies. Alternatively/in addition, the deficit might be explained by low gas metallicities ($Z<0.1 Z_{\odot}$).

Cosmic metal enrichment is one of the key physical processes regulating galaxy formation and the evolution of the intergalactic medium (IGM). However, determining the metal content of the most distant galaxies has proven so far almost impossible; also, absorption line experiments at $z\sim6$ become increasingly difficult because of instrumental limitations and the paucity of background quasars. With the advent of ALMA, far-infrared emission lines provide a novel tool to study early metal enrichment. Among these, the [CII] line at 157.74 $\mu$m is the most luminous line emitted by the interstellar medium of galaxies. It can also resonant scatter CMB photons inducing characteristic intensity fluctuations ($\Delta I/I_{CMB}$) near the peak of the CMB spectrum, thus allowing to probe the low-density IGM. We compute both [CII] galaxy emission and metal-induced CMB fluctuations at $z\sim 6$ by using Adaptive Mesh Refinement cosmological hydrodynamical simulations and produce mock observations to be directly compared with ALMA BAND6 data ($\nu_{obs}\sim 272$ GHz). The [CII] line flux is correlated with $M_{UV}$ as $\log(F_{peak}/\mu{\rm Jy})=-27.205-2.253\,M_{UV}-0.038\,M_{UV}^2$. Such relation is in very good agreement with recent ALMA observations (e.g. Maiolino et al. 2015; Capak et al. 2015) of $M_{UV}<-20$ galaxies. We predict that a $M_{UV}=-19$ ($M_{UV}=-18$) galaxy can be detected at $4\sigma$ in $\simeq40$ (2000) hours, respectively. CMB resonant scattering can produce $\simeq\pm 0.1\,\mu$Jy/beam emission/absorptions features that are very challenging to be detected with current facilities. The best strategy to detect these signals consists in the stacking of deep ALMA observations pointing fields with known $M_{UV}\simeq-19$ galaxies. This would allow to simultaneously detect both [CII] emission from galactic reionization sources and CMB fluctuations produced by $z\sim6$ metals.

The intensity mapping of the [CII] 157.7 $\rm \mu$m fine-structure emission line represents an ideal experiment to probe star formation activity in galaxies, especially in those that are too faint to be individually detected. Here, we investigate the feasibility of such an experiment for $z > 5$ galaxies. We construct the $L_{\rm CII} - M_{\rm h}$ relation from observations and simulations, then generate mock [CII] intensity maps by applying this relation to halo catalogs built from large scale N-body simulations. Maps of the extragalactic far-infrared (FIR) continuum, referred to as "foreground", and CO rotational transition lines and [CI] fine-structure lines referred to as "contamination", are produced as well. We find that, at 316 GHz (corresponding to $z_{\rm CII} = 5$), the mean intensities of the extragalactic FIR continuum, [CII] signal, all CO lines from $J=1$ to 13 and two [CI] lines are $\sim 3\times10^5$ Jy sr$^{-1}$, $\sim 1200$ Jy sr$^{-1}$, $\sim 800$ Jy sr$^{-1}$ and $\sim 100$ Jy sr$^{-1}$, respectively. We discuss a method that allows us to subtract the FIR continuum foreground by removing a spectrally smooth component from each line of sight, and to suppress the CO/[CI] contamination by discarding pixels that are bright in contamination emission. The $z > 5$ [CII] signal comes mainly from halos in the mass range $10^{11-12} \,M_\odot$; as this mass range is narrow, intensity mapping is an ideal experiment to investigate these early galaxies. In principle such signal is accessible to a ground-based telescope with a 6 m aperture, 150 K system temperature, a $128\times128$ pixels FIR camera in 5000 hr total integration time, however it is difficult to perform such an experiment by using currently available telescopes.

We study the Initial Mass Function (IMF) and host halo properties of Intermediate Mass Black Holes (IMBH, 10^{4-6} Msun) formed inside metal-free, UV illuminated atomic cooling haloes (virial temperature T_vir > 10^4 K) either via the direct collapse of the gas or via an intermediate Super Massive Star (SMS) stage. We achieve this goal in three steps: (a) we derive the gas accretion rate for a proto-SMS to undergo General Relativity instability and produce a direct collapse black hole (DCBH) or to enter the ZAMS and later collapse into a IMBH; (b) we use merger-tree simulations to select atomic cooling halos in which either a DCBH or SMS can form and grow, accounting for metal enrichment and major mergers that halt the growth of the proto-SMS by gas fragmentation. We derive the properties of the host halos and the mass distribution of black holes at this stage, and dub it the "Birth Mass Function"; (c) we follow the further growth of the DCBH due to accretion of leftover gas in the parent halo and compute the final IMBH mass.We consider two extreme cases in which minihalos (T_vir < 10^4 K) can (fertile) or cannot (sterile) form stars and pollute their gas leading to a different IMBH IMF. In the (fiducial) fertile case the IMF is bimodal extending over a broad range of masses, M= (0.5-20)x10^5 Msun, and the DCBH accretion phase lasts from 10 to 100 Myr. If minihalos are sterile, the IMF spans the narrower mass range M= (1-2.8)x10^6 Msun, and the DCBH accretion phase is more extended (70-120 Myr). We conclude that a good seeding prescription is to populate halos (a) of mass 7.5 < log (M_h/Msun) < 8, (b) in the redshift range 8 < z < 17, (c) with IMBH in the mass range 4.75 < log (M_BH/Msun) < 6.25.

By combining cosmological simulations with Frontier Fields project lens models we find that, in the most optimistic case, galaxies as faint as $m \approx 33 - 34$ (AB magnitude at $1.6 \rm \mu m$) can be detected in the Frontier Fields. Such faint galaxies are hosted by dark matter halos of mass $\sim10^9 M_\odot$ and dominate the ionizing photon budget over currently observed bright galaxies, thus allowing for the first time the investigation of the dominant reionization sources. In addition, the observed number of these galaxies can be used to constrain the role of feedback in suppressing star formation in small halos: for example, if galaxy formation is suppressed in halos with circular velocity $v_c < 50$ km s$^{-1}$, galaxies fainter than $m=31$ should not be detected in the FFs.

It has been proposed that the first, intermediate-mass ($\approx 10^{5-6}~M_\odot$) black holes might form through direct collapse of unpolluted gas in atomic-cooling halos exposed to a strong Lyman-Werner (LW) or near-infrared (NIR) radiation. As these systems are expected to be Compton-thick, photons above 13.6 eV are largely absorbed and re-processed into lower energy bands. It follows that direct collapse black holes (DCBHs) are very bright in the LW/NIR bands, typically outshining small high-redshift galaxies by more than 10 times. Once the first DCBHs form, they then trigger a runaway process of further DCBH formation, producing a sudden rise in their cosmic mass density. The universe enters the "DCBH era" at $z \approx 20$ when a large fraction of atomic-cooling halos are experiencing DCBH formation. By combining the clustering properties of the radiation sources with Monte Carlo simulations we show that in this scenario the DCBH mass density rises from $\sim 5$~$M_\odot$ Mpc$^{-3}$ at $z\sim 30$ to the peak value $\sim5\times10^5 M_\odot$ Mpc$^{-3}$ at $z \sim 14$ in our fiducial model. However, the abundance of \textit{active} (accreting) DCBHs drops after $z \sim 14$, as gas in the potential formation sites (unpolluted halos with virial temperature slightly above $10^4$~K) is photoevaporated. This effect almost completely suppresses DCBH formation after $z\sim 13$. The DCBH formation era lasts only $\approx 150$ Myr, but it might crucially provide the seeds of the supermassive black holes (SMBHs) powering $z\sim6$ quasars.

In this paper we investigate the nature and distribution of large neutral regions during the late epoch of reionization. In the "bubble model" of reionization, the mass distribution of large ionized regions ("bubbles") during the early stage of reionization is obtained by using the excursion set model, where the ionization of a region corresponds to the first up-crossing of a barrier by random trajectories. We generalize this idea, and develop a method to predict the distribution of large scale neutral regions during the late stage of reionization, taking into account the ionizing background after the percolation of HII regions. The large scale neutral regions which we call "neutral islands" are not individual galaxies or minihalos, but larger regions where fewer galaxies formed and hence ionized later, and they are identified in the excursion set model with the first down-crossings of the island barrier. Assuming that the consumption rate of ionizing background photons is proportional to the surface area of the neutral islands, we obtained the size distribution of the neutral islands. We also take the "bubbles-in-island" effect into account by considering the conditional probability of up-crossing a bubble barrier after down-crossing the island barrier. We find that this effect is very important. An additional barrier is set to avoid islands being percolated through. We find that there is a characteristic scale for the neutral islands, while the small islands are rapidly swallowed up by the ionizing background, this characteristic scale does not change much as the reionization proceeds.

Angular fluctuations of the Near InfraRed Background (NIRB) intensity are observed up to scales $\simlt 1^{\ensuremath{^{\circ}}}$. Their interpretation is challenging as even after removing the contribution from detected sources, the residual signal is $>10$ times higher than expected from distant galaxies below the detection limit and first stars. We propose here a novel interpretation in which early, intermediate mass, accreting direct collapse black holes (DCBH), which are too faint to be detected individually in current surveys, could explain the observed fluctuations. We find that a population of highly obscured ($N_{\rm H}\simgt 10^{25} \rm cm^{-2}$) DCBHs formed in metal-free halos with virial temperature $10^4$ K at $z\simgt 12$, can explain the observed level $\approx 10^{-3}$ (nW m$^{-2}$ sr$^{-1})^2$ of the 3.6 and 4.5 $\mu$m fluctuations on scales $>100''$. The signal on smaller scales is instead produced by undetected galaxies at low and intermediate redshifts. Albeit Compton-thick, at scales $\theta> 100''$ DCBHs produce a CXB (0.5-2 keV)-NIRB ($4.5 \rm \mu m$) cross-correlation signal of $\simeq 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ nW m$^{-2}$ sr$^{-1}$ slightly dependent on the specific value of the absorbing gas column ($N_{\rm H} \approx 10^{25} \rm cm^{-2}$) adopted and in agreement with the recent measurements by \cite{2012arXiv1210.5302C}. At smaller scales the cross-correlation is dominated by the emission of high-mass X-ray binaries (HMXB) hosted by the same low-$z$, undetected galaxies accounting for small scale NIRB fluctuations. These results outline the great potential of the NIRB as a tool to investigate the nature of the first galaxies and black holes.

Several independent measurements have confirmed the existence of fluctuations ($\delta F_{\rm obs}\approx 0.1 \rm nW/m^{2}/sr$ at $3.6 \rm \mu m$) up to degree angular scales in the source-subtracted Near InfraRed Background (NIRB) whose origin is unknown. By combining high resolution cosmological N-body/hydrodynamical simulations with an analytical model, and by matching galaxy Luminosity Functions (LFs) and the constraints on reionization simultaneously, we predict the NIRB absolute flux and fluctuation amplitude produced by high-$z$ ($z > 5$) galaxies (some of which harboring Pop III stars, shown to provide a negligible contribution). This strategy also allows us to make an empirical determination of the evolution of ionizing photon escape fraction: we find $f_{\rm esc} = 1$ at $z \ge 11$, decreasing to $\approx 0.05$ at $z = 5$. In the wavelength range $1.0-4.5 \rm \mu m$, the predicted cumulative flux is $F =0.2-0.04 \rm nW/m^2/sr$. However, we find that the radiation from high-$z$ galaxies (including those undetected by current surveys) is insufficient to explain the amplitude of the observed fluctuations: at $l=2000$, the fluctuation level due to $z > 5$ galaxies is $\delta F = 0.01-0.002 \rm nW/m^2/sr$, with a relative wavelength-independent amplitude $\delta F/F = 4$%. The source of the missing power remains unknown. This might indicate that an unknown component/foreground, with a clustering signal very similar to that of high-$z$ galaxies, dominates the source-subtracted NIRB fluctuation signal.

In a semi-numerical model of reionization, the evolution of ionization fraction is simulated approximately by the ionizing photon to baryon ratio criterion. In this paper we incorporate a semi-analytical model of galaxy formation based on the Millennium II N-body simulation into the semi-numerical modeling of reionization. The semi-analytical model is used to predict the production of ionizing photons, then we use the semi-numerical method to model the reionization process. Such an approach allows more detailed modeling of the reionization, and also connects observations of galaxies at low and high redshifts to the reionization history. The galaxy formation model we use was designed to match the low-$z$ observations, and it also fits the high redshift luminosity function reasonably well, but its prediction on the star formation falls below the observed value, and we find that it also underpredicts the stellar ionizing photon production rate, hence the reionization can not be completed at $z \sim 6$ without taking into account some other potential sources of ionization photons. We also considered simple modifications of the model with more top heavy initial mass functions (IMF), with which the reionization can occur at earlier epochs. The incorporation of the semi-analytical model may also affect the topology of the HI regions during the EoR, and the neutral regions produced by our simulations with the semi-analytical model appeared less poriferous than the simple halo-based models.

Compared with the cold dark matter (CDM) model, in the warm dark matter (WDM) model formation of small scale structure is suppressed. It is often thought that this would delay the reionization of the intergalactic medium (IGM), as the star formation rate during the epoch of reionization (EOR) would be lowered. However, during the later stage of the EOR, a large portion of the ionizing photons are consumed by recombination inside the minihalos, where the gas has higher density and recombination rates than the gas in the IGM. The suppression of small scale structure would therefore reduce the recombination rate, and could potentially shorten the reionization process. This effect is investigated here by using the analytical "bubble model" of reionization. We find that in some cases, though the initiation of the EOR is delayed in the WDM model, its completion could be even earlier than the CDM case, but the effect is generally small. We obtain limits on the WDM particles mass for different reionization redshifts.

It has been proposed that during the formation of the first generation stars there might be a "dark star" phase in which the power of the star comes from dark matter annihilation. The adiabatic contraction process to form the dark star would result in a highly concentrated density profile of the host halo at the same time, which may give enhanced indirect detection signals of dark matter. In this work we investigate the extragalactic $\gamma$-ray background from dark matter annihilation with such a dark star formation scenario, and employ the isotropic $\gamma$-ray data from Fermi-LAT to constrain the model parameters of dark matter. The results suffer from large uncertainties of both the formation rate of the first generation stars and the subsequent evolution effects of the host halos of the dark stars. We find, in the most optimistic case for $\gamma$-ray production via dark matter annihilation, the expected extragalactic $\gamma$-ray flux will be enhanced by 1-2 orders of magnitude. In such a case, the annihilation cross section of the supersymmetric dark matter can be constrained to the thermal production level, and the leptonic dark matter model which is proposed to explain the positron/electron excesses can be well excluded. Conversely, if the positron/electron excesses are of a dark matter annihilation origin, then the early Universe environment is such that no dark star can form.

The cosmic electron and positron excesses have been explained as possible dark matter (DM) annihilation products. In this work we investigate the possible effects of such a DM annihilation scenario during the evolution history of the Universe. We first calculate the extragalactic $\gamma$-ray background (EGRB), which is produced through the final state radiation of DM annihilation to charged leptons and the inverse Compton scattering between electrons/positrons and the cosmic microwave background. The DM halo profile and the minimal halo mass, which are not yet well determined from the current N-body simulations, are constrained by the EGRB data from EGRET and Fermi telescopes. Then we discuss the impact of such leptonic DM models on cosmic evolution, such as the reionization and heating of intergalactic medium, neutral Hydrogen 21 cm signal and suppression of structure formation. We show that the impact on the Hydrogen 21 cm signal might show interesting signatures of DM annihilation, but the influence on star formation is not remarkable. Future observations of the 21 cm signals could be used to place new constraints on the properties of DM.

We study the statistical properties of the cosmological 21 cm signal from both the intergalactic medium (IGM) and minihalos, using a reionization simulation that includes a self--consistent treatment of minihalo photoevaporation. We consider two models for minihalo formation and three typical thermal states of the IGM -- heating purely by ionization, heating from both ionizing and $\rm Ly\alpha$ photons, and a maximal "strong heating" model. We find that the signal from the IGM is almost always dominant over that from minihalos. In our calculation, the differential brightness temperature, ${\delta}T_b,$ of minihalos is never larger than 2 mK. Although there are indeed some differences in the signals from the minihalos and from the IGM, even with the planned generation of radio telescopes it will be unfeasible to detect them. However, minihalos significantly affect the ionization state of the IGM and the corresponding 21 cm flux.