Gravitational waves can provide an accurate measurement of the luminosity distance to the source, but cannot provide the source redshift unless the degeneracy between mass and redshift can be broken. This makes it essential to infer the redshift of the source independently to measure the expansion history of the Universe. We show that by exploiting the clustering scale of the gravitational wave sources with galaxies of known redshift, we can infer the expansion history from redshift unknown gravitational wave sources. Using gravitational wave sources with unknown redshift that are detectable from the network of gravitational wave detectors with Advanced LIGO design sensitivity, we will be able to obtain accurate and precise measurements of the local Hubble constant, the expansion history of the universe, and the gravitational wave bias parameter, which captures the distribution of gravitational wave sources with respect to the redshift tracer distribution. This technique is not only limited to the low redshift gravitational wave sources, but will be also applicable to the high redshift gravitational wave sources detectable from Laser Interferometer Space Antenna (LISA), Cosmic Explorer (CE), and Einstein Telescope (ET). Moreover, this method is also going to be applicable to samples of supernovae and fast radio bursts with unknown or photometric redshifts.

Measurements of time delays between multiple quasar images produced by strong lensing are reaching a sensitivity that makes them a promising cosmological probe. Future surveys will provide significantly more measurements, reaching unprecedented depth in redshift, making strong lensing time delay (SLTD) observations competitive with other background probes. We forecast constraints on the nature of dark energy from upcoming SLTD surveys, simulating future catalogues with different numbers of lenses distributed up to redshift $z\sim 1$ and focusing on cosmological parameters such as the Hubble constant $H_0$ and parametrisations of the dark energy equation of state. We also explore the impact of our ability to precisely model the lens mass profile and its environment, on the forecasted constraints. We find that in the most optimistic cases, SLTD will constrain $H_0$ at the level of $\sim 0.1\%$, while the CPL equation of state parameters, $w_0$ and $w_a$, can be determined with errors $\sigma_{w_0}\sim 0.05$ and $\sigma_{w_a}\sim 0.3$, respectively. Furthermore, we investigate the bias introduced when a wrong cosmological model is assumed for the analysis. We find that the value of $H_0$ could be biased up to $10 \sigma$, assuming a perfect knowledge of the lens profile, when a $\Lambda$CDM model is used to analyse data that really belong to a $w$CDM cosmology with $w=-0.9$. Based on these findings, we identify a consistency check of the assumed cosmological model in future SLTD surveys, by splitting the dataset in several redshift bins. Depending on the characteristics of the survey, this could provide a smoking gun for dark energy.

We reconstruct the Equation of State of Dark Energy (EoS) from current data using a non-parametric approach where, rather than assuming a specific time evolution of this function, we bin it in time. We treat the transition between the bins with two different methods, i.e. a smoothed step function and a Gaussian Process reconstruction, investigating whether or not the two approaches lead to compatible results. Additionally, we include in the reconstruction procedure a correlation between the values of the EoS at different times in the form of a theoretical prior that takes into account a set of viability and stability requirements that one can impose on models alternative to $\Lambda$CDM. In such case, we necessarily specialize to broad, but specific classes of alternative models, i.e. Quintessence and Horndeski gravity. We use data coming from CMB, Supernovae and BAO surveys. We find an overall agreement between the different reconstruction methods used; with both approaches, we find a time dependence of the mean of the reconstruction, with different trends depending on the class of model studied. The constant EoS predicted by the $\Lambda$CDM model falls anyway within the $1\sigma$ bounds of our analysis.

The expansion history of the Universe reconstructed from a combination of recent data indicates a preference for a changing Dark Energy (DE) density. Moreover, the DE density appears to be increasing with cosmic time, with its equation of state being below -1 on average, and possibly crossing the so-called phantom divide. Scalar-tensor theories, in which the scalar field mediates a force between matter particles, offer a natural framework in which the effective DE equation of state can be less than -1 and cross the phantom barrier. We consider the generalized Brans-Dicke (GBD) class of scalar-tensor theories and reconstruct their Lagrangian given the effective DE density extracted from recent data. Then, given the reconstructed Lagrangian, we solve for the linear perturbations and investigate the characteristic signatures of these reconstructed GBD in the cosmological observables, such as the cosmic microwave background (CMB) anisotropy, the galaxy number counts, and their cross-correlations. In particular, we demonstrate that the Integrated Sachs-Wolfe (ISW) effect probed by the cross-correlation of CMB with the matter distribution can rule out scalar-tensor theories as the explanation of the observed DE dynamics independently from the laboratory and solar system fifth force constraints.

The tachyonic instability is associated with the unboundedness of the Hamiltonian from below and results in an unstable low-$k$ regime. In the cosmological exploration of modified gravity, it is seldom taken into account, with more focus given to the popular no-ghost and no-gradient conditions. The latter though are intrinsically high-$k$ statements. Here we combine all three conditions into a full set of requirements that we show to guarantee stability on the whole range of cosmological scales. We then explore the impact of the different conditions on the parameter space of scalar-tensor gravity, with particular emphasis on the no-tachyon one. We focus on Horndeski gravity and also consider separately the two subclasses of $f(R)$ and Generalized Brans Dicke theories. We identify several interesting features, for instance in the parameter space of designer $f(R)$ on a $w$CDM background, shedding light on previous findings. When looking at the phenomenological functions $\Sigma$ and $\mu$, associated to the weak lensing and clustering potential respectively, we find that in the case of Generalized Brans Dicke the no-tachyon condition clearly cuts models with $\mu\,,\,\Sigma>1$. This effect is less prevalent in the Horndeski case due to the larger amount of free functions in the theory.

We present a major upgrade of MGCAMB, a patch for the Einstein-Boltzmann solver CAMB used for phenomenological tests of general relativity against cosmological datasets. This new version is compatible with the latest CosmoMC code and includes a consistent implementation of massive neutrinos and dynamical dark energy. The code has been restructured to make it easier to modify to fit the custom needs of specific users. To illustrate the capabilities of the code, we present joint constraints on the modified growth, massive neutrinos and the dark energy equation of state from the latest cosmological observations, including the recent galaxy counts and weak lensing measurements from the Dark Energy Survey, and find a good consistency with the {\Lambda}CDM model.

Ongoing and upcoming cosmological surveys will significantly improve our ability to probe the equation of state of dark energy, $w_{\rm DE}$, and the phenomenology of Large Scale Structure. They will allow us to constrain deviations from the $\Lambda$CDM predictions for the relations between the matter density contrast and the weak lensing and the Newtonian potential, described by the functions $\Sigma$ and $\mu$, respectively. In this work, we derive the theoretical prior for the joint covariance of $w_{\rm DE}$, $\Sigma$ and $\mu$, expected in general scalar-tensor theories with second order equations of motion (Horndeski gravity), focusing on their time-dependence at certain representative scales. We employ Monte-Carlo methods to generate large ensembles of statistically independent Horndeski models, focusing on those that are physically viable and in broad agreement with local tests of gravity, the observed cosmic expansion history and the measurement of the speed of gravitational waves from a binary neutron star merger. We identify several interesting features and trends in the distribution functions of $w_{\rm DE}$, $\Sigma$ and $\mu$, as well as in their covariances; we confirm the high degree of correlation between $\Sigma$ and $\mu$ in scalar-tensor theories. The derived prior covariance matrices will allow us to reconstruct jointly $w_{\rm DE}$, $\Sigma$ and $\mu$ in a non-parametric way.

The splashback radius $r_\mathrm{sp}$ has been identified in cosmological $N$-body simulations as an important scale associated with gravitational collapse and the phase-space distribution of recently accreted material. We employ a semi-analytical approach to study the spherical collapse of dark matter haloes in symmetron gravity and provide insights into how the phenomenology of splashback is affected. The symmetron is a scalar-tensor theory of gravity which exhibits a screening mechanism whereby higher-density regions are screened from the effects of a fifth force. In this model, we find that, as over-densities grow over cosmic time, the inner region becomes heavily screened. In particular, we identify a sector of the parameter space for which material currently sitting at $r_\mathrm{sp}$ has followed, during the collapse, the formation of this screened region. As a result, we find that for this part of the parameter space the splashback radius is maximally affected by the symmetron force and we predict changes in $r_\mathrm{sp}$ up to around $10\%$ compared to its General Relativity value. Because this margin is within the precision of present splashback experiments, we expect this feature to soon provide constraints for symmetron gravity on previously unexplored scales.

