On the scale of galactic haloes, the distribution of matter in the cosmos is affected by energetic, non-gravitational processes; so-called baryonic feedback. A lack of knowledge about the details of how feedback processes redistribute matter is a source of uncertainty for weak-lensing surveys, which accurately probe the clustering of matter in the Universe over a wide range of scales. We develop a cosmology-dependent model for the matter distribution that simultaneously accounts for the clustering of dark matter, gas and stars. We inform our model by comparing it to power spectra measured from the BAHAMAS suite of hydrodynamical simulations. As well as considering matter power spectra, we also consider spectra involving the electron-pressure field, which directly relates to the thermal Sunyaev-Zel'dovich (tSZ) effect. We fit parameters in our model so that it can simultaneously model both matter and pressure data and such that the distribution of gas as inferred from tSZ has influence on the matter spectrum predicted by our model. We present two variants; one that matches the feedback-induced suppression seen in the matter-matter power spectrum at the per-cent level and a second that matches the matter-matter data slightly less well (~2 per cent), but that is able to simultaneously model the matter-electron pressure spectrum at the ~15 per-cent level. We envisage our models being used to simultaneously learn about cosmological parameters and the strength of baryonic feedback using a combination of tSZ and lensing auto- and cross-correlation data.

We evaluate the ability of Convolutional Neural Networks (CNNs) to predict galaxy cluster masses in the BAHAMAS hydrodynamical simulations. We train four separate single-channel networks using: stellar mass, soft X-ray flux, bolometric X-ray flux, and the Compton $y$ parameter as observational tracers, respectively. Our training set consists of $\sim$6400 synthetic cluster images generated from the simulation, while an additional $\sim$1600 images form a test set. We also train a `multi-channel' CNN by combining the four observational tracers. We utilize $\texttt{Keras}$ with a $\texttt{Tensorflow}$ backend to train the network, and all four converge within 2000 epochs. The cluster masses predicted from these networks are evaluated using the average fractional difference between predicted cluster mass and true cluster mass. The resulting predictions are especially precise for halo masses in the range $10^{13.25}M_{\odot}<M<10^{14.5}M_{\odot}$, where all five networks produce mean mass biases of order $\approx 1\%$ with a scatter on the mean bias of $\approx 0.5\%$. The network trained with Compton $y$ parameter maps yields the most precise predictions. We interpret the network's behaviour using two diagnostic tests to determine which features are used to predict cluster mass. The CNN trained with stellar mass images detect galaxies (not surprisingly), while CNNs trained with gas-based tracers utilise the shape of the signal to estimate cluster mass.

The location of a galaxy cluster's centroid is typically derived from observations of the galactic and/or gas component of the cluster, but these typically deviate from the true centre. This can produce bias when observations are combined to study average cluster properties. Using data from the BAHAMAS cosmological hydrodynamic simulations we study this bias in both two and three dimensions for 2000 clusters over the $10^{13} - 10^{15} ~\mathrm{M_{\odot}}$ mass range. We quantify and model the offset distributions between observationally-motivated centres and the `true' centre of the cluster, which is taken to be the most gravitationally bound particle measured in the simulation. We fit the cumulative distribution function of offsets with an exponential distribution and a Gamma distribution fit well with most of the centroid definitions. The galaxy-based centres can be seen to be divided into a mis-centred group and a well-centred group, with the well-centred group making up about $60\%$ of all the clusters. Gas-based centres are overall less scattered than galaxy-based centres. We also find a cluster-mass dependence of the offset distribution of gas-based centres, with generally larger offsets for smaller mass clusters. We then measure cluster density profiles centred at each choice of the centres and fit them with empirical models. Stacked, mis-centred density profiles fit to the Navarro-Frenk-White dark-matter profile and Komatsu-Seljak gas profile show that recovered shape and size parameters can significantly deviate from the true values. For the galaxy-based centres, this can lead to cluster masses being underestimated by up to $10\%$.

