The nature of the time correlated noise component (the 1/f noise) of single dish radio telescopes is critical to the detectability of the HI signal in intensity mapping experiments. In this paper, we present the 1/f noise properties of the MeerKAT receiver system using South Celestial Pole (SCP) tracking data. We estimate both the temporal power spectrum density and the 2D power spectrum density for each of the antennas and polarizations. We apply Singular Value Decomposition (SVD) to the dataset and show that, by removing the strongest components, the 1/f noise can be drastically reduced, indicating that it is highly correlated in frequency. Without SVD mode subtraction, the knee frequency over a $20\,$MHz integration is higher than $0.1\,\rm Hz$; with just $2$~mode subtraction, the knee frequency is reduced to $\sim 3\times 10^{-3}\,{\rm Hz}$, indicating that the system induced 1/f-type variations are well under the thermal noise fluctuations over a few hundred seconds time scales. The 2D power spectrum shows that the 1/f-type variations are restricted to a small region in the time-frequency space, either with long wavelength correlations in frequency or in time. This gives a wide range of cosmological scales where the 21cm signal can be measured without further need to calibrate the gain time fluctuations. Finally, we demonstrate that a simple power spectrum parameterization is sufficient to describe the data and provide fitting parameters for both the 1D and 2D power spectrum.

The efficient classification of different types of supernova is one of the most important problems for observational cosmology. However, spectroscopic confirmation of most objects in upcoming photometric surveys, such as the The Rubin Observatory Legacy Survey of Space and Time (LSST), will be unfeasible. The development of automated classification processes based on photometry has thus become crucial. In this paper we investigate the performance of machine learning (ML) classification on the final cosmological constraints using simulated lightcurves from The Supernova Photometric Classification Challenge, released in 2010. We study the use of different feature sets for the lightcurves and many different ML pipelines based on either decision tree ensembles or automated search processes. To construct the final catalogs we propose a threshold selection method, by employing a Bias-Variance tradeoff. This is a very robust and efficient way to minimize the Mean Squared Error. With this method we were able to get very strong cosmological constraints, which allowed us to keep ~75% of the total information in the type Ia SNe when using the SALT2 feature set and ~33% for the other cases (based on either the Newling model or on standard wavelet decomposition).

We investigate the prospects of measuring the cosmic magnification effect by cross-correlating neutral hydrogen intensity mapping (HI IM) maps with background optical galaxies. We forecast the signal-to-noise ratio for HI IM data from SKA1-MID and HIRAX, combined with LSST photometric galaxy samples. We find that, thanks to their different resolutions, SKA1-MID and HIRAX are highly complementary in such an analysis. We predict that SKA1-MID can achieve a detection with a signal-to-noise ratio of $\sim 15$ on a multipole range of $\ell \lesssim 200$, while HIRAX can reach a signal-to-noise ratio of $\sim 50$ on $200 < \ell < 2000$. We conclude that measurements of the cosmic magnification signal will be possible on a wide redshift range with foreground HI intensity maps up to $z \lesssim 2$, while optimal results are obtained when $0.6 \lesssim z \lesssim 1.3$. Finally, we perform a signal to noise analysis that shows how these measurements can constrain the HI parameters across a wide redshift range.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

Upcoming 21cm surveys with the SKA1-LOW telescope will enable imaging of the neutral hydrogen distribution on cosmological scales in the early Universe. These surveys are expected to generate huge imaging datasets that will encode more information than the power spectrum. This provides an alternative unique way to constrain the reionization history, which might break the degeneracy in the power spectral analysis. Using Convolutional Neural Networks (CNN), we create a fast estimator of the neutral fraction from the 21cm maps that are produced by our large semi-numerical simulation. Our estimator is able to efficiently recover the neutral fraction ($x_{\rm HI}$) at several redshifts with a high accuracy of 99\% as quantified by the coefficient of determination $R^{2}$. Adding the instrumental effects from the SKA design slightly increases the loss function, but nevertheless we are still able to recover the neutral fraction with a similar high accuracy of 98\%, which is only 1 per cent less. While a weak dependence on redshift is observed, the accuracy increases rapidly with decreasing neutral fraction. This is due to the fact that the instrumental noise increases towards high redshift where the Universe is highly neutral. Our results show the promise of directly using 21cm-tomography to constrain the reionization history in a model independent way, complementing similar efforts, such as those of the optical depth measurements from the Cosmic Microwave Background (CMB) observations by {\it Planck}.

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from bright astrophysical foregrounds. If not properly calibrated, each antenna element's complex gain can impart spectral structure on otherwise spectrally-smooth foregrounds, complicating that statistical separation. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as redundant-baseline calibration resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions---both in data and in simulations---and present strategies for mitigating that structure.

We present the radio luminosity function (RLF) of optically-selected quasars below 1~mJy, constructed by applying a Bayesian-fitting stacking technique to objects well below the nominal radio flux-density limit. We test the technique using simulated data, confirming that we can reconstruct the RLF over three orders of magnitude below the typical $5\sigma$ detection threshold. We apply our method to 1.4-GHz flux-densities from the Faint Images of the Radio Sky at Twenty-cm survey (FIRST), extracted at the positions of optical quasars from the Sloan Digital Sky Survey (SDSS) over seven redshift bins up to $z=2.15$ {and measure the RLF down to two orders of magnitude below the FIRST detection threshold}. In the lowest redshift bin ($0.2<z<0.45$), we find that our measured RLF agrees well with deeper data from the literature. The RLF for the radio-loud quasars flattens below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 25.5$ and becomes steeper again below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 24.8$, where radio-quiet quasars start to emerge. The radio luminosity where radio-quiet quasars emerge coincides with the luminosity where star-forming galaxies are expected to start to dominate the radio source counts. This implies that there could be a significant contribution from star formation in the host galaxies, but additional data is required to investigate this further. The higher-redshift bins show a similar behaviour as for the lowest-$z$ bin, implying that the same physical process may be responsible.

We introduce the "Asymmetric Radiative Transfer In Shells Technique" (ARTIST), a new method for photon propagation on large scales that explicitly conserves photons, propagates photons at the speed of light, approximately accounts for photon directionality, and closely reproduces results of more detailed radiative transfer (RT) codes. Crucially, it is computationally fast enough to evolve the large cosmological volumes required to predict the 21cm power spectrum on scales that will be probed by future experiments targeting the Epoch of Reionisation (EoR). Most semi-numerical models aimed at predicting the EoR 21cm signal make use of an excursion set formalism (ESF) approach, which achieves computational viability by compromising on photon conservation, constraining ionised regions to be spherical by construction, and not accounting for light-travel time. By implementing our RT method within the semi-numerical code SimFast21, we show that ARTIST predicts a significantly different evolution for the EoR ionisation field compared to the code's native ESF. In particular, ARTIST predicts a more gradual evolution of the volume-averaged ionisation fraction, and up to an order-of-magnitude difference in the ionisation power, depending on the physical parameters assumed. Its application to large-scale EoR simulations will therefore allow more physically-motivated constraints to be obtained for key EoR parameters, such as the escape fraction.

We present a Bayesian Stacking technique to directly measure the HI mass function (HIMF) and its evolution with redshift using galaxies formally below the nominal detection threshold. We generate galaxy samples over several sky areas given an assumed HIMF described by a Schechter function and simulate the HI emission lines with different levels of background noise to test the technique. We use Multinest to constrain the parameters of the HIMF in a broad redshift bin, demonstrating that the HIMF can be accurately reconstructed, using the simulated spectral cube far below the HI mass limit determined by the $5\sigma$ flux-density limit, i.e. down to $M_{\rm HI} = 10^{7.5}$ M$_{\odot}$ over the redshift range $0 < z < 0.55$ for this particular simulation, with a noise level similar to that expected for the MIGHTEE survey. We also find that the constraints on the parameters of the Schechter function, $\phi_{\star}$, $M_\star$ and $\alpha$ can be reliably fit, becoming tighter as the background noise decreases as expected, although the constraints on the redshift evolution are not significantly affected. All the parameters become better constrained as the survey area increases. In summary, we provide an optimal method for estimating the HI mass at cosmological distances that allows us to constrain the HI mass function below the detection threshold in forthcoming HI surveys. This study is a first step towards the measurement of the HIMF at high ($z>0.1$) redshifts.

Accurate simulations of the 21-cm signal from the Epoch of Reionization require the generation of maps at specific frequencies given the values of astrophysical and cosmological parameters. The peculiar velocities of the hydrogen atoms producing the 21-cm radiation result in a shift in the observed frequency of the 21-cm radiation and changes the amplitude of the signal itself. This is not an effect we can remove but instead needs to be accurately modelled to ensure we infer the correct physical parameters from an observation. We describe the full calculation of the distortion of the 21-cm signal, and propose a new code that integrates the 21-cm intensity along the line of sight for each individual light cone pixel to fully capture the intensity contributions from multiple redshifts. This algorithm naturally deals with the typical divergences found in standard approaches, allowing for large optical depths and 21-cm absorption events at high redshifts. We find the new method results in up to a 60% decrease in power on the largest scales on the sky, and an increase of over 80% on the smallest scales on the sky. We find that the new implementation of the light cone results in a longer tail of bright temperatures in the brightness temperature distribution, as a result of the successful circumventing of a previous cap that had to be implemented to avoid a divergence in the brightness temperature. We conclude that this full treatment of the evolution of the light cone pixel can be an important effect.

Global (i.e. sky-averaged) $21$~cm signal experiments can measure the evolution of the universe from the Cosmic Dawn to the Epoch of Reionization. These measurements are challenged by the presence of bright foreground emission that can be separated from the cosmological signal if its spectrum is smooth. This assumption fails in the case of single polarization antennas as they measure linearly polarized foreground emission - which is inevitably Faraday rotated through the interstellar medium. We investigate the impact of Galactic polarized foregrounds on the extraction of the global 21~cm signal through realistic sky and dipole simulations both in a low frequency band from $50$ to $100$~MHz, where a 21~cm absorption profile is expected, and in a higher frequency band ($100-200$~MHz). We find that the presence of a polarized contaminant with complex frequency structure can bias the amplitude and the shape of the reconstructed signal parameters in both bands. We investigate if polarized foregrounds can explain the unexpected $21$~cm Cosmic Dawn signal recently reported by the EDGES collaboration. We find that unaccounted polarized foreground contamination can produce an enhanced and distorted $21$~cm absorption trough similar to the anomalous profile reported by Bowman et al. (2018), and whose amplitude is in mild tension with the assumed input Gaussian profile (at $\sim 1.5 \sigma$ level). Moreover, we note that, under the hypothesis of contamination from polarized foreground, the amplitude of the reconstructed EDGES signal can be overestimated by around $30\%$, mitigating the requirement for an explanation based on exotic physics.

After a decade of design and construction, South Africa's SKA-MID precursor MeerKAT has begun its science operations. To make full use of the widefield capability of the array, it is imperative that we have an accurate model of the primary beam of its antennas. We have used an available L-band full-polarization astro-holographic observation and electromagnetic simulation to create sparse representations of the beam using principal components and Zernike polynomials. The spectral behaviour of the spatial coefficients has been modelled using discrete cosine transform. We have provided the Zernike-based model over a diameter of 10 degrees in an associated software tool that can be useful for direction dependent calibration and imaging. The model is more accurate for the diagonal elements of the beam Jones matrix and at lower frequencies. As we get more accurate beam measurements and simulations in the future, especially for the cross-polarization patterns, our pipeline can be used to create more accurate sparse representations of MeerKAT beam.

The angular power spectrum is a gauge-independent observable that is in principle the natural tool for analysing galaxy number counts. In practice, the problem is that the computational requirements for next-generation spectroscopic surveys such as Euclid and the Square Kilometre Array are currently unfeasible. We propose a new method to save computational time for spectroscopic angular power spectra. This hybrid method is modelled on the Fourier power spectrum approach of treating relatively thick redshift bins (redshift width ~0.1) as separate surveys. In the hybrid method, each thick bin is further subdivided into thin bins (redshift width ~0.01); all the correlations within each thick bin are computed, while cross-bin correlations beyond the thick bins are neglected. Constraints on cosmological parameters from the hybrid method are comparable to those from the standard galaxy power spectrum analysis - but they have the advantage that cosmic evolution, wide-angle and lensing effects are naturally included, while no Alcock-Paczynski correction is needed. The hybrid method delivers much tighter constraints than a 2D tomographic approach that is typical for photometric surveys, which considers only thick bins and the correlations between them. Furthermore, for standard cosmological parameters our method is not biased by neglecting the effects of lensing on number counts, while the tomographic method is strongly biased.

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.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

The near-infrared background between 0.5 $\mu$m to 2 $\mu$m contains a wealth of information related to radiative processes in the universe. Infrared background anisotropies encode the redshift-weighted total emission over cosmic history, including any spatially diffuse and extended contributions. The anisotropy power spectrum is dominated by undetected galaxies at small angular scales and diffuse background of Galactic emission at large angular scales. In addition to these known sources, the infrared background also arises from intra-halo light (IHL) at $z < 3$ associated with tidally-stripped stars during galaxy mergers. Moreover, it contains information on the very first galaxies from the epoch of reionization (EoR). The EoR signal has a spectral energy distribution (SED) that goes to zero near optical wavelengths due to Lyman absorption, while other signals have spectra that vary smoothly with frequency. Due to differences in SEDs and spatial clustering, these components may be separated in a multi-wavelength-fluctuation experiment. To study the extent to which EoR fluctuations can be separated in the presence of IHL, extra-galactic and Galactic foregrounds, we develop a maximum likelihood technique that incorporates a full covariance matrix among all the frequencies at different angular scales. We apply this technique to simulated deep imaging data over a 2$\times$100 deg$^2$ sky area from 0.75 $\mu$m to 5 $\mu$m in 9 bands and find that such a "frequency tomography" can successfully reconstruct both the amplitude and spectral shape for representative EoR, IHL and the foreground signals.

We use full sky simulations, including the effects of foreground contamination and removal, to explore multi-tracer synergies between a SKA-like 21cm intensity mapping survey and a LSST-like photometric galaxy redshift survey. In particular we study ratios of auto and cross-correlations between the two tracers as estimators of the ratio of their biases, a quantity that should benefit considerably from the cosmic variance cancellation of the multi-tracer approach. We show how well we should be able to measure the bias ratio on very large scales (down to $\ell \sim 3$), which is crucial to measure primordial non-Gaussianity and general relativistic effects on large scale structure. We find that, in the absence of foregrounds but with realistic noise levels of such surveys, the multi-tracer estimators are able to improve on the sensitivity of a cosmic-variance contaminated measurement by a factor of $2-4$. When foregrounds are included, estimators using the 21cm auto-correlation become biased. However, we show that cross-correlation estimators are immune to this and do not incur in any significant penalty in terms of sensitivity from discarding the auto-correlation data. However, the loss of long-wavelength radial modes caused by foreground removal in combination with the low redshift resolution of photometric surveys, reduces the sensitivity of the multi-tracer estimator, albeit still better than the cosmic variance contaminated scenario even in the noise free case. Finally we explore different alternative avenues to avoid this problem.

We study synergies between HI 21cm and H$\alpha$ intensity map observations, focusing on SKA1-like and SPHEREx-like surveys. We forecast how well such a combination can measure features in the angular power spectrum on the largest scales, that arise from primordial non-Gaussianity and from general relativistic effects. For the first time we consider Doppler, Sachs-Wolfe and integrated SW effects separately. We confirm that the single-tracer surveys on their own cannot detect general relativistic effects and can constrain the non-Gaussianity parameter $f_{\rm NL}$ only slightly better than Planck. Using the multi-tracer technique, constraints on $f_{\rm NL}$ can be pushed down to $\sim1$. Amongst the general relativistic effects, the Doppler term is detectable with the multi-tracer. The Sachs-Wolfe terms and the integrated SW effect are still not detectable.