We revisit the cosmology of Covariant Galileon gravity in view of the most recent cosmological data sets, including weak lensing. As a higher derivative theory, Covariant Galileon models do not have a $\Lambda$CDM limit and predict a very different structure formation pattern compared with the standard $\Lambda$CDM scenario. Previous cosmological analyses suggest that this model is marginally disfavoured, yet can not be completely ruled out. In this work we use a more recent and extended combination of data, and we allow for more freedom in the cosmology, by including a massive neutrino sector with three different mass hierarchies. We use the Planck measurements of Cosmic Microwave Background temperature and polarization; Baryonic Acoustic Oscillations measurements by BOSS DR12; local measurements of $H_0$; the joint light-curve analysis supernovae sample; and, for the first time, weak gravitational lensing from the KiDS collaboration. We find, that in order to provide a reasonable fit, a non-zero neutrino mass is indeed necessary, but we do not report any sizable difference among the three neutrino hierarchies. Finally, the comparison of the Bayesian Evidence to the $\Lambda$CDM one shows that in all the cases considered, Covariant Galileon models are statistically ruled out by cosmological data.

We propose a simple, nonlocal modification to general relativity (GR) on large scales, which provides a model of late-time cosmic acceleration in the absence of the cosmological constant and with the same number of free parameters as in standard cosmology. The model is motivated by adding to the gravity sector an extra spin-2 field interacting nonlocally with the physical metric coupled to matter. The form of the nonlocal interaction is inspired by the simplest form of the Deser-Woodard (DW) model, $\alpha R\frac{1}{\Box}R$, with one of the Ricci scalars being replaced by a constant $m^{2}$, and gravity is therefore modified in the infrared by adding a simple term of the form $m^2\frac{1}{\Box}R$ to the Einstein-Hilbert term. We study cosmic expansion histories, and demonstrate that the new model can provide background expansions consistent with observations if $m$ is of the order of the Hubble expansion rate today, in contrast to the simple DW model with no viable cosmology. The model is best fit by $w_0\sim-1.075$ and $w_a\sim0.045$. We also compare the cosmology of the model to that of Maggiore and Mancarella (MM), $m^2R\frac{1}{\Box^2}R$, and demonstrate that the viable cosmic histories follow the standard-model evolution more closely compared to the MM model. We further demonstrate that the proposed model possesses the same number of physical degrees of freedom as in GR. Finally, we discuss the appearance of ghosts in the local formulation of the model, and argue that they are unphysical and harmless to the theory, keeping the physical degrees of freedom healthy.

Phenomenological functions $\Sigma$ and $\mu$, also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$, are commonly used to parameterize modifications of the growth of large-scale structure in alternative theories of gravity. We study the values these functions can take in Horndeski theories, i.e. the class of scalar-tensor theories with second order equations of motion. We restrict our attention to models that are in a broad agreement with tests of gravity and the observed cosmic expansion history. In particular, we require the speed of gravity to be equal to the speed of light today, as required by the recent detection of gravitational waves and electromagnetic emission from a binary neutron star merger. We examine the correlations between the values of $\Sigma$ and $\mu$ analytically within the quasi-static approximation, and numerically, by sampling the space of allowed solutions. We confirm that the conjecture made in [Pogosian:2016pwr], that $(\Sigma-1)(\mu -1) \ge 0$ in viable Horndeski theories, holds very well. Along with that, we check the validity of the quasi-static approximation within different corners of Horndeski theory. Our results show that, even with the tight bound on the present day speed of gravitational waves, there is room within Horndeski theories for non-trivial signatures of modified gravity at the level of linear perturbations.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

We investigate the impact of general conditions of theoretical stability and cosmological viability on dynamical dark energy models. As a powerful example, we study whether minimally coupled, single field Quintessence models that are safe from ghost instabilities, can source the CPL expansion history recently shown to be mildly favored by a combination of CMB (Planck) and Weak Lensing (KiDS) data. We find that in their most conservative form, the theoretical conditions impact the analysis in such a way that smooth single field Quintessence becomes significantly disfavored with respect to the standard LCDM cosmological model. This is due to the fact that these conditions cut a significant portion of the (w0;wa) parameter space for CPL, in particular eliminating the region that would be favored by weak lensing data. Within the scenario of a smooth dynamical dark energy parametrized with CPL, weak lensing data favors a region that would require multiple fields to ensure gravitational stability.

EFTCAMB/EFTCosmoMC are publicly available patches to the CAMB/CosmoMC codes implementing the effective field theory approach to single scalar field dark energy and modified gravity models. With the present numerical notes we provide a guide to the technical details of the code. Moreover we reproduce, as they appear in the code, the complete set of the modified equations and the expressions for all the other relevant quantities used to construct these patches. We submit these notes to the arXiv to grant full and permanent access to this material which provides very useful guidance to the numerical implementation of the EFT framework. We will update this set of notes when relevant modifications to the EFTCAMB/EFTCosmoMC codes will be released. The present version is based on the version of EFTCAMB/EFTCosmoMC Sep17.

Constraining the Dark Energy (DE) equation of state, w, is one of the primary science goals of ongoing and future cosmological surveys. In practice, with imperfect data and incomplete redshift coverage, this requires making assumptions about the evolution of w with redshift z. These assumptions can be manifested in a choice of a specific parametric form, which can potentially bias the outcome, or else one can reconstruct w(z) non-parametrically, by specifying a prior covariance matrix that correlates values of w at different redshifts. In this work, we derive the theoretical prior covariance for the effective DE equation of state predicted by general scalar-tensor theories with second order equations of motion (Horndeski theories). This is achieved by generating a large ensemble of possible scalar-tensor theories using a Monte Carlo methodology, including the application of physical viability conditions. We also separately consider the special sub-case of the minimally coupled scalar field, or quintessence. The prior shows a preference for tracking behaviors in the most general case. Given the covariance matrix, theoretical priors on parameters of any specific parametrization of w(z) can also be readily derived by projection.

Phenomenological functions $\Sigma$ and $\mu$ (also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$) are commonly used to parameterize possible modifications of the Poisson equation relating the matter density contrast to the lensing and the Newtonian potentials, respectively. They will be well constrained by future surveys of large scale structure. But what would the implications of measuring particular values of these functions be for modified gravity theories? We ask this question in the context of general Horndeski class of single field scalar-tensor theories with second order equations of motion. We find several consistency conditions that make it possible to rule out broad classes of theories based on measurements of $\Sigma$ and $\mu$ that are independent of their parametric forms. For instance, a measurement of $\Sigma \ne 1$ would rule out all models with a canonical form of kinetic energy, while finding $\Sigma-1$ and $\mu-1$ to be of opposite sign would strongly disfavour the entire class of Horndeski models. We separately examine the large and the small scale limits, the possibility of scale-dependence, and the consistency with bounds on the speed of gravitational waves. We identify sub-classes of Horndeski theories that can be ruled out based on the measured difference between $\Sigma$ and $\mu$.

We present a generalization of the effective field theory (EFT) formalism for dark energy and modified gravity models to include operators with higher order spatial derivatives. This allows the extension of the EFT framework to a wider class of gravity theories such as Horava gravity. We present the corresponding extended action, both in the EFT and the Arnowitt-Deser-Misner (ADM) formalism, and proceed to work out a convenient mapping between the two, providing a self contained and general procedure to translate a given model of gravity into the EFT language at the basis of the Einstein-Boltzmann solver EFTCAMB. Putting this mapping at work, we illustrate, for several interesting models of dark energy and modified gravity, how to express them in the ADM notation and then map them into the EFT formalism. We also provide for the first time, the full mapping of GLPV models into the EFT framework. We next perform a thorough analysis of the physical stability of the generalized EFT action, in absence of matter components. We work out viability conditions that correspond to the absence of ghosts and modes that propagate with a negative speed of sound in the scalar and tensor sector, as well as the absence of tachyonic modes in the scalar sector. Finally, we extend and generalize the phenomenological basis in terms of $\alpha$-functions introduced to parametrize Horndeski models, to cover all theories with higher order spatial derivatives included in our extended action. We elaborate on the impact of the additional functions on physical quantities, such as the kinetic term and the speeds of propagation for scalar and tensor modes.

