We study specific star formation rate (sSFR) and gas profiles of star forming and green valley galaxies in the Simba cosmological hydrodynamic simulation. Star-forming galaxy half-light radii ($R_{\rm half}$) at $z=0$ and their evolution ($\propto(1+z)^{-0.78}$) agree with observations. Passive galaxy $R_{\rm half}$ agree with observations at high redshift, but by $z=0$ are too large, owing to numerical heating. We compare Simba $z=0$ sSFR radial profiles for star forming and green valley galaxies to observations. Simba shows strong central depressions in star formation rate (SFR), sSFR, and gas fraction in green valley galaxies and massive star-forming systems, qualitatively as observed, owing to black hole X-ray feedback, which pushes central gas outwards. Turning off X-ray feedback leads to centrally peaked sSFR profiles as in other simulations. In conflict with observations, Simba yields green valley galaxies with strongly dropping sSFR profiles beyond $\gtrsim R_{\rm half}$, regardless of AGN feedback. The central depression owes to lowering molecular gas content; the drop in the outskirts owes to reduced star formation efficiency. Simba's satellites have higher central sSFR and lower outskirts sSFR than centrals, in qualitative agreement with observations. At $z=2$ Simba does not show central depressions in massive star-forming galaxies, suggesting Simba's X-ray feedback should be more active at high-$z$. High resolution tests indicate central sSFR suppression is not sensitive to numerical resolution. Reproducing the central sSFR depression in $z=0$ green valley galaxies represents a unique success of Simba. The remaining discrepancies highlight the importance of SFR and gas profiles in constraining quenching mechanisms.

We provide a procedure for identification of dominating compressible and Alfvenic MHD modes or isotropic turbulence in synchrotron emission polarization maps of Galactic objects. The results for the region of North Galactic Pole, Orion molecular cloud complex and the star-forming complex Cygnus X are presented.

We present measurements of the two-dimensional genus of the SDSS-III BOSS catalogs to constrain cosmological parameters governing the shape of the matter power spectrum. The BOSS data are divided into twelve concentric shells over the redshift range $0.2 < z < 0.6$, and we extract the genus from the projected two-dimensional galaxy density fields. We compare the genus amplitudes to their Gaussian expectation values, exploiting the fact that this quantity is relatively insensitive to non-linear gravitational collapse. The genus amplitude provides a measure of the shape of the linear matter power spectrum, and is principally sensitive to $\Omega_{\rm c}h^{2}$ and scalar spectral index $n_{\rm s}$. A strong negative degeneracy between $\Omega_{\rm c}h^{2}$ and $n_{\rm s}$ is observed, as both can increase small scale power by shifting the peak and tilting the power spectrum respectively. We place a constraint on the particular combination $n_{\rm s}^{3/2} \Omega_{\rm c}h^{2}$ -- we find $n_{\rm s}^{3/2} \Omega_{\rm c}h^{2} = 0.1121 \pm 0.0043$ after combining the LOWZ and CMASS data sets, assuming a flat $\Lambda$CDM cosmology. This result is practically insensitive to reasonable variations of the power spectrum amplitude and linear galaxy bias. Our results are consistent with the Planck best fit $n_{\rm s}^{3/2}\Omega_{\rm c}h^{2} = 0.1139 \pm 0.0009$.

We study the morphology of the cosmic microwave background temperature and polarization fields using the shape and alignment parameters, $\beta$ and $\alpha$, that are constructed from the contour Minkowski tensor. The primary goal of our paper is to understand the effect of weak gravitational lensing on the morphology of the CMB fields. In order to isolate different physical effects that can be potentially confused with the effect of lensing, we first study the effect of varying the cosmology on $\alpha$ and $\beta$, and show that they are relatively insensitive to variation of cosmological parameters. Next we analyze the signatures of hemispherical anisotropy, and show that information of such anisotropy in $\alpha$ gets washed out at small angular scales and become pronounced only at large angular scales. For $\beta$ we find characteristic distortions which vary with the field threshold. We then study the effect of weak gravitational lensing using simulations of lensed temperature and $E$ and $B$ modes. We quantify the distortion induced in the fields across different angular scales. We find that lensing makes structures of all fields increasingly more anisotropic as we probe down to smaller scales. We find distinct behaviour of morphological distortions as a function of threshold for the different fields. The effect is small for temperature and $E$ mode, while it is significantly large for $B$ mode. Further, we find that lensing does not induce statistical anisotropy, as expected from the isotropic distribution of large scale structure of matter. We expect that the results obtained in this work will provide insights on the reconstruction of the lensing potential.