We study the new mechanism of the axion production suggested recently in [1,2]. This mechanism is based on the so-called Axion Quark Nugget (AQN) dark matter model, which was originally invented to explain the similarity of the dark and visible cosmological matter densities. We perform numerical simulations to evaluate the axion flux on the Earth's surface. We examine annual and daily modulations, which have been studied previously and are known to occur for any type of dark matter. We also discuss a novel type of short time enhancements which are unique to the AQN model: the statistical fluctuations and burst-like amplification, both of which can drastically amplify the axion signal, up to a factor $\sim10^2-10^3$ for a very short period of time. The present work studies the AQN-induced axions within the mass window $10^{-6}{\rm\,eV}\lesssim m_a\lesssim10^{-3}\rm\,eV$ with typical velocities $\langle v_a\rangle\sim0.6c$. We also comment on the broadband detection strategy to search for such relativistic axions by studying the daily and annual time modulations as well as random burst-like amplifications.

We construct the mean thermal Sunyaev-Zel'dovich (tSZ) Comptonization y profile around Luminous Red Galaxies (LRGs) in the redshift range 0.16 < z < 0.47 from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) using the Planck y map. The mean central tSZ signal for the full sample is y ~ 1.8 * 10^(-7) and we detect tSZ emission out to ~30 arcmin, which is well beyond the 10 arcmin angular resolution of the y map and well beyond the virial radii of the LRGs. We compare the measured profile with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations. This comparison agrees well for models that include feedback from active galactic nuclei (AGN), but not with hydrodynamic models without this energetic feedback mechanism. This suggests that an additional heating mechanism is required over SNe feedback and star formation to explain the y data profile. We also compare our results with predictions based on the halo model with a universal pressure profile (UPP) giving the y signal. The predicted profile is consistent with the data, but only if we account for the clustering of haloes via a two-halo term and if halo masses are estimated using the mean stellar-to-halo mass (SHM) relation of Coupon et al. (2015) or Wang et al.(2016) estimated from gravitational lensing measurements. We also discuss the importance of scatter in the SHM relation on the model predictions.

We present a tomographic cosmic shear analysis of the Kilo-Degree Survey (KiDS) combined with the VISTA Kilo-Degree Infrared Galaxy Survey (VIKING). This is the first time that a full optical to near-infrared data set has been used for a wide-field cosmological weak lensing experiment. This unprecedented data, spanning $450~$deg$^2$, allows us to improve significantly the estimation of photometric redshifts, such that we are able to include robustly higher-redshift sources for the lensing measurement, and - most importantly - solidify our knowledge of the redshift distributions of the sources. Based on a flat $\Lambda$CDM model we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.737_{-0.036}^{+0.040}$ in a blind analysis from cosmic shear alone. The tension between KiDS cosmic shear and the Planck-Legacy CMB measurements remains in this systematically more robust analysis, with $S_8$ differing by $2.3\sigma$. This result is insensitive to changes in the priors on nuisance parameters for intrinsic alignment, baryon feedback, and neutrino mass. KiDS shear measurements are calibrated with a new, more realistic set of image simulations and no significant B-modes are detected in the survey, indicating that systematic errors are under control. When calibrating our redshift distributions by assuming the 30-band COSMOS-2015 photometric redshifts are correct (following the Dark Energy Survey and the Hyper Suprime-Cam Survey), we find the tension with Planck is alleviated. The robust determination of source redshift distributions remains one of the most challenging aspects for future cosmic shear surveys.

We search the Planck data for a thermal Sunyaev-Zel'dovich (tSZ) signal due to gas filaments between pairs of Luminous Red Galaxies (LRG's) taken from the Sloan Digital Sky Survey Data Release 12 (SDSS/DR12). We identify $\sim$260,000 LRG pairs in the DR12 catalog that lie within 6-10 $h^{-1} \mathrm{Mpc}$ of each other in tangential direction and within 6 $h^{-1} \mathrm{Mpc}$ in radial direction. We stack pairs by rotating and scaling the angular positions of each LRG so they lie on a common reference frame, then we subtract a circularly symmetric halo from each member of the pair to search for a residual signal between the pair members. We find a statistically significant (5.3$\sigma$) signal between LRG pairs in the stacked data with a magnitude $\Delta y = (1.31 \pm 0.25) \times 10^{-8}$. The uncertainty is estimated from two Monte Carlo null tests which also establish the reliability of our analysis. Assuming a simple, isothermal, cylindrical filament model of electron over-density with a radial density profile proportional to $r_c/r$ (as determined from simulations), where $r$ is the perpendicular distance from the cylinder axis and $r_c$ is the core radius of the density profile, we constrain the product of over-density and filament temperature to be $\delta_c \times (T_{\rm e}/10^7 \, {\rm K}) \times (r_c/0.5h^{-1} \, {\rm Mpc}) = 2.7 \pm 0.5$. To our knowledge, this is the first detection of filamentary gas at over-densities typical of cosmological large-scale structure. We compare our result to the BAHAMAS suite of cosmological hydrodynamic simulations (McCarthy et al. 2017) and find a slightly lower, but marginally consistent Comptonization excess, $\Delta y = (0.84 \pm 0.24) \times 10^{-8}$.

