The Square Kilometer Array is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z>6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centered on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z=7-10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at redshifts z=10 or below. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z~10, due to the very small number of line-of-sight interlopers expected at similar redshifts.

The Square Kilometre Array (SKA) is a planned large radio interferometer designed to operate over a wide range of frequencies, and with an order of magnitude greater sensitivity and survey speed than any current radio telescope. The SKA will address many important topics in astronomy, ranging from planet formation to distant galaxies. However, in this work, we consider the perspective of the SKA as a facility for studying physics. We review four areas in which the SKA is expected to make major contributions to our understanding of fundamental physics: cosmic dawn and reionisation; gravity and gravitational radiation; cosmology and dark energy; and dark matter and astroparticle physics. These discussions demonstrate that the SKA will be a spectacular physics machine, which will provide many new breakthroughs and novel insights on matter, energy and spacetime.

Several investigations of the X-ray variability of active galactic nuclei (AGN) using the normalised excess variance (${\sigma^2_{\rm NXS}}$) parameter have shown that variability has a strong anti-correlation with black hole mass ($M_{\rm BH}$) and X-ray luminosity ($L_{\rm X}$). In this study we confirm these previous correlations and find no evidence of a redshift evolution. Using observations from XMM-Newton, we determine the ${\sigma^2_{\rm NXS}}$ and $L_{\rm X}$ for a sample of 1091 AGN drawn from the XMM-Newton Cluster Survey (XCS) - making this the largest study of X-ray spectral properties of AGNs. We created light-curves in three time-scales; 10 ks, 20 ks and 40 ks and used these to derive scaling relations between ${\sigma^2_{\rm NXS}}$, $L_{\rm X}$ (2.0-10 keV range) and literature estimates of $M_{\rm BH}$ from reverberation mapping. We confirm the anti-correlation between $M_{\rm BH}$ and ${\sigma^2_{\rm NXS}}$ and find a positive correlation between $M_{\rm BH}$ and $L_{\rm X}$. The use of ${\sigma^2_{\rm NXS}}$ is practical only for pointed observations where the observation time is tens of kiloseconds. For much shorter observations one cannot accurately quantify variability to estimate $M_{\rm BH}$. Here we describe a method to derive $L_{\rm X}$ from short duration observations and used these results as an estimate for $M_{\rm BH}$. We find that it is possible to estimate $L_{\rm X}$ from observations of just a few hundred seconds and that when correlated with $M_{\rm BH}$, the relation is statistically similar to the relation of $M_{\rm BH}$-$L_{\rm X}$ derived from a spectroscopic analysis of full XMM observations. This method may be particularly useful to the eROSITA mission, an all-sky survey, which will detect $>$10$^{6}$ AGN.

The 4MOST Cosmology Redshift Survey (CRS) will perform stringent cosmological tests via spectroscopic clustering measurements that will complement the best lensing, cosmic microwave background and other surveys in the southern hemisphere. The combination of carefully selected samples of bright galaxies, luminous red galaxies, emission-line galaxies and quasars, totalling about 8 million objects over the redshift range $z = 0.15$ to $3.5$, will allow definitive tests of gravitational physics. Many key science questions will be addressed by combining CRS spectra of these targets with data from current or future facilities such as the Large Synoptic Survey Telescope, the Square Kilometre Array and the Euclid mission.