We present a Bayesian full-spectral-fitting analysis of 75 massive ($M_* > 10^{10.3} M_\odot$) UVJ-selected galaxies at redshifts of $1.0 < z < 1.3$, combining extremely deep rest-frame ultraviolet spectroscopy from VANDELS with multi-wavelength photometry. By the use of a sophisticated physical plus systematic uncertainties model, constructed within the Bagpipes code, we place strong constraints on the star-formation histories (SFHs) of individual objects. We firstly constrain the stellar mass vs stellar age relationship, finding a steep trend towards earlier average formation with increasing stellar mass of $1.48^{+0.34}_{-0.39}$ Gyr per decade in mass, although this shows signs of flattening at $M_* > 10^{11} M_\odot$. We show that this is consistent with other spectroscopic studies from $0 < z < 2$. This relationship places strong constraints on the AGN-feedback models used in cosmological simulations. We demonstrate that, although the relationships predicted by Simba and IllustrisTNG agree well with observations at $z=0.1$, they are too shallow at $z=1$, predicting an evolution of $<0.5$ Gyr per decade in mass. Secondly, we consider the connections between green-valley, post-starburst and quiescent galaxies, using our inferred SFH shapes and the distributions of galaxy physical properties on the UVJ diagram. The majority of our lowest-mass galaxies ($M_* \sim 10^{10.5} M_\odot$) are consistent with formation in recent ($z<2$), intense starburst events, with timescales of $\lesssim500$ Myr. A second class of objects experience extended star-formation epochs before rapidly quenching, passing through both green-valley and post-starburst phases. The most massive galaxies in our sample are extreme systems: already old by $z=1$, they formed at $z\sim5$ and quenched by $z=3$. However, we find evidence for their continued evolution through both AGN and rejuvenated star-formation activity.

We present the ensemble expectation values for the translation invariant, rank-2 Minkowski tensors in three-dimensions, for a linearly redshift space distorted Gaussian random field. The Minkowski tensors $W^{0,2}_{1}$, $W^{0,2}_{2}$ are sensitive to global anisotropic signals present within a field, and by extracting these statistics from the low redshift matter density one can place constraints on the redshift space distortion parameter $\beta = f/b$. We begin by reviewing the calculation of the ensemble expectation values $\langle W^{0,2}_{1} \rangle$, $\langle W^{0,2}_{2} \rangle $ for isotropic, Gaussian random fields, then consider how these results are modified by the presence of a linearly anisotropic signal. Under the assumption that all fields remain Gaussian, we calculate the anisotropic correction due to redshift space distortion in a coordinate system aligned with the line of sight, finding inequality between the diagonal elements of $\langle W^{0,2}_{1} \rangle $, $\langle W^{0,2}_{2} \rangle $. The ratio of diagonal elements of these matrices provides a set of statistics that are sensitive only to the redshift space distortion parameter $\beta$. We estimate the Fisher information that can be extracted from the Minkowski tensors, and find $W^{0,2}_{1}$ is more sensitive to $\beta$ than $W^{0,2}_{2}$, and a measurement of $W^{0,2}_{1}$ accurate to $\sim 1\%$ can yield a $\sim 4\%$ constraint on $\beta$. Finally, we discuss the difference between using the matrix elements of the Minkowski tensors directly against measuring the eigenvalues. For the purposes of cosmological parameter estimation we advocate the use of the matrix elements, to avoid spurious anisotropic signals that can be generated by the eigenvalue decomposition.

