In a recent work, we emphasized that an excess in tri-lepton events plus missing energy observed by the ATLAS experiment at the LHC could be interpreted as a signal of low energy supersymmetry. In such a scenario the lightest neutralino mass is approximately $m_\chi \simeq 60$ GeV and the direct Dark Matter detection cross section is naturally below the current bound. In this work we present simple extensions of this scenario that lead to an explanation of the gamma ray excess at the center of the galaxy observed by Fermi-LAT, as well as the anti-proton excess observed by AMS-02. These extensions include the addition of a small CP violating phase in the neutralino sector or the addition of a light CP-odd Higgs scalar. Our study is of special relevance in view of a recent analysis that casts doubt on the previously accepted preference for mili-second pulsars as the origin of the galactic center excess.

The nearby cluster Abell 1795 is used as a testbed to examine whether hot gas in cluster galaxies is stripped by the ram pressure of the intracluster medium (ICM). The expected X-ray emission in and around Abell 1795 galaxies is likely dominated by the ICM, low-mass X-ray binaries, active galactic nuclei, and hot gas halos. In order to constrain these components, we use archival Chandra X-ray Observatory and Sloan Digital Sky Survey (SDSS) observations of Abell 1795 and identify 58 massive (M_star>10^10 M_sun) spectroscopic cluster members within 5 arcmin of the Chandra optical axis. X-ray images at 0.5-1.5 keV and 4-8 keV were created for each cluster member and then stacked into two clustercentric radius bins: inner (0.25<R/R_500<1) and outer (1<R/R_500<2.5). Surface brightness profiles of inner and outer cluster members are fit using Markov chain Monte Carlo sampling in order to generate model parameters and measure the 0.5-1.5 keV luminosities of each model component. Leveraging effective total Chandra exposure times of 3.4 and 1.7 Msec for inner and outer cluster members, respectively, we report the detection of hot gas halos, in a statistical sense, around outer cluster members. Outer members have 0.5-1.5 keV hot halo luminosities (L_X = 8.1(-3.5/+5)x10^39 erg/s) that are six times larger than the upper limit for inner cluster members (L_X < 1.3x10^39 erg/s). This result suggests that the ICM is removing hot gas from the halos of Abell 1795 members as they fall into the cluster.

We analyze a low energy effective model of Dark Matter in which the thermal relic density is provided by a singlet Majorana fermion which interacts with the Higgs fields via higher dimensional operators. Direct detection signatures may be reduced if blind spot solutions exist, which naturally appear in models with extended Higgs sectors. Explicit mass terms for the Majorana fermion can be forbidden by a $Z_3$ symmetry, which in addition leads to a reduction of the number of higher dimensional operators. Moreover, a weak scale mass for the Majorana fermion is naturally obtained from the vacuum expectation value of a scalar singlet field. The proper relic density may be obtained by the $s$-channel interchange of Higgs and gauge bosons, with the longitudinal mode of the $Z$ boson (the neutral Goldstone mode) playing a relevant role in the annihilation process. This model shares many properties with the Next-to-Minimal Supersymmetric extension of the Standard Model (NMSSM) with light singlinos and heavy scalar and gauge superpartners. In order to test the validity of the low energy effective field theory, we compare its predictions with those of the ultraviolet complete NMSSM. Extending our framework to include $Z_3$ neutral Majorana fermions, analogous to the bino in the NMSSM, we find the appearance of a new bino-singlino well tempered Dark Matter region.

We explore 7.5 billion years of evolution in the star formation activity of massive ($M_{\star}>10^{10.1}\,M_{\odot}$) cluster galaxies using a sample of 25 clusters over $0.15<z<1$ from the Cluster Lensing And Supernova survey with Hubble and 11 clusters over $1<z<1.5$ from the IRAC Shallow Cluster Survey. Galaxy morphologies are determined visually using high-resolution Hubble Space Telescope images. Using the spectral energy distribution fitting code CIGALE, we measure star formation rates, stellar masses, and 4000 \AA\ break strengths. The latter are used to separate quiescent and star-forming galaxies (SFGs). From $z\sim1.3$ to $z\sim0.2$, the specific star formation rate (sSFR) of cluster SFGs and quiescent galaxies decreases by factors of three and four, respectively. Over the same redshift range, the sSFR of the entire cluster population declines by a factor of 11, from $0.48\pm0.06\;\mathrm{Gyr}^{-1}$ to $0.043\pm0.009\;\mathrm{Gyr}^{-1}$. This strong overall sSFR evolution is driven by the growth of the quiescent population over time; the fraction of quiescent cluster galaxies increases from $28^{+8}_{-19}\%$ to $88^{+5}_{-4}\%$ over $z\sim1.3\rightarrow0.2$. The majority of the growth occurs at $z\gtrsim0.9$, where the quiescent fraction increases by 0.41. While the sSFR of the majority of star-forming cluster galaxies is at the level of the field, a small subset of cluster SFGs have low field-relative star formation activity, suggestive of long-timescale quenching. The large increase in the fraction of quiescent galaxies above $z\sim0.9$, coupled with the field-level sSFRs of cluster SFGs, suggests that higher redshift cluster galaxies are likely being quenched quickly. Assessing those timescales will require more accurate stellar population ages and star formation histories.