We consider Horava gravity within the framework of the effective field theory (EFT) of dark energy and modified gravity. We work out a complete mapping of the theory into the EFT language for an action including all the operators which are relevant for linear perturbations with up to sixth order spatial derivatives. We then employ an updated version of the EFTCAMB/EFTCosmoMC package to study the cosmology of the low-energy limit of Horava gravity and place constraints on its parameters using several cosmological data sets. In particular we use cosmic microwave background (CMB) temperature-temperature and lensing power spectra by Planck 2013, WMAP low-l polarization spectra, WiggleZ galaxy power spectrum, local Hubble measurements, Supernovae data from SNLS, SDSS and HST and the baryon acoustic oscillations measurements from BOSS, SDSS and 6dFGS. We get improved upper bounds, with respect to those from Big Bang Nucleosynthesis, on the deviation of the cosmological gravitational constant from the local Newtonian one. At the level of the background phenomenology, we find a relevant rescaling of the Hubble rate at all epoch, which has a strong impact on the cosmological observables; at the level of perturbations, we discuss in details all the relevant effects on the observables and find that in general the quasi-static approximation is not safe to describe the evolution of perturbations. Overall we find that the effects of the modifications induced by the low-energy Horava gravity action are quite dramatic and current data place tight bounds on the theory parameters.

We investigate the impact of modified theories of gravity on the kinetic Sunyaev-Zeldovich (kSZ) effect of the cosmic microwave background. We focus on a specific class of $f(R)$ models of gravity and compare their predictions for the kSZ power spectrum to that of the $\Lambda$CDM model. We use a publicly available modified version of Halofit to properly include the nonlinear matter power spectrum of $f(R)$ in the modeling of the kSZ signal. We find that the well-known modifications of the growth rate of structure in $f(R)$ can indeed induce sizable changes in the kSZ signal, which are more significant than the changes induced by modifications of the expansion history. We discuss prospects of using the kSZ signal as a complementary probe of modified gravity, giving an overview of assumptions and possible caveats in the modeling.

We show how to fully map a specific model of modified gravity into the Einstein-Boltzmann solver EFTCAMB. This approach consists in few steps and allows to obtain the cosmological phenomenology of a model with minimal effort. We discuss all these steps, from the solution of the dynamical equations for the cosmological background of the model to the use of the mapping relations to cast the model into the effective field theory language and use the latter to solve for perturbations. We choose the Hu-Sawicki f(R) model of gravity as our working example. After solving the background and performing the mapping, we interface the algorithm with EFTCAMB and take advantage of the effective field theory framework to integrate the full dynamics of linear perturbations, returning all quantities needed to accurately compare the model with observations. We discuss some observational signatures of this model, focusing on the linear growth of cosmic structures. In particular we present the behavior of $f\sigma_8$ and $E_G$ that, unlike the $\Lambda$CDM scenario, are generally scale dependent in addition to redshift dependent. Finally, we study the observational implications of the model by comparing its cosmological predictions to the Planck 2015 data, including CMB lensing, the WiggleZ galaxy survey and the CFHTLenS weak lensing survey measurements. We find that while WiggleZ data favor a non-vanishing value of the Hu-Sawicki model parameter, $\log_{10}(-f^0_{R})$, and consequently a large value of $\sigma_8$, CFHTLenS drags the estimate of $\log_{10}(-f^0_{R})$ back to the $\Lambda$CDM limit.

We use measurements from the Planck satellite mission and galaxy redshift surveys over the last decade to test three of the basic assumptions of the standard model of cosmology, $\Lambda$CDM: the spatial curvature of the universe, the nature of dark energy and the laws of gravity on large scales. We obtain improved constraints on several scenarios that violate one or more of these assumptions. We measure $w_0=-0.94\pm0.17$ (18\% measurement) and $1+w_a=1.16\pm0.36$ (31\% measurement) for models with a time-dependent equation of state, which is an improvement over current best constraints \citep{Aubourg2014}. In the context of modified gravity, we consider popular scalar tensor models as well as a parametrization of the growth factor. In the case of one-parameter $f(R)$ gravity models with a $\Lambda$CDM background, we constrain $B_0 < 1.36 \times 10^{-5} $ (1$\sigma$ C.L.), which is an improvement by a factor of 4 on the current best \citep{XU2015}. We provide the very first constraint on the coupling parameters of general scalar-tensor theory and stringent constraint on the only free coupling parameter of Chameleon models. We also derive constraints on extended Chameleon models, improving the constraint on the coupling by a factor of 6 on the current best \citep{Hojjati2011} . We also measure $\gamma = 0.612 \pm 0.072$ (11.7\% measurement) for growth index parametrization. We improve all the current constraints by combining results from various galaxy redshift surveys in a coherent way, which includes a careful treatment of scale-dependence introduced by modified gravity.

We introduce EFTCAMB/EFTCosmoMC as publicly available patches to the commonly used CAMB/CosmoMC codes. We briefly describe the structure of the codes, their applicability and main features. To illustrate the use of these patches, we obtain constraints on parametrized pure EFT and designer $f(R)$ models, both on $\Lambda$CDM and $w$CDM background expansion histories, using data from Planck temperature and lensing potential spectra, WMAP low-$\ell$ polarization spectra (WP), and baryon acoustic oscillation (BAO). Upon inspecting theoretical stability of the models on the given background, we find non-trivial parameter spaces that we translate into viability priors. We use different combinations of data sets to show their individual effects on cosmological and model parameters. Our data analysis results show that, depending on the adopted data sets, in the $w$CDM background case this viability priors could dominate the marginalized posterior distributions. Interestingly, with Planck+WP+BAO+lensing data, in $f(R)$ gravity models, we get very strong constraints on the constant dark energy equation of state, $w_0\in(-1,-0.9997)\,\,(95% {\rm C.L.})$.

We revisit the degeneracy between massive neutrinos and generalized theories of gravity in the framework of effective field theory of cosmic acceleration. In particular we consider f(R) theories and a class of non-minimally coupled models parametrized via a coupling to gravity which is linear in the scale factor. In the former case, we find a slightly lower degeneracy with respect to what found so far in the literature, due to the fact that we implement exact designer f(R) models and evolve the full linear dynamics of perturbations. As a consequence, our bounds are slightly tighter on the f(R) parameter but looser on the summed neutrino mass. We also set a new upper bound on the Compton wavelength parameter ${\rm Log}_{10}B_0<-4.1$ at 95% C.L. with fixed summed neutrino mass ($\sum m_{\nu}=0.06$ eV) in f(R) gravity with the combined data sets from cosmic microwave background temperature and lensing power spectra of Planck collaboration as well as galaxy power spectrum from WiggleZ dark energy survey. We do not observe a sizable degeneracy between massive neutrinos and modified gravity in the linear parametrization of non-minimally gravitational coupling model. The analysis is performed with an updated version of the EFTCAMB/EFTCosmoMC package, which is now publicly available and extends the first version of the code with the consistent inclusion of massive neutrinos, tensor modes, several alternative background histories and designer quintessence models.

In general relativity gravitational waves propagate at the speed of light, however in alternative theories of gravity that might not be the case. We study the effects of a modified speed of gravity, $c_T^2$, on the B-modes of the Cosmic Microwave Background (CMB) anisotropy in polarisation. We find that a departure from the light speed value would leave a characteristic imprint on the BB spectrum part induced by tensors, manifesting as a shift in the angular scale of its peaks. We derive constraints by using the available {\it Planck} and BICEP2 datasets showing how $c_T^2$ can be measured, albeit obtaining weak constraints due to the overall poor accuracy of the current BB power spectrum measurements. The present constraint corresponds to $c_T^2 = 1.30 \pm 0.79$ and $c_T^2< 2.85$ at $95%$ C.L. by assuming a power law primordial tensor power spectrum and $c_T^2<2.33$ at $95%$ C.L. if the running of the spectral index is allowed. We derive forecasts for the next generation CMB satellites, which we find capable of tightly constraining $c_T^2$ at percent level, comparable with bounds from binary pulsar measurements, largely due to the absence of degeneracy with other cosmological parameters.

