Next generation telescopes such as the James Webb Space Telescope (JWST) and the Wide Field Infrared Survey Telescope (WFIRST) will enable us to study the first billion years of our Universe in unprecedented detail. In this work we use the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework, that couples galaxy formation and reionization (for a wide range of reionization feedback models), to estimate the cosmic variance expected in the UV Luminosity Function (UV LF) and the Stellar Mass Function (SMF) in JWST surveys. We find that different reionization scenarios play a minor role in the cosmic variance. Most of the cosmic variance is completely driven by the underlying density field and increases above 100% for M$_{UV} \sim$ -17.5 (-20) at z = 12 (6) for the JADES-deep survey (the deep JWST Advanced Extragalactic Survey with an area of 46 arcmin$^2$); the cosmic variance decreases with an increasing survey area roughly independently of redshift. We find that the faint-end slope of the Lyman Break Galaxies (LBGs) UV LF becomes increasingly shallower with increasing reionization feedback and show how JWST observations will be able to distinguish between different models of reionization feedback at z > 9, even accounting for cosmic variance. We also show the environments (in terms of density and ionization fields) of Lyman Break Galaxies during the EoR. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas.

We introduce a new self-consistent model of galaxy evolution and reionization, ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS), which couples a state-of-the-art N-body simulation with the semi-analytical galaxy evolution DELPHI and the semi-numerical reionization scheme CIFOG. ASTRAEUS includes all the key processes of galaxy formation and evolution (including accretion, mergers, supernova and radiative feedback) and follows the time and spatial evolution of the ionized regions in the intergalactic medium (IGM). Importantly, it explores different radiative feedback models that cover the physically plausible parameter space, ranging from a weak and delayed to a strong and immediate reduction of gas mass available for star formation. From our simulation suite that covers the different radiative feedback prescriptions and ionization topologies, we find that radiative feedback continuously reduces star formation in galaxies with $M_h<10^{9.5}M_{\odot}$ upon local reionization; larger mass halos are unaffected even for the strongest and immediate radiative feedback cases during reionization. For this reason, the ionization topologies of different radiative feedback scenarios differ only on scales smaller than $1-2$Mpc, and significant deviations are only found when physical parameters (e.g. the escape fraction of ionizing photons) are altered based on galactic properties. Finally, we find observables (the ultra-violet luminosity function, stellar mass function, reionization histories and ionization topologies) are hardly affected by the choice of the used stellar population synthesis models that either model single stars or binaries.

This white paper highlights the crucial and urgent synergies required between WFIRST, Subaru Hyper Suprime-Cam or other >25m-class telescopes galaxy observations and SKA 21cm measurements to constrain the nature of reionization (ionization history and topology) and its sources.

The 21-cm power spectrum (PS) has been shown to be a powerful discriminant of reionization and cosmic dawn astrophysical parameters. However, the 21-cm tomographic signal is highly non-Gaussian. Therefore there is additional information which is wasted if only the PS is used for parameter recovery. Here we showcase astrophysical parameter recovery directly from 21-cm images, using deep learning with convolutional neural networks (CNN). Using a database of 2D images taken from 10,000 21-cm lightcones (each generated from different cosmological initial conditions), we show that a CNN is able to recover parameters describing the first galaxies: (i) Tvir , their minimum host halo virial temperatures (or masses) capable of hosting efficient star formation; (ii) {\zeta} , their typical ionizing efficiencies; (iii) LX/SFR , their typical soft-band X-ray luminosity to star formation rate; and (iv) E0 , the minimum X-ray energy capable of escaping the galaxy into the IGM. For most of their allowed ranges, log Tvir and log LX/SFR are recovered with < 1% uncertainty, while {\zeta} and E0 are recovered within 10% uncertainty. Our results are roughly comparable to the accuracy obtained from Monte Carlo Markov Chain sampling of the PS with 21CMMC for the two mock observations analyzed previously, although we caution that we do not yet include noise and foreground contaminants in this proof-of-concept study.