We present accurate measurements of the linear, quadratic, and cubic local bias of dark matter halos, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength overdensity. This can be seen as an exact implementation of the peak-background split argument. We compare the results with the linear and quadratic bias measured from the halo-matter power spectrum and bispectrum, and find good agreement. On the other hand, the standard peak-background split applied to the Sheth & Tormen (1999) and Tinker et al. (2008) halo mass functions matches the measured linear bias parameter only at the level of 10%. The prediction from the excursion set-peaks approach performs much better, which can be attributed to the stochastic moving barrier employed in the excursion set-peaks prediction. We also provide convenient fitting formulas for the nonlinear bias parameters $b_2(b_1)$ and $b_3(b_1)$, which work well over a range of redshifts.

We report on the first measurement of the three-point function with the position-dependent correlation function from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 10 CMASS sample. This new observable measures the correlation between two-point functions of galaxy pairs within different subvolumes, $\hat{\xi}({\rm r},{\rm r}_L)$, where ${\rm r}_L$ is the location of a subvolume, and the corresponding mean overdensities, $\bar{\delta}({\rm r}_L)$. This correlation, which we call the "integrated three-point function", $i\zeta(r)=\langle\hat{\xi}({\rm r},{\rm r}_L)\bar{\delta}({\rm r}_L)\rangle$, measures a three-point function of two short- and one long-wavelength modes, and is generated by nonlinear gravitational evolution and possibly also by the physics of inflation. The $i\zeta(r)$ measured from the BOSS data lies within the scatter of those from the mock galaxy catalogs in redshift space, yielding a ten-percent-level determination of the amplitude of $i\zeta(r)$. The tree-level perturbation theory in redshift space predicts how this amplitude depends on the linear and quadratic nonlinear galaxy bias parameters ($b_1$ and $b_2$), as well as on the amplitude and linear growth rate of matter fluctuations ($\sigma_8$ and $f$). Combining $i\zeta(r)$ with the constraints on $b_1\sigma_8$ and $f\sigma_8$ from the global two-point correlation function and that on $\sigma_8$ from the weak lensing signal of BOSS galaxies, we measure $b_2=0.41\pm0.41$ (68% C.L.) assuming standard perturbation theory at the tree level and the local bias model.

We examine the cosmological implications of the measurements of the linear growth rate of cosmological structure obtained in a companion paper from the power spectrum and bispectrum monopoles of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey Data, Release 11, CMASS galaxies. This measurement was of $f^{0.43}\sigma_8$, where $\sigma_8$ is the amplitude of dark matter density fluctuations, and $f$ is the linear growth rate, at the effective redshift of the survey, $z_{\rm eff}=0.57$. In conjunction with Cosmic Microwave Background (CMB) data, interesting constraints can be placed on models with non-standard neutrino properties and models where gravity deviates from general relativity on cosmological scales. In particular, the sum of the masses of the three species of the neutrinos is constrained to $m_\nu<0.49\,{\rm eV}$ (at 95\% confidence level) when the $f^{0.43}\sigma_8$ measurement is combined with state-of-the-art CMB measurements. Allowing the effective number of neutrinos to vary as a free parameter does not significantly change these results. When we combine the measurement of $f^{0.43}\sigma_8$ with the complementary measurement of $f\sigma_8$ from the monopole and quadrupole of the two-point correlation function we are able to obtain an independent measurements of $f$ and $\sigma_8$. We obtain $f=0.63\pm0.16$ and $\sigma_8=0.710\pm0.086$ (68\% confidence level). This is the first time when these parameters have been able to be measured independently using the redshift-space power spectrum and bispectrum measurements from galaxy clustering data only.

We analyse the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 sample, which consists of $690 827$ galaxies in the redshift range $0.43 < z < 0.70$ and has a sky coverage of $8 498$ deg$^2$ corresponding to an effective volume of $\sim6\,\rm{Gpc}^3$. We fit the Fourier space statistics, the power spectrum and bispectrum monopoles to measure the linear and quadratic bias parameters, $b_1$ and $b_2$, for a non-linear non-local bias model, the growth of structure parameter $f$ and the amplitude of dark matter density fluctuations parametrised by $\sigma_8$. We obtain $b_1(z_{\rm eff})^{1.40}\sigma_8(z_{\rm eff})=1.672\pm 0.060$ and $b_2^{0.30}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.579\pm0.082$ at the effective redshift of the survey, $z_{\rm eff}=0.57$. The main cosmological result is the constraint on the combination $f^{0.43}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.582\pm0.084$, which is complementary to $f\sigma_8$ constraints obtained from 2-point redshift space distortion analyses. A less conservative analysis yields $f^{0.43}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.584\pm0.051$. We ensure that our result is robust by performing detailed systematic tests using a large suite of survey galaxy mock catalogs and N-body simulations. The constraints on $f^{0.43}\sigma_8$ are useful for setting additional constrains on neutrino mass, gravity, curvature as well as the number of neutrino species from galaxy surveys analyses (as presented in a companion paper).