We implement the effective field theory (EFT) approach to dark energy and modified gravity in the public Einstein-Boltzmann solver CAMB. The resulting code, which we dub EFTCAMB, is a powerful and versatile tool that can be used for several objectives. It can be employed to evolve the full dynamics of linear scalar perturbations in any given single field dark energy or modified gravity model, once the latter is mapped into the EFT formalism. It offers a numerical implementation of EFT as a model-independent framework to test gravity on cosmological scales. EFTCAMB has a built-in check for the fulfillment of general stability conditions such as the absence of ghost and superluminal propagation of perturbations. It handles phantom-divide crossing models and does not contain any quasi-static approximation, but rather evolves the full dynamics of perturbations on linear scales. As we will show, the latter is an important feature in view of the accuracy and scale range of upcoming surveys. We show the reliability and applicability of our code by evolving the dynamics of linear perturbations and extracting predictions for power spectra in several models. In particular we perform a thorough analysis of f(R) theories, comparing our outputs with those of an existing code for LCDM backgrounds, and finding an agreement that can reach 0.1% for models with a Compton wavelength consistent with current cosmological data. We then showcase the flexibility of our code studying two different scenarios. First we produce new results for designer f(R) models with a time-varying dark energy equation of state. Second, we extract predictions for linear observables in some parametrized EFT models with a phantom-divide crossing equation of state for dark energy.

The effective field theory (EFT) of dark energy relies on three functions of time to describe the background dynamics. The viability of these functions is investigated here by means of a thorough dynamical analysis. While the system is underdetermined, and one can always find a set of functions reproducing any expansion history, we are able to determine general compatibility conditions for these functions by requiring a viable background cosmology. In particular, we identify a set of variables that allows us to transform the non-autonomous system of equations into an infinite-dimensional one characterized by a significant recursive structure. We then analyze several autonomous sub-systems, obtained truncating the original one at increasingly higher dimension, that correspond to increasingly general models of dark energy and modified gravity. Furthermore, we exploit the recursive nature of the system to draw some general conclusions on the different cosmologies that can be recovered within the EFT formalism and the corresponding compatibility requirements for the EFT functions. The machinery that we set up serves different purposes. It offers a general scheme for performing dynamical analysis of dark energy and modified gravity models within the model independent framework of EFT; the general results, obtained with this technique, can be projected into specific models, as we show in one example. It also can be used to determine appropriate ansatze for the three EFT background functions when studying the dynamics of cosmological perturbations in the context of large scale structure tests of gravity.

At linear order in cosmological perturbations, departures from the growth in the cosmological standard model can be quantified in terms of two functions of redshift $z$ and Fourier number $k$. Previous studies have performed principal component forecasts for several choices of these two functions, based on expected capabilities of upcoming large structure surveys. It is typically found that there will be many well-constrained degrees of freedom. However, not all and, probably most, of these degrees of freedom were physical if the parametrization had allowed for an arbitrary $k$-dependence. In this paper, we restrict the $k$-dependence to that allowed in local theories of gravity under the quasi-static approximation, i.e. ratios of polynomials in $k$, and identify the best constrained features in the ($z$,$k$)-dependence of the commonly considered functions $\mu$ and $\gamma$ as measured by an LSST-like weak lensing survey. We estimate the uncertainty in the measurements of the eigenmodes of modified growth. We find that imposing the theoretical prior on $k$-dependence reduces the number of degrees of freedom and the covariance between parameters. On the other hand, imaging surveys like LSST are not as sensitive to the $z$-dependence as they are to the $k$-dependence of the modified growth functions. This trade off provides us with, more or less, the same number of well-constrained eigenmodes (with respect to our prior) as found before, but now these modes are physical.

The next generation of large scale surveys will not only measure cosmological parameters within the framework of General Relativity, but will also allow for precision tests of the framework itself. At the order of linear perturbations, departures from the growth in the LCDM model can be quantified in terms of two functions of time and Fourier number k. We argue that in local theories of gravity, in the quasi-static approximation, these functions must be ratios of polynomials in k, with the numerator of one function being equal to the denominator of the other. Moreover, the polynomials are even and of second degree in practically all viable models considered today. This means that, without significant loss of generality, one can use data to constraint only five functions of a single variable, instead of two functions of two variables. Furthermore, since the five functions are expected to be slowly varying, one can fit them to data in a non-parametric way with the aid of an explicit smoothness prior. We discuss practical application of this parametrization to forecasts and fits.

Gamma-Ray Bursts (GRBs) have been proposed as a leading candidate for acceleration of ultra high-energy cosmic rays, which would be accompanied by emission of TeV neutrinos produced in proton-photon interactions during acceleration in the GRB fireball. Two analyses using data from two years of the IceCube detector produced no evidence for this neutrino emission, placing strong constraints on models of neutrino and cosmic-ray production in these sources.

We combine recent measurements of Cosmic Microwave Background Anisotropies, Supernovae luminosity distances and Baryonic Acoustic Oscillations to derive constraints on the dark energy equation of state w in the redshift range 0<z<2, using a principal components approach. We find no significant deviations from the expectations of a cosmological constant. However, combining the datasets we find slight indication for w<-1 at low redshift, thus highlighting how these datasets prefer a non-constant w. Nevertheless the cosmological constant is still in agreement with these observations, while we find that some classes of alternative models may be in tension with the inferred w(z) behaviour.

We examine the accuracy of the quasi-static approximation and a parametric form commonly employed for evolving linear perturbations in f(R) gravity theories and placing cosmological constraints. In particular, we analyze the nature and the importance of the near horizon effects that are often neglected. We find that for viable models, with a small present value of the scalaron field, such corrections are entirely negligible with no impact on observables. We also find that the one-parameter form, commonly used to represent the modified equations for linear perturbations in f(R), leads to theoretical systematic errors that are relatively small and, therefore, is adequate for placing constraint on f(R) models.

The first dedicated search for ultra-high energy (UHE) tau neutrinos of astrophysical origin was performed using the IceCube detector in its 22-string configuration with an instrumented volume of roughly 0.25 km^3. The search also had sensitivity to UHE electron and muon neutrinos. After application of all selection criteria to approximately 200 live-days of data, we expect a background of 0.60 +/- 0.19 (stat.) $^{+0.56}_{-0.58}$ (syst.) events and observe three events, which after inspection emerge as being compatible with background but are kept in the final sample. Therefore, we set an upper limit on neutrinos of all-flavors from UHE astrophysical sources at 90% CL of $E^{2} \Phi(\nu_{x}) < 16.3 * 10^-8 GeV cm^-2 sr^-1 s^-1 over an estimated primary neutrino energy range of 340 TeV to 200 PeV.

We forecast the constraints on modified theories of gravity from the cosmic microwave background (CMB) anisotropies bispectrum that arises from correlations between lensing and the Integrated Sachs-Wolfe effect. In models of modified gravity the evolution of the metric potentials is generally altered and the contribution to the CMB bispectrum signal can differ significantly from the one expected in the standard cosmological model.We adopt a parametrised approach and focus on three different classes of models: Linder's growth index, Chameleon-type models and f(R) theories. We show that the constraints on the parameters of the models will significantly improve with future CMB bispectrum measurements.

In this paper we present the results of searches for periodic neutrino emission from a catalog of binary systems. Such modulation, observed in the photon flux, would be caused by the geometry of these systems. In the analysis, the period is fixed by these photon observations, while the phase and duration of the neutrino emission are treated as free parameters to be fit with the data. If the emission occurs during ~20% or less of the total period, this analysis achieves better sensitivity than a time-integrated analysis. We use the IceCube data taken from May 31, 2007 to April 5, 2008 with its 22-string configuration, and from April 5, 2008 to May 20, 2009 with its 40-string configuration. No evidence for neutrino emission is found, with the strongest excess occurring for Cygnus X-3 at 2.1 sigma significance after accounting for trials. Neutrino flux upper limits for both periodic and time-integrated emission are provided.

