The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

The galaxy power spectrum is one of the central quantities in cosmology. It contains information about the primordial inflationary process, the matter clustering, the baryon-photon interaction, the effects of gravity, the galaxy-matter bias, the cosmic expansion, the peculiar velocity field, etc.. Most of this information is however difficult to extract without assuming a specific cosmological model, for instance $\Lambda$CDM and standard gravity. In this paper we explore instead how much information can be obtained that is independent of the cosmological model, both at background and linear perturbation level. We determine the full set of model-independent statistics that can be constructed by combining two redshift bins and two distinct tracers. We focus in particular on the statistics $r(k,z_1,z_2)$, defined as the ratio of $f\sigma_8(z)$ at two redshift shells, and we show how to estimate it with a Fisher matrix approach. Finally, we forecast the constraints on $r$ that can be achieved by future galaxy surveys, and compare it with the standard single-tracer result. We find that $r$ can be measured with a precision from 3 to 11%, depending on the survey. Using two tracers, we find improvements in the constraints up to a factor of two.

We use the improved IllustrisTNG300 magneto-hydrodynamical cosmological simulation to revisit the effect that secondary halo bias has on the clustering of the central galaxy population. With a side length of 205 $h^{-1}$Mpc and significant improvements on the sub-grid model with respect to the previous Illustris boxes, IllustrisTNG300 allows us to explore the dependencies of galaxy clustering over a large cosmological volume and wide halo-mass range. We show, at high statistical significance, that the halo assembly bias signal (i.e., the secondary dependence of halo bias on halo formation redshift) manifests itself on the clustering of the central galaxy population when this is split by stellar mass, colour, specific star formation rate, and surface density. A significant detection is also obtained for galaxy size: at fixed halo mass, larger central galaxies are more tightly clustered than smaller central galaxies in haloes of mass M$_{\rm vir} \lesssim 10^{12.5}$ $h^{-1}$M$_{\odot}$. This effect, however, seems to be uncorrelated with halo formation time, unlike the rest of the secondary dependencies analysed. We also explore the transmission of the halo spin bias signal, i.e., the secondary dependence of halo bias on halo spin. Although galaxy spin retains little information about the total spin of the halo, the correlation is enough to produce a significant galaxy spin bias signal. We discuss possible ways to probe the spin bias effects with observations.

The next generation of galaxy surveys will allow us to test one of the most fundamental assumptions of the standard cosmology, i.e., that gravity is governed by the general theory of relativity (GR). In this paper we investigate the ability of the Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) to constrain GR and its extensions. Based on the J-PAS information on clustering and gravitational lensing, we perform a Fisher matrix forecast on the effective Newton constant, $\mu$, and the gravitational slip parameter, $\eta$, whose deviations from unity would indicate a breakdown of GR. Similar analysis is also performed for the DESI and Euclid surveys and compared to J-PAS with two configurations providing different areas, namely an initial expectation with 4000 $\mathrm{deg}^2$ and the future best case scenario with 8500 $\mathrm{deg}^2$. We show that J-PAS will be able to measure the parameters $\mu$ and $\eta$ at a sensitivity of $2\% - 7\%$, and will provide the best constraints in the interval $z = 0.3 - 0.6$, thanks to the large number of ELGs detectable in that redshift range. We also discuss the constraining power of J-PAS for dark energy models with a time-dependent equation-of-state parameter of the type $w(a)=w_0+w_a(1-a)$, obtaining $\Delta w_0=0.058$ and $\Delta w_a=0.24$ for the absolute errors of the dark energy parameters.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

We use mock galaxy data from the VIMOS Public Extragalactic Redshift Survey (VIPERS) to test the performance of the Multi-Tracer Optimal Estimator (MTOE) of Abramo et al. as a tool to measure the monopoles of the power spectra of multiple tracers of the large-scale structure, $P^{(0)}_\alpha(k)$. We show that MTOE provides more accurate measurements than the standard technique of Feldman, Kaiser & Peacock (FKP), independently of the tracer-selection strategy adopted, on both small and large scales. The largest improvements on individual $P^{(0)}_\alpha(k)$ are obtained using a colour-magnitude selection on small scales, due to MTOE being naturally better equipped to deal with shot noise: we report an average error reduction with respect to FKP of $\sim$ 40$\%$ at $k \, [h$ Mpc$^{-1}]\gtrsim 0.3$. On large scales ($k[h$ Mpc$^{-1}]\lesssim0.1$), the gain in accuracy resulting from cosmic-variance cancellation is $\sim$ 10$\%$ for the ratios of $P^{(0)}_\alpha(k)$. We have carried out a Monte-Carlo Markov Chain analysis to determine the impact of these gains on several quantities derived from $P^{(0)}_\alpha(k)$. If we push the measurement to scales $0.3 < k \, [h$ Mpc$^{-1}]< 0.5$, the average improvements are $\sim$ 30 $\%$ for the amplitudes of the monopoles, $\sim$ 70 $\%$ for the monopole ratios, and $\sim$ 20 $\%$ for the galaxy biases. Our results highlight the potential of MTOE to shed light upon the physics that operate both on large and small cosmological scales. The effect of MTOE on cosmological constraints using VIPERS data will be addressed in a separate paper.