We show that in a certain, angle-averaged squeezed limit, the $N$-point function of matter is related to the response of the matter power spectrum to a long-wavelength density perturbation, $P^{-1}d^nP(k|\delta_L)/d\delta_L^n|_{\delta_L=0}$, with $n=N-2$. By performing N-body simulations with a homogeneous overdensity superimposed on a flat Friedmann-Robertson-Lema\^itre-Walker (FRLW) universe using the \emph{separate universe} approach, we obtain measurements of the nonlinear matter power spectrum response up to $n=3$, which is equivalent to measuring the fully nonlinear matter $3-$ to $5-$point function in this squeezed limit. The sub-percent to few percent accuracy of those measurements is unprecedented. We then test the hypothesis that nonlinear $N$-point functions at a given time are a function of the linear power spectrum at that time, which is predicted by standard perturbation theory (SPT) and its variants that are based on the ideal pressureless fluid equations. Specifically, we compare the responses computed from the separate universe simulations and simulations with a rescaled initial (linear) power spectrum amplitude. We find discrepancies of 10\% at $k\simeq 0.2 - 0.5 \,h\,{\rm Mpc}^{-1}$ for $5-$ to $3-$point functions at $z=0$. The discrepancy occurs at higher wavenumbers at $z=2$. Thus, SPT and its variants, carried out to arbitrarily high order, are guaranteed to fail to describe matter $N$-point functions ($N>2$) around that scale.

We measure the star formation rates (SFRs) of massive ($M_{\star}>10^{10.1}M_{\odot}$) early-type galaxies (ETGs) in a sample of 11 high-redshift ($1.0 < z < 1.5$) galaxy clusters drawn from the IRAC Shallow Cluster Survey (ISCS). We identify ETGs visually from Hubble Space Telescope imaging and select likely cluster members as having either an appropriate spectroscopic redshift or red sequence color. Mid-infrared SFRs are measured using Spitzer 24 $\mu$m data for isolated cluster galaxies for which contamination by neighbors, and active galactic nuclei, can be ruled out. Cluster ETGs show enhanced specific star formation rates (sSFRs) compared to cluster galaxies in the local Universe, but have sSFRs more than four times lower than that of field ETGs at $1 < z < 1.5$. Relative to the late-type cluster population, isolated ETGs show substantially quenched mean SFRs, yet still contribute 12% of the overall star formation activity measured in $1 < z < 1.5$ clusters. We find that new ETGs are likely being formed in ISCS clusters; the fraction of cluster galaxies identified as ETGs increases from 34% to 56% from $z \sim 1.5 \rightarrow 1.25$. While the fraction of cluster ETGs that are highly star-forming ($\textrm{SFR}\geq26\ M_{\odot}$ yr$^{-1}$) drops from 27% to 10% over the same period, their sSFRs are roughly constant. All these factors taken together suggest that, particularly at $z\gtrsim1.25$, the events that created these distant cluster ETGs$-$likely mergers, at least among the most massive$-$were both recent and gas-rich.

We present a phenomenological modification of the standard perturbation theory prediction for the bispectrum in redshift space that allows us to extend the model to mildly non-linear scales over a wide range of redshifts, $z\leq1.5$. We find that we can describe the bispectrum of dark matter particles with $\sim5%$ accuracy for $k_i\lesssim0.10\,h/{\rm Mpc}$ at $z=0$, for $k_i\lesssim0.15\,h/{\rm Mpc}$ at $z=0.5$, for $k_i\lesssim0.17\,h/{\rm Mpc}$ at $z=1.0$ and for $k_i\lesssim0.20\,h/{\rm Mpc}$ at $z=1.5$. We also test that the fitting formula is able to describe with similar accuracy the bispectrum of cosmologies with different $\Omega_m$, in the range $0.2\lesssim \Omega_m \lesssim 0.4$, and consequently with different values of the logarithmic grow rate $f$ at $z=0$, $0.4\lesssim f(z=0) \lesssim 0.6$. We apply this new formula to recover the bias parameters, $f$ and $\sigma_8$, by combining the redshift space power spectrum monopole and quadrupole with the bispectrum monopole for both dark matter particles and haloes. We find that the combination of these three statistics can break the degeneracy between $b_1$, $f$ and $\sigma_8$. For dark matter particles the new model can be used to recover $f$ and $\sigma_8$ with $\sim1%$ accuracy. For dark matter haloes we find that $f$ and $\sigma_8$ present larger systematic shifts, $\sim10%$. The systematic offsets arise because of limitations in the modelling of the interplay between bias and redshift space distortions, and represent a limitation as the statistical errors of forthcoming surveys reach this level. Conveniently, we find that these residual systematics are mitigated for combinations of parameters. The improvement on the modelling of the bispectrum presented in this paper will be useful for extracting information from current and future galaxy surveys. [abridged]

In this article we compare the halo mass function predicted by the excursion set theory with a drifting diffusive barrier against the results of N-body simulations for several cosmological models. This includes the standard LCDM case for a large range of halo masses, models with different types of primordial non-Gaussianity, and the Ratra-Peebles quintessence model of Dark Energy. We show that in all those cosmological scenarios, the abundance of dark matter halos can be described by a drifting diffusive barrier, where the two parameters describing the barrier have physical content. In the case of the Gaussian LCDM, the statistics are precise enough to actually predict those parameters at different redshifts from the initial conditions. Furthermore, we found that the stochasticity in the barrier is nonnegligible making the simple deterministic spherical collapse model a bad approximation even at very high halo masses. We also show that using the standard excursion set approach with a barrier inspired by peak patches leads to inconsistent predictions of the halo mass function.

