Future Square Kilometre Array (SKA) surveys are expected to generate huge datasets of 21cm maps on cosmological scales from the Epoch of Reionization (EoR). We assess the viability of exploiting machine learning techniques, namely, convolutional neural networks (CNN), to simultaneously estimate the astrophysical and cosmological parameters from 21cm maps from semi-numerical simulations. We further convert the simulated 21cm maps into SKA-like mock maps using the detailed SKA antennae distribution, thermal noise and a recipe for foreground cleaning. We successfully design two CNN architectures (VGGNet-like and ResNet-like) that are both efficiently able to extract simultaneously three astrophysical parameters, namely the photon escape fraction (f$_{\rm esc}$), the ionizing emissivity power dependence on halo mass ($C_{\rm ion}$) and the ionizing emissivity redshift evolution index ($D_{\rm ion}$), and three cosmological parameters, namely the matter density parameter ($\Omega_{m}$), the dimensionless Hubble constant ($h$), and the matter fluctuation amplitude ($\sigma_{8}$), from 21cm maps at several redshifts. With the presence of noise from SKA, our designed CNNs are still able to recover these astrophysical and cosmological parameters with great accuracy ($R^{2} > 92\%$), improving to $R^{2} > 99\%$ towards low redshift and low neutral fraction values. Our results show that future 21cm observations can play a key role to break degeneracy between models and tightly constrain the astrophysical and cosmological parameters, using only few frequency channels.

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}.

We examine the properties of damped Lyman-$\alpha$ absorbers (DLAs) emerging from a single set of cosmological initial conditions in two state-of-the-art cosmological hydrodynamic simulations: {\sc Simba} and {\sc Technicolor Dawn}. The former includes star formation and black hole feedback treatments that yield a good match with low-redshift galaxy properties, while the latter uses multi-frequency radiative transfer to model an inhomogeneous ultraviolet background (UVB) self-consistently and is calibrated to match the Thomson scattering optical depth, UVB amplitude, and Ly-$\alpha$ forest mean transmission at $z>5$. Both simulations are in reasonable agreement with the measured stellar mass and star formation rate functions at $z\geq 3$, and both reproduce the observed neutral hydrogen cosmological mass density, $\Omega_{\rm HI}(z)$. However, the DLA abundance and metallicity distribution are sensitive to the galactic outflows' feedback and the UVB amplitude. Adopting a strong UVB and/or slow outflows under-produces the observed DLA abundance, but yields broad agreement with the observed DLA metallicity distribution. By contrast, faster outflows eject metals to larger distances, yielding more metal-rich DLAs whose observational selection may be more sensitive to dust bias. The DLA metallicity distribution in models adopting an ${\rm H}_2$-regulated star formation recipe includes a tail extending to $[M/H] \ll -3$, lower than any DLA observed to date, owing to curtailed star formation in low-metallicity galaxies. Our results show that DLA observations play an imporant role in constraining key physical ingredients in galaxy formation models, complementing traditional ensemble statistics such as the stellar mass and star formation rate functions.

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.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Active Galactic Nuclei (AGN) and star-forming galaxies are leading candidates for being the luminous sources that reionized our Universe. Next-generation 21cm surveys are promising to break degeneracies between a broad range of reionization models, hence revealing the nature of the source population. While many current efforts are focused on a measurement of the 21cm power spectrum, some surveys will also image the 21cm field during reionization. This provides further information with which to determine the nature of reionizing sources. We create a Convolutional Neural Network (CNN) that is efficiently able to distinguish between 21cm maps that are produced by AGN versus galaxies scenarios with an accuracy of 92-100%, depending on redshift and neutral fraction range. An exception to this is when our Universe is highly ionized, since the source models give near-identical 21cm maps in that case. When adding thermal noise from typical 21cm experiments, the classification accuracy depends strongly on the effectiveness of foreground removal. Our results show that if foregrounds can be removed reasonably well, SKA, HERA and LOFAR should be able to discriminate between source models with greater accuracy at a fixed redshift. Only future SKA 21cm surveys are promising to break the degeneracies in the power spectral analysis.