The equipartition of magnetic and thermal energy in the interstellar medium (ISM) indicates the magneto-hydrodynamic nature of the interstellar turbulence, which can be decomposed into three wave modes: Alfv\'en, fast and slow magnetosonic modes\cite{Herlofson50,LG01,CL03}. Even for studies performed in the case of subsonic turbulence, the magnetosonic modes have been revealed to be more effective than the Alfv\'enic modes for processes such as cosmic ray (CR) transport and acceleration \cite{YL04,Lynn14}. The multiphase nature of ISM and diversity of driving mechanisms give rise to spatial variation of turbulence properties. Nevertheless, the employed model of magneto-hydrodynamic turbulence is often oversimplified being assumed to be only Alfv\'enic or even hydrodynamic due to a lack of observational evidence. Here we report the employment of our novel method, the signature from polarization analysis (SPA), on unveiling the plasma modes in interstellar turbulence. Its application leads to the first discovery of magnetosonic modes in the Cygnus X region. It is found that the magnetosonic modes overlap to a high degree with Fermi cocoon, consistent with theoretical expectations. Moreover, through comparison with the spectrum at other wavelengths, the plasma modes of turbulence are unveiled in active star formation zones. The SPA casts light on the plasma modes composition of the Galactic turbulence, and marks the onset of a new era in the study of interstellar turbulence and accordingly our understandings of relevant processes including cosmic ray transport and star formation.

We use morphological descriptors, Betti numbers and Contour Minkowski Tensor (CMT) on 21cm brightness temperature excursion sets, to study the ionization and heating history of the intergalactic medium (IGM) during and before the Epoch of Reionization (EoR). The ratio of eigenvalues of the CMT denoted by $\beta$, gives shape information while it's trace gives the contour length of holes and connected regions. We simulate the matter density, neutral hydrogen fraction, spin temperature and brightness temperature field using the publicly available code 21cmFAST in a redshift range of $z=20.22$ to $z=6$. We study the redshift evolution of three quantities - the Betti number counts $N_{con,hole}$, the characteristic size $r^{ch}_{con,hole}$ and shape anisotropy parameter $\beta^{ch}_{con,hole}$ of connected regions and holes for these fields and investigate the different physical origins of their evolution. We make a qualitative comparison of different models of heating and ionization during the EoR. We obtain different regimes of morphological evolution of brightness temperature, depending upon how the shapes and sizes of connected regions and holes change with redshift for different astrophysical settings affecting the ionization and heating history of the IGM during and before the EoR. We find that the morphology of the brightness temperature field traces the morphology of ionized regions below a certain redshift value depending upon the model, where $\Delta r^{ch}_{hole}<10 \%$ and $\Delta \beta^{ch}_{hole}<1 \%$ relative to the $x_{HI}$ field. This difference decreases with redshift. Therefore, the ionization history of the IGM can be reconstructed using the morphological description of $\delta T_b$ in real space.

We introduce the Simba simulations, the next generation of the Mufasa cosmological galaxy formation simulations run with Gizmo's meshless finite mass hydrodynamics. Simba includes updates to Mufasa's sub-resolution star formation and feedback prescriptions, and introduces black hole growth via the torque-limited accretion model of Angl\'es-Alc\'azar et al. (2017) from cold gas and Bondi accretion from hot gas, along with black hole feedback via kinetic bipolar outflows and X-ray energy. Ejection velocities are taken to be ~10^3 km/s at high Eddington ratios, increasing to ~8000 km/s at Eddington ratios below 2%, with a constant momentum input of 20L/c. Simba further includes an on-the-fly dust production, growth, and destruction model. Our Simba run with (100 Mpc/h)^3 and 1024^3 gas elements reproduces numerous observables, including galaxy stellar mass functions at z=0-6, the stellar mass--star formation rate main sequence, HI and H2 fractions, the mass-metallicity relation at z=0 and z=2, star-forming galaxy sizes, hot gas fractions in massive halos, and z=0 galaxy dust properties. However, Simba also yields an insufficiently sharp truncation of the z=0 mass function, and too-large sizes for low-mass quenched galaxies. We show that Simba's jet feedback is primarily responsible for quenching massive galaxies.

We apply the Minkowski tensor statistics to three dimensional Gaussian random fields. Minkowski tensors contain information regarding the orientation and shape of excursion sets, that is not present in the scalar Minkowski functionals. They can be used to quantify globally preferred directions, and additionally provide information on the mean shape of subsets of a field. This makes them ideal statistics to measure the anisotropic signal generated by redshift space distortion in the low redshift matter density field. We review the definition of the Minkowski tensor statistics in three dimensions, focusing on two coordinate invariant quantities $W^{0,2}_{1}$ and $W^{0,2}_{2}$. We calculate the ensemble average of these $3 \times 3$ matrices for an isotropic Gaussian random field, finding that they are proportional to products of the identity matrix and a corresponding scalar Minkowski functional. We show how to numerically reconstruct $W^{0,2}_{1}$ and $W^{0,2}_{2}$ from discretely sampled fields and apply our algorithm to isotropic Gaussian fields generated from a linear $\Lambda$CDM matter power spectrum. We then introduce anisotropy by applying a linear redshift space distortion operator to the matter density field, and find that both $W^{0,2}_{1}$ and $W^{0,2}_{2}$ exhibit a distinct signal characterised by inequality between their diagonal components. We discuss the physical origin of this signal and how it can be used to constrain the redshift space distortion parameter $\Upsilon \equiv f/b$.