In this paper we present a thorough discussion about the photometric redshift (photo-z) performance of the Southern Photometric Local Universe Survey (S-PLUS). This survey combines a 7 narrow + 5 broad passband filter system, with a typical photometric-depth of r$\sim$21 AB. For this exercise, we utilize the Data Release 1 (DR1), corresponding to 336 deg$^{2}$ from the Stripe-82 region. We rely on the \texttt{BPZ2} code to compute our estimates, using a new library of SED models, which includes additional templates for quiescent galaxies. When compared to a spectroscopic redshift control sample of $\sim$100k galaxies, we find a precision of $\sigma_{z}<$0.8\%, $<$2.0\% or $<$3.0\% for galaxies with magnitudes r$<$17, $<$19 and $<$21, respectively. A precision of 0.6\% is attained for galaxies with the highest \texttt{Odds} values. These estimates have a negligible bias and a fraction of catastrophic outliers inferior to 1\%. We identify a redshift window (i.e., 0.26$<z<$0.32) where our estimates double their precision, due to the simultaneous detection of two emission-lines in two distinct narrow-bands; representing a window opportunity to conduct statistical studies such as luminosity functions. We forecast a total of $\sim$2M, $\sim$16M and $\sim$32M galaxies in the S-PLUS survey with a photo-z precision of $\sigma_{z}<$1.0\%, $<$2.0\% and $<$2.5\% after observing 8000 $deg^{2}$. We also derive redshift Probability Density Functions, proving their reliability encoding redshift uncertainties and their potential recovering the $n(z)$ of galaxies at $z<0.4$, with an unprecedented precision for a photometric survey in the southern hemisphere.

We develop a new, simple, fast and parameter-free method to construct dark matter halo catalogues. This method requires as inputs only the linear matter power spectrum and the threshold density for halo formation in linear theory. It directly uses excursion set ideas and Lagrangian perturbation theory to produce halo catalogues with the correct abundance, large scale power spectrum, bispectrum and velocity field. These halo catalogues can be used for the fast construction of mock galaxy catalogues, allowing for the evaluation of covariance matrices for multiple observables. Because of its robustness and predictive nature, this method can be easily adapted to produce catalogues with e.g. primordial non-Gaussianities, modified theories of gravity and non-standard dark energy models, enabling detailed studies of these models in the context of next-generation surveys. We implement this method in a C code, and present numerical comparisons with theoretical predictions as well as full N-body cosmological simulations.

We estimate the constraining power of J-PAS for parameters of an interacting dark energy cosmology. The survey is expected to map several millions of luminous red galaxies, emission line galaxies and quasars in an area of thousands of square degrees in the northern sky with precise photometric redshift measurements. Forecasts for the DESI and Euclid surveys are also evaluated and compared to J-PAS. With the Fisher matrix approach, we find that J-PAS can place constraints on the interaction parameter comparable to those from DESI, with an absolute uncertainty of about $0.02$, when the interaction term is proportional to the dark matter energy density, and almost as good, of about $0.01$, when the interaction is proportional to the dark energy density. For the equation of state of dark energy, the constraints from J-PAS are slightly better in the two cases (uncertainties $0.04$ - $0.05$ against $0.05$ - $0.07$ around the fiducial value $-1$). Both surveys stay behind Euclid but follow it closely, imposing comparable constraints in all specific cases considered.

We show how to obtain constraints on $\beta=f/b$, the ratio of the matter growth rate and the bias that quantifies the linear redshift-space distortions, that are independent of the cosmological model, using multiple tracers of large-scale structure. For a single tracer the uncertainties on $\beta$ are constrained by the uncertainties in the amplitude and shape of the power spectrum, which is limited by cosmic variance. However, for two or more tracers this limit does not apply, since taking the ratio of power spectra cosmic variance cancels out, and in the linear (Kaiser) approximation one measures directly the quantity $(1+ \beta_1 \mu^2)^2/(1+ \beta_2 \mu^2)^2$, where $\mu$ is the angle of a given mode with the line of sight. We provide analytic formulae for the Fisher matrix for one and two tracers, and quantify the signal-to-noise ratio needed to make effective use of the multiple-tracer technique. We also forecast the errors on $\beta$ for a survey like Euclid.

Halo bias is the main link between the matter distribution and dark matter halos. In its simplest form, halo bias is determined by halo mass, but there are known additional dependencies on other halo properties which are of consequence for accurate modeling of galaxy clustering. Here we present the most precise measurement of these secondary-bias dependencies on halo age, concentration, and spin, for a wide range of halo masses spanning from 10$^{10.7}$ to 10$^{14.7}$ $h^{-1}$ M$_{\odot}$. At the high-mass end, we find no strong evidence of assembly bias for masses above M$_{vir}$ $\sim10^{14}$ $h^{-1}$ M$_{\odot}$. Secondary bias exists, however, for halo concentration and spin, up to cluster-size halos, in agreement with previous findings. For halo spin, we report, for the first time, two different regimes: above M$_{vir}\sim$10$^{11.5}$ $h^{-1}$ M$_{\odot}$, halos with larger values of spin have larger bias, at fixed mass, with the effect reaching almost a factor 2. This trend reverses below this characteristic mass. In addition to these results, we test, for the first time, the performance of a multi-tracer method for the determination of the relative bias between different subsets of halos. We show that this method increases significantly the signal-to-noise of the secondary-bias measurement as compared to a traditional approach. This analysis serves as the basis for follow-up applications of our multi-tracer method to real data.

We investigate the relationship between the turnaround radius ($R_t$) and the virial mass ($M_v$) of cosmic structures in the context of $\Lambda$CDM model and in an $f(R)$ model of modified gravity -- namely, the Hu-Sawicki model. The $R_t$ is the distance from the center of the cosmic structure to the shell that is detaching from the Hubble flow at a given time, while the $M_v$ is defined, for this work, as the mass enclosed within the volume where the density is $200$ times the background density. We consider that gravitationally bound astrophysical systems follow a Navarro-Frenk-White density profile, while beyond the virial radius ($R_v$) the profile is approximated by the 2-halo term of the matter correlation function. By combining them together with the information drawn from solving the spherical collapse for the structures, we are able to connect two observables: the $R_t$ and the $M_v$. We show that, in $\Lambda$CDM, the turnaround mass ($M_t$) at $z=0$ is related to the $M_v$ of that same structure by $M_t \simeq 3.07 \, M_v$, while in terms of the radii we have that $R_t \simeq 3.7 \, R_v$ (for $M_v$ of $10^{13} \, h^{-1} \, M_\odot$). In the $f(R)$ model, we have $M_t \simeq 3.43 \, M_v$ and $R_t \simeq 4.1 \, R_v$, for $|f_{R0}|=10^{-6}$ and the same mass scale. Therefore, the difference between $\Lambda$CDM and $f(R)$ in terms of these observable relations is of order $\sim 10-20\%$ even for a relatively mild strength of the modification of gravity ($|f_{R0}|=10^{-6}$). For the $R_t$ itself we find a difference of $\sim 9\%$ between the weakly modification in gravity considered in this work ($|f_{R0}|=10^{-6}$) and $\Lambda$CDM for a mass of $10^{13} \, h^{-1} \, M_\odot$. Once observations allow precisions of this order or better in measurements $R_t$, as well as the $M_v$, these quantities will become powerful tests of modified gravity.