Groups and clusters of galaxies are a current focus of astronomical research owing to their role in determining the environmental effects on galaxies and the constraints they provide to cosmology. The eROSITA X-ray telescope on board the Spectrum Roentgen Gamma observatory will be launched in 2019 and will have completed eight scans of the full sky when 4MOST starts operating. The experiment will detect groups and clusters of galaxies through X-ray emission from the hot intergalactic medium. The purpose of the 4MOST eROSITA Galaxy Cluster Redshift Survey is to provide spectroscopic redshifts of the optical counterparts to the X-ray emission from 40,000 groups and clusters of galaxies so as to perform dynamical estimates of the total mass and to measure the properties of the member galaxies. The survey aims to obtain precise redshift measurements of the photometrically identified brightest cluster galaxies at redshift $z > 0.7$. At lower redshifts ($z < 0.7$), the programme aims to sample over 15 member galaxies per cluster and enable dynamical mass measurements to calibrate the clusters for cosmological experiments. At $z < 0.2$, eROSITA will also detect X-ray emission from galaxy groups and filaments. 4MOST spectroscopic data from the survey will be used for optical identification of galaxy groups down to eROSITA's mass detection limits of $10^{13} M_\odot$, as well as the detection of the largest filaments for pioneering studies of their X-ray emission.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Future large-scale spectroscopic astronomical surveys, e.g. Euclid, will enable the compilation of vast new catalogues of clusters and voids in the galaxy distribution. By combining the constraining power of both cluster and void number counts, such surveys could place stringent simultaneous limits on the sum of neutrino masses $M_\nu$ and the dark energy equation of state $w(z) = w_0 + w_a z/(1+z)$. For minimal normal-hierarchy neutrino masses, we forecast that Euclid clusters + voids ideally could reach uncertainties $\sigma(M_\nu) \lesssim 15$ meV, $\sigma(w_0) \lesssim~0.02$, $\sigma(w_a) \lesssim 0.07$, independent of other data. Such precision is competitive with expectations for e.g. galaxy clustering and weak lensing in future cosmological surveys, and could reject an inverted neutrino mass hierarchy at $\gtrsim 99\%$ confidence.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

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

Modern scientific cosmology pushes the boundaries of knowledge and the knowable. This is prompting questions on the nature of scientific knowledge. A central issue is what defines a 'good' model. When addressing global properties of the Universe or its initial state this becomes a particularly pressing issue. How to assess the probability of the Universe as a whole is empirically ambiguous, since we can examine only part of a single realisation of the system under investigation: at some point, data will run out. We review the basics of applying Bayesian statistical explanation to the Universe as a whole. We argue that a conventional Bayesian approach to model inference generally fails in such circumstances, and cannot resolve, e.g., the so-called 'measure problem' in inflationary cosmology. Implicit and non-empirical valuations inevitably enter model assessment in these cases. This undermines the possibility to perform Bayesian model comparison. One must therefore either stay silent, or pursue a more general form of systematic and rational model assessment. We outline a generalised axiological Bayesian model inference framework, based on mathematical lattices. This extends inference based on empirical data (evidence) to additionally consider the properties of model structure (elegance) and model possibility space (beneficence). We propose this as a natural and theoretically well-motivated framework for introducing an explicit, rational approach to theoretical model prejudice and inference beyond data.

We characterize the X-ray luminosity--temperature ($L_{\rm X}-T$) relation using a sample of 353 clusters and groups of galaxies with temperatures in excess of 1 keV, spanning the redshift range $0.1 < z < 0.6$, the largest ever assembled for this purpose. All systems are part of the ${\it XMM-Newton}$ Cluster Survey (XCS), and have also been independently identified in Sloan Digital Sky Survey (SDSS) data using the redMaPPer algorithm. We allow for redshift evolution of the normalisation and intrinsic scatter of the $L_{\rm X}-T$ relation, as well as, for the first time, the possibility of a temperature-dependent change-point in the exponent of such relation. However, we do not find strong statistical support for deviations from the usual modelling of the $L_{\rm X}-T$ relation as a single power-law, where the normalisation evolves self-similarly and the scatter remains constant with time. Nevertheless, assuming {\it a priori} the existence of the type of deviations considered, then faster evolution than the self-similar expectation for the normalisation of the $L_{\rm X}-T$ relation is favoured, as well as a decrease with redshift in the scatter about the $L_{\rm X}-T$ relation. Further, the preferred location for a change-point is then close to 2 keV, possibly marking the transition between the group and cluster regimes. Our results also indicate an increase in the power-law exponent of the $L_{\rm X}-T$ relation when moving from the group to the cluster regime, and faster evolution in the former with respect to the later, driving the temperature-dependent change-point towards higher values with redshift.