We apply the Minkowski Tensor statistics to two dimensional slices of the three dimensional density field. The Minkowski Tensors are a set of functions that are sensitive to directionally dependent signals in the data, and furthermore can be used to quantify the mean shape of density peaks. We begin by introducing our algorithm for constructing bounding perimeters around subsets of a two dimensional field, and reviewing the definition of Minkowski Tensors. Focusing on the translational invariant statistic $W^{1,1}_{2}$ - a $2 \times 2$ matrix - we calculate its eigenvalues for both the entire excursion set ($\Lambda_{1},\Lambda_{2}$) and for individual connected regions and holes within the set ($\lambda_{1},\lambda_{2}$). The ratio of eigenvalues $\Lambda_{2}/\Lambda_{1}$ informs us of the presence of global anisotropies in the data, and $\langle \lambda_{2}/\lambda_{1} \rangle$ is a measure of the mean shape of peaks and troughs in the density field. We study these quantities for a Gaussian field, then consider how they are modified by the effect of gravitational collapse using the latest Horizon Run 4 cosmological simulation. We find $\Lambda_{1,2}$ are essentially independent of gravitational collapse, as the process maintains statistical isotropy. However, the mean shape of peaks is modified significantly - overdensities become relatively more circular compared to underdensities of the same area. When applying the statistic to a redshift space distorted density field, we find a significant signal in the eigenvalues $\Lambda_{1,2}$, suggesting that they can be used to probe the large-scale velocity field.

Self tuning is one of the few methods for dynamically cancelling a large cosmological constant and yet giving an accelerating universe. Its drawback is that it tends to screen all sources of energy density, including matter. We develop a model that tempers the self tuning so the dynamical scalar field still cancels an arbitrary cosmological constant, including the vacuum energy through any high energy phase transitions, without affecting the matter fields. The scalar-tensor gravitational action is simple, related to cubic Horndeski gravity, with a nonlinear derivative interaction plus a tadpole term. Applying shift symmetry and using the property of degeneracy of the field equations we find families of functions that admit de Sitter solutions with expansion rates that are independent of the magnitude of the cosmological constant and preserve radiation and matter dominated phases. That is, the method can deliver a standard cosmic history including current acceleration, despite the presence of a Planck scale cosmological constant.

We study the topology of the matter density field in two dimensional slices, and consider how we can use the amplitude $A$ of the genus for cosmological parameter estimation. Using the latest Horizon Run 4 simulation data, we calculate the genus of the smoothed density field constructed from lightcone mock galaxy catalogs. Information can be extracted from the amplitude of the genus by considering both its redshift evolution and magnitude. The constancy of the genus amplitude with redshift can be used as a standard population, from which we derive constraints on the equation of state of dark energy $w_{\rm de}$ - by measuring $A$ at $z \sim 0.1$ and $z \sim 1$, we can place an order $\Delta w_{\rm de} \sim {\cal O}(15\%)$ constraint on $w_{\rm de}$. By comparing $A$ to its Gaussian expectation value we can potentially derive an additional stringent constraint on the matter density $\Delta \Omega_{\rm mat} \sim 0.01$. We discuss the primary sources of contamination associated with the two measurements - redshift space distortion and shot noise. With accurate knowledge of galaxy bias, we can successfully remove the effect of redshift space distortion, and the combined effect of shot noise and non-linear gravitational evolution is suppressed by smoothing over suitably large scales $R_{\rm G} \ge 15 {\rm Mpc}/h$. Without knowledge of the bias, we discuss how joint measurements of the two and three dimensional genus can be used to constrain the growth factor $\beta = f/b$. The method can be applied optimally to redshift slices of a galaxy distribution generated using the drop-off technique.