We investigate the turnaround radius in the spherical collapse model, both in General Relativity and in modified gravity, in particular $f(R)$ scenarios. The phases of spherical collapse are marked by the density contrast in the instant of turnaround $\delta_t$, and by the linear density contrast in the moment of collapse, $\delta_c$. We find that the effective mass of the extra scalar degree of freedom which arises in modified gravity models has an impact on $\delta_t$ of up to $\sim10\%$, and that $\delta_c$ can increase by $\sim1.0\%$. We also compute the turnaround radius, $R_t$, which in modified gravity models can increase by up to $\sim 6\%$ at $z \simeq 0$.

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag.

We present the results of prompt optical follow-up of the electromagnetic counterpart of the gravitational-wave event GW170817 by the Transient Optical Robotic Observatory of the South Collaboration (TOROS). We detected highly significant dimming in the light curves of the counterpart (Delta g=0.17+-0.03 mag, Delta r=0.14+-0.02 mag, Delta i=0.10 +- 0.03 mag) over the course of only 80 minutes of observations obtained ~35 hr after the trigger with the T80-South telescope. A second epoch of observations, obtained ~59 hr after the event with the EABA 1.5m telescope, confirms the fast fading nature of the transient. The observed colors of the counterpart suggest that this event was a "blue kilonova" relatively free of lanthanides.

We present ELDAR, a new method that exploits the potential of medium- and narrow-band filter surveys to securely identify active galactic nuclei (AGN) and determine their redshifts. Our methodology improves on traditional approaches by looking for AGN emission lines expected to be identified against the continuum, thanks to the width of the filters. To assess its performance, we apply ELDAR to the data of the ALHAMBRA survey, which covered an effective area of $2.38\,{\rm deg}^2$ with 20 contiguous medium-band optical filters down to F814W$\simeq 24.5$. Using two different configurations of ELDAR in which we require the detection of at least 2 and 3 emission lines, respectively, we extract two catalogues of type-I AGN. The first is composed of 585 sources ($79\,\%$ of them spectroscopically-unknown) down to F814W$=22.5$ at $z_{\rm phot}>1$, which corresponds to a surface density of $209\,{\rm deg}^{-2}$. In the second, the 494 selected sources ($83\,\%$ of them spectroscopically-unknown) reach F814W$=23$ at $z_{\rm phot}>1.5$, for a corresponding number density of $176\,{\rm deg}^{-2}$. Then, using samples of spectroscopically-known AGN in the ALHAMBRA fields, for the two catalogues we estimate a completeness of $73\,\%$ and $67\,\%$, and a redshift precision of $1.01\,\%$ and $0.86\,\%$ (with outliers fractions of $8.1\,\%$ and $5.8\,\%$). At $z>2$, where our selection performs best, we reach $85\,\%$ and $77\,\%$ completeness and we find no contamination from galaxies.

We study the detectability of large-scale velocity effects on galaxy clustering, by simulating galaxy surveys and combining the clustering of different types of tracers of large-scale structure. We employ a set of lognormal mocks that simulate a $20.000$ deg$^2$ near-complete survey up to $z=0.8$, in which each galaxy mock traces the spatial distribution of dark matter of that mock with a realistic bias prescription. We find that the ratios of the monopoles of the power spectra of different types of tracers carry most of the information that can be extracted from a multi-tracer analysis. In particular, we show that with a multi-tracer technique it will be possible to detect velocity effects with $\gtrsim 3 \sigma$. Finally, we investigate the degeneracy of these effects with the (local) non-Gaussianity parameter $f_{\rm NL}$, and how large-scale velocity contributions could be mistaken for the signatures of primordial non-Gaussianity.

In the Halo Model, galaxies are hosted by dark matter halos, while the halos themselves are biased tracers of the underlying matter distribution. Measurements of galaxy correlation functions include contributions both from galaxies in different halos, and from galaxies in the same halo (the so-called 1-halo terms). We show that, for highly biased tracers, the 1-halo term of the power spectrum obeys a steep scaling relation in terms of bias. We also show that the 1-halo term of the trispectrum has a steep scaling with bias. The steepness of these scaling relations is such that the 1-halo terms can become key contributions to the $n$-point correlation functions, even at large scales. We interpret these results through analytical arguments and semi-analytical calculations in terms of the statistical properties of halos.

We present optimal quadratic estimators for the Fourier analysis of cosmological surveys that detect several different types of tracers of large-scale structure. Our estimators can be used to simultaneously fit the matter power spectrum and the biases of the tracers - as well as redshift-space distortions (RSDs), non-Gaussianities (NGs), or any other effects that are manifested through differences between the clusterings of distinct species of tracers. Our estimators reduce to the one by Feldman, Kaiser & Peacock (ApJ 1994, FKP) in the case of a survey consisting of a single species of tracer. We show that the multi-tracer estimators are unbiased, and that their covariance is given by the inverse of the multi-tracer Fisher matrix (Abramo, MNRAS 2013; Abramo & Leonard, MNRAS 2013). When the biases, RSDs and NGs are fixed to their fiducial values, and one is only interested in measuring the underlying power spectrum, our estimators are projected into the estimator found by Percival, Verde & Peacock (MNRAS 2003). We have tested our estimators on simple (lognormal) simulated galaxy maps, and we show that it performs as expected, being either equivalent or superior to the FKP method in all cases we analyzed. Finally, we have shown how to extend the multi-tracer technique to include the 1-halo term of the power spectrum.