Combining galaxy cluster and void abundances is a novel, powerful way to constrain deviations from General Relativity and the $\Lambda$CDM model. For a flat $w$CDM model with growth of large-scale structure parameterized by the redshift-dependent growth index $\gamma(z) = \gamma_0 + \gamma_1 z/(1+z)$ of linear matter perturbations, combining void and cluster abundances in future surveys with Euclid and the 4-metre Multi-Object Spectroscopic Telescope (4MOST) could improve the Figure of Merit for ($w, \gamma_0, \gamma_1$) by a factor of 20 compared to individual abundances. In an ideal case, improvement on current cosmological data is a Figure of Merit factor $600$ or more.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

We measure the evolution of the velocity dispersion--temperature ($\sigma_{\rm v}$--$T_{\rm X}$) relation up to $z = 1$ using a sample of 38 galaxy clusters drawn from the \textit{XMM} Cluster Survey. This work improves upon previous studies by the use of a homogeneous cluster sample and in terms of the number of high redshift clusters included. We present here new redshift and velocity dispersion measurements for 12 $z > 0.5$ clusters observed with the GMOS instruments on the Gemini telescopes. Using an orthogonal regression method, we find that the slope of the relation is steeper than that expected if clusters were self-similar, and that the evolution of the normalisation is slightly negative, but not significantly different from zero ($\sigma_{\rm v} \propto T^{0.86 \pm 0.14} E(z)^{-0.37 \pm 0.33}$). We verify our results by applying our methods to cosmological hydrodynamical simulations. The lack of evolution seen in our data is consistent with simulations that include both feedback and radiative cooling.

We describe updates to the \redmapper{} algorithm, a photometric red-sequence cluster finder specifically designed for large photometric surveys. The updated algorithm is applied to $150\,\mathrm{deg}^2$ of Science Verification (SV) data from the Dark Energy Survey (DES), and to the Sloan Digital Sky Survey (SDSS) DR8 photometric data set. The DES SV catalog is locally volume limited, and contains 786 clusters with richness $\lambda>20$ (roughly equivalent to $M_{\rm{500c}}\gtrsim10^{14}\,h_{70}^{-1}\,M_{\odot}$) and $0.2<z<0.9$. The DR8 catalog consists of 26311 clusters with $0.08<z<0.6$, with a sharply increasing richness threshold as a function of redshift for $z\gtrsim 0.35$. The photometric redshift performance of both catalogs is shown to be excellent, with photometric redshift uncertainties controlled at the $\sigma_z/(1+z)\sim 0.01$ level for $z\lesssim0.7$, rising to $\sim0.02$ at $z\sim0.9$ in DES SV. We make use of \emph{Chandra} and \emph{XMM} X-ray and South Pole Telescope Sunyaev-Zeldovich data to show that the centering performance and mass--richness scatter are consistent with expectations based on prior runs of \redmapper{} on SDSS data. We also show how the \redmapper{} \photoz{} and richness estimates are relatively insensitive to imperfect star/galaxy separation and small-scale star masks.

Combining galaxy cluster and void abundances breaks the degeneracy between mean matter density $\Omega_{\rm m}$ and power spectrum normalization $\sigma_8$. In a first for voids, we constrain $\Omega_{\rm m} = 0.21 \pm 0.10$ and $\sigma_8 = 0.95 \pm 0.21$ for a flat $\Lambda$CDM universe, using extreme-value statistics on the claimed largest cluster and void. The Planck-consistent results detect dark energy with two objects, independently of other dark energy probes. Cluster-void studies also offer complementarity in scale, density, and non-linearity - of particular interest for testing modified-gravity models.

In this work we model the quintessence potential in a Taylor series expansion, up to second order, around the present-day value of the scalar field. The field is evolved in a thawing regime assuming zero initial velocity. We use the latest data from the Planck satellite, baryonic acoustic oscillations observations from the Sloan Digital Sky Survey, and Supernovae luminosity distance information from Union2.1 to constrain our models parameters, and also include perturbation growth data from the WiggleZ, BOSS and the 6dF surveys. The supernova data provide the strongest individual constraint on the potential parameters. We show that the growth data performance is competitive with the other datasets in constraining the dark energy parameters we introduce. We also conclude that the combined constraints we obtain for our model parameters, when compared to previous works of nearly a decade ago, have shown only modest improvement, even with new growth of structure data added to previously-existent types of data.