We advocate for the idea that there is a fundamentally new mechanism for axion production on Earth, as recently suggested in Fischer et al. (2018) and Liang & Zhitnitsky (2018). We specifically focus on production of axions within Earth, with low velocities such that they will be trapped in the gravitational field. Our computations are based on the so-called Axion Quark Nugget (AQN) dark matter model, which was originally invented to explain the similarity of the dark and visible cosmological matter densities. This occurs in the model irrespective of the axion mass $m_\mathrm{a}$ or initial misalignment angle $\theta_0$. Annihilation of antimatter AQNs with visible matter inevitably produce axions when AQNs hit Earth. The emission rate of axions with velocities below escape velocity is very tiny compared to the overall emission, however these axions will be accumulated over the 4.5 billion year life time of the Earth, which greatly enhances the discovery potential. We perform numerical simulations with a realistically modeled incoming AQN velocity and mass distribution, and explore how AQNs interact as they travel through the interior of the Earth. We use this to estimate the axion flux on the surface of the Earth, the velocity-spectral features of trapped axions, the typical annihilation pattern of AQN, and the density profile of the axion halo around the Earth. Knowledge of these properties is necessary to make predictions for the observability of trapped axions using CAST, ADMX, MADMAX, CULTASK, ORPHEUS, ARIADNE, CASPEr, ABRACADABRA, QUAX, DM Radio.

We present a general method to compute the nonlinear matter power spectrum for dark energy and modified gravity scenarios with percent-level accuracy. By adopting the halo model and nonlinear perturbation theory, we predict the reaction of a $\Lambda$CDM matter power spectrum to the physics of an extended cosmological parameter space. By comparing our predictions to $N$-body simulations we demonstrate that with no-free parameters we can recover the nonlinear matter power spectrum for a wide range of different $w_0$-$w_a$ dark energy models to better than 1% accuracy out to $k \approx 1 \, h \, {\rm Mpc}^{-1}$. We obtain a similar performance for both DGP and $f(R)$ gravity, with the nonlinear matter power spectrum predicted to better than 3% accuracy over the same range of scales. When including direct measurements of the halo mass function from the simulations, this accuracy improves to 1%. With a single suite of standard $\Lambda$CDM $N$-body simulations, our methodology provides a direct route to constrain a wide range of non-standard extensions to the concordance cosmology in the high signal-to-noise nonlinear regime.

Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0<z<3$ to elucidate the nature of dark matter and dark energy. In this redshift range, galaxy formation is known to have a significant impact on observables such as two-point correlations of galaxy shapes and positions, altering their amplitude and scale dependence beyond the expected statistical uncertainty of upcoming experiments at separations under 10 Mpc. Successful extraction of information in such a regime thus requires, at the very least, unbiased models for the impact of galaxy formation on the matter distribution, and can benefit from complementary observational priors. This work reviews the current state of the art in the modelling of baryons for cosmology, from numerical methods to approximate analytical prescriptions, and makes recommendations for studies in the next decade, including a discussion of potential probe combinations that can help constrain the role of baryons in cosmological studies. We focus, in particular, on the modelling of the matter power spectrum, $P(k,z)$, as a function of scale and redshift, and of the observables derived from this quantity. This work is the result of a workshop held at the University of Oxford in November of 2018.

