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

In light of the recent BICEP2 B-mode polarization detection, which implies a large inflationary tensor-to-scalar ratio r_{0.05}=0.2^{+0.07}_{-0.05}, we re-examine the evidence for an extra sterile massive neutrino, originally invoked to account for the tension between the cosmic microwave background (CMB) temperature power spectrum and local measurements of the expansion rate H0 and cosmological structure. With only the standard active neutrinos and power-law scalar spectra, this detection is in tension with the upper limit of r<0.11 (95% confidence) from the lack of a corresponding low multipole excess in the temperature anisotropy from gravitational waves. An extra sterile species with the same energy density as is needed to reconcile the CMB data with H0 measurements can also alleviate this new tension. By combining data from the Planck and ACT/SPT temperature spectra, WMAP9 polarization, H_0, baryon acoustic oscillation and local cluster abundance measurements with BICEP2 data, we find the joint evidence for a sterile massive neutrino increases to DeltaNeff=0.98\pm 0.26 for the effective number and ms= 0.52\pm 0.13 eV for the effective mass or 3.8 sigma and 4 sigma evidence respectively. We caution the reader that these results correspond to a joint statistical evidence and, in addition, astrophysical systematic errors in the clusters and H0 measurements, and small-scale CMB data could weaken our conclusions.

The use of galaxy clusters as cosmological probes hinges on our ability to measure their masses accurately and with high precision. Hydrostatic mass is one of the most common methods for estimating the masses of individual galaxy clusters, which suffer from biases due to departures from hydrostatic equilibrium. Using a large, mass-limited sample of massive galaxy clusters from a high-resolution hydrodynamical cosmological simulation, in this work we show that in addition to turbulent and bulk gas velocities, acceleration of gas introduces biases in the hydrostatic mass estimate of galaxy clusters. In unrelaxed clusters, the acceleration bias is comparable to the bias due to non-thermal pressure associated with merger-induced turbulent and bulk gas motions. In relaxed clusters, the mean mass bias due to acceleration is small (<3%), but the scatter in the mass bias can be reduced by accounting for gas acceleration. Additionally, this acceleration bias is greater in the outskirts of higher redshift clusters where mergers are more frequent and clusters are accreting more rapidly. Since gas acceleration cannot be observed directly, it introduces an irreducible bias for hydrostatic mass estimates. This acceleration bias places limits on how well we can recover cluster masses from future X-ray and microwave observations. We discuss implications for cluster mass estimates based on X-ray, Sunyaev-Zeldovich effect, and gravitational lensing observations and their impact on cluster cosmology.

Current measurements of the low and high redshift Universe are in tension if we restrict ourselves to the standard six parameter model of flat $\Lambda$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, $H_0$, 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 dataset combination, with statistical errors for clusters, a model with a massive sterile neutrino shows 3.5$\sigma$ evidence for a non-zero 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 1) 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; 2) the results we find are at some variance with analyses that do not include cluster measurements; and 3) some tension remains among the datasets even when new neutrino physics is included.

We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

In the hierarchical structure formation model, clusters of galaxies form through a sequence of mergers and continuous mass accretion, which generate significant random gas motions especially in their outskirts where material is actively accreting. Non-thermal pressure provided by the internal gas motions affects the thermodynamic structure of the X-ray emitting intracluster plasma and introduces biases in the physical interpretation of X-ray and Sunyaev-Zeldovich effect observations. However, we know very little about the nature of gas motions in galaxy clusters. The ASTRO-H X-ray mission, scheduled to launch in 2015, will have a calorimeter capable of measuring gas motions in galaxy clusters at the level of <100 km/s. In this work, we predict the level of merger-induced gas motions expected in the Lambda-CDM model using hydrodynamical simulations of galaxy cluster formation. We show that the gas velocity dispersion is larger in more massive clusters, but exhibits a large scatter. We show that systems with large gas motions are morphologically disturbed, while early forming, relaxed groups show a smaller level of gas motions. By analyzing mock ASTRO-H observations of simulated clusters, we show that such observations can accurately measure the gas velocity dispersion out to the outskirts of nearby relaxed galaxy clusters. ASTRO-H analysis of merging clusters, on the other hand, requires multi-component spectral fitting and enables unique studies of substructures in galaxy clusters by measuring both the peculiar velocities and the velocity dispersion of gas within individual sub-clusters.