We report on a measurement of the cosmic ray energy spectrum with the IceTop air shower array, the surface component of the IceCube Neutrino Observatory at the South Pole. The data used in this analysis were taken between June and October, 2007, with 26 surface stations operational at that time, corresponding to about one third of the final array. The fiducial area used in this analysis was 0.122 km^2. The analysis investigated the energy spectrum from 1 to 100 PeV measured for three different zenith angle ranges between 0{\deg} and 46{\deg}. Because of the isotropy of cosmic rays in this energy range the spectra from all zenith angle intervals have to agree. The cosmic-ray energy spectrum was determined under different assumptions on the primary mass composition. Good agreement of spectra in the three zenith angle ranges was found for the assumption of pure proton and a simple two-component model. For zenith angles {\theta} < 30{\deg}, where the mass dependence is smallest, the knee in the cosmic ray energy spectrum was observed between 3.5 and 4.32 PeV, depending on composition assumption. Spectral indices above the knee range from -3.08 to -3.11 depending on primary mass composition assumption. Moreover, an indication of a flattening of the spectrum above 22 PeV were observed.

The next generation of weak lensing surveys will trace the evolution of matter perturbations and gravitational potentials from the matter dominated epoch until today. Along with constraining the dynamics of dark energy, they will probe the relations between matter overdensities, local curvature, and the Newtonian potential. We work with two functions of time and scale to account for any modifications of these relations in the linear regime from those in the LCDM model. We perform a Principal Component Analysis (PCA) to find the eigenmodes and eigenvalues of these functions for surveys like DES and LSST. This paper builds on and significantly extends the PCA analysis of Zhao et al. (2009) in several ways. In particular, we consider the impact of some of the systematic effects expected in weak lensing surveys. We also present the PCA in terms of other choices of the two functions needed to parameterize modified growth on linear scales, and discuss their merits. We analyze the degeneracy between the modified growth functions and other cosmological parameters, paying special attention to the effective equation of state w(z). Finally, we demonstrate the utility of the PCA as an efficient data compression stage which enables one to easily derive constraints on parameters of specific models without recalculating Fisher matrices from scratch.

A search for an excess of muon-neutrinos from dark matter annihilations in the Sun has been performed with the AMANDA-II neutrino telescope using data collected in 812 days of livetime between 2001 and 2006 and 149 days of livetime collected with the AMANDA-II and the 40-string configuration of IceCube during 2008 and early 2009. No excess over the expected atmospheric neutrino background has been observed. We combine these results with the previously published IceCube limits obtained with data taken during 2007 to obtain a total livetime of 1065 days. We provide an upper limit at 90% confidence level on the annihilation rate of captured neutralinos in the Sun, as well as the corresponding muon flux limit at the Earth, both as functions of the neutralino mass in the range 50 GeV-5000 GeV. We also derive a limit on the neutralino-proton spin-dependent and spin-independent cross section. The limits presented here improve the previous results obtained by the collaboration between a factor of two and five, as well as extending the neutralino masses probed down to 50 GeV. The spin-dependent cross section limits are the most stringent so far for neutralino masses above 200 GeV, and well below direct search results in the mass range from 50 GeV to 5 TeV.

Context. Transient neutrino sources such as Gamma-Ray Bursts (GRBs) and Supernovae (SNe) are hypothesized to emit bursts of high-energy neutrinos on a time-scale of \lesssim 100 s. While GRB neutrinos would be produced in high relativistic jets, core-collapse SNe might host soft-relativistic jets, which become stalled in the outer layers of the progenitor star leading to an efficient production of high-energy neutrinos. Aims. To increase the sensitivity to these neutrinos and identify their sources, a low-threshold optical follow-up program for neutrino multiplets detected with the IceCube observatory has been implemented. Methods. If a neutrino multiplet, i.e. two or more neutrinos from the same direction within 100 s, is found by IceCube a trigger is sent to the Robotic Optical Transient Search Experiment, ROTSE. The 4 ROTSE telescopes immediately start an observation program of the corresponding region of the sky in order to detect an optical counterpart to the neutrino events. Results. No statistically significant excess in the rate of neutrino multiplets has been observed and furthermore no coincidence with an optical counterpart was found. Conclusion. The search allows, for the first time, to set stringent limits on current models predicting a high-energy neutrino flux from soft relativistic hadronic jets in core-collapse SNe. We conclude that a sub-population of SNe with typical Lorentz boost factor and jet energy of 10 and 3\times10^{51} erg, respectively, does not exceed 4.2% at 90% confidence.

We report on a search for extremely-high energy neutrinos with energies greater than $10^6$ GeV using the data taken with the IceCube detector at the South Pole. The data was collected between April 2008 and May 2009 with the half completed IceCube array. The absence of signal candidate events in the sample of 333.5 days of livetime significantly improves model independent limit from previous searches and allows to place a limit on the diffuse flux of cosmic neutrinos with an $E^{-2}$ spectrum in the energy range $2.0 \times 10^{6}$ $-$ $6.3 \times 10^{9}$ GeV to a level of $E^2 \phi \leq 3.6 \times 10^{-8}$ ${\rm GeV cm^{-2} sec^{-1}sr^{-1}}$.

The detection of acoustic signals from ultra-high energy neutrino interactions is a promising method to measure the tiny flux of cosmogenic neutrinos expected on Earth. The energy threshold for this process depends strongly on the absolute noise level in the target material. The South Pole Acoustic Test Setup (SPATS), deployed in the upper part of four boreholes of the IceCube Neutrino Observatory, has monitored the noise in Antarctic ice at the geographic South Pole for more than two years down to 500 m depth. The noise is very stable and Gaussian distributed. Lacking an in-situ calibration up to now, laboratory measurements have been used to estimate the absolute noise level in the 10 to 50 kHz frequency range to be smaller than 20 mPa. Using a threshold trigger, sensors of the South Pole Acoustic Test Setup registered acoustic pulse-like events in the IceCube detector volume and its vicinity. Acoustic signals from refreezing IceCube holes and from anthropogenic sources have been used to localize acoustic events. Monte Carlo simulations of sound propagating from the established sources to the SPATS sensors have allowed to check corresponding model expectations. An upper limit on the neutrino flux at energies $E_\nu > 10^{11}$ GeV is derived from acoustic data taken over eight months.

This paper presents searches for flaring sources of neutrinos using the IceCube neutrino telescope. For the first time, a search is performed over the entire parameter space of energy, direction and time looking for neutrino flares of 20 microseconds to a year duration from astrophysical sources among the atmospheric neutrino and muon backgrounds. Searches which integrate over time are less sensitive to flares because they are affected by a larger background of atmospheric neutrinos and muons that can be reduced by the time constraint. Flaring sources considered here, such as active galactic nuclei, soft gamma-ray repeaters and gamma-ray bursts, are promising candidate neutrino emitters. We used mainly data taken between April 5, 2008 and May 20, 2009 by a partially completed configuration of IceCube with 40 strings. For the presented searches an unbinned maximum likelihood method is used with a time-dependent prior to test several different source hypotheses. An "untriggered" search covers any possible time-dependent emission from sources not correlated to any other observation using other astrophysical messengers such as photons. Moreover, a similar time scan is performed for a predefined catalogue of sources that exhibit intense photon flares. Searches triggered by multi-wavelength information on flares from blazars and soft gamma-ray repeaters are performed using the 40 string data and also the data taken by the previous configuration of 22 strings in operation between May 31, 2007 and April 5, 2008. Flares for which extensive and continuous monitoring is available from Fermi-LAT and SWIFT and flares detected by imaging Cherenkov telescopes with shorter time-scale monitoring are considered. The results from all searches are compatible with a fluctuation of the background.