Next-generation 21cm observations will enable imaging of reionization on very large scales. These images will contain more astrophysical and cosmological information than the power spectrum, and hence providing an alternative way to constrain the contribution of different reionizing sources populations to cosmic reionization. Using Convolutional Neural Networks, we present a simple network architecture that is sufficient to discriminate between Galaxy-dominated versus AGN-dominated models, even in the presence of simulated noise from different experiments such as the HERA and SKA.

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.

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.

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.

Bimetric theory describes gravitational interactions in the presence of an extra spin-2 field. Previous work has suggested that its cosmological solutions are generically plagued by instabilities. We show that by taking the Planck mass for the second metric, $M_f$, to be small, these instabilities can be pushed back to unobservably early times. In this limit, the theory approaches general relativity with an effective cosmological constant which is, remarkably, determined by the spin-2 interaction scale. This provides a late-time expansion history which is extremely close to $\Lambda$CDM, but with a technically-natural value for the cosmological constant. We find $M_f$ should be no larger than the electroweak scale in order for cosmological perturbations to be stable by big-bang nucleosynthesis. We further show that in this limit the helicity-0 mode is no longer strongly-coupled at low energy scales.

We consider the semi-classical dynamics of a free massive scalar field in a homogeneous and isotropic cosmological spacetime. The scalar field is quantized using the polymer quantization method assuming that it is described by a gaussian coherent state. For quadratic potentials, the semi-classical equations of motion yield a universe that has an early "polymer inflation" phase which is generic and almost exactly de Sitter, followed by a epoch of slow-roll inflation. We compute polymer corrections to the slow roll formalism, and discuss the probability of inflation in this model using a physical Hamiltonian arising from time gauge fixing. We also show how in this model, it is possible to obtain a significant amount of slow-roll inflation from sub-Planckain initial data, hence circumventing some of the criticisms of standard scenarios. These results show the extent to which a quantum gravity motivated quantization method affects early universe dynamics.

We obtain the general cosmological evolution equations for a classically consistent theory of bimetric gravity. Their analytic solutions are demonstrated to generically allow for a cosmic evolution starting out from a matter dominated FLRW universe while relaxing towards a de Sitter (anti-de Sitter) phase at late cosmic time. In particular, we examine a subclass of models which contain solutions that are able to reproduce the expansion history of the cosmic concordance model inspite of the nonlinear couplings of the two metrics. This is demonstrated explicitly by fitting these models to observational data from Type Ia supernovae, Cosmic Microwave Background and Baryon Acoustic Oscillations.

The young open star cluster M25 (IC 4725) is located in the direction of the galactic center in a crowded region, near much irregular absorption features on Sagittarius arm. This cluster has some difficult observing problems due to its southern location. The mass data available in the literature have been gathered to reinvestigate this cluster using most photometric tools to determine its main photometric parameters. This system is found to be at a distance of 600 pc from the sun, of -52.82 pc from the galactic plane; and the median age is 9.45 x 10^(7) yr. More than 220 stars with mean reddening of 0.50 mag and absorption of 1.62 mag are found within the cluster.

There are good indications that fundamental physics gives rise to a modified space-momentum uncertainty relation that implies the existence of a minimum length scale. We implement this idea in the scalar field theory that describes density perturbations in flat Robertson-Walker space-time. This leads to a non-linear time-dependent dispersion relation that encodes the effects of Planck scale physics in the inflationary epoch. Unruh type dispersion relations naturally emerge in this approach, while unbounded ones are excluded by the minimum length principle. We also find red-shift induced modifications of the field theory, due to the reduction of degrees of freedom at high energies, that tend to dampen the fluctuations at trans-Planckian momenta. In the specific example considered, this feature helps determine the initial state of the fluctuations, leading to a flat power spectrum.

We report on CCD optical (B and V passbands) and near IR (J and K bands) observations in the region of the old open cluster NGC 2141. By combining the two sets of photometry (500 stars in common) we derive new estimates of the cluster fundamental parameters. We confirm that the cluster is 2.5 Gyrs old, but, with respect to previous investigations, we obtain a slightly larger reddening (E(B-V)=0.40), and a slightly shorter distance (3.8 kpc) from the Sun. Finally we present Luminosity Function (LF) in the V band, which is another age indicator. We provide a good fit for the age range inferred from isochrones by assuming the Kroupa et al. (1993) IMF up to $M_V$=5.0. We interpret the disagreemt at fainter magnitudes as an evidence of mass segregation.