We study the characteristics of a narrow band type Ia supernova survey through simulations based on the upcoming Javalambre Physics of the accelerating universe Astrophysical Survey (J-PAS). This unique survey has the capabilities of obtaining distances, redshifts, and the SN type from a single experiment thereby circumventing the challenges faced by the resource-intensive spectroscopic follow-up observations. We analyse the flux measurements signal-to-noise ratio and bias, the supernova typing performance, the ability to recover light curve parameters given by the SALT2 model, the photometric redshift precision from type Ia supernova light curves and the effects of systematic errors on the data. We show that such a survey is not only feasible but may yield large type Ia supernova samples (up to 250 supernovae at $z<0.5$ per month of search) with low core collapse contamination ($\sim 1.5$ per cent), good precision on the SALT2 parameters (average $\sigma_{m_B}=0.063$, $\sigma_{x_1}=0.47$ and $\sigma_c=0.040$) and on the distance modulus (average $\sigma_{\mu}=0.16$, assuming an intrinsic scatter $\sigma_{\mathrm{int}}=0.14$), with identified systematic uncertainties $\sigma_{\mathrm{sys}}\lesssim 0.10 \sigma_{\mathrm{stat}}$. Moreover, the filters are narrow enough to detect most spectral features and obtain excellent photometric redshift precision of $\sigma_z=0.005$, apart from $\sim$ 2 per cent of outliers. We also present a few strategies for optimising the survey's outcome. Together with the detailed host galaxy information, narrow band surveys can be very valuable for the study of supernova rates, spectral feature relations, intrinsic colour variations and correlations between supernova and host galaxy properties, all of which are important information for supernova cosmological applications.

The Javalambre-Physics of the Accelerated Universe Astrophysical Survey (J-PAS) is a narrow band, very wide field Cosmological Survey to be carried out from the Javalambre Observatory in Spain with a purpose-built, dedicated 2.5m telescope and a 4.7 sq.deg. camera with 1.2Gpix. Starting in late 2015, J-PAS will observe 8500sq.deg. of Northern Sky and measure $0.003(1+z)$ photo-z for $9\times10^7$ LRG and ELG galaxies plus several million QSOs, sampling an effective volume of $\sim 14$ Gpc$^3$ up to $z=1.3$ and becoming the first radial BAO experiment to reach Stage IV. J-PAS will detect $7\times 10^5$ galaxy clusters and groups, setting constrains on Dark Energy which rival those obtained from its BAO measurements. Thanks to the superb characteristics of the site (seeing ~0.7 arcsec), J-PAS is expected to obtain a deep, sub-arcsec image of the Northern sky, which combined with its unique photo-z precision will produce one of the most powerful cosmological lensing surveys before the arrival of Euclid. J-PAS unprecedented spectral time domain information will enable a self-contained SN survey that, without the need for external spectroscopic follow-up, will detect, classify and measure $\sigma_z\sim 0.5\%$ redshifts for $\sim 4000$ SNeIa and $\sim 900$ core-collapse SNe. The key to the J-PAS potential is its innovative approach: a contiguous system of 54 filters with $145\AA$ width, placed $100\AA$ apart over a multi-degree FoV is a powerful "redshift machine", with the survey speed of a 4000 multiplexing low resolution spectrograph, but many times cheaper and much faster to build. The J-PAS camera is equivalent to a 4.7 sq.deg. "IFU" and it will produce a time-resolved, 3D image of the Northern Sky with a very wide range of Astrophysical applications in Galaxy Evolution, the nearby Universe and the study of resolved stellar populations.

The fluctuations in the temperature and polarization of the cosmic microwave background are described by a hierarchy of Boltzmann equations. In its integral form, this Boltzmann hierarchy can be converted from the usual Fourier-space base into a position-space and causal description. We show that probability densities for random flights play a key role in this description. The integral system can be treated as a perturbative series in the number of steps of the random flights, and the properties of random flight probabilities impose constraints on the domains of dependence. We show that, as a result of these domains, a Fourier-Bessel decomposition can be employed in order to calculate the random flight probability densities. We also illustrate how the H-theorem applies to the cosmic microwave background: by using analytical formulae for the asymptotic limits of these probability densities, we show that, as the photon distribution approaches a state of equilibrium, both the temperature anisotropies and the net polarization must vanish.

We present a formal derivation of the key equations governing gravitational lensing in arbitrary space-times, starting from the basic properties of Jacobi fields and their expressions in terms of the exponential map. A careful analysis of Jacobi fields and Jacobi classes near the origin of a light beam determines the nature of the singular behavior of the optical deformation matrix. We also show that potential problems that could arise from this singularity do not invalidate the conclusions of the original argument presented by Seitz, Schneider & Ehlers (1994).

Galaxy surveys that map multiple species of tracers of large-scale structure can improve the constraints on some cosmological parameters far beyond the limits imposed by a simplistic interpretation of cosmic variance. This enhancement derives from comparing the relative clustering between different tracers of large-scale structure. We present a simple but fully generic expression for the Fisher information matrix of surveys with any (discrete) number of tracers, and show that the enhancement of the constraints on bias-sensitive parameters are a straightforward consequence of this multi-tracer Fisher matrix. In fact, the relative clustering amplitudes between tracers are eigenvectors of this multi-tracer Fisher matrix. The diagonalized multi-tracer Fisher matrix clearly shows that while the effective volume is bounded by the physical volume of the survey, the relational information between species is unbounded. As an application, we study the expected enhancements in the constraints of realistic surveys that aim at mapping several different types of tracers of large-scale structure. The gain obtained by combining multiple tracers is highest at low redshifts, and in one particular scenario we analyzed, the enhancement can be as large as a factor of ~3 for the accuracy in the determination of the redshift distortion parameter, and a factor ~5 for the local non-Gaussianity parameter. Radial and angular distance determinations from the baryonic features in the power spectrum may also benefit from the multi-tracer approach.