The IceCube Neutrino Observatory is a 1 km$^{3}$ detector currently taking data at the South Pole. One of the main strategies used to look for astrophysical neutrinos with IceCube is the search for a diffuse flux of high-energy neutrinos from unresolved sources. A hard energy spectrum of neutrinos from isotropically distributed astrophysical sources could manifest itself as a detectable signal that may be differentiated from the atmospheric neutrino background by spectral measurement. This analysis uses data from the IceCube detector collected in its half completed configuration which operated between April 2008 and May 2009 to search for a diffuse flux of astrophysical muon neutrinos. A total of 12,877 upward going candidate neutrino events have been selected for this analysis. No evidence for a diffuse flux of astrophysical muon neutrinos was found in the data set leading to a 90 percent C.L. upper limit on the normalization of an $E^{-2}$ astrophysical $\nu_{\mu}$ flux of $8.9 \times 10^{-9} \ \mathrm{GeV \ cm^{-2} \ s^{-1} \ sr^{-1}}$. The analysis is sensitive in the energy range between $35 \ \mathrm{TeV} - 7 \ \mathrm{PeV}$. The 12,877 candidate neutrino events are consistent with atmospheric muon neutrinos measured from 332 GeV to 84 TeV and no evidence for a prompt component to the atmospheric neutrino spectrum is found.

This paper describes the response of the IceCube neutrino telescope located at the geographic South Pole to outbursts of MeV neutrinos from the core collapse of nearby massive stars. IceCube was completed in December 2010 forming a lattice of 5160 photomultiplier tubes that monitor a volume of ~ 1 cubic kilometer in the deep Antarctic ice for particle induced photons. The telescope was designed to detect neutrinos with energies greater than 100 GeV. Owing to subfreezing ice temperatures, the photomultiplier dark noise rates are particularly low. Hence IceCube can also detect large numbers of MeV neutrinos by observing a collective rise in all photomultiplier rates on top of the dark noise. With 2 ms timing resolution, IceCube can detect subtle features in the temporal development of the supernova neutrino burst. For a supernova at the galactic center, its sensitivity matches that of a background-free megaton-scale supernova search experiment. The sensitivity decreases to 20 standard deviations at the galactic edge (30 kpc) and 6 standard deviations at the Large Magellanic Cloud (50 kpc). IceCube is sending triggers from potential supernovae to the Supernova Early Warning System. The sensitivity to neutrino properties such as the neutrino hierarchy is discussed, as well as the possibility to detect the neutronization burst, a short outbreak of electron neutrinos released by electron capture on protons soon after collapse. Tantalizing signatures, such as the formation of a quark star or a black hole as well as the characteristics of shock waves, are investigated to illustrate IceCube's capability for supernova detection.

Between May 2009 and May 2010, the IceCube neutrino detector at the South Pole recorded 32 billion muons generated in air showers produced by cosmic rays with a median energy of 20 TeV. With a data set of this size, it is possible to probe the southern sky for per-mille anisotropy on all angular scales in the arrival direction distribution of cosmic rays. Applying a power spectrum analysis to the relative intensity map of the cosmic ray flux in the southern hemisphere, we show that the arrival direction distribution is not isotropic, but shows significant structure on several angular scales. In addition to previously reported large-scale structure in the form of a strong dipole and quadrupole, the data show small-scale structure on scales between 15 degrees and 30 degrees. The skymap exhibits several localized regions of significant excess and deficit in cosmic ray intensity. The relative intensity of the smaller-scale structures is about a factor of 5 weaker than that of the dipole and quadrupole structure. The most significant structure, an excess localized at right ascension 122.4 degrees and declination -47.4 degrees, extends over at least 20 degrees in right ascension and has a post-trials significance of 5.3 sigma. The origin of this anisotropy is still unknown.

In this paper we report the first observation in the Southern hemisphere of an energy dependence in the Galactic cosmic ray anisotropy up to a few hundred TeV. This measurement was performed using cosmic ray induced muons recorded by the partially deployed IceCube observatory between May 2009 and May 2010. The data include a total of 33$\times 10^{9}$ muon events with a median angular resolution of $\sim3^{\circ}$ degrees. A sky map of the relative intensity in arrival direction over the Southern celestial sky is presented for cosmic ray median energies of 20 and 400 TeV. The same large-scale anisotropy observed at median energies around 20 TeV is not present at 400 TeV. Instead, the high energy skymap shows a different anisotropy structure including a deficit with a post-trial significance of -6.3$\sigma$. This anisotropy reveals a new feature of the Galactic cosmic ray distribution, which must be incorporated into theories of the origin and propagation of cosmic rays.

I study the profile of the Chameleon field around a radially pulsating mass. Focusing on the case in which the background (static) Chameleon profile exhibits a thin-shell, I add small perturbations to the source in the form of time-dependent radial pulsations. It is found that the Chameleon field inherits a time-dependence, there is a resultant scalar radiation from the region of the source and the metric outside the spherically symmetric mass is not static. This has several interesting and potentially testable consequences.

We present the results for a search of high-energy muon neutrinos with the IceCube detector in coincidence with the Crab nebula flare reported on September 2010 by various experiments. Due to the unusual flaring state of the otherwise steady source we performed a prompt analysis of the 79-string configuration data to search for neutrinos that might be emitted along with the observed gamma-rays. We performed two different and complementary data selections of neutrino events in the time window of 10 days around the flare. One event selection is optimized for discovery of E^-2 neutrino spectrum typical of 1st order Fermi acceleration. A similar event selection has also been applied to the 40-string data to derive the time-integrated limits to the neutrino emission from the Crab. The other event selection was optimized for discovery of neutrino spectra with softer spectral index and TeV energy cut-offs as observed for various galactic sources in gamma-rays. The 90% CL best upper limits on the Crab flux during the 10 day flare are 4.73 x 10^-11 cm-2 s-1 TeV-1 for an E^-2 neutrino spectrum and 2.50 x 10^-10 cm-2 s-1 TeV-1 for a softer neutrino spectra of E-2.7, as indicated by Fermi measurements during the flare. IceCube has also set a time-integrated limit on the neutrino emission of the Crab using 375.5 days of livetime of the 40-string configuration data. This limit is compared to existing models of neutrino production from the Crab and its impact on astrophysical parameters is discussed. The most optimistic predictions of some models are already rejected by the IceCube neutrino telescope with more than 90% CL.

The next generation of weak lensing surveys will trace the growth of large scale perturbations through a sequence of epochs, offering an opportunity to test General Relativity (GR) on cosmological scales. We review in detail the parametrization used in MGCAMB to describe the modified growth expected in alternative theories of gravity and generalized dark energy models. We highlight its advantages and examine several theoretical aspects. In particular, we show that the same set of equations can be consistently used on super-horizon and sub-horizon linear scales. We also emphasize the sensitivity of data to scale-dependent features in the growth pattern, and propose using Principal Component Analysis to converge on a practical set of parameters which is most likely to detect departures from GR. The connection with other parametrizations is also discussed.

IceCube has become the first neutrino telescope with a sensitivity below the TeV neutrino flux predicted from gamma-ray bursts if GRBs are responsible for the observed cosmic-ray flux above $10^{18}$ eV. Two separate analyses using the half-complete IceCube detector, one a dedicated search for neutrinos from $p \gamma$-interactions in the prompt phase of the GRB fireball, and the other a generic search for any neutrino emission from these sources over a wide range of energies and emission times, produced no evidence for neutrino emission, excluding prevailing models at 90% confidence.

SN 2008D, a core collapse supernova at a distance of 27 Mpc, was serendipitously discovered by the Swift satellite through an associated X-ray flash. Core collapse supernovae have been observed in association with long gamma-ray bursts and X-ray flashes and a physical connection is widely assumed. This connection could imply that some core collapse supernovae possess mildly relativistic jets in which high-energy neutrinos are produced through proton-proton collisions. The predicted neutrino spectra would be detectable by Cherenkov neutrino detectors like IceCube. A search for a neutrino signal in temporal and spatial correlation with the observed X-ray flash of SN 2008D was conducted using data taken in 2007-2008 with 22 strings of the IceCube detector. Events were selected based on a boosted decision tree classifier trained with simulated signal and experimental background data. The classifier was optimized to the position and a "soft jet" neutrino spectrum assumed for SN 2008D. Using three search windows placed around the X-ray peak, emission time scales from 100 - 10000 s were probed. No events passing the cuts were observed in agreement with the signal expectation of 0.13 events. Upper limits on the muon neutrino flux from core collapse supernovae were derived for different emission time scales and the principal model parameters were constrained.

