
A significant challenge facing photometric surveys for cosmological purposes
is the need to produce reliable redshift estimates. The estimation of
photometric redshifts (photozs) has been consolidated as the standard strategy
to bypass the high production costs and incompleteness of spectroscopic
redshift samples. Trainingbased photoz methods require the preparation of a
highquality list of spectroscopic redshifts, which needs to be constantly
updated. The photoz training, validation, and estimation must be performed in
a consistent and reproducible way in order to accomplish the scientific
requirements. To meet this purpose, we developed an integrated webbased data
interface that not only provides the framework to carry out the above steps in
a systematic way, enabling the ease testing and comparison of different
algorithms, but also addresses the processing requirements by parallelizing the
calculation in a transparent way for the user. This framework called the
Science Portal (hereafter Portal) was developed in the context the Dark Energy
Survey (DES) to facilitate scientific analysis. In this paper, we show how the
Portal can provide a reliable environment to access vast data sets, provide
validation algorithms and metrics, even in the case of multiple photozs
methods. It is possible to maintain the provenance between the steps of a chain
of workflows while ensuring reproducibility of the results. We illustrate how
the Portal can be used to provide photoz estimates using the DES first year
(Y1A1) data. While the DES collaboration is still developing techniques to
obtain more precise photozs, having a structured framework like the one
presented here is critical for the systematic vetting of DES algorithmic
improvements and the consistent production of photozs in the future DES
releases.

We present a novel approach for creating scienceready catalogs through a
software infrastructure developed for the Dark Energy Survey (DES). We
integrate the data products released by the DES Data Management and additional
products created by the DES collaboration in an environment known as DES
Science Portal. Each step involved in the creation of a scienceready catalog
is recorded in a relational database and can be recovered at any time. We
describe how the DES Science Portal automates the creation and characterization
of lightweight catalogs for DES Year 1 Annual Release, and show its flexibility
in creating multiple catalogs with different inputs and configurations.
Finally, we discuss the advantages of this infrastructure for large surveys
such as DES and the Large Synoptic Survey Telescope. The capability of creating
scienceready catalogs efficiently and with full control of the inputs and
configurations used is an important asset for supporting science analysis using
data from large astronomical surveys.

Approximate methods to full Nbody simulations provide a fast and accurate
solution to the development of mock catalogues for the modeling of galaxy
clustering observables. In this paper we extend ICECOLA (Izard et al. 2016),
based on an optimized implementation of the approximate COLA method, to produce
weak lensing maps and halo catalogues in the light cone using an integrated and
self consistent approach. We show that despite the approximate dynamics, the
catalogues thus produced enable an accurate modeling of weak lensing
observables one decade beyond the characteristic scale where the growth becomes
nonlinear. In particular, we compare ICECOLA to the MICEGC Nbody simulation
for some fiducial cases representative of upcoming surveys and find that, for
sources at redshift $z=1$, their convergence power spectra agree to within one
percent up to high multipoles (i.e., of order $1000$). The corresponding shear
two point functions, $\xi_{+}$ and $\xi_{}$, yield similar accuracy down to
$2$ and $20$ arcmin respectively, while tangential shear around a $z=0.5$ lens
sample is accurate down to $4$ arcmin. We show that such accuracy is stable
against an increased angular resolution of the weak lensing maps. Hence, this
opens the possibility of using approximate methods for the joint modeling of
galaxy clustering and weak lensing observables and their covariance in ongoing
and future galaxy surveys.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

We present the discovery and preliminary characterization of a
gravitationally lensed quasar with a source redshift $z_{s}=2.74$ and image
separation of $2.9"$ lensed by a foreground $z_{l}=0.40$ elliptical galaxy.
Since the images of gravitationally lensed quasars are the superposition of
multiple point sources and a foreground lensing galaxy, we have developed a
morphology independent multiwavelength approach to the photometric selection
of lensed quasar candidates based on Gaussian Mixture Models (GMM) supervised
machine learning. Using this technique and $gi$ multicolour photometric
observations from the Dark Energy Survey (DES), near IR $JK$ photometry from
the VISTA Hemisphere Survey (VHS) and WISE mid IR photometry, we have
identified a candidate system with two catalogue components with $i_{AB}=18.61$
and $i_{AB}=20.44$ comprised of an elliptical galaxy and two blue point
sources. Spectroscopic followup with NTT and the use of an archival AAT
spectrum show that the point sources can be identified as a lensed quasar with
an emission line redshift of $z=2.739\pm0.003$ and a foreground early type
galaxy with $z=0.400\pm0.002$. We model the system as a single isothermal
ellipsoid and find the Einstein radius $\theta_E \sim 1.47"$, enclosed mass
$M_{enc} \sim 4 \times 10^{11}$M$_{\odot}$ and a time delay of $\sim$52 days.
The relatively wide separation, month scale time delay duration and high
redshift make this an ideal system for constraining the expansion rate beyond a
redshift of 1.