We employ the rank-2 {\em contour} Minkowski Tensor in two dimensions to probe length and time scales of ionized bubbles during the epoch of reionization. We demonstrate that the eigenvalues of this tensor provide excellent probes of the distribution of the sizes of ionized bubbles, and from it the characteristic bubble sizes, at different redshifts. We show that ionized bubbles are not circular, and hence not spherical in three dimensions, as is often assumed for simplified analytic arguments. We quantify their shape anisotropy by using the ratio of the two eigenvalues. The shape parameter provides the characteristic time epochs when bubble mergers begin and end. Our method will be very useful to reconstruct the reionization history using data of the brightness temperature field.

We generalize the translation invariant tensor-valued Minkowski Functionals which are defined on two-dimensional flat space to the unit sphere. We apply them to level sets of random fields. The contours enclosing boundaries of level sets of random fields give a spatial distribution of random smooth closed curves. We obtain analytic expressions for the ensemble expectation values for the matrix elements of the tensor-valued Minkowski Functionals for isotropic Gaussian and Rayleigh fields. We elucidate the way in which the elements of the tensor Minkowski Functionals encode information about the nature and statistical isotropy (or departure from isotropy) of the field. We then implement our method to compute the tensor-valued Minkowski Functionals numerically and demonstrate how they encode statistical anisotropy and departure from Gaussianity by applying the method to maps of the Galactic foreground emissions from the PLANCK data.

We study the two-dimensional topology of the galactic distribution when projected onto two-dimensional spherical shells. Using the latest Horizon Run 4 simulation data, we construct the genus of the two-dimensional field and consider how this statistic is affected by late-time nonlinear effects -- principally gravitational collapse and redshift space distortion (RSD). We also consider systematic and numerical artifacts such as shot noise, galaxy bias, and finite pixel effects. We model the systematics using a Hermite polynomial expansion and perform a comprehensive analysis of known effects on the two-dimensional genus, with a view toward using the statistic for cosmological parameter estimation. We find that the finite pixel effect is dominated by an amplitude drop and can be made less than $1\%$ by adopting pixels smaller than $1/3$ of the angular smoothing length. Nonlinear gravitational evolution introduces time-dependent coefficients of the zeroth, first, and second Hermite polynomials, but the genus amplitude changes by less than $1\%$ between $z=1$ and $z=0$ for smoothing scales $R_{\rm G} > 9 {\rm Mpc/h}$. Non-zero terms are measured up to third order in the Hermite polynomial expansion when studying RSD. Differences in shapes of the genus curves in real and redshift space are small when we adopt thick redshift shells, but the amplitude change remains a significant $\sim {\cal O}(10\%)$ effect. The combined effects of galaxy biasing and shot noise produce systematic effects up to the second Hermite polynomial. It is shown that, when sampling, the use of galaxy mass cuts significantly reduces the effect of shot noise relative to random sampling.

We study features in the bispectrum of the primordial curvature perturbation correlated with the reconstructed primordial power spectrum from the observed cosmic microwave background temperature data. We first show how the bispectrum can be completely specified in terms of the power spectrum and its first two derivatives, valid for any configuration of interest. Then using a model-independent reconstruction of the primordial power spectrum from the Planck angular power spectrum of temperature anisotropies, we compute the bispectrum in different triangular configurations. We find that in the squeezed limit at k ~ 0.06/Mpc and k ~ 0.014/Mpc there are marginal 2sigma deviations from the standard featureless bispectrum, which meanwhile is consistent with the reconstructed bispectrum in the equilateral configuration.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We test the isotropy of the local Universe using low redshift supernova data from various catalogs and the non-parametric method of smoothed residuals. Using a recently developed catalog which combines supernova data from various surveys, we show that the isotropic hypothesis of a Universe with zero velocity perturbation can be rejected with moderate significance, with $p$-value $\sim 0.07$ out to redshift $z < 0.045$. We estimate the direction of maximal anisotropy on the sky for various pre-existing catalogs and show that it remains relatively unaffected by the light curve fitting procedure. However the recovered direction is biased by the underlying distribution of data points on the sky. We estimate both the uncertainty and bias in the direction by creating mock data containing a randomly oriented bulk flow and using our method to reconstruct its direction. We conclude that the inhomogeneous nature of the data introduces a directional bias in galactic latitude of approximately $|\Delta b_{\rm max}| \sim 18^{\circ}$ for the supernova catalog considered in this work, and after correcting for this effect we infer the direction of maximum anisotropy as $(b,\ell) = (20^{\circ}, 276^{\circ}) \pm (12^{\circ},29^{\circ})$ in galactic coordinates. Finally we compare the anisotropic signal in the data to mock realisations in which large scale velocity perturbations are consistently accounted for at the level of linear perturbation theory. We show that including the effect of the velocity perturbation in our mock catalogs degrades the significance of the anisotropy considerably, with $p$-value increasing to $\sim 0.29$. One can conclude from our analysis that there is a moderate deviation from isotropy in the supernova data, but the signal is consistent with a large scale bulk velocity expected within $\Lambda$CDM.

