We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.Comment: 15 pages, 13 figures, to be submitted to MNRAS. For high-resolution figures, please see http://www-theorie.physik.unizh.ch/~agertz

We study the importance of baryonic physics on predictions of the matter power spectrum as it is relevant for forthcoming weak lensing surveys. We quantify the impact of baryonic physics using a set of three cosmological numerical simulations. Each simulation has the same initial density field, but models a different set of physical processes. The first simulation evolves the density field using gravity alone, the second includes non-radiative gasdynamics, and the third includes radiative heating and cooling of baryons, star formation, and supernova feedback. We find that baryonic processes alter predictions for the matter power spectrum significantly relative to models that include only gravitational interactions. Our results imply that future weak lensing experiments such as LSST and SNAP will be very sensitive to the poorly-understood physics governing the nonlinear evolution of the baryonic component of the universe. The net effect is significantly larger in the case of the model with cooling and star formation, in which case our results imply that contemporary surveys such as the CFHT Wide survey may also be sensitive to baryonic processes. In particular, this effect could be important for forecasts of the constraining power of future surveys if information from scales larger than l ~ 1000 is included in the analysis. We find that deviations are caused primarily by the rearrangement of matter within individual dark matter halos relative to the gravity-only case, rather than a large-scale rearrangement of matter. Consequently, we propose a simple model, based on the phenomenological halo model of dark matter clustering, for baryonic effects that can be used to aid in the interpretation of forthcoming weak lensing data.Comment: 14 pages, 8 figures. Submitted to Ap

Cosmological parameter estimation is entering a new era. Large collaborations need to coordinate high-stakes analyses using multiple methods; furthermore such analyses have grown in complexity due to sophisticated models of cosmology and systematic uncertainties. In this paper we argue that modularity is the key to addressing these challenges: calculations should be broken up into interchangeable modular units with inputs and outputs clearly defined. We present a new framework for cosmological parameter estimation, CosmoSIS, designed to connect together, share, and advance development of inference tools across the community. We describe the modules already available in CosmoSIS, including camb, Planck, cosmic shear calculations, and a suite of samplers. We illustrate it using demonstration code that you can run out-of-the-box with the installer available at http://bitbucket.org/joezuntz/cosmosis.

Current measurements of the low and high redshift Universe are in tension if we restrict ourselves to the standard six-parameter model of flat ΛCDM. This tension has two parts. First, the Planck satellite data suggest a higher normalization of matter perturbations than local measurements of galaxy clusters. Second, the expansion rate of the Universe today, H0, derived from local distance-redshift measurements is significantly higher than that inferred using the acoustic scale in galaxy surveys and the Planck data as a standard ruler. The addition of a sterile neutrino species changes the acoustic scale and brings the two into agreement; meanwhile, adding mass to the active neutrinos or to a sterile neutrino can suppress the growth of structure, bringing the cluster data into better concordance as well. For our fiducial data set combination, with statistical errors for clusters, a model with a massive sterile neutrino shows 3.5σ evidence for a nonzero mass and an even stronger rejection of the minimal model. A model with massive active neutrinos and a massless sterile neutrino is similarly preferred. An eV-scale sterile neutrino mass--of interest for short baseline and reactor anomalies--is well within the allowed range. We caution that (i) unknown astrophysical systematic errors in any of the data sets could weaken this conclusion, but they would need to be several times the known errors to eliminate the tensions entirely; (ii) the results we find are at some variance with analyses that do not include cluster measurements; and (iii) some tension remains among the data sets even when new neutrino physics is included.