We present a new method to measure the redshiftdependent galaxy bias by
combining information from the galaxy density field and the weak lensing field.
This method is based on Amara et al. (2012), where they use the galaxy density
field to construct a biasweighted convergence field kg. The main difference
between Amara et al. (2012) and our new implementation is that here we present
another way to measure galaxy bias using tomography instead of bias
parameterizations. The correlation between kg and the true lensing field k
allows us to measure galaxy bias using different zerolag correlations, such as
<kgk>/<kk> or <kgkg>/<kgk>. Our method measures the linear bias factor on
linear scales under the assumption of no stochasticity between galaxies and
matter. We use the MICE simulation to measure the linear galaxy bias for a
fluxlimited sample (i < 22.5) in tomographic redshift bins using this method.
This paper is the first that studies the accuracy and systematic uncertainties
associated with the implementation of the method, and the regime where it is
consistent with the linear galaxy bias defined by projected 2point correlation
functions (2PCF). We find that our method is consistent with linear bias at the
percent level for scales larger than 30 arcmin, while nonlinearities appear at
smaller scales. This measurement is a good complement to other measurements of
bias, since it does not depend strongly on sigma8 as the 2PCF measurements. We
apply this method to the Dark Energy Survey Science Verification data in a
followup paper.

Host galaxy identification is a crucial step for modern supernova (SN)
surveys such as the Dark Energy Survey (DES) and the Large Synoptic Survey
Telescope (LSST), which will discover SNe by the thousands. Spectroscopic
resources are limited, so in the absence of realtime SN spectra these surveys
must rely on host galaxy spectra to obtain accurate redshifts for the Hubble
diagram and to improve photometric classification of SNe. In addition, SN
luminosities are known to correlate with hostgalaxy properties. Therefore,
reliable identification of host galaxies is essential for cosmology and SN
science. We simulate SN events and their locations within their host galaxies
to develop and test methods for matching SNe to their hosts. We use both real
and simulated galaxy catalog data from the Advanced Camera for Surveys General
Catalog and MICECATv2.0, respectively. We also incorporate "hostless" SNe
residing in undetected faint hosts into our analysis, with an assumed hostless
rate of 5%. Our fully automated algorithm is run on catalog data and matches
SNe to their hosts with 91% accuracy. We find that including a machine learning
component, run after the initial matching algorithm, improves the accuracy
(purity) of the matching to 97% with a 2% cost in efficiency (true positive
rate). Although the exact results are dependent on the details of the survey
and the galaxy catalogs used, the method of identifying host galaxies we
outline here can be applied to any transient survey.

We present cosmological constraints from the Dark Energy Survey (DES) using a
combined analysis of angular clustering of red galaxies and their
crosscorrelation with weak gravitational lensing of background galaxies. We
use a 139 square degree contiguous patch of DES data from the Science
Verification (SV) period of observations. Using large scale measurements, we
constrain the matter density of the Universe as Omega_m = 0.31 +/ 0.09 and the
clustering amplitude of the matter power spectrum as sigma_8 = 0.74 +/ 0.13
after marginalizing over seven nuisance parameters and three additional
cosmological parameters. This translates into S_8 = sigma_8(Omega_m/0.3)^{0.16}
= 0.74 +/ 0.12 for our fiducial lens redshift bin at 0.35 <z< 0.5, while S_8 =
0.78 +/ 0.09 using two bins over the range 0.2 <z< 0.5. We study the
robustness of the results under changes in the data vectors, modelling and
systematics treatment, including photometric redshift and shear calibration
uncertainties, and find consistency in the derived cosmological parameters. We
show that our results are consistent with previous cosmological analyses from
DES and other data sets and conclude with a joint analysis of DES angular
clustering and galaxygalaxy lensing with Planck CMB data, Baryon Accoustic
Oscillations and Supernova type Ia measurements.