We present a detailed investigation of the impact of astrophysical processes on the shape and amplitude of the kinetic Sunyaev-Zel'dovich (kSZ) power spectrum from the post-reionization epoch. This is achieved by constructing a new model of the kSZ power spectrum which we calibrate to the results of hydrodynamic simulations. By construction, our method accounts for all relevant density and velocity modes and so is unaffected by the limited box size of our simulations. We find that radiative cooling and star-formation can reduce the amplitude of the kSZ power spectrum by up to 33%, or 1 uK^2 at ell = 3000. This is driven by a decrease in the mean gas density in groups and clusters due to the conversion of gas into stars. Variations in the redshifts at which helium reionization occurs can effect the amplitude by a similar fraction, while current constraints on cosmological parameters (namely sigma_8) translate to a further +-15% uncertainty on the kSZ power spectrum. We demonstrate how the models presented in this work can be constrained -- reducing the astrophysical uncertainty on the kSZ signal -- by measuring the redshift dependence of the signal via kSZ tomography. Finally, we discuss how the results of this work can help constrain the duration of reionization via measurements of the kinetic SZ signal sourced by inhomogeneous (or patchy) reionization.

In this work, we examine the effects of mergers on the hydrostatic mass estimate of galaxy clusters using high-resolution Eulerian cosmological simulations. We utilize merger trees to isolate the last merger for each cluster in our sample and follow the time evolution of the hydrostatic mass bias as the systems relax. We find that during a merger, a shock propagates outward from the parent cluster, resulting in an overestimate in the hydrostatic mass bias. After the merger, as a cluster relaxes, the bias in hydrostatic mass estimate decreases but remains at a level of -5-10% with 15-20% scatter within r500. We also investigate the post-merger evolution of the pressure support from bulk motions, a dominant cause of this residual mass bias. At r500, the contribution from random motions peaks at 30% of the total pressure during the merger and quickly decays to \sim 10-15% as a cluster relaxes. Additionally, we use a measure of the random motion pressure to correct the hydrostatic mass estimate. We discover that 4 Gyr after mergers, the direct effects of the merger event on the hydrostatic mass bias have become negligible. Thereafter, the mass bias is primarily due to residual bulk motions in the gas which are not accounted for in the hydrostatic equilibrium equation. We present a hydrostatic mass bias correction method that can recover the unbiased cluster mass for relaxed clusters with 9% scatter at r500 and 11% scatter in the outskirts, within r200.

As emphasized by previous studies, proper treatment of the density fluctuation on the fundamental scale of a cosmological simulation volume - the "DC mode" - is critical for accurate modeling of spatial correlations on scales > 10% of simulation box size. We provide further illustration of the effects of the DC mode on the abundance of halos in small boxes and show that it is straightforward to incorporate this mode in cosmological codes that use the "supercomoving" variables. The equations governing evolution of dark matter and baryons recast with these variables are particularly simple and include the expansion factor, and hence the effect of the DC mode, explicitly only in the Poisson equation.

Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

We present high-resolution cosmological hydrodynamic simulations of three galaxy clusters employing a two-temperature model for the intracluster medium. We show that electron temperatures in cluster outskirts are significantly lower than the mean gas temperature, because Coulomb collisions are insufficient to keep electrons and ions in thermal equilibrium. This deviation is larger in more massive and less relaxed systems, ranging from 5% in relaxed clusters to 30% for clusters undergoing major mergers. The presence of non-equilibrium electrons leads to significant suppression of the SZE signal at large cluster-centric radius. The suppression of the electron pressure also leads to an underestimate of the hydrostatic mass. Merger-driven, internal shocks may also generate significant populations of non-equilibrium electrons in the cluster core, leading to a 5% bias on the integrated SZ mass proxy during cluster mergers.