Using the science verification data of the Dark Energy Survey (DES) for a new sample of 106 X-Ray selected clusters and groups, we study the stellar mass growth of Bright Central Galaxies (BCGs) since redshift 1.2. Compared with the expectation in a semi-analytical model applied to the Millennium Simulation, the observed BCGs become under-massive/under-luminous with decreasing redshift. We incorporate the uncertainties associated with cluster mass, redshift, and BCG stellar mass measurements into analysis of a redshift-dependent BCG-cluster mass relation, $m_{*}\propto(\frac{M_{200}}{1.5\times 10^{14}M_{\odot}})^{0.24\pm 0.08}(1+z)^{-0.19\pm0.34}$, and compare the observed relation to the model prediction. We estimate the average growth rate since $z = 1.0$ for BCGs hosted by clusters of $M_{200, z}=10^{13.8}M_{\odot}$, at $z=1.0$: $m_{*, BCG}$ appears to have grown by $0.13\pm0.11$ dex, in tension at $\sim 2.5 \sigma$ significance level with the $0.40$ dex growth rate expected from the semi-analytic model. We show that the buildup of extended intra-cluster light after $z=1.0$ may alleviate this tension in BCG growth rates.

Modified gravity has garnered interest as a backstop against dark matter and dark energy (DE). As one possible modification, the graviton can become massive, which introduces a new scalar field - here with a Galileon-type symmetry. The field can lead to a nontrivial equation of state (EOS) of DE which is density-and-scale-dependent. Tension between Type Ia supernovae and Planck could be reduced. In voids the scalar field dramatically alters the EOS of DE, induces a soon-observable gravitational slip between the two metric potentials, and develops a topological defect (domain wall) due to a nontrivial vacuum structure for the field.

This paper presents both the result of a search for fossil systems (FSs) within the XMM Cluster Survey and the Sloan Digital Sky Survey and the results of a study of the stellar mass assembly and stellar populations of their fossil galaxies. In total, 17 groups and clusters are identified at z < 0.25 with large magnitude gaps between the first and fourth brightest galaxies. All the information necessary to classify these systems as fossils is provided. For both groups and clusters, the total and fractional luminosity of the brightest galaxy is positively correlated with the magnitude gap. The brightest galaxies in FSs (called fossil galaxies) have stellar populations and star formation histories which are similar to normal brightest cluster galaxies (BCGs). However, at fixed group/cluster mass, the stellar masses of the fossil galaxies are larger compared to normal BCGs, a fact that holds true over a wide range of group/cluster masses. Moreover, the fossil galaxies are found to contain a significant fraction of the total optical luminosity of the group/cluster within 0.5R200, as much as 85%, compared to the non-fossils, which can have as little as 10%. Our results suggest that FSs formed early and in the highest density regions of the universe and that fossil galaxies represent the end products of galaxy mergers in groups and clusters. The online FS catalog can be found at http://www.astro.ljmu.ac.uk/~xcs/Harrison2012/XCSFSCat.html.

We present a list of 15 clusters of galaxies, serendipitously detected by the XMM Cluster Survey (XCS), that have a high probability of detection by the Planck satellite. Three of them already appear in the Planck Early Sunyaev-Zel'dovich (ESZ) catalogue. The estimation of the Planck detection probability assumes the flat Lambda cold dark matter (LambdaCDM) cosmology most compatible with 7-year Wilkinson Microwave Anisotropy Probe (WMAP7) data. It takes into account the XCS selection function and Planck sensitivity, as well as the covariance of the cluster X-ray luminosity, temperature, and integrated comptonization parameter, as a function of cluster mass and redshift, determined by the Millennium Gas Simulations. We also characterize the properties of the galaxy clusters in the final data release of the XCS that we expect Planck will have detected by the end of its extended mission. Finally, we briefly discuss possible joint applications of the XCS and Planck data.