The large-scale statistics of observables such as the galaxy density are chiefly determined by their dependence on the local coarse-grained matter density. This dependence can be measured directly and efficiently in N-body simulations by using the fact that a uniform density perturbation with respect to some fiducial background cosmology is equivalent to modifying the background and including curvature, i.e., by simulating a "separate universe". We derive this mapping to fully non-linear order, and provide a step-by-step description of how to perform and analyse the separate universe simulations. This technique can be applied to a wide range of observables. As an example, we calculate the response of the non-linear matter power spectrum to long-wavelength density perturbations, which corresponds to the angle-averaged squeezed limit of the matter bispectrum and higher $n$-point functions. Using only a modest simulation volume, we obtain results with percent-level precision over a wide range of scales.

We investigate the cosmological implications of the latest growth of structure measurement from the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 with particular focus on the sum of the neutrino masses, $\sum m_{\nu}$. We examine the robustness of the cosmological constraints from the Baryon Acoustic Oscillation (BAO) scale, the Alcock-Paczynski effect and redshift-space distortions ($D_V/r_s$, $F_{\rm AP}$, $f\sigma_8$) of \citet{Beutler:2013b}, when introducing a neutrino mass in the power spectrum template. We then discuss how the neutrino mass relaxes discrepancies between the Cosmic Microwave Background (CMB) and other low-redshift measurements within $\Lambda$CDM. Combining our cosmological constraints with WMAP9 yields $\sum m_{\nu} = 0.36\pm0.14\,$eV ($68\%$ c.l.), which represents a $2.6\sigma$ preference for non-zero neutrino mass. The significance can be increased to $3.3\sigma$ when including weak lensing results and other BAO constraints, yielding $\sum m_{\nu} = 0.35\pm0.10\,$eV ($68\%$ c.l.). However, combining CMASS with Planck data reduces the preference for neutrino mass to $\sim 2\sigma$. When removing the CMB lensing effect in the Planck temperature power spectrum (by marginalising over $A_{\rm L}$), we see shifts of $\sim 1\sigma$ in $\sigma_8$ and $\Omega_m$, which have a significant effect on the neutrino mass constraints. In case of CMASS plus Planck without the $A_{\rm L}$-lensing signal, we find a preference for a neutrino mass of $\sum m_{\nu} = 0.34\pm 0.14\,$eV ($68\%$ c.l.), in excellent agreement with the WMAP9+CMASS value. The constraint can be tightened to $3.4\sigma$ yielding $\sum m_{\nu} = 0.36\pm 0.10\,$eV ($68\%$ c.l.) when weak lensing data and other BAO constraints are included.