We show that a large-area imaging survey using narrow-band filters could detect quasars in sufficiently high number densities, and with more than sufficient accuracy in their photometric redshifts, to turn them into suitable tracers of large-scale structure. If a narrow-band optical survey can detect objects as faint as i=23, it could reach volumetric number densities as high as 10^{-4} h^3 Mpc^{-3} (comoving) at z~1.5 . Such a catalog would lead to precision measurements of the power spectrum up to z~3-4. We also show that it is possible to employ quasars to measure baryon acoustic oscillations at high redshifts, where the uncertainties from redshift distortions and nonlinearities are much smaller than at z<1. As a concrete example we study the future impact of J-PAS, which is a narrow-band imaging survey in the optical over 1/5 of the unobscured sky with 42 filters of ~100 A full-width at half-maximum. We show that J-PAS will be able to take advantage of the broad emission lines of quasars to deliver excellent photometric redshifts, \sigma_{z}~0.002(1+z), for millions of objects.

Starting from the Fisher matrix for counts in cells, I derive the full Fisher matrix for surveys of multiple tracers of large-scale structure. The key assumption is that the inverse of the covariance of the galaxy counts is given by the naive matrix inverse of the covariance in a mixed position-space and Fourier-space basis. I then compute the Fisher matrix for the power spectrum in bins of the three-dimensional wavenumber k; the Fisher matrix for functions of position x (or redshift z) such as the linear bias of the tracers and/or the growth function; and the cross-terms of the Fisher matrix that expresses the correlations between estimations of the power spectrum and estimations of the bias. When the bias and growth function are fully specified, and the Fourier-space bins are large enough that the covariance between them can be neglected, the Fisher matrix for the power spectrum reduces to the widely used result that was first derived by Feldman, Kaiser and Peacock (1994). Assuming isotropy, an exact calculation of the Fisher matrix can be performed in the case of a constant-density, volume-limited survey. I then show how the exact Fisher matrix in the general case can be obtained in terms of a series of volume-limited surveys.

Fossil systems are defined to be X-ray bright galaxy groups with a 2-magnitude difference between their two brightest galaxies within half the projected virial radius,and represent an interesting extreme of the population of galaxy agglomerations.However,the physical conditions and processes leading to their formation are still poorly constrained.We compare the outskirts of fossil systems with that of normal groups to understand whether environmental conditions play a significant role in their formation.We study galaxy groups in both,numerical simulations and observations.We use a variety of statistical tools including the spatial cross-correlation function and the local density parameter \Delta_5 to probe differences in the density and structure of the environments of normal and fossil systems in the Millennium simulation.We find that the number density of galaxies surrounding fossil systems evolves from greater than that observed around normal systems at z=0.69, to lower than the normal systems by z=0.Both fossil and normal systems exhibit an increment in their otherwise radially declining local density measure (\Delta_5) at distances of order 2.5r_{vir} from the system centre.We show that this increment is more noticeable for fossil systems than normal systems and demonstrate that this difference is linked to the earlier formation epoch of fossil groups.Despite the importance of the assembly time, we show that the environment is different for fossil and non-fossil systems with similar masses and formation times along their evolution.We also confirm that the physical characteristics identified in the Millennium simulation can also be detected in SDSS observations.Our results confirm the commonly held belief that fossil systems assembled earlier than normal systems but also show that the surroundings of fossil groups could be responsible for the formation of their large magnitude gap.

We investigate the effect of an interaction between dark energy and dark matter upon the dynamics of galaxy clusters. This effect is computed through the Layser-Irvine equation, which describes how an astrophysical system reaches virial equilibrium and was modified to include the dark interactions. Using observational data from almost 100 purportedly relaxed galaxy clusters we put constraints on the strength of the couplings in the dark sector. We compare our results with those from other observations and find that a positive (in the sense of energy flow from dark energy to dark matter) non vanishing interaction is consistent with the data within several standard deviations.

We review the basic hypotheses which motivate the statistical framework used to analyze the cosmic microwave background, and how that framework can be enlarged as we relax those hypotheses. In particular, we try to separate as much as possible the questions of gaussianity, homogeneity and isotropy from each other. We focus both on isotropic estimators of non-gaussianity as well as statistically anisotropic estimators of gaussianity, giving particular emphasis on their signatures and the enhanced "cosmic variances" that become increasingly important as our putative Universe becomes less symmetric. After reviewing the formalism behind some simple model-independent tests, we discuss how these tests can be applied to CMB data when searching for large scale "anomalies"

This paper makes two points. First, we show that the line-of-sight solution to cosmic microwave anisotropies in Fourier space, even though formally defined for arbitrarily large wavelengths, leads to position-space solutions which only depend on the sources of anisotropies inside the past light-cone of the observer. This happens order by order in a series expansion in powers of the visibility $\gamma=e^{-\mu}$, where $\mu$ is the optical depth to Thompson scattering. We show that the CMB anisotropies are regulated by spacetime window functions which have support only inside the past light-cone of the point of observation. Second, we show that the Fourier-Bessel expansion of the physical fields (including the temperature and polarization momenta) is an alternative to the usual Fourier basis as a framework to compute the anisotropies. In that expansion, for each multipole $l$ there is a discrete tower of momenta $k_{i,l}$ (not a continuum) which can affect physical observables, with the smallest momenta being $k_{1,l} ~ l$. The Fourier-Bessel modes take into account precisely the information from the sources of anisotropies that propagates from the initial value surface to the point of observation - no more, no less. We also show that the physical observables (the temperature and polarization maps), and hence the angular power spectra, are unaffected by that choice of basis. This implies that the Fourier-Bessel expansion is the optimal scheme with which one can compute CMB anisotropies. (Abridged)