We use Chandra X-ray and Spitzer infrared observations to explore the AGN and starburst populations of XMMXCS J2215.9-1738 at z=1.46, one of the most distant spectroscopically confirmed galaxy clusters known. The high resolution X-ray imaging reveals that the cluster emission is contaminated by point sources that were not resolved in XMM observations of the system, and have the effect of hardening the spectrum, leading to the previously reported temperature for this system being overestimated. From a joint spectroscopic analysis of the Chandra and XMM data, the cluster is found to have temperature T=4.1_-0.9^+0.6 keV and luminosity L_X=(2.92_-0.35^+0.24)x10^44 erg/s extrapolated to a radius of 2 Mpc. As a result of this revised analysis, the cluster is found to lie on the sigma_v-T relation, but the cluster remains less luminous than would be expected from self-similar evolution of the local L_X-T relation. Two of the newly discovered X-ray AGN are cluster members, while a third object, which is also a prominent 24 micron source, is found to have properties consistent with it being a high redshift, highly obscured object in the background. We find a total of eight >5 sigma 24 micron sources associated with cluster members (four spectroscopically confirmed, and four selected using photometric redshifts), and one additional 24 micron source with two possible optical/near-IR counterparts that may be associated with the cluster. Examining the IRAC colors of these sources, we find one object is likely to be an AGN. Assuming that the other 24 micron sources are powered by star formation, their infrared luminosities imply star formation rates ~100 M_sun/yr. We find that three of these sources are located at projected distances of <250 kpc from the cluster center, suggesting that a large amount of star formation may be taking place in the cluster core, in contrast to clusters at low redshift.

We forecast the constraints on the values of sigma_8, Omega_m, and cluster scaling relation parameters which we expect to obtain from the XMM Cluster Survey (XCS). We assume a flat Lambda-CDM Universe and perform a Monte Carlo Markov Chain analysis of the evolution of the number density of galaxy clusters that takes into account a detailed simulated selection function. Comparing our current observed number of clusters shows good agreement with predictions. We determine the expected degradation of the constraints as a result of self-calibrating the luminosity-temperature relation (with scatter), including temperature measurement errors, and relying on photometric methods for the estimation of galaxy cluster redshifts. We examine the effects of systematic errors in scaling relation and measurement error assumptions. Using only (T,z) self-calibration, we expect to measure Omega_m to +-0.03 (and Omega_Lambda to the same accuracy assuming flatness), and sigma_8 to +-0.05, also constraining the normalization and slope of the luminosity-temperature relation to +-6 and +-13 per cent (at 1sigma) respectively in the process. Self-calibration fails to jointly constrain the scatter and redshift evolution of the luminosity-temperature relation significantly. Additional archival and/or follow-up data will improve on this. We do not expect measurement errors or imperfect knowledge of their distribution to degrade constraints significantly. Scaling-relation systematics can easily lead to cosmological constraints 2sigma or more away from the fiducial model. Our treatment is the first exact treatment to this level of detail, and introduces a new `smoothed ML' estimate of expected constraints.

We update and extend our previous work reconstructing the potential of a quintessence field from current observational data. We extend the cosmological dataset to include new supernova data, plus information from the cosmic microwave background and from baryon acoustic oscillations. We extend the modelling by considering Pad\'e approximant expansions as well as Taylor series, and by using observations to assess the viability of the tracker hypothesis. We find that parameter constraints have improved by a factor of two, with a strengthening of the preference of the cosmological constant over evolving quintessence models. Present data show some signs, though inconclusive, of favouring tracker models over non-tracker models under our assumptions.

We describe an algorithm which directly determines the quintessence potential from observational data, without using an equation of state parametrisation. The strategy is to numerically determine observational quantities as a function of the expansion coefficients of the quintessence potential, which are then constrained using a likelihood approach. We further impose a model selection criterion, the Bayesian Information Criterion, to determine the appropriate level of the potential expansion. In addition to the potential parameters, the present-day quintessence field velocity is kept as a free parameter. Our investigation contains unusual model types, including a scalar field moving on a flat potential, or in an uphill direction, and is general enough to permit oscillating quintessence field models. We apply our method to the `gold' Type Ia supernovae sample of Riess et al. (2004), confirming the pure cosmological constant model as the best description of current supernovae luminosity-redshift data. Our method is optimal for extracting quintessence parameters from future data.