Dwarf galaxies are ideal laboratories to study the physics of the interstellar medium (ISM). Emission lines have been widely used to this aim. Retrieving the full information encoded in the spectra is therefore essential. This can be efficiently and reliably done using Machine Learning (ML) algorithms. Here, we apply the ML code GAME to MUSE (Multi Unit Spectroscopic Explorer) and PMAS (Potsdam Multi Aperture Spectrophotometer) Integral Field Unit (IFU) observations of two nearby Blue Compact Galaxies (BCGs): Henize 2-10 and IZw18. We derive spatially resolved maps of several key ISM physical properties. We find that both galaxies show a remarkably uniform metallicity distribution. Henize 2-10 is a star forming dominated galaxy, with a Star Formation Rate (SFR) of about 1.2 M$_{\odot}$ yr$^{-1}$. Henize 2-10 features dense and dusty ($A_V$ up to 5-7 mag) star forming central sites. We find IZw18 to be very metal-poor ($Z$= 1/20 $Z_{\odot}$). IZw18 has a strong interstellar radiation field, with a large ionization parameter. We also use models of PopIII stars spectral energy distribution as a possible ionizing source for the HeII $\lambda$4686 emission detected in the IZw18 NW component. We find that PopIII stars could provide a significant contribution to the line intensity. The upper limit to the PopIII star formation is 52% of the total IZw18 SFR.

We present an updated, optimized version of GAME (GAlaxy Machine learning for Emission lines), a code designed to infer key interstellar medium physical properties from emission line intensities of UV/optical/far infrared galaxy spectra. The improvements concern: (a) an enlarged spectral library including Pop III stars; (b) the inclusion of spectral noise in the training procedure, and (c) an accurate evaluation of uncertainties. We extensively validate the optimized code and compare its performance against empirical methods and other available emission line codes (PYQZ and HII-CHI-MISTRY) on a sample of 62 SDSS stacked galaxy spectra and 75 observed HII regions. Very good agreement is found for metallicity. However, ionization parameters derived by GAME tend to be higher. We show that this is due to the use of too limited libraries in the other codes. The main advantages of GAME are the simultaneous use of all the measured spectral lines, and the extremely short computational times. We finally discuss the code potential and limitations.

We present a new approach based on Supervised Machine Learning (SML) algorithms to infer key physical properties of galaxies (density, metallicity, column density and ionization parameter) from their emission line spectra. We introduce a numerical code (called GAME, GAlaxy Machine learning for Emission lines) implementing this method and test it extensively. GAME delivers excellent predictive performances, especially for estimates of metallicity and column densities. We compare GAME with the most widely used diagnostics (e.g. R$_{23}$, [NII]$\lambda$6584 / H$\alpha$ indicators) showing that it provides much better accuracy and wider applicability range. GAME is particularly suitable for use in combination with Integral Field Unit (IFU) spectroscopy, both for rest-frame optical/UV nebular lines and far-infrared/sub-mm lines arising from Photo-Dissociation Regions. Finally, GAME can also be applied to the analysis of synthetic galaxy maps built from numerical simulations.

High Energy Neutrino Astronomy has been revealed by a sudden change in the flavor composition around maximal energies since three years of recording in ice km detector. However these discover didn' t led to any clear promised Neutrino Astronomy land yet. No correlation with hard gamma sources was found. Moreover the astrophysical spectra expected at Fermi value, seem to converge to a softer value, also required to avoid any Glashow resonant neutrino signal. We suggest a main solution within a composite flux ruled by prompt atmospheric neutrinos. Nevertheless the recent discover of twentyone crossing muons at hundreds TeVs had shown a first narrow overlapped doublet and a correlated track with a peculiar hard UHECR event clustering pointing to a well known microjet in bynary precession; these first connections hint for a non negligible astrophysical component making neutrino astronomy anyway already more than a hope.

In this paper we review all the up-to-date Ultra High Energy Cosmic Rays (UHECR) events reported by AUGER, by Telescope Array (TA) and by AGASA in common coordinate maps. We also confirm our earliest (2008-2013) model, where UHECR are mostly made by light nuclei (namely He, Be, B), which explains the Virgo absence and confirms M82 as the main source for North TA Hot Spot. Many more sources, such as NGC253 and several galactic ones, are possible candidates for most of the 376 UHECR events. Several correlated map, already considered in recent years, are then reported to show all the events, with their statistical correlation values.

The variabile star Delta Scorpii is in conjunction with the Sun at the end of November each year. We studied its magnitude by averaging the observations of 28 Nov - 1 Dec from 1996 to 2013 using the coronograph LASCO C3 on-board the SOHO Satellite and we extended of four years, i.e. 25 % of the total light curve, back to 1996, with respect to the present AAVSO dataset on this star. The 0.2 magnitude scatters of the single measurements have been studied and the sources of such disturbances are vignetting and diffraction patterns from the coronograph. The new data collected on Delta Scorpii show its minimum at mv=2.5 magnitudes for 1996 and 1997, confirming the values observed during the minimum of 2009, and the main periodicity of 11 years in the stellar variability.