We present a Fisher information study of the statistical impact of galaxy bias and selection effects on the estimation of key cosmological parameters from galaxy redshift surveys; in particular, the angular diameter distance, Hubble parameter, and linear growth rate at a given redshift, the cold dark matter density, and the tilt and running of the primordial power spectrum. The line-of-sight-dependent selection contributions we include here are known to exist in real galaxy samples. We determine the maximum wavenumber included in the analysis by requiring that the next-order corrections to the galaxy power spectrum or bispectrum, treated here at next-to-leading and leading order, respectively, produce shifts of $\lesssim 0.25\sigma$ on each of the six cosmological parameters. With the galaxy power spectrum alone, selection effects can deteriorate the constraints severely, especially on the linear growth rate. Adding the galaxy bispectrum helps break parameter degeneracies significantly. We find that a joint power spectrum-bispectrum analysis of a Euclid-like survey can still measure the linear growth rate to 10% precision after complete marginalization over selection bias. We also discuss systematic parameter shifts arising from ignoring selection effects and/or other bias parameters, and emphasize that it is necessary to either control selection effects at the percent level or marginalize over them. We obtain similar results for the Roman Space Telescope and HETDEX.

We measure power spectrum response functions in the presence of baryonic physical processes using separate universe simulations with the IllustrisTNG galaxy formation model. The response functions describe how the small-scale power spectrum reacts to long-wavelength perturbations and they can be efficiently measured with the separate universe technique by absorbing the effects of the long modes into a modified cosmology. Specifically, we focus on the total first-order matter power spectrum response to an isotropic density fluctuation $R_1(k,z)$, which is fully determined by the logarithmic derivative of the nonlinear matter power spectrum ${\rm dln}P_m(k,z)/{\rm dln}k$ and the growth-only response function $G_1(k,z)$. We find that $G_1(k,z)$ is not affected by the baryonic physical processes in the simulations at redshifts $z < 3$ and on all scales probed ($k \lesssim 15h/{\rm Mpc}$, i.e. length scales $\gtrsim 0.4 {\rm Mpc}/h$). In practice, this implies that the power spectrum fully specifies the baryonic dependence of its response function. Assuming an idealized lensing survey setup, we evaluate numerically the baryonic impact on the squeezed-lensing bispectrum and the lensing super-sample power spectrum covariance, which are given in terms of responses. Our results show that these higher-order lensing statistics can display varying levels of sensitivity to baryonic effects compared to the power spectrum, with the squeezed-bispectrum being the least sensitive. We also show that ignoring baryonic effects on lensing covariances slightly overestimates the error budget (and is therefore conservative from the point of view of parameter error bars) and likely has negligible impact on parameter biases in inference analyses.

We study the impact that baryon-CDM relative density perturbations $\delta_{bc}$ have on galaxy formation using cosmological simulations with the IllustrisTNG model. These isocurvature (non-adiabatic) perturbations can be induced primordially, if multiple fields are present during inflation, and are generated before baryon-photon decoupling when baryons did not comove with CDM. The presence of long-wavelength $\delta_{bc}$ perturbations in our simulations is mimicked by modifying the ratios of the cosmic densities of baryons $\Omega_b$ and CDM $\Omega_c$, at fixed total matter density $\Omega_m$. We measure the corresponding galaxy bias parameter $b_{\delta}^{bc}$ as the response of galaxy abundances to $\delta_{bc}$. When selecting by total host halo mass, $b_{\delta}^{bc}$ is negative and it decreases with mass and redshift. Stellar-mass selected simulated galaxies show a weaker or even the opposite trend because of the competing effects of $\delta_{bc}$ on the halo mass function and stellar-to-halo-mass relations. We show that simple modeling of the latter two effects describes $b_{\delta}^{bc}$ for stellar-mass-selected objects well. We find $b_{\delta}^{bc} =0.6$ for $M_* = 10^{11}\ M_{\odot}/h$ and $z=0.5$, which is representative of BOSS DR12 galaxies. For $\delta_{bc}$ modes generated by baryon-photon interactions, we estimate the impact on the DR12 power spectrum to be below $1\%$, and shifts on inferred distance and growth rate parameters should not exceed $0.1\%$.

We study the impact that large-scale perturbations of (i) the matter density and (ii) the primordial gravitational potential with local primordial non-Gaussianity (PNG) have on galaxy formation using the IllustrisTNG model. We focus on the linear galaxy bias $b_1$ and the coefficient $b_\phi$ of the scale-dependent bias induced by PNG, which describe the response of galaxy number counts to these two types of perturbations, respectively. We perform our study using separate universe simulations, in which the effect of the perturbations is mimicked by changes to the cosmological parameters: modified cosmic matter density for $b_1$ and modified amplitude $\mathcal{A}_s$ of the primordial scalar power spectrum for $b_\phi$. We find that the widely used universality relation $b_\phi = 2\delta_c(b_1 - 1)$ is a poor description of the bias of haloes and galaxies selected by stellar mass $M_*$, which is instead described better by $b_\phi(M_*) = 2\delta_c(b_1(M_*) - p)$ with $p \in [0.4, 0.7]$. This is explained by the different impact that matter overdensities and local PNG have on the median stellar-to-halo-mass relation. A simple model of this impact allows us to describe the stellar mass dependence of $b_1$ and $b_\phi$ fairly well. Our results also show a nontrivial relation between $b_1$ and $b_\phi$ for galaxies selected by color and black hole mass accretion rate. Our results provide refined priors on $b_\phi$ for local PNG constraints and forecasts using galaxy clustering. Given that the widely used universality relation underpredicts $b_\phi(M_*)$, existing analyses may underestimate the true constraining power on local PNG.

We derive, using functional methods and the bias expansion, the conditional likelihood for observing a specific tracer field given an underlying matter field. This likelihood is necessary for Bayesian-inference methods. If we neglect all stochastic terms apart from the ones appearing in the auto two-point function of tracers, we recover the result of Schmidt et al., 2018. We then rigorously derive the corrections to this result, such as those coming from a non-Gaussian stochasticity (which include the stochastic corrections to the tracer bispectrum) and higher-derivative terms. We discuss how these corrections can affect current applications of Bayesian inference. We comment on possible extensions to our result, with a particular eye towards primordial non-Gaussianity. This work puts on solid theoretical grounds the EFT-based approach to Bayesian forward modeling.

We present updates on the cosmology inference using the effective field theory (EFT) likelihood presented previously in Schmidt et al., 2018, Elsner et al., 2019 [1,2]. Specifically, we add a cutoff to the initial conditions that serve as starting point for the matter forward model. We show that this cutoff, which was not employed in any previous related work, is important to regularize loop integrals that otherwise involve small-scale, non-perturbative modes. We then present results on the inferred value of the linear power spectrum normalization $\sigma_{8}$ from rest-frame halo catalogs using both second- and third-order bias expansions, imposing uniform priors on all bias parameters. Due to the perfect bias-$\sigma_{8}$ degeneracy at linear order, constraints on $\sigma_{8}$ rely entirely on nonlinear information. The results show the expected convergence behavior when lowering the cutoff in wavenumber, $\Lambda$. When including modes up to $k \leq \Lambda = 0.1\,h\,{\rm Mpc}^{-1}$ in the second-order case, $\sigma_{8}$ is recovered to within $\lesssim 6\,\%$ for a range of halo masses and redshifts. The systematic bias shrinks to $4\,\%$ or less for the third-order bias expansion on the same range of scales. Together with additional evidence we provide, this shows that the residual mismatch in $\sigma_{8}$ can be attributed to higher-order bias contributions. We conclude that the EFT likelihood is able to infer unbiased cosmological constraints, within expected theoretical systematic errors, from physical biased tracers on quasilinear scales

In the EFT of biased tracers the noise field $\varepsilon_g$ is not exactly uncorrelated with the nonlinear matter field $\delta$. Its correlation with $\delta$ is effectively captured by adding stochasticities to each bias coefficient. We show that if these stochastic fields are Gaussian (the impact of their non-Gaussianity being subleading on quasi-linear scales anyway) it is possible to resum exactly their effect on the conditional likelihood ${\cal P}[\delta_g|\delta]$ to observe a galaxy field $\delta_g$ given an underlying $\delta$. This resummation allows to take them into account in EFT-based approaches to Bayesian forward modeling. We stress that the resulting corrections to a purely Gaussian conditional likelihood with white-noise covariance are the most relevant on scales where the EFT is under control: they are more important than any non-Gaussianity of the noise $\varepsilon_g$.

We present anisotropic ``separate universe'' simulations which modify the N-body code Gadget-4 in order to represent a large-scale tidal field through an anisotropic expansion factor. These simulations are used to measure the linear, quasi-linear and nonlinear response of the matter power spectrum to a spatially uniform trace-free tidal field up to wavenumber $k = 10 h \text{Mpc}^{-1}$. Together with the response to a large-scale overdensity measured in previous work, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that the response amplitude does not approach zero on small scales in physical coordinates, but rather a constant value, $R_K\approx 0.5$ for $k \geq 3 h\text{Mpc}^{-1}$ at $z=0$. This shows that even the inner regions of haloes are affected by the large-scale tidal field. We also measure directly the alignment of halo shapes with the tidal field, finding a clear signal which increases with halo mass.

Core-collapse supernovae can condense large masses of dust post-explosion. However, sputtering and grain-grain collisions during the subsequent passage of the dust through the reverse shock can potentially destroy a significant fraction of the newly formed dust before it can reach the interstellar medium. Here we show that in oxygen-rich supernova remnants like Cassiopeia A the penetration and trapping within silicate grains of the same impinging ions of oxygen, silicon and magnesium that are responsible for grain surface sputtering can significantly reduce the net loss of grain material. We model conditions representative of dusty clumps (density contrast $\chi=100$) passing through the reverse shock in the oxygen-rich Cassiopeia A remnant and find that, compared to cases where the effect is neglected, as well as facilitating the formation of grains larger than those that had originally condensed, ion trapping increases the surviving masses of silicate dust by factors of up to two to four, depending on initial grain radii. For higher density contrasts ($\chi\gtrsim180$), we find that the effect of gas accretion on the surface of dust grains surpasses ion trapping, and the survival rate increases to ${\sim}55 \%$ of the initial dust mass for $\chi=256$.

We investigate the potential of the galaxy power spectrum to constrain compensated isocurvature perturbations (CIPs), primordial fluctuations in the baryon density that are compensated by fluctuations in CDM density to ensure an unperturbed total matter density. We show that CIPs contribute to the galaxy overdensity at linear order, and if they are close to scale-invariant, their effects are nearly perfectly degenerate with the local PNG parameter $f_{\rm nl}$ if they correlate with the adiabatic perturbations. This degeneracy can however be broken by analyzing multiple galaxy samples with different bias parameters, or by taking CMB priors on $f_{\rm nl}$ into account. Parametrizing the amplitude of the CIP power spectrum as $P_{\sigma\sigma} = A_{\rm CIP}^2P_{\mathcal{R}\mathcal{R}}$ (where $P_{\mathcal{R}\mathcal{R}}$ is the adiabatic power spectrum) we find, for a number of fiducial galaxy samples in a simplified forecast setup, that constraints on $A_{\rm CIP}$, relative to those on $f_{\rm nl}$, of order $\sigma_{A_{\rm CIP}}/\sigma_{f_{\rm nl}} \approx 1-2$ are achievable for CIPs correlated with adiabatic perturbations, and $\sigma_{A_{\rm CIP}}/\sigma_{f_{\rm nl}} \approx 5$ for the uncorrelated case. These values are independent of survey volume, and suggest that current galaxy data are already able to improve significantly on the tightest existing constraints on CIPs from the CMB. Future galaxy surveys that aim to achieve $\sigma_{f_{\rm nl}} \sim 1$ have the potential to place even stronger bounds on CIPs.

The effective-field-theory (EFT) approach to the clustering of galaxies and other biased tracers allows for an isolation of the cosmological information that is protected by symmetries, in particular the equivalence principle, and thus is robust to the complicated dynamics of dark matter, gas, and stars on small scales. All existing implementations proceed by making predictions for the lowest-order $n$-point functions of biased tracers, as well as their covariance, and comparing with measurements. Recently, we presented an EFT-based expression for the conditional probability of the density field of a biased tracer given the matter density field, which in principle combines information from arbitrarily high order $n$-point functions. Here, we report results based on this likelihood by applying it to halo catalogs in real space, specifically on the inference of the power spectrum normalization $\sigma_8$. We include bias terms up to second order as well as the leading higher-derivative term. For a cutoff value of $\Lambda = 0.1 h\,{\rm Mpc}^{-1}$, we recover the ground-truth value of $\sigma_8$ to within 95% CL for different halo samples and redshifts. We discuss possible sources for the remaining systematic bias in $\sigma_8$ as well as future developments.

