The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

Delensing is an increasingly important technique to reverse the gravitational lensing of the cosmic microwave background (CMB) and thus reveal primordial signals the lensing may obscure. We present a first demonstration of delensing on Planck temperature maps using the cosmic infrared background (CIB). Reversing the lensing deflections in Planck CMB temperature maps using a linear combination of the 545 and 857GHz maps as a lensing tracer, we find that the lensing effects in the temperature power spectrum are reduced in a manner consistent with theoretical expectations. In particular, the characteristic sharpening of the acoustic peaks of the temperature power spectrum resulting from successful delensing is detected at a significance of 16$\rm{\sigma}$, with an amplitude of $A_{\rm{delens}} = 1.12 \pm 0.07$ relative to the expected value of unity. This first demonstration on data of CIB delensing, and of delensing techniques in general, is significant because lensing removal will soon be essential for achieving high-precision constraints on inflationary B-mode polarization.

We present new measurements of the power spectra of the cosmic infrared background (CIB) anisotropies using the Planck 2015 full-mission HFI data at 353, 545, and 857 GHz over 20000 square degrees. We use techniques similar to those applied for the cosmological analysis of Planck, subtracting dust emission at the power spectrum level. Our analysis gives stable solutions for the CIB power spectra with increasing sky coverage up to about 50% of the sky. These spectra agree well with Hi cleaned spectra from Planck measured on much smaller areas of sky with low Galactic dust emission. At 545 and 857 GHz our CIB spectra agree well with those measured from Herschel data. We find that the CIB spectra at l > 500 are well fitted by a power-law model for the clustered CIB, with a shallow index {\gamma}^cib = 0.53\pm0.02. This is consistent with the CIB results at 217 GHz from the cosmological parameter analysis of Planck. We show that a linear combination of the 545 and 857 GHz Planck maps is dominated by CIB fluctuations at multipoles l > 300.

By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present a significant detection of baryons participating in large-scale bulk flows around central galaxies (CGs) at redshift $z\approx 0.1$. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the CGC (Central Galaxy Catalogue) samples extracted from Sloan Digital Sky Survey (DR7) data. For the foreground-cleaned maps, we find $1.8$-$2.5\sigma$ detections of the kSZ signal, which are consistent with the kSZ evidence found in individual Planck raw frequency maps, although lower than found in the WMAP-9yr W band ($3.3\sigma$). We further reconstruct the peculiar velocity field from the CG density field, and compute for the first time the cross-correlation function between kSZ temperature fluctuations and estimates of CG radial peculiar velocities. This correlation function yields a $3.0$-$3.7$$\sigma$ detection of the peculiar motion of extended gas on Mpc scales, in flows correlated up to distances of 80-100 $h^{-1}$ Mpc. Both the pairwise momentum estimates and kSZ temperature-velocity field correlation find evidence for kSZ signatures out to apertures of 8 arcmin and beyond, corresponding to a physical radius of $> 1$ Mpc, more than twice the mean virial radius of halos. This is consistent with the predictions from hydro simulations that most of the baryons are outside the virialized halos. We fit a simple model, in which the temperature-velocity cross-correlation is proportional to the signal seen in a semi-analytic model built upon N-body simulations, and interpret the proportionality constant as an "effective" optical depth to Thomson scattering. We find $\tau_T=(1.4\pm0.5)\times 10^{-4}$; the simplest interpretation of this measurement is that much of the gas is in a diffuse phase, which contributes little signal to X-ray or thermal SZ observations.