Self-annihilating or decaying dark matter in the Galactic halo might produce high energy neutrinos detectable with neutrino telescopes. We have conducted a search for such a signal using 276 days of data from the IceCube 22-string configuration detector acquired during 2007 and 2008. The effect of halo model choice in the extracted limit is reduced by performing a search that considers the outer halo region and not the Galactic Center. We constrain any large scale neutrino anisotropy and are able to set a limit on the dark matter self-annihilation cross section of <\sigma_{A}v> \simeq 10^{-22} cm^3/s for WIMP masses above 1 TeV, assuming a monochromatic neutrino line spectrum.

We report on the first search for atmospheric and for diffuse astrophysical neutrino-induced showers (cascades) in the IceCube detector using 257 days of data collected in the year 2007-2008 with 22 strings active. A total of 14 events with energies above 16 TeV remained after event selections in the diffuse analysis, with an expected total background contribution of $8.3\pm 3.6$. At 90% confidence we set an upper limit of $E^2\Phi_{90%CL}<3.6\times10^{-7} GeV \cdot cm^{-2} \cdot s^{-1}\cdot sr^{-1} $ on the diffuse flux of neutrinos of all flavors in the energy range between 24 TeV and 6.6 PeV assuming that $\Phi \propto E^{-2}$ and that the flavor composition of the $\nu_e : \nu_\mu : \nu_\tau$ flux is $1 : 1 : 1$ at the Earth. The atmospheric neutrino analysis was optimized for lower energies. A total of 12 events were observed with energies above 5 TeV. The observed number of events is consistent with the expected background, within the uncertainties.

A measurement of the atmospheric muon neutrino energy spectrum from 100 GeV to 400 TeV was performed using a data sample of about 18,000 up-going atmospheric muon neutrino events in IceCube. Boosted decision trees were used for event selection to reject mis-reconstructed atmospheric muons and obtain a sample of up-going muon neutrino events. Background contamination in the final event sample is less than one percent. This is the first measurement of atmospheric neutrinos up to 400 TeV, and is fundamental to understanding the impact of this neutrino background on astrophysical neutrino observations with IceCube. The measured spectrum is consistent with predictions for the atmospheric muon neutrino plus muon antineutrino flux.

A search for sidereal modulation in the flux of atmospheric muon neutrinos in IceCube was performed. Such a signal could be an indication of Lorentz-violating physics. Neutrino oscillation models, derivable from extensions to the Standard Model, allow for neutrino oscillations that depend on the neutrino's direction of propagation. No such direction-dependent variation was found. A discrete Fourier transform method was used to constrain the Lorentz and CPT-violating coefficients in one of these models. Due to the unique high energy reach of IceCube, it was possible to improve constraints on certain Lorentz-violating oscillations by three orders of magnitude with respect to limits set by other experiments.

We report on the results of the search for extremely-high energy (EHE) neutrinos with energies above $10^7$ GeV obtained with the partially ($\sim$30%) constructed IceCube in 2007. From the absence of signal events in the sample of 242.1 days of effective livetime, we derive a 90% C.L. model independent differential upper limit based on the number of signal events per energy decade at $E^2 \phi_{\nu_e+\nu_\mu+\nu_\tau}\simeq 1.4 \times 10^{-6}$ GeV cm$^{-2}$ sec$^{-1}$ sr$^{-1}$ for neutrinos in the energy range from $3\times10^7$ to $3\times10^9$ GeV.

We test General Relativity (GR) using current cosmological data: the cosmic microwave background (CMB) from WMAP5 (Komatsu et al. 2009), the integrated Sachs-Wolfe (ISW) effect from the cross-correlation of the CMB with six galaxy catalogs (Giannantonio et al. 2008), a compilation of supernovae Type Ia (SNe) including the latest SDSS SNe (Kessler et al. 2009), and part of the weak lensing (WL) data from CFHTLS (Fu et al. 2008, Kilbinger et al. 2009) that probe linear and mildly non-linear scales. We first test a model where the effective Newton's constant, mu, and the ratio of the two gravitational potentials, eta, transit from the GR value to another constant at late times; in this case, we find that standard GR is fully consistent with the combined data. The strongest constraint comes from the ISW effect which would arise from this gravitational transition; the observed ISW signal imposes a tight constraint on a combination of mu and eta that characterizes the lensing potential. Next, we consider four pixels in time and space for each function mu and eta, and perform a Principal Component Analysis (PCA) finding that seven of the resulting eight eigenmodes are consistent with GR within the errors. Only one eigenmode shows a 2-sigma deviation from the GR prediction, which is likely to be due to a systematic effect. However, the detection of such a deviation demonstrates the power of our time- and scale-dependent PCA methodology when combining observations of structure formation and expansion history to test GR.

IceCube Collaboration Contributions to the 2009 International Cosmic Ray Conference

We study the effects of modified theories of gravity on the cosmic microwave background (CMB) anisotropies power spectrum, and in particular on its large scales, where the integrated Sachs-Wolfe (ISW) effect is important. Starting with a general parametrisation, we then specialise to f(R) theories and theories with Yukawa-type interactions between dark matter particles. In these models, the evolution of the metric potentials is altered, and the contribution to the ISW effect can differ significantly from that in the standard model of cosmology. We proceed to compare these predictions with observational data for the CMB and the ISW, performing a full Monte Carlo Markov chain (MCMC) analysis. In the case of f(R) theories, the result is an upper limit on the lengthscale associated to the extra scalar degree of freedom characterising these theories. With the addition of data from the Hubble diagram of Type Ia supernovae, we obtain an upper limit on the lengthscale of the theory of B_0 < 0.4, or correspondingly \lambda_1 < 1900 Mpc/h at 95% c.l. improving previous CMB constraints. For Yukawa-type models we get a bound on the coupling 0.75 < \beta_1 < 1.25 at the 95% c.l. We also discuss the implications of the assumed priors on the estimation of modified gravity parameters, showing that a marginally less conservative choice improves the f(R) constraints to \lambda_1 < 1400 Mpc/h, corresponding to B_0 < 0.2 at 95% c.l.

Over 5,000 PMTs are being deployed at the South Pole to compose the IceCube neutrino observatory. Many are placed deep in the ice to detect Cherenkov light emitted by the products of high-energy neutrino interactions, and others are frozen into tanks on the surface to detect particles from atmospheric cosmic ray showers. IceCube is using the 10-inch diameter R7081-02 made by Hamamatsu Photonics. This paper describes the laboratory characterization and calibration of these PMTs before deployment. PMTs were illuminated with pulses ranging from single photons to saturation level. Parameterizations are given for the single photoelectron charge spectrum and the saturation behavior. Time resolution, late pulses and afterpulses are characterized. Because the PMTs are relatively large, the cathode sensitivity uniformity was measured. The absolute photon detection efficiency was calibrated using Rayleigh-scattered photons from a nitrogen laser. Measured characteristics are discussed in the context of their relevance to IceCube event reconstruction and simulation efforts.

We study the evolution of linear perturbations in metric f(R) models of gravity and identify a potentially observable characteristic scale-dependent pattern in the behavior of cosmological structures. While at the background level viable f(R) models must closely mimic LCDM, the differences in their prediction for the growth of large scale structures can be sufficiently large to be seen with future weak lensing surveys. While working in the Jordan frame, we perform an analytical study of the growth of structures in the Einstein frame, demonstrating the equivalence of the dynamics in the two frames. We also provide a physical interpretation of the results in terms of the dynamics of an effective dark energy fluid with a non-zero shear. We find that the growth of structure in f(R) is enhanced, but that there are no small scale instabilities associated with the additional attractive "fifth force". We then briefly consider some recently proposed observational tests of modified gravity and their utility for detecting the f(R) pattern of structure growth.

We consider contributions to non-Gaussianity of the Cosmic Microwave Background (CMB) from remnants of post-inflationary phase transitions in the very early universe. Such signatures can optimistically be used to discover evidence of new particle physics through cosmological observations. More conservatively they provide a potential obstacle to extracting information about the non-Gaussian nature of primordial density fluctuations from any detection in the CMB. We demonstrate this explicitly by computing the bispectrum from global textures, which may be generated from a wide class of particle physics models.

