We present a combined radio/X-ray analysis of the poorly studied galaxy cluster Abell 2495 (z=0.07923) based on new EVLA and Chandra data. We also analyze and discuss Halpha emission and optical continuum data retrieved from the literature. We find an offset of 6 kpc between the cluster BCG (MCG+02-58-021) and the peak of the X-ray emission, suggesting that the cooling process is not taking place on the central galaxy nucleus. We propose that sloshing of the ICM could be responsible for this separation. Furthermore, we detect a second, 4 kpc offset between the peak of the Halpha emission and that of the X-ray emission. Optical images highlight the presence of a dust filament extending up to 6 kpc in the cluster BCG, and allow us to estimate a dust mass within the central 7 kpc of 1.7e+5 Msun. Exploiting the dust to gas ratio and the L_Halpha-M_mol relation, we argue that a significant amount (up to 10^9 Msun) of molecular gas should be present in the BCG of this cluster. We also investigate the presence of ICM depressions, finding two putative systems of cavities; the inner pair is characterized by t age = 18 Myr and P cav = 1.2e+43 erg/s, the outer one by t age = 53 Myr and P cav = 5.6e+42 erg/s. Their age difference appears to be consistent with the free-fall time of the central cooling gas and with the offset timescale estimated with the Halpha kinematic data, suggesting that sloshing is likely playing a key role in this environment. Furthermore, the cavities' power analysis shows that the AGN energy injection is able to sustain the feedback cycle, despite cooling being offset from the BCG nucleus.

We present new ALMA CO(2--1) observations of two well studied group-centered elliptical galaxies: NGC~4636 and NGC~5846. In addition, we include a revised analysis of Cycle 0 ALMA observations of the central galaxy in the NGC~5044 group that has been previously published. We find evidence that molecular gas, in the form of off-center orbiting clouds, is a common presence in bright group-centered galaxies (BGG). CO line widths are $\gtrsim 10$ times broader than Galactic molecular clouds, and using the reference Milky Way $X_{CO}$, the total molecular mass ranges from as low as $2.6\times 10^5 M_\odot$ in NGC~4636 to $6.1\times 10^7 M_\odot$ in NGC~5044. With these parameters the virial parameters of the molecular structures is $\gg 1$. Complementary observations of NGC~5846 and NGC~4636 using the ALMA Compact Array (ACA) do not exhibit any detection of a CO diffuse component at the sensitivity level achieved by current exposures. The origin of the detected molecular features is still uncertain, but these ALMA observations suggest that they are the end product of the hot gas cooling process and not the result of merger events. Some of the molecular clouds are associated with dust features as revealed by HST dust extinction maps suggesting that these clouds formed from dust-enhanced cooling. The global nonlinear condensation may be triggered via the chaotic turbulent field or buoyant uplift. The large virial parameter of the molecular structures and correlation with the warm ($10^3 - 10^5 K$)/hot ($\ge10^6$) phase velocity dispersion provide evidence that they are unbound giant molecular associations drifting in the turbulent field, consistently with numerical predictions of the chaotic cold accretion process. Alternatively, the observed large CO line widths may be generated by molecular gas flowing out from cloud surfaces due to heating by the local hot gas atmosphere.

We use high-resolution Herschel-PACS data of 2 nearby elliptical galaxies, IC1459 & NGC2768 to characterize their dust and stellar content. IC1459 & NGC2768 have an unusually large amount of dust for elliptical galaxies (1-3 x 10^5 Msun), this dust is also not distributed along the stellar content. Using data from GALEX (ultraviolet) to PACS (far-infrared), we analyze the spectral energy distribution (SED) of these galaxies with CIGALEMC as a function of the projected position, binning images in 7.2" pixels. From this analysis, we derive maps of SED parameters, such as the metallicity, the stellar mass, the fraction of young star and the dust mass. The larger amount of dust in FIR maps seems related in our model to a larger fraction of young stars which can reach up to 4% in the dustier area. The young stellar population is fitted as a recent (~ 0.5 Gyr) short burst of star formation for both galaxies. The metallicities, which are fairly large at the center of both galaxies, decrease with the radial distance with fairly steep gradient for elliptical galaxies.

We explore the properties of early-type galaxies (ETGs), including ellipticals (E) and lenticulars (S0), in rich environments such as clusters of galaxies (Virgo and Coma). The L_24/L_K distribution of ETGs in both Virgo and Coma clusters shows that some S0s have a much larger L_24/L_K ratio (0.5 to ~2 dex) than the bulk of the ETG population. This could be interpreted as an enhanced star formation rate in these lenticulars. We compare the optical colors of galaxies in these two clusters and investigate the nature of these sources with a large L24/L_K ratio by looking at their spatial distribution within the cluster, by analyzing their optical spectra and by looking at their optical colors compared to late-types. We obtain 10 Coma and 3 Virgo early-type sources with larger L24/L_K ratios than the bulk of their population. We call these sources Mid-Infrared Enhanced Galaxies (MIEGs). In Coma, they are mostly located in the South-West part of the cluster where a substructure is falling onto the main cluster. MIEGs present lower g-r color than the rest of the ETG sample, because of a blue continuum. We interpret the excess L24/L_K ratio as evidence for an enhanced star-formation induced as a consequence of their infall into the main cluster.

We compute the properties of a sample of 221 local early-type galaxies with a spectral energy distribution (SED) modeling software, CIGALEMC. Concentrating on the star forming activity and dust contents, we derive parameters such as the specific star formation rate, the dust luminosity, dust mass and temperature. 52 % of our sample is composed of elliptical (E) galaxies and 48 % of lenticular (S0) galaxies. We find a larger proportion of S0 galaxies among galaxies with a large specific star formation rate (sSFR) and large specific dust emission. The stronger activity of S0 galaxies is confirmed by larger dust masses. We investigate the relative proportion of active galactic nuclei (AGN) and star-forming (SF) galaxies in our sample using spectroscopic SDSS data and near-infrared selection techniques, and find a larger proportion of AGN dominated galaxy in the S0 sample than the E one. This could corroborate a scenario where blue galaxies evolve into red ellipticals by passing through a S0 AGN active period while quenching its star formation. Finally, we find a good agreement comparing our estimates with color indicators.

The old, red stars which constitute the bulges of galaxies, and the massive black holes at their centres, are the relics of a period in cosmic history when galaxies formed stars at remarkable rates and active galactic nuclei (AGN) shone brightly from accretion onto black holes. It is widely suspected, but unproven, that the tight correlation in mass of the black hole and stellar components results from the AGN quenching the surrounding star formation as it approaches its peak luminosity. X-rays trace emission from AGN unambiguously, while powerful star-forming galaxies are usually dust-obscured and are brightest at infrared to submillimetre wavelengths. Here we report observations in the submillimetre and X-ray which show that rapid star formation was common in the host galaxies of AGN when the Universe was 2-6 Gyrs old, but that the most vigorous star formation is not observed around black holes above an X-ray luminosity of 10^44 erg/s. This suppression of star formation in the host galaxies of powerful AGN is a key prediction of models in which the AGN drives a powerful outflow, expelling the interstellar medium of its host galaxy and transforming the galaxy's properties in a brief period of cosmic time.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

We present measurements of the auto- and cross-frequency power spectra of the cosmic infrared background (CIB) at 250, 350, and 500um (1200, 860, and 600 GHz) from observations totaling ~ 70 deg^2 made with the SPIRE instrument aboard the Herschel Space Observatory. We measure a fractional anisotropy dI / I = 14 +- 4%, detecting signatures arising from the clustering of dusty star-forming galaxies in both the linear (2-halo) and non-linear (1-halo) regimes; and that the transition from the 2- to 1-halo terms, below which power originates predominantly from multiple galaxies within dark matter halos, occurs at k_theta ~ 0.1 - 0.12 arcmin^-1 (l ~ 2160 - 2380), from 250 to 500um. New to this paper is clear evidence of a dependence of the Poisson and 1-halo power on the flux-cut level of masked sources --- suggesting that some fraction of the more luminous sources occupy more massive halos as satellites, or are possibly close pairs. We measure the cross-correlation power spectra between bands, finding that bands which are farthest apart are the least correlated, as well as hints of a reduction in the correlation between bands when resolved sources are more aggressively masked. In the second part of the paper we attempt to interpret the measurements in the framework of the halo model. With the aim of fitting simultaneously with one model the power spectra, number counts, and absolute CIB level in all bands, we find that this is achievable by invoking a luminosity-mass relationship, such that the luminosity-to-mass ratio peaks at a particular halo mass scale and declines towards lower and higher mass halos. Our best-fit model finds that the halo mass which is most efficient at hosting star formation in the redshift range of peak star-forming activity, z ~ 1-3, is log(M_peak/M_sun) ~ 12.1 +- 0.5, and that the minimum halo mass to host infrared galaxies is log(M_min/M_sun) ~ 10.1 +- 0.6.