We use Planck data to detect the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) effect and the infrared emission from the galaxies that make up the the cosmic infrared background (CIB). We first perform a stacking analysis towards Planck-confirmed galaxy clusters. We detect infrared emission produced by dusty galaxies inside these clusters and demonstrate that the infrared emission is about 50% more extended than the tSZ effect. Modelling the emission with a Navarro--Frenk--White profile, we find that the radial profile concentration parameter is $c_{500} = 1.00^{+0.18}_{-0.15}$. This indicates that infrared galaxies in the outskirts of clusters have higher infrared flux than cluster-core galaxies. We also study the cross-correlation between tSZ and CIB anisotropies, following three alternative approaches based on power spectrum analyses: (i) using a catalogue of confirmed clusters detected in Planck data; (ii) using an all-sky tSZ map built from Planck frequency maps; and (iii) using cross-spectra between Planck frequency maps. With the three different methods, we detect the tSZ-CIB cross-power spectrum at significance levels of (i) 6 $\sigma$, (ii) 3 $\sigma$, and (iii) 4 $\sigma$. We model the tSZ-CIB cross-correlation signature and compare predictions with the measurements. The amplitude of the cross-correlation relative to the fiducial model is $A_{\rm tSZ-CIB}= 1.2\pm0.3$. This result is consistent with predictions for the tSZ-CIB cross-correlation assuming the best-fit cosmological model from Planck 2015 results along with the tSZ and CIB scaling relations.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

Using \Planck\ data combined with the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we address the study of peculiar motions by searching for evidence of the kinetic Sunyaev-Zeldovich effect (kSZ). By implementing various filters designed to extract the kSZ generated at the positions of the clusters, we obtain consistent constraints on the radial peculiar velocity average, root mean square (rms), and local bulk flow amplitude at different depths. For the whole cluster sample of average redshift 0.18, the measured average radial peculiar velocity with respect to the cosmic microwave background (CMB) radiation at that redshift, i.e., the kSZ monopole, amounts to $72 \pm 60$ km s$^{-1}$. This constitutes less than 1% of the relative Hubble velocity of the cluster sample with respect to our local CMB frame. While the linear $\Lambda$CDM prediction for the typical cluster radial velocity rms at $z=0.15$ is close to 230km s$^{-1}$, the upper limit imposed by \Planck\ data on the cluster subsample corresponds to 800 km s$^{-1}$ at 95% confidence level, i.e., about three times higher. \Planck\ data also set strong constraints on the local bulk flow in volumes centred on the Local Group. There is no detection of bulk flow as measured in any comoving sphere extending to the maximum redshift covered by the cluster sample. A blind search for bulk flows in this sample has an upper limit of 254 km s$^{-1}$ (95% confidence level) dominated by CMB confusion and instrumental noise, indicating that the Universe is largely homogeneous on Gpc scales. In this context, in conjunction with supernova observations, \Planck\ is able to rule out a large class of inhomogeneous void models as alternatives to dark energy or modified gravity. The \Planck\ constraints on peculiar velocities and bulk flows are thus consistent with the $\Lambda$CDM scenario.

The presence of massive neutrinos has a characteristic impact on the growth of large scale structures such as galaxy clusters. We forecast on the capability of the number count and power spectrum measured from the ongoing and future Sunyaev-Zeldovich (SZ) cluster surveys, combined with cosmic microwave background (CMB) observation to constrain the total neutrino mass $\mnu$ in a flat $\Lambda$CDM cosmology. We adopt self-calibration for the mass-observable scaling relation, and evaluate constraints for the South Pole Telescope normal and with polarization (SPT, SPTPol), Planck, and Atacama Cosmology Telescope Polarization (ACTPol) surveys. We find that a sample of $\approx1000$ clusters obtained from the Planck cluster survey plus extra information from CMB lensing extraction could tighten the current upper bound on the sum of neutrino masses to $\sigma_{\mnu}=0.17$ eV at 68% C.L. Our analysis shows that cluster number counts and power spectrum provide complementary constraints and as a result they help reducing the error bars on $\mnu$ by a factor of 4-8 when both probes are combined. We also show that the main strength of cluster measurements in constraining $\mnu$ is when good control of cluster systematics is available. When applying a weak prior on the mass-observable relations, which can be at reach in the upcoming cluster surveys, we obtain $\sigma_{\mnu}=0.48$ eV using cluster only probes and, more interestingly, $\sigma_{\mnu}=0.08$ eV using cluster + CMB which corresponds to a $S/N\approx4$ detection for $\mnu\ge0.3$ eV. We analyze and discuss the degeneracies of $\mnu$ with other parameters and investigate the sensitivity of neutrino mass constraints with various surveys specifications.