We use the Sloan Digital Sky Survey to investigate the properties of massive elliptical galaxies in the local Universe (z\leq0.08) that have unusually blue optical colors. Through careful inspection, we distinguish elliptical from non-elliptical morphologies among a large sample of similarly blue galaxies with high central light concentrations (c_r\geq2.6). These blue ellipticals comprise 3.7 per cent of all c_r\geq2.6 galaxies with stellar masses between 10^10 and 10^11 h^{-2} {\rm M}_{\sun}. Using published fiber spectra diagnostics, we identify a unique subset of 172 non-star-forming ellipticals with distinctly blue urz colors and young (< 3 Gyr) light-weighted stellar ages. These recently quenched ellipticals (RQEs) have a number density of 2.7-4.7\times 10^{-5}\,h^3\,{\rm Mpc}^{-3} and sufficient numbers above 2.5\times10^{10} h^{-2} {\rm M}_{\sun} to account for more than half of the expected quiescent growth at late cosmic time assuming this phase lasts 0.5 Gyr. RQEs have properties that are consistent with a recent merger origin (i.e., they are strong `first-generation' elliptical candidates), yet few involved a starburst strong enough to produce an E+A signature. The preferred environment of RQEs (90 per cent reside at the centers of < 3\times 10^{12}\,h^{-1}{\rm M}_{\sun} groups) agrees well with the `small group scale' predicted for maximally efficient spiral merging onto their halo center and rules out satellite-specific quenching processes. The high incidence of Seyfert and LINER activity in RQEs and their plausible descendents may heat the atmospheres of small host halos sufficiently to maintain quenching.

The influence of large-scale density fluctuations on structure formation on small scales is described by the three-point correlation function (bispectrum) in the so-called "squeezed configurations," in which one wavenumber, say $k_3$, is much smaller than the other two, i.e., $k_3\ll k_1\approx k_2$. This bispectrum is generated by non-linear gravitational evolution and possibly also by inflationary physics. In this paper, we use this fact to show that the bispectrum in the squeezed configurations can be measured without employing three-point function estimators. Specifically, we use the "position-dependent power spectrum," i.e., the power spectrum measured in smaller subvolumes of the survey (or simulation box), and correlate it with the mean overdensity of the corresponding subvolume. This correlation directly measures an integral of the bispectrum dominated by the squeezed configurations. Measuring this correlation is only slightly more complex than measuring the power spectrum itself, and sidesteps the considerable complexity of the full bispectrum estimation. We use cosmological $N$-body simulations of collisionless particles with Gaussian initial conditions to show that the measured correlation between the position-dependent power spectrum and the long-wavelength overdensity agrees with the theoretical expectation. The position-dependent power spectrum thus provides a new, efficient, and promising way to measure the squeezed-limit bispectrum from large-scale structure observations such as galaxy redshift surveys.

In this paper, we study an extension of the standard model with a vector-like generation of leptons. This model provides a viable dark matter candidate and a possibility to enhance the Higgs decay rate into a pair of photons. We evaluate constraints from electroweak precision tests and from vacuum stability, and find that the latter provide an upper limit on the lepton Yukawa couplings. A large enhancement of the Higgs di-photon rate can therefore only be obtained when the mass of the lightest charged lepton is close to the LEP limit. The relic density constraint suggests a co-annihilation scenario with a small mass difference between the lightest charged and neutral leptons, which also weakens the LEP limit on the lightest charged lepton mass and allows for larger Higgs di-photon decay rates. Cross sections for direct detection of the dark matter candidate are calculated, and prospects for detecting the new particles at the LHC are discussed briefly.

Many of the most exciting questions in astrophysics and cosmology, including the majority of observational probes of dark energy, rely on an understanding of the nonlinear regime of structure formation. In order to fully exploit the information available from this regime and to extract cosmological constraints, accurate theoretical predictions are needed. Currently such predictions can only be obtained from costly, precision numerical simulations. This paper is the third in a series aimed at constructing an accurate calibration of the nonlinear mass power spectrum on Mpc scales for a wide range of currently viable cosmological models, including dark energy. The first two papers addressed the numerical challenges, and the scheme by which an interpolator was built from a carefully chosen set of cosmological models. In this paper we introduce the "Coyote Univers"' simulation suite which comprises nearly 1,000 N-body simulations at different force and mass resolutions, spanning 38 wCDM cosmologies. This large simulation suite enables us to construct a prediction scheme, or emulator, for the nonlinear matter power spectrum accurate at the percent level out to k~1 h/Mpc. We describe the construction of the emulator, explain the tests performed to ensure its accuracy, and discuss how the central ideas may be extended to a wider range of cosmological models and applications. A power spectrum emulator code is released publicly as part of this paper.

We investigate the accuracy of Eulerian perturbation theory for describing the matter and galaxy power spectra in real and redshift space in light of future observational probes for precision cosmology. Comparing the analytical results with a large suite of N-body simulations (160 independent boxes of 13.8 (Gpc/h)^3 volume each, which are publicly available), we find that re-summing terms in the standard perturbative approach predicts the real-space matter power spectrum with an accuracy of < 2% for k < 0.20 h/Mpc at redshifts z < 1.5. This is obtained following the widespread technique of writing the resummed propagator in terms of 1-loop contributions. We show that the accuracy of this scheme increases by considering higher-order terms in the resummed propagator. By combining resummed perturbation theories with several models for the mappings from real to redshift space discussed in the literature, the multipoles of the dark-matter power spectrum can be described with sub-percent deviations from N-body results for k < 0.15h/Mpc at z < 1. As a consequence, the logarithmic growth rate, f, can be recovered with sub-percent accuracy on these scales. Extending the models to massive dark-matter haloes in redshift space, our results describe the monopole term from N-body data within 2% accuracy for scales k < 0.15 h/Mpc at z < 0.5; here f can be recovered within < 5% when the halo bias is known. We conclude that these techniques are suitable to extract cosmological information from future galaxy surveys.

We implement Slow Roll Reconstruction -- an optimal solution to the inverse problem for inflationary cosmology -- within ModeCode, a publicly available solver for the inflationary dynamics. We obtain up-to-date constraints on the reconstructed inflationary potential, derived from the WMAP 7-year dataset and South Pole Telescope observations, combined with large scale structure data derived from SDSS Data Release 7. Using ModeCode in conjunction with the MultiNest sampler, we compute Bayesian evidence for the reconstructed potential at each order in the truncated slow roll hierarchy. We find that the data are well-described by the first two slow roll parameters, \epsilon and \eta, and that there is no need to include a nontrivial \xi parameter.

We present measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5\sigma in both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7\sigma. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This {\it beat coupling} effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated {\it average} density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this {\it local average effect} and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only $\sim 10 %$ of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations. For the variances, we find excellent agreement with the analytic expressions for $k < 0.2 h$Mpc$^{-1}$ at $z=0.5$, while the correlation coefficients agree to beyond $k=0.4 h$Mpc$^{-1}$. As expected, the range of agreement increases towards higher redshift and decreases slightly towards $z=0$. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix. [abridged]

We have performed cosmological N-body simulations which include the effect of the masses of the individual neutrino species. The simulations were aimed at studying the effect of different neutrino hierarchies on the matter power spectrum. Compared to the linear theory predictions, we find that non-linearities enhance the effect of hierarchy on the matter power spectrum at mildly non-linear scales. The difference between the different hierarchies is about 0.5% for a sum of neutrino masses of 0.1eV. Albeit this is a small effect, it is potentially measurable from upcoming surveys. In combination with neutrinoless double-beta decay experiments, this opens up the possibility of using the sky to determine if neutrinos are Majorana or Dirac fermions.

In this paper we present an improved fitting formula for the dark matter bispectrum motivated by the previous phenomenological approach of Scoccimarro & Couchman (2001). We use a set of LCDM simulations to calibrate the fitting parameters in the k-range of 0.03 h/Mpc<k<0.4 h/Mpc and in the redshift range of 0<z<1.5. This new proposed fit describes well the BAO-features although it was not designed to. The deviation between the simulations output and our analytic prediction is typically less than 5% and in the worst case is never above 10%. We envision that this new analytic fitting formula will be very useful in providing reliable predictions for the non-linear dark matter bispectrum for LCDM models.

We report new mid-eclipse times of the two close binaries NSVS14256825 and HS0705+6700, harboring an sdB primary and a low-mass main-sequence secondary. Both objects display clear variations in the measured orbital period, which can be explained by the action of a third object orbiting the binary. If this interpretation is correct, the third object in NSVS14256825 is a giant planet with a mass of roughly 12 M_Jup. For HS0705+6700, we provide evidence that strengthens the case for the suggested periodic nature of the eclipse time variation and reduces the uncertainties in the parameters of the brown dwarf implied by that model. The derived period is 8.4 yr and the mass is 31 M_Jup, if the orbit is coplanar with the binary. This research is part of the PlanetFinders project, an ongoing collaboration between professional astronomers and student groups at high schools.

We present N-body simulations for generic non-Gaussian initial conditions with the aim of exploring and modelling the scale-dependent halo bias. This effect is evident at very large scales requiring large simulation boxes. In addition, the previously available prescription to implement generic non-Gaussian initial conditions has been improved to keep under control higher-order terms which were spoiling the power spectrum on large scales. We pay particular attention to the differences between physical, inflation-motivated primordial bispectra and their factorizable templates, and to the operational definition of the non-Gaussian halo bias (which has both a scale-dependent and an approximately scale-independent contributions). We find that analytic predictions for both the non-Gaussian halo mass function and halo bias work well once a calibration factor (which was introduced before) is calibrated on simulations. The halo bias remains therefore an extremely promising tool to probe primordial non-Gaussianity and thus to give insights into the physical mechanism that generated the primordial perturbations. The simulation outputs and tables of the analytic predictions will be made publicly available via the non-Gaussian comparison project web site http://icc.ub.edu/~liciaverde/NGSCP.html

We study a limit of the nearly-Peccei-Quinn-symmetric Next-to-Minimal Supersymmetric Standard Model possessing novel Higgs and dark matter (DM) properties. In this scenario, there naturally co-exist three light singlet-like particles: a scalar, a pseudoscalar, and a singlino-like DM candidate, all with masses of order 0.1-10 GeV. The decay of a Standard Model-like Higgs boson to pairs of the light scalars or pseudoscalars is generically suppressed, avoiding constraints from collider searches for these channels. For a certain parameter window annihilation into the light pseudoscalar and exchange of the light scalar with nucleons allow the singlino to achieve the correct relic density and a large direct detection cross section consistent with the CoGeNT and DAMA/LIBRA preferred region simultaneously. This parameter space is consistent with experimental constraints from LEP, the Tevatron, and Upsilon- and flavor physics.

This paper describes the Hubble Space Telescope imaging data products and data reduction procedures for the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). This survey is designed to document the evolution of galaxies and black holes at $z\sim1.5-8$, and to study Type Ia SNe beyond $z>1.5$. Five premier multi-wavelength sky regions are selected, each with extensive multiwavelength observations. The primary CANDELS data consist of imaging obtained in the Wide Field Camera 3 / infrared channel (WFC3/IR) and UVIS channel, along with the Advanced Camera for Surveys (ACS). The CANDELS/Deep survey covers \sim125 square arcminutes within GOODS-N and GOODS-S, while the remainder consists of the CANDELS/Wide survey, achieving a total of \sim800 square arcminutes across GOODS and three additional fields (EGS, COSMOS, and UDS). We summarize the observational aspects of the survey as motivated by the scientific goals and present a detailed description of the data reduction procedures and products from the survey. Our data reduction methods utilize the most up to date calibration files and image combination procedures. We have paid special attention to correcting a range of instrumental effects, including CTE degradation for ACS, removal of electronic bias-striping present in ACS data after SM4, and persistence effects and other artifacts in WFC3/IR. For each field, we release mosaics for individual epochs and eventual mosaics containing data from all epochs combined, to facilitate photometric variability studies and the deepest possible photometry. A more detailed overview of the science goals and observational design of the survey are presented in a companion paper.