We present a new simulation technique in which any chosen mode $\mathbf{k}$ of the density contrast field can be amplified by an amplitude $\Delta$. These amplified-mode simulations allow us to study the response of the halo density field to a long-wavelength mode other than the DC mode. In this sense they are a generalization of the separate-universe simulations to finite-wavelength modes. In particular, we use these simulations to obtain robust measurements of the first higher-derivative bias of dark matter halos $b_{\nabla^2 \delta}$. We find a negative bias at all mass considered, roughly following the $-R_L^2(M)$ relation, the Lagrangian radius of halos squared, as naively expected. We compare our results with those obtained from a fit to the 1-loop halo-matter power spectrum, as well as with the recent results from Abidi and Baldauf (2018), and to the prediction from the peak theory.

We present a general perturbative effective field theory (EFT) description of galaxy shape correlations, which are commonly known as intrinsic alignments. This rigorous approach extends current analytical modelling strategies in that it only relies on the equivalence principle. We present our results in terms of three-dimensional statistics for two- and three-point functions of both galaxy shapes and number counts. In case of the two-point function, we recover the well-known linear alignment result at leading order, but also present the full next-to-leading order expressions. In case of the three-point function we present leading order results for all the auto- and cross-correlations of galaxy shapes and densities. We use a spherical tensor basis to decompose the tensor perturbations in different helicity modes, which allows us to make use of isotropy and parity properties in the correlators. Combined with the results on projection presented in a forthcoming companion paper, our framework is directly applicable to accounting for intrinsic alignment contamination in weak lensing surveys, and to extracting cosmological information from intrinsic alignments.

Observations of the ejecta of core-collapse supernovae have shown that dust grains form in over-dense gas clumps in the expanding ejecta. The clumps are later subject to the passage of the reverse shock and a significant amount of the newly formed dust material can be destroyed due to the high temperatures and high velocities in the post-shock gas. To determine dust survival rates, we have performed a set of hydrodynamic simulations using the grid-based code AstroBEAR in order to model a shock wave interacting with a clump of gas and dust. Afterwards, dust motions and dust destruction rates are computed using our newly developed external, post-processing code Paperboats, which includes gas and plasma drag, grain charging, kinematic and thermal sputtering as well as grain-grain collisions. We have determined dust survival rates for the oxygen-rich supernova remnant Cassiopeia A as a function of initial grain sizes, dust materials and clump gas densities.

The reverse shock in the ejecta of core-collapse supernovae is potentially able to destroy newly formed dust material. In order to determine dust survival rates, we have performed a set of hydrodynamic simulations using the grid-based code AstroBEAR in order to model a shock wave interacting with clumpy supernova ejecta. Dust motions and destruction rates were computed using our newly developed external, post-processing code Paperboats, which includes gas drag, grain charging, sputtering and grain-grain collisions. We have determined dust destruction rates for the oxygen-rich supernova remnant Cassiopeia A as a function of initial grain sizes and clump gas density. We found that up to 30 % of the carbon dust mass is able to survive the passage of the reverse shock if the initial grain size distribution is narrow with radii around ~10 - 50 nm for high gas densities, or with radii around ~0.5 - 1.5 ${\mu}$m for low and medium gas densities. Silicate grains with initial radii around 10 - 30 nm show survival rates of up to 40 % for medium and high density contrasts, while silicate material with micron sized distributions is mostly destroyed. For both materials, the surviving dust mass is rearranged into a new size distribution that can be approximated by two components: a power-law distribution of small grains and a log-normal distribution of grains having the same size range as the initial distribution. Our results show that grain-grain collisions and sputtering are synergistic and that grain-grain collisions can play a crucial role in determining the surviving dust budget in supernova remnants.

Europa's surface exhibits morphological features which, associated with a low crater density, might be interpreted to have formed as a result of recent cryovolcanic activity. In particular, the morphology of smooth deposits covering parts of the surface, and their relationship to the surrounding terrains, suggest that they result from liquid extrusions. Furthermore, recent literature suggests that the emplacement of liquid-related features, such as double ridges, lenticulae and chaos could result from the presence of liquid reservoirs beneath the surface. We model the ascent of liquid water through a fracture or a pipe-like conduit from a subsurface reservoir to Europa\textquoteright s surface and calculate the eruption time-scale and the total volume extruded during the eruption, as a function of the reservoir volume and depth. We also estimate the freezing time of a subsurface reservoir necessary to trigger an eruption. Our model is derived for pure liquid water and for a briny mixture outlined by Kargel (1991): 81 wt% H$_{2}$O + 16 wt% MgSO$_{4}$ + 3 wt% Na$_{2}$SO$_{4}$. Considering compositional data for salt impurities for Europa, we discuss the effect of MgSO$_{4}$ and Na$_{2}$SO$_{4}$ on the cryomagma freezing time-scale and ascent. For plausible reservoir volumes and depths in the range of $\mathrm{10^{6}\:m^{3}\leq V\leq10^{10}\:m^{3}}$ and $\mathrm{1\:km\leq H\leq10\:km}$ respectively, the total extruded cryolava volume ranges from $10^{3}\,\mathrm{m^{3}}$ to $10^{8}\,\mathrm{m^{3}}$ and the duration of the eruptions varies from few minutes to few tens of hours. The freezing time-scale of the cryomagma reservoirs varies with cryomagma composition and the temperature gradient in the ice shell: from a few days to a thousand years for pure water cryomagma, and from a few months to a 10$^{4}$ years for briny cryomagma.

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

The galaxy bispectrum contains a wealth of information about the early universe, gravity, as well as astrophysics such as galaxy bias. In this paper, we study the parity-odd part of the galaxy bispectrum which is hitherto unexplored. In the standard cosmological model, the odd-parity bispectrum is generated by galaxy velocities through redshift-space distortions. While small in the case of General Relativity coupled with smooth dark energy, the signal could be larger in modified gravity scenarios. Thus, apart from being a very useful consistency test of measurements of galaxy clustering, the odd bispectrum offers a novel avenue for searching for new physics.

The fact that the spatial nonlocality of galaxy formation is controlled by some short length scale like the Lagrangian radius is the cornerstone of the bias expansion for large-scale-structure tracers. However, the first sources of ionizing radiation between $z\approx 15$ and $z\approx 6$ are expected to have significant effects on the formation of galaxies we observe at lower redshift, at least on low-mass galaxies. These radiative-transfer effects introduce a new scale in the clustering of galaxies, i.e. the finite distance which ionizing radiation travels until it reaches a given galaxy. This mean free path can be very large, of order $100\,h^{-1}\,{\rm Mpc}$. Consequently, higher-derivative terms in the bias expansion could turn out to be non-negligible even on these scales: treating them perturbatively would lead to a massive loss in predictivity and, for example, could spoil the determination of the BAO feature or constraints on the neutrino mass. Here, we investigate under what assumptions an explicit non-perturbative model of radiative-transfer effects can maintain the robustness of large-scale galaxy clustering as a cosmological probe.

The local galaxy bias formalism relies on the energy constraint equation at the formation time to relate the metric perturbation to the matter density contrast. In the Newtonian approximation, this relationship is linear, which allows us to specify the initial galaxy density as a function of local physical operators. In general relativity however, the relationship is intrinsically nonlinear and a modulation of the short-wavelength mode by the long-wavelength mode might be expected. We describe in detail how to obtain local coordinates where the coupling of the long- to the short-wavelength modes is removed through a change of coordinates (in the absence of the primordial non-Gaussianity). We derive the general-relativistic correction to the galaxy bias expansion at second order. The correction does not come from the modulation of small-scale clustering by the long-wavelength mode; instead, it arises from distortions of the volume element by the long-wavelength mode and it does not lead to new bias parameters.

Current technology is not able to map the topography of rocky exoplanets, simply because the objects are too faint and far away to resolve them. Nevertheless, indirect effect of topography should be soon observable thanks to photometry techniques, and the possibility of detecting specular reflections. In addition, topography may have a strong effect on Earth-like exoplanet climates because oceans and mountains affect the distribution of clouds \citep{Houze2012}. Also topography is critical for evaluating surface habitability \citep{Dohm2015}. We propose here a general statistical theory to describe and generate realistic synthetic topographies of rocky exoplanetary bodies. In the solar system, we have examined the best-known bodies: the Earth, Moon, Mars and Mercury. It turns out that despite their differences, they all can be described by multifractral statistics, although with different parameters. Assuming that this property is universal, we propose here a model to simulate 2D spherical random field that mimics a rocky planetary body in a stellar system. We also propose to apply this model to estimate the statistics of oceans and continents to help to better assess the habitability of distant worlds.

This review presents a comprehensive overview of galaxy bias, that is, the statistical relation between the distribution of galaxies and matter. We focus on large scales where cosmic density fields are quasi-linear. On these scales, the clustering of galaxies can be described by a perturbative bias expansion, and the complicated physics of galaxy formation is absorbed by a finite set of coefficients of the expansion, called bias parameters. The review begins with a detailed derivation of this very important result, which forms the basis of the rigorous perturbative description of galaxy clustering, under the assumptions of General Relativity and Gaussian, adiabatic initial conditions. Key components of the bias expansion are all leading local gravitational observables, which include the matter density but also tidal fields and their time derivatives. We hence expand the definition of local bias to encompass all these contributions. This derivation is followed by a presentation of the peak-background split in its general form, which elucidates the physical meaning of the bias parameters, and a detailed description of the connection between bias parameters and galaxy statistics. We then review the excursion-set formalism and peak theory which provide predictions for the values of the bias parameters. In the remainder of the review, we consider the generalizations of galaxy bias required in the presence of various types of cosmological physics that go beyond pressureless matter with adiabatic, Gaussian initial conditions: primordial non-Gaussianity, massive neutrinos, baryon-CDM isocurvature perturbations, dark energy, and modified gravity. Finally, we discuss how the description of galaxy bias in the galaxies' rest frame is related to clustering statistics measured from the observed angular positions and redshifts in actual galaxy catalogs.

We present the complete expression for the next-to-leading (1-loop) order galaxy power spectrum and the leading-order galaxy bispectrum in redshift space in the general bias expansion, or equivalently the effective field theory of biased tracers. We consistently include all line-of-sight dependent selection effects. These are degenerate with many, but not all, of the redshift-space distortion contributions, and have not been consistently derived before. Moreover, we show that, in the framework of effective field theory, a consistent bias expansion in redshift space must include these selection contributions. Physical arguments about the tracer sample considered and its observational selection have to be used to justify neglecting the selection contributions. In summary, the next-to-leading order galaxy power spectrum and leading-order galaxy bispectrum in the general bias expansion are described by 22 parameters, which reduces to 11 parameters if selection effects can be neglected. All contributions to the power spectrum can be written in terms of 28 independent loop integrals.

Conventional approaches to cosmology inference from galaxy redshift surveys are based on n-point functions, which are under rigorous perturbative control on sufficiently large scales. Here, we present an alternative approach, which employs a likelihood at the level of the galaxy density field. By integrating out small-scale modes based on effective-field theory arguments, we prove that this likelihood is under perturbative control if certain specific conditions are met. We further show that the information captured by this likelihood is equivalent to the combination of the next-to-leading order galaxy power spectrum, leading-order bispectrum, and BAO reconstruction. Combined with MCMC sampling and MAP optimization techniques, our results allow for fully Bayesian cosmology inference from large-scale structure that is under perturbative control. We illustrate this via a first demonstration of unbiased cosmology inference from nonlinear large-scale structure using this likelihood. In particular, we show unbiased estimates of the power spectrum normalization $\sigma_{8}$ from a catalog of simulated dark matter halos, where nonlinear information is crucial in breaking the $b_{1} - \sigma_{8}$ degeneracy.