We search for planar deviations of statistical isotropy in the Wilkinson Microwave Anisotropy Probe (WMAP) data by applying a recently introduced angular-planar statistics both to full-sky and to masked temperature maps, including in our analysis the effect of the residual foreground contamination and systematics in the foreground removing process as sources of error. We confirm earlier findings that full-sky maps exhibit anomalies at the planar ($l$) and angular ($\ell$) scales $(l,\ell)=(2,5),(4,7),$ and $(6,8)$, which seem to be due to unremoved foregrounds since this features are present in the full-sky map but not in the masked maps. On the other hand, our test detects slightly anomalous results at the scales $(l,\ell)=(10,8)$ and $(2,9)$ in the masked maps but not in the full-sky one, indicating that the foreground cleaning procedure (used to generate the full-sky map) could not only be creating false anomalies but also hiding existing ones. We also find a significant trace of an anomaly in the full-sky map at the scale $(l,\ell)=(10,5)$, which is still present when we consider galactic cuts of 18.3% and 28.4%. As regards the quadrupole ($\ell=2$), we find a coherent over-modulation over the whole celestial sphere, for all full-sky and cut-sky maps. Overall, our results seem to indicate that current CMB maps derived from WMAP data do not show significant signs of anisotropies, as measured by our angular-planar estimator. However, we have detected a curious coherence of planar modulations at angular scales of the order of the galaxy's plane, which may be an indication of residual contaminations in the full- and cut-sky maps.

Gaussianity and statistical isotropy of the Universe are modern cosmology's minimal set of hypotheses. In this work we introduce a new statistical test to detect observational deviations from this minimal set. By defining the temperature correlation function over the whole celestial sphere, we are able to independently quantify both angular and planar dependence (modulations) of the CMB temperature power spectrum over different slices of this sphere. Given that planar dependence leads to further modulations of the usual angular power spectrum $C_l$, this test can potentially reveal richer structures in the morphology of the primordial temperature field. We have also constructed an unbiased estimator for this angular-planar power spectrum which naturally generalizes the estimator for the usual $C_l$'s. With the help of a chi-square analysis, we have used this estimator to search for observational deviations of statistical isotropy in WMAP's 5 year release data set (ILC5), where we found only slight anomalies on the angular scales $l=7$ and $l=8$. Since this angular-planar statistic is model-independent, it is ideal to employ in searches of statistical anisotropy (e.g., contaminations from the galactic plane) and to characterize non-Gaussianities.

We investigate the influence of an interaction between dark energy and dark matter upon the dynamics of galaxy clusters. We obtain the general Layser-Irvine equation in the presence of interactions, and find how, in that case, the virial theorem stands corrected. Using optical, X-ray and weak lensing data from 33 relaxed galaxy clusters, we put constraints on the strength of the coupling in the dark sector. We find that this coupling is small but positive, indicating that dark energy can be decaying into dark matter.

All models of dynamical dark energy possess fluctuations, which affect the number of galaxy clusters in the Universe. We have studied the impact of dark energy clustering on the number of clusters using a generalization of the spherical collapse model and the Press-Schechter formalism. Our statistical analysis is performed in a 7-parameter space using the Fisher matrix method, for several hypothetical Sunyaev-Zel'dovich and weak lensing (shear maps) surveys. In some scenarios, the impact of these fluctuations is large enough that their effect could already be detected by existing instruments such as the South Pole Telescope, when its data is combined with WMAP and SDSS. Future observations could go much further and probe the nature of dark energy by distinguishing between different models on the basis of their perturbations, not only their expansion histories.

The abundance and distribution of collapsed objects such as galaxy clusters will become an important tool to investigate the nature of dark energy and dark matter. Number counts of very massive objects are sensitive not only to the equation of state of dark energy, which parametrizes the smooth component of its pressure, but also to the sound speed of dark energy as well, which determines the amount of pressure in inhomogeneous and collapsed structures. Since the evolution of these structures must be followed well into the nonlinear regime, and a fully relativistic framework for this regime does not exist yet, we compare two approximate schemes: the widely used spherical collapse model, and the pseudo-Newtonian approach. We show that both approximation schemes convey identical equations for the density contrast, when the pressure perturbation of dark energy is parametrized in terms of an effective sound speed. We also make a comparison of these approximate approaches to general relativity in the linearized regime, which lends some support to the approximations.

We discuss the intriguing possibility that dark energy may change its equation of state in situations where large dark energy fluctuations are present. We show indications of this dynamical mutation in some generic models of dark energy.

We study non-linear structure formation in the presence of dark energy. The influence of dark energy on the growth of large-scale cosmological structures is exerted both through its background effect on the expansion rate, and through its perturbations as well. In order to compute the rate of formation of massive objects we employ the Spherical Collapse formalism, which we generalize to include fluids with pressure. We show that the resulting non-linear evolution equations are identical to the ones obtained in the Pseudo-Newtonian approach to cosmological perturbations, in the regime where an equation of state serves to describe both the background pressure relative to density, and the pressure perturbations relative to the density perturbations as well. We then consider a wide range of constant and time-dependent equations of state (including phantom models) parametrized in a standard way, and study their impact on the non-linear growth of structure. The main effect is the formation of dark energy structure associated with the dark matter halo: non-phantom equations of state induce the formation of a dark energy halo, damping the growth of structures; phantom models, on the other hand, generate dark energy voids, enhancing structure growth. Finally, we employ the Press-Schechter formalism to compute how dark energy affects the number of massive objects as a function of redshift.

By demanding that a bounce is nonsingular and that perturbations are well-behaved at all times, we narrow the scope of possible models with one degree of freedom that can describe a bounce in the absence of spatial curvature. We compute the general properties of the transfer matrix of perturbations through the bounce, and show that spectral distortions of the Bardeen potential $\Phi$ are generically produced only for the small wavelengths, although the spectrum of long wavelength curvature perturbations produced in a contracting phase gets propagated unaffected through such a bounce.

We describe how a survey of the polarization of the cosmic microwave background induced by Compton scattering in galaxy clusters can be used to make a full spatial reconstruction of the primordial (z~1089) matter distribution inside our surface of last scattering. This "polarization tomography" can yield a spatial map of the initial state of the Universe just as gravitational collapse was beginning to drive structure formation. We present a transparent method and simple formulas from which one can compute the 3D primordial map in real and in Fourier space, given a 3D map of the polarization due to galaxy clusters. The advantage of the real space reconstruction is that it is free from the statistical uncertainties which are inherent in the Fourier space reconstruction. We discuss how noise, partial sky covering and depth of the survey can affect the results.