Next generation galaxy surveys demand the development of massive ensembles of
galaxy mocks to model the observables and their covariances, what is
computationally prohibitive using $N$body simulations. COLA is a novel method
designed to make this feasible by following an approximate dynamics but with up
to 3 orders of magnitude speedups when compared to an exact $N$body. In this
paper we investigate the optimization of the code parameters in the compromise
between computational cost and recovered accuracy in observables such as
twopoint clustering and halo abundance. We benchmark those observables with a
stateoftheart $N$body run, the MICE Grand Challenge simulation (MICEGC).
We find that using 40 time steps linearly spaced since $z_i \sim 20$, and a
force mesh resolution three times finer than that of the number of particles,
yields a matter power spectrum within $1\%$ for $k \lesssim 1\,h {\rm
Mpc}^{1}$ and a halo mass function within $5\%$ of those in the $N$body. In
turn the halo bias is accurate within $2\%$ for $k \lesssim 0.7\,h {\rm
Mpc}^{1}$ whereas, in redshift space, the halo monopole and quadrupole are
within $4\%$ for $k \lesssim 0.4\,h {\rm Mpc}^{1}$. These results hold for a
broad range in redshift ($0 < z < 1$) and for all halo mass bins investigated
($M > 10^{12.5} \, h^{1} \, {\rm M_{\odot}}$). To bring accuracy in clustering
to one percent level we study various methods that recalibrate halo masses
and/or velocities. We thus propose an optimized choice of COLA code parameters
as a powerful tool to optimally exploit future galaxy surveys.

We present spectroscopic confirmation of two new lensed quasars via data
obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been
selected from the Dark Energy Survey (DES) and WISE based on their multiband
photometry and extended morphology in DES images. Images of DES J01155244 show
two blue point sources at either side of a red galaxy. Our longslit data
confirm that both point sources are images of the same quasar at $z_{s}=1.64.$
The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is
in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two
blue components are visible in the DES and SDSS images. The SDSS fiber spectrum
shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II
and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground
lens galaxy. The longslit Magellan spectra show that the blue components are
resolved images of the same quasar. The Einstein Radius is $0.68$"
corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other
candidates were observed and rejected, two being lowredshift pairs of
starburst galaxies, and one being a quasar behind a blue star. These first
confirmation results provide an important empirical validation of the
datamining and modelbased selection that is being applied to the entire DES
dataset.