The mass distribution of halos, as specified by the halo mass function, is a key input for several cosmological probes. The sizes of $N$-body simulations are now such that, for the most part, results need no longer be statistics-limited, but are still subject to various systematic uncertainties. We investigate and discuss some of the reasons for these differences. Quantifying error sources and compensating for them as appropriate, we carry out a high-statistics study of dark matter halos from 67 $N$-body simulations to investigate the mass function and its evolution for a reference $\Lambda$CDM cosmology and for a set of $w$CDM cosmologies. For the reference $\Lambda$CDM cosmology (close to WMAP5), we quantify the breaking of universality in the form of the mass function as a function of redshift, finding an evolution of as much as 10% away from the universal form between redshifts $z=0$ and $z=2$. For cosmologies very close to this reference we provide a fitting formula to our results for the (evolving) $\Lambda$CDM mass function over a mass range of $6\cdot 10^{11}-3\cdot 10^{15}$ M$_{\odot}$ to an estimated accuracy of about 2%. The set of $w$CDM cosmologies is taken from the Coyote Universe simulation suite. The mass functions from this suite (which includes a $\Lambda$CDM cosmology and others with $w\simeq-1$) are described by the fitting formula for the reference $\Lambda$CDM case at an accuracy level of 10%, but with clear systematic deviations. We argue that, as a consequence, fitting formulae based on a universal form for the mass function may have limited utility in high precision cosmological applications.