Accurate inference of cosmology from weak lensing shear requires an accurate shear power spectrum covariance matrix. Here, we investigate this accuracy requirement and quantify the relative importance of the Gaussian (G), super-sample covariance (SSC) and connected non-Gaussian (cNG) contributions to the covariance. Specifically, we forecast cosmological parameter constraints for future wide-field surveys and study how different covariance matrix components affect parameter bounds. Our main result is that the cNG term represents only a small and potentially negligible contribution to statistical parameter errors: the errors obtained using the G+SSC subset are within $\lesssim 5\%$ of those obtained with the full G+SSC+cNG matrix for a Euclid-like survey. This result also holds for the shear two-point correlation function, variations in survey specifications and for different analytical prescriptions of the cNG term. The cNG term is that which is often tackled using numerically expensive ensembles of survey realizations. Our results suggest however that the accuracy of analytical or approximate numerical methods to compute the cNG term is likely to be sufficient for cosmic shear inference from the next generation of surveys.

We present direct measurements of cubic bias parameters of dark matter halos from the halo-matter-matter-matter trispectrum. We measure this statistic efficiently by cross-correlating the halo field measured in N-body simulations with specific third-order nonlocal transformations of the initial density field in the same simulation. Together with the recent Abidi & Baldauf (2018), these are the first measurements of halo bias using the four-point function that have been reported to date. We also obtain constraints on the quadratic bias parameters. For all individual cubic parameters involving the tidal field $\mathcal{K}_{ij}$, we find broad consistency with the prediction of the Lagrangian local-in-matter-density ansatz, with some indications of a positive Lagrangian coefficient $b_{\rm td}^L$ multiplying the time derivative of $\mathcal{K}_{ij}$. For the quadratic tidal bias ($b_{K^2}$), we obtain a significant detection of a negative Lagrangian tidal bias.

Long-short mode coupling during inflation, encoded in the squeezed bispectrum of curvature perturbations, induces a dependence of the local, small-scale power spectrum on long-wavelength perturbations, leading to a scale-dependent halo bias. While this scale dependence is absent in the large-scale limit for single-field inflation models that satisfy the consistency relation, certain models such as resonant non-Gaussianity show a peculiar behavior on intermediate scales. We reconsider the predictions for the halo bias in this model by working in Conformal Fermi Coordinates, which isolate the physical effects of long-wavelength perturbations on short-scale physics. We find that the bias oscillates with scale with an envelope similar to that of equilateral non-Gaussianity. Moreover, the bias shows a peculiar modulation with the halo mass. Unfortunately, we find that upcoming surveys will be unable to detect the signal because of its very small amplitude. We also discuss non-Gaussianity due to interactions between the inflaton and massive fields: our results for the bias agree with those in the literature.

We introduce a cosmological model based on the normal branch of DGP braneworld gravity with a smooth dark energy component on the brane. The expansion history in this model is identical to LambdaCDM, thus evading all geometric constraints on the DGP cross-over scale r_c. This model can serve as a first approximation to more general braneworld models whose cosmological solutions have not been obtained yet. We study the formation of large scale structure in this model in the linear and non-linear regime using N-body simulations for different values of r_c. The simulations use the code presented in (F.S., arXiv:0905.0858) and solve the full non-linear equation for the brane-bending mode in conjunction with the usual gravitational dynamics. The brane-bending mode is attractive rather than repulsive in the DGP normal branch, hence the sign of the modified gravity effects is reversed compared to those presented in arXiv:0905.0858. We compare the simulation results with those of ordinary LambdaCDM simulations run using the same code and initial conditions. We find that the matter power spectrum in this model shows a characteristic enhancement peaking at k ~ 0.7 h/Mpc. We also find that the abundance of massive halos is significantly enhanced. Other results presented here include the density profiles of dark matter halos, and signatures of the brane-bending mode self-interactions (Vainshtein mechanism) in the simulations. Independently of the expansion history, these results can be used to place constraints on the DGP model and future generalizations through their effects on the growth of cosmological structure.

The Martian climate is governed by an annual cycle, that results in the condensation of CO$_{2}$ ice during winter, up to a meter thick at the pole and thousands of kilometers in extension. Water and dust may be trapped during the condensation and freed during the sublimation. In addition, ice may be translucent or granular depending on the deposition process (snow vs direct condensation), annealing efficiency, and dust sinking process. The determination of ice translucency is of particular interest to confirm or reject the cold jet model (also known as Kieffer model). This work is focused on the dune field of Richardson Crater in which strong interactions between the water, dust and CO$_{2}$ cycles are observed. We analyzed CRISM hyperspectral images in the near IR using radiative transfer model inversion. We demonstrate that among the states of CO$_{2}$ ice, the translucent state is observed most frequently. The monitoring of surface characteristics shows a decrease in the thickness of the ice during the spring consistently with climate models simulations. We estimate a very low dust content of a few ppmv into the CO$_{2}$ ice, consistent with the formation scenario of cold jets. The water impurities is around 0.1\%v, almost stable during the spring, suggesting a water escape from the surface of subliming CO$_{2}$ ice layer. The water ice grain size varies in a range 1 to 50 microns. From these results, we propose the following new mechanism of small water ice grain suspension: as a cold jet occurs, water ice grains of various sizes are lifted from the surface. These jets happen during daytime, when the general upward gas flux from the subliming CO$_{2}$ ice layer is strong enough to carry the smaller grains, while the bigger fall back on the CO$_{2}$ ice layer. The smaller water grains are carried away and integrated to the general atmospheric circulation.

The Hapke model has been widely used to describe the photometrical behavior of planetary surface through the Bi-directional Reflectance Distribution Function (BRDF), but the uncertainties about retrieved parameters has been difficult to handle so far. A recent study proposed to estimate the uncertainties using a Bayesian approach (Schmidt et al., Icarus 2015). In the present article, we first propose an improved numerical implementation to speed up the uncertainties estimation. Then, we conduct two synthetic studies about photometric measurements in order to analyze the influence of observation geometry: First, we introduce the concept of "efficiency" of a set of geometries to sample the photometric behavior. A set of angular sampling elements (noted as geometry) is efficient if the retrieved Hapke parameters are close to the expected ones. We compared different geometries and found that the principal plane with high incidence is the most efficient geometry among the tested ones. In particular, such geometries are better than poorly sampled full BRDF. Second, we test the analysis scheme of a collection of photometric data acquired from various locations in order to answer the question: are these locations photometrically homogeneous or not? For instance, this question arises when combining data from an entire planetary body, where each spatial position is sampled at a single geometry. We tested the ability of the Bayesian method to decipher two situations, in the presence of noise: (i) a photometrically homogeneous surface (all observations with the same photometric behavior), or (ii) an heterogeneous surface with two distinct photometrical properties (half observations with photometric behavior 1, other half with photometric behavior 2). We show that the naive interpretation of the results provided by Bayesian method is not able to solve this problem.

Topography is the expression of both internal and external processes of a planetary body. Thus hypsometry (the study of topography) is a way to decipher the dynamic of a planet. For that purpose, the statistics of height and slopes may be described by different tools, at local and global scale. We propose here to use the multifractal approach to describe fields of topography. This theory both encompass height and slopes and other statistical moment of the field, tacking into account the scale invariance. Contrary to the widely used fractal formalism, multifractal is able to describe the intermittency of the topography field. As we commonly observe the juxtapostion of rough and smooth at given scale, the multifractal framework seems to be appropriate for hypsometric studies. Here we analyze the data at global scale of the Earth, Mars, Mercury and the Moon and find that the statistics are in good agreement with the multifractal theory for scale larger than 10km. Surprisingly, the analysis shows that all bodies have the same fractal behavior for scale smaller than 10km. We hypothesized that dynamic topography of the mantle may be the explanation at large scale, whereas the smaller scales behavior may be related to elastic thickness.

In this paper we carry out anisotropic "separate universe" simulations by including a large-scale tidal field in the N-body code \textsc{gadget}4 using an anisotropic expansion factor $A_{ij}$. We use the code in a pure \textit{particle-mesh} (PM) mode to simulate the evolution of 16 realizations of an initial density field with and without a large-scale tidal field, which are then used to measure the \textit{response function} describing how the tidal field influences structure formation in the linear and non-linear regimes. Together with the previously measured response to a large scale overdensity, this completely describes the nonlinear matter bispectrum in the squeezed limit. We find that, contrary to the density response, the tidal response never significantly exceeds the large-scale perturbation-theory prediction even on nonlinear scales for the redshift range we discuss. We develop a simple halo model that takes into account the effect of the tidal field and compare it with our direct measurement from the anisotropic N-body simulations.

We derive the complete super-sample covariance (SSC) of the matter and weak lensing convergence power spectra using the power spectrum response formalism to accurately describe the coupling of super- to sub-survey modes. The SSC term is completely characterized by the survey window function, the nonlinear matter power spectrum and the full first-order nonlinear power spectrum response function, which describes the response to super-survey density and tidal field perturbations. Generalized separate universe simulations can efficiently measure these responses in the nonlinear regime of structure formation, which is necessary for lensing applications. We derive the lensing SSC formulae for two cases: one under the Limber and flat-sky approximations, and a more general one that goes beyond the Limber approximation in the super-survey mode and is valid for curved sky applications. Quantitatively, we find that for sky fractions $f_{\rm sky} \approx 0.3$ and a single source redshift at $z_S=1$, the use of the flat-sky and Limber approximation underestimates the total SSC contribution by $\approx 10\%$. The contribution from super-survey tidal fields to the lensing SSC, which has not been included in cosmological analyses so far, is shown to represent about $5\%$ of the total lensing covariance on multipoles $\ell_1,\ell_2 \gtrsim 300$. The SSC is the dominant off-diagonal contribution to the total lensing covariance, making it appropriate to include these tidal terms and beyond flat-sky/Limber corrections in cosmic shear analyses.

Active dark flows known as recurring slope lineae have been observed on the warmest slopes of equatorial Mars. The morphology, composition and seasonality of the lineae suggest a role of liquid water in their formation. However, internal and atmospheric sources of water appear to be insufficient to sustain the observed slope activity. Experimental evidence suggests that under the low atmospheric pressure at the surface of Mars, gas can flow upwards through porous Martian soil due to thermal creep under surface regions heated by the Sun, and disturb small particles. Here we present numerical simulations to demonstrate that such a dry process involving the pumping of rarefied gas in the Martian soil due to temperature contrasts can explain the formation of the recurring slope lineae. In our simulations, solar irradiation followed by shadow significantly reduces the angle of repose due to the resulting temporary temperature gradients over shaded terrain, and leads to flow at intermediate slope angles. The simulated flow locations are consistent with observed recurring slope lineae that initiate in rough and bouldered terrains with local shadows over the soil. We suggest that this dry avalanche process can explain the formation of the recurring slope lineae on Mars without requiring liquid water or CO2 frost activity.

Modified gravity models which include an additional propagating degree of freedom are typically studied in the quasi-static limit, where the propagation is neglected, and the wave equation of the field is replaced with a Poisson-type equation. Recently, it has been proposed that, in the context of models with symmetron- or chameleon-type screening, scalar waves from astrophysical or cosmological events could have a significant effect on the screening of the Solar System, and hence invalidate these models. Here, we quantitatively investigate the impact of scalar waves by solving the full field equation linearised in the wave amplitude. In the symmetron case, we find that the quantitative effect of waves is generally negligible, even for the largest amplitudes of waves that are physically expected. In order to spoil the screening in the Solar System, a significant amount of wave energy would have to be focused on the Solar System by arranging the sources in a spherical shell centred on Earth. In the chameleon case, we are able to rule out any significant effects of propagating waves on Solar System tests.

We delineate the conditions under which the consistency relation for the non-Gaussian bias and the universality of the halo mass function hold in the context of microscopic Lagrangian descriptions of halos. The former is valid provided that the collapse barrier depends only on the physical fields (instead of fields normalized by their variance for example) and explicitly includes the effect of {\it all} physical fields such as the tidal shear. The latter holds provided that the response of the halo number density to a long-wavelength density fluctuation is equivalent to the response induced by shifting the spherical collapse threshold. Our results apply to any Lagrangian halo bias prescription. Effective "moving" barriers, which are ubiquitous in the literature, do not generally satisfy the consistency relation. Microscopic barriers including the tidal shear lead to two additional, second-order Lagrangian bias parameters which ensure that the consistency relation is satisfied. We provide analytic expressions for them.