The Core Cosmology Library (CCL) provides routines to compute basic cosmological observables to a high degree of accuracy, which have been verified with an extensive suite of validation tests. Predictions are provided for many cosmological quantities, including distances, angular power spectra, correlation functions, halo bias and the halo mass function through state-of-the-art modeling prescriptions available in the literature. Fiducial specifications for the expected galaxy distributions for the Large Synoptic Survey Telescope (LSST) are also included, together with the capability of computing redshift distributions for a user-defined photometric redshift model. A rigorous validation procedure, based on comparisons between CCL and independent software packages, allows us to establish a well-defined numerical accuracy for each predicted quantity. As a result, predictions for correlation functions of galaxy clustering, galaxy-galaxy lensing and cosmic shear are demonstrated to be within a fraction of the expected statistical uncertainty of the observables for the models and in the range of scales of interest to LSST. CCL is an open source software package written in C, with a python interface and publicly available at https://github.com/LSSTDESC/CCL.

I examine differences in non-linear structure formation between cosmological models that share a $z=0$ linear power spectrum in both shape and amplitude, but that differ via their growth history. $N$-body simulations of these models display an approximately identical large-scale-structure skeleton, but reveal deeply non-linear differences in the demographics and properties of haloes. I investigate to what extent the spherical-collapse model can help in understanding these differences, in both real and redshift space. I discuss how this is difficult to do if one attempts to identify haloes directly, because in that case one is subject to the vagaries of halo finding algorithms. However, I demonstrate that the halo model of structure formation provides an accurate non-linear response in the power spectrum, but only if results from spherical collapse that include formation hysteresis are properly incorporated. I comment on how this fact can be used to provide per cent level accurate matter power spectrum predictions for dark energy models for $k\leq5\,h\mathrm{Mpc}^{-1}$ by using the halo model as a correction to accurate $\Lambda$CDM simulations. In the appendix I provide some fitting functions for the linear-collapse threshold ($\delta_\mathrm{c}$) and virialized overdensity ($\Delta_\mathrm{v}$) that are valid for a wide range of dark energy models. I also make my spherical-collapse code available at https://github.com/alexander-mead/collapse

Constraints on cosmological parameters from large-scale structure have traditionally been obtained from two-point statistics. However, non-linear structure formation renders these statistics insufficient in capturing the full information content available, necessitating the measurement of higher-order moments to recover information which would otherwise be lost. We construct quantities based on non-linear and non-local transformations of weakly non-Gaussian fields that Gaussianize the full multivariate distribution at a given order in perturbation theory. Our approach does not require a model of the fields themselves and takes as input only the first few polyspectra, which could be modelled or measured from simulations or data, making our method particularly suited to observables lacking a robust perturbative description such as the weak-lensing shear. We apply our method to simulated density fields, finding a significantly reduced bispectrum and an enhanced correlation with the initial field. We demonstrate that our method reconstructs a large proportion of the linear baryon acoustic oscillations, improving the information content over the raw field by 35%. We apply the transform to toy 21cm intensity maps, showing that our method still performs well in the presence of complications such as redshift-space distortions, beam smoothing, pixel noise, and foreground subtraction. We discuss how this method might provide a route to constructing a perturbative model of the fully non-Gaussian multivariate likelihood function.

We perform a combined analysis of cosmic shear tomography, galaxy-galaxy lensing tomography, and redshift-space multipole power spectra (monopole and quadrupole) using 450 deg$^2$ of imaging data by the Kilo Degree Survey (KiDS) overlapping with two spectroscopic surveys: the 2-degree Field Lensing Survey (2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS). We restrict the galaxy-galaxy lensing and multipole power spectrum measurements to the overlapping regions with KiDS, and self-consistently compute the full covariance between the different observables using a large suite of $N$-body simulations. We methodically analyze different combinations of the observables, finding that galaxy-galaxy lensing measurements are particularly useful in improving the constraint on the intrinsic alignment amplitude (by 30%, positive at $3.5\sigma$ in the fiducial data analysis), while the multipole power spectra are useful in tightening the constraints along the lensing degeneracy direction (e.g. factor of two stronger matter density constraint in the fiducial analysis). The fully combined constraint on $S_8 \equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3} = 0.742 \pm 0.035$, which is an improvement by 20% compared to KiDS alone, corresponds to a $2.6\sigma$ discordance with Planck, and is not significantly affected by fitting to a more conservative set of scales. Given the tightening of the parameter space, we are unable to resolve the discordance with an extended cosmology that is simultaneously favored in a model selection sense, including the sum of neutrino masses, curvature, evolving dark energy, and modified gravity. The complementarity of our observables allows for constraints on modified gravity degrees of freedom that are not simultaneously bounded with either probe alone, and up to a factor of three improvement in the $S_8$ constraint in the extended cosmology compared to KiDS alone.