Comparing clustering of differently biased tracers of the dark matter distribution offers the opportunity to reduce the cosmic variance error in the measurement of certain cosmological parameters. We develop a formalism that includes bias non-linearities and stochasticity. Our formalism is general enough that can be used to optimise survey design and tracers selection and optimally split (or combine) tracers to minimise the error on the cosmologically interesting quantities. Our approach generalises the one presented by McDonald & Seljak (2009) of circumventing sample variance in the measurement of $f\equiv d \ln D/d\ln a$. We analyse how the bias, the noise, the non-linearity and stochasticity affect the measurements of $Df$ and explore in which signal-to-noise regime it is significantly advantageous to split a galaxy sample in two differently-biased tracers. We use N-body simulations to find realistic values for the parameters describing the bias properties of dark matter haloes of different masses and their number density. We find that, even if dark matter haloes could be used as tracers and selected in an idealised way, for realistic haloes, the sample variance limit can be reduced only by up to a factor $\sigma_{2tr}/\sigma_{1tr}\simeq 0.6$. This would still correspond to the gain from a three times larger survey volume if the two tracers were not to be split. Before any practical application one should bear in mind that these findings apply to dark matter haloes as tracers, while realistic surveys would select galaxies: the galaxy-host halo relation is likely to introduce extra stochasticity, which may reduce the gain further.

Near-future cosmological observations targeted at investigations of dark energy pose stringent requirements on the accuracy of theoretical predictions for the clustering of matter. Currently, N-body simulations comprise the only viable approach to this problem. In this paper we demonstrate that N-body simulations can indeed be sufficiently controlled to fulfill these requirements for the needs of ongoing and near-future weak lensing surveys. By performing a large suite of cosmological simulation comparison and convergence tests we show that results for the nonlinear matter power spectrum can be obtained at 1% accuracy out to k~1 h/Mpc. The key components of these high accuracy simulations are: precise initial conditions, very large simulation volumes, sufficient mass resolution, and accurate time stepping. This paper is the first in a series of three, with the final aim to provide a high-accuracy prediction scheme for the nonlinear matter power spectrum.

The power spectrum of density fluctuations is a foundational source of cosmological information. Precision cosmological probes targeted primarily at investigations of dark energy require accurate theoretical determinations of the power spectrum in the nonlinear regime. To exploit the observational power of future cosmological surveys, accuracy demands on the theory are at the one percent level or better. Numerical simulations are currently the only way to produce sufficiently error-controlled predictions for the power spectrum. The very high computational cost of (precision) N-body simulations is a major obstacle to obtaining predictions in the nonlinear regime, while scanning over cosmological parameters. Near-future observations, however, are likely to provide a meaningful constraint only on constant dark energy equation of state 'wCDM' cosmologies. In this paper we demonstrate that a limited set of only 37 cosmological models -- the "Coyote Universe" suite -- can be used to predict the nonlinear matter power spectrum at the required accuracy over a prior parameter range set by cosmic microwave background observations. This paper is the second in a series of three, with the final aim to provide a high-accuracy prediction scheme for the nonlinear matter power spectrum for wCDM cosmologies.

We address the issue of setting up generic non-Gaussian initial conditions for N-body simulations. We consider inflationary-motivated primordial non-Gaussianity where the perturbations in the Bardeen potential are given by a dominant Gaussian part plus a non-Gaussian part specified by its bispectrum. The approach we explore here is suitable for any bispectrum, i.e. it does not have to be of the so-called separable or factorizable form. The procedure of generating a non-Gaussian field with a given bispectrum (and a given power spectrum for the Gaussian component) is not univocal, and care must be taken so that higher-order corrections do not leave a too large signature on the power spectrum. This is so far a limiting factor of our approach. We then run N-body simulations for the most popular inflationary-motivated non-Gaussian shapes. The halo mass function and the non-linear power spectrum agree with theoretical analytical approximations proposed in the literature, even if they were so far developed and tested only for a particular shape (the local one). We plan to make the simulations outputs available to the community via the non-Gaussian simulations comparison project web site http://icc.ub.edu/~liciaverde/NGSCP.html.