We present a calculation of the matter power spectrum covariance matrix ${\rm Cov}(\bf{k}_1,\bf{k}_2)$ that uses power spectrum responses to accurately describe the coupling between large- and small-scale modes beyond the perturbative regime. These response functions can be measured with (small-volume) N-body simulations, which is why the response contributions to the covariance remain valid and predictive at all orders in perturbation theory. A novel and key step presented here is the use of responses to compute loop contributions with soft loop momenta, which extends the application of the response approach beyond that of previously considered squeezed $n$-point functions. The calculation presented here does not involve any fitting parameters. When including response-type terms up to 1-loop order in perturbation theory, we find that our calculation captures the bulk of the total covariance as estimated from simulations up to values of $k_1, k_2 \sim 1\ h/{\rm Mpc}$. Moreover, the prediction is guaranteed to be accurate whenever the softer mode is sufficiently linear, ${\rm min}\{k_1,k_2\} \lesssim 0.08\ h/{\rm Mpc}$. We identify and discuss straightforward improvements in the context of the response approach, which are expected to further increase the accuracy of the calculation presented here.

Reionization, the only phase transition in the Universe since recombination, is a key event in the cosmic history of baryonic matter. We derive, in the context of the large-scale bias expansion, the imprints of the epoch of reionization in the large-scale distribution of galaxies, and identify two contributions of particular importance. First, the Compton scattering of CMB photons off the free electrons lead to a drag force on the baryon fluid. This drag induces a relative velocity between baryons and CDM which is of the same order of magnitude as the primordially-induced relative velocity, and enters in the evolution of the relative velocity as calculated by Boltzmann codes. This leads to a unique contribution to galaxy bias involving the matter velocity squared. The second important effect is a modulation of the galaxy density by the ionizing radiation field through radiative-transfer effects, which is captured in the bias expansion by so-called higher-derivative terms. We constrain both of these imprints using the power spectrum of the BOSS DR12 galaxy sample. While they do not lead to a shift in the baryon acoustic oscillation scale, including these terms is important for unbiased cosmology constraints from the shape of the galaxy power spectrum.

Since the discovery of the accelerated expansion of the Universe, the constraints on the equation of state $w_\text{DE}$ of dark energy, the stress-energy component responsible for the acceleration, have tightened significantly. These constraints generally assume an equation of state that is slowly varying in time. We argue that there is good theoretical motivation to consider "monodromic" scenarios with periodic modulations of the dark energy potential. We provide a simple parametrization of such models, and show that these leave room for significant, periodic departures of $w_\text{DE}$ from -1. Moreover, simple models with non-standard kinetic term result in interesting large-scale structure phenomenology beyond that of standard slow-roll dark energy. All these scenarios are best constrained in a dedicated search, as current analyses average over relatively wide redshift ranges.

The squeezed-limit bispectrum, which is generated by nonlinear gravitational evolution as well as inflationary physics, measures the correlation of three wavenumbers, in the configuration where one wavenumber is much smaller than the other two. Since the squeezed-limit bispectrum encodes the impact of a large-scale fluctuation on the small-scale power spectrum, it can be understood as how the small-scale power spectrum "responds" to the large-scale fluctuation. Viewed in this way, the squeezed-limit bispectrum can be calculated using the response approach even in the cases which do not submit to perturbative treatment. To illustrate this point, we apply this approach to the cross-correlation between the large-scale quasar density field and small-scale Lyman-$\alpha$ forest flux power spectrum. In particular, using separate universe simulations which implement changes in the large-scale density, velocity gradient, and primordial power spectrum amplitude, we measure how the Lyman-$\alpha$ forest flux power spectrum responds to the local, long-wavelength quasar overdensity, and equivalently their squeezed-limit bispectrum. We perform a Fisher forecast for the ability of future experiments to constrain local non-Gaussianity using the bispectrum of quasars and the Lyman-$\alpha$ forest. Combining with quasar and Lyman-$\alpha$ forest power spectra to constrain the biases, we find that for DESI the expected $1-\sigma$ constraint is ${\rm err}[f_{\rm NL}]\sim60$. Ability for DESI to measure $f_{\rm NL}$ through this channel is limited primarily by the aliasing and instrumental noise of the Lyman-$\alpha$ forest flux power spectrum. The combination of response approach and separate universe simulations provides a novel technique to explore the constraints from the squeezed-limit bispectrum between different observables.

We measure the large-scale bias of dark matter halos in simulations with non-Gaussian initial conditions of the local type, and compare this bias to the response of the mass function to a change in the primordial amplitude of fluctuations. The two are found to be consistent, as expected from physical arguments, for three halo-finder algorithms which use different Spherical Overdensity (SO) and Friends-of-Friends (FoF) methods. On the other hand, we find that the commonly used prediction for universal mass functions, that the scale-dependent bias is proportional to the first-order Gaussian Lagrangian bias, does not yield a good agreement with the measurements. For all halo finders, high-mass halos show a non-Gaussian bias suppressed by 10-15% relative to the universal mass function prediction. For SO halos, this deviation changes sign at low masses, where the non-Gaussian bias becomes larger than the universal prediction.

Peculiar velocities of objects in the nearby universe are correlated due to the gravitational pull of large-scale structure. By measuring these velocities, we have a unique opportunity to test the cosmological model at the lowest redshifts. We perform this test, using current data to constrain the amplitude of the "signal" covariance matrix describing the velocities and their correlations. We consider a new, well-calibrated "Supercal" set of low-redshift SNe Ia as well as a set of distances derived from the fundamental plane relation of 6dFGS galaxies. Analyzing the SN and galaxy data separately, both results are consistent with the peculiar velocity signal of our fiducial $\Lambda$CDM model, ruling out the noise-only model with zero peculiar velocities at greater than $7\sigma$ (SNe) and $8\sigma$ (galaxies). When the two data sets are combined appropriately, the precision of the test increases slightly, resulting in a constraint on the signal amplitude of $A = 1.05_{-0.21}^{+0.25}$, where $A = 1$ corresponds to our fiducial model. Equivalently, we report an 11% measurement of the product of the growth rate and amplitude of mass fluctuations evaluated at $z_\text{eff} = 0.02$, $f \sigma_8 = 0.428_{-0.045}^{+0.048}$, valid for our fiducial $\Lambda$CDM model. We explore the robustness of the results to a number of conceivable variations in the analysis and find that individual variations shift the preferred signal amplitude by less than ${\sim}0.5\sigma$. We briefly discuss our Supercal SN Ia results in comparison with our previous results using the JLA compilation.

We introduce a rigorous definition of general power-spectrum responses as resummed vertices with two hard and $n$ soft momenta in cosmological perturbation theory. These responses measure the impact of long-wavelength perturbations on the local small-scale power spectrum. The kinematic structure of the responses (i.e., their angular dependence) can be decomposed unambiguously through a "bias" expansion of the local power spectrum, with a fixed number of physical response coefficients, which are only a function of the hard wavenumber $k$. Further, the responses up to $n$-th order completely describe the $(n+2)$-point function in the squeezed limit, i.e. with two hard and $n$ soft modes, which one can use to derive the response coefficients. This generalizes previous results, which relate the angle-averaged squeezed limit to isotropic response coefficients. We derive the complete expression of first- and second-order responses at leading order in perturbation theory, and present extrapolations to nonlinear scales based on simulation measurements of the isotropic response coefficients. As an application, we use these results to predict the non-Gaussian part of the angle-averaged matter power spectrum covariance ${\rm Cov}^{\rm NG}_{\ell = 0}(k_1,k_2)$, in the limit where one of the modes, say $k_2$, is much smaller than the other. Without any free parameters, our model results are in very good agreement with simulations for $k_2 \lesssim 0.06\ h/{\rm Mpc}$, and for any $k_1 \gtrsim 2 k_2$. The well-defined kinematic structure of the power spectrum response also permits a quick evaluation of the angular dependence of the covariance matrix. While we focus on the matter density field, the formalism presented here can be generalized to generic tracers such as galaxies.

We present precise measurements of the assembly bias of dark matter halos, i.e. the dependence of halo bias on other properties than the mass, using curved "separate universe" N-body simulations which effectively incorporate an infinite-wavelength matter overdensity into the background density. This method measures the LIMD (local-in-matter-density) bias parameters $b_n$ in the large-scale limit. We focus on the dependence of the first two Eulerian biases $b^E_1$ and $b^E_2$ on four halo properties: the concentration, spin, mass accretion rate, and ellipticity. We quantitatively compare our results with previous works in which assembly bias was measured on fairly small scales. Despite this difference, our findings are in good agreement with previous results. We also look at the joint dependence of bias on two halo properties in addition to the mass. Finally, using the excursion set peaks model, we attempt to shed new insights on how assembly bias arises in this analytical model.

Density perturbations in cosmology, i.e. spherically symmetric adiabatic perturbations of a Friedmann-Lema\^itre-Robertson-Walker (FLRW) spacetime, are locally exactly equivalent to a different FLRW solution, as long as their wavelength is much larger than the sound horizon of all fluid components. This fact is known as the "separate universe" paradigm. However, no such relation is known for anisotropic adiabatic perturbations, which correspond to an FLRW spacetime with large-scale tidal fields. Here, we provide a closed, fully relativistic set of evolutionary equations for the nonlinear evolution of such modes, based on the conformal Fermi (CFC) frame. We show explicitly that the tidal effects are encoded by the Weyl tensor, and are hence entirely different from an anisotropic Bianchi I spacetime, where the anisotropy is sourced by the Ricci tensor. In order to close the system, certain higher derivative terms have to be dropped. We show that this approximation is equivalent to the local tidal approximation of Hui and Bertschinger (1995). We also show that this very simple set of equations matches the exact evolution of the density field at second order, but fails at third and higher order. This provides a useful, easy-to-use framework for computing the fully relativistic growth of structure at second order.

Gravity is a non-linear theory, and hence, barring cancellations, the initial super-horizon perturbations produced by inflation must contain some minimum amount of mode coupling, or primordial non-Gaussianity. In single-field slow-roll models, where this lower bound is saturated, non-Gaussianity is controlled by two observables: the tensor-to-scalar ratio, which is uncertain by more than fifty orders of magnitude; and the scalar spectral index, or tilt, which is relatively well measured. It is well known that to leading and next-to-leading order in derivatives, the contributions proportional to the tilt disappear from any local observable, and suspicion has been raised that this might happen to all orders, allowing for an arbitrarily low amount of primordial non-Gaussianity. Employing Conformal Fermi Coordinates, we show explicitly that this is not the case. Instead, a contribution of order the tilt appears in local observables. In summary, the floor of physical primordial non-Gaussianity in our Universe has a squeezed-limit scaling of $k_\ell^2/k_s^2$, similar to equilateral and orthogonal shapes, and a dimensionless amplitude of order $0.1\times(n_\mathrm{s}-1)$.

Correlations between intrinsic galaxy shapes on large-scales arise due to the effect of the tidal field of the large-scale structure. Anisotropic primordial non-Gaussianity induces a distinct scale-dependent imprint in these tidal alignments on large scales. Motivated by the observational finding that the alignment strength of luminous red galaxies depends on how galaxy shapes are measured, we study the use of two different shape estimators as a multi-tracer probe of intrinsic alignments. We show, by means of a Fisher analysis, that this technique promises a significant improvement on anisotropic non-Gaussianity constraints over a single-tracer method. For future weak lensing surveys, the uncertainty in the anisotropic non-Gaussianity parameter, $A_2$, is forecast to be $\sigma(A_2)\approx 50$, $\sim 40\%$ smaller than currently available constraints from the bispectrum of the Cosmic Microwave Background. This corresponds to an improvement of a factor of $4-5$ over the uncertainty from a single-tracer analysis.

We present high signal-to-noise galaxy-galaxy lensing measurements of the BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS and CS82. We compare this signal with predictions from mock catalogs trained to match observables including the stellar mass function and the projected and two dimensional clustering of CMASS. We show that the clustering of CMASS, together with standard models of the galaxy-halo connection, robustly predicts a lensing signal that is 20-40% larger than observed. Detailed tests show that our results are robust to a variety of systematic effects. Lowering the value of $S_{\rm 8}=\sigma_{\rm 8} \sqrt{\Omega_{\rm m}/0.3}$ compared to Planck2015 reconciles the lensing with clustering. However, given the scale of our measurement ($r<10$ $h^{-1}$ Mpc), other effects may also be at play and need to be taken into consideration. We explore the impact of baryon physics, assembly bias, massive neutrinos, and modifications to general relativity on $\Delta\Sigma$ and show that several of these effects may be non-negligible given the precision of our measurement. Disentangling cosmological effects from the details of the galaxy-halo connection, the effects of baryons, and massive neutrinos, is the next challenge facing joint lensing and clustering analyses. This is especially true in the context of large galaxy samples from Baryon Acoustic Oscillation surveys with precise measurements but complex selection functions.