We search for self tuning solutions to the Einstein-scalar field equations for the simplest class of `Fab-Four' models with constant potentials. We first review the conditions under which self tuning occurs in a cosmological spacetime, and by introducing a small modification to the original theory - introducing the second and third Galileon terms - show how one can obtain de Sitter states where the expansion rate is independent of the vacuum energy. We then consider whether the same self tuning mechanism can persist in a spherically symmetric inhomogeneous spacetime. We show that there are no asymptotically flat solutions to the field equations in which the vacuum energy is screened, other than the trivial one (Minkowski space). We then consider the possibility of constructing Schwarzschild de Sitter spacetimes for the modified Fab Four plus Galileon theory. We argue that the only model that can successfully screen the vacuum energy in both an FLRW and Schwarzschild de Sitter spacetime is one containing `John' $\sim G^{\mu}{}_{\nu} \partial_{\mu}\phi\partial^{\nu}\phi$ and a canonical kinetic term $\sim \partial_{\alpha}\phi \partial^{\alpha}\phi$. This behaviour was first observed in (Babichev&Charmousis,2013). The screening mechanism, which requires redundancy of the scalar field equation in the `vacuum', fails for the `Paul' term in an inhomogeneous spacetime.

We test the isotropy of the local distribution of galaxies using the 2MASS extended source catalogue. By decomposing the full sky survey into distinct patches and using a combination of photometric and spectroscopic redshift data, we use both parametric and non-parametric methods to obtain the shape of the luminosity function in each patch. We use the shape of the luminosity function to test the statistical isotropy of the underlying galaxy distribution. The parametric estimator shows some evidence of a hemispherical asymmetry in the north/south Galactic plane. However the non-parametric estimator exhibits no significant anisotropy, with the galaxy distribution being consistent with the assumption of isotropy in all regions considered. The parametric asymmetry is attributed to the relatively poor fit of the functional form to the underlying data. When using the non-parametric estimator, we do find a dipole in the shape of the luminosity function, with maximal deviation from isotropy at galactic coordinate $(b,l)=(30^{\circ},315^{\circ})$. However we can ascribe no strong statistical significance to this observation.

We discuss some details regarding the method of smoothed residuals, which has recently been used to search for anisotropic signals in low-redshift distance measurements (Supernovae). In this short note we focus on some details regarding the implementation of the method, particularly the issue of effectively detecting signals in data that are inhomogeneously distributed on the sky. Using simulated data, we argue that the original method proposed in Colin et al. [1] will not detect spurious signals due to incomplete sky coverage, and that introducing additional Gaussian weighting to the statistic as in [2] can hinder its ability to detect a signal. Issues related to the width of the Gaussian smoothing are also discussed.

We compute the critical density of collapse for spherically symmetric overdensities in a class of f(R) modified gravity models. For the first time we evolve the Einstein, scalar field and non-linear fluid equations, making the minimal simplifying assumptions that the metric potentials and scalar field remain quasi-static throughout the collapse. Initially evolving a top hat profile, we find that the density threshold for collapse depends significantly on the initial conditions imposed, specifically the choice of size and shape. By imposing `natural' initial conditions, we obtain a fitting function for the spherical collapse delta_c as a function of collapse redshift, mass of the overdensity and f_{R0}, the background scalar field value at z=0. By extending delta_c into drifting and diffusing barrier within the context of excursion set theory, we obtain a realistic mass function that might be used to confront this class of scalar-tensor models with observations of dark matter halos. The proposed analytic formula for the halo mass function was tested against Monte Carlo random walks for a wide class of moving barriers and can therefore be applied to other modified gravity theories.

