Recent analyses of the cosmic microwave background (CMB) and the Lyman-alpha forest indicate a mild preference for a deviation from a power law primordial matter power spectrum (a so-called negative `running'). We use an extension to the BAHAMAS suite of cosmological hydrodynamic simulations to explore the effects that a running scalar spectral index has on large-scale structure (LSS), using Planck CMB constraints to initialize the simulations. We focus on 5 key statistics: i) the non-linear matter power spectrum ii) the halo mass function; iii) the halo two-point auto correlation function; iv) total mass halo density profiles; and v) the halo concentration-mass relation. In terms of the matter power spectrum, we find that a running scalar spectral index affects all k-scales examined in this study, with a negative (positive) running leading to an amplification (suppression) of power. These effects should be easily detectable with upcoming surveys such as LSST and Euclid. In the mass range sampled, a positive running leads to an increase in the mass of galaxy groups and clusters, with the favoured negative running leading to a decrease in mass of lower-mass (M <~ 10^13 M_solar) halos, but an increase for the most massive (M >~ 10^13 M_solar) halos. Changes in the mass are generally confined to 5-10% which, while not insignificant, cannot by itself reconcile the claimed tension between the primary CMB and cluster number counts. We find that running does not significantly affect the shapes of density profiles of matched halos, changing only their amplitude. Finally, we demonstrate that the observed effects on LSS due to a running scalar spectral index are separable from those of baryonic effects to typically a few percent precision.

In this work we consider the impact of spatially-uniform but time-varying dark energy (or `dynamical dark energy', DDE) on large-scale structure in a spatially flat universe, using large cosmological hydrodynamical simulations that form part of the BAHAMAS project. As DDE changes the expansion history of the universe, it impacts the growth of structure. We explore variations in DDE that are constrained to be consistent with the cosmic microwave background. We find that DDE can affect the clustering of matter and haloes at the ~10% level (suppressing it for so-called `freezing' models, while enhancing it for `thawing' models), which should be distinguishable with upcoming large-scale structure surveys. DDE cosmologies can also enhance or suppress the halo mass function (with respect to LCDM) over a wide range of halo masses. The internal properties of haloes are minimally affected by changes in DDE, however. Finally, we show that the impact of baryons and associated feedback processes is largely independent of the change in cosmology and that these processes can be modelled separately to typically better than a few percent accuracy

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

We introduce a simple analytic model of galaxy formation that links the growth of dark matter haloes in a cosmological background to the build-up of stellar mass within them. The model aims to identify the physical processes that drive the galaxy-halo co-evolution through cosmic time. The model restricts the role of baryonic astrophysics to setting the relation between galaxies and their haloes. Using this approach, galaxy properties can be directly predicted from the growth of their host dark matter haloes. We explore models in which the effective star formation efficiency within haloes is a function of mass (or virial temperature) and independent of time. Despite its simplicity, the model reproduces self-consistently the shape and evolution of the cosmic star formation rate density, the specific star formation rate of galaxies, and the galaxy stellar mass function, both at the present time and at high redshifts. By systematically varying the effective star formation efficiency in the model, we explore the emergence of the characteristic shape of the galaxy stellar mass function. The origin of the observed double Schechter function at low redshifts is naturally explained by two efficiency regimes in the stellar to halo mass relation, namely, a stellar feedback regulated stage, and a supermassive black hole regulated stage. By providing a set of analytic differential equations, the model can be easily extended and inverted, allowing the roles and impact of astrophysics and cosmology to be explored and understood.

We investigate the effect of the accelerated expansion of the Universe due to a cosmological constant, $\Lambda$, on the cosmic star formation rate. We utilise hydrodynamical simulations from the EAGLE suite, comparing a $\Lambda$CDM Universe to an Einstein-de Sitter model with $\Lambda=0$. Despite the differences in the rate of growth of structure, we find that dark energy, at its observed value, has negligible impact on star formation in the Universe. We study these effects beyond the present day by allowing the simulations to run forward into the future ($t>13.8$ Gyr). We show that the impact of $\Lambda$ becomes significant only when the Universe has already produced most of its stellar mass, only decreasing the total co-moving density of stars ever formed by ${\approx}15\%$. We develop a simple analytic model for the cosmic star formation rate that captures the suppression due to a cosmological constant. The main reason for the similarity between the models is that feedback from accreting black holes dramatically reduces the cosmic star formation at late times. Interestingly, simulations without feedback from accreting black holes predict an upturn in the cosmic star formation rate for $t>15$ Gyr due to the rejuvenation of massive ($ > 10^{11} \mathrm{M}_{\odot}$) galaxies. We briefly discuss the implication of the weak dependence of the cosmic star formation on $\Lambda$ in the context of the anthropic principle.

Models of the very early universe, including inflationary models, are argued to produce varying universe domains with different values of fundamental constants and cosmic parameters. Using the cosmological hydrodynamical simulation code from the eagle collaboration, we investigate the effect of the cosmological constant on the formation of galaxies and stars. We simulate universes with values of the cosmological constant ranging from Lambda = 0 to Lambda_0 = 300, where Lambda_0 is the value of the cosmological constant in our Universe. Because the global star formation rate in our Universe peaks at t = 3.5 Gyr, before the onset of accelerating expansion, increases in Lambda of even an order of magnitude have only a small effect on the star formation history and efficiency of the universe. We use our simulations to predict the observed value of the cosmological constant, given a measure of the multiverse. Whether the cosmological constant is successfully predicted depends crucially on the measure. The impact of the cosmological constant on the formation of structure in the universe does not seem to be a sharp enough function of Lambda to explain its observed value alone.

We estimate the expected event rate of gravitational wave signals from mergers of supermassive black holes that could be resolved by a space-based interferometer, such as the Evolved Laser Interferometer Space Antenna (eLISA), utilising the reference cosmological hydrodynamical simulation from the EAGLE suite. These simulations assume a $\Lambda$CDM cosmogony with state-of-the-art subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. They have been shown to reproduce the observed galaxy population with unprecedented fidelity. We combine the merger rates of supermassive black holes in EAGLE with the latest phenomenological waveform models to calculate the gravitational waves signals from the intrinsic parameters of the merging black holes. The EAGLE models predict $\sim2$ detections per year by a gravitational wave detector such as eLISA. We find that these signals are largely dominated by mergers between seed mass black holes merging at redshifts between $z\sim2$ and $z\sim1$. In order to investigate the dependence on the assumed black hole seed mass, we introduce an additional model with a black hole seed mass an order of magnitude smaller than in our reference model. We also consider a variation of the reference model where a prescription for the expected delays in the black hole merger timescale has been included after their host galaxies merge. We find that the merger rate is similar in all models, but that the initial black hole seed mass could be distinguished through their detected gravitational waveforms. Hence, the characteristic gravitational wave signals detected by eLISA will provide profound insight into the origin of supermassive black holes and the initial mass distribution of black hole seeds.