Future weak lensing surveys will map the evolution of matter perturbations and gravitational potentials, yielding a new test of general relativity on cosmic scales. They will probe the relations between matter overdensities, local curvature, and the Newtonian potential. These relations can be modified in alternative gravity theories or by the effects of massive neutrinos or exotic dark energy fluids. We introduce two functions of time and scale which account for any such modifications in the linear regime. We use a principal component analysis to find the eigenmodes of these functions that cosmological data will constrain. The number of constrained modes gives a model-independent forecast of how many parameters describing deviations from general relativity could be constrained, along with $w(z)$. The modes' scale and time dependence tell us which theoretical models will be better tested.

We constrain the maximum flux from extragalactic neutrino point sources by using diffuse neutrino flux limits. We show that the maximum flux from extragalactic point sources is E^2(dN/dE) < 1.4 x 10^-9 (L_nu/2x10^43 erg/s)^1/3 GeV cm-^2 s^-1 from individual point sources with average neutrino luminosity per decade, L_nu. It depends only slightly on factors such as the inhomogeneous matter density distribution in the local universe, the luminosity distribution, and the assumed spectral index. The derived constraints are at least one order of magnitude below the current experimental limits from direct searches. Significant constraints are also derived on the number density of neutrino sources and on the total neutrino power density.

Theoretical approaches to explaining the observed acceleration of the universe are reviewed. We briefly discuss the evidence for cosmic acceleration, and the implications for standard General Relativity coupled to conventional sources of energy-momentum. We then address three broad methods of addressing an accelerating universe: the introduction of a cosmological constant, its problems and origins; the possibility of dark energy, and the associated challenges for fundamental physics; and the option that an infrared modification of general relativity may be responsible for the large-scale behavior of the universe.

Viable models of modified gravity designed to produce cosmic acceleration at the current epoch, closely mimic the $\Lambda$CDM model at the level of background cosmology. However, this degeneracy is generically broken at the level of linear perturbations, where the modifications induce a peculiar scale-dependent pattern. A similar pattern is expected in models of coupled dark energy. I present the main results published in Pogosian and Silvestri [1], on the growth of structure in f(R) theories of gravity, and in Zhao et al. [2], on the potential of upcoming and future tomographic surveys to detect departures from the growth of cosmic structure expected within General Relativity with a cosmological constant.

The peculiar velocity of the Local Group, reconstructed from inhomogeneities in the local density field, differs in direction and magnitude from the velocity inferred from the Cosmic Microwave Background dipole. We investigate whether generalized theories of gravity, which predict a modified growth of perturbations, are able to alleviate this discrepancy. We introduce a general formalism for calculating the real-space peculiar-velocity field for modified gravity and theories with interactions in the dark sector. For different classes of theories - scalar tensor and higher-dimensional gravity - we reconstruct the Local Group peculiar velocity using groups of galaxies identified in the 2MASS Redshift Survey. We show that, for realistic parameters, modifications to General Relativity cannot account for the angular discrepancy between the reconstructed Local Group velocity and the dipole in the Cosmic Microwave Background.

In alternative theories of gravity, designed to produce cosmic acceleration at the current epoch, the growth of large scale structure can be modified. We study the potential of upcoming and future tomographic surveys such as DES and LSST, with the aid of CMB and supernovae data, to detect departures from the growth of cosmic structure expected within General Relativity. We employ parametric forms to quantify the potential time- and scale-dependent variation of the effective gravitational constant, and the differences between the two Newtonian potentials. We then apply the Fisher matrix technique to forecast the errors on the modified growth parameters from galaxy clustering, weak lensing, CMB, and their cross-correlations across multiple photometric redshift bins. We find that even with conservative assumptions about the data, DES will produce non-trivial constraints on modified growth, and that LSST will do significantly better.

We consider predictions for structure formation from modifications to general relativity in which the Einstein-Hilbert action is replaced by a general function of the Ricci scalar. We work without fixing a gauge, as well as in explicit popular coordinate choices, appropriate for the modification of existing cosmological code. We present the framework in a comprehensive and practical form that can be directly compared to standard perturbation analyses. By considering the full evolution equations, we resolve perceived instabilities previously suggested, and instead find a suppression of perturbations. This result presents significant challenges for agreement with current cosmological structure formation observations. The findings apply to a broad range of forms of f(R) for which the modification becomes important at low curvatures, disfavoring them in comparison with the LCDM scenario. As such, these results provide a powerful method to rule out a wide class of modified gravity models aimed at providing an alternative explanation to the dark energy problem.

The Antarctic Muon And Neutrino Detector Array (AMANDA) has been taking data since 2000 and its data acquisition system was upgraded in January 2003 to read out the complete digitized waveforms from the buried Photomultipliers (PMTs) using Transient Waveform Recorders (TWR). This system currently runs in parallel with the standard AMANDA data acquisition system. Once AMANDA is incorporated into the 1 km^3 detector IceCube, only the TWR system will be kept. We report results from a first atmospheric neutrino analysis on data collected in 2003 with TWR. Good agreement in event rate and angular distribution verify the performance of the TWR system. A search of the northern hemisphere for localized event clusters shows no statistically significant excess, thus a flux limit is calculated, which is in full agreement with previous results based on the standard AMANDA data acquisition system. We also update the status of a search for diffusely distributed neutrinos with ultra high energy (UHE) using data collected by the TWR system.

One way to account for the acceleration of the universe is to modify general relativity, rather than introducing dark energy. Typically, such modifications introduce new degrees of freedom. It is interesting to consider models with no new degrees of freedom, but with a modified dependence on the conventional energy-momentum tensor; the Palatini formulation of $f(R)$ theories is one example. Such theories offer an interesting testing ground for investigations of cosmological modified gravity. In this paper we study the evolution of structure in these ``modified-source gravity'' theories. In the linear regime, density perturbations exhibit scale dependent runaway growth at late times and, in particular, a mode of a given wavenumber goes nonlinear at a higher redshift than in the standard $\Lambda$CDM model. We discuss the implications of this behavior and why there are reasons to expect that the growth will be cut off in the nonlinear regime. Assuming that this holds in a full nonlinear analysis, we briefly describe how upcoming measurements may probe the differences between the modified theory and the standard $\Lambda$CDM model.

We report new limits on cosmic neutrino fluxes from the test flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment, which completed an 18.4 day flight of a prototype long-duration balloon payload, called ANITA-lite, in early 2004. We search for impulsive events that could be associated with ultra-high energy neutrino interactions in the ice, and derive limits that constrain several models for ultra-high energy neutrino fluxes. We rule out the long-standing Z-burst model as the source for the ultra-high energy cosmic rays.

The Antarctic Impulsive Transient Antenna (ANITA) experiment is being developed to search for ultra-high-energy (UHE) neutrino interactions ($>3\times10^{18}$ eV) in the Antarctic ice cap. A neutrino interaction in the ice will produce a radio pulse by the means of the Askaryan effect. The large radio transparency of ice allows for such a pulse to be recorded by a cluster of balloon-borne antennas. The details of the ANITA instrument, now in a construction phase, and the science we hope to achieve is discussed. In order to prepare for the main mission, we have flown ANITA-lite during the 2003/04 austral season. ANITA-lite consisted of two quad-ridge horn antennas and a prototype RF (radio frequency) triggering and recording system. Here we present the results of an impulsive RF background survey of Antarctica, as well as proof-of-principle gain, tracking, and timing calibrations conducted by observing solar radio emissions and calibration radio-pulses. A preliminary UHE neutrino flux limit based on ANITA-lite data is also presented.

We describe a new experiment to search for neutrinos with energies above 3 x 10^18 eV based on the observation of short duration radio pulses that are emitted from neutrino-initiated cascades. The primary objective of the ANtarctic Impulse Transient Antenna (ANITA) mission is to measure the flux of Greisen-Zatsepin-Kuzmin (GZK) neutrinos and search for neutrinos from Active Galactic Nuclei (AGN). We present first results obtained from the successful launch of a 2-antenna prototype instrument (called ANITA-lite) that circled Antarctica for 18 days during the 03/04 Antarctic campaign and show preliminary results from attenuation length studies of electromagnetic waves at radio frequencies in Antarctic ice. The ANITA detector is funded by NASA, and the first flight is scheduled for December 2006.