We investigate the large scale anomalies in the angular distribution of the cosmic microwave background radiation as measured by WMAP using several tests. These tests, based on the multipole vector expansion, measure correlations between the phases of the multipoles as expressed by the directions of the multipole vectors and their associated normal planes. We have computed the probability distribution functions for 46 such tests, for the multipoles l=2-5. We confirm earlier findings that point to a high level of alignment between l=2 (quadrupole) and l=3 (octopole), but with our tests we do not find significant planarity in the octopole. In addition, we have found other possible anomalies in the alignment between the octopole and the l=4 (hexadecupole) components, as well as in the planarity of l=4 and l=5. We introduce the notion of a total likelihood to estimate the relevance of the low-multipoles tests of non-gaussianity. We show that, as a result of these tests, the CMB maps which are most widely used for cosmological analysis lie within the ~ 10% of randomly generated maps with lowest likelihoods.

We discuss how an extended foreground of the cosmic microwave background (CMB) can account for the anomalies in the low multipoles of the CMB anisotropies. The distortion needed to account for the anomalies is consistent with a cold spot with the spatial geometry of the Local Supercluster (LSC) and a temperature quadrupole of order DeltaT_2^2 ~ 50 microK^2. If this hypothetic foreground is subtracted from the CMB data, the amplitude of the quadrupole (l=2) is substantially increased, and the statistically improbable alignment of the quadrupole with the octopole (l=3) is substantially weakened, increasing dramatically the likelihood of the "cleaned" maps. By placing the foreground on random locations and then computing the likelihood of the cleaned maps we can estimate the most likely place for this foreground. Although the 1-year WMAP data clearly points the location of this hypothetical foreground to the LSC or its specular image (i.e., the vicinity of the poles of the cosmic dipole axis), the three-year data seems to point to these locations as well as the north ecliptic pole. We show that this is consistent with the symmetries of the cosmic quadrupole. We also discuss a possible mechanism that could have generated this foreground: the thermal Sunyaev-Zeldovich effect caused by hot electrons in the LSC. We argue that the temperature and density of the hot gas which are necessary to generate such an effect, though in the upper end of the expected range of values, are marginally consistent with present observations of the X-ray background of spectral distortions of the CMB.

We show that phantom dark energy, if it is described by a K-essence theory, has three fundamental problems: first, its hamiltonian is unbounded from below. Second, classical stability precludes the equation of state from crossing the ``Lambda-barrier'', $w_\Lambda=-1$. Finally, both the equation of state and the sound speed are unbounded -- the first, from below, the second, from above -- if the kinetic term is not bounded by dynamics.

We study the CMB constraints on two Dark Energy models described by scalar fields with different Lagrangians, namely a Klein-Gordon and a Born-Infeld field. The speed of sound of field fluctuations are different in these two theories, and therefore the predictions for CMB and structure formation are different. Employing the WMAP data on CMB, we make a likelihood analysis on a grid of theoretical models. We constrain the parameters of the models and compute the probability distribution functions for the equation of state. We show that the effect of the different sound speeds affects the low multipoles of CMB anisotropies, but is at most marginal for the class of models studied here.