We perform a validation of estimates of the growth rate of structure, described by the parameter combination $f\sigma_8$, in modified gravity cosmologies. We consider an analysis pipeline based on the redshift-space distortion modelling of the clustering wedges statistic of the galaxy correlation function and apply it to mock catalogues of $\Lambda{\rm CDM}$ and the normal branch of DGP cosmologies. We employ a halo occupation distribution approach to construct our mocks, which we ensure resemble the CMASS sample from BOSS in terms of the total galaxy number density and large scale amplitude of the power spectrum monopole. We show that the clustering wedges model successfully recovers the true growth rate difference between DGP and $\Lambda{\rm CDM}$, even for cases with over 40\% enhancement in $f\sigma_8$ compared to $\Lambda{\rm CDM}$. The unbiased performance of the clustering wedges model allows us to use the growth rate values estimated from the BOSS DR12 data to constrain the cross-over scale $r_c$ of DGP gravity to $\left[r_cH_0\right]^{-1} < 0.97$ ($2\sigma$) or $r_c > 3090\ {\rm Mpc}/h$, cutting into the interesting region of parameter space with $r_c \sim H_0^{-1}$ using constraints from the growth of structure alone.

Pre-recombination acoustic oscillations induce non-adiabatic perturbations between baryons and dark matter, corresponding to a constant relative-density $\delta_{bc}$ and decaying relative-velocity perturbation $\vec{v}_{bc}$. Due to their significant large-scale correlations and prominent baryon acoustic oscillation (BAO) features, these modes are potentially important for the use of the BAO as standard ruler. We present a complete treatment of the effects of the baryon-CDM perturbations on galaxy clustering in the context of a rigorous perturbative bias expansion. The leading effects are proportional to $\delta_{bc}$ and $\theta_{bc} = \partial_i v_{bc}^i$. We estimate the magnitude of these terms through the excursion set approach. The contribution from $v_{bc}^2$, which has attracted significant attention recently, contributes at subleading (1-loop) order. The relative-density contribution $\delta_{bc}$ is expected to be by far the largest contribution. We also point out contributions to the galaxy velocity bias, the largest of which is simply $v_{bc}$, leading to a term $\propto \mu^2\theta_{bc}$ in the redshift-space galaxy power spectrum $P_g^s(k,\mu)$. Complete expressions of the galaxy power spectrum at 1-loop order are given, which contain several new terms.

Understanding the reason for the observed accelerated expansion of the Universe represents one of the fundamental open questions in physics. In cosmology, a classification has emerged among physical models for the acceleration, distinguishing between Dark Energy and Modified Gravity. In this review, we give a brief overview of models in both categories as well as their phenomenology and characteristic observable signatures in cosmology. We also introduce a rigorous distinction between Dark Energy and Modified Gravity based on the strong and weak equivalence principles.

The origin of the martian methane is still poorly understood. A plausible explanation is that methane could have been produced either by hydrothermal alteration of basaltic crust or by serpentinization of ultramafic rocks producing hydrogen and reducing crustal carbon into methane. Once formed, methane storage on Mars is commonly associated with the presence of hidden clathrate reservoirs. Here, we alternatively suggest that chabazite and clinoptilolite, which belong to the family of zeolites, may form a plausible storage reservoir of methane in the martian subsurface. Because of the existence of many volcanic terrains, zeolites are expected to be widespread on Mars and their Global Equivalent Layer may range up to more than $\sim$1 km, according to the most optimistic estimates. If the martian methane present in chabazite and clinoptilolite is directly sourced from an abiotic source in the subsurface, the destabilization of a localized layer of a few millimeters per year may be sufficient to explain the current observations. The sporadic release of methane from these zeolites requires that they also remained isolated from the atmosphere during its evolution. The methane release over the ages could be due to several mechanisms such as impacts, seismic activity or erosion. If the methane outgassing from excavated chabazite and/or clinoptilolite prevails on Mars, then the presence of these zeolites around Gale Crater could explain the variation of methane level observed by Mars Science Laboratory.

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 show that intrinsic (not lensing-induced) correlations between galaxy shapes offer a new probe of primordial non-Gaussianity and inflationary physics which is complementary to galaxy number counts. Specifically, intrinsic alignment correlations are sensitive to an anisotropic squeezed limit bispectrum of the primordial perturbations. Such a feature arises in solid inflation, as well as more broadly in the presence of light higher spin fields during inflation (as pointed out recently by Arkani-Hamed and Maldacena). We present a derivation of the all-sky two-point correlations of intrinsic shapes and number counts in the presence of non-Gaussianity with general angular dependence, and show that a quadrupolar (spin-2) anisotropy leads to the analog in galaxy shapes of the well-known scale-dependent bias induced in number counts by isotropic (spin-0) non-Gaussianity. Moreover, in presence of non-zero anisotropic non-Gaussianity, the quadrupole of galaxy shapes becomes sensitive to far superhorizon modes. These effects come about because long-wavelength modes induce a local anisotropy in the initial power spectrum, with which galaxies will correlate. We forecast that future imaging surveys could provide constraints on the amplitude of anisotropic non-Gaussianity that are comparable to those from the Cosmic Microwave Background (CMB). These are complementary as they probe different physical scales. The constraints, however, depend on the sensitivity of galaxy shapes to the initial conditions which we only roughly estimate from observed tidal alignments.

The halo model is a theoretically and empirically well-motivated framework for predicting the statistics of the nonlinear matter distribution in the Universe. However, current incarnations of the halo model suffer from two major deficiencies: $(i)$ they do not enforce the stress-energy conservation of matter; $(ii)$ they are not guaranteed to recover exact perturbation theory results on large scales. Here, we provide a formulation of the halo model ("EHM") that remedies both drawbacks in a consistent way, while attempting to maintain the predictivity of the approach. In the formulation presented here, mass and momentum conservation are guaranteed on large scales, and results of perturbation theory and the effective field theory can in principle be matched to any desired order on large scales. We find that a key ingredient in the halo model power spectrum is the halo stochasticity covariance, which has been studied to a much lesser extent than other ingredients such as mass function, bias, and profiles of halos. As written here, this approach still does not describe the transition regime between perturbation theory and halo scales realistically, which is left as an open problem. We also show explicitly that, when implemented consistently, halo model predictions do not depend on any properties of low-mass halos that are smaller than the scales of interest.