ABRIGED Herschel/SPIRE has provided confusion limited maps of deep fields at 250, 350, and 500um, as part of the HerMES survey. Due to confusion, only a small fraction of the Cosmic Infrared Background can be resolved into individually-detected sources. Our goal is to produce deep galaxy number counts and redshift distributions below the confusion limit, which we then use to place strong constraints on the origins of the cosmic infrared background and on models of galaxy evolution. We individually extracted the bright SPIRE with a method using the positions, the flux densities, and the redshifts of the 24um sources as a prior, and derived the number counts and redshift distributions of the bright SPIRE sources. For fainter SPIRE sources, we reconstructed the number counts and the redshift distribution below the confusion limit using the deep 24um catalogs associated with photometric redshift and information provided by the stacking of these sources into the deep SPIRE maps. Finally, by integrating all these counts, we studied the contribution of the galaxies to the CIB as a function of their flux density and redshift. Through stacking, we managed to reconstruct the source counts per redshift slice down to ~2 mJy in the three SPIRE bands, which lies about a factor 10 below the 5sigma confusion limit. None of the pre-existing population models are able to reproduce our results at better than 3sigma. Finally, we extrapolate our counts to zero flux density in order to derive an estimate of the total contribution of galaxies to the CIB, finding 10.1, 6.5, and 2.8 nW/m2/sr at 250, 350, and 500um, respectively. These values agree well with FIRAS absolute measurements, suggesting our number counts and their extrapolation are sufficient to explain the CIB. Finally, combining our results with other works, we estimate the energy budget contained in the CIB between 8 and 1000um: 26 nW/m2/sr.

The Hubble morphological sequence from early to late galaxies corresponds to an increasing rate of specific star formation. The Hubble sequence also follows a banana-shaped correlation between 24 and 70 micron luminosities, both normalized with the K-band luminosity. We show that this correlation is significantly tightened if galaxies with central AGN emission are removed, but the cosmic scatter of elliptical galaxies in both 24 and 70 micron luminosities remains significant along the correlation. We find that the 24 micron variation among ellipticals correlates with stellar metallicity, reflecting emission from hot dust in winds from asymptotic giant branch stars of varying metallicity. Infrared surface brightness variations in elliptical galaxies indicate that the K - 24 color profile is U-shaped for reasons that are unclear. In some elliptical galaxies cold interstellar dust emitting at 70 and 160 microns may arise from recent gas-rich mergers. However, we argue that most of the large range of 70 micron luminosity in elliptical galaxies is due to dust transported from galactic cores by feedback events in (currently IR-quiet) active galactic nuclei. Cooler dusty gas naturally accumulates in the cores of elliptical galaxies due to dust-cooled local stellar mass loss and may accrete onto the central black hole, releasing energy. AGN-heated gas can transport dust in cores 5-10 kpc out into the hot gas atmospheres where it radiates extended 70 micron emission but is eventually destroyed by sputtering. This, and some modest star formation, defines a cycle of dust creation and destruction. Elliptical galaxies evidently undergo large transient excursions in the banana plot in times comparable to the sputtering time or AGN duty cycle, 10 Myrs. Normally regarded as passive, elliptical galaxies are the most active galaxies in the IR color-color correlation.

We have constructed an extended halo model (EHM) which relates the total stellar mass and star-formation rate (SFR) to halo mass (M_h). An empirical relation between the distribution functions of total stellar mass of galaxies and host halo mass, tuned to match the spatial density of galaxies over 0<z<2 and the clustering properties at z~0, is extended to include two different scenarios describing the variation of SFR on M_h. We also present new measurements of the redshift evolution of the average SFR for star-forming galaxies of different stellar mass up to z=2, using data from the Herschel Multi-tiered Extragalactic Survey (HerMES) for infrared-bright galaxies. Combining the EHM with the halo accretion histories from numerical simulations, we trace the stellar mass growth and star-formation history in halos spanning a range of masses. We find that: (1) The intensity of the star-forming activity in halos in the probed mass range has steadily decreased from z~2 to 0; (2) At a given epoch, halos in the mass range between a few times 10^{11} M_Sun and a few times 10^{12} M_Sun are the most efficient at hosting star formation; (3) The peak of SFR density shifts to lower mass halos over time; (4) Galaxies that are forming stars most actively at z~2 evolve into quiescent galaxies in today's group environments, strongly supporting previous claims that the most powerful starbursts at z~2 are progenitors of today's elliptical galaxies.

We present a measurement of the angular power spectrum of the cosmic far-infrared background (CFIRB) anisotropies in one of the extragalactic fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) at 250, 350 and 500 \mu m bands. Consistent with recent measurements of the CFIRB power spectrum in Herschel-SPIRE maps, we confirm the existence of a clear one-halo term of galaxy clustering on arcminute angular scales with large-scale two-halo term of clustering at 30 arcminutes to angular scales of a few degrees. The power spectrum at the largest angular scales, especially at 250 \mu m, is contaminated by the Galactic cirrus. The angular power spectrum is modeled using a conditional luminosity function approach to describe the spatial distribution of unresolved galaxies that make up the bulk of the CFIRB. Integrating over the dusty galaxy population responsible for the background anisotropies, we find that the cosmic abundance of dust, relative to the critical density, to be between \Omega_dust=10^{-6} and 8 x 10^{-6} in the redshift range z ~ 0-3. This dust abundance is consistent with estimates of the dust content in the Universe using quasar reddening and magnification measurements in the SDSS.

We present a list of 13 candidate gravitationally lensed submillimeter galaxies (SMGs) from 95 square degrees of the Herschel Multi-tiered Extragalactic Survey, a surface density of 0.14\pm0.04deg^{-2}. The selected sources have 500um flux densities (S_500) greater than 100mJy. Gravitational lensing is confirmed by follow-up observations in 9 of the 13 systems (70%), and the lensing status of the four remaining sources is undetermined. We also present a supplementary sample of 29 (0.31\pm0.06deg^{-2}) gravitationally lensed SMG candidates with S_500=80--100mJy, which are expected to contain a higher fraction of interlopers than the primary candidates. The number counts of the candidate lensed galaxies are consistent with a simple statistical model of the lensing rate, which uses a foreground matter distribution, the intrinsic SMG number counts, and an assumed SMG redshift distribution. The model predicts that 32--74% of our S_500>100mJy candidates are strongly gravitationally lensed (mu>2), with the brightest sources being the most robust; this is consistent with the observational data. Our statistical model also predicts that, on average, lensed galaxies with S_500=100mJy are magnified by factors of ~9, with apparently brighter galaxies having progressively higher average magnification, due to the shape of the intrinsic number counts. 65% of the sources are expected to have intrinsic 500micron flux densities less than 30mJy. Thus, samples of strongly gravitationally lensed SMGs, such as those presented here, probe below the nominal Herschel detection limit at 500 micron. They are good targets for the detailed study of the physical conditions in distant dusty, star-forming galaxies, due to the lensing magnification, which can lead to spatial resolutions of ~0.01" in the source plane.