We perform a Fisher matrix analysis to forecast the capability of ongoing and future Sunyaev-Zeldovich cluster surveys in constraining the deviations from Gaussian distribution of primordial density perturbations. We use the constraining power of the cluster number counts and clustering properties to forecast limits on the $\fnl$ parameter. The primordial non-Gaussianity effects on the mass function and halo bias are considered. We adopt self-calibration for the mass-observable scaling relation, and evaluate constraints for the SPT, Planck, CCAT--like, SPTPol and ACTPol surveys. We show that the scale-dependence of halo bias induced by the local NG provides strong constraints on $\fnl$, while the results from number count are two orders of magnitude worse. When combining information from number counts and power spectrum, the \planck\ cluster catalog provides the tightest constraint with $\sigma_{\fnl}=7$ (68% C.L.) even for relatively conservative assumptions on the expected cluster yields and systematics. This value is a factor of 2 smaller than the $1\sigma$ error as measured by WMAP CMB measurements, and comparable to what expected from Planck. We find that the results are mildly sensitive to the mass threshold of the surveys, but strongly depend on the survey coverage: a full-sky survey like Planck is more favorable because it can probe longer wavelengths modes which are most sensitive to NG effects. In addition, the constraints are largely insensitive to priors on nuisance parameters as they are mainly driven by the power spectrum probe which has a mild dependence on the mass-observable relations.

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

We have looked for bulk motions of galaxy clusters in the WMAP~7 year data. We isolate the kinetic Sunyaev-Zeldovich (SZ) signal by filtering the WMAP Q, V and W band maps with multi-frequency matched filters, that utilize the spatial properties of the kinetic SZ signal to optimize detection. We try two filters: a filter that has no spectral dependence, and a filter that utilizes the spectral properties of the kinetic and thermal SZ signals to remove the thermal SZ bias. We measure the monopole and dipole spherical harmonic coefficients of the kinetic SZ signal, as well as the $\ell=2-5$ modes, at the locations of 736 ROSAT observed galaxy clusters. We find no significant power in the kinetic SZ signal at these multipoles with either filter, consistent with the $\Lambda$CDM prediction. Our limits are a factor of ~ 3 more sensitive than the claimed bulk flow detection of~\citet{2009ApJ...691.1479K}. Using simulations we estimate that in maps filtered by our matched filter with no spectral dependence there is a thermal SZ dipole that would be mistakenly measured as a bulk motion of $\sim \! 2000-4000$ km/s. For the WMAP data the signal to noise ratio obtained with the unbiased filter is almost an order of magnitude lower.

We predict the performance of the Planck satellite in determining the bulk flow through kinetic Sunyaev-Zeldovich (kSZ) measurements. As velocity tracers, we use ROSAT All-Sky Survey (RASS) clusters as well as expected cluster catalogs from the upcoming missions Planck and eRosita (All-Sky Survey: EASS). We implement a semi-analytical approach to simulate realistic Planck maps as well as Planck and eRosita cluster catalogs. We adopt an unbiased kinetic SZ filter (UF) and matched filter (MF) to maximize the cluster kSZ signal to noise ratio. We find that the use of Planck CMB maps in conjunction with the currently existing ROSAT cluster sample improves current upper limits on the bulk flow determination by a factor \sim 5 (\sim 10) when using the MF (UF). The accuracy of bulk flow measurement increases with the depth and abundance of the cluster sample: for an input bulk velocity of 500 km/s, the UF recovered velocity errors decrease from 94 km/s for RASS, to 73 km/s for Planck and to 24 km/s for EASS; while the systematic bias decreases from 44% for RASS, 5% for Planck, to 0% for EASS. The 95% upper limit for the recovered bulk flow direction \Delta {\alpha} ranges between 4 \circ and 60 \circ depending on cluster sample and adopted filter. The kSZ dipole determination is mainly limited by the effects of thermal SZ (tSZ) emission in all cases but the one of EASS clusters analyzed with the unbiased filter. This fact makes the UF preferable to the MF when analyzing Planck maps.