We investigate the efficiency of screening mechanisms in the hybrid metric-Palatini gravity. The value of the field is computed around spherical bodies embedded in a background of constant density. We find a thin shell condition for the field depending on the background field value. In order to quantify how the thin shell effect is relevant, we analyze how it behaves in the neighborhood of different astrophysical objects (planets, moons or stars). We find that the condition is very well satisfied except only for some peculiar objects. Furthermore we establish bounds on the model using data from solar system experiments such as the spectral deviation measured by the Cassini mission and the stability of the Earth-Moon system, which gives the best constraint to date on $f(R)$ theories. These bounds contribute to fix the range of viable hybrid gravity models.

According to the Cosmological Principle, the matter distribution on very large scales should have a kinematic dipole that is aligned with that of the CMB. We determine the dipole anisotropy in the number counts of two all-sky surveys of radio galaxies. For the first time, this analysis is presented for the TGSS survey, allowing us to check consistency of the radio dipole at low and high frequencies by comparing the results with the well-known NVSS survey. We match the flux thresholds of the catalogues, with flux limits chosen to minimise systematics, and adopt a strict masking scheme. We find dipole directions that are in good agreement with each other and with the CMB dipole. In order to compare the amplitude of the dipoles with theoretical predictions, we produce sets of lognormal realisations. Our realisations include the theoretical kinematic dipole, galaxy clustering, Poisson noise, simulated redshift distributions which fit the NVSS and TGSS source counts, and errors in flux calibration. The measured dipole for NVSS is $\sim\!2$ times larger than predicted by the mock data. For TGSS, the dipole is almost $\sim\! 5$ times larger than predicted, even after checking for completeness and taking account of errors in source fluxes and in flux calibration. Further work is required to understand the nature of the systematics that are the likely cause of the anomalously large TGSS dipole amplitude.

We analyze a class of generalized inflationary models proposed in Ref. [1], known as $\beta$-exponential inflation. We show that this kind of potential can arise in the context of brane cosmology, where the field describing the size of the extra-dimension is interpreted as the inflaton. We discuss the observational viability of this class of model in light of the latest Cosmic Microwave Background (CMB) data from the Planck Collaboration through a Bayesian analysis, and impose tight constraints on the model parameters. We find that the CMB data alone prefer \textit{weakly} the minimal standard model ($\Lambda$CDM) over the $\beta$-exponential inflation. However, when current local measurements of the Hubble parameter, $H_0$, are considered, the $\beta$-inflation model is {\it moderately} preferred over the $\Lambda$CDM cosmology, making the study of this class of inflationary models interesting in the context of the current $H_0$ tension.

The detection of the redshifted cosmological $21$~cm line signal requires the removal of the Galactic and extragalactic foreground emission, which is orders of magnitude brighter anywhere in the sky. Foreground cleaning methods currently used are efficient in removing spectrally smooth components. However, they struggle in the presence of not spectrally smooth contamination that is, therefore, potentially the most dangerous one. An example of this is the polarized synchrotron emission, which is Faraday rotated by the interstellar medium and leaks into total intensity due to instrumental imperfections. In this work we present new full-sky simulations of this polarized synchrotron emission in the $50-200$~MHz range, obtained from the observed properties of diffuse polarized emission at low frequencies. The simulated polarized maps are made publicly available, aiming to provide more realistic templates to simulate the effect of instrumental leakage and the effectiveness of foreground separation techniques.

Recent precision cosmological parameter constraints imply that the spatial curvature of the Universe is essentially dynamically negligible - but only if relatively strong assumptions are made about the equation of state of dark energy (DE). When these assumptions are relaxed, strong degeneracies arise that make it hard to disentangle DE and curvature, degrading the constraints. We show that forthcoming 21cm intensity mapping experiments such as HIRAX are ideally designed to carry out model-independent curvature measurements, as they can measure the clustering signal at high redshift with sufficient precision to break many of the degeneracies. We consider two different model-independent methods, based on `avoiding' the DE-dominated regime and non-parametric modelling of the DE equation of state respectively. Our forecasts show that HIRAX will be able to improve upon current model-independent constraints by around an order of magnitude, reaching percent-level accuracy even when an arbitrary DE equation of state is assumed. In the same model-independent analysis, the sample variance limit for a similar survey is another order of magnitude better.

Neutral hydrogen (HI) will soon be the dark matter tracer observed over the largest volumes of Universe thanks to the 21 cm intensity mapping technique. To unveil cosmological information it is indispensable to understand the HI distribution with respect to dark matter. Using a full one-loop derivation of the power spectrum of HI, we show that higher order corrections change the amplitude and shape of the power spectrum on typical cosmological (linear) scales. These effects go beyond the expected dark matter non-linear corrections and include non-linearities in the way the HI signal traces dark matter. We show that, on linear scales at z = 1, the HI bias drops by up to 15% in both real and redshift space, which results in underpredicting the mass of the halos in which HI lies. Non-linear corrections give rise to a significant scale dependence when redshift space distortions arise, in particular on the scale range of the baryonic acoustic oscillations (BAO). There is a factor of 5 difference between the linear and full HI power spectra over the full BAO scale range, which will modify the ratios between the peaks. This effect will also be seen in other types of survey and it will be essential to take it into account in future experiments in order to match the expectations of precision cosmology.

Recent results have suggested that active galactic nuclei (AGN) could provide enough photons to reionise the Universe. We assess the viability of this scenario using a semi-numerical framework for modeling reionisation, to which we add a quasar contribution by constructing a Quasar Halo Occupation Distribution (QHOD) based on Giallongo et al. observations. Assuming a constant QHOD, we find that an AGN-only model cannot simultaneously match observations of the optical depth $\tau_e$, neutral fraction, and ionising emissivity. Such a model predicts $\tau_e$ too low by $\sim 2\sigma$ relative to Planck constraints, and reionises the Universe at $z\lesssim 5$. Arbitrarily increasing the AGN emissivity to match these results yields a strong mismatch with the observed ionising emissivity at $z\sim 5$. If we instead assume a redshift-independent AGN luminosity function yielding an emissivity evolution like that assumed in Madau & Haardt model, then we can match $\tau_e$ albeit with late reionisation, however such evolution is inconsistent with observations at $z\sim 4-6$ and poorly motivated physically. These results arise because AGN are more biased towards massive halos than typical reionising galaxies, resulting in stronger clustering and later formation times. AGN-dominated models produce larger ionising bubbles that are reflected in $\sim\times 2$ more 21cm power on all scales. A model with equal parts galaxies and AGN contribution is still (barely) consistent with observations, but could be distinguished using next-generation 21cm experiments HERA and SKA-low. We conclude that, even with recent claims of more faint AGN than previously thought, AGN are highly unlikely to dominate the ionising photon budget for reionisation.

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 discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).

Multi-wavelength astronomical studies brings a wealth of science within reach. One way to achieve a cross-wavelength analysis is via `stacking', i.e. combining precise positional information from an image at one wavelength with data from one at another wavelength in order to extract source-flux distributions and other derived quantities. For the first time we extend stacking to include the effects of confusion. We develop our algorithm in a fully Bayesian framework and apply it to the Square Kilometre Array Design Study (SKADS) simulation in order to extract galaxy number counts. Previous studies have shown that recovered source counts are highly biased high when confusion is non-negligible. With this new method, source counts are returned correctly. We also describe a novel estimator for quantifying the impact of confusion on stacking analyses. This method is an essential step in exploiting scientific return for upcoming deep radio surveys, e.g. MIGHTEE on MeerKAT.

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

We propose a new method to test modified gravity theories, taking advantage of the available data on extrasolar planets. We computed the deviations from the Kepler third law and use that to constrain gravity theories beyond General Relativity. We investigate gravity models which incorporate three screening mechanisms: the Chameleon, the Symmetron and the Vainshtein. We find that data from exoplanets orbits are very sensitive to the screening mechanisms putting strong constraints in the parameter space for the Chameleon models and the Symmetron, complementary and competitive to other methods, like interferometers and solar system. With the constraints on Vainshtein we are able to work beyond the hypothesis that the crossover scale is of the same order of magnitude than the Hubble radius $r_c \sim H_0^{-1}$, which makes the screening work automatically, testing how strong this hypothesis is and the viability of other scales.

The recent low value of Planck (2016) integrated optical depth to Thomson scattering suggests that the reionization occurred fairly suddenly, disfavoring extended reionization scenarios. This will have a significant impact on the 21cm power spectrum. Using a semi-numerical framework, we improve our model from Hassan et al. (2016) to include time-integrated ionisation and recombination effects, and find that this leads to more sudden reionisation. It also yields larger HII bubbles which leads to an order of magnitude more 21cm power on large scales, while suppressing the small scale ionization power. Local fluctuations in the neutral hydrogen density play the dominant role in boosting the 21cm power spectrum on large scales, while recombinations are subdominant. We use a Monte Carlo Markov Chain approach to constrain our model to observations of the star formation rate functions at z = 6,7,8 from Bouwens et al. (2015), the Planck (2016) optical depth measurements, and the Becker & Bolton (2013) ionising emissivity data at z~5. We then use this constrained model to perform 21cm forecasting for LOFAR, HERA, and SKA in order to determine how well such data can characterise the sources driving reionisation. We find that the 21cm power spectrum alone can somewhat constrain the halo mass dependence of ionising sources, the photon escape fraction and ionising amplitude, but combining the 21cm data with other current observations enables us to separately constrain all these parameters. Our framework illustrates how 21cm data can play a key role in understanding the sources and topology of reionisation as observations improve.

The objective of this paper is to present the main characteristics of a wide-field MFAA precursor that we envisage to be built at the SKA site in South Africa. Known as MANTIS (the Mid-Frequency Aperture Array Transient and Intensity-Mapping System), this ambitious instrument will represent the next logical step towards the MFAA based SKA telescope. The goal is to use innovative aperture array technology at cm wavelengths, in order to demonstrate the feasibility of deploying huge collecting areas at modest construction and operational cost. Such a transformative step is required in order to continue the exponential progress in radio telescope performance, and to make the ambitious scale of the SKA Phase 2 a realistic near-time proposition.

Intensity mapping is now becoming a useful tool to study the large-scale structure of the universe through spatial variations in the integrated emission from galaxies and the intergalactic medium. We study intensity mapping of the H-alpha 6563AA, [OIII]5007AA, [OII]3727AA and H-beta 4861AA lines at 0.8<z<5.2. The mean intensities of these four emission lines are estimated using the observed luminosity functions (LFs), cosmological simulations, and the star formation rate density (SFRD) derived from observations at z<5. We calculate the intensity power spectra and consider the foreground contamination of other lines at lower redshifts. We use the proposed NASA small explorer SPHEREx (the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) as a case study for the detectability of the intensity power spectra of the four emission lines. We also investigate the cross correlation with the 21-cm line probed by CHIME (the Canadian Hydrogen Intensity Mapping Experiment), Tianlai experiment and SKA (the Square Kilometer Array) at 0.8<z<2.4. We find both the auto and cross power spectra can be well measured for the H-alpha, [OIII] and [OII] lines at z<3, while it is more challenging for the H-beta line. Finally, we estimate the constraint on the SFRD from intensity mapping, and find we can reach accuracy higher than 7% at z<4, which is better than usual measurements using the LFs of galaxies.

We study the dispersion relation for parallel propagating ion-cyclotron (IC) waves in a dusty plasma, considering situations where the velocity dispersion along perpendicular direction is greater than along the parallel direction, and considering the use of product-bi-kappa (PBK) velocity distributions for the plasma particles. The results obtained by numerical solution of the dispersion relation, in a case with isotropic Maxwellian distributions for electrons and PBK distribution for ions, show the occurrence of the electromagnetic ion-cyclotron instability (EMIC), and show that the decrease in the kappa indexes of the PBK ion distribution leads to significant increase in the magnitude of the growth rates and in the range of wavenumber for which the instability occurs. On the other hand, for anisotropic Maxwellian distribution for ions and PBK distribution for electrons, the decrease of the kappa index in the PBK electron distribution contributes to reduce the growth rate of the EMIC instability, but the reduction effect is less pronounced than the increase obtained with ion PBK distribution with the same $\kappa$ index. The results obtained also show that, as a general rule, the presence of a dust population contributes to reduce the instability in magnitude of the growth rates and range, but that in the case of PBK ion distribution with small kappa indexes the instability may continue to occur for dust populations which would eliminate completely the instability in the case of bi-Maxwellian ion distributions. It has also been seen that the anisotropy due to the kappa indexes in the ion PBK distribution is not so efficient in producing the EMIC instability as the ratio of perpendicular and parallel ion temperatures, for equivalent value of the effective temperature.

It is usually assumed that we will need to wait until next-generation surveys like Euclid, LSST and SKA, in order to improve on the current best constraints on primordial non-Gaussianity from the Planck experiment. We show that two contemporary surveys, with the SKA precursor MeerKAT and the Dark Energy Survey (DES), can be combined using the multi-tracer technique to deliver an accuracy on measurement of $f_{\rm NL}$ that is up to three times better than Planck.

We present a systematic study of the intensity mapping technique using updated models for the different emission lines from galaxies and identify which ones are more promising for cosmological studies of the post reionization epoch. We consider the emission of ${\rm Ly\alpha}$, ${\rm H\alpha}$, H$\beta$, optical and infrared oxygen lines, nitrogen lines, CII and the CO rotational lines. We then identify that ${\rm Ly\alpha}$, ${\rm H\alpha}$, OII, CII and the lowest rotational CO lines are the best candidates to be used as IM probes. These lines form a complementary set of probes of the galaxies emission spectra. We then use reasonable experimental setups from current, planned or proposed experiments to access the detectability of the power spectrum of each emission line. Intensity mapping of ${\rm Ly\alpha}$ emission from $z=2$ to 3 will be possible in the near future with HETDEX, while far-infrared lines require new dedicated experiments. We also show that the proposed SPHEREx satellite can use OII and ${\rm H\alpha}$ IM to study the large-scale distribution of matter in intermediate redshifts of 1 to 4. We found that submilimeter experiments with bolometers can have similar performances at intermediate redshifts using CII and CO(3-2).