In order to attract and retain excellent researchers and diverse individuals in astrophysics, we recommend action be taken in several key areas impacting young scientists: (1) Maintain balance between large collaborations and individual projects through distribution of funding; encourage public releases of observational and simulation data for use by a broader community. (2) Improve the involvement of women, particularly at leading institutions. (3) Address the critical shortage of child care options and design reasonable profession-wide parental leave policies. (4) Streamline the job application and hiring process. We summarize our reasons for bringing these areas to the attention of the committee, and we suggest several practical steps that can be taken to address them.

Groups and clusters of galaxies occupy a special position in the hierarchy of large-scale cosmic structures because they are the largest and the most massive (from ~10^13 Msun to over 10^15 Msun) objects in the universe that have had time to undergo gravitational collapse. The large masses of clusters imply that their contents have been accreted from regions of ~8-40 comoving Mpc in size and should thus be representative of the mean matter content of the universe. During the next decade sensitive multi-wavelength observations should be able to map the radial distributions of all main mass components (stars, cold, warm, and hot gas and total mass) at z<~ 1 out to the virial radius. At the same time, comparative studies of real and simulated cluster samples sould allow us to use clusters as veritable astrophysical laboratories for studying galaxy formation, as well as testing our theoretical models of structure formation and underlying assumptions about fundamental physics governing the universe.

Cosmological tests based on cluster counts require accurate calibration of the space density of massive halos, but most calibrations to date have ignored complex gas physics associated with halo baryons. We explore the sensitivity of the halo mass function to baryon physics using two pairs of gas-dynamic simulations that are likely to bracket the true behavior. Each pair consists of a baseline model involving only gravity and shock heating, and a refined physics model aimed at reproducing the observed scaling of the hot, intracluster gas phase. One pair consists of billion-particle re-simulations of the original 500 Mpc/h Millennium Simulation of Springel et al. (2005), run with the SPH code Gadget-2 and using a refined physics treatment approximated by preheating (PH) at high redshift. The other pair are high-resolution simulations from the adaptive-mesh refinement code ART, for which the refined treatment includes cooling, star formation, and supernova feedback (CSF). We find that, although the mass functions of the gravity-only (GO) treatments are consistent with the recent calibration of Tinker et al. (2008), both pairs of simulations with refined baryon physics show significant deviations. Relative to the GO case, the masses of ~10^{14} Msun/h halos in the PH and CSF treatments are shifted by averages of -15 \pm 1 percent and +12 \pm 5 percent, respectively. These mass shifts cause ~ 30% deviations in number density relative to the Tinker function, significantly larger than the 5% statistical uncertainty of that calibration.

Numerical studies indicate that uncertainties in the treatment of baryonic physics can affect predictions for shear power spectra at a level that is significant for forthcoming surveys such as DES, SNAP, and LSST. Correspondingly, we show that baryonic effects can significantly bias dark energy parameter measurements. Eliminating such biases by neglecting information in multipoles beyond several hundred leads to weaker parameter constraints by a factor of approximately 2 to 3 compared with using information out to multipoles of several thousand. Fortunately, the same numerical studies that explore the influence of baryons indicate that they primarily affect power spectra by altering halo structure through the relation between halo mass and mean effective halo concentration. We explore the ability of future weak lensing surveys to constrain both the internal structures of halos and the properties of the dark energy simultaneously as a first step toward self calibrating for the physics of baryons. This greatly reduces parameter biases and no parameter constraint is degraded by more than 40% in the case of LSST or 30% in the cases of SNAP or DES. Modest prior knowledge of the halo concentration relation greatly improves even these forecasts. Additionally, we find that these surveys can constrain effective halo concentrations near m~10^14 Msun/h and z~0.2 to better than 10% with shear power spectra alone. These results suggest that inferring dark energy parameters with measurements of shear power spectra can be made robust to baryonic effects and may simultaneously be competitive with other methods to inform models of galaxy formation. (Abridged)

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.

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.