We test extensions to the standard cosmological model with weak gravitational lensing tomography using 450 deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). In these extended cosmologies, which include massive neutrinos, nonzero curvature, evolving dark energy, modified gravity, and running of the scalar spectral index, we also examine the discordance between KiDS and cosmic microwave background measurements from Planck. The discordance between the two datasets is largely unaffected by a more conservative treatment of the lensing systematics and the removal of angular scales most sensitive to nonlinear physics. The only extended cosmology that simultaneously alleviates the discordance with Planck and is at least moderately favored by the data includes evolving dark energy with a time-dependent equation of state (in the form of the $w_0-w_a$ parameterization). In this model, the respective $S_8 = \sigma_8 \sqrt{\Omega_{\rm m}/0.3}$ constraints agree at the $1\sigma$ level, and there is `substantial concordance' between the KiDS and Planck datasets when accounting for the full parameter space. Moreover, the Planck constraint on the Hubble constant is wider than in LCDM and in agreement with the Riess et al. (2016) direct measurement of $H_0$. The dark energy model is moderately favored as compared to LCDM when combining the KiDS and Planck measurements, and remains moderately favored after including an informative prior on the Hubble constant. In both of these scenarios, marginalized constraints in the $w_0-w_a$ plane are discrepant with a cosmological constant at the $3\sigma$ level. Moreover, KiDS constrains the sum of neutrino masses to 4.0 eV (95% CL), finds no preference for time or scale dependent modifications to the metric potentials, and is consistent with flatness and no running of the spectral index. The analysis code is public at https://github.com/sjoudaki/kids450

We investigate the impact of astrophysical systematics on cosmic shear cosmological parameter constraints from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), and the concordance with cosmic microwave background measurements by Planck. We present updated CFHTLenS cosmic shear tomography measurements extended to degree scales using a covariance calibrated by a new suite of N-body simulations. We analyze these measurements with a new model fitting pipeline, accounting for key systematic uncertainties arising from intrinsic galaxy alignments, baryonic effects in the nonlinear matter power spectrum, and photometric redshift uncertainties. We examine the impact of the systematic degrees of freedom on the cosmological parameter constraints, both independently and jointly. When the systematic uncertainties are considered independently, the intrinsic alignment amplitude is the only degree of freedom that is substantially preferred by the data. When the systematic uncertainties are considered jointly, there is no consistently strong preference in favor of the more complex models. We quantify the level of concordance between the CFHTLenS and Planck datasets by employing two distinct data concordance tests, grounded in Bayesian evidence and information theory. We find that the two data concordance tests largely agree with one another, and that the level of concordance between the CFHTLenS and Planck datasets is sensitive to the exact details of the systematic uncertainties included in our analysis, ranging from decisive discordance to substantial concordance as the treatment of the systematic uncertainties becomes more conservative. The least conservative scenario is the one most favored by the cosmic shear data, but it is also the one that shows the greatest degree of discordance with Planck. The data and analysis code are public at https://github.com/sjoudaki/cfhtlens_revisited

We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). For a flat $\Lambda$CDM cosmology with a prior on $H_0$ that encompasses the most recent direct measurements, we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.745\pm0.039$. This result is in good agreement with other low redshift probes of large scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A $2.3$-$\sigma$ tension in $S_8$ and `substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved `self-calibrating' version of $lens$fit validated using an extensive suite of image simulations. Four-band $ugri$ photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov Chains are available at http://kids.strw.leidenuniv.nl.

We discuss 76 large amplitude transients (Delta-m>1.5) occurring in the nuclei of galaxies, nearly all with no previously known Active Galactic Nucleus (AGN). They have been discovered as part of the Pan-STARRS1 (PS1) 3pi survey, by comparison with SDSS photometry a decade earlier, and then monitored with the Liverpool Telescope, and studied spectroscopically with the William Herschel Telescope (WHT). Based on colours, light curve shape, and spectra, these transients fall into four groups. A few are misclassified stars or objects of unknown type. Some are red/fast transients and are known or likely nuclear supernovae. A few are either radio sources or erratic variables and so likely blazars. However the majority (~66%) are blue and evolve slowly, on a timescale of years. Spectroscopy shows them to be AGN at z~ 0.3 - 1.4, which must have brightened since the SDSS photometry by around an order of magnitude. It is likely that these objects were in fact AGN a decade ago, but too weak to be recognised by SDSS; they could then be classed as "hypervariable" AGN. By searching the SDSS Stripe 82 quasar database, we find 15 similar objects. We discuss several possible explanations for these slow blue hypervariables - (i) unusually luminous tidal disruption events; (ii) extinction events; (iii) changes in accretion state; and (iv) large amplitude microlensing by stars in foreground galaxies. A mixture of explanations (iii) and (iv) seems most likely. Both hold promise of considerable new insight into the AGN phenomenon.

We present an optimized variant of the halo model, designed to produce accurate matter power spectra well into the non-linear regime for a wide range of cosmological models. To do this, we introduce physically motivated free parameters into the halo-model formalism and fit these to data from high-resolution $N$-body simulations. For a variety of $\Lambda$CDM and $w$CDM models the halo-model power is accurate to $\simeq 5$ per cent for $k\leq 10h\,\mathrm{Mpc}$ and $z\leq 2$. An advantage of our new halo model is that it can be adapted to account for the effects of baryonic feedback on the power spectrum. We demonstrate this by fitting the halo model to power spectra from the OWLS hydrodynamical simulation suite via parameters that govern halo internal structure. We are able to fit all feedback models investigated at the 5 per cent level using only two free parameters, and we place limits on the range of these halo parameters for feedback models investigated by the OWLS simulations. Accurate predictions to high-$k$ are vital for weak lensing surveys, and these halo parameters could be considered nuisance parameters to marginalize over in future analyses to mitigate uncertainty regarding the details of feedback. Finally, we investigate how lensing observables predicted by our model compare to those from simulations and from HALOFIT for a range of $k$-cuts and feedback models and quantify the angular scales at which these effects become important. Code to calculate power spectra from the model presented in this paper can be found at https://github.com/alexander-mead/hmcode

We present an accurate non-linear matter power spectrum prediction scheme for a variety of extensions to the standard cosmological paradigm, which uses the tuned halo model previously developed in Mead (2015b). We consider dark energy models that are both minimally and non-minimally coupled, massive neutrinos and modified gravitational forces with chameleon and Vainshtein screening mechanisms. In all cases we compare halo-model power spectra to measurements from high-resolution simulations. We show that the tuned halo model method can predict the non-linear matter power spectrum measured from simulations of parameterised $w(a)$ dark energy models at the few per cent level for $k<10\,h\mathrm{Mpc}^{-1}$, and we present theoretically motivated extensions to cover non-minimally coupled scalar fields, massive neutrinos and Vainshtein screened modified gravity models that result in few per cent accurate power spectra for $k<10\,h\mathrm{Mpc}^{-1}$. For chameleon screened models we achieve only 10 per cent accuracy for the same range of scales. Finally, we use our halo model to investigate degeneracies between different extensions to the standard cosmological model, finding that the impact of baryonic feedback on the non-linear matter power spectrum can be considered independently of modified gravity or massive neutrino extensions. In contrast, considering the impact of modified gravity and massive neutrinos independently results in biased estimates of power at the level of 5 per cent at scales $k>0.5\,h\mathrm{Mpc}^{-1}$. An updated version of our publicly available HMcode can be found at https://github.com/alexander-mead/HMcode

We present deep radio continuum observations of the star-forming core of the Serpens South Infrared Dark Cloud with the Karl G. Jansky Very Large Array (VLA). Observations were conducted in two bands centered at 7.25 GHz (4.14 cm) and 4.75 GHz (6.31 cm) with an rms of 8.5 and 11.1 microJy/beam, respectively. We also use 2MASS, Spitzer and Herschel data to put our radio observations in the context of young stellar populations characterized by near and far infrared observations. Within a 5 arcmin x 5 arcmin region of interest around the central cluster, we detect roughly eighteen radio sources, seven of which we determine are protostellar in nature due to their radio spectral indices and their association with infrared sources. We find evidence for a previously undetected embedded Class 0 protostar and reaffirm Class 0 protostellar classifications determined by previous millimeter wavelength continuum studies. We use our infrared data to derive mid-infrared luminosities for three of our protostellar sources and find relative agreement between the known YSO radio luminosity vs bolometric luminosity correlation. Lastly, we marginally detect an additional six radio sources at the 2-3 sigma level that lie within two arcseconds of infrared YSO candidates, providing motivation for higher sensitivity studies to clarify the nature of these sources and further probe embedded and/or low luminosity YSOs in Serpens South.

Viable modifications of gravity that may produce cosmic acceleration need to be screened in high-density regions such as the Solar System, where general relativity is well tested. Screening mechanisms also prevent strong anomalies in the large-scale structure and limit the constraints that can be inferred on these gravity models from cosmology. We find that by suppressing the contribution of the screened high-density regions in the matter power spectrum, allowing a greater contribution of unscreened low densities, modified gravity models can be more readily discriminated from the concordance cosmology. Moreover, by variation of density thresholds, degeneracies with other effects may be dealt with more adequately. Specializing to chameleon gravity as a worked example for screening in modified gravity, employing N-body simulations of f(R) models and the halo model of chameleon theories, we demonstrate the effectiveness of this method. We find that a percent-level measurement of the clipped power at k < 0.3 h/Mpc can yield constraints on chameleon models that are more stringent than what is inferred from Solar System tests or distance indicators in unscreened dwarf galaxies. Finally, we verify that our method is also applicable to the Vainshtein mechanism.

We develop and test an algorithm to rescale a simulated dark-matter particle distribution or halo catalogue from a standard gravity model to that of a modified gravity model. This method is based on that of Angulo & White but with some additional ingredients to account for (i) scale-dependent growth of linear density perturbations and (ii) screening mechanisms that are generic features of viable modified gravity models. We attempt to keep the method as general as possible, so that it may plausibly be applied to a wide range of modified theories, although tests against simulations are restricted to a subclass of $f(R)$ models at this stage. We show that rescaling allows the power spectrum of matter to be reproduced at the $\sim 3$ per cent level in both real and redshift space up to $k=0.1h\,\mathrm{Mpc}^{-1}$ if we change the box size and alter the particle displacement field; this limit can be extended to $k=1h\,\mathrm{Mpc}^{-1}$ if we additionally alter halo internal structure. We simultaneously develop an algorithm that can be applied directly to a halo catalogue, in which case the halo mass function and clustering can be reproduced at the $\sim 5$ per cent level. Finally we investigate the clustering of halo particle distributions, generated from rescaled halo catalogues, and find that a similar accuracy can be reached.

We present and test a method for modifying the catalogue of dark matter haloes produced from a given cosmological simulation, so that it resembles the result of a simulation with an entirely different set of parameters. This extends the method of Angulo & White (2010), which rescales the full particle distribution from a simulation. Working directly with the halo catalogue offers an advantage in speed, and also allows modifications of the internal structure of the haloes to account for nonlinear differences between cosmologies. Our method can be used directly on a halo catalogue in a self contained manner without any additional information about the overall density field: although the large-scale displacement field is required by the method, this can be inferred from the halo catalogue alone. We show proof of concept of our method by rescaling a matter-only simulation with no baryon acoustic oscillation (BAO) features to a more standard LCDM model containing a cosmological constant and a BAO signal. In conjunction with the halo occupation approach, this method provides a basis for the rapid generation of mock galaxy samples spanning a wide range of cosmological parameters.

We discuss the extension to redshift space of a rescaling algorithm, designed to alter the effective cosmology of a pre-existing simulated particle distribution or catalogue of dark matter haloes. The rescaling approach was initially developed by Angulo & White and was adapted and applied to halo catalogues in real space in our previous work. This algorithm requires no information other than the initial and target cosmological parameters, and it contains no tuned parameters. It is shown here that the rescaling method also works well in redshift space, and that the rescaled simulations can reproduce the growth rate of cosmological density fluctuations appropriate for the target cosmology. Even when rescaling a grossly non-standard model with Lambda=0 and zero baryons, the redshift-space power spectrum of standard LCDM can be reproduced to about 5% error for k<0.2h Mpc^-1. The ratio of quadrupole-to-monopole power spectra remains correct to the same tolerance up to k=1h Mpc^-1, provided that the input halo catalogue contains measured internal velocity dispersions.

Recently, the Kepler Space Telescope has detected several planets in orbit around a close binary star system. These so-called circumbinary planets will experience non-trivial spatial and temporal distributions of radiative flux on their surfaces, with features not seen in their single-star orbiting counterparts. Earthlike circumbinary planets inhabited by photosynthetic organisms will be forced to adapt to these unusual flux patterns. We map the flux received by putative Earthlike planets (as a function of surface latitude/longitude and time) orbiting the binary star systems Kepler-16 and Kepler-47, two star systems which already boast circumbinary exoplanet detections. The longitudinal and latitudinal distribution of flux is sensitive to the centre of mass motion of the binary, and the relative orbital phases of the binary and planet. Total eclipses of the secondary by the primary, as well as partial eclipses of the primary by the secondary add an extra forcing term to the system. We also find that the patterns of darkness on the surface are equally unique. Beyond the planet's polar circles, the surface spends a significantly longer time in darkness than latitudes around the equator, due to the stars' motions delaying the first sunrise of spring (or hastening the last sunset of autumn). In the case of Kepler-47, we also find a weak longitudinal dependence for darkness, but this effect tends to average out if considered over many orbits. In the light of these flux and darkness patterns, we consider and discuss the prospects and challenges for photosynthetic organisms, using terrestrial analogues as a guide.

We present optical spectra and light curves for three hydrogen-poor super-luminous supernovae followed by the Public ESO Spectroscopic Survey of Transient Objects (PESSTO). Time series spectroscopy from a few days after maximum light to 100 days later shows them to be fairly typical of this class, with spectra dominated by Ca II, Mg II, Fe II and Si II, which evolve slowly over most of the post-peak photospheric phase. We determine bolometric light curves and apply simple fitting tools, based on the diffusion of energy input by magnetar spin-down, 56Ni decay, and collision of the ejecta with an opaque circumstellar shell. We investigate how the heterogeneous light curves of our sample (combined with others from the literature) can help to constrain the possible mechanisms behind these events. We have followed these events to beyond 100-200 days after peak, to disentangle host galaxy light from fading supernova flux and to differentiate between the models, which predict diverse behaviour at this phase. Models powered by radioactivity require unrealistic parameters to reproduce the observed light curves, as found by previous studies. Both magnetar heating and circumstellar interaction still appear to be viable candidates. A large diversity is emerging in observed tail-phase luminosities, with magnetar models failing in some cases to predict the rapid drop in flux. This would suggest either that magnetars are not responsible, or that the X-ray flux from the magnetar wind is not fully trapped. The light curve of one object shows a distinct re-brightening at around 100d after maximum light. We argue that this could result either from multiple shells of circumstellar material, or from a magnetar ionisation front breaking out of the ejecta.

Photosynthetic life requires sufficient photosynthetically active radiation (PAR) to metabolise. On Earth, plant behaviour, physiology and metabolism are sculpted around the night-day cycle by an endogenous biological circadian clock. The evolution of life was influenced by the Earth-Sun orbital dynamic, which generates the photo-environment incident on the planetary surface. In this work the unusual photo-environment of an Earth-like planet (ELP) in 3:2 spin orbit resonance is explored. Photo-environments on the ELP are longitudinally differentiated, in addition to differentiations relating to latitude and depth (for aquatic organisms) which are familiar on Earth. The light environment on such a planet could be compatible with Earth's photosynthetic life although the threat of atmospheric freeze-out and prolonged periods of darkness would present significant challenges. We emphasise the relationship between the evolution of life on a planetary body with its orbital dynamics.