Cluster abundances are oddly insensitive to canonical early dark energy. Early dark energy with sound speed equal to the speed of light cannot be distinguished from a quintessence model with the equivalent expansion history for $z<2$ but negligible early dark energy density, despite the different early growth rate. However, cold early dark energy, with a sound speed much smaller than the speed of light, can give a detectable signature. Combining cluster abundances with cosmic microwave background power spectra can determine the early dark energy fraction to 0.3 % and distinguish a true sound speed of 0.1 from 1 at 99 % confidence. We project constraints on early dark energy from the Euclid cluster survey, as well as the Dark Energy Survey, using both current and projected Planck CMB data, and assess the impact of cluster mass systematics. We also quantify the importance of dark energy perturbations, and the role of sound speed during a crossing of $w=-1$.

Wide area cosmological surveys enable investigation of whether dark energy properties are the same in different directions on the sky. Cosmic microwave background observations strongly restrict any dynamical effects from anisotropy, in an integrated sense. For more local constraints we compute limits from simulated distance measurements for various distributions of survey fields in a Bianchi I anisotropic universe. We then consider the effects of fitting for line of sight properties where isotropic dynamics is assumed (testing the accuracy through simulations) and compare sensitivities of observational probes for anisotropies, from astrophysical systematics as well as dark energy. We also point out some interesting features of anisotropic expansion in Bianchi I cosmology.

Galileon gravity is a robust theoretical alternative to general relativity with a cosmological constant for explaining cosmic acceleration, with interesting properties such as having second order field equations and a shift symmetry. While either its predictions for the cosmic expansion or growth histories can approach standard \Lambda CDM, we demonstrate the incompatibility of both doing so simultaneously. Already current observational constraints can severely disfavor an entire class of Galileon gravity models that do not couple directly to matter, ruling them out as an alternative to \Lambda CDM.

We investigate circumstances under which one can generalize Horndeski's most general scalar-tensor theory of gravity. Specifically we demonstrate that a nonlinear combination of purely kinetic gravity terms can give rise to an accelerating universe without the addition of extra propagating degrees of freedom on cosmological backgrounds, and exhibit self tuning to bring a large cosmological constant under control. This nonlinear approach leads to new properties that may be instructive for exploring the behaviors of gravity.