The Square Kilometre Array (SKA) will produce spectroscopic surveys of tens to hundreds of millions of HI galaxies, eventually covering 30,000 sq. deg. and reaching out to redshift z~2. The huge volumes probed by the SKA will allow for some of the best constraints on primordial non-Gaussianity, based on measurements of the large-scale power spectrum. We investigate various observational set-ups for HI galaxy redshift surveys, compatible with the SKA Phase 1 and Phase 2 (full SKA) configurations. We use the corresponding number counts and bias for each survey from realistic simulations and derive the magnification bias and the evolution of source counts directly from these. For the first time, we produce forecasts that fully include the general relativistic effects on the galaxy number counts. These corrections to the standard analysis become important on very large scales, where the signal of primordial non-Gaussianity grows strongest. Our results show that, for the full survey, the non-Gaussianity parameter fNL can be constrained down to an accuracy of 1.54. This improves the current limit set by the Planck satellite by a factor of five, using a completely different approach.

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a new 400-800MHz radio interferometer under development for deployment in South Africa. HIRAX will comprise 1024 six meter parabolic dishes on a compact grid and will map most of the southern sky over the course of four years. HIRAX has two primary science goals: to constrain Dark Energy and measure structure at high redshift, and to study radio transients and pulsars. HIRAX will observe unresolved sources of neutral hydrogen via their redshifted 21-cm emission line (`hydrogen intensity mapping'). The resulting maps of large-scale structure at redshifts 0.8-2.5 will be used to measure Baryon Acoustic Oscillations (BAO). HIRAX will improve upon current BAO measurements from galaxy surveys by observing a larger cosmological volume (larger in both survey area and redshift range) and by measuring BAO at higher redshift when the expansion of the universe transitioned to Dark Energy domination. HIRAX will complement CHIME, a hydrogen intensity mapping experiment in the Northern Hemisphere, by completing the sky coverage in the same redshift range. HIRAX's location in the Southern Hemisphere also allows a variety of cross-correlation measurements with large-scale structure surveys at many wavelengths. Daily maps of a few thousand square degrees of the Southern Hemisphere, encompassing much of the Milky Way galaxy, will also open new opportunities for discovering and monitoring radio transients. The HIRAX correlator will have the ability to rapidly and eXperimentciently detect transient events. This new data will shed light on the poorly understood nature of fast radio bursts (FRBs), enable pulsar monitoring to enhance long-wavelength gravitational wave searches, and provide a rich data set for new radio transient phenomena searches. This paper discusses the HIRAX instrument, science goals, and current status.

It is well known that the probability distribution function (PDF) of galaxy density contrast is approximately lognormal; whether the PDF of mass fluctuations derived from weak lensing convergence (kappa_WL) is lognormal is less well established. We derive PDFs of the galaxy and projected matter density distributions via the Counts in Cells (CiC) method. We use maps of galaxies and weak lensing convergence produced from the Dark Energy Survey (DES) Science Verification data over 139 deg^2. We test whether the underlying density contrast is well described by a lognormal distribution for the galaxies, the convergence and their joint PDF. We confirm that the galaxy density contrast distribution is well modeled by a lognormal PDF convolved with Poisson noise at angular scales from 10-40 arcmin (corresponding to physical scales of 3-10 Mpc). We note that as kappa_WL is a weighted sum of the mass fluctuations along the line of sight, its PDF is expected to be only approximately lognormal. We find that the kappa_WL distribution is well modeled by a lognormal PDF convolved with Gaussian shape noise at scales between 10 and 20 arcmin, with a best-fit chi^2/DOF of 1.11 compared to 1.84 for a Gaussian model, corresponding to p-values 0.35 and 0.07 respectively, at a scale of 10 arcmin. Above 20 arcmin a simple Gaussian model is sufficient. The joint PDF is also reasonably fitted by a bivariate lognormal. As a consistency check we compare the variances derived from the lognormal modelling with those directly measured via CiC. Our methods are validated against maps from the MICE Grand Challenge N-body simulation.

Intensity mapping of the neutral hydrogen brightness temperature promises to provide a three-dimensional view of the universe on very large scales. Nonlinear effects are typically thought to alter only the small-scale power, but we show how they may bias the extraction of cosmological information contained in the power spectrum on ultra-large scales. For linear perturbations to remain valid on large scales, we need to renormalize perturbations at higher order. In the case of intensity mapping, the second-order contribution to clustering from weak lensing dominates the nonlinear contribution at high redshift. Renormalization modifies the mean brightness temperature and therefore the evolution bias. It also introduces a term that mimics white noise. These effects may influence forecasting analysis on ultra-large scales.

We have investigated structure formation in the $\gamma$ gravity $f(R)$ model with {\it N}-body simulations. The $\gamma$ gravity model is a proposal which, unlike other viable $f(R)$ models, not only changes the gravitational dynamics, but can in principle also have signatures at the background level that are different from those obtained in $\Lambda$CDM (Cosmological constant, Cold Dark Matter). The aim of this paper is to study the nonlinear regime of the model in the case where, at late times, the background differs from $\Lambda$CDM. We quantify the signatures produced on the power spectrum, the halo mass function, and the density and velocity profiles. To appreciate the features of the model, we have compared it to $\Lambda$CDM and the Hu-Sawicki $f(R)$ models. For the considered set of parameters we find that the screening mechanism is ineffective, which gives rise to deviations in the halo mass function that disagree with observations. This does not rule out the model per se, but requires choices of parameters such that $|f_{R0}|$ is much smaller, which would imply that its cosmic expansion history cannot be distinguished from $\Lambda$CDM at the background level.

We use type Ia supernovae (SN Ia) data in combination with recent baryonic acoustic oscillations (BAO) and cosmic microwave background (CMB) observations to constrain a kink-like parametrization of the deceleration parameter ($q$). This $q$-parametrization can be written in terms of the initial ($q_i$) and present ($q_0$) values of the deceleration parameter, the redshift of the cosmic transition from deceleration to acceleration ($z_t$) and the redshift width of such transition ($\tau$). By assuming a flat space geometry, $q_i=1/2$ and adopting a likelihood approach to deal with the SN Ia data we obtain, at the 68% confidence level (C.L.), that: $z_t=0.56^{+0.13}_{-0.10}$, $\tau=0.47^{+0.16}_{-0.20}$ and $q_0=-0.31^{+0.11}_{-0.11}$ when we combine BAO/CMB observations with SN Ia data processed with the MLCS2k2 light-curve fitter. When in this combination we use the SALT2 fitter we get instead, at the same C.L.: $z_t=0.64^{+0.13}_{-0.07}$, $\tau=0.36^{+0.11}_{-0.17}$ and $q_0=-0.53^{+0.17}_{-0.13}$. Our results indicate, with a quite general and model independent approach, that MLCS2k2 favors Dvali-Gabadadze-Porrati-like cosmological models, while SALT2 favors $\Lambda$CDM-like ones. Progress in determining the transition redshift and/or the present value of the deceleration parameter depends crucially on solving the issue of the difference obtained when using these two light-curve fitters.

We revisit the kink-like parametrization of the deceleration parameter q(z), which considers a transition, at $z_t$, from cosmic deceleration to acceleration. In this parametrization the initial, at $z \gg z_t$, value of the q-parameter is $q_i$, its final, z=-1, value is $q_f$ and the duration of the transition is parametrized by $\tau$. By assuming a flat space geometry we obtain constraints on the free parameters of the model using recent data from type Ia supernovae, BAO, CMB and the Hubble parameter H(z). The use of H(z) data introduces an explicit dependence of the combined likelihood on the present value of the Hubble parameter $H_0$, allowing us to explore the influence of different priors when marginalizing over this parameter. We also study the importance of the CMB information in the results by considering data from WMAP7, WMAP9 and Planck. We show that the contours and best fit do not depend much on the different CMB data used and that the considered new BAO data is responsible for most of the improvement in the results. Assuming a flat space geometry, $q_i=1/2$ and expressing the present value of the deceleration parameter $q_0$ as a function of the other three free parameters, we obtain $z_t=0.67^{+0.10}_{-0.08}$, $\tau=0.26^{+0.14}_{-0.10}$ and $q_0=-0.48^{+0.11}_{-0.13}$, at 68\% of confidence level, with an uniform prior over $H_0$. If in addition we fix $q_f=-1$, as in flat $\Lambda$CDM, DGP and Chaplygin quartessence that are special models described by our parametrization, we get $z_t=0.66^{+0.03}_{-0.04}$, $\tau=0.33^{+0.04}_{-0.04}$ and $q_0=-0.54^{+0.05}_{-0.07}$, in excellent agreement with flat $\Lambda$CDM for which $\tau=1/3$. We also obtain for flat wCDM, another dark energy model described by our parametrization, the constraint on the equation of state parameter -1.22 < w < -0.78 at more than 99% confidence level.

Cosmic Dawn Intensity Mapper is a "Probe Class" mission concept for reionization studies of the universe. It will be capable of spectroscopic imaging observations between 0.7 to 6-7 microns in the near-Infrared. The primary observational objective is pioneering observations of spectral emission lines of interest throughout the cosmic history, but especially from the first generation of distant, faint galaxies when the universe was less than 800 million years old. With spectro-imaging capabilities, using a set of linear variable filters (LVFs), CDIM will produce a three-dimensional tomographic view of the epoch of reionization (EoR). CDIM will also study galaxy formation over more than 90% of the cosmic history and will move the astronomical community from broad-band astronomical imaging to low-resolution (R=200-300) spectro-imaging of the universe.

We explore the impact of incorporating physically motivated ionisation and recombination rates on the history and topology of cosmic reionisation, by incorporating inputs from small-volume hydrodynamic simulations into a semi-numerical code, SimFast21, that evolves reionisation on large scales. We employ radiative hydrodynamic simulations to parameterize the ionisation rate Rion and recombination rate Rrec as functions of halo mass, overdensity and redshift. We find that Rion is super-linearly dependent on halo mass (Rion ~ Mh^1.41), in contrast to previous assumptions. We implement these scalings into SimFast21 to identify the ionized regions. We tune our models to be consistent with recent observations of the optical depth, ionizing emissivity, and neutral fraction by the end of reionisation. We require an average photon escape fraction fesc=0.04 within ~ 0.5 cMpc cells, independent of halo mass or redshift, to simultaneously match these data. We present predictions for the 21cm power spectrum, and show that it is converged with respect to simulation volume. We find that introducing superlinearly mass-dependent ionisations increases the duration of reionisation and boosts the small-scale 21cm power by ~ 2-3 at intermediate phases of reionisation. Introducing inhomogeneous recombinations reduces ionised bubble sizes and suppresses large-scale 21cm power by ~ 2-3. Moreover, gas clumping on sub-cell scales has a minimal effect on the 21cm power, indicating that robust predictions do not depend on the behaviour of kpc-scale structures. The superlinear ionisations significantly increase the median halo mass scale for ionising photon output to >10^10 Mo, giving greater hope for detecting most of ionising sources with next-generation facilities. These results highlight the importance of more accurately treating ionising sources and recombinations for modeling reionisation and its 21cm signal.

(Abridged) We investigate the intrinsic shapes of low-luminosity galaxies in the central 300 kpc of the Virgo cluster using deep imaging obtained as part of the NGVS. We build a sample of nearly 300 red-sequence cluster members in the yet unexplored $-14 < M_{g} < -8$ magnitude range. The observed distribution of apparent axis ratios is then fit by families of triaxial models with normally-distributed intrinsic ellipticities and triaxialities. We develop a Bayesian framework to explore the posterior distribution of the model parameters, which allows us to work directly on discrete data, and to account for individual, surface brightness-dependent axis ratio uncertainties. For this population we infer a mean intrinsic ellipticity E=0.43, and a mean triaxiality T=0.16. This implies that faint Virgo galaxies are best described as a family of thick, nearly oblate spheroids with mean intrinsic axis ratios 1:0.94:0.57. We additionally attempt a study of the intrinsic shapes of Local Group satellites of similar luminosities. For the LG population we infer a slightly larger mean intrinsic ellipticity E=0.51, and the paucity of objects with round apparent shapes translates into more triaxial mean shapes, 1:0.76:0.49. We finally compare the intrinsic shapes of NGVS low-mass galaxies with samples of more massive quiescent systems, and with field, star-forming galaxies of similar luminosities. We find that the intrinsic flattening in this low-luminosity regime is almost independent of the environment in which the galaxy resides--but there is a hint that objects may be slightly rounder in denser environments. The comparable flattening distributions of low-luminosity galaxies that have experienced very different degrees of environmental effects suggests that internal processes are the main drivers of galaxy structure at low masses--with external mechanisms playing a secondary role.

We investigate the cross-correlation between the spatial distribution of Lyman Break Galaxies (LBGs) and the 21cm intensity mapping signal at $z\sim[3-5]$. At these redshifts, galactic feedback is supposed to only marginally affect the matter power spectrum, and the neutral hydrogen distribution is independently constrained by quasar spectra. Using a high resolution N-body simulation, populated with neutral hydrogen a posteriori, we forecast for the expected LBG-21cm cross-spectrum and its error for a 21cm field observed by the Square Kilometre Array (SKA1-LOW and SKA1-MID), combined with a spectroscopic LBG survey with the same volume. The cross power can be detected with a signal-to-noise ratio (SNR) up to ~10 times higher (and down to ~4 times smaller scales) than the 21cm auto-spectrum for this set-up, with the SNR depending only very weakly on redshift and the LBG population. We also show that while both the 21cm auto- and LBG-21cm cross-spectra can be reliably recovered after the cleaning of smooth-spectrum foreground contamination, only the cross-power is robust to problematic non-smooth foregrounds like polarized synchrotron emission.

Future surveys of large-scale structure will be able to measure perturbations on the scale of the cosmological horizon, and so could potentially probe a number of novel relativistic effects that are negligibly small on sub-horizon scales. These effects leave distinctive signatures in the power spectra of clustering observables and, if measurable, would open a new window on relativistic cosmology. We quantify the size and detectability of the effects for the most relevant future large-scale structure experiments: spectroscopic and photometric galaxy redshift surveys, intensity mapping surveys of neutral hydrogen, and radio continuum surveys. Our forecasts show that next-generation experiments, reaching out to redshifts $z\simeq 4$, will not be able to detect previously-undetected general-relativistic effects by using individual tracers of the density field, although the contribution of weak lensing magnification on large scales should be clearly detectable. We also perform a rigorous joint forecast for the detection of primordial non-Gaussianity through the excess power it produces in the clustering of biased tracers on large scales, finding that uncertainties of $\sigma(f_{\rm NL})\sim 1-2$ should be achievable. We study the level of degeneracy of these large-scale effects with several tracer-dependent nuisance parameters, quantifying the minimal priors on the latter that are needed for an optimal measurement of the former. Finally, we discuss the systematic effects that must be mitigated to achieve this level of sensitivity, and some alternative approaches that should help to improve the constraints. The computational tools developed to carry out this study, which requires the full-sky computation of the theoretical angular power spectra for $\mathcal{O}(100)$ redshift bins, as well as realistic models of the luminosity function, are publicly available.

In this paper, we use a semi-analytical approach to analyze the global structure of the phase space of the planar planetary 3/1 mean-motion resonance, in cases where the outer planet is more massive than its inner companion. We show that the resonant dynamics can be described using only two fundamental parameters, the total angular momentum and the scaling parameter. The topology of the Hamiltonian function describing the resonant behaviour is studied on the representative planes that allows us to investigate a large domain of the phase space of the three-body problem without time-expensive numerical integrations of the equations of motion, and without any restriction on the magnitude of the planetary eccentricities. The families of the well known Apsidal Corotation Resonances (ACR) parameterized by the planetary mass ratio are obtained and their stability is analyzed. The main dynamical features in the domains around ACR are also investigated in detail by means of spectral analysis techniques, which allow us to detect the regions of different regimes of motion of resonant systems. The construction of dynamical maps for various values of the total angular momentum shows the evolution of domains of stable motion with the eccentricities, identifying possible configurations suitable for exoplanetary systems.

We present the results of the first test plates of the extended Baryon Oscillation Spectroscopic Survey. This paper focuses on the emission line galaxies (ELG) population targetted from the Dark Energy Survey (DES) photometry. We analyse the success rate, efficiency, redshift distribution, and clustering properties of the targets. From the 9000 spectroscopic redshifts targetted, 4600 have been selected from the DES photometry. The total success rate for redshifts between 0.6 and 1.2 is 71\% and 68\% respectively for a bright and faint, on average more distant, samples including redshifts measured from a single strong emission line. We find a mean redshift of 0.8 and 0.87, with 15 and 13\% of unknown redshifts respectively for the bright and faint samples. In the redshift range 0.6<z<1.2, for the most secure spectroscopic redshifts, the mean redshift for the bright and faint sample is 0.85 and 0.9 respectively. Star contamination is lower than 2\%. We measure a galaxy bias averaged on scales of 1 and 10~Mpc/h of 1.72 \pm 0.1 for the bright sample and of 1.78 \pm 0.12 for the faint sample. The error on the galaxy bias have been obtained propagating the errors in the correlation function to the fitted parameters. This redshift evolution for the galaxy bias is in agreement with theoretical expectations for a galaxy population with MB-5\log h < -21.0. We note that biasing is derived from the galaxy clustering relative to a model for the mass fluctuations. We investigate the quality of the DES photometric redshifts and find that the outlier fraction can be reduced using a comparison between template fitting and neural network, or using a random forest algorithm.

Next-generation cosmological surveys will probe ever larger volumes of the Universe, including the largest scales, near and beyond the horizon. On these scales, the galaxy power spectrum carries signatures of local primordial non-Gaussianity (PNG) and horizon-scale general relativistic (GR) effects. However, cosmic variance limits the detection of horizon-scale effects. Combining different surveys via the multi-tracer method allows us to reduce the effect down cosmic variance. This method benefits from large bias differences between two tracers of the underlying dark matter distribution, which suggests a multi-wavelength combination of large volume surveys that are planned on a similar timescale. We show that the combination of two contemporaneous surveys, a large neutral hydrogen intensity mapping survey in SKA Phase1 and a Euclid-like photometric survey, will provide unprecedented constraints on PNG as well as detection of the GR effects. We forecast that the error on local PNG will break through the cosmic variance limit on cosmic microwave background surveys and achieve $\sigma(f_{NL})\simeq1.4-0.5$, depending on assumed priors, bias, and sky coverage. GR effects are more robust to changes in the assumed fiducial model, and we forecast that they can be detected with a signal-to-noise of about $14$.

We report the detection of a rare transiting brown dwarf with a mass of 59 M_Jup and radius of 1.1 R_Jup around the metal-rich, [Fe/H] = +0.44, G9V star CoRoT-33. The orbit is eccentric (e = 0.07) with a period of 5.82 d. The companion, CoRoT-33b, is thus a new member in the so-called brown dwarf desert. The orbital period is within 3% to a 3:2 resonance with the rotational period of the star. CoRoT-33b may be an important test case for tidal evolution studies. The true frequency of brown dwarfs close to their host stars (P < 10 d) is estimated to be approximately 0.2% which is about six times smaller than the frequency of hot Jupiters in the same period range. We suspect that the frequency of brown dwarfs declines faster with decreasing period than that of giant planets.

This paper deals with the application of the creep tide theory (Ferraz-Mello, Cel. Mech. Dyn. Astron. vol. 116, 109, 2013) to the study of the rotation of stars hosting massive close-in planets. The stars have nearly the same tidal relaxation factors as gaseous planets and the evolution of their rotation is similar to that of close-in hot Jupiters: they tidally evolve towards a stationary solution. However, stellar rotation may also be affected by stellar wind braking. Thus, while the rotation of a quiet host star evolves towards a stationary attractor with a frequency ($1+6e^2$) times the orbital mean-motion of the companion, the continuous loss of angular momentum in an active star displaces the stationary solution towards slower values: Active host stars with big close-in companions tend to have rotational periods larger than the orbital periods of their companions. The study of some hypothetical examples shows that because of tidal evolution, the rules of gyrochronology cannot be used to estimate the age of one system with a large close-in companion, no matter if the star is quiet or active, if the current semi-major axis of the companion is smaller than 0.03--0.04 AU. Details on the evolution of the systems: CoRoT LRc06E21637, CoRoT-27, Kepler-75, CoRoT-2, CoRoT-18, CoRoT-14 and on hypothetical systems with planets of mass 1--4 M_Jup in orbit around a star similar to the Sun are given.

Context. We present the discovery of two transiting extrasolar planets by the satellite CoRoT. Aims. We aim at a characterization of the planetary bulk parameters, which allow us to further investigate the formation and evolution of the planetary systems and the main properties of the host stars. Methods. We used the transit light curve to characterize the planetary parameters relative to the stellar parameters. The analysis of HARPS spectra established the planetary nature of the detections, providing their masses. Further photometric and spectroscopic ground-based observations provided stellar parameters (log g,Teff,v sin i) to characterize the host stars. Our model takes the geometry of the transit to constrain the stellar density into account, which when linked to stellar evolutionary models, determines the bulk parameters of the star. Because of the asymmetric shape of the light curve of one of the planets, we had to include the possibility in our model that the stellar surface was not strictly spherical. Results. We present the planetary parameters of CoRoT-28b, a Jupiter-sized planet (mass 0.484+/-0.087MJup; radius 0.955+/-0.066RJup) orbiting an evolved star with an orbital period of 5.208 51 +/- 0.000 38 days, and CoRoT-29b, another Jupiter-sized planet (mass 0.85 +/- 0.20MJup; radius 0.90 +/- 0.16RJup) orbiting an oblate star with an orbital period of 2.850 570 +/- 0.000 006 days. The reason behind the asymmetry of the transit shape is not understood at this point. Conclusions. These two new planetary systems have very interesting properties and deserve further study, particularly in the case of the star CoRoT-29.

Non-Gaussianity in the primordial fluctuations that seeded structure formation produces a signal in the galaxy power spectrum on very large scales. This signal contains vital information about the primordial Universe, but it is very challenging to extract, because of cosmic variance and large-scale systematics - especially after the Planck experiment has already ruled out a large amplitude for the signal. Cosmic variance and experimental systematics can be alleviated by the multi-tracer method. Here we address another systematic - introduced by not using the correct relativistic analysis of the power spectrum on very large scales. In order to reduce the errors on fNL, we need to include measurements on the largest possible scales. Failure to include the relativistic effects on these scales can introduce significant bias in the best-fit value of fNL from future galaxy surveys.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

We produce simulations of emission of the atomic CII line in large sky fields in order to determine the current prospects for mapping this line during the high redshift Epoch of Reionization. We estimate the CII line intensity, redshift evolution and spatial fluctuations using observational relations between CII emission and the SFR in a galaxy for the frequency range of 200 GHz to 300 GHz. We obtained a frequency averaged intensity of CII emission of ${\rm I_{\rm CII}=(4 \pm 2)\times10^{2}\, Jy\, \rm sr^{-1}}$ in the redshift range $z\, \sim\, 5.3\, -\, 8.5$. Observations of CII emission in this frequency range will suffer contamination from emission lines at lower redshifts, in particular from the CO rotation lines. For the relevant frequency range we estimated the CO contamination (originated in emission from galaxies at $z\, <\, 2.5$), using simulations, to be ${\rm I_{\rm CO} \approx 10^{3}\, Jy \, sr^{-1}}$ and independently confirmed the result based in observational relations. We generated maps as a function of angle and frequency using detailed simulations of the CII and CO emission across several redshifts in order to properly take into account the observational pipeline and light cone effects. In order to reduce the foreground contamination we found that we should mask galaxies below redshifts $\sim 2.5$ with a CO flux in one of the CO(J:2-1) to CO(J:6-5) lines higher than ${\rm 5\times 10^{-22}\, W\ m^{-2}}$ or a AB magnitude lower than ${\rm m_{\rm K}\, =\, 22}$. We estimate that the additional continuum contamination is of the order of ${\rm 10^{5}\, Jy\, sr^{-1}}$. It is also considered the possibility of cross correlating foreground lines with galaxies in order to probe the intensity of the foregrounds.

We present a framework for forecasting cosmological constraints from future neutral hydrogen intensity mapping experiments at low to intermediate redshifts. In the process, we establish a simple way of comparing such surveys with optical galaxy redshift surveys. We explore a wide range of experimental configurations and assess how well a number of cosmological observables (the expansion rate, growth rate, and angular diameter distance) and parameters (the densities of dark energy and dark matter, spatial curvature, the dark energy equation of state, etc.) will be measured by an extensive roster of upcoming experiments. A number of potential contaminants and systematic effects are also studied in detail. The overall picture is encouraging -- if autocorrelation calibration can be controlled to a sufficient level, Phase I of the SKA should be able to constrain the dark energy equation of state about as well as a DETF Stage IV galaxy redshift survey like Euclid, in roughly the same timeframe.

The Square Kilometre Array (SKA) will conduct the biggest spectroscopic galaxy survey ever, by detecting the 21cm emission line of neutral hydrogen (HI) from around a billion galaxies over 3/4 of the sky, out to a redshift of z~2. This will allow the redshift-space matter power spectrum, and corresponding dark energy observables, to be measured with unprecedented precision. In this paper, we present an improved model of the HI galaxy number counts and bias from semi-analytic simulations, and use it to calculate the expected yield of HI galaxies from surveys with a variety of Phase 1 and 2 SKA configurations. We illustrate the relative performance of the different surveys by forecasting errors on the radial and transverse scales of the baryon acoustic oscillation (BAO) feature, finding that the full "billion galaxy survey" with SKA2 will deliver the largest dark energy figure of merit of any current or future large-scale structure survey.

Measuring radio source counts is critical for characterizing new extragalactic populations, brings a wealth of science within reach and will inform forecasts for SKA and its pathfinders. Yet there is currently great debate (and few measurements) about the behaviour of the 1.4-GHz counts in the microJy regime. One way to push the counts to these levels is via 'stacking', the covariance of a map with a catalogue at higher resolution and (often) a different wavelength. For the first time, we cast stacking in a fully bayesian framework, applying it to (i) the SKADS simulation and (ii) VLA data stacked at the positions of sources from the VIDEO survey. In the former case, the algorithm recovers the counts correctly when applied to the catalogue, but is biased high when confusion comes into play. This needs to be accounted for in the analysis of data from any relatively-low-resolution SKA pathfinders. For the latter case, the observed radio source counts remain flat below the 5-sigma level of 85 microJy as far as 40 microJy, then fall off earlier than the flux hinted at by the SKADS simulations and a recent P(D) analysis (which is the only other measurement from the literature at these flux-density levels, itself extrapolated in frequency). Division into galaxy type via spectral-energy distribution reveals that normal spiral galaxies dominate the counts at these fluxes.

The Planck CMB experiment has delivered the best constraints so far on primordial non-Gaussianity, ruling out early-Universe models of inflation that generate large non-Gaussianity. Although small improvements in the CMB constraints are expected, the next frontier of precision will come from future large-scale surveys of the galaxy distribution. The advantage of such surveys is that they can measure many more modes than the CMB -- in particular, forthcoming radio surveys with the SKA will cover huge volumes. Radio continuum surveys deliver the largest volumes, but with the disadvantage of no redshift information. In order to mitigate this, we use two additional observables. First, the integrated Sachs-Wolfe effect -- the cross-correlation of the radio number counts with the CMB temperature anisotropies -- helps to reduce systematics on the large scales that are sensitive to non-Gaussianity. Second, optical data allows for cross-identification in order to gain some redshift information. We show that, while the single redshift bin case can provide a sigma(fNL) ~ 20, and is therefore not competitive with current and future constraints on non-Gaussianity, a tomographic analysis could improve the constraints by an order of magnitude, even with only two redshift bins. A huge improvement is provided by the addition of high-redshift sources, so having cross-ID for high-z galaxies and an even higher-z radio tail is key to enabling very precise measurements of fNL. Our results show that SKA continuum surveys could provide constraints competitive with CMB and forthcoming optical surveys, potentially allowing a measurement of sigma(fNL) ~ 1 to be made. Moreover, these measurements would act as a useful check of results obtained with other probes at other redshift ranges with other methods.

We present DES13S2cmm, the first spectroscopically-confirmed superluminous supernova (SLSN) from the Dark Energy Survey (DES). We briefly discuss the data and search algorithm used to find this event in the first year of DES operations, and outline the spectroscopic data obtained from the European Southern Observatory (ESO) Very Large Telescope to confirm its redshift (z = 0.663 +/- 0.001 based on the host-galaxy emission lines) and likely spectral type (type I). Using this redshift, we find M_U_peak = -21.05 +0.10 -0.09 for the peak, rest-frame U-band absolute magnitude, and find DES13S2cmm to be located in a faint, low metallicity (sub-solar), low stellar-mass host galaxy (log(M/M_sun) = 9.3 +/- 0.3); consistent with what is seen for other SLSNe-I. We compare the bolometric light curve of DES13S2cmm to fourteen similarly well-observed SLSNe-I in the literature and find it possesses one of the slowest declining tails (beyond +30 days rest frame past peak), and is the faintest at peak. Moreover, we find the bolometric light curves of all SLSNe-I studied herein possess a dispersion of only 0.2-0.3 magnitudes between +25 and +30 days after peak (rest frame) depending on redshift range studied; this could be important for 'standardising' such supernovae, as is done with the more common type Ia. We fit the bolometric light curve of DES13S2cmm with two competing models for SLSNe-I - the radioactive decay of 56Ni, and a magnetar - and find that while the magnetar is formally a better fit, neither model provides a compelling match to the data. Although we are unable to conclusively differentiate between these two physical models for this particular SLSN-I, further DES observations of more SLSNe-I should break this degeneracy, especially if the light curves of SLSNe-I can be observed beyond 100 days in the rest frame of the supernova.

We present new calculations of unified line profiles for hydrogen perturbed by collisions with protons. We report on new calculations of the potential energies and dipole moments which allow the evaluation of profiles for the lines of the Lyman series up to Lyman$\delta$ and the Balmer series up to Balmer10. Unified calculations only existed for the lines Lyman$\alpha$ to Lyman$\gamma$ and Balmer$\alpha$ including the H$_2^+$ quasi-molecule. These data are available as online material accompanying this paper and should be included in atmosphere models, in place of the Stark effect of protons, since the quasi-molecular contributions cause not only satellites, but large asymmetries that are unaccounted for in models that assume Stark broadening of electrons and protons are equal.

We present the science enabled by cross-correlations of the SKA1-LOW 21-cm EoR surveys with other line mapping programs. In particular, we identify and investigate potential synergies with planned programs, such as the line intensity mapping of redshifted CO rotational lines, [CII] and Ly-$\alpha$ emissions during reionization. We briefly describe how these tracers of the star-formation rate at $z \sim 8$ can be modeled jointly before forecasting their auto- and cross-power spectra measurements with the nominal 21cm EoR survey. The use of multiple line tracers would be invaluable to validate and enrich our understanding of the EoR.

The exceptional sensitivity of the SKA will allow observations of the Cosmic Dawn and Epoch of Reionization (CD/EoR) in unprecedented detail, both spectrally and spatially. This wealth of information is buried under Galactic and extragalactic foregrounds, which must be removed accurately and precisely in order to reveal the cosmological signal. This problem has been addressed already for the previous generation of radio telescopes, but the application to SKA is different in many aspects. In this chapter we summarise the contributions to the field of foreground removal in the context of high redshift and high sensitivity 21-cm measurements. We use a state-of-the-art simulation of the SKA Phase 1 observations complete with cosmological signal, foregrounds and frequency-dependent instrumental effects to test both parametric and non-parametric foreground removal methods. We compare the recovered cosmological signal using several different statistics and explore one of the most exciting possibilities with the SKA --- imaging of the ionized bubbles. We find that with current methods it is possible to remove the foregrounds with great accuracy and to get impressive power spectra and images of the cosmological signal. The frequency-dependent PSF of the instrument complicates this recovery, so we resort to splitting the observation bandwidth into smaller segments, each of a common resolution. If the foregrounds are allowed a random variation from the smooth power law along the line of sight, methods exploiting the smoothness of foregrounds or a parametrization of their behaviour are challenged much more than non-parametric ones. However, we show that correction techniques can be implemented to restore the performances of parametric approaches, as long as the first-order approximation of a power law stands.

The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will have many advantages. For instance, the growth of perturbations is well understood for those modes, since it falls fully within the linear regime. Also, such scales are unaffected by the poorly understood feedback of baryonic physics. On ultra-large cosmic scales, two key effects become significant: primordial non-Gaussianity and relativistic corrections to cosmological observables. Moreover, if late-time acceleration is driven not by dark energy but by modifications to general relativity, then such modifications should become apparent near and above the horizon scale. As a result, the SKA is forecast to deliver transformational constraints on non-Gaussianity and to probe gravity on super-horizon scales for the first time.

Employing a nonparametric approach of the principal component analysis (PCA), we forecast the future constraint on the equation of state $w(z)$ of dark energy, and on the effective Newton constant $\mu(k,z)$, which parameterise the effect of modified gravity, using the planned SKA HI galaxy survey. Combining with the simulated data of Planck and Dark Energy Survey (DES), we find that SKA Phase 1 (SKA1) and SKA Phase 2 (SKA2) can well constrain $3$ and $5$ eigenmodes of $w(z)$ respectively. The errors of the best measured modes can be reduced to 0.04 and 0.023 for SKA1 and SKA2 respectively, making it possible to probe dark energy dynamics. On the other hand, SKA1 and SKA2 can constrain $7$ and $20$ eigenmodes of $\mu(k,z)$ respectively within 10\% sensitivity level. Furthermore, 2 and 7 modes can be constrained within sub percent level using SKA1 and SKA2 respectively. This is a significant improvement compared to the combined datasets without SKA.

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 6 - 27 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2.

In this chapter we provide an overview of the current status of the simulations and modelling of the Cosmic Dawn and Epoch of Reionization. We discuss the modelling requirements as dictated by the characteristic scales of the problem and the SKA instrumental properties and the planned survey parameters. Current simulations include most of the relevant physical processes. They can follow the full nonlinear dynamics and are now reaching the required scale and dynamic range, although small-scale physics still needs to be included at sub-grid level. However, despite a significant progress in developing novel numerical methods for efficient utilization of current hardware they remain quite computationally expensive. In response, a number of alternative approaches, particularly semi-analytical/semi-numerical methods, have been developed. While necessarily more approximate, if appropriately constructed and calibrated on simulations they could be used to quickly explore the vast parameter space available. Further work is still required on including some physical processes in both simulations and semi-analytical modelling. This hybrid approach of fast, approximate modelling calibrated on numerical simulations can then be used to construct large libraries of reionization models for reliable interpretation of the observational data.

HI intensity mapping (IM) is a novel technique capable of mapping the large-scale structure of the Universe in three dimensions and delivering exquisite constraints on cosmology, by using HI as a biased tracer of the dark matter density field. This is achieved by measuring the intensity of the redshifted 21cm line over the sky in a range of redshifts without the requirement to resolve individual galaxies. In this chapter, we investigate the potential of SKA1 to deliver HI intensity maps over a broad range of frequencies and a substantial fraction of the sky. By pinning down the baryon acoustic oscillation and redshift space distortion features in the matter power spectrum -- thus determining the expansion and growth history of the Universe -- these surveys can provide powerful tests of dark energy models and modifications to General Relativity. They can also be used to probe physics on extremely large scales, where precise measurements of spatial curvature and primordial non-Gaussianity can be used to test inflation; on small scales, by measuring the sum of neutrino masses; and at high redshifts where non-standard evolution models can be probed. We discuss the impact of foregrounds as well as various instrumental and survey design parameters on the achievable constraints. In particular we analyse the feasibility of using the SKA1 autocorrelations to probe the large-scale signal.

The Square Kilometer Array (SKA) has the potential to produce galaxy redshift surveys which will be competitive with other state of the art cosmological experiments in the next decade. In this chapter we summarise what capabilities the first and the second phases of the SKA will be able to achieve in its current state of design. We summarise the different cosmological experiments which are outlined in further detail in other chapters of this Science Book. The SKA will be able to produce competitive Baryonic Oscillation (BAOs) measurements in both its phases. The first phase of the SKA will provide similar measurements as optical and IR experiments with completely different systematic effects whereas the second phase being transformational in terms of its statistical power. The SKA will produce very accurate Redshift Space Distortions (RSD) measurements, being superior to other experiments at lower redshifts, due to the large number of galaxies. Cross correlations of the galaxy redshift data from the SKA with radio continuum surveys and optical surveys will provide extremely good calibration of photometric redshifts as well as extremely good bounds on modifications of gravity. Basing on a Principle Component Analysis (PCA) approach, we find that the SKA will be able to provide competitive constraints on dark energy and modified gravity models. Due to the large area covered the SKA it will be a transformational experiment in measuring physics from the largest scales such as non-Gaussian signals from $\textrm{f}_{\textrm{nl}}$. Finally, the SKA might produce the first real time measurement of the redshift drift. The SKA will be a transformational machine for cosmology as it grows from an early Phase 1 to its full power.

The new frontier of cosmology will be led by three-dimensional surveys of the large-scale structure of the Universe. Based on its all-sky surveys and redshift depth, the SKA is destined to revolutionize cosmology, in combination with future optical/ infrared surveys such as Euclid and LSST. Furthermore, we will not have to wait for the full deployment of the SKA in order to see transformational science. In the first phase of deployment (SKA1), all-sky HI intensity mapping surveys and all-sky continuum surveys are forecast to be at the forefront on the major questions of cosmology. We give a broad overview of the major contributions predicted for the SKA. The SKA will not only deliver precision cosmology -- it will also probe the foundations of the standard model and open the door to new discoveries on large-scale features of the Universe.

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

This chapter describes the assumed specifications and sensitivities for HI galaxy surveys with SKA1 and SKA2. It addresses the expected galaxy number densities based on available simulations as well as the clustering bias over the underlying dark matter. It is shown that a SKA1 HI galaxy survey should be able to find around $5\times 10^6$ galaxies over 5,000 deg$^2$ (up to $z\sim 0.8$), while SKA2 should find $\sim 10^9$ galaxies over 30,000 deg$^2$ (up to $z\sim 2.5$). The numbers presented here have been used throughout the cosmology chapters for forecasting.

The Square Kilometre Array (SKA) will offer an unprecedented view onto the early Universe, using interferometric observations of the redshifted 21cm line. The 21cm line probes the thermal and ionization state of the cosmic gas, which is governed by the birth and evolution of the first structures in our Universe. Here we show how the evolution of the 21cm signal will allow us to study when the first generations of galaxies appeared, what were their properties, and what was the structure of the intergalactic medium. We highlight qualitative trends which will offer robust insights into the early Universe.

The imprint of baryon acoustic oscillations (BAO) in large-scale structure can be used as a standard ruler for mapping out the cosmic expansion history, and hence for testing cosmological models. In this article we briefly describe the scientific background to the BAO technique, and forecast the potential of the Phase 1 and 2 SKA telescopes to perform BAO surveys using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys. We find that a 25,000 sq. deg. intensity mapping survey on a Phase 1 array will preferentially constrain the radial BAO, providing a highly competitive 2% constraint on the expansion rate at z ~ 2. A 30,000 sq. deg. galaxy redshift survey on SKA2 will outperform all other planned experiments for z < 1.4.

Radio continuum surveys have, in the past, been of restricted use in cosmology. Most studies have concentrated on cross-correlations with the cosmic microwave background to detect the integrated Sachs-Wolfe effect, due to the large sky areas that can be surveyed. As we move into the SKA era, radio continuum surveys will have sufficient source density and sky area to play a major role in cosmology on the largest scales. In this chapter we summarise the experiments that can be carried out with the SKA as it is built up through the coming decade. We show that the SKA can play a unique role in constraining the non-Gaussianity parameter to \sigma(f_NL) ~ 1, and provide a unique handle on the systematics that inhibit weak lensing surveys. The SKA will also provide the necessary data to test the isotropy of the Universe at redshifts of order unity and thus evaluate the robustness of the cosmological principle.Thus, SKA continuum surveys will turn radio observations into a central probe of cosmological research in the coming decades.

The peculiar motion of galaxies can be a particularly sensitive probe of gravitational collapse. As such, it can be used to measure the dynamics of dark matter and dark energy as well the nature of the gravitational laws at play on cosmological scales. Peculiar motions manifest themselves as an overall anisotropy in the measured clustering signal as a function of the angle to the line-of-sight, known as redshift-space distortion (RSD). Limiting factors in this measurement include our ability to model non-linear galaxy motions on small scales and the complexities of galaxy bias. The anisotropy in the measured clustering pattern in redshift-space is also driven by the unknown distance factors at the redshift in question, the Alcock-Paczynski distortion. This weakens growth rate measurements, but permits an extra geometric probe of the Hubble expansion rate. In this chapter we will briefly describe the scientific background to the RSD technique, and forecast the potential of the SKA phase 1 and the SKA2 to measure the growth rate using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys.

21cm intensity mapping experiments aim to observe the diffuse neutral hydrogen (HI) distribution on large scales which traces the Cosmic structure. The Square Kilometre Array (SKA) will have the capacity to measure the 21cm signal over a large fraction of the sky. However, the redshifted 21cm signal in the respective frequencies is faint compared to the Galactic foregrounds produced by synchrotron and free-free electron emission. In this article, we review selected foreground subtraction methods suggested to effectively separate the 21cm signal from the foregrounds with intensity mapping simulations or data. We simulate an intensity mapping experiment feasible with SKA phase 1 including extragalactic and Galactic foregrounds. We give an example of the residuals of the foreground subtraction with a independent component analysis and show that the angular power spectrum is recovered within the statistical errors on most scales. Additionally, the scale of the Baryon Acoustic Oscillations is shown to be unaffected by foreground subtraction.

We present new 2 cm and 6 cm maps of H2CO, radio recombination lines, and the radio continuum in the W51 star forming complex acquired with Arecibo and the Green Bank Telescope at ~50" resolution. We use H2CO absorption to determine the relative line-of-sight positions of molecular and ionized gas. We measure gas densities using the H2CO densitometer, including continuous measurements of the dense gas mass fraction (DGMF) over the range $10^4$ cm$^{-3}$ < n(H$_2$) < $10^6$ cm$^{-3}$- this is the first time a dense gas mass fraction has been measured over a range of densities with a single data set. The DGMF in W51A is high,f >~70% above $n>10^4$ cm$^{-3}$, while it is low, f<20%, in W51 B. We did not detect any H2CO emission throughout the W51 GMC; all gas dense enough to emit under normal conditions is in front of bright continuum sources and therefore is seen in absorption instead. The data set has been made public at http://dx.doi.org/10.7910/DVN/26818. Conclusions. (1) The dense gas fraction in the W51 A and B clouds shows that W51 A will continue to form stars vigorously, while star formation has mostly ended in W51 B. The lack of dense, star-forming gas around W51 C indicates that collect-and-collapse is not acting or is inefficient in W51. (2) Ongoing high-mass star formation is correlated with n ~ 1x10$^5$ cm$^{-3}$ gas. Gas with n > 10$^4$ cm$^{-3}$ is weakly correlated with low and moderate mass star formation, but does not strongly correlate with high-mass star formation. (3) The nondetection of H$_2$CO emission implies that the emission detected in other galaxies, e.g. Arp 220, comes from high-density gas that is not directly affiliated with already-formed massive stars. Either the non-star-forming ISM of these galaxies is very dense, implying the star formation density threshold is higher, or H II regions have their emission suppressed.

Astronomy at or below the 'survey threshold' has expanded significantly since the publication of the original 'Science with the Square Kilometer Array' in 1999 and its update in 2004. The techniques in this regime may be broadly (but far from exclusively) defined as 'confusion' or 'P(D)' analyses (analyses of one-point statistics), and 'stacking', accounting for the flux-density distribution of noise-limited images co-added at the positions of objects detected/isolated in a different waveband. Here we discuss the relevant issues, present some examples of recent analyses, and consider some of the consequences for the design and use of surveys with the SKA and its pathfinders.

The creation of cold dark matter cosmology model is studied beyond the linear perturbation level. The skewness is explicitly computed and the results are compared to those from the $\Lambda$CDM model. It is explicitly shown that both models have the same signature for the skewness and cannot be distinguished by using this observable.

It has been suggested that HD 200964 is the first exoplanetary system with two Jovian planets evolving in the 4/3 mean- motion resonance. Previous scenarios to simulate the formation of two giant planets in the stable 4/3 resonance configuration have failed. Moreover, the orbital parameters available in the literature point out an unstable configuration of the planetary pair. The purpose of this paper is i) to determine the orbits of the planets from the RV measurements and update the value of the stellar mass (1.57 M), ii) to analyse the stability of the planetary evolution in the vicinity and inside the 4/3 MMR, and iii) to elaborate a possible scenario for the formation of systems in the 4/3 MMR. The results of the formation simulations are able to very closely reproduce the 4/3 resonant dynamics of the best-fit config- uration obtained in this paper. Moreover, the confidence interval of the fit matches well with the very narrow stable region of the 4/3 mean-motion resonance. The formation process of the HD 200964 system is very sensitive to the planetary masses and protoplanetary disk parameters. Only a thin, flat disk allows the embryo-sized planets to reach the 4/3 resonant configuration. The stable evolution of the resonant planets is also sensitive to the mass of the central star, because of overlapping high-order resonances inside the 4/3 resonance. Regardless of the very narrow domain of stable motion, the confidence interval of our fit closely matches the stability area.

We make use of a large set of fast simulations of an intensity mapping experiment with characteristics similar to those expected of the Square Kilometre Array (SKA) in order to study the viability and limits of blind foreground subtraction techniques. In particular, we consider different approaches: polynomial fitting, principal component analysis (PCA) and independent component analysis (ICA). We review the motivations and algorithms for the three methods, and show that they can all be described, using the same mathematical framework, as different approaches to the blind source separation problem. We study the efficiency of foreground subtraction both in the angular and radial (frequency) directions, as well as the dependence of this efficiency on different instrumental and modelling parameters. For well-behaved foregrounds and instrumental effects we find that foreground subtraction can be successful to a reasonable level on most scales of interest. We also quantify the effect that the cleaning has on the recovered signal and power spectra. Interestingly, we find that the three methods yield quantitatively similar results, with PCA and ICA being almost equivalent.

We present a code to generate mock observations of 21 cm intensity mapping experiments. The emphasis of the code is on reducing the computational cost of running a full-blown simulation, trading computational time for accuracy. The code can be used to generate independent realizations of the cosmological signal and foregrounds, which are necessary, for instance, in order to obtain realistic forecasts for future intensity mapping experiments. The code is able to reproduce the correct angular and radial clustering pattern for the cosmological signal, including redshift-space distortions, lightcone evolution and bias. Furthermore, it is possible to simulate a variety of foregrounds, including the potentially problematic polarized synchrotron emission.

We present the discovery of a candidate multiply-transiting system, the first one found in the CoRoT mission. Two transit-like features with periods of 5.11 and 11.76d are detected in the CoRoT light curve, around a main sequence K1V star of r=15.1. If the features are due to transiting planets around the same star, these would correspond to objects of 3.7$\pm$0.4 and 5.0$\pm$0.5 R_earth respectively. Several radial velocities serve to provide an upper limit of 5.7 M_earth for the 5.11~d signal, and to tentatively measure a mass of 28$^{+11}_{-11}$ M_earth for the object transiting with a 11.76~d period. These measurements imply low density objects, with a significant gaseous envelope. The detailed analysis of the photometric and spectroscopic data serve to estimate the probability that the observations are caused by transiting Neptune-sized planets as $>$26$\times$ higher than a blend scenario involving only one transiting planet, and $>$900$\times$ higher than a scenario involving two blends and no planets. The radial velocities show a long term modulation that might be attributed to a 1.5 M_jup planet orbiting at 1.8~A.U. from the host, but more data are required to determine the precise orbital parameters of this companion.

We investigate the creation of cold dark matter (CCDM) cosmology as an alternative to explain the cosmic acceleration. Particular attention is given to the evolution of density perturbations and constraints coming from recent observations. By assuming negligible effective sound speed we compare CCDM predictions with redshift-space-distortion based f(z) sigma_8(z) measurements. We identify a subtle issue associated with which contribution in the density contrast should be used in this test and then show that the CCDM results are the same as those obtained with LambdaCDM. These results are then contrasted with the ones obtained at the background level. For the background tests we have used type Ia supernovae data (Union 2.1 compilation) in combination with baryonic acoustic oscillations and cosmic microwave background observations and also measurements of the Hubble parameter at different redshifts. As a consequence of the studies we have performed at both the background and perturbation levels, we explicitly show that CCDM is observationally degenerate with respect to LambdaCDM (dark degeneracy). The need to overcome the lack of a fundamental microscopic basis for the CCDM is the major challenge for this kind of model.

The intensity mapping of Ly-alpha emission during the epoch of reionization (EoR) will be contaminated by foreground emission lines from lower redshifts. We calculate the mean intensity and power spectrum of Ly-alpha emission at z~7, and estimate the uncertainties according to the relevant astrophysical processes. We find that the low-redshift emission lines from 6563 A H-alpha, 5007 A [OIII] and 3727 A [OII] will be the strong contaminants on the observed Ly-alpha power spectrum. We make use of both the star formation rate (SFR) and luminosity functions (LF) to estimate the mean intensity and power spectra of the three foreground lines at z~0.5 for H-alpha, z~0.9 for [OIII] and z~1.6 for [OII], as they will contaminate the Ly-alpha emission at z~7. The [OII] line is found to be the strongest. We analyze the masking of the bright survey pixels with a foreground line above some line intensity threshold as a way to reduce the contamination in the intensity mapping survey. We find that the foreground contamination can be neglected if we remove the pixels with fluxes above 1.4x10^-20 W/m^2.

We present a large-scale clustering analysis of radio galaxies in the Very Large Array (VLA) Faint Images of the Radio Sky at Twenty-cm (FIRST) survey over the Galaxy And Mass Assembly (GAMA) survey area, limited to S1.4 GHz >1 mJy with spectroscopic and photometric redshift limits up to r < 19.8 and r < 22 mag, respectively. For the GAMA spectroscopic matches, we present the redshift-space and projected correlation functions, the latter of which yielding a correlation length r0 ~ 8.2 Mpc/h and linear bias of ~1.9 at z ~ 0.34. Furthermore, we use the angular two-point correlation function w({\theta}) to determine spatial clustering properties at higher redshifts. We find r0 to increase from ~6 to ~14 Mpc/h between z = 0.3 and z = 1.55, with the corresponding bias increasing from ~2 to ~10 over the same range. Our results are consistent with the bias prescription implemented in the SKADS simulations at low redshift, but exceed these predictions at z > 1. This is indicative of an increasing (rather than fixed) halo mass and/or AGN fraction at higher redshifts or a larger typical halo mass for the more abundant FRI sources.

We explore the use of different radio galaxy populations as tracers of different mass halos and therefore, with different bias properties, to constrain primordial non-Gaussianity of the local type. We perform a Fisher matrix analysis based on the predicted auto and cross angular power spectra of these populations, using simulated redshift distributions as a function of detection flux and the evolution of the bias for the different galaxy types (Star forming galaxies, Starburst galaxies, Radio-Quiet Quasars, FRI and FRII AGN galaxies). We show that such a multi-tracer analysis greatly improves the information on non-Gaussianity by drastically reducing the cosmic variance contribution to the overall error budget. By using this method applied to future surveys, we predict a constraint of sigma_fnl=3.6 on the local non-Gaussian parameter for a galaxy detection flux limit of 10 \muJy and sigma_fnl=2.2 for 1 \muJy. We show that this significantly improves on the constraints obtained when using the whole undifferentiated populations (sigma_fnl=48 for 10 \muJy and sigma_fnl=12 for 1 \muJy). We conclude that continuum radio surveys alone have the potential to constrain primordial non-Gaussianity to an accuracy at least a factor of two better than the present constraints obtained with Planck data on the CMB bispectrum, opening a window to obtain sigma_fnl~1 with the Square Kilometer Array.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

The large-scale structure of the Universe supplies crucial information about the physical processes at play at early times. Unresolved maps of the intensity of 21 cm emission from neutral hydrogen HI at redshifts z~1-5 are the best hope of accessing the ultralarge-scale information, directly related to the early Universe. A purpose-built HI intensity experiment may be used to detect the large scale effects of primordial non-Gaussianity, placing stringent bounds on different models of inflation. We argue that it may be possible to place tight constraints on the non-Gaussianity parameter f_NL, with an error close to ~1.

We present a statistical method based on a maximum likelihood approach to constrain the number counts of extragalactic sources below the nominal flux-density limit of continuum imaging surveys. We extract flux densities from a radio map using positional information from an auxiliary catalogue and show that we can model the number counts of this undetected population down to flux density levels well below the detection threshold of the radio survey. We demonstrate the capabilities that our method will have with future generation wide-area radio surveys by performing simulations over various sky areas. We show that it is possible to accurately constrain the number counts of the simulated distribution down to one tenth of the flux noise rms with just a sky area of 100 deg$^2$. We then test the application of our method using data from the Faint Images of the Radio Sky at Twenty-Centimeters survey (FIRST). We extract flux densities from the FIRST map, sensitive to 150 $\mu$Jy/beam (1 $\sigma$), using the positional information from a catalogue in the same field, also acquired at the same frequency, sensitive to 12 $\mu$Jy/beam (1 $\sigma$). Implementing our method, with known source positions, we are able to recover the right differential number counts of the noise-dominated FIRST map fluxes down to a flux density level which is one-tenth the FIRST detection threshold.

One of the most promising ways to study the epoch of reionization (EoR) is through radio observations of the redshifted 21-cm line emission from neutral hydrogen. These observations are complicated by the fact that the mapping of redshifts to line-of-sight positions is distorted by the peculiar velocities of the gas. Such distortions can be a source of error if they are not properly understood, but they also encode information about cosmology and astrophysics. We study the effects of redshift space distortions on the power spectrum of 21-cm radiation from the EoR using large scale $N$-body and radiative transfer simulations. We quantify the anisotropy introduced in the 21-cm power spectrum by redshift space distortions and show how it evolves as reionization progresses and how it relates to the underlying physics. We go on to study the effects of redshift space distortions on LOFAR observations, taking instrument noise and foreground subtraction into account. We find that LOFAR should be able to directly observe the power spectrum anisotropy due to redshift space distortions at spatial scales around $k \sim 0.1$ Mpc$^{-1}$ after $\gtrsim$ 1000 hours of integration time. At larger scales, sample errors become a limiting factor, while at smaller scales detector noise and foregrounds make the extraction of the signal problematic. Finally, we show how the astrophysical information contained in the evolution of the anisotropy of the 21-cm power spectrum can be extracted from LOFAR observations, and how it can be used to distinguish between different reionization scenarios.

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

We calculate the absolute intensity and anisotropies of the Lyman-alpha radiation field present during the epoch of reionization. We consider emission from both galaxies and the intergalactic medium (IGM) and take into account the main contributions to the production of Lyman-alpha photons: recombinations, collisions, continuum emission from galaxies and scattering of Lyman-n photons in the IGM. We find that the emission from individual galaxies dominates over the IGM with a total Lyman-alpha intensity (times frequency) of about (1.43-3.57)x10^{-8} erg s^{-1} cm^{-2} sr^{-1} at a redshift of 7. This intensity level is low so it is unlikely that the Lyman-\alpha background during reionization can be established by an experiment aiming at an absolute background light measurement. Instead we consider Lyman-\alpha intensity mapping with the aim of measuring the anisotropy power spectrum which has rms fluctuations at the level of 1 x 10^{-16} [erg s^[-1} cm^{-2} sr^{-1}]^2 at a few Mpc scales. These anisotropies could be measured with a spectrometer at near-IR wavelengths from 0.9 to 1.4 \mu m with fields in the order of 0.5 to 1 sq. degrees. We recommend that existing ground-based programs using narrow band filters also pursue intensity fluctuations to study statistics on the spatial distribution of faint Lyman-\alpha emitters. We also discuss the cross-correlation signal with 21 cm experiments that probe HI in the IGM during reionization. A dedicated sub-orbital or space-based Lyman-\alpha intensity mapping experiment could provide a viable complimentary approach to probe reionization, when compared to 21 cm experiments, and is likely within experimental reach.

Molecular hydrogen is now understood to be the main coolant of the primordial gas clouds leading to the formation of the very first stars and galaxies. The line emissions associated with molecular hydrogen should then be a good tracer of the matter distribution at the onset of reionization of the universe. Here we propose intensity mapping of H2 line emission in rest-frame mid-infrared wavelengths to map out the spatial distribution of gas at redshifts z > 10. We calculate the expected mean intensity and clustering power spectrum for several H2 lines. We find that the 0-0S(3) rotational line at a rest wavelength of 9.66 microns is the brightest line over the redshift range of 10 to 30 with an intensity of about 5 to 10 Jy/sr at z~15. To reduce astrophysical and instrumental systematics, we propose the cross-correlation between multiple lines of the H2 rotational and vibrational line emission spectrum. Our estimates of the intensity can be used as a guidance in planning instruments for future mid-IR spectroscopy missions such as SPICA.

In this paper, we explore how the forthcoming generation of large-scale radio continuum surveys, with the inclusion of some degree of redshift information, can constrain cosmological parameters. By cross-matching these radio surveys with shallow optical to near-infrared surveys, we can essentially separate the source distribution into a low- and a high-redshift sample, thus providing a constraint on the evolution of cosmological parameters such as those related to dark energy. We examine two radio surveys, the Evolutionary Map of the Universe (EMU) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). A crucial advantage is their combined potential to provide a deep, full-sky survey. The surveys used for the cross-identifications are SkyMapper and SDSS, for the southern and northern skies, respectively. We concentrate on the galaxy clustering angular power spectrum as our benchmark observable, and find that the possibility of including such low redshift information yields major improvements in the determination of cosmological parameters. With this approach, and provided a good knowledge of the galaxy bias evolution, we are able to put strict constraints on the dark energy parameters, i.e. w_0=-0.9+/-0.041 and w_a=-0.24+/-0.13, with type Ia supernovae and CMB priors (with a one-parameter bias in this case); this corresponds to a Figure of Merit (FoM) > 600, which is twice better than what is obtained by using only the cross-identified sources and greater than four time better than the case without any redshift information at all.

A fraction of the extragalactic near-infrared (near-IR) background light involves redshifted photons from the ultraviolet (UV) emission from galaxies present during reionization at redshifts above 6. The absolute intensity and the anisotropies of the near-IR background provide an observational probe of the first-light galaxies and their spatial distribution. We estimate the extragalactic background light intensity during reionization by accounting for the stellar and nebular emission from first-light galaxies. We require the UV photon density from these galaxies to generate a reionization history that is consistent with the optical depth to electron scattering from cosmic microwave background measurements. We also require the bright-end luminosity function of galaxies in our models to reproduce the measured Lyman drop-out luminosity functions at redshifts of 6 to 8. The absolute intensity is about 0.1 to 0.3 nW m$^{-2}$ sr$^{-1}$ at the peak of its spectrum at $\sim$ 1.1 $\mu$m. We also discuss the anisotropy power spectrum of the near-IR background using a halo model to describe the galaxy distribution. We compare our predictions for the anisotropy power spectrum to existing measurements from deep near-IR imaging data from {\it Spitzer}/IRAC, {\it Hubble}/NICMOS, and {\it AKARI}. The predicted rms fluctuations at tens of arcminute angular scales are roughly an order of magnitude smaller than the existing measurements. While strong arguments have been made that the measured fluctuations do not have an origin involving faint low-redshift galaxies, we find that the existing measurements are also incompatible with an origin during the era of reionization. The measured near-IR background anisotropies remain unexplained and could be associated with an unidentified non-astrophysical origin.

We present theoretical calculations of the quasi-molecular line profiles for the Lyman (Lya, Lyb, Lyg, Lyd) and Balmer (Ha, Hb, Hg, Hd, He, H8, H9, and H10) series, perturbed by collisions with protons. In all calculations we include the dependence of the dipole moments as a function of internuclear distance during the collision. The broadening from ion collisions must be added to the normal electron Stark broadening.

We have done a new analysis of the available observations for the GJ581 exoplanetary system. Today this system is controversial due to choices that can be done in the orbital determination. The main ones are the ocurrence of aliases and the additional bodies - the planets f and g - announced in Vogt et al. 2010. Any dynamical study of exoplanets requires the good knowledge of the orbital elements and the investigations involving the planet g are particularly interesting, since this body would lie in the Habitable Zone (HZ) of the star GJ581. This region,for this system, is very attractive of the dynamical point of view due to several resonances of two and three bodies present there. In this work, we investigate the conditions under which the planet g may exist. We stress the fact that the planet g is intimately related with the orbital elements of the planet d; more precisely, we conclude that it is not possible to disconnect its existence from the determination of the eccentricity of the planet d. Concerning the planet f, we have found one solution with period $\approx 450$ days, but we are judicious about any affirmation concernig this body because its signal is in the threshold of detection and the high period is in a spectral region where the ocorruence of aliases is very common. Besides, we outline some dynamical features of the habitable zone with the dynamical map and point out the role played by some resonances laying there.

The atomic CII fine-structure line is one of the brightest lines in a typical star-forming galaxy spectrum with a luminosity ~ 0.1% to 1% of the bolometric luminosity. It is potentially a reliable tracer of the dense gas distribution at high redshifts and could provide an additional probe to the era of reionization. By taking into account of the spontaneous, stimulated and collisional emission of the CII line, we calculate the spin temperature and the mean intensity as a function of the redshift. When averaged over a cosmologically large volume, we find that the CII emission from ionized carbon in individual galaxies is larger than the signal generated by carbon in the intergalactic medium (IGM). Assuming that the CII luminosity is proportional to the carbon mass in dark matter halos, we also compute the power spectrum of the CII line intensity at various redshifts. In order to avoid the contamination from CO rotational lines at low redshift when targeting a CII survey at high redshifts, we propose the cross-correlation of CII and 21-cm line emission from high redshifts. To explore the detectability of the CII signal from reionization, we also evaluate the expected errors on the CII power spectrum and CII-21 cm cross power spectrum based on the design of the future milimeter surveys. We note that the CII-21 cm cross power spectrum contains interesting features that captures physics during reionization, including the ionized bubble sizes and the mean ionization fraction, which are challenging to measure from 21-cm data alone. We propose an instrumental concept for the reionization CII experiment targeting the frequency range of $\sim$ 200 to 300 GHz with 1, 3 and 10 meter apertures and a bolometric spectrometer array with 64 independent spectral pixels with about 20,000 bolometers.

The large-scale structure of the Universe can be mapped with unresolved intensity fluctuations of the 21 cm line. The power spectrum of the intensity fluctuations has been proposed as a probe of the baryon acoustic oscillations at low to moderate redshifts with interferometric experiments now under consideration. We discuss the contamination to the low-redshift 21 cm intensity power spectrum generated by the 18 cm OH line since the intensity fluctuations of the OH line generated at a slightly higher redshift contribute to the intensity fluctuations observed in an experiment. We assume the OH megamaser luminosity is correlated with the star formation rate, and use the simulation to estimate the OH signal and the spatial anisotropies. We also use a semi-analytic simulation to predict the 21 cm power spectrum. At z=1 to 3, we find that the OH contamination could reach 0.1% to 1% of the 21 cm rms fluctuations at the scale of the first peak of the baryon acoustic oscillation. When z>3 the OH signal declines quickly, so that the contamination on the 21-cm becomes negligible at high redshifts.

Primordial non-Gaussianity is a crucial test of inflationary cosmology. We consider the impact of non-Gaussianity on the ionization power spectrum from 21 cm emission during the epoch of reionization. We focus on the power spectrum on large scales at redshifts of 7 to 8 and explore the expected constraint on the local non-Gaussianity parameter f_NL for current and next-generation 21 cm experiments. We show that experiments such as SKA and MWA could measure f_NL values of order 10. This can be improved by an order of magnitude with a fast-Fourier transform telescope like Omniscope.

We present forecasts for constraints on cosmological models which can be obtained by forthcoming radio continuum surveys: the wide surveys with the LOw Frequency ARray (LOFAR), Australian Square Kilometre Array Pathfinder (ASKAP) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). We use simulated catalogues appropriate to the planned surveys to predict measurements obtained with the source auto-correlation, the cross-correlation between radio sources and CMB maps (the Integrated Sachs-Wolfe effect), the cross-correlation of radio sources with foreground objects due to cosmic magnification, and a joint analysis together with the CMB power spectrum and supernovae. We show that near future radio surveys will bring complementary measurements to other experiments, probing different cosmological volumes, and having different systematics. Our results show that the unprecedented sky coverage of these surveys combined should provide the most significant measurement yet of the Integrated Sachs-Wolfe effect. In addition, we show that using the ISW effect will significantly tighten constraints on modified gravity parameters, while the best measurements of dark energy models will come from galaxy auto-correlation function analyses. Using the combination of EMU and WODAN to provide a full sky survey, it will be possible to measure the dark energy parameters with an uncertainty of \{$\sigma (w_0) = 0.05$, $\sigma (w_a) = 0.12$\} and the modified gravity parameters \{$\sigma (\eta_0) = 0.10$, $\sigma (\mu_0) = 0.05$\}, assuming Planck CMB+SN(current data) priors. Finally, we show that radio surveys would detect a primordial non-Gaussianity of $f_{\rm NL}$ = 8 at 1-$\sigma$ and we briefly discuss other promising probes.

This investigation uses the excellent HARPS radial velocity measurements of CoRoT-7 to re-determine the planet masses and to explore techniques able to determine mass and elements of planets discovered around active stars when the relative variation of the radial velocity due to the star activity cannot be considered as just noise and can exceed the variation due to the planets. The main technique used here is a self-consistent version of the high-pass filter used by Queloz et al. (2009) in the first mass determination of CoRoT-7b and CoRoT-7c. The results are compared to those given by two alternative techniques: (1) The approach proposed by Hatzes et al. (2010) using only those nights in which 2 or 3 observations were done; (2) A pure Fourier analysis. In all cases, the eccentricities are taken equal to zero as indicated by the study of the tidal evolution of the system; the periods are also kept fixed at the values given by Queloz et al. Only the observations done in the time interval BJD 2,454,847 - 873 are used because they include many nights with multiple observations; otherwise it is not possible to separate the effects of the rotation fourth harmonic (5.91d = Prot/4) from the alias of the orbital period of CoRoT-7b (0.853585 d). The results of the various approaches are combined to give for the planet masses the values 8.0 \pm 1.2 MEarth for CoRoT-7b and 13.6 \pm 1.4 MEarth for CoRoT 7c. An estimation of the variation of the radial velocity of the star due to its activity is also given.The results obtained with 3 different approaches agree to give masses larger than those in previous determinations. From the existing internal structure models they indicate that CoRoT-7b is a much denser super-Earth. The bulk density is 11 \pm 3.5 g.cm-3 . CoRoT-7b may be rocky with a large iron core.

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

Observations of anisotropies in the brightness temperature of the 21 cm line of neutral hydrogen from the period before reionization would shed light on the dawn of the first stars and galaxies. In this paper, we use large-scale semi-numerical simulations to analyse the imprint on the 21 cm signal of spatial fluctuations in the Lyman-alpha flux arising from the clustering of the first galaxies. We show that an experiment such as the Square Kilometer Array (SKA) can probe this signal at the onset of reionization, giving us important information about the UV emission spectra of the first stars and characterizing their host galaxies. SKA-pathfinders with ~ 10% of the full collecting area should be capable of making a statistical detection of the 21 cm power spectrum at redshifts z < 20 (corresponding to frequencies $\nu$ > 67 MHz). We then show that the SKA should be able to measure the three dimensional power spectrum as a function of the angle with the line of sight and discuss the use of the redshift space distortions as a way to separate out the different components of the 21 cm power spectrum. We demonstrate that, at least on large scales where the Lyman-alpha fluctuations are linear, they can be used as a model independent way to extract the power spectra due to these Lyman-alpha fluctuations.

We propose observations of the molecular gas distribution during the era of reionization. At $z \sim $ 6 to 8, $^{12}$CO$(J=1-0)$ line intensity results in a mean brightness temperature of about 0.5 $\mu$K with rms fluctuations of 0.1 $\mu$K at 1 to 10 Mpc spatial scales, corresponding to 30 arcminute angular scales. This intensity fluctuations can be mapped with an interferometer, similar to existing and planned 21-cm background experiments, but operating at $\sim$ 12 to 17 GHz. We discuss the feasibility of detecting the cross-correlation between HI and CO molecular gas since such a cross-correlation has the advantage that it will be independent of systematics and most foregrounds in each of the 21-cm and CO$(1-0)$ line experiments. Additional instruments tuned to higher-order transitions of the CO molecule or an instrument operating with high spectral resolution at mm-wavelengths targeting 158 $\mu$m CII could further improve the reionization studies with molecular gas. The combined 21-cm and CO line observations has the potential to establish the relative distribution of gas in the inter-galactic medium and molecular gas that are clumped in individual first-light galaxies that are closely connected to the formation of massive stars in these galaxies.

The temperature fluctuations and polarization of the Cosmic Microwave Background (CMB) are now a well-known probe of the Universe at an infant age of 400,000 years. During the transit to us from the surface of last scattering, the CMB photons are expected to undergo modifications induced by the intervening large-scale structure. Among the expected secondary effects is the weak gravitational lensing of the CMB by the foreground dark matter distribution. We derive a quadratic estimator that uses the non-Gaussianities generated by the lensing effect at the four-point function level to extract the power spectrum of lensing potential fluctuations integrated out to z ~ 1100 with peak contributions from potential fluctuations at z of 2 to 3. Using WMAP 7-year temperature maps, we report the first direct constraints of this lensing potential power spectrum and find that it has an amplitude of A_L = 0.96 +/- 0.60, 1.06 +/- 0.69 and 0.97 +/- 0.47 using the W, V and W+V bands, respectively.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

While limited to low spatial resolution, the next generation low-frequency radio interferometers that target 21 cm observations during the era of reionization and prior will have instantaneous fields-of-view that are many tens of square degrees on the sky. Predictions related to various statistical measurements of the 21 cm brightness temperature must then be pursued with numerical simulations of reionization with correspondingly large volume box sizes, of order 1000 Mpc on one side. We pursue a semi-numerical scheme to simulate the 21 cm signal during and prior to Reionization by extending a hybrid approach where simulations are performed by first laying down the linear dark matter density field, accounting for the non-linear evolution of the density field based on second-order linear perturbation theory as specified by the Zel'dovich approximation, and then specifying the location and mass of collapsed dark matter halos using the excursion-set formalism. The location of ionizing sources and the time evolving distribution of ionization field is also specified using an excursion-set algorithm. We account for the brightness temperature evolution through the coupling between spin and gas temperature due to collisions, radiative coupling in the presence of Lyman-alpha photons and heating of the intergalactic medium, such as due to a background of X-ray photons. The hybrid simulation method we present is capable of producing the required large volume simulations with adequate resolution in a reasonable time so a large number of realizations can be obtained with variations in assumptions related to astrophysics and background cosmology that govern the 21 cm signal.

We estimate the conditions for detectability of two planets in a 2/1 mean-motion resonance from radial velocity data, as a function of their masses, number of observations and the signal-to-noise ratio. Even for a data set of the order of 100 observations and standard deviations of the order of a few meters per second, we find that Jovian-size resonant planets are difficult to detect if the masses of the planets differ by a factor larger than $\sim 4$. This is consistent with the present population of real exosystems in the 2/1 commensurability, most of which have resonant pairs with similar minimum masses, and could indicate that many other resonant systems exist, but are presently beyond the detectability limit. Furthermore, we analyze the error distribution in masses and orbital elements of orbital fits from synthetic data sets for resonant planets in the 2/1 commensurability. For various mass ratios and number of data points we find that the eccentricity of the outer planet is systematically over estimated, although the inner planet's eccentricity suffers a much smaller effect. If the initial conditions correspond to small amplitude oscillations around stable apsidal corotation resonances (ACR), the amplitudes estimated from the orbital fits are biased toward larger amplitudes, in accordance to results found in real resonant extrasolar systems.

We measure several properties of the reionization process and the corresponding low-frequency 21-cm signal associated with the neutral hydrogen distribution, using a large volume, high resolution simulation of cosmic reionization. The brightness temperature of the 21-cm signal is derived by post-processing this numerical simulation with a semi-analytical prescription. Our study extends to high redshifts (z ~ 25) where, in addition to collisional coupling, our post-processed simulations take into account the inhomogeneities in the heating of the neutral gas by X-rays and the effect of an inhomogeneous Lya radiation field. Unlike the well-studied case where spin temperature is assumed to be significantly greater than the temperature of the cosmic microwave background due to uniform heating of the gas by X-rays, spatial fluctuations in both the Lya radiation field and X-ray intensity impact predictions related to the brightness temperature at z > 10, during the early stages of reionization and gas heating. The statistics of the 21-cm signal from our simulation are then compared to existing analytical models in the literature and we find that these analytical models provide a reasonably accurate description of the 21-cm power spectrum at z < 10. Such an agreement is useful since analytical models are better suited to quickly explore the full astrophysical and cosmological parameter space relevant for future 21-cm surveys. We find, nevertheless, non-negligible differences that can be attributed to differences in the inhomogeneous X-ray heating and Lya coupling at z > 10 and, with upcoming interferometric data, these differences in return can provide a way to better understand the astrophysical processes during reionization.

A number of recent numerical investigations concluded that the remnants of rare structures formed at very high redshift, such as the very first stars and bright redshift z~6 QSOs, are preferentially located at the center of the most massive galaxy clusters at redshift z=0. In this paper we readdress this question using a combination of cosmological simulations of structure formation and extended Press-Schechter formalism and we show that the typical remnants of Population III stars are instead more likely to be found in a group environment, that is in dark matter halos of mass ~2x10^{13} h^{-1}M_sun. Similarly, the descendants of the brightest z~6 QSOs are expected to be in medium-sized clusters (mass of a few 10^{14} h^{-1}M_sun), rather than in the most massive superclusters (M>10^{15} h^{-1}M_sun) found within the typical 1 Gpc^3 cosmic volume where a bright z~6 QSO lives. The origin of past claims that the most massive clusters preferentially host these remnants is rooted in the numerical method used to initialize their numerical simulations: Only a small region of the cosmological volume of interest was simulated with sufficient resolution to identify low-mass halos at early times, and this region was chosen to host the most massive halo in the cosmological volume at late times. The conclusion that the earliest structures formed in the entire cosmological volume evolve into the most massive halo at late times was thus arrived at by construction. We demonstrate that, to the contrary, the first structures to form in a cosmological region evolve into relatively typical objects at later times. We propose alternative numerical methods for simulating the earliest structures in cosmological volumes.

Star-forming galaxies constitute the majority of galaxies with stellar masses >10^10 M_Sun/h^2 at z~2 and dominate the star-formation rate density of the Universe at this early epoch. It is thus critical to understand their origins, evolution, and connection to the underlying dark matter distribution. To this end, we identify the dark matter halos (including subhalos) that are likely to contain star-forming galaxies at z~2 (z2SFGs) within a large dissipationless cosmological simulation and then use halo merger histories to follow the evolution of z2SFG descendants to z~1 and z~0. The evolved halos at these epochs are then confronted with an array of observational data in order to uncover the likely descendants of z2SFGs. Though the evolved halos have clustering strengths comparable to red galaxies at z~1 and z~0, we find that the bulk of z2SFGs do not evolve into red galaxies, at either epoch. This conclusion is based primarily on the fact that the space density of z2SFGs is much higher than that of lower redshift red galaxies, even when accounting for the merging of z2SFG descendants, which decreases the number density of z2SFG descendants by at most a factor of two by z~0. Of the ~50% of z2SFGs that survive to z~0, ~70% reside at the center of z~0 dark matter halos with M>10^12 M_Sun/h. Halo occupation modeling of z~0 galaxies suggests that such halos are occupied by galaxies with M_r<-20.5, implying that these z2SFGs evolve into ``typical'' ~ L* galaxies today, including our own Galaxy. The remaining ~30% become satellite galaxies by z~0, and comparison to halo occupation modeling suggests that they are rather faint, with M_r<-19.5. (ABRIDGED)

We discuss the quantization of a class of relativistic fluid models defined in terms of one real and two complex conjugate potentials with values on a K\"{a}hler manifold, and parametrized by the K\"{a}hler potential $K(z,\bar{z})$ and a real number $\lambda$. In the hamiltonian formulation, the canonical conjugate momenta of the potentials are subjected to second class constraints which allow us to apply the symplectic projector method in order to find the physical degrees of freedom and the physical hamiltonian. We construct the quantum theory for that class of models by employing the canonical quantization methods. We also show that a semiclassical theory in which the K\"{a}hler and the complex potential are not quantized has a highly degenerate vacuum. Also, we define and compute the quantum topological number (quantum linking number) operator which has non-vanishing contributions from the K\"{a}hler and complex potentials only. Finally, we show that the vacuum and the states formed by tensoring the number operators eigenstates have zero linking number and show that linear combinations of the tensored number operators eigenstates which have the form of entangled states have non-zero linking number.

We discuss new observational constraints on the abundance of faint high redshift Lyman-alpha emitters secured from a deep Keck near-infrared spectroscopic survey which utilizes the strong magnification provided by lensing galaxy clusters at intermediate redshift (z=0.2-0.5). In each of 9 clusters we have undertaken a systematic `blind' search for line emission with NIRSPEC in the J-band within carefully-selected regions which offer very high magnifications (>10-50x) for background sources with redshifts z~10. The high magnification enables the detection of emission at unprecedented limits (10^41 - 10^42 erg/s), much fainter than those of conventional narrow band imaging and other spectroscopic searches. As the comoving volumes probed are small, our survey is designed to address the important question of whether intrinsically feeble star forming galaxies could provide the dominant ionizing flux at z~10. Our survey has yielded six promising (>5 sigma) candidate Lyman-alpha emitters which lie between z=8.7 and z=10.2. While the predicted reionization photon budget depends upon a large number of physical assumptions, our first glimpse at the z~10 Universe suggests that low luminosity star-forming galaxies contribute a significant proportion of the UV photons necessary for cosmic reionization. [Abridged]

We have observed again two stars inside the ZZ Ceti instability strip that were previously classified as not-observed-to-vary (NOV) by Mukadam et al. (2004) and found them to be low-amplitude variables. Some evidence points to a pure ZZ Ceti instability strip; other evidence contests it. The two stars previously classified as NOV have Sloan Digital Sky Survey (SDSS) spectroscopic effective temperatures that place them inside the ZZ Ceti instability strip, and they were "contaminating" the strip as constant stars, which could indicate that the instability strip was no longer a simple evolutionary stage. A pure instability strip indicates that pulsation is a normal phase which all DAs must go through. We used effective temperatures derived from SDSS optical spectra by comparing them with model atmospheres to look for pulsators through time-resolved photometry and stars previously classified as NOV. Our new results indicate, but do not prove, a pure instability strip, because there are still other NOV stars that need to be observed again. Additionally, we have discovered five other ZZ Ceti stars based on their effective temperatures.

We study the prospects for extracting cosmological and astrophysical parameters from the low radio frequency 21-cm background due to the spin-flip transition of neutral Hydrogen during and prior to the reionization of the Universe. We make use of the angular power spectrum of 21-cm anisotropies, which exists due to inhomogeneities in the neutral Hydrogen density field, the gas temperature field, the gas velocity field, and the spatial distribution of the Lyman-$\alpha$ intensity field associated with first luminous sources that emit UV photons. We extract parameters that describe both the underlying mass power spectrum and the global cosmology, as well as a set of simplified astrophysical parameters that connect fluctuations in the dark matter to those that govern 21-cm fluctuations. We also marginalize over a model for the foregrounds at low radio frequencies. In this general description, we find large degeneracies between cosmological parameters and the astrophysical parameters, though such degeneracies are reduced when strong assumptions are made with respect to the spin temperature relative to the CMB temperature or when complicated sources of anisotropy in the brightness temperature are ignored. Some of the degeneracies between cosmological and astrophysical parameters are broken when 21-cm anisotropy measurements are combined with information from the CMB, such as the temperature and the polarization measurements with Planck. While the overall improvement on the cosmological parameter estimates is not significant when measurements from first-generation interferometers are combined with Planck, such a combination can measure astrophysical parameters such as the ionization fraction in several redshift bins with reasonable accuracy.

Twenty-one centimeter tomography is emerging as a powerful tool to explore the end of the cosmic dark ages and the reionization epoch, but it will only be as good as our ability to accurately model and remove astrophysical foreground contamination. Previous treatments of this problem have focused on the angular structure of the signal and foregrounds and what can be achieved with limited spectral resolution (bandwidths in the 1 MHz range). In this paper we introduce and evaluate a ``blind'' method to extract the multifrequency 21cm signal by taking advantage of the smooth frequency structure of the Galactic and extragalactic foregrounds. We find that 21 cm tomography is typically limited by foregrounds on scales $k\ll 1h/$Mpc and limited by noise on scales $k\gg 1h/$Mpc, provided that the experimental bandwidth can be made substantially smaller than 0.1 MHz. Our results show that this approach is quite promising even for scenarios with rather extreme contamination from point sources and diffuse Galactic emission, which bodes well for upcoming experiments such as LOFAR, MWA, PAST, and SKA.

We have determined the total numbers and specific frequencies of blue, metal-poor globular clusters (GCs) in eight spiral and early-type galaxies. These data, along with five measurements from the literature, show a trend of increasing blue GC specific frequency with increasing mass of the host galaxy. The increase is not accounted for in a simple galaxy formation model in which ellipticals and their GC systems are formed by the merger of typical spiral galaxies. The data appear broadly consistent with hierarchical formation scenarios in which metal-poor GCs are formed over a finite period in the early Universe during the initial stages of galaxy assembly. In this picture, the observed trend is related to biasing, in the sense that the more massive galaxies of today began assembling earlier and therefore formed relatively more GCs during this early epoch of metal-poor GC formation. We discuss how comparisons of the observed specific frequency of metal-poor GCs with model calculations can constrain the formation redshift of these objects.

We study the prospects for extracting detailed statistical properties of the neutral Hydrogen distribution during the era of reionization using the brightness temperature fluctuations from redshifted 21 cm line emission. Detection of this signal is complicated by contamination from foreground sources such as diffuse Galactic synchrotron and free-free emission at low radio frequencies, extragalactic free-free emission from ionized regions and radio point sources. We model these foregrounds to determine the extent to which 21 cm fluctuations can be detected with upcoming experiments. We find that not only the level of correlation from one frequency to another, but also the functional form of the foreground correlations has a substantial impact on foreground removal. We calculate how well the angular power spectra of the 21cm fluctuations can be determined. We also show that the large-scale bias of the neutral Hydrogen gas distribution with respect to the density field, can be determined with high precision and used to distinguish between different reionization histories.

We report the detection of a z~7 galaxy strongly lensed by the massive galaxy cluster Abell 2218 (z=0.175) at 3.6 and 4.5 um using the Spitzer Observatory and at 1.1 um using the Hubble Space Telescope. The new data indicate a refined photometric redshift in the range of 6.6-6.8 depending on the presence of Ly-alpha emission. The spectral energy distribution is consistent with having a significant Balmer break, suggesting that the galaxy is in the poststarburst stage with an age of at least ~50 Myr and quite possibly a few hundred Myr. This suggests the possibility that a mature stellar population is already in place at such a high redshift. Compared with typical Lyman break galaxies at z~3-4, the stellar mass is an order of magnitude smaller (~10^{9} Msun), but the specific star formation rate (star formation rate/M_{star}) is similarly large (> 10^{-9} yr^{-1}), indicating equally vigorous star-forming activity.

The consequences of taking the generalized Chaplygin gas as the dark energy constituent of the Universe on the gravitational waves are studied and the spectrum obtained from this model, for the flat case, is analyzed. Besides its importance for the study of the primordial Universe, the gravitational waves represent an additional perspective (besides the CMB temperature and polarization anisotropies) to evaluate the consistence of the different dark energy models and establish better constraints to their parameters. The analysis presented here takes this fact into consideration to open one more perspective of verification of the generalized Chapligin gas model applicability. Nine particular cases are compared: one where no dark energy is present; two that simulate the $\Lambda$-CDM model; two where the gas acts like the traditional Chaplygin gas; and four where the dark energy is the generalized Chaplygin gas. The different spectra permit to distinguish the $\Lambda$-CDM and the Chaplygin gas scenarios.

We discuss the observational properties of a remarkably faint triply-imaged galaxy revealed in a deep z-band ACS observation of the lensing cluster A2218. A well-constrained mass model for the cluster, which incorporates the outcome of recent Keck spectroscopic campaigns, suggests that the triple system arises via a z>6 source viewed at high magnification (~x25). Optical and IR photometry from HST and Keck confirms the lensing hypothesis and suggests a significant discontinuity occurs in the spectral energy distribution within 9250-9850A. If this break is associated with Gunn-Peterson absorption from neutral hydrogen,6.6<z<7.1 is inferred. Deep Keck spectroscopy conducted using both optical and IR spectrographs fails to reveal any prominent emission lines in this region. However, an IR stellar continuum is detected whose decline below 9800A suggests a spectroscopic redshift of z~7. Regardless of the precise redshift, the source is remarkably compact (<1 kpc) and faint (z=28.0) yet is undergoing vigorous star formation at a rate ~2.6 Mo/yr. An intriguing property is the steep slope of the UV continuum implied by the photometry which may suggest that the source is representative of an early population of galaxies responsible for cosmic reionization. Independent verification of these results is highly desirable but our attempts highlight the difficulty of studying such sources with present facilities and the challenges faced in pushing back the frontiers of the observable universe beyond z~6.5.

We derive the luminosity function of high-redshift Lyman alpha emitting sources from a deep, blind, spectroscopic survey that utilized strong-lensing magnification by intermediate-redshift clusters of galaxies. We observed carefully selected regions near 9 clusters, consistent with magnification factors generally greater than 10 for the redshift range 4.5<z<6.7. Eleven emission-line candidates were located in the range 2.2<z<5.6 whose identification we justify as Lyman alpha, in most cases via further spectroscopic observations. The selection function we constructed for our survey takes into account our varying intrinsic Lyman alpha line sensitivity as a function of wavelength and sky position. By virtue of the strong magnification factor, we provide constraints on the Lyman alpha luminosity function to unprecedented limits of 10^40 erg/s, corresponding to a star-formation rate of 0.01 Msun/yr. Our cumulative z=5 Lyman alpha luminosity function is consistent with a power law form, n(>L) proportional to L^-1 over 10^41 to 10^42.5 erg/s. When combined with the results of other surveys, limited at higher luminosities, our results suggest evidence for the suppression of star formation in low-mass halos, as predicted in popular models of galaxy formation.

We present different methods used to identify high redshift (z>5) objects in the high-magnification regions of lensing galaxy clusters, taking advantage of very well constrained lensing models. The research procedures are explained and discussed. The detection of emission lines in the optical/NIR spectra, such as Lyman-alpha, allows us to determine the redshift of these sources. Thanks to the lensing magnification, it is possible to identify and to study more distant or intrinsically fainter objects with respect to standard field surveys.

Lyman-alpha emission from high-redshift galaxies may be a powerful probe of the ionization history of the IGM at z>6: the observed Lyman-alpha emission line is sensitive to the neutral fraction of IGM hydrogen in the range 0.1-1. We present calculations of observed Lyman-alpha emission lines from z>6 galaxies, illustrating the effect of varying the many free parameters associated with the emitting galaxy, its halo, and the IGM around the galaxy. In particular, we use a dynamic model of the IGM that includes the effect of IGM infall toward the emitting galaxy. Galactic winds may play a crucial role in determining observed Lyman-alpha line fluxes. We compare our model predictions with observations of two z=6.5 galaxies and conclude that, if galactic winds are allowed for, existing observations place no constraint on the neutral fraction of the IGM at z=6.5. Future space-based observations will constrain the importance of galactic winds; if winds are unimportant for the observed z=6.5 galaxies, our models suggest that the IGM neutral fraction at z=6.5 is <~0.1.

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

We present an observational test of the spectrum of the ionizing background at z~5; the test is sufficiently sensitive to determine whether Pop II or Pop III stars are the dominant source of ionizing radiation. The ionizing background at z~5 may reflect the nature of the sources responsible for the final overlap phase of reionization. We find that rest-frame extreme-UV HeI absorption will be detectable in deep spectral observations of some rare z~5 quasars; the ratio of HeI to HI absorption reflects the shape of the ionizing background in the photon energy range between 13.6 and 24.6 eV. Most z~5 quasars have too much HI absorption along their line of sight for HeI absorption to be observed. However, based on current measurements of HI absorber statistics, we use Monte Carlo simulations to demonstrate that the Sloan Digital Sky Survey (SDSS) will discover a sufficient number of z~5 quasars to turn up a quasar suitable for measuring HeI absorption (and we illustrate a selection method to identify that quasar). From simulated observations of a suitable z~5 quasar with a 10-meter telescope, we show that a constraint on the spectral slope of the ionizing background at that redshift can be obtained.

We describe different methods for estimating the bispectrum of Cosmic Microwave Background data. In particular we construct a minimum variance estimator for the flat-sky limit and compare results with previously-studied frequentist methods. Application to the MAXIMA dataset shows consistency with primordial Gaussianity. Weak quadratic non-Gaussianity is characterised by a tunable parameter $f_{NL}$, corresponding to non-Gaussianity at a level $\sim 10^{-5}f_{NL}$ (ratio of non-Gaussian to Gaussian terms), and we find limits of $|f_{NL}|<950$ for the minimum-variance estimator and $|f_{NL}|<1650$ for the usual frequentist estimator. These are the tightest limits on primordial non-Gaussianity which include the full effects of the radiation transfer function.

We calculate the contribution to the cosmic infrared background from very massive metal-free stars at high redshift. We explore two plausible star-formation models and two limiting cases for the reprocessing of the ionizing stellar emission. We find that Population III stars may contribute significantly to the cosmic near-infrared background if the following conditions are met: (i) The first stars were massive, with M > ~100 M_sun. (ii) Molecular hydrogen can cool baryons in low-mass haloes. (iii) Pop III star formation is ongoing, and not shut off through negative feedback effects. (iv) Virialized haloes form stars at about 40 per cent efficiency up to the redshift of reionization, z~7. (v) The escape fraction of the ionizing radiation into the intergalactic medium is small. (vi) Nearly all of the stars end up in massive black holes without contributing to the metal enrichment of the Universe.

We use the measurement of the cosmic microwave background taken during the MAXIMA-1 flight to estimate the bispectrum of cosmological perturbations. We propose an estimator for the bispectrum that is appropriate in the flat sky approximation, apply it to the MAXIMA-1 data and evaluate errors using bootstrap methods. We compare the estimated value with what would be expected if the sky signal were Gaussian and find that it is indeed consistent, with a $\chi^2$ per degree of freedom of approximately unity. This measurement places constraints on models of inflation.

We discuss the physical nature of a remarkably faint pair of Lyman alpha-emitting images discovered close to the giant cD galaxy in the lensing cluster Abell 2218 (z=0.18) during a systematic survey for highly-magnified star-forming galaxies beyond z=5. A well-constrained mass model suggests the pair arises via a gravitationally-lensed source viewed at high magnification. Keck spectroscopy confirms the lensing hypothesis and implies the unlensed source is a very faint (I~30) compact (<150 pc) and isolated object at z=5.576 whose optical emission is substantially contained within the Lyman alpha emission line; no stellar continuum is detectable. The available data suggest the source is a promising candidate for an isolated ~10^6 solar mass system seen producing its first generation of stars close to the epoch of reionization.

We reexamine Wald's no-hair theorem for global anisotropy in the brane world scenarios. We derive a set of sufficient conditions which must be satisfied by the brane matter and bulk metric so that a homogeneous and anisotropic brane asymptotically evolves to a de Sitter spacetime in the presence of a positive cosmological constant on the brane. We discuss the violations of these sufficient conditions and we show that a negative nonlocal energy density or the presence of strong anisotropic stress (i.e., a magnetic field) may lead the brane to collapse. We discuss the generality of these conditions.

Adams et al. (1997b) have noted that according to our current understanding of the unification of fundamental interactions, there should have been phase transitions associated with spontaneous symmetry breaking {\em during} the inflationary era. This may have resulted in the breaking of scale-invariance of the primordial density perturbation for brief periods. A possible such feature was identified in the power spectrum of galaxy clustering in the APM survey at the scale $k \sim 0.1 h$ Mpc^{-1} and it was shown that the secondary acoustic peaks in the power spectrum of the CMB anisotropy should consequently be suppressed. We demonstrate that this prediction is confirmed by the recent Boomerang and Maxima observations, which favour a step-like spectral feature in the range $k \sim (0.06-0.6)h$ Mpc^{-1}, independently of the similar previous indication from the APM data. Such a spectral break enables an excellent fit to both APM and CMB data with a baryon density consistent with the BBN value. It also allows the possibility of a matter-dominated universe with zero cosmological constant, which we show can now account for even the evolution of the abundance of rich clusters.

In the light of the recent Boomerang and Maxima observations of the CMB which show an anomalously low second acoustic peak, we reexamine the prediction by Adams et al (1997) that this would be the consequence of a 'step' in the primordial spectrum induced by a spontaneous symmetry breaking phase transition during primordial inflation. We demonstrate that a deviation from scale-invariance around $k\sim0.1h$~Mpc$^{-1}$ can simultaneously explain both the feature identified earlier in the APM galaxy power spectrum as well the recent CMB anisotropy data, with a baryon density consistent with the BBN value. Such a break also allows a good fit to the data on cluster abundances even for a critical density matter-dominated universe with zero cosmological constant.