We present a comparison of major methodologies of fast generating mock halo
or galaxy catalogues. The comparison is done for twopoint and the threepoint
clustering statistics. The reference catalogues are drawn from the BigMultiDark
Nbody simulation. Both friendoffriends (including distinct halos only) and
spherical overdensity (including distinct halos and subhalos) catalogs have
been used with the typical number density of a largevolume galaxy surveys. We
demonstrate that a proper biasing model is essential for reproducing the power
spectrum at quasilinear and even smaller scales. With respect to various
clustering statistics a methodology based on perturbation theory and a
realistic biasing model leads to very good agreement with Nbody simulations.
However, for the quadrupole of the correlation function or the power spectrum,
only the method based on semiNbody simulation could reach high accuracy (1%
level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full Nbody solutions
will remain indispensable to produce reference catalogues. Nevertheless, we
have demonstrated that the far more efficient approximate solvers can reach a
few percent accuracy in terms of clustering statistics at the scales
interesting for the largescale structure analysis after calibration with a few
reference Nbody calculations. This makes them useful for massive production
aimed at covariance studies, to scan large parameter spaces, and to estimate
uncertainties in data analysis techniques, such as baryon acoustic oscillation
reconstruction, redshift distortion measurements, etc.

Measurements of the linear growth factor $D$ at different redshifts $z$ are
key to distinguish among cosmological models. One can estimate the derivative
$dD(z)/d\ln(1+z)$ from redshift space measurements of the 3D anisotropic galaxy
twopoint correlation $\xi(z)$, but the degeneracy of its transverse (or
projected) component with galaxy bias $b$, i.e. $\xi_{\perp}(z) \propto\ D^2(z)
b^2(z)$, introduces large errors in the growth measurement. Here we present a
comparison between two methods which break this degeneracy by combining second
and thirdorder statistics. One uses the shape of the reduced threepoint
correlation and the other a combination of thirdorder one and twopoint
cumulants. These methods use the fact that, for Gaussian initial conditions and
scales larger than $20$ $h^{1}$Mpc, the reduced thirdorder matter
correlations are independent of redshift (and therefore of the growth factor)
while the thirdorder galaxy correlations depend on $b$. We use matter and halo
catalogs from the MICEGC simulation to test how well we can recover $b(z)$ and
therefore $D(z)$ with these methods in 3D real space. We also present a new
approach, which enables us to measure $D$ directly from the redshift evolution
of second and thirdorder galaxy correlations without the need of modelling
matter correlations. For haloes with masses lower than $10^{14}$
$h^{1}$M$_\odot$, we find $10%$ deviations between the different estimates of
$D$, which are comparable to current observational errors. At higher masses we
find larger differences that can probably be attributed to the breakdown of the
bias model and nonPoissonian shot noise.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to
ESA in May 2013 as a largeclass mission for investigating within the framework
of the ESA Cosmic Vision program a set of important scientific questions that
require high resolution, high sensitivity, fullsky observations of the sky
emission at wavelengths ranging from millimeterwave to the farinfrared.
PRISM's main objective is to explore the distant universe, probing cosmic
history from very early times until now as well as the structures, distribution
of matter, and velocity flows throughout our Hubble volume. PRISM will survey
the full sky in a large number of frequency bands in both intensity and
polarization and will measure the absolute spectrum of sky emission more than
three orders of magnitude better than COBE FIRAS. The aim of this Extended
White Paper is to provide a more detailed overview of the highlights of the new
science that will be made possible by PRISM

We study cell count moments up to fifth order of the distributions of haloes,
of halo substructures as a proxy for galaxies, and of mass in the context of
the halo model and compare theoretical predictions to the results of numerical
simulations. On scales larger than the size of the largest cluster, we present
a simple point cluster model in which results depend only on clustercluster
correlations and on the distribution of the number of objects within a cluster,
or cluster occupancy. The point cluster model leads to expressions for moments
of galaxy counts in which the volumeaveraged moments on large scales approach
those of the halo distribution and on smaller scales exhibit hierarchical
clustering with amplitudes $S_k$ determined by moments of the occupancy
distribution. In this limit, the halo model predictions are purely
combinatoric, and have no dependence on halo profile, concentration parameter,
or potential asphericity. The full halo model introduces only two additional
effects: on large scales, haloes of different mass have different clustering
strengths, introducing relative bias parameters; and on the smallest scales,
halo structure is resolved and details of the halo profile become important,
introducing shapedependent form factors. Because of differences between
discrete and continuous statistics, the hierarchical amplitudes for galaxies
and for mass behave differently on small scales even if galaxy number is
exactly proportional to mass, a difference that is not necessarily well
described in terms of bias.

We introduce a new set of large Nbody runs, the MICE simulations, that
provide a unique combination of very large cosmological volumes with good mass
resolution. They follow the gravitational evolution of ~ 8.5 billion particles
(2048^3) in volumes covering up to 450 (Gpc/h)^3. Our main goal is to
accurately model and calibrate basic cosmological probes that will be used by
upcoming astronomical surveys. Here we take advantage of the very large volumes
of MICE to make a robust sampling of the highmass tail of the halo mass
function (MF). We discuss and avoid possible systematic effects in our study,
and do a detailed analysis of different error estimators. We find that
available fits to the local abundance of halos (Warren et al. (2006)) match
well the abundance in MICE up to M ~ 10^{14}\Msun, but significantly deviate
for larger masses, underestimating the mass function by 10% (30%) at M = 3.16 x
10^{14}\Msun (10^{15}\Msun). Similarly, the widely used Sheth & Tormen (1999)
fit, if extrapolated to high redshift assuming universality, leads to an
underestimation of the cluster abundance by 30%, 20% and 15% at z=0, 0.5, 1 for
M ~ [7  2.5  0.8] x 10^{14}\Msun respectively ($\nu = \delta_c/\sigma ~ 3$).
We provide a recalibration of the halo MF valid over 5 orders of magnitude in
mass, 10^{10} < M/(\Msun) < 10^{15}, that accurately describes its redshift
evolution up to z=1. We explore the impact of this recalibration on the
determination of darkenergy, and conclude that using available fits may
systematically bias the estimate of w by as much as 50% for mediumdepth (z <=
1) surveys. MICE halo catalogues are publicly available at
http://www.ice.cat/mice

How did the universe evolve? The fine angular scale (l>1000) temperature and
polarization anisotropies in the CMB are a Rosetta stone for understanding the
evolution of the universe. Through detailed measurements one may address
everything from the physics of the birth of the universe to the history of star
formation and the process by which galaxies formed. One may in addition track
the evolution of the dark energy and discover the net neutrino mass.
We are at the dawn of a new era in which hundreds of square degrees of sky
can be mapped with arcminute resolution and sensitivities measured in
microKelvin. Acquiring these data requires the use of special purpose
telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and
the South Pole Telescope (SPT). These new telescopes are outfitted with a new
generation of custom mmwave kilopixel arrays. Additional instruments are in
the planning stages.

Crosscorrelation between the CMB and largescale structure is a powerful
probe of darkenergy and gravity on the largest physical scales. We introduce a
novel estimator, the CMBvelocity correlation, that has most of his power on
large scales and that, at low redshift, delivers up to factor of two higher
signaltonoise ratio than the recently detected CMBdark matter density
correlation expected from the Integrated SachsWolfe effect. We propose to use
a combination of peculiar velocities measured from supernovae type Ia and
kinetic SunyaevZeldovich cluster surveys to reveal this signal and forecast
darkenergy constraints that can be achieved with future surveys. We stress
that low redshift peculiar velocity measurements should be exploited with
complementary deeper largescale structure surveys for precision cosmology.

Constraining cosmological parameters from measurements of the Integrated
SachsWolfe effect requires developing robust and accurate methods for
computing statistical errors in the crosscorrelation between maps. This paper
presents a detailed comparison of such error estimation applied to the case of
crosscorrelation of Cosmic Microwave Background (CMB) and largescale
structure data. We compare theoretical models for error estimation with
montecarlo simulations where both the galaxy and the CMB maps vary around a
fiducial autocorrelation and crosscorrelation model which agrees well with
the current concordance LCDM cosmology. Our analysis compares estimators both
in harmonic and configuration (or real) space, quantifies the accuracy of the
error analysis and discuss the impact of partial sky survey area and the choice
of input fiducial model on darkenergy constraints. We show that purely
analytic approaches yield accurate errors even in surveys that cover only 10%
of the sky and that parameter constraints strongly depend on the fiducial model
employed. Alternatively, we discuss the advantages and limitations of error
estimators that can be directly applied to data. In particular, we show that
errors and covariances from the JackKnife method agree well with the
theoretical approaches and simulations. We also introduce a novel method in
real space that is computationally efficient and can be applied to real data
and realistic survey geometries. Finally, we present a number of new findings
and prescriptions that can be useful for analysis of real data and forecasts,
and present a critical summary of the analyses done to date.

The Dark Energy Survey (DES) will use a new imaging camera on the Blanco 4m
telescope at CTIO to image 5000 square degrees of sky in the South Galactic Cap
in four optical bands, and to carry out repeat imaging over a smaller area to
identify and measure lightcurves of Type Ia supernovae. The main imaging area
overlaps the planned SunyaevZel'dovich survey of the South Pole Telescope. The
idea behind DES is to use four distinct and largely independent methods to
probe the properties of dark energy: baryon oscillations of the power spectrum,
abundance and spatial distribution of clusters, weak gravitational lensing, and
Type Ia supernovae. This white paper outlines, in broad terms, some of the
theoretical issues associated with the first three of these probes (the issues
for supernovae are mostly different in character), and with the general task of
characterizing dark energy and distinguishing it from alternative explanations
for cosmic acceleration. A companion white paper discusses the kind of
numerical simulations and other theoretical tools that will be needed to
address the these issues and to create mock catalogs that allow endtoend
tests of analysis procedures. Although we have been thinking about these
problems in the specific context of DES, many of them are also relevant to
other planned dark energy studies.

We measure the angular power spectrum of the WMAP firstyear temperature
anisotropy maps. We use SpICE (Spatially Inhomogeneous Correlation Estimator)
to estimate Cl's for multipoles l=2900 from all possible crosscorrelation
channels. Except for the mapmaking stage, our measurements provide an
independent analysis of that by Hinshaw etal (2003). Despite the different
methods used, there is virtually no difference between the two measurements for
l < 700 ; the highest l's are still compatible within 1sigma errors. We use a
novel intrabin variance method to constrain Cl errors in a model independent
way. When applied to WMAP data, the intrabin variance estimator yields
diagonal errors 10% larger than those reported by the WMAP team for 100 < l <
450. This translates into a 2.4 sigma detection of systematics since no
difference is expected between the SpICE and the WMAP team estimator window
functions in this multipole range. With our measurement of the Cl's and errors,
we get chi^2/d.o.f. = 1.042 for a bestfit LCDM model, which has a 14%
probability, whereas the WMAP team obtained chi^2/d.o.f. = 1.066, which has a
5% probability. We assess the impact of our results on cosmological parameters
using Markov Chain Monte Carlo simulations. From WMAP data alone, assuming
spatially flat power law LCDM models, we obtain the reionization optical depth
tau = 0.145 +/ 0.067, spectral index n_s = 0.99 +/ 0.04, Hubble constant h =
0.67 +/ 0.05, baryon density Omega_b h^2 = 0.0218 +/ 0.0014, cold dark matter
density Omega_{cdm} h^2 = 0.122 +/ 0.018, and sigma_8 = 0.92 +/ 0.12,
consistent with a reionization redshift z_{re} = 16 +/ 5 (68% CL).

Models with late time cosmic acceleration, such as the Lambdadominated CDM
model, predict a freeze out for the growth of linear gravitational potential at
moderate redshift z<1, what can be observed as temperature anisotropies in the
CMB: the so called integrated SachsWolfe (ISW) effect. We present a direct
measurement of the ISW effect based on the angular crosscorrelation function,
w_{TG}, of CMB temperature anisotropies and darkmatter fluctuations traced by
galaxies. We crosscorrelate the firstyear WMAP data in combination with the
APM Galaxy survey. On the largest scales, theta = 410 deg, we detect a
nonvanishing crosscorrelation at 98.8 % significance level, with a 1sigma
error of w_{TG} = 0.35 +/ 0.14 microK, which favors large values of
Omega_Lambda \simeq 0.8 for flat FRW models. On smaller scales, theta < 1deg,
the correlations disappear. This is contrary to what would be expected from the
ISW effect, but the absence of correlations may be simply explained if the ISW
signal was being cancelled by anticorrelations arising from the thermal
SunyaevZeldovich (SZ) effect.

We present results for the crosscorrelation between the WMAP 1styear cosmic
microwave background (CMB) anisotropy data and optical galaxy surveys: the APM
and SDSS DR1 catalogs. Our measurement of a positive CMBgalaxy correlation on
large angles (\theta > 4 deg) yields significant detections of the Integrated
SachsWolfe (ISW) effect and provides a new estimate of darkenergy in the
universe, \Omega_\Lambda=0.690.86 (2 \sigma range). In addition, the
correlated signal on small angles (\theta<1 deg) reveals the imprint left by
hot intracluster gas in the CMB photons: the thermal SunyaevZeldovich (SZ)
effect

We present a crosscorrelation analysis of the WMAP cosmic microwave
background (CMB) temperature anisotropies and the SDSS galaxy density
fluctuations. We find significant detections of the angular CMBgalaxy
correlation for both the flux limited galaxy sample (z~0.3) and the high
redshift (z ~ 0.5) color selected sample. The signal is compatible with that
expected from the integrated SachsWolfe (ISW) effect at large angles (\theta >
3deg) and the SunyaevZeldovich (SZ) effect at small scales (\theta < 1 deg).
The detected correlation at lowz is in good agreement with a previous analysis
using the APM survey (z~0.15). The combined analysis of all 3 samples yields a
total significance better than 3 sigma for ISW and about 2.7 \sigma for SZ,
with a Compton parameter y~10^(6). For a given flat LCDM model, the ISW effect
depends both on the value of \Omega_\Lambda and the galaxy bias b. To break
this degeneracy, we estimate the bias using the ratio between the galaxy and
mass autocorrelation functions in each sample. With our bias estimation, all
samples consistently favor a best fit darkenergy dominated model:
\Omega_\Lambda ~ 0.8, with a 2 \sigma error \Omega_\Lambda=0.690.86.

We present a statistical analysis of the interstellar medium (ISM)
polarization from the largest compilation available of starlight data, which
comprises ~ 5500 stars. The measured correlation between the mean polarization
degree and extinction indicates that ISM dust grains are not fully aligned with
the uniform component of the largescale Galactic magnetic field. Moreover, we
estimate the ratio of the uniform to the random planeofthesky components of
the magnetic field to be B_u/B_r = 0.8. From the analysis of starlight
polarization degree and position angle we find that the magnetic field broadly
follows Galactic structures on largescales. On the other hand, the angular
power spectrum C_l of the polarization degree for Galactic plane data is found
to be consistent with a powerlaw, C_l ~ l^{1.5} (where l = 180 deg/\theta is
the multipole order), for angular scales \theta > 10 arcmin. We argue that this
data set can be used to estimate diffuse polarized emission at microwave
frequencies.