Galileon gravity offers a robust gravitational theory for explaining cosmic acceleration, having a rich phenomenology of testable behaviors. We explore three classes of Galileon models -- standard uncoupled, and linearly or derivatively coupled to matter -- investigating the expansion history with particular attention to early time and late time attractors, as well as the linear perturbations. From the relativistic and nonrelativistic Poisson equations we calculate the generalizations of the gravitational strength (Newton's constant), deriving its early and late time behavior. By scanning through the parameters we derive distributions of the gravitational strength at various epochs and trace the paths of gravity in its evolution. Using ghost-free and stability criteria we restrict the allowed parameter space, finding in particular that the linear and derivative coupled models are severely constrained by classical instabilities in the early universe.

We study $F(R)$ modified gravity models which are capable of driving the accelerating epoch of the Universe at the present time whilst not destroying the standard Big Bang and inflationary cosmology. Recent studies have shown that a weak curvature singularity with $|R|\to\infty$ can arise generically in viable $F(R)$ models of present dark energy (DE) signaling an internal incompleteness of these models. In this work we study how this problem is cured by adding a quadratic correction with a sufficiently small coefficient to the $F(R)$ function at large curvatures. At the same time, this correction eliminates two more serious problems of previously constructed viable $F(R)$ DE models: unboundedness of the mass of a scalar particle (scalaron) arising in $F(R)$ gravity and the scalaron overabundance problem. Such carefully constructed models can also yield both an early time inflationary epoch and a late time de Sitter phase with vastly different values of $R$. The reheating epoch in these combined models of primordial and present dark energy is completely different from that of the old $R + R^{2}/6M^{2}$ inflationary model, mainly due to the fact that values of the effective gravitational constant at low and intermediate curvatures are different for positive and negative $R$. This changes the number of e-folds during the observable part of inflation that results in a different value of the primordial power spectrum index.

We examine general physical parameterisations for viable gravitational models in the $f(R)$ framework. This is related to the mass of an additional scalar field, called the scalaron, that is introduced by the theories. Using a simple parameterisation for the scalaron mass $M(a)$ we show there is an exact correspondence between the model and popular parameterisations of the modified Poisson equation $\mu(a,k)$ and the ratio of the Newtonian potentials $\eta(a,k)$. However, by comparing the aforementioned model against other viable scalaron theories we highlight that the common form of $\mu(a,k)$ and $\eta(a,k)$ in the literature does not accurately represent $f(R)$ behaviour. We subsequently construct an improved description for the scalaron mass (and therefore $\mu(a,k)$ and $\eta(a,k)$) which captures their essential features and has benefits derived from a more physical origin. We study the scalaron's observational signatures and show the modification to the background Friedmann equation and CMB power spectrum to be small. We also investigate its effects in the linear and non linear matter power spectrum--where the signatures are evident--thus giving particular importance to weak lensing as a probe of these models. Using this new form, we demonstrate how the next generation Euclid survey will constrain these theories and its complementarity to current solar system tests. In the most optimistic case Euclid, together with a Planck prior, can constrain a fiducial scalaron mass $M_{0} = 9.4 \times 10^{-30}{\rm eV}$ at the $\sim 20 %$ level. However, the decay rate of the scalaron mass, with fiducial value $\nu = 1.5$, can be constrained to $\sim 3%$ uncertainty.

We study the evolution of density perturbations for a class of $f(R)$ models which closely mimic $\Lambda$CDM background cosmology. Using the quasi-static approximation, and the fact that these models are equivalent to scalar-tensor gravity, we write the modified Friedmann and cosmological perturbation equations in terms of the mass $M$ of the scalar field. Using the perturbation equations, we then derive an analytic expression for the growth parameter $\gamma$ in terms of $M$, and use our result to reconstruct the linear matter power spectrum. We find that the power spectrum at $z \sim 0$ is characterized by a tilt relative to its General Relativistic form, with increased power on small scales. We discuss how one has to modify the standard, constant $\gamma$ prescription in order to study structure formation for this class of models. Since $\gamma$ is now scale and time dependent, both the amplitude and transfer function associated with the linear matter power spectrum will be modified. We suggest a simple parameterization for the mass of the scalar field, which allows us to calculate the matter power spectrum for a broad class of $f(R)$ models.

If the equation of state of dark energy is anisotropic there will be additional quadrupole anisotropy in the cosmic microwave background induced by the time dependent anisotropic stress quantified in terms of $\Delta w$. Assuming that the entire amplitude of the observed quadrupole is due to this anisotropy, we conservatively impose a limit of $|\Delta w| < 2.1\times 10^{-4}$ for any value of $w\ge -1$ assuming that $\Omega_{\rm m}<0.5$. This is considerably tighter than that which comes from SNe. Stronger limits, upto a factor of 10, are possible for specific values of $\Omega_{\rm m}$ and $w$. Since we assume this component is uncorrelated with the stochastic component from inflation, we find that both the expectation value and the sample variance are increased. There no improvement in the likelihood of an anomalously low quadrupole as suggested by previous work on an elliptical universe.

We study cosmological expansion in F(R) gravity using the trace of the field equations. High frequency oscillations in the Ricci scalar, whose amplitude increase as one evolves backward in time, have been predicted in recent works. We show that the approximations used to derive this result very quickly breakdown in any realistic model due to the non-linear nature of the underlying problem. Using a combination of numerical and semi-analytic techniques, we study a range of models which are otherwise devoid of known pathologies. We find that high frequency asymmetric oscillations and a singularity at finite time appear to be present for a wide range of initial conditions. We show that this singularity can be avoided with a certain range of initial conditions, which we find by evolving the models forwards in time. In addition we show that the oscillations in the Ricci scalar are highly suppressed in the Hubble parameter and scale factor.

Modified gravity models are subject to a number of consistency requirements which restrict the form that the function $F(R)$ can take. We study a particular class of $F(R)$ functions which satisfy various constraints that have been found in the literature. These models have a late time accelerating epoch, and an acceptable matter era. We calculate the Friedmann equation for our models, and show that in order to satisfy the constraints we impose, they must mimic General Relativity plus $\Lambda$ throughout the cosmic history, with exponentially suppressed corrections. We also find that the free parameters in our model must be fine tuned to obtain an acceptable late time accelerating phase. We discuss the generality of this conclusion.