We systematically study the effects of varying the starting redshift z_i for cosmological simulations in the highly non-linear regime. Our primary focus lies with the (individual) properties of dark matter halos -- namely the mass, spin, triaxiality, and concentration -- where we find that even substantial variations in z_i leave only a small imprint, at least for the probed mass range M \in [10^{10}, 10^{13}] Msun/h and when investigated at redshift z=0. We further compare simulations started by using the standard Zel'dovich approximation to runs based upon initial conditions produced with second order Lagrangian perturbation theory. Here we observe the same phenomenon, i.e. that differences in the studied (internal) properties of dark matter haloes are practically undetectable. These findings are (for the probed mass range) in agreement with other work in the literature. We therefore conclude that the commonly used technique for setting up cosmological simulations leads to stable results at redshift z=0 for the mass, the spin parameter, the triaxiality, and the concentration of dark matter haloes.

The measurement of the scale of the baryon acoustic oscillations (BAO) in the galaxy power spectrum as a function of redshift is a promising method to constrain the equation-of-state parameter of the dark energy w. To measure the scale of the BAO precisely, a substantial volume of space must be surveyed. We test whether light-cone effects are important and whether the scaling relations used to compensate for an incorrect reference cosmology are in this case sufficiently accurate. We investigate the degeneracies in the cosmological parameters and the benefits of using the two-dimensional anisotropic power spectrum. Finally, we estimate the uncertainty with which w can be measured by proposed surveys at redshifts of about z=3 and z=1, respectively. In the simulated survey we find that light-cone effects are small and that the simple scaling relations used to correct for the cosmological distortion work fairly well even for large survey volumes. The analysis of the two-dimensional anisotropic power spectra enables an independent determination to be made of the apparent scale of the BAO, perpendicular and parallel to the line of sight. This is essential for two-parameter w-models, such as the redshift-dependent dark energy model w=w_0+(1-a)w_a. Using Planck priors for the matter and baryon density and Delta(H_0)=5% for the Hubble constant, we estimate that the BAO measurements of future surveys around z=3 and z=1 will be able to constrain, independently of other cosmological probes, a constant w to ~ 12% and ~ 11% (68% c.l.), respectively.

Models of Universal Extra Dimensions (UED) at the TeV scale lead to the presence of Kaluza Klein (KK) excitations of the ordinary fermions and bosons of the Standard Model that may be observed at hadron and lepton colliders. A conserved discrete symmetry, KK-parity, ensures the stability of the lightest KK particle (LKP), which, if neutral, becomes a good dark matter particle. It has been recently shown that for a certain range of masses of the LKP a relic density consistent with the experimentally observed one may be obtained. These works, however, ignore the impact of KK graviton production at early times. Whether the G^1 is the LKP or not, the G^n tower thus produced can decay to the LKP, and depending on the reheating temperature, may lead to a modification of the relic density. In this article, we show that this effect may lead to a relevant modification of the range of KK masses consistent with the observed relic density. Additionally, if evidence for UED is observed experimentally, we find a stringent upper limit on the reheating temperature depending on the mass of the LKP observed.

The MareNostrum Universe is one of the biggest SPH cosmological simulations done so far. It contains more than 2 billion particles (2 times 1024^3) in a 500 Mpc/h cubic volume. This simulation has been performed on the MareNostrum supercomputer at the Barcelona Supercomputer Center. We have obtained more than 0.5 million halos with masses greater than a typical Milky Way galaxy halo. We report results about the halo mass function, the shapes of dark matter and gas distributions in halos, the baryonic fraction in galaxy clusters and groups, baryon oscillations in the dark matter and the halo power spectra as well as the distribution and evolution of the gas fraction at large scales.

We study the generation of primordial magnetic fields, coherent over cosmologically interesting scales, by gravitational creation of charged scalar particles during the reheating period. We show that magnetic fields consistent with those detected by observation may obtained if the particle mean life \tau_s is in the range 10^{-14} sec \leq \tau_s \leq 10{-7} sec. We apply this mechanism to minimal gauge mediated supersymmetry-breaking models, in the case in which the lightest stau \tilde\tau_1 is the next-to-lightest supersymmetric particle. We show that, for a large range of phenomenologically acceptable values of the supersymmetry-breaking scale \sqrt{F}, the generated primordial magnetic field can be strong enough to seed the galactic dynamo.

We study the Higgs-boson mass spectrum of the minimal supersymmetric standard model, in which the tree-level CP invariance of the Higgs potential is broken explicitly by loop effects of soft-CP-violating Yukawa interactions related to scalar quarks of the third generation. The analysis is performed by considering the CP-non-invariant renormalization-group improved effective potential through next-to-leading order that includes leading logarithms due to two-loop Yukawa and QCD corrections. We find that the three neutral Higgs particles predicted by the theory may strongly mix with one another, thereby significantly modifying their tree-level couplings to fermions and to the $W^\pm$ and Z bosons. We analyze the phenomenological consequences of such a minimal supersymmetric scenario of explicit CP violation on the production rates of the lightest Higgs particle, and discuss strategies for its potential discovery at high-energy colliders.