We presented the results of an analysis of four XMM-Newton observations of the starburst galaxy IC342 taken over a four-year span from 2001 to 2005, with an emphasis on investigating the long-term flux and spectral variability of the X-ray point sources. We detected a total of 61 X-ray sources within 35' $\times$ 30' of the galaxy down to a luminosity of (1-2)$\times$1037 erg s-1 depending on the local background. We found that 39 of the 61 detected sources showed long-term variability, in which 26 of them were classified as X-ray transients. We also found 19 sources exhibiting variations in hardness ratios or undergoing spectral transitions among observations, and were identified as spectral variables. In particular, 8 of the identified X-ray transients showed spectral variability in addition to flux variability. The diverse patterns of variability observed is indicative of a population of X-ray binaries. We used X-ray colors, flux and spectral variability, and in some cases the optical or radio counterparts to classify the detected X-ray sources into several stellar populations. We identified a total of 11 foreground stars, 1 supersoft sources (SSS), 3 quasisoft sources (QSS), and 2 supernova remnants (SNR). The identified SSS/QSS are located near or on the spiral arms, associate with young stellar populations; the 2 SNR are very close to the starburst nucleus where current star formation activities are dominated. We also discovered a spectral change in the nuclear source of IC342 for the first time by a series of X-ray spectrum analysis.

We presented the result of a high resolution (FWHM~0.5'') 12 ks Chandra HRC-I observation of the starburst galaxy IC342 taken on 2 April 2006. We identified 23 X-ray sources within the central 30' x 30' region of IC342. Our HRC-I observation resolved the historical Ultraluminous X-ray sources (ULX), X3, near the nucleus into 2 sources, namely C12 and C13, for the first time. The brighter source C12, with L(0.08-10keV)=(6.66\pm0.45)\times10^{38}ergs^-1, was spatially extended (~82 pc x 127 pc). From the astrometric registration of the X-ray image, C12 was at R.A.=03h:46m:48.43s, decl.=+68d05m47.45s, and was closer to the nucleus than C13. Thus we concluded that source was not an ULX and must instead be associated with the nucleus. The fainter source C13, with L(0.08-10keV)=(5.1\pm1.4) x 10^{37}ergs^-1 was consistent with a point source and located $6.51'' at P.A. 240 degree of C12. We also analyzed astrometrically corrected optical Hubble Space Telescope and radio Very Large Array images, a comparison with the X-ray image showed similarities in their morphologies. Regions of star formation within the central region of IC342 were clearly visible in HST H alpha image and this was the region where 3 optical star clusters and correspondingly our detected X-ray source C12 were observed. We found that a predicted X-ray emission from starburst was very close to the observed X-ray luminosity of C12, suggesting that nuclear X-ray emission in IC342 was dominated by starburst. Furthermore, we discussed the possibility of AGN in the nucleus of IC342. Although our data was not enough to give a firm existence of an AGN, it could not be discarded.

We report the identification of a recurrent ultraluminous X-ray source (ULX), a highly absorbed X-ray source (possibly a background AGN), and a young supernova remnant near the center of the starburst galaxy M82. From a series of Chandra observations taken from 1999 to 2005, we found that the transient ULX first appeared in 1999 October. The source turned off in 2000 January, but later reappeared and has been active since then. The X-ray luminosity of this source varies from below the detection level (~2.5e38 erg/s) to its active state in between ~7e39 erg/s and 1.3e40 erg/s (in the 0.5-10 keV energy band) and shows unusual spectral changes. The X-ray spectra of some Chandra observations are best fitted with an absorbed power-law model with photon index ranging from 1.3 to 1.7. These spectra are similar to those of Galactic black hole binary candidates seen in the low/hard state except that a very hard spectrum was seen in one of the observations. By comparing with near infrared images taken with the Hubble Space Telescope, the ULX is found to be located within a young star cluster. Radio imaging indicates that it is associated with a H II region. We suggest that the ULX is likely to be a > 100 solar mass intermediate-mass black hole in the low/hard state. In addition to the transient ULX, we also found a highly absorbed hard X-ray source which is likely to be an AGN and an ultraluminous X-ray emitting young supernova remnant which may be related to a 100-year old gamma-ray burst event, within 2 arcsec of the transient ULX.