We show that the thermal Sunyaev-Zeldovich effect caused by hot electrons in the Local Supercluster (LSC) can explain the abnormal quadrupole and octopole of the cosmic microwave background (CMB) that were measured by WMAP and COBE. The distortion needed to account for the low observed quadrupole is a spot in the direction of the LSC with a temperature decrease of order \Delta T \approx - 7 \mu K for \nu ~ 20 -- 90 Ghz photons. The temperature and density of the hot gas which can generate such an effect are consistent with observations of the X-ray background. If this hypothetic foreground is subtracted from the WMAP data, we find that the amplitude of the quadrupole (l=2) is substantially increased, and that the ``planarity'' of both the quadrupole and the octopole (l=3) are weakened. For smaller scales the effect decays and, at least in our simplified model, it does not affect the angular power spectrum at l>10. Moreover, since the Sunyaev-Zeldovich effect increases the temperature of photons with frequencies above 218 GHz, observations sensitive in that range (such as PLANCK's HFI) will be able to confirm whether the LSC indeed affects the CMB.

We investigate tachyon dynamics with an inverse power-law potential $V(\phi) \propto \phi^{-\alpha}$. We find global attractors of the dynamics leading to a dust behaviour for $\alpha > 2$ and to an accelerating universe for $0 < \alpha \le 2$. We study linear cosmological perturbations and we show that metric fluctuations are constant on large scales in both cases. In presence of an additional perfect fluid, the tachyon with this potential behaves as dust or dark energy.

We show that Einstein's gravity coupled to a non-minimally coupled scalar field is stable even for high values of the scalar field, when the sign of the Einstein-Hilbert action is reversed. We also discuss inflationary solutions and a possible new mechanism of reheating.

We demonstrate the existence of a secular back-reaction on inflation using a simple scalar model. The model consists of a massless, minimally coupled scalar with a quartic self-interaction which is a spectator to $\Lambda$-driven inflation. To avoid problems with coincident propagators, and to make the scalars interact more like gravitons, we impose a covariant normal ordering prescription which has the effect of removing tadpole graphs. This version of the theory exhibits a secular slowing at three loop order due to interactions between virtual infrared scalars which are ripped apart by the inflating background. The effect is quantified using an invariant observable and all orders bounds are given. We also argue that, although stochastic effects can have either sign, the slowing mechanism is superimposed upon them.

We define a scalar measure of the local expansion rate based on how astronomers determine the Hubble constant. Our observable is the inverse conformal d'Alembertian acting on a unit ``standard candle.'' Because this quantity is an integral over the past lightcone of the observation point it provides a manifestly causal and covariant technique for averaging over small fluctuations. For an exactly homogeneous and isotropic spacetime our scalar gives minus one half times the inverse square of the Hubble parameter. Our proposal is that it be assigned this meaning generally and that it be employed to decide the issue of whether or not there is a significant quantum gravitational back-reaction on inflation. Several techniques are discussed for promoting the scalar to a full invariant by giving a geometrical description for the point of observation. We work out an explicit formalism for evaluating the invariant in perturbation theory. The results for two simple models are presented in subsequent papers.

We use an invariant operator to study the quantum gravitational back-reaction to scalar perturbations during chaotic inflation. Our operator is the inverse covariant d'Alembertian expressed as a function of the local value of the inflaton. In the slow roll approximation this observable gives $-1/(2 H^2)$ for an arbitrary homogeneous and isotropic geometry, hence it is a good candidate for measuring the local expansion rate even when the spacetime is not perfectly homogeneous and isotropic. Corrections quadratic in the scalar creation and annihilation operators of the initial value surface are included using the slow-roll and long wavelength approximations. The result is that all terms which could produce a significant secular back-reaction cancel from the operator, before one even takes its expectation value. Although it is not relevant to the current study, we also develop a formalism for using stochastic samples to study back-reaction.

We investigate the evolution of cosmological perturbations in scenarios with a quintessence scalar field, both analytically and numerically. In the tracking regime for quintessence, we find the long wavelength solutions for quintessential perturbations. We discuss the possibility of isocurvature modes generated by the quintessence sector and their impact on observations.

Inflation leads us to expect a spectrum of gravitational waves (tensor perturbations) extending to wavelengths much bigger than the present observable horizon. Although these gravity waves are not directly observable, the energy density that they contribute grows in importance during the radiation- and dust-dominated ages of the universe. We show that the back reaction of tensor perturbations during matter domination is limited from above, since gravitational waves of wavelength $\lambda$ have a share of the total energy density $\Delta \rho(\lambda)/\rho$ during matter domination that is at most equal to the share of the total energy density that they had when the mode $\lambda$ exited the Hubble radius $H^{-1}$ during inflation. This work is to be contrasted to that of Sahni, who analyzed the energy density of gravity waves only insofar as their wavelengths are smaller than $H^{-1}$. Such a cut-off in the spectral energy of gravity waves leads to the breakdown of energy conservation, and we show that this anomaly is eliminated simply by taking into account the energy density and pressure of long wavelength gravitational waves as well as short wavelength ones.

We consider quantum mechanical corrections to a homogeneous, isotropic and spatially flat geometry whose scale factor expands classically as a general power of the co-moving time. The effects of both gravitons and the scalar inflaton are computed at one loop using the manifestly causal formalism of Schwinger with the Feynman rules recently developed by Iliopoulos {\it et al.} We find no significant effect, in marked contrast with the result obtained by Mukhanov {\it et al.} for chaotic inflation based on a quadratic potential. By applying the canonical technique of Mukhanov {\it et al.} to the exponential potentials of power law inflation, we show that the two methods produce the same results, within the approximations employed, for these backgrounds. We therefore conclude that the shape of the inflaton potential can have an enormous impact on the one loop back-reaction.

We extend, for the case of a general scalar potential, the inflaton-graviton Feynman rules recently developed by Iliopoulos {\it et al.} As an application we compute the leading term, for late co-moving times, of the one loop back reaction on the expansion rate for $V(\phi) = \frac12 m^2 \phi^2$. This is expressed as the logarithmic time derivative of the scale factor in the coordinate system for which the expectation value of the metric has the form: $<0 | g_{\mu\nu}({\bar t},{\vec x}) | 0 > dx^{\mu} dx^{\nu} = - d{\bar t}^2 + a^2({\bar t}) d{\vec x} \cdot d{\vec x}$. This quantity should be a gauge independent observable. Our result for it agrees exactly with that inferred from the effect previously computed by Mukhanov {\it et al.} using canonical quantization. It is significant that the two calculations were made with completely different schemes for fixing the gauge, and that our computation was done using the standard formalism of covariant quantization. This should settle some of the issues recently raised by Unruh.

We derive the implications for anisotropies in the cosmic microwave background following from a model of inflation in which a bare cosmological constant is gradually screened by an infrared process in quantum gravity. The model predicts that the amplitude of scalar perturbations is $A_S = (2.0 \pm .2) \times 10^{-5}$, that the tensor-to-scalar ratio is $r \approx 1.7 \times 10^{-3}$, and that the scalar and tensor spectral indices are $n \approx .97$ and $n_T \approx -2.8 \times 10^{-4}$, respectively. By comparing the model's power spectrum with the COBE 4-year RMS quadrupole, the mass scale of inflation is determined to be $M = (.72 \pm .03) \times 10^{16}~{\rm GeV}$. At this scale the model produces about $10^8$ e-foldings of inflation, so another prediction is $\Omega = 1$.

We study the effective energy-momentum tensor (EMT) for cosmological perturbations and formulate the gravitational back-reaction problem in a gauge invariant manner. We analyze the explicit expressions for the EMT in the cases of scalar metric fluctuations and of gravitational waves and derive the resulting equations of state. The formalism is applied to investigate the back-reaction effects in chaotic inflation. We find that for long wavelength scalar and tensor perturbations, the effective energy density is negative and thus counteracts any pre-existing cosmological constant. For scalar perturbations during an epoch of inflation, the equation of state is de Sitter-like.

The back reaction of gravitational perturbations in a homogeneous background is determined by an effective energy-momentum tensor quadratic in the perturbations. We show that this nonlinear feedback effect is important in the case of long wavelength scalar perturbations in inflationary universe models. We also show how to solve an old problem concerning the gauge dependence of the effective energy-momentum tensor of perturbations.

We derive the effective energy-momentum tensor for cosmological perturbations and prove its gauge-invariance. The result is applied to study the influence of perturbations on the behaviour of the Friedmann background in inflationary Universe scenarios. We found that the back reaction of cosmological perturbations on the background can become important already at energies below the self-reproduction scale.