We report new, sensitive observations of two z ~ 3-3.5 FIR-luminous radio galaxies, 6C1909+72 and B3J2330+3927, in 12CO J=1-0 with the Karl Jansky VLA and at 100-500um using Herschel, alongside new and archival 12CO J=4-3 observations from IRAM PdBI. We introduce a new colour-colour diagnostic plot to constrain the redshifts of several distant, dusty galaxies in our target fields. A bright SMG near 6C1909+72 likely shares the same node or filament as the signpost AGN, but it is not detected in CO despite ~20,000 km/s of velocity coverage. Also in the 6C1909+72 field, a large, red dust feature spanning ~500 kpc is aligned with the radio jet. We suggest several processes by which metal-rich material may have been transported, favouring a collimated outflow reminiscent of the jet-oriented metal enrichment seen in local cluster environments. Our interferometric imaging reveals a gas-rich companion to B3J2330+3927; indeed, all bar one of the eight z >~ 2 radio galaxies (or companions) detected in CO provide some evidence that starburst activity in radio-loud AGN at high redshift is driven by the interaction of two or more gas-rich systems in which a significant mass of stars has already formed, rather than via steady accretion of cold gas from the cosmic web. We find that the CO Tb ratios in radio-loud AGN host galaxies are significantly higher than those seen in similarly intense starbursts where AGN activity is less pronounced. Our most extreme example, where L'(CO4-3)/L'(CO1-0) > 2.7, provides evidence that significant energy is being deposited rapidly into the molecular gas via X-rays and/or mechanical (`quasar-mode') feedback from the AGN, leading to a high degree of turbulence globally and a low optical depth in 12CO - feedback that may lead to the cessation of star formation on a timescale commensurate with that of the jet activity, <~10 Myr.

We describe the infrared properties of sources detected over ~36 deg^2 of sky in the GAMA 15-hr equatorial field, using data from both the Herschel Astrophysical Terahertz Large-Area Survey (H-ATLAS) and Wide-field Infrared Survey (WISE). With 5-sigma point-source depths of 34 and 0.048 mJy at 250 micron and 3.4 micron, respectively, we are able to identify 50.6% of the H-ATLAS sources in the WISE survey, corresponding to a surface density of ~630 deg^{-2}. Approximately two-thirds of these sources have measured spectroscopic or optical/near-IR photometric redshifts of z<1. For sources with spectroscopic redshifts at z<0.3, we find a linear correlation between the infrared luminosity at 3.4 micron and that at 250 micron, with +-50% scatter over ~1.5 orders of magnitude in luminosity, ~10^9 - 10^{10.5} L_sun. By contrast, the matched sources without previously measured redshifts (r>~20.5) have 250-350 micron flux density ratios that suggest either high-redshift galaxies (z>~1.5) or optically faint low-redshift galaxies with unusually low temperatures (T<~20). Their small 3.4-250 micron flux ratios favor a high-redshift galaxy population, as only the most actively star-forming galaxies at low redshift (e.g., Arp 220) exhibit comparable flux density ratios. Furthermore, we find a relatively large AGN fraction (~30%) in a 12 micron flux-limited subsample of H-ATLAS sources, also consistent with there being a significant population of high-redshift sources in the no-redshift sample.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

While the arcminute-scale Cosmic Microwave Background (CMB) anisotropies are due to secondary effects, point sources dominate the total anisotropy power spectrum. At high frequencies the point sources are primarily in the form of dusty, star-forming galaxies. Both Herschel and Planck have recently measured the anisotropy power spectrum of cosmic infrared background (CIB) generated by dusty, star-forming galaxies from degree to sub-arcminute angular scales, including the non-linear clustering of these galaxies at multipoles of 3000 to 6000 relevant to CMB secondary anisotropy studies. We scale the CIB angular power spectra to CMB frequencies and interpret the combined WMAP-7 year and arcminute-scale Atacama Cosmology Telescope (ACT) and South Pole Telescope (SPT) CMB power spectra measurements to constrain the Sunyaev-Zel'dovich (SZ) effects. Allowing the CIB clustering amplitude to vary, we constrain the amplitudes of thermal and kinetic SZ power spectra at 150 GHz.

We investigate the potential of submm-mm and submm-mm-radio photometric redshifts using a sample of mm-selected sources as seen at 250, 350 and 500 {\mu}m by the SPIRE instrument on Herschel. From a sample of 63 previously identified mm-sources with reliable radio identifications in the GOODS-N and Lockman Hole North fields 46 (73 per cent) are found to have detections in at least one SPIRE band. We explore the observed submm/mm colour evolution with redshift, finding that the colours of mm-sources are adequately described by a modified blackbody with constant optical depth {\tau} = ({\nu}/{\nu}0)^{\beta} where {\beta} = +1.8 and {\nu}0 = c/100 {\mu}m. We find a tight correlation between dust temperature and IR luminosity. Using a single model of the dust temperature and IR luminosity relation we derive photometric redshift estimates for the 46 SPIRE detected mm-sources. Testing against the 22 sources with known spectroscopic, or good quality optical/near-IR photometric, redshifts we find submm/mm photometric redshifts offer a redshift accuracy of |z|/(1+z) = 0.16 (< |z| >= 0.51). Including constraints from the radio-far IR correlation the accuracy is improved to |z|/(1 + z) = 0.15 (< |z| >= 0.45). We estimate the redshift distribution of mm-selected sources finding a significant excess at z > 3 when compared to ~ 850 {\mu}m selected samples.

We introduce a fast Markov Chain Monte Carlo (MCMC) exploration of the astrophysical parameter space using a modified version of the publicly available code CIGALE (Code Investigating GALaxy emission). The original CIGALE builds a grid of theoretical Spectral Energy Distribution (SED) models and fits to photometric fluxes from Ultraviolet (UV) to Infrared (IR) to put contraints on parameters related to both formation and evolution of galaxies. Such a grid-based method can lead to a long and challenging parameter extraction since the computation time increases exponentially with the number of parameters considered and results can be dependent on the density of sampling points, which must be chosen in advance for each parameter. Markov Chain Monte Carlo methods, on the other hand, scale approximately linearly with the number of parameters, allowing a faster and more accurate exploration of the parameter space by using a smaller number of efficiently chosen samples. We test our MCMC version of the code CIGALE (called CIGALEMC) with simulated data. After checking the ability of the code to retrieve the input parameters used to build the mock sample, we fit theoretical SEDs to real data from the well known and studied SINGS sample. We discuss constraints on the parameters and show the advantages of our MCMC sampling method in terms of accuracy of the results and optimization of CPU time.

The Sunyaev-Zel'dovich (SZ) effect has a distinct spectral signature that allows its separation from fluctuations in the cosmic microwave background (CMB) and foregrounds. Using CMB anisotropies measured in Wilkinson Microwave Anisotropy Probe's five-year maps, we constrain the SZ fluctuations at large, degree angular scales corresponding to multipoles in the range from 10 to 400. We provide upper bounds on SZ fluctuations at multipoles greater than 50, and find evidence for a hemispherically asymmetric signal at ten degrees angular scales. The amplitude of the detected signal cannot be easily explained with the allowed number density and temperature of electrons in the Galactic halo. We have failed to explain the excess signal as a residual from known Galactic foregrounds or instrumental uncertainties such as 1/f-noise.

We present the results of a gravitational lensing analysis of the bright $\zs=2.957$ sub-millimeter galaxy (SMG), HERMES J105751.1+573027 found in {\it Herschel}/SPIRE Science Demonstration Phase data from the Herschel Multi-tiered Extragalactic Survey (HerMES) project. The high resolution imaging available in optical and Near-IR channels, along with CO emission obtained with the Plateau de Bure Interferometer, allow us to precisely estimate the intrinsic source extension and hence estimate the total lensing magnification to be $\mu=10.9\pm 0.7$. We measure the half-light radius $R_{\rm eff}$ of the source in the rest-frame Near-UV and $V$ bands that characterize the unobscured light coming from stars and find $R_{\rm eff,*}= [2.0 \pm 0.1]$ kpc, in good agreement with recent studies on the Submillimeter Galaxy population. This lens model is also used to estimate the size of the gas distribution ($R_{\rm eff,gas}= [1.1\pm0.5]$) kpc by mapping back in the source plane the CO (J=5-4) transition line emission. The lens modeling yields a relatively large Einstein radius $R_{\rm Ein}= 4\farcs10 \pm 0\farcs02$, corresponding to a deflector velocity dispersion of [$483\pm 16] \,\kms$. This shows that HERMES J105751.1+573027 is lensed by a {\it galaxy group-size} dark matter halo at redshift $\zl\sim 0.6$. The projected dark matter contribution largely dominates the mass budget within the Einstein radius with $f_{\rm dm}(<R_{\rm Ein})\sim 80%$. This fraction reduces to $f_{\rm dm}(<R_{\rm eff,G1}\simeq 4.5\kpc)\sim 47%$ within the effective radius of the main deflecting galaxy of stellar mass $M_{\rm *,G1}=[8.5\pm 1.6] \times 10^{11}\msun$. At this smaller scale the dark matter fraction is consistent with results already found for massive lensing ellipticals at $z\sim0.2$ from the SLACS survey.

We report the discovery of a bright ($f(250\mum) > 400$ mJy), multiply-lensed submillimeter galaxy \obj\ in {\it Herschel}/SPIRE Science Demonstration Phase data from the HerMES project. Interferometric 880\mum\ Submillimeter Array observations resolve at least four images with a large separation of $\sim 9\arcsec$. A high-resolution adaptive optics $K_p$ image with Keck/NIRC2 clearly shows strong lensing arcs. Follow-up spectroscopy gives a redshift of $z=2.9575$, and the lensing model gives a total magnification of $\mu \sim 11 \pm 1$. The large image separation allows us to study the multi-wavelength spectral energy distribution (SED) of the lensed source unobscured by the central lensing mass. The far-IR/millimeter-wave SED is well described by a modified blackbody fit with an unusually warm dust temperature, $88 \pm 3$ K. We derive a lensing-corrected total IR luminosity of $(1.43 \pm 0.09) \times 10^{13}\, \mathrm{L}_{\odot}$, implying a star formation rate of $\sim 2500\, \mathrm{M}_{\odot}\, \mathrm{yr}^{-1}$. However, models primarily developed from brighter galaxies selected at longer wavelengths are a poor fit to the full optical-to-millimeter SED. A number of other strongly lensed systems have already been discovered in early {\it Herschel} data, and many more are expected as additional data are collected.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

The extragalactic background light at far-infrared wavelengths originates from optically-faint, dusty, star-forming galaxies in the universe with star-formation rates at the level of a few hundred solar masses per year. Due to the relatively poor spatial resolution of far-infrared telescopes, the faint sub-millimetre galaxies are challenging to study individually. Instead, their average properties can be studied using statistics such as the angular power spectrum of the background intensity variations. A previous attempt at measuring this power spectrum resulted in the suggestion that the clustering amplitude is below the level computed with a simple ansatz based on a halo model. Here we report a clear detection of the excess clustering over the linear prediction at arcminute angular scales in the power spectrum of brightness fluctuations at 250, 350, and 500 microns. From this excess, we find that sub-millimetre galaxies are located in dark matter halos with a minimum mass of log[M_min/M_sun ]= 11.5^+0.7_-0.2 at 350 microns. This minimum dark matter halo mass corresponds to the most efficient mass scale for star formation in the universe, and is lower than that predicted by semi-analytical models for galaxy formation.

The temperature fluctuations and polarization of the Cosmic Microwave Background (CMB) are now a well-known probe of the Universe at an infant age of 400,000 years. During the transit to us from the surface of last scattering, the CMB photons are expected to undergo modifications induced by the intervening large-scale structure. Among the expected secondary effects is the weak gravitational lensing of the CMB by the foreground dark matter distribution. We derive a quadratic estimator that uses the non-Gaussianities generated by the lensing effect at the four-point function level to extract the power spectrum of lensing potential fluctuations integrated out to z ~ 1100 with peak contributions from potential fluctuations at z of 2 to 3. Using WMAP 7-year temperature maps, we report the first direct constraints of this lensing potential power spectrum and find that it has an amplitude of A_L = 0.96 +/- 0.60, 1.06 +/- 0.69 and 0.97 +/- 0.47 using the W, V and W+V bands, respectively.

We study the Spectral Energy Distribution (SED) and the power spectrum of Galactic cirrus emission observed in the 14 deg^2 Science Demonstration Phase field of the Herschel-ATLAS using Herschel and IRAS data from 100 to 500 um. We compare the SPIRE 250, 350 and 500um maps with IRAS 100um emission, binned in 6' pixels. We assume a modified black-body SED with dust emissivity parameter beta (F ~ lambda^(-beta)) and a single dust temperature T_d, and find that the dust temperature and emissivity index varies over the science demonstration field as 10< T_rm < 25 K and 1 < beta< 4. The latter values are somewhat higher than the range of beta often quoted in the literature (1< beta< 2). We estimate the mean values of these parameters to be T_d=19.0 +/- 2.4 K and beta = 1.4 +/- 0.4. In regions of bright cirrus emission, we find that the dust has similar temperatures with T_d = 18.0 +/- 2.5 K, and similar values of beta, ranging from 1.4 +- 0.5 to 1.9+/- 0.5. We show that T_d and beta associated with diffuse cirrus emission are anti-correlated and can be described by the relationship: beta(T_d) = NT_d^alpha with [N=116+/-38, alpha=-1.4+/1 0.1]. The strong correlation found in this analysis is not just limited to high density clumps of cirrus emission as seen in previous studies, but is also seen in diffuse cirrus in low density regions. To provide an independent measure of $T_{\rm d}$ and $\beta$, we obtain the angular power spectrum of the cirrus emission in the {\it IRAS} and SPIRE maps, which is consistent with a power spectrum of the form P(k)=P_0(k/k_0)^gamma where gamma = \^a H R2.6+/-m 0.2 for scales of 50-200' in the SPIRE maps. The cirrus rms fluctuation amplitude at angular scales of 100' is consistent with a modified blackbody SED with T_d = 20.1+/- 0.9 K and beta = 1.3+/- 0.2, in agreement with the values obtained above.

We have analysed the rest-frame far infrared (FIR) properties of a sample of massive (Mstar > 10^11Msun) galaxies at 2<z<3 in the GOODS (Great Observatories Origins Deep Survey) North field using the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. To conduct this analysis we take advantage of the data from the HerMES key program. The sample comprises 45 massive galaxies with structural parameters characterised with HST NICMOS-3. We study detections at submm Herschel bands, together with Spitzer 24{\mu}m data, as a function of the morphological type, mass and size. We find that 26/45 sources are detected at MIPS-24{\mu}m and 15/45 (all MIPS-24{\mu}m detections) are detected at SPIRE-250{\mu}m, with disk-like galaxies more easily detected. We derive star formation rates (SFR) and specific star formation rates (sSFR) by fitting the spectral energy distribution (SED) of our sources, taking into account non-detections for SPIRE and systematic effects for MIPS derived quantities. We find that the mean SFR for the spheroidal galaxies (50-100 Msun*yr^-1) is substantially (a factor ~ 3) lower than the mean value presented by disk-like galaxies (250-300 Msun*yr^-1).

We describe the production and verification of sky maps of the five SPIRE fields observed as part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demonstration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that preserve extended diffuse emission on the sky while exploiting the repeated observations of the same region of the sky with many detectors in multiple scan directions to minimize residual instrument noise. We specify here the SHIM algorithm and outline the various tests that were performed to determine and characterize the quality of the maps and verify that the astrometry, point source flux and power on all relevant angular scales meets the needs of the HerMES science goals. These include multiple jackknife tests, determination of the map transfer function and detailed examination of the power spectra of both sky and jackknife maps. The map transfer function is approximately unity on scales from one arcminute to one degree. Final maps (v1.0), including multiple jackknives, as well as the SHIM pipeline, have been used by the HerMES team for the production of SDP papers.

Dusty, star forming galaxies contribute to a bright, currently unresolved cosmic far-infrared background. Deep Herschel-SPIRE images designed to detect and characterize the galaxies that comprise this background are highly confused, such that the bulk lies below the classical confusion limit. We analyze three fields from the HerMES programme in all three SPIRE bands (250, 350, and 500 microns); parameterized galaxy number count models are derived to a depth of ~2 mJy/beam, approximately 4 times the depth of previous analyses at these wavelengths, using a P(D) (probability of deflection) approach for comparison to theoretical number count models. Our fits account for 64, 60, and 43 per cent of the far-infrared background in the three bands. The number counts are consistent with those based on individually detected SPIRE sources, but generally inconsistent with most galaxy number counts models, which generically overpredict the number of bright galaxies and are not as steep as the P(D)-derived number counts. Clear evidence is found for a break in the slope of the differential number counts at low flux densities. Systematic effects in the P(D) analysis are explored. We find that the effects of clustering have a small impact on the data, and the largest identified systematic error arises from uncertainties in the SPIRE beam.

We present first results of a study of the submillimetre (rest frame far-infrared) properties of z~3 Lyman Break Galaxies (LBGs) and their lower-redshift counterparts BX/BM galaxies, based on Herschel-SPIRE observations of the Northern field of the Great Observatories Origins Deep Survey (GOODS-N). We use stacking analysis to determine the properties of LBGs well below the current limit of the survey. Although LBGs are not detected individually, stacking the infrared luminous LBGs (those detected with Spitzer at 24 microns yields a statistically significant submm detection with mean flux <S_{250}>= 5.9+/-1.4 mJy confirming the power of SPIRE in detecting UV-selected high-redshift galaxies at submillimetre wavelengths. In comparison, the Spitzer 24 microns detected BX/BM galaxies appear fainter with a stacked value of <S_{250}> = 2.7 +/-0.8 mJy. By fitting the Spectral Energy Distributions (SEDs) we derive median infrared luminosities, L_{IR}, of 2.8x10^{12} Lsun and 1.5x10^{11} Lsun for z~3 LBGs and BX/BMs, respectively. We find that $L_{IR} estimates derived from present measurements are in good agreement with those based on UV data for z~2 BX/BM galaxies, unlike the case for z~3 infrared luminous LBGs where the UV underestimates the true $L_{IR}. Although sample selection effects may influence this result we suggest that differences in physical properties (such as morphologies, dust distribution and extent of star-forming regions) between z ~3 LBGs and z~2 BX/BMs may also play a significant role.

The Spectral and Photometric Imaging Receiver (SPIRE) on Herschel has been carrying out deep extragalactic surveys, one of whose aims is to establish spectral energy distributions (SED)s of individual galaxies spanning the infrared/submillimeter (IR/SMM) wavelength region. We report observations of the (IR/SMM) emission from the Lockman North field (LN) and Great Observatories Origins Deep Survey field North (GOODS-N). Because galaxy images in the wavelength range covered by Herschel generally represent a blend with contributions from neighboring galaxies, we present sets of galaxies in each field especially free of blending at 250, 350, and 500 microns. We identify the cumulative emission of these galaxies and the fraction of the far infrared cosmic background radiation they contribute. Our surveys reveal a number of highly luminous galaxies at redshift z ~< 3 and a novel relationship between infrared and visible emission that shows a dependence on luminosity and redshift.

We present the cross-identification and source photometry techniques used to process Herschel SPIRE imaging taken as part of the Herschel Multi-Tiered Extragalactic Survey (HerMES). Cross-identifications are performed in map-space so as to minimise source blending effects. We make use of a combination of linear inversion and model selection techniques to produce reliable cross-identification catalogues based on Spitzer MIPS 24 micron source positions. Testing on simulations and real Herschel observations show that this approach gives robust results for even the faintest sources S250~10 mJy. We apply our new technique to HerMES SPIRE observations taken as part of the science demostration phase of Herschel. For our real SPIRE observations we show that, for bright unconfused sources, our flux density estimates are in good agreement with those produced via more traditional point source detection methods (SussExtractor; Savage & Oliver et al. 2006) by Smith et al. 2010. When compared to the measured number density of sources in the SPIRE bands, we show that our method allows the recovery of a larger fraction of faint sources than these traditional methods. However this completeness is heavily dependant on the relative depth of the existing 24 micron catalogues and SPIRE imaging. Using our deepest multi-wavelength dataset in GOODS-N, we estimate that the use of shallow 24 micron in our other fields introduces an incompleteness at faint levels of between 20-40 per cent at 250 micron.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

We outline the expected constraints on non-Gaussianity from the cosmic microwave background (CMB) with current and future experiments, focusing on both the third (f_{NL}) and fourth-order (g_{NL} and \tau_{NL}) amplitudes of the local configuration or non-Gaussianity. The experimental focus is the skewness (two-to-one) and kurtosis (two-to-two and three-to-one) power spectra from weighted maps. In adition to a measurement of \tau_{NL} and g_{NL} with WMAP 5-year data, our study provides the first forecasts for future constraints on g_{NL}. We describe how these statistics can be corrected for the mask and cut-sky through a window function, bypassing the need to compute linear terms that were introduced for the previous-generation non-Gaussianity statistics, such as the skewness estimator. We discus the ratio A_{NL} = \tau_{NL}/(6f_{NL}/5)^2 as an additional test of single-field inflationary models and discuss the physical significance of each statistic. Using these estimators with WMAP 5-Year V+W-band data out to l_{max}=600 we constrain the cubic order non-Gaussianity parameters \tau_{NL}, and g_{NL} and find -7.4 < g_{NL}/10^5 < 8.2 and -0.6 < \tau_{NL}/10^4 < 3.3 improving the previous COBE-based limit on \tau_{NL} < 10^8 nearly four orders of magnitude with WMAP.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We consider corrections to the underlying cosmology due to secondary contributions from weak gravitational lensing, the integrated Sachs-Wolfe effect, and the Sunyaev-Zel'dovich effect contained in the trispectrum. We incorporate these additional contributions to the covariance of a binned angular power spectrum of temperature anisotropies in the analysis of current and prospective data sets. Although recent experiments such as ACBAR and CBI are not particularly sensitive to these additional non-Gaussian effects, the interpretation of Planck and CMBPol anisotropy spectra will require an accounting of non-Gaussian covariance leading to a degradation in cosmological parameter estimates by up to 20% and 30%, respectively.

The reliability of infrared (IR) and ultraviolet (UV) emissions to measure star formation rates in galaxies is investigated for a large sample of galaxies observed with the SPIRE and PACS instruments on Herschel as part of the HerMES project. We build flux-limited 250 micron samples of sources at redshift z<1, cross-matched with the Spitzer/MIPS and GALEX catalogues. About 60 % of the Herschel sources are detected in UV. The total IR luminosities, L_IR, of the sources are estimated using a SED-fitting code that fits to fluxes between 24 and 500 micron. Dust attenuation is discussed on the basis of commonly-used diagnostics: the L_IR/L_UV ratio and the slope, beta, of the UV continuum. A mean dust attenuation A_UV of ~ 3 mag is measured in the samples. L_IR/L_UV is found to correlate with L_IR. Galaxies with L_IR > 10 ^{11} L_sun and 0.5< z<1 exhibit a mean dust attenuation A_UV about 0.7 mag lower than that found for their local counterparts, although with a large dispersion. Our galaxy samples span a large range of beta and L_IR/L_UV values which, for the most part, are distributed between the ranges defined by the relations found locally for starburst and normal star-forming galaxies. As a consequence the recipe commonly applied to local starbursts is found to overestimate the dust attenuation correction in our galaxy sample by a factor ~2-3 .

While limited to low spatial resolution, the next generation low-frequency radio interferometers that target 21 cm observations during the era of reionization and prior will have instantaneous fields-of-view that are many tens of square degrees on the sky. Predictions related to various statistical measurements of the 21 cm brightness temperature must then be pursued with numerical simulations of reionization with correspondingly large volume box sizes, of order 1000 Mpc on one side. We pursue a semi-numerical scheme to simulate the 21 cm signal during and prior to Reionization by extending a hybrid approach where simulations are performed by first laying down the linear dark matter density field, accounting for the non-linear evolution of the density field based on second-order linear perturbation theory as specified by the Zel'dovich approximation, and then specifying the location and mass of collapsed dark matter halos using the excursion-set formalism. The location of ionizing sources and the time evolving distribution of ionization field is also specified using an excursion-set algorithm. We account for the brightness temperature evolution through the coupling between spin and gas temperature due to collisions, radiative coupling in the presence of Lyman-alpha photons and heating of the intergalactic medium, such as due to a background of X-ray photons. The hybrid simulation method we present is capable of producing the required large volume simulations with adequate resolution in a reasonable time so a large number of realizations can be obtained with variations in assumptions related to astrophysics and background cosmology that govern the 21 cm signal.

We measure the angular correlation function, w(theta), from 0.5 to 30 arcminutes of detected sources in two wide fields of the Herschel Multi-tiered Extragalactic Survey (HerMES). Our measurements are consistent with the expected clustering shape from a population of sources that trace the dark matter density field, including non-linear clustering at arcminute angular scales arising from multiple sources that occupy the same dark matter halos. By making use of the halo model to connect the spatial clustering of sources to the dark matter halo distribution, we estimate source bias and halo occupation number for dusty sub-mm galaxies at z ~ 2. We find that sub-mm galaxies with 250 micron flux densities above 30 mJy reside in dark matter halos with mass above (5\pm4) x 10^12 M_sun, while (14\pm8)% of such sources appear as satellites in more massive halos.

Aims: We aim to demonstrate that the Herschel ATLAS (H-ATLAS) is suitable for a blind and unbiased survey for debris disks by identifying candidate debris disks associated with main sequence stars in the initial science demonstration field of the survey. We show that H-ATLAS reveals a population of far-infrared/sub-mm sources that are associated with stars or star-like objects on the SDSS main-sequence locus. We validate our approach by comparing the properties of the most likely candidate disks to those of the known population. Methods: We use a photometric selection technique to identify main sequence stars in the SDSS DR7 catalogue and a Bayesian Likelihood Ratio method to identify H-ATLAS catalogue sources associated with these main sequence stars. Following this photometric selection we apply distance cuts to identify the most likely candidate debris disks and rule out the presence of contaminating galaxies using UKIDSS LAS K-band images. Results: We identify 78 H-ATLAS sources associated with SDSS point sources on the main-sequence locus, of which two are the most likely debris disk candidates: H-ATLAS J090315.8 and H-ATLAS J090240.2. We show that they are plausible candidates by comparing their properties to the known population of debris disks. Our initial results indicate that bright debris disks are rare, with only 2 candidates identified in a search sample of 851 stars. We also show that H-ATLAS can derive useful upper limits for debris disks associated with Hipparcos stars in the field and outline the future prospects for our debris disk search programme.

We have determined the luminosity function of 250um-selected galaxies detected in the ~14 sq.deg science demonstration region of the Herschel-ATLAS project out to a redshift of z=0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250um luminosity density out to z=0.2 where it is 3.6+1.4-0.9 times higher than the local value.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Local luminosity functions are fundamental benchmarks for high-redshift galaxy formation and evolution studies as well as for models describing these processes. Determining the local luminosity function in the submillimeter range can help to better constrain in particular the bolometric luminosity density in the local Universe, and Herschel offers the first opportunity to do so in an unbiased way by imaging large sky areas at several submillimeter wavelengths. We present the first Herschel measurement of the submillimeter 0<z<0.2 local luminosity function and infrared bolometric (8-1000 $\mu$m) local luminosity density based on SPIRE data from the HerMES Herschel Key Program over 14.7 deg^2. Flux measurements in the three SPIRE channels at 250, 350 and 500 \mum are combined with Spitzer photometry and archival data. We fit the observed optical-to-submillimeter spectral energy distribution of SPIRE sources and use the 1/V_{max} estimator to provide the first constraints on the monochromatic 250, 350 and 500 \mum as well as on the infrared bolometric (8-1000 \mum) local luminosity function based on Herschel data. We compare our results with modeling predictions and find a slightly more abundant local submillimeter population than predicted by a number of models. Our measurement of the infrared bolometric (8-1000 \mum) local luminosity function suggests a flat slope at low luminosity, and the inferred local luminosity density, 1.31_-0.21^+0.24 x 10^8 Lsun Mpc^-3, is consistent with the range of values reported in recent literature.

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, and GOODS-N fields in three photometric bands at 250, 350 and 500 micrometers. We compare these with expectations from the literature and discuss associated uncertainties and biases in the SPIRE data. We identify a 500 micrometer flux limited selection of sources from the HerMES point source catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 1<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S_250/S_350<0.8 at fluxes above 50 mJy as observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies.

Aims. The Herschel-ATLAS survey (H-ATLAS) will be the largest area survey to be undertaken by the Herschel Space Observatory. It will cover 550 sq. deg. of extragalactic sky at wavelengths of 100, 160, 250, 350 and 500 microns when completed, reaching flux limits (5 sigma) from 32 to 145mJy. We here present galaxy number counts obtained for SPIRE observations of the first ~14 sq. deg. observed at 250, 350 and 500 microns. Methods. Number counts are a fundamental tool in constraining models of galaxy evolution. We use source catalogs extracted from the H-ATLAS maps as the basis for such an analysis. Correction factors for completeness and flux boosting are derived by applying our extraction method to model catalogs and then applied to the raw observational counts. Results. We find a steep rise in the number counts at flux levels of 100-200mJy in all three SPIRE bands, consistent with results from BLAST. The counts are compared to a range of galaxy evolution models. None of the current models is an ideal fit to the data but all ascribe the steep rise to a population of luminous, rapidly evolving dusty galaxies at moderate to high redshift.

We present a derivation of the star formation rate per comoving volume of quasar host galaxies, derived from stacking analyses of far-infrared to mm-wave photometry of quasars with redshifts 0<z<6 and absolute I-band magnitudes -22>I_AB>-32. We use the science demonstration observations of the first ~16 deg^2 from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) in which there are 240 quasars from the Sloan Digital Sky Survey (SDSS) and a further 171 from the 2dF-SDSS LRG and QSO (2SLAQ) survey. We supplement this data with a compilation of data from IRAS, ISO, Spitzer, SCUBA and MAMBO. H-ATLAS alone statistically detects the quasars in its survey area at >5sigma at 250, 350 and 500um. From the compilation as a whole we find striking evidence of downsizing in quasar host galaxy formation: low-luminosity quasars with absolute magnitudes in the range -22>I_AB>-24 have a comoving star formation rate (derived from 100um rest-frame luminosities) peaking between redshifts of 1 and 2, while high-luminosity quasars with I_AB<-26 have a maximum contribution to the star formation density at z~3. The volume-averaged star formation rate of -22>I_AB>-24 quasars evolves as (1+z)^{2.3 +/- 0.7} at z<2, but the evolution at higher luminosities is much faster reaching (1+z)^{10 +/- 1} at -26>I_AB>-28. We tentatively interpret this as a combination of a declining major merger rate with time and gas consumption reducing fuel for both black hole accretion and star formation.

We present colour-colour diagrams of detected sources in the Herschel-ATLAS Science Demonstration Field from 100 to 500 microns using both PACS and SPIRE. We fit isothermal modified black bodies to the spectral energy distribution (SED) to extract the dust temperature of sources with counterparts in Galaxy And Mass Assembly (GAMA) or SDSS surveys with either a spectroscopic or a photometric redshift. For a subsample of 330 sources detected in at least three FIR bands with a significance greater than 3 $\sigma$, we find an average dust temperature of $(28 \pm 8)$K. For sources with no known redshift, we populate the colour-colour diagram with a large number of SEDs generated with a broad range of dust temperatures and emissivity parameters, and compare to colours of observed sources to establish the redshift distribution of this sample. For another subsample of 1686 sources with fluxes above 35 mJy at 350 microns and detected at 250 and 500 microns with a significance greater than 3$\sigma$, we find an average redshift of $2.2 \pm 0.6$.

We present measurements of the angular correlation function of galaxies selected from the first field of the H-ATLAS survey. Careful removal of the background from galactic cirrus is essential, and currently dominates the uncertainty in our measurements. For our 250 micron-selected sample we detect no significant clustering, consistent with the expectation that the 250 micron-selected sources are mostly normal galaxies at z<~ 1. For our 350 micron and 500 micron-selected samples we detect relatively strong clustering with correlation amplitudes A of 0.2 and 1.2 at 1', but with relatively large uncertainties. For samples which preferentially select high redshift galaxies at z~2-3 we detect significant strong clustering, leading to an estimate of r_0 ~ 7-11 h^{-1} Mpc. The slope of our clustering measurements is very steep, delta~2. The measurements are consistent with the idea that sub-mm sources consist of a low redshift population of normal galaxies and a high redshift population of highly clustered star-bursting galaxies.

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts at <100 mJy. We have directly resolved 15% of the infrared extra-galactic background at the wavelength near where it peaks.

We have carried out two extremely deep surveys with SPIRE, one of the two cameras on Herschel, at 250 microns, close to the peak of the far-infrared background. We have used the results to investigate the evolution of the rest-frame 250-micron luminosity function out to z=2. We find evidence for strong evolution out to a redshift of around 1 but evidence for at most weak evolution beyond this redshift. Our results suggest that a significant part of the stars and metals in the Universe today were formed at z<1.4 in spiral galaxies.

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.

We use Herschel PACS and SPIRE observations of the edge-on spiral galaxy UGC 4754, taken as part of the H-ATLAS SDP observations, to investigate the dust energy balance in this galaxy. We build detailed SKIRT radiative models based on SDSS and UKIDSS maps and use these models to predict the far-infrared emission. We find that our radiative transfer model underestimates the observed FIR emission by a factor of two to three. Similar discrepancies have been found for other edge-on spiral galaxies based on IRAS, ISO, and SCUBA data. Thanks to the good sampling of the SED at FIR wavelengths, we can rule out an underestimation of the FIR emissivity as the cause for this discrepancy. Instead we support highly obscured star formation that contributes little to the optical extinction as a more probable explanation.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

We measure two higher-order power spectra involving weighted cubic and squared temperature anisotropy maps from WMAP 5-year data to study the trispectrum generated by primordial non-Gaussianity. Using these measurements combined with Gaussian and noise simulations, we constrain the cubic order non-Gaussianity parameters \tau_{NL}, and g_{NL}. With V+W-band data out to l_{max}=600, we find -7.4 < g_{\rm NL}/10^5 < 8.2 and -0.6 < \tau_{\rm NL}/10^4 < 3.3 improving the previous COBE-based limit on \tau_{\rm NL} < 10^8 nearly four orders of magnitude with WMAP. We find that the ratio of trispectrum to bispectrum amplitude as captured by the ratio of \tau_{\rm Nl}/(6f_{\rm NL}/5)^2 ranges from -3 to 21 at the 95% confidence level.

In this report we discuss the impact of polarized foregrounds on a future CMBPol satellite mission. We review our current knowledge of Galactic polarized emission at microwave frequencies, including synchrotron and thermal dust emission. We use existing data and our understanding of the physical behavior of the sources of foreground emission to generate sky templates, and start to assess how well primordial gravitational wave signals can be separated from foreground contaminants for a CMBPol mission. At the estimated foreground minimum of ~100 GHz, the polarized foregrounds are expected to be lower than a primordial polarization signal with tensor-to-scalar ratio r=0.01, in a small patch (~1%) of the sky known to have low Galactic emission. Over 75% of the sky we expect the foreground amplitude to exceed the primordial signal by about a factor of eight at the foreground minimum and on scales of two degrees. Only on the largest scales does the polarized foreground amplitude exceed the primordial signal by a larger factor of about 20. The prospects for detecting an r=0.01 signal including degree-scale measurements appear promising, with 5 sigma_r ~0.003 forecast from multiple methods. A mission that observes a range of scales offers better prospects from the foregrounds perspective than one targeting only the lowest few multipoles. We begin to explore how optimizing the composition of frequency channels in the focal plane can maximize our ability to perform component separation, with a range of typically 40 < nu < 300 GHz preferred for ten channels. Foreground cleaning methods are already in place to tackle a CMBPol mission data set, and further investigation of the optimization and detectability of the primordial signal will be useful for mission design.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

We measure the skewness power spectrum of the Cosmic Microwave Background (CMB) anisotropies optimized for a detection of the secondary bispectrum generated by the correlation of the CMB lensing potential with integrated Sachs-Wolfe effect and the Sunyaev-Zel'dovich effect. The covariance of our measurements is generated by Monte-Carlo simulations of Gaussian CMB fields with noise properties consistent with Wilkinson Microwave Anisotropy Probe (WMAP) 5-year data. When interpreting multi-frequency measurements we also take into account the confusion resulting from unresolved radio point sources. We analyze Q, V and W-band WMAP 5-year raw and foreground-cleaned maps using the KQ75 mask out to l_max=600. We find no significant evidence for a non-zero non-Gaussian signal from the lensing-secondary correlation in all three bands and we constrain the overall amplitude of the cross power spectrum between CMB lensing potential and the sum of SZ and ISW fluctuations to be 0.42 \pm 0.86 and 1.19 \pm 0.86 in combined V and W-band raw and foreground-cleaned maps provided by the WMAP team, respectively. The point source amplitude at the bispectrum level measured with this skewness power spectrum is higher than previous measurements of point source non-Gaussianity. We also consider an analysis where we also account for the primordial non-Gaussianity in addition to lensing-secondary bispectrum and point sources. The focus of this paper is on secondary anisotropies. Consequently the estimator is not optimised for primordial non-Gaussianity and the limit we find on local non-Gaussianity from the foreground-cleaned V+W maps is f_NL = -13 \pm 62, when marginalized over point sources and lensing-ISW/SZ constributions to the total bispectrum.

The Sunyaev-Zel'dovich (SZ) effect is the inverse Compton-scattering of cosmic microwave background (CMB) photons by hot electrons in the intervening gas throughout the universe. The effect has a distinct spectral signature that allows its separation from other signals in multifrequency CMB datasets. Using CMB anisotropies measured at three frequencies by the BOOMERanG 2003 flight we constrain SZ fluctuations in the 10 arcmin to 1 deg angular range. Propagating errors and potential systematic effects through simulations, we obtain an overall upper limit of 15.3 uK (2 sigma) for rms SZ fluctuations in a broad bin between multipoles of of 250 and 1200 at the Rayleigh-Jeans (RJ) end of the spectrum. When combined with other CMB anisotropy and SZ measurements, we find that the local universe normalization of the density perturbations is sigma-8(SZ) < 0.96 at the 95% confidence level, consistent with sigma-8 determined from primordial perturbations.

We constrain the primordial non-Gaussianity parameter of the local model f_{NL} using the skewness power spectrum associated with the two-to-one cumulant correlator of cosmic microwave background temperature anisotropies. This bispectrum-related power spectrum was constructed after weighting the temperature map with the appropriate window functions to form an estimator that probes the multipolar dependence of the underlying bispectrum associated with the primordial non-Gaussianity. We also estimate a separate skewness power spectrum sensitive more strongly to unresolved point sources. When compared to previous attempts at measuring the primordial non-Gaussianity with WMAP data, our estimators have the main advantage that we do not collapse information to a single number. When model fitting the two-to-one skewness power spectrum we make use of bispectra generated by the primordial non-Gaussianity, radio point sources, and lensing-secondary correlation. We analyze Q, V and W-band WMAP 5-year data using the KQ75 mask out to l_{max}=600. Using V and W-band data and marginalizing over model parameters related to point sources and lensing-secondary bispectrum, our overall and preferred constraint on f_{NL} is 11.0\pm23.7 at the 68% confidence level (-36.4 < f_{NL} < 58.4 at 95% confidence). We find no evidence for a non-zero value of f_{NL} even marginally at the 1\sigma level.

Measurements of Cosmic Microwave Background (CMB) anisotropy have served as the best experimental probe of the early universe to date. The inflationary paradigm, inspired in part by the extreme isotropy of the CMB, is now a cornerstone in modern cosmology. Inflation has passed a series of rigorous experimental tests, but we still do not understand the physical mechanism or energy scale behind inflation. A general prediction of inflation and one that can provide certain insights into inflationary physics is a background of primordial gravitational waves. These perturbations leave a distinct signature in the CMB B-modes of polarization. The EPIC (Experimental Probe of Inflationary Cosmology) study team has investigated several CMB polarization mission concepts to carry out a definitive measurement of the inflationary B-mode polarization spectrum. In this report we study a mission with an aperture intermediate between the two missions discussed in our previous report. EPIC-IM's increased aperture allows access to a broader science case than the small EPIC-Low Cost mission. In addition to the search for inflationary gravitational waves, the increase aperture allows us to mine the scale polarization and lensing shear polarization signals down to cosmological limits, so that we extract virtually all the cosmological information available from the CMB. In addition, a modest number of channels operating at higher frequencies allows for an all-sky measurement of polarized Galactic dust, which will provide a rich dataset for Galactic science related to magnetic fields. Using a combination of a large sensitivity focal plane with a new optical design, and an efficient 4K mechanical cooler, EPIC-IM realizes higher sensitivity than EPIC-Comprehensive Science mission.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

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 mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

We apply spherical needlets to the Wilkinson Microwave Anisotropy Probe 5-year cosmic microwave background (CMB) dataset, to search for imprints of non-isotropic features in the CMB sky. We use the needlets localization properties to resolve peculiar features in the CMB sky and to study how these features contribute to the anisotropy power spectrum of the CMB. In addition to the now well-known "cold spot" of the CMB map in the southern hemisphere, we also find two hot spots at greater than 99% confidence level, again in the southern hemisphere and closer to the Galactic plane. While the cold spot contributes to the anisotropy power spectrum in the multipoles between l=6 to l=33, the hot spots are found to be dominating the anisotropy power in the range between l=6 and l=18. Masking both the cold and the two hot spots results in a reduction by about 15% in the amplitude of the angular power spectrum of CMB around l=10. The resulting changes to the cosmological parameters when the power spectrum is estimated masking these features (in addition to the WMAP team's KQ85 mask) are within the 1$\sigma$ errors published with the WMAP mask only. We also study the asymmetry between the angular power spectra evaluated on the northern and southern hemispheres. When the features detected by needlets are masked, we find that the difference in the power, measured in terms of the anisotropy variance between l=4 and l=18, is reduced by a factor 2. We make available a mask related to needlet features for more detailed studies on asymmetries in the CMB anisotropy sky.

We measure several properties of the reionization process and the corresponding low-frequency 21-cm signal associated with the neutral hydrogen distribution, using a large volume, high resolution simulation of cosmic reionization. The brightness temperature of the 21-cm signal is derived by post-processing this numerical simulation with a semi-analytical prescription. Our study extends to high redshifts (z ~ 25) where, in addition to collisional coupling, our post-processed simulations take into account the inhomogeneities in the heating of the neutral gas by X-rays and the effect of an inhomogeneous Lya radiation field. Unlike the well-studied case where spin temperature is assumed to be significantly greater than the temperature of the cosmic microwave background due to uniform heating of the gas by X-rays, spatial fluctuations in both the Lya radiation field and X-ray intensity impact predictions related to the brightness temperature at z > 10, during the early stages of reionization and gas heating. The statistics of the 21-cm signal from our simulation are then compared to existing analytical models in the literature and we find that these analytical models provide a reasonably accurate description of the 21-cm power spectrum at z < 10. Such an agreement is useful since analytical models are better suited to quickly explore the full astrophysical and cosmological parameter space relevant for future 21-cm surveys. We find, nevertheless, non-negligible differences that can be attributed to differences in the inhomogeneous X-ray heating and Lya coupling at z > 10 and, with upcoming interferometric data, these differences in return can provide a way to better understand the astrophysical processes during reionization.

Recent work suggests that Type Ia supernovae (SNe) are composed of two distinct populations: prompt and delayed. By explicitly incorporating properties of host galaxies, it may be possible to target and eliminate systematic differences between these two putative populations. However, any resulting {\em post}-calibration shift in luminosity between the components will cause a redshift-dependent systematic shift in the Hubble diagram. Utilizing an existing sample of 192 SNe Ia, we find that the average luminosity difference between prompt and delayed SNe is constrained to be $(4.5 \pm 8.9)%$. If the absolute difference between the two populations is 0.025 mag, and this is ignored when fitting for cosmological parameters, then the dark energy equation of state (EOS) determined from a sample of 2300 SNe Ia is biased at $\sim1\sigma$. By incorporating the possibility of a two-population systematic, this bias can be eliminated. However, assuming no prior on the strength of the two-population effect, the uncertainty in the best-fit EOS is increased by a factor of 2.5, when compared to the equivalent sample with no underlying two-population systematic. To avoid introducing a bias in the EOS parameters, or significantly degrading the measurement accuracy, it is necessary to control the post-calibration luminosity difference between prompt and delayed SN populations to better than 0.025 mag.

Our ignorance of the dark energy is generally described by a two-parameter equation of state. In these approaches a particular {\it ad hoc} functional form is assumed, and only two independent parameters are incorporated. We propose a model-independent, multi-parameter approach to fitting the dark energy, and show that next-generation surveys will constrain the equation of state in three or more independent redshift bins to better than 10%. Future knowledge of the dark energy will surpass two numbers (e.g., [$w_0$,$w_1$] or [$w_0$,$w_a$]), and we propose a more flexible approach to the analysis of present and future data.

The faint radio point sources that are unresolved in cosmic microwave background (CMB) anisotropy maps are likely to be a biased tracer of the large-scale structure dark matter distribution. While the shot-noise contribution to the angular power spectrum of unresolved radio point sources is included either when optimally constructing the CMB angular power spectrum, as with WMAP data, or when extracting cosmological parameters, we suggest that clustering part of the point source power spectrum should also be included. This is especially necessary at high frequencies above 150 GHz, where the clustering of far-IR sources is expected to dominate the shot-noise level of the angular power spectrum at tens of arcminute angular scales of both radio and sub-mm sources. We make an estimate of source clustering of unresolved radio sources in both WMAP and ACBAR, and marginalize over the amplitude of source clustering in each CMB data set when model fitting for cosmological parameters. For the combination of WMAP 5-year data and ACBAR, we find that the spectral index changes from the value of $0.963 \pm 0.014$ to $0.959 \pm 0.014$ (at 68% c.l.) when the clustering power spectrum of point sources is included in model fits. While we find that the differences are marginal with and without source clustering in current data, it may be necessary to account for source clustering with future datasets such as Planck, especially to properly model fit anisotropies at arcminute angular scales. If clustering is not accounted and point sources are modeled with a shot-noise only out to $l \sim 2000$, the spectral index will be biased by about 1.5$\sigma$.

The gravitational magnification and demagnification of Type Ia supernovae (SNe) modify their positions on the Hubble diagram, shifting the distance estimates from the underlying luminosity-distance relation. This can introduce a systematic uncertainty in the dark energy equation of state (EOS) estimated from SNe, although this systematic is expected to average away for sufficiently large data sets. Using mock SN samples over the redshift range $0 < z \leq 1.7$ we quantify the lensing bias. We find that the bias on the dark energy EOS is less than half a percent for large datasets ($\gtrsim$ 2,000 SNe). However, if highly magnified events (SNe deviating by more than 2.5$\sigma$) are systematically removed from the analysis, the bias increases to $\sim$ 0.8%. Given that the EOS parameters measured from such a sample have a 1$\sigma$ uncertainty of 10%, the systematic bias related to lensing in SN data out to $z \sim 1.7$ can be safely ignored in future cosmological measurements.

Dusty, starforming galaxies and active galactic nuclei that contribute to the integrated background intensity at far-infrared wavelengths trace the large-scale structure. Below the point source detection limit, correlations in the large-scale structure lead to clustered anisotropies in the unresolved component of the far-infrared background (FIRB). The angular power spectrum of the FIRB anisotropies could be measured in large-area surveys with the Spectral and Photometric Imaging Receiver (SPIRE) on the upcoming Herschel observatory. To study statistical properties of these anisotropies, the confusion from foreground Galactic dust emission needs to be reduced even in the ``cleanest'' regions of the sky.The multi-frequency coverage of SPIRE allows the foreground dust to be partly separated from the extragalactic background composed of dusty starforming galaxies as well as faint normal galaxies. The separation improves for fields with sizes greater than a few hundred square degrees and when combined with Planck data. We show that an area of about $\sim$ 400 degrees$^2$ observed for about 1000 hours with Herschel-SPIRE and complemented by Planck provides maximal information on the anisotropy power spectrum. We discuss the scientific studies that can be done with measurements of the unresolved FIRB anisotropies including a determination of the large scale bias and the small-scale halo occupation distribution of FIRB sources with fluxes below the point-source detection level.

B-modes of Cosmic Microwave Background (CMB) polarization can be created by a primordial gravitational wave background. If this background was created by Inflation, then the amplitude of the polarization signal is proportional the energy density of the universe during inflation. The primordial signal will be contaminated by polarized foregrounds including dust and synchrotron emission within the galaxy. In light of the WMAP polarization maps, we consider the ability of several hypothetical CMB polarization experiments to separate primordial CMB B-mode signal from galactic foregrounds. We also study the optimization of a CMB experiment with a fixed number of detectors in the focal plane to determine how the detectors should be distributed in different frequency bands to minimize foreground confusion. We show that the optimal configuration requires observations in at least 5 channels spread over the frequency range between 30 GHz and 500 GHz with substantial coverage around 150 GHz. If a low-resolution space experiment using 1000 detectors to reach a noise level of about 1000 nK$^2$ concentrates on roughly 66% of the sky with the least foreground contamination the minimum detectable level of the tensor-to-scalar ratio would be about 0.002 at the 99% confidence level for an optical depth of 0.1.

Archeops is a balloon--borne experiment widely inspired by the Planck satellite and by its High Frequency Instrument (HFI). It is mainly dedicated to measure the Cosmic Microwave Background (CMB) temperature anisotropies at high angular resolution (about 12 arcminutes) over a large fraction of the sky (around 30 %) in the millimetre and submillimetre range at 143, 217, 353 and 545 GHz. Further, the Archeops 353 GHz channel consists of three pairs of polarized sensitive bolometers designed to detect the polarized diffuse emission of Galactic dust. We present in this paper the update of the instrumental setup as well as the inflight performance for the last Archeops flight campaign in February 2002 from Kiruna (Sweden). We also describe the processing and analysis of the Archeops time ordered data for that campaign which lead to the measurement of the CMB anisotropies power spectrum in the multipole range l=10-700 (Benoit et al. 2003a, Tristram et al. 2005) and to the first measurement of the dust polarized emission at large angular scales and its polarized. We present maps of 30 % of the sky of the Galactic emission, including the Galactic plane, in the four Archeops channels at 143, 217, 353 and 545 GHz and maps of the CMB anisotropies at 143 and 217 GHz. These are the firstever available sub--degree resolution maps in the millimetre and submillimetre range of the large angular-scales Galactic dust diffuse emission and CMB temperature anisotropies respectively.

Compton scattering of Cosmic Microwave Background (CMB) photons on galaxy cluster electrons produces a linear polarization, which contains some information on the local quadrupole at the cluster location. We use N-body simulations to create, for the first time, maps of this polarization signal. We then look at the different properties of the polarization with respect to the cluster position and redshift.

We present improved results on the measurement of the angular power spectrum of the Cosmic Microwave Background (CMB) temperature anisotropies using the data from the last Archeops flight. This refined analysis is obtained by using the 6 most sensitive photometric pixels in the CMB bands centered at 143 and 217 GHz and 20% of the sky, mostly clear of foregrounds. Using two different cross-correlation methods, we obtain very similar results for the angular power spectrum. Consistency checks are performed to test the robustness of these results paying particular attention to the foreground contamination level which remains well below the statistical uncertainties. The multipole range from l=10 to l=700 is covered with 25 bins, confirming strong evidence for a plateau at large angular scales (the Sachs-Wolfe plateau) followed by two acoustic peaks centered around l=220 and l=550 respectively. These data provide an independent confirmation, obtained at different frequencies, of the WMAP first year results.

We present the first measurement of temperature and polarization angular power spectra of the diffuse emission of Galactic dust at 353 GHz as seen by Archeops on 20% of the sky. The temperature angular power spectrum is compatible with that provided by the extrapolation to 353 GHz of IRAS and DIRBE maps using \cite{fds} model number 8. For Galactic latitudes $|b| \geq 5$ deg we report a 4 sigma detection of large scale ($3\leq \ell \leq 8$) temperature-polarization cross-correlation $(\ell+1)C_\ell^{TE}/2\pi = 76\pm 21 \mu\rm{K_{RJ}}^2$ and set upper limits to the $E$ and $B$ modes at $11 \mu\rm{K_{RJ}}^2$. For Galactic latitudes $|b| \geq 10$ deg, on the same angular scales, we report a 2 sigma detection of temperature-polarization cross-correlation $(\ell+1)C_\ell^{TE}/2\pi = 24\pm 13 \mu\rm{K_{RJ}}^2$. These results are then extrapolated to 100 GHz to estimate the contamination in CMB measurements by polarized diffuse Galactic dust emission. The $TE$ signal is then $1.7\pm0.5$ and $0.5\pm0.3 \mu\rm{K^2_{CMB}}$ for $|b| \geq 5$ and 10 deg. respectively. The upper limit on $E$ and $B$ becomes $0.2 \mu\rm{K^2_{CMB}} (2\sigma)$. If polarized dust emission at higher Galactic latitude cuts is similar to the one we report here, then dust polarized radiation will be a major foreground for determining the polarization power spectra of the CMB at high frequencies above 100 GHz.

Several recent papers have studied lensing of the CMB by large-scale structures, which probes the projected matter distribution from $z=10^3$ to $z\simeq 0$. This interest is motivated in part by upcoming high resolution, high sensitivity CMB experiments, such as APEX/SZ, ACT, SPT or Planck, which should be sensitive to lensing. In this paper we examine the reconstruction of the large-scale dark matter distribution from lensed CMB temperature anisotropies. We go beyond previous work in using numerical simulations to include higher order, non-Gaussian effects and study how well the quadratic estimator of \cite{Hu01a} is able to recover the input field. We also study contamination by kinetic Sunyaev-Zel'dovich signals, which is spectrally indistinguishable from lensed CMB anisotropies. We finish by estimating the sensitivity of the previously cited experiments.

We investigate what the lensing information contained in high resolution, low noise CMB temperature maps can teach us about cluster mass profiles. We create lensing fields and Sunyaev-Zel'dovich effect maps from N-body simulations and apply them to primary CMB anisotropies modeled as a Gaussian random field. We examine the success of several techniques of cluster mass reconstruction using CMB lensing information, and make an estimate of the observational requirements necessary to achieve a satisfactory result.

We present a method designed to estimate the noise power spectrum in the time domain for CMB experiments. The noise power spectrum is extracted from the time ordered data avoiding the contamination coming from sky signal and accounting the pixellisation of the signal and the projection of the noise when making intermediate sky projections. This method is simple to implement and relies on Monte-Carlo simulations, it runs on a simple desk computer. We also propose a trick for filtering data before making coadded maps in order to avoid ringing due to the presence of signal in the timelines. These algorithms were succesfully tested on Archeops data.

We present the first determination of the Galactic polarized emission at 353 GHz by Archeops. The data were taken during the Arctic night of February 7, 2002 after the balloon--borne instrument was launched by CNES from the Swedish Esrange base near Kiruna. In addition to the 143 GHz and 217 GHz frequency bands dedicated to CMB studies, Archeops had one 545 GHz and six 353 GHz bolometers mounted in three polarization sensitive pairs that were used for Galactic foreground studies. We present maps of the I, Q, U Stokes parameters over 17% of the sky and with a 13 arcmin resolution at 353 GHz (850 microns). They show a significant Galactic large scale polarized emission coherent on the longitude ranges [100, 120] and [180, 200] deg. with a degree of polarization at the level of 4-5%, in agreement with expectations from starlight polarization measurements. Some regions in the Galactic plane (Gem OB1, Cassiopeia) show an even stronger degree of polarization in the range 10-20%. Those findings provide strong evidence for a powerful grain alignment mechanism throughout the interstellar medium and a coherent magnetic field coplanar to the Galactic plane. This magnetic field pervades even some dense clouds. Extrapolated to high Galactic latitude, these results indicate that interstellar dust polarized emission is the major foreground for PLANCK-HFI CMB polarization measurement.

We analyze the cosmological constraints that Archeops places on adiabatic cold dark matter models with passive power-law initial fluctuations. Because its angular power spectrum has small bins in l and large l coverage down to COBE scales, Archeops provides a precise determination of the first acoustic peak in terms of position at multipole l_peak=220 +- 6, height and width. An analysis of Archeops data in combination with other CMB datasets constrains the baryon content of the Universe, Omega(b)h^2 = 0.022 (+0.003,-0.004), compatible with Big-Bang nucleosynthesis and with a similar accuracy. Using cosmological priors obtainedfrom recent non-CMB data leads to yet tighter constraints on the total density, e.g. Omega(tot)=1.00 (+0.03,-0.02) using the HST determination of the Hubble constant. An excellent absolute calibration consistency is found between Archeops and other CMB experiments, as well as with the previously quoted best fit model.The spectral index n is measured to be 1.04 (+0.10,-0.12) when the optical depth to reionization, tau, is allowed to vary as a free parameter, and 0.96 (+0.03,-0.04) when tau is fixed to zero, both in good agreement with inflation.

We present a determination by the Archeops experiment of the angular power spectrum of the cosmic microwave background anisotropy in 16 bins over the multipole range l=15-350. Archeops was conceived as a precursor of the Planck HFI instrument by using the same optical design and the same technology for the detectors and their cooling. Archeops is a balloon-borne instrument consisting of a 1.5 m aperture diameter telescope and an array of 21 photometers maintained at ~100 mK that are operating in 4 frequency bands centered at 143, 217, 353 and 545 GHz. The data were taken during the Arctic night of February 7, 2002 after the instrument was launched by CNES from Esrange base (Sweden). The entire data cover ~ 30% of the sky.This first analysis was obtained with a small subset of the dataset using the most sensitive photometer in each CMB band (143 and 217 GHz) and 12.6% of the sky at galactic latitudes above 30 degrees where the foreground contamination is measured to be negligible. The large sky coverage and medium resolution (better than 15 arcminutes) provide for the first time a high signal-to-noise ratio determination of the power spectrum over angular scales that include both the first acoustic peak and scales probed by COBE/DMR. With a binning of Delta(l)=7 to 25 the error bars are dominated by sample variance for l below 200. A companion paper details the cosmological implications.

Archeops, a balloon-borne experiment, will provide a measurement of CMB anistropies from large to small angular scale thanks to its large sky coverage (30%), its high angular resolution (10 arcminutes), and its high signal-to-noise ratio due to high sensitivity 100 mK cooled bolometers. We will therefore be able to put strong constraints on the value of the cosmological parameters. Archeops flew already twice, once in Sicily for a technical flight in 1999, and once from Sweden for its first scientific flight in January 2001. I describe here Archeops'main characteristics, the preliminary results from the scientific flight, the expected precision of this flight for power spectrum measurements, and perspectives for the next flights this winter.

Archeops is a balloon-borne instrument dedicated to measuring cosmic microwave background (CMB) temperature anisotropies at high angular resolution (8 arcminutes) over a large fraction (25%) of the sky in the millimetre domain. Based on Planck High Frequency Instrument (HFI) technology, cooled bolometers (0.1 K) scan the sky in total power mode with large circles at constant elevation. During the course of a 24-hour Arctic-night balloon flight, Archeops will observe a complete annulus on the sky in four frequency bands centered at 143, 217, 353 and 545 GHz with an expected sensitivity to CMB fluctuations of \~100muK for each of the 90 thousand 20 arcminute average pixels. We describe the instrument and its performance obtained during a test flight from Trapani (Sicily) to Spain in July 1999.