In a previous paper (Chassefi\`ere et al., Icarus 223, 878-891, 2013), we have shown that most volcanic sulfur released to early Mars atmosphere could have been trapped in the cryosphere under the form of CO2-SO2 clathrates. Huge amounts of sulfur, up to the equivalent of a ~1 bar atmosphere of SO2, would have been stored in the Noachian cryosphere, then massively released to the atmosphere during Hesperian due to rapidly decreasing CO2 pressure. It would have resulted in the formation of the large sulfate deposits observed mainly in Hesperian terrains, whereas no or little sulfates are found at the Noachian. In the present paper, we first clarify some aspects of our previous work. We discuss the possibility of a smaller cooling effect of sulfur particles, or even of a net warming effect. We point out the fact that CO2-SO2 clathrates formed through a progressive enrichment of a preexisting reservoir of CO2 clathrates and discuss processes potentially involved in the slow formation of a SO2-rich upper cryosphere. We show that episodes of sudden destabilization at the Hesperian may generate 1000 ppmv of SO2 in the atmosphere and contribute to maintaining the surface temperature above the water freezing point.

The supersonic motion of gravitating objects through a gaseous medium constitutes a classical problem in theoretical astrophysics. Its application covers a broad range of objects and scales from planets up to galaxies. Especially the dynamical friction, caused by the forming wake behind the object, plays an important role for the dynamics of the system. To calculate the dynamical friction, standard formulae, based on linear theory are often used. It is our goal to check the general validity of these formulae and provide suitable expressions for the dynamical friction acting on the moving object, based on the basic physical parameters of the problem. We perform sequences of high resolution numerical studies of rigid bodies moving supersonically through a homogeneous medium, and calculate the total drag acting on the object, which is the sum of gravitational and hydro drag. We study cases without gravity with purely hydrodynamical drag, as well as gravitating objects. From the final equilibrium state of the simulations, we compute for gravitating objects the dynamical friction by direct numerical integration of the gravitational pull acting on the embedded object. The numerical experiments confirm the known scaling laws for the dependence of the dynamical friction on the basic physical parameters as derived in earlier semi-analytical studies. As a new important result we find that the shock's stand-off distance is revealed as the minimum spatial interaction scale of dynamical friction. Below this radius, the gas settles into a hydrostatic state, which causes no net gravitational pull onto the moving body. Finally, we derive an analytic estimate for the stand-off distance that can be used to calculate the dynamical friction force.

We revisit the effect of peculiar velocities on low-redshift type Ia supernovae. Velocities introduce an additional guaranteed source of correlations between supernova magnitudes that should be considered in all analyses of nearby supernova samples but has largely been neglected in the past. Applying a likelihood analysis to the latest compilation of nearby supernovae, we find no evidence for the presence of these correlations, although, given the significant noise, the data is also consistent with the correlations predicted for the standard $\Lambda$CDM model. We then consider the dipolar component of the velocity correlations -- the frequently studied "bulk velocity" -- and explicitly demonstrate that including the velocity correlations in the data covariance matrix is crucial for drawing correct and unambiguous conclusions about the bulk flow. In particular, current supernova data is consistent with no excess bulk flow on top of what is expected in $\Lambda$CDM and effectively captured by the covariance. We further clarify the nature of the apparent bulk flow that is inferred when the velocity covariance is ignored. We show that a significant fraction of this quantity is expected to be noise bias due to uncertainties in supernova magnitudes and not any physical peculiar motion.

We present a systematic study of galaxy biasing in the presence of primordial non-Gaussianity. For a large class of non-Gaussian initial conditions, we define a general bias expansion and prove that it is closed under renormalization, thereby showing that the basis of operators in the expansion is complete. We then study the effects of primordial non-Gaussianity on the statistics of galaxies. We show that the equivalence principle enforces a relation between the scale-dependent bias in the galaxy power spectrum and that in the dipolar part of the bispectrum. This provides a powerful consistency check to confirm the primordial origin of any observed scale-dependent bias. Finally, we also discuss the imprints of anisotropic non-Gaussianity as motivated by recent studies of higher-spin fields during inflation.

Disformal theories of gravity are scalar-tensor theories where the scalar couples derivatively to matter via the Jordan frame metric. These models have recently attracted interest in the cosmological context since they admit accelerating solutions. We derive the solution for a static isolated mass in generic disformal gravity theories and transform it into the parameterised post-Newtonian form. This allows us to investigate constraints placed on such theories by local tests of gravity. The tightest constraints come from preferred-frame effects due to the motion of the Solar System with respect to the evolving cosmological background field. The constraints we obtain improve upon the previous solar system constraints by two orders of magnitude, and constrain the scale of the disformal coupling for generic models to $\mathcal{M} \gtrsim 100$ eV. These constraints render all disformal effects irrelevant for cosmology.

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

The abundance of massive galaxy clusters is a powerful probe of departures from General Relativity (GR) on cosmic scales. Despite current stringent constraints placed by stellar and galactic tests, on larger scales alternative theories of gravity such as $f(R)$ can still work as effective theories. Here we present constraints on two popular models of $f(R)$, Hu-Sawicki and "designer", derived from a fully self-consistent analysis of current samples of X-ray selected clusters and accounting for all the covariances between cosmological and astrophysical parameters. Using cluster number counts in combination with recent data from the cosmic microwave background (CMB) and the CMB lensing potential generated by large scale structures, as well as with other cosmological constraints on the background expansion history and its mean matter density, we obtain the upper bounds $\log_{10}|f_{R0}| < 4.79$ and $\log_{10}B_0 < 3.75$ at the 95.4 per cent confidence level, for the Hu-Sawicki (with $n=1$) and designer models, respectively. The robustness of our results derives from high quality cluster growth data for the most massive clusters known out to redshifts $z \sim 0.5$, a tight control of systematic uncertainties including an accurate and precise mass calibration from weak gravitational lensing data, and the use of the full shape of the halo mass function over the mass range of our data.

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.

Hapke proposed a convenient and widely used analytical model to describe the spectro-photometry of granular materials. Using a compilation of the published data, Hapke (2012, Icarus, 221, 1079-1083) recently studied the relationship of b and c for natural examples and proposed the hockey stick relation (excluding b>0.5 and c>0.5). For the moment, there is no theoretical explanation for this relationship. One goal of this article is to study a possible bias due to the retrieval method. We expand here an innovative Bayesian inversion method in order to study into detail the uncertainties of retrieved parameters. On Emission Phase Function (EPF) data, we demonstrate that the uncertainties of the retrieved parameters follow the same hockey stick relation, suggesting that this relation is due to the fact that b and c are coupled parameters in the Hapke model instead of a natural phenomena. Nevertheless, the data used in the Hapke (2012) compilation generally are full Bidirectional Reflectance Diffusion Function (BRDF) that are shown not to be subject to this artifact. Moreover, the Bayesian method is a good tool to test if the sampling geometry is sufficient to constrain the parameters (single scattering albedo, surface roughness, b, c, opposition effect). We performed sensitivity tests by mimicking various surface scattering properties and various single image-like/disk resolved image, EPF-like and BRDF-like geometric sampling conditions. The second goal of this article is to estimate the favorable geometric conditions for an accurate estimation of photometric parameters in order to provide new constraints for future observation campaigns and instrumentations.

We show that dust absorption in disk galaxies leads to a color- and orientation-dependent centroid shift which is expected to be observable in multi-band imaging surveys. This centroid shift is an interesting new probe which contains astrophysically and cosmologically relevant information: it can be used to probe the dust content of a large sample of galaxies, and to reduce the shape noise due to inclination of disk galaxies for weak lensing shear. Specifically, we find that data sets comparable to CFHTLenS, the Dark Energy Survey (DES) or the Hyper Suprime-Cam (HSC) survey should provide a dust measurement for several hundred galaxies per square degree. Conversely, given knowledge of the dust optical depth, this technique will significantly lower the shape noise for the brightest galaxies in the sample (signal-to-noise greater than a few hundred), thereby increasing their relative importance for the weak lensing shear measurement.

(abridged version) The separate universe conjecture states that in General Relativity a density perturbation behaves locally (i.e. on scales much smaller than the wavelength of the mode) as a separate universe with different background density and curvature. We prove this conjecture for a spherical compensated tophat density perturbation of arbitrary amplitude and radius in $\Lambda$CDM. We then use Conformal Fermi Coordinates to generalize this result to scalar perturbations of arbitrary configuration and scale. In this case, the separate universe conjecture holds for the isotropic part of the perturbations. The anisotropic part on the other hand is exactly captured by a tidal field in the Newtonian form. We show that the separate universe picture is restricted to scales larger than the sound horizons of all fluid components. We then derive an expression for the locally measured matter bispectrum induced by a long-wavelength mode of arbitrary wavelength. We show that nonlinear gravitational dynamics does not generate observable contributions that scale like local-type non-Gaussianity $f_{\rm NL}^{\rm loc}$, and hence does not contribute to a scale-dependent galaxy bias $\Delta b \propto k^{-2}$ on large scales; rather, the locally measurable long-short mode coupling assumes a form essentially identical to subhorizon perturbation theory results, once the long-mode density perturbation is replaced by the synchronous-comoving gauge density perturbation. Apparent $f_{\rm NL}^{\rm loc}$-type contributions arise through projection effects on photon propagation, which depend on the specific large-scale structure tracer and observable considered, and are in principle distinguishable from the local mode coupling induced by gravity. We conclude that any observation of $f_{\rm NL}^{\rm loc}$ beyond these projection effects signals a departure from standard single-clock inflation.

Fermi Normal Coordinates (FNC) are a useful frame for isolating the locally observable, physical effects of a long-wavelength spacetime perturbation. Their cosmological application, however, is hampered by the fact that they are only valid on scales much smaller than the horizon. We introduce a generalization that we call Conformal Fermi Coordinates (CFC). CFC preserve all the advantages of FNC, but in addition are valid outside the horizon. They allow us to calculate the coupling of long- and short-wavelength modes on all scales larger than the sound horizon of the cosmological fluid, starting from the epoch of inflation until today, by removing the complications of the second order Einstein equations to a large extent, and eliminating all gauge ambiguities. As an application, we present a calculation of the effect of long-wavelength tensor modes on small scale density fluctuations. We recover previous results, but clarify the physical content of the individual contributions in terms of locally measurable effects and "projection" terms.

The analysis of the surface texture from the particle (grain size, shape and internal structure) to its organization (surface roughness) provides information on the geological processes. CRISM multi-angular observations (varied emission angles) allow to characterize the surface scattering behavior which depends on the composition but also the material physical properties (e.g., grain size, shape, internal structure, the surface roughness). After an atmospheric correction by the Multi-angle Approach for Retrieval of the Surface Reflectance from CRISM Observations, the surface reflectances at different geometries are analyzed by inverting the Hapke photometric model depending on the single scattering albedo, the 2-term phase function, the macroscopic roughness and the 2-term opposition effects. Surface photometric maps are created to observe the spatial variations of surface scattering properties as a function of geological units at the CRISM spatial resolution (200m/pixel). An application at the Mars Exploration Rover (MER) landing sites located at Gusev Crater and Meridiani Planum where orbital and in situ observations are available, is presented. Complementary orbital observations (e.g. CRISM spectra, THermal EMission Imaging System, High Resolution Imaging Science Experiment images) are used for interpreting the estimated Hapke photometric parameters in terms of physical properties. The in situ observations are used as ground truth to validate the interpretations. Varied scattering properties are observed inside a CRISM observation (5x10km) suggesting that the surfaces are controlled by local geological processes (e.g. volcanic resurfacing, aeolian and impact processes) rather than regional or global. Consistent results with the in situ observations are observed thus validating the approach and the use of photometry for the characterization of Martian surface physical properties.

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 study the effect of time evolution on galaxy bias. We argue that at any order in perturbations, the galaxy density contrast can be expressed in terms of a finite set of locally measurable operators made of spatial and temporal derivatives of the Newtonian potential. This is checked in an explicit third order calculation. There is a systematic way to derive a basis for these operators. This basis spans a larger space than the expansion in gravitational and velocity potentials usually employed, although new operators only appear at fourth order. The basis is argued to be closed under renormalization. Most of the arguments also apply to the structure of the counter-terms in the effective theory of large-scale structure.

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.

The detection of B-modes in the Cosmic Microwave Background (CMB) polarization by the BICEP2 experiment, if interpreted as evidence for a primordial gravitational wave background, has enormous ramifications for cosmology and physics. It is crucial to test this hypothesis with independent measurements. A gravitational wave background leads to B-modes in galaxy shape correlations (shear) both through lensing and tidal alignment effects. Since the systematics and foregrounds of galaxy shapes and CMB polarization are entirely different, a detection of a cross-correlation between the two observables would provide conclusive proof for the existence of a primordial gravitational wave background. We find that upcoming weak lensing surveys will be able to detect the cross-correlation between B-modes of the CMB and galaxy shapes. However, this detection is not sufficient to confirm or falsify the hypothesis of a primordial origin for CMB B-mode polarization.

We review recent studies that rigorously define several key observables of the large-scale structure of the Universe in a general relativistic context. Specifically, we consider i) redshift perturbation of cosmic clock events; ii) distortion of cosmic rulers, including weak lensing shear and magnification; iii) observed number density of tracers of the large-scale structure. We provide covariant and gauge-invariant expressions of these observables. Our expressions are given for a linearly perturbed flat Friedmann-Robertson-Walker metric including scalar, vector, and tensor metric perturbations. While we restrict ourselves to linear order in perturbation theory, the approach can be straightforwardly generalized to higher order.

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.

The leading locally observable effect of a long-wavelength metric perturbation corresponds to a tidal field. We derive the tidal field induced by scalar, vector, and tensor perturbations, and use second order perturbation theory to calculate the effect on the locally measured small-scale density fluctuations. For sub-horizon scalar perturbations, we recover the standard perturbation theory result ($F_2$ kernel). For tensor modes of wavenumber $k_L$, we find that effects persist for $k_L\tau \gg 1$, i.e. even long after the gravitational wave has entered the horizon and redshifted away, i.e. it is a "fossil" effect. We then use these results, combined with the "ruler perturbations" of arXiv:1204.3625, to predict the observed distortion of the small-scale matter correlation function induced by a long-wavelength tensor mode. We also estimate the observed signal in the B mode of the cosmic shear from a gravitational wave background, including both tidal (intrinsic alignment) and projection (lensing) effects. The non-vanishing tidal effect in the $k_L\tau \gg 1$ limit significantly increases the intrinsic alignment contribution to shear B modes, especially at low redshifts $z \lesssim 2$.

The B-(curl-)mode of the correlation of galaxy ellipticities (shear) can be used to detect a stochastic gravitational wave background, such as that predicted by inflation. In this paper, we derive the tensor mode contributions to shear from both gravitational lensing and intrinsic alignments, using the gauge-invariant, full-sky results of arXiv:1204.3625. We find that the intrinsic alignment contribution, calculated using the linear alignment model, is larger than the lensing contribution by an order of magnitude or more, if the alignment strength for tensor modes is of the same order as for scalar modes. This contribution also extends to higher multipoles. These results make the prospects for probing tensor modes using galaxy surveys less pessimistic than previously thought, though still very challenging.

The Planetary Fourier Spectrometer (PFS) onboard Mars Express (MEx) is the instrument with the highest spectral resolution observing Mars from orbit since January 2004. It permits studying the atmospheric structure, major and minor compounds. The present time version of the calibration is limited by the effects of mechanical vibration, currently not corrected. We proposed here a new approach to correct for the vibrations based on semi-blind deconvolution of the measurements. This new approach shows that a correction can be done efficiently with 85% reduction of the artefacts, in a equivalent manner to the stacking of 10 spectra. Our strategy is not fully automatic due to the dependence on some regularisation parameters. It may be applied on the complete PFS dataset, correcting the large-scale perturbation due to microvibrations for each spectrum independently. This approach is validated on actual PFS data of Short Wavelength Channel (SWC), perturbed by microvibrations. A coherence check can be performed and also validate our approach. Unfortunately, the coherence check can be done only on the first 310 orbits of MEx only, until the laser line has been switch off. More generally, this work may apply to numerically "deshake" Fourier Transform Spectrometer (FTS), widely used in space experiments or in the laboratory.

The squeezed limit of the three-point function of cosmological perturbations is a powerful discriminant of different models of the early Universe. We present a conceptually simple and complete framework to relate any primordial bispectrum in this limit to late time observables, such as the CMB temperature bispectrum and the scale-dependent halo bias. We employ a series of convenient coordinate transformations to capture the leading non-linear effects of cosmological perturbation theory on these observables. This makes crucial use of Fermi Normal Coordinates and their conformal generalization, which we introduce here and discuss in detail. As an example, we apply our formalism to standard slow-roll single-field inflation. We show explicitly that Maldacena's results for the squeezed limits of the scalar bispectrum [proportional to (ns-1) in comoving gauge] and the tensor-scalar-scalar bispectrum lead to no deviations from a Gaussian universe, except for projection effects. In particular, the primordial contributions to the squeezed CMB bispectrum and scale dependent halo bias vanish, and there are no primordial "fossil" correlations between long-wavelength tensor perturbations and small-scale perturbations. The contributions to observed correlations are then only due to projection effects such as gravitational lensing and redshift perturbations.

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

The quantity and quality of cosmic structure observations have greatly accelerated in recent years. Further leaps forward will be facilitated by imminent projects, which will enable us to map the evolution of dark and baryonic matter density fluctuations over cosmic history. The way that these fluctuations vary over space and time is sensitive to the nature of dark matter and dark energy. Dark energy and gravity both affect how rapidly structure grows; the greater the acceleration, the more suppressed the growth of structure, while the greater the gravity, the more enhanced the growth. While distance measurements also constrain dark energy, the comparison of growth and distance data tests whether General Relativity describes the laws of physics accurately on large scales. Modified gravity models are able to reproduce the distance measurements but at the cost of altering the growth of structure (these signatures are described in more detail in the accompanying paper on Novel Probes of Gravity and Dark Energy). Upcoming surveys will exploit these differences to determine whether the acceleration of the Universe is due to dark energy or to modified gravity. To realize this potential, both wide field imaging and spectroscopic redshift surveys play crucial roles. Projects including DES, eBOSS, DESI, PFS, LSST, Euclid, and WFIRST are in line to map more than a 1000 cubic-billion-light-year volume of the Universe. These will map the cosmic structure growth rate to 1% in the redshift range 0<z<2, over the last 3/4 of the age of the Universe.

We present a derivation of two-point correlations of general tracers in the peak-background split (PBS) framework by way of a rigorous definition of the PBS argument. Our expressions only depend on connected matter correlators and "renormalized" bias parameters with clear physical interpretation, and are independent of any coarse-graining scale. This result should be contrasted with the naive expression derived from a local bias expansion of the tracer number density with respect to the matter density perturbation \delta_L coarse-grained on a scale R_L. In the latter case, the predicted tracer correlation function receives contributions of order <\delta_L^n> at each perturbative order n, whereas, in our formalism, these are absorbed in the PBS bias parameters at all orders. Further, this approach naturally predicts both a scale-dependent bias ~ k^2 such as found for peaks of the density field, and the scale-dependent bias induced by primordial non-Gaussianity in the initial conditions. The only assumption made about the tracers is that their abundance at a given position depends solely on the matter distribution within a finite region around that position.

Primordial non-Gaussianity, in particular the coupling of modes with widely different wavelengths, can have a strong impact on the large-scale clustering of tracers through a scale-dependent bias with respect to matter. We demonstrate that the standard derivation of this non-Gaussian scale-dependent bias is in general valid only in the extreme squeezed limit of the primordial bispectrum, i.e. for clustering over very large scales. We further show how the treatment can be generalized to describe the scale-dependent bias on smaller scales, without making any assumptions on the nature of tracers apart from a dependence on the small-scale fluctuations within a finite region. If the leading scale-dependent bias \Delta b \propto k^{\alpha}, then the first subleading term will scale as k^{\alpha+2}. This correction typically becomes relevant as one considers clustering over scales k >~ 0.01 h Mpc^{-1}.

The comparison between dynamical mass and lensing mass provides a targeted test for a wide range of modified gravity models. In our previous paper we showed, through numerical simulations, that the measurement of the line-of-sight velocity dispersion around stacked massive clusters whose lensing masses are known allows for stringent constraints on modified gravity on scales of 2 - 15 Mpc/h. In this work we develop a semi-analytical approach based on the halo model to describe the phase-space distribution and the line-of-sight velocity dispersion for different tracers. The model distinguishes contributions from the halo pairwise velocity and the virial velocity within halos. We also discuss observational complications, in particular the contribution from Hubble flow, and show how our model can incorporate these complications. We then incorporate the effects of modified gravity (specifically, f(R) and braneworld models), and show that the model predictions are in excellent agreement with modified gravity simulations. More broadly, the phase-space distribution provides a sensitive test of our understanding of hierarchical structure formation when confronted with observations via this model.

We derive general covariant expressions for the six independent observable modes of distortion of ideal standard rulers in a perturbed Friedmann-Robertson-Walker spacetime. Our expressions are gauge-invariant and valid on the full sky. These six modes are most naturally classified in terms of their rotational properties on the sphere, yielding two scalars, two vector (spin-1), and two tensor (spin-2) components. One scalar corresponds to the magnification, while the spin-2 components correspond to the shear. The vector components allow for a polar/axial decomposition analogous to the E/B-decomposition for the shear. Scalar modes do not contribute to the axial (B-)vector, opening a new avenue to probing tensor modes. Our results apply, but are not limited to, the distortion of correlation functions (of the CMB, 21cm emission, or galaxies) as well as to weak lensing shear and magnification, all of which can be seen as methods relying on "standard rulers".

In a perturbed Universe, comoving tracers on a two-dimensional surface of constant observed redshift are at different proper time since the Big Bang. For tracers whose age is known independently, one can measure these perturbations of the proper time. Examples of such sources include cosmic events which only happen during a short period of cosmic history, as well as evolving standard candles and standard rulers. In this paper we derive a general gauge-invariant linear expression for this perturbation in terms of space-time perturbations. As an example, we show that the observed temperature perturbations of the cosmic microwave background (CMB) on large scales are exactly given by these proper time perturbations. Together with the six ruler perturbations derived in Schmidt and Jeong (2012), this completes the set of independent observables which can be measured with standard rulers and candles.

In massive gravity, galileon, and braneworld explanations of cosmic acceleration, force modifications are screened by nonlinear derivative self-interactions of the scalar field mediating that force. Interactions between the field of a central body ("A") and an orbiting body ("B") imply that body B does not move as a test body in the field of body A if the orbit is smaller than the Vainshtein radius of body B. We find through numerical solutions of the joint field at the position of B that the A-field Laplacian is nearly perfectly screened by the B self-field, whereas first derivative or net forces are reduced in a manner that scales with the mass ratio of the bodies as (M_B/M_A)^{3/5}. The latter causes mass-dependent reductions in the universal perihelion precession rate, with deviations for the Earth-Moon system at the ~4% level. In spite of universal coupling, which preserves the microscopic equivalence principle, the motion of macroscopic screened bodies depends on their mass providing in principle a means for testing the Vainshtein mechanism.

The present article proposes an approach to analyze the photometric properties of the surface materials from multi-angle observations acquired by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) on-board the Mars Reconnaissance Orbiter. We estimate photometric parameters using Hapke model in a Bayesian inversion framework. This work also represents a validation of the atmospheric correction provided by the Multi-angle Approach for Retrieval of Surface Reflectance from CRISM Observations (MARS-ReCO) proposed in the companion article.The latter algorithm retrieves photometric curves of surface materials in reflectance units after removing the aerosol contribution. This validation is done by comparing the estimated photometric parameters to those obtained from in situ measurements by Panoramic Camera instrument at the Mars Exploration Rover (MER)-Spirit and MER-Opportunity landing sites. Consistent photometric parameters with those from in situ measurements are found, demonstrating that MARS-ReCO gives access to accurate surface reflectance. Moreover the assumption of a non-Lambertian surface as included in MARS-ReCO is shown to be significantly more precise to estimate surface photometric properties from space in comparison to methods based on a Lambertian surface assumption. In the future, the presented method will allow us to map from orbit the surface bidirectional reflectance and the related photometric parameters in order to characterize the Martian surface.

We consider the prospects for measuring the pairwise kinetic Sunyaev-Zel'dovich (kSZ) signal from galaxy clusters discovered in large photometric surveys such as the Dark Energy Survey (DES). We project that the DES cluster sample will, in conjunction with existing mm-wave data from the South Pole Telescope (SPT), yield a detection of the pairwise kSZ signal at the 8-13 sigma level, with sensitivity peaking for clusters separated by ~100 Mpc distances. A next-generation version of SPT would allow for a 18-30 sigma detection and would be limited by variance from the kSZ signal itself and residual thermal Sunyaev-Zel'dovich (tSZ) signal. Throughout our analysis we assume photometric redshift errors, which wash out the signal for clusters separated by <~50 Mpc; a spectroscopic survey of the DES sample would recover this signal and allow for a 26-43 sigma detection, and would again be limited by kSZ/tSZ variance. Assuming a standard model of structure formation, these high-precision measurements of the pairwise kSZ signal will yield detailed information on the gas content of the galaxy clusters. Alternatively, if the gas can be sufficiently characterized by other means (e.g. using tSZ, X-ray, or weak lensing), then the relative velocities of the galaxy clusters can be isolated, thereby providing a precision measurement of gravity on 100 Mpc scales. We briefly consider the utility of these measurements for constraining theories of modified gravity.

Non-Gaussianity in the inflationary perturbations can couple observable scales to modes of much longer wavelength (even superhorizon), leaving as a signature a large-angle modulation of the observed cosmic microwave background (CMB) power spectrum. This provides an alternative origin for a power asymmetry which is otherwise often ascribed to a breaking of statistical isotropy. The non-Gaussian modulation effect can be significant even for typical ~10^{-5} perturbations, while respecting current constraints on non-Gaussianity, if the squeezed limit of the bispectrum is sufficiently infrared divergent. Just such a strongly infrared-divergent bispectrum has been claimed for inflation models with a non-Bunch-Davies initial state, for instance. Upper limits on the observed CMB power asymmetry place stringent constraints on the duration of inflation in such models.

Observed angular positions and redshifts of large-scale structure tracers such as galaxies are affected by gravitational waves through volume distortion and magnification effects. Thus, a gravitational wave background can in principle be probed through clustering statistics of large-scale structure. We calculate the observed angular clustering of galaxies in the presence of a gravitational wave background at linear order including all relativistic effects. For a scale-invariant spectrum of gravitational waves, the effects are most significant at the smallest multipoles (2 <= l <= 5), but typically suppressed by six or more orders of magnitude with respect to scalar contributions for currently allowed amplitudes of the inflationary gravitational wave background. We also discuss the most relevant second-order terms, corresponding to the distortion of tracer correlation functions by gravitational waves. These provide a natural application of the approach recently developed in arXiv:1204.3625.

We use Luminous Red Galaxies from the Sloan Digital Sky Survey II to test the cosmological structure growth in two alternatives to the standard LCDM+GR cosmological model. We compare observed three-dimensional clustering in SDSS DR7 with theoretical predictions for the standard vanilla LCDM+GR model, Unified Dark Matter cosmologies and the normal branch DGP. In computing the expected correlations in UDM cosmologies, we derive a parameterized formula for the growth factor in these models. For our analysis we apply the methodology tested in Raccanelli et al. 2010 and use the measurements of Samushia et al. 2011, that account for survey geometry, non-linear and wide-angle effects and the distribution of pair orientation. We show that the estimate of the growth rate is potentially degenerate with wide-angle effects, meaning that extremely accurate measurements of the growth rate on large scales will need to take such effects into account. We use measurements of the zeroth and second order moments of the correlation function from SDSS DR7 data and the Large Suite of Dark Matter Simulations, and perform a likelihood analysis to constrain the parameters of the models. Using information on the clustering up to r_max = 120 Mpc/h, and after marginalizing over the bias, we find, for UDM models, a speed of sound < 6.1e-4, and, for the nDGP model, a cross-over scale r_c > 340 Mpc, at 95% confidence level.

In General Relativity, the average velocity field of dark matter around galaxy clusters is uniquely determined by the mass profile. The latter can be measured through weak lensing. We propose a new method of measuring the velocity field (phase space density) by stacking redshifts of surrounding galaxies from a spectroscopic sample. In combination with lensing, this yields a direct test of gravity on scales of 1-30 Mpc. Using N-body simulations, we show that this method can improve upon current constraints on f(R) and DGP model parameters by several orders of magnitude when applied to upcoming imaging and redshift surveys.

We show how the peak-background split can be generalized to predict the effect of non-local primordial non-Gaussianity on the clustering of halos. Our approach is applicable to arbitrary primordial bispectra. We show that the scale-dependence of halo clustering predicted in the peak-background split (PBS) agrees with that of the local-biasing model on large scales. On smaller scales, k >~ 0.01 h/Mpc, the predictions diverge, a consequence of the assumption of separation of scales in the peak-background split. Even on large scales, PBS and local biasing do not generally agree on the amplitude of the effect outside of the high-peak limit. The scale dependence of the biasing - the effect that provides strong constraints to the local-model bispectrum - is far weaker for the equilateral and self-ordering-scalar-field models of non-Gaussianity. The bias scale dependence for the orthogonal and folded models is weaker than in the local model (~ 1/k), but likely still strong enough to be constraining. We show that departures from scale-invariance of the primordial power spectrum may lead to order-unity corrections, relative to predictions made assuming scale-invariance - to the non-Gaussian bias in some of these non-local models for non-Gaussianity. An Appendix shows that a non-local model can produce the local-model bispectrum, a mathematical curiosity we uncovered in the course of this investigation.

We present a new test of gravitational interactions at the r\sim(0.2-20)Mpc scale, around the virial radius of dark matter halos measured through cluster-galaxy lensing of maxBCG clusters from the Sloan Digital Sky Survey (SDSS). We employ predictions from self-consistent simulations of f(R) gravity to find an upper bound on the background field amplitude of f_R0<3.5x10^-3 at the 1D-marginalized 95% confidence level. As a model-independent assessment of the constraining power of cluster profiles measured through weak gravitational lensing, we also constrain the amplitude F_0 of a phenomenological modification based on the profile enhancement induced by f(R) gravity when not including effects from the increased cluster abundance in f(R). In both scenarios, dark-matter-only simulations of the concordance model corresponding to f_R0=0 and F_0=0 are consistent with the lensing measurements, i.e., at the 68% and 95% confidence level, respectively.

Several recent studies have shown how to properly calculate the observed clustering of galaxies in a relativistic context, and uncovered corrections to the Newtonian calculation that become significant on scales near the horizon. Here, we retrace these calculations and show that, on scales approaching the horizon, the observed galaxy power spectrum depends strongly on which gauge is assumed to relate the intrinsic fluctuations in galaxy density to matter perturbations through a linear bias relation. Starting from simple physical assumptions, we derive a gauge-invariant expression relating galaxy density perturbations to matter density perturbations on large scales, and show that it reduces to a linear bias relation in synchronous-comoving gauge, corroborating an assumption made in several recent papers. We evaluate the resulting observed galaxy power spectrum, and show that it leads to corrections similar to an effective non-Gaussian bias corresponding to a local (effective) fNL < 0.5. This number can serve as a guideline as to which surveys need to take into account relativistic effects. We also discuss the scale-dependent bias induced by primordial non-Gaussianity in the relativistic context, which again is simplest in synchronous-comoving gauge.

We investigate the constraining power of current and future Sunyaev-Zeldovich cluster surveys on the f(R) gravity model. We use a Fisher matrix approach, adopt self-calibration for the mass- observable scaling relation, and evaluate constraints for the SPT, Planck, SPTPol and ACTPol surveys. The modified gravity effects on the mass function, halo bias, matter power spectrum, and mass-observable relation are taken into account. We show that, relying on number counts only, the Planck cluster catalog is expected to reduce current upper limits by about a factor of four, to {\sigma}fR0 = 3 {\times} 10-5 (68% confidence level). Adding the cluster power spectrum further improves the constraints to {\sigma}fR0 = 10-5 for SPT and Planck, and {\sigma}fR0 = 3 {\times} 10-6 for SPTPol, pushing cluster constraints significantly beyond the limit where number counts have no constraining power due to the chameleon screening mechanism. Further, the combination of both observables breaks degeneracies, especially with the expansion history (effective dark energy density and equation of state). The constraints are only mildly worsened by the use of self-calibration but depend strongly on the mass threshold of the cluster samples.

This article presents a comparative study of three different types of estimators used for supervised linear unmixing of two MEx/OMEGA hyperspectral cubes. The algorithms take into account the constraints of the abundance fractions, in order to get physically interpretable results. Abundance maps show that the Bayesian maximum a posteriori probability (MAP) estimator proposed in Themelis and Rontogiannis (2008) outperforms the other two schemes, offering a compromise between complexity and estimation performance. Thus, the MAP estimator is a candidate algorithm to perform ice and minerals detection on large hyperspectral datasets.

Weak-gravitational-lensing distortions to the intensity pattern of 21-cm radiation from the dark ages can be decomposed geometrically into curl and curl-free components. Lensing by primordial gravitational waves induces a curl component, while the contribution from lensing by density fluctuations is strongly suppressed. Angular fluctuations in the 21-cm background extend to very small angular scales, and measurements at different frequencies probe different shells in redshift space. There is thus a huge trove of information with which to reconstruct the curl component of the lensing field, allowing tensor-to-scalar ratios conceivably as small as r ~ 10^{-9} - far smaller than those currently accessible - to be probed.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

In this paper we analyze a suite of cosmological simulations of modified gravitational action f(R) models, where cosmic acceleration is induced by a scalar field that acts as a fifth force on all forms of matter. In particular, we focus on the bispectrum of the dark matter density field on mildly non-linear scales. For models with the same initial power spectrum, the dark matter bispectrum shows significant differences for cases where the final dark matter power spectrum also differs. Given the different dependence on bias of the galaxy power spectrum and bispectrum, bispectrum measurements can close the loophole of galaxy bias hiding differences in the power spectrum. Alternatively, changes in the initial power spectrum can also hide differences. By constructing LCDM models with very similar final non-linear power spectra, we show that the differences in the bispectrum are reduced (<4%) and are comparable with differences in the imperfectly matched power spectra. These results indicate that the bispectrum depends mainly on the power spectrum and less sensitively on the gravitational signatures of the f(R) model. This weak dependence of the matter bispectrum on gravity makes it useful for breaking degeneracies associated with galaxy bias, even for models beyond general relativity.

Weak lensing is commonly measured using shear through galaxy ellipticities, or using the effect of magnification bias on galaxy number densities. Here, we report on the first detection of weak lensing magnification with a new, independent technique using the distribution of galaxy sizes and magnitudes. These data come for free in galaxy surveys designed for measuring shear. We present the magnification estimator and apply it to an X-ray selected sample of galaxy groups in the COSMOS HST survey. The measurement of the projected surface density \Sigma(r) is consistent with the shear measurements within the uncertainties, and has roughly 40% of the signal-to-noise of the latter. We discuss systematic issues and challenges to realizing the potential of this new probe of weak lensing.

Recent results of N-body simulations have shown that current theoretical models are not able to correctly predict the amplitude of the scale-dependent halo bias induced by primordial non-Gaussianity, for models going beyond the simplest, local quadratic case. Motivated by these discrepancies, we carefully examine three theoretical approaches based on (1) the statistics of thresholded regions, (2) a peak-background split method based on separation of scales, and (3) a peak-background split method using the conditional mass function. We first demonstrate that the statistics of thresholded regions, which is shown to be equivalent at leading order to a local bias expansion, cannot explain the mass-dependent deviation between theory and N-body simulations. In the two formulations of the peak-background split on the other hand, we identify an important, but previously overlooked, correction to the non-Gaussian bias that strongly depends on halo mass. This new term is in general significant for any primordial non-Gaussianity going beyond the simplest local fNL model. In a separate paper, we compare these new theoretical predictions with N-body simulations, showing good agreement for all simulated types of non-Gaussianity.

The large-scale clustering of galaxies can serve as a probe of primordial non-Gaussianity in the Universe competitive with the anisotropies of the CMB. Here, we present improved theoretical predictions which include an important, previously overlooked correction to the bias. We demonstrate that the new predictions are able to reproduce the results of N-body simulations, explaining the significant departures seen from previous theoretical results. These refined predictions open the way to accurate constraints on primordial physics with large-scale structure surveys.

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.

Many models of inflation predict oscillatory features in the bispectrum of primordial fluctuations. Since it has been shown that primordial non-Gaussianity can lead to a scale-dependent halo bias, we investigate the effect of oscillations in the three-point function on the clustering of dark-matter halos. Interestingly, we find that features in the inflaton potential such as oscillations or sharp steps get imprinted in the mass dependence of the non-Gaussian halo bias. In this paper, we focus on models displaying a sharp feature in the inflaton potential as well as Resonant non-Gaussianity. In both cases, we find a strong scale dependence for the non-Gaussian halo bias with a slope similar to that of the local model. In the resonant case, we find that the non-Gaussian bias oscillates with halo mass, a novel feature that is unique to this type of models. In the case of a sharp feature in the inflaton potential, we find that the clustering of halos is enhanced at the mass scale corresponding to the Fourier mode that exited the horizon when the inflaton was crossing the feature in the potential. Both of these are new effects that open the possibility of characterizing the inflationary potential with large-scale-structure surveys. We briefly discuss the prospects for detecting these non-Gaussian effects.

Ongoing and planned enhancements of the Pierre Auger Observatory

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory

Technical reports on operations and monitoring of the Pierre Auger Observatory

Astrophysical observations of ultra-high-energy cosmic rays with the Pierre Auger Observatory

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers.

The Pierre Auger Collaboration has reported evidence for anisotropy in the distribution of arrival directions of the cosmic rays with energies $E>E_{th}=5.5\times 10^{19}$ eV. These show a correlation with the distribution of nearby extragalactic objects, including an apparent excess around the direction of Centaurus A. If the particles responsible for these excesses at $E>E_{th}$ are heavy nuclei with charge $Z$, the proton component of the sources should lead to excesses in the same regions at energies $E/Z$. We here report the lack of anisotropies in these directions at energies above $E_{th}/Z$ (for illustrative values of $Z=6,\ 13,\ 26$). If the anisotropies above $E_{th}$ are due to nuclei with charge $Z$, and under reasonable assumptions about the acceleration process, these observations imply stringent constraints on the allowed proton fraction at the lower energies.

Estimators for weak lensing observables such as shear and convergence generally have non-linear corrections, which, in principle, make weak lensing power spectra sensitive to primordial non-Gaussianity. In this paper, we quantitatively evaluate these contributions for weak lensing auto- and cross-correlation power spectra, and show that they are strongly suppressed by projection effects. This is a consequence of the central limit theorem, which suppresses departures from Gaussianity when the projection reaches over several correlation lengths of the density field, L_P~55 [Mpc/h]. Furthermore, the typical scales that contribute to projected bispectra are generally smaller than those that contribute to projected power spectra. Both of these effects are not specific to lensing, and thus affect the statistics of non-linear tracers (e.g., peaks) of any projected density field. Thus, the clustering of biased tracers of the three-dimensional density field is generically more sensitive to non-Gaussianity than observables constructed from projected density fields.

We present the results of searches for dipolar-type anisotropies in different energy ranges above $2.5\times 10^{17}$ eV with the surface detector array of the Pierre Auger Observatory, reporting on both the phase and the amplitude measurements of the first harmonic modulation in the right-ascension distribution. Upper limits on the amplitudes are obtained, which provide the most stringent bounds at present, being below 2% at 99% $C.L.$ for EeV energies. We also compare our results to those of previous experiments as well as with some theoretical expectations.

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs "radio-hybrid" measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluoresence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request.

In this study we assess the presence, nature and properties of ices - in particular water ice - that occur within these spots using HIRISE and CRISM observations, as well as the LMD Global Climate Model. Our studies focus on Richardson crater (72{\deg}S, 179{\deg}E) and cover southern spring and summer (LS 175{\deg} - 17 341{\deg}). Three units have been identified of these spots: dark core, gray ring and bright halo. Each unit show characteristic changes as the season progress. In winter, the whole area is covered by CO2 ice with H2O ice contamination. Dark spots form during late winter and early spring. During spring, the dark spots are located in a 10 cm thick depression compared to the surrounding bright ice-rich layer. They are spectrally characterized by weak CO2 ice signatures that probably result from spatial mixing of CO2 ice rich and ice free regions within pixels, and from mixing of surface signatures due to aerosols scattering. The bright halo shaped by winds shows stronger CO2 absorptions than the average ice covered terrain, which is consistent with a formation process involving CO2 re-condensation. According to spectral, morphological and modeling considerations, the gray ring is composed of a thin layer of a few tens of {\mu}m of water ice. Two sources/processes could participate to the enrichment of water ice in the gray ring unit: (i) water ice condensation at the surface in early fall (prior to the condensation of a CO2 rich winter layer) or during winter time (due to cold trapping of the CO2 layer); (ii) ejection of dust grains surrounded by water ice by the geyser activity responsible for the dark spot. In any case, water ice remains longer in the gray ring unit after the complete sublimation of the CO2. Finally, we also looked for liquid water in the near-IR CRISM spectra using linear unmixing modeling but found no conclusive evidence for it.

As one of the most powerful probes of cosmological structure formation, the abundance of massive galaxy clusters is a sensitive probe of modifications to gravity on cosmological scales. In this paper, we present results from N-body simulations of a general class of f(R) models, which self-consistently solve the non-linear field equation for the enhanced forces. Within this class we vary the amplitude of the field, which controls the range of the enhanced gravitational forces, both at the present epoch and as a function of redshift. Most models in the literature can be mapped onto the parameter space of this class. Focusing on the abundance of massive dark matter halos, we compare the simulation results to a simple spherical collapse model. Current constraints lie in the large-field regime, where the chameleon mechanism is not important. In this regime, the spherical collapse model works equally well for a wide range of models and can serve as a model-independent tool for placing constraints on f(R) gravity from cluster abundance. Using these results, we show how constraints from the observed local abundance of X-ray clusters on a specific f(R) model can be mapped onto other members of this general class of models.