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.

Quantitative modeling of the spectro-spatial distributions of energetic electrons and protons in galactic halos is needed in order to determine their interactions with the local plasma and radiation fields, and also to estimate their residual spectral densities in intracluster and intergalactic environments. We develop a semi-analytic approach for calculating the particle distributions in the halo based on a detailed diffusion model for particle propagation from acceleration sites and interactions in the galactic disk. Important overall normalization of our models is based on results from detailed modeling in the Galactic disk with the GALPROP code. This provides the essential input for determining particle distributions in the outer disk, which are used as source terms for calculating the distributions in the extensive halo for a range of values of key parameters affecting energy losses and propagation mode. Our modeling approach is applied to the two edge-on star-forming galaxies NGC 4631 and NGC 4666, for which recent mapping of radio emission in the inner halo provides the required overall normalization. We predict the levels and spatial profiles of radio, X-ray, and gamma-ray emission in the halos of these galaxies. Our quantitative modeling enables us to estimate the total calorimetric efficiencies of electrons and protons in star-forming galaxies, and to predict their residual spectral distributions in the outer halo and intergalactic space.

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.

As energetic particles diffuse out of radio and star-forming galaxies (SFGs), their intracluster density builds up to a level that could account for a substantial part or all the emission from a radio halo. We calculate the particle time-dependent, spectro-spatial distributions from a solution of a diffusion equation with radio galaxies as sources of electrons, and SFGs as sources of both electrons and protons. Whereas strong radio galaxies are typically found in the cluster (e.g., Coma) core, the fraction of SFGs increases with distance from the cluster center. Scaling particle escape rates from their sources to the reasonably well determined Galactic rates, and for realistic gas density and magnetic field spatial profiles, we find that predicted spectra and spatial profiles of radio emission from primary and secondary electrons are roughly consistent with those deduced from current measurements of the Coma halo (after subtraction of emission from the relic Coma A). Nonthermal X-ray emission is predicted to be mostly by Compton scattering of electrons from radio galaxies off the CMB, whereas $\gamma$-ray emission is primarily from the decay of neutral pions produced in interactions of protons from SFGs with protons in intracluster gas.

Hot gas in filamentary structures induces CMB aniostropy through the SZ effect. Guided by results from N-body simulations, we model the morphology and gas properties of filamentary gas and determine the power spectrum of the anisotropy. Our treatment suggests that power levels can be an appreciable fraction of the cluster contribution at multipoles $\ell\lesssim 1500$. Its spatially irregular morphology and larger characteristic angular scales can help to distinguish this SZ signature from that of clusters. In addition to intrinsic interest in this most extended SZ signal as a probe of filaments, its impact on cosmological parameter estimation should also be assessed. We find that filament `noise' can potentially bias determination of $A_s$, $n_s$, and $w$ (the normalization of the primordial power spectrum, the scalar index, and the dark energy equation of state parameter, respectively) by more than the nominal statistical uncertainty in Planck SZ survey data. More generally, when inferred from future optimal cosmic-variance-limited CMB experiments, we find that virtually all parameters will be biased by more than the nominal statistical uncertainty estimated for these next generation CMB experiments.

We report statistical results for dark matter (DM) velocity anisotropy, \beta, from a sample of some 6000 cluster-size halos (at redshift zero) identified in a \Lambda CDM hydrodynamical adaptive mesh refinement simulation performed with the Enzo code. These include profiles of \beta\ in clusters with different masses, relaxation states, and at several redshifts, modeled both as spherical and triaxial DM configurations. Specifically, although we find a large scatter in the DM velocity anisotropy profiles of different halos (across elliptical shells extending to at least ~$1.5 r_{vir}$), universal patterns are found when these are averaged over halo mass, redshift, and relaxation stage. These are characterized by a very small velocity anisotropy at the halo center, increasing outward to about 0.27 and leveling off at about $0.2 r_{vir}$. Indirect measurements of the DM velocity anisotropy fall on the upper end of the theoretically expected range. Though measured indirectly, the estimations are derived by using two different surrogate measurements - X-ray and galaxy dynamics. Current estimates of the DM velocity anisotropy are based on very small cluster sample. Increasing this sample will allow testing theoretical predictions, including the speculation that the decay of DM particles results in a large velocity boost. We also find, in accord with previous works, that halos are triaxial and likely to be more prolate when unrelaxed, whereas relaxed halos are more likely to be oblate. Our analysis does not indicate that there is significant correlation (found in some previous studies) between the radial density slope, \gamma, and \beta\ at large radii, $0.3 r_{vir} < r < r_{vir}$.

Extensive surveys of galaxy clusters motivate us to assess the likelihood of cluster-cluster lensing (CCL), namely, gravitational-lensing of a background cluster by a foreground cluster. We briefly describe the characteristics of CCLs in optical, X-ray and SZ measurements, and calculate their predicted numbers for $\Lambda$CDM parameters and a viable range of cluster mass functions and their uncertainties. The predicted number of CCLs in the strong-lensing regime varies from several ($<10$) to as high as a few dozen, depending mainly on whether lensing triaxiality bias is accounted for, through the c-M relation. A much larger number is predicted when taking into account also CCL in the weak-lensing regime. In addition to few previously suggested CCLs, we report a detection of a possible CCL in A383, where background candidate high-$z$ structures are magnified, as seen in deep Subaru observations.

The growth of structure in the universe begins at the time of radiation-matter equality, which corresponds to energy scales of $\sim 0.4 eV$. All tracers of dark matter evolution are expected to be sensitive to neutrino masses on this and smaller scales. Here we explore the possibility of using cluster number counts and power spectrum obtained from ongoing SZ surveys to constrain neutrino masses. Specifically, we forecast the capability of ongoing measurements with the PLANCK satellite and the ground-based SPT experiment, as well as measurements with the proposed EPIC satellite, to set interesting bounds on neutrino masses from their respective SZ surveys. We also consider an ACT-like CMB experiment that covers only a few hundred ${\rm deg^{2}}$ also to explore the tradeoff between the survey area and sensitivity and what effect this may have on inferred neutrino masses. We find that for such an experiment a shallow survey is preferable over a deep and low-noise scanning scheme. We also find that projected results from the PLANCK SZ survey can, in principle, be used to determine the total neutrino mass with a ($1\sigma$) uncertainty of $0.28 eV$, if the detection limit of a cluster is set at the $5\sigma$ significance level. This is twice as large as the limits expected from PLANCK CMB lensing measurements. The corresponding limits from the SPT and EPIC surveys are $\sim 0.44 eV$ and $\sim 0.12 eV$, respectively. Mapping an area of 200 deg$^{2}$, ACT measurements are predicted to attain a $1\sigma$ uncertainty of 0.61 eV; expanding the observed area to 4,000 deg$^{2}$ will decrease the uncertainty to 0.36 eV.

We identify 13 sets of multiply-lensed galaxies around MACS J0717.5+3745 ($z=0.546$), outlining a very large tangential critical curve of major axis $\sim2.8\arcmin$, filling the field of HST/ACS. The equivalent circular Einstein radius is $\theta_{e}= 55 \pm 3\arcsec$ (at an estimated source redshift of $z_{s}\sim2.5$), corresponding to $r_e\simeq 350\pm 20 kpc$ at the cluster redshift, nearly three times greater than that of A1689 ($r_e\simeq 140 kpc$ for $z_{s}=2.5$). The mass enclosed by this critical curve is very large, $7.4\pm 0.5 \times 10^{14}M_{\odot}$ and only weakly model dependent, with a relatively shallow mass profile within $r<250 kpc$, reflecting the unrelaxed appearance of this cluster. This shallow profile generates a much higher level of magnification than the well known relaxed lensing clusters of higher concentration, so that the area of sky exceeding a magnification of $>10\times$, is $\simeq 3.5\sq\arcmin$ for sources with $z\simeq 8$, making MACS J0717.5+3745 a compelling target for accessing faint objects at high redshift. We calculate that only one such cluster, with $\theta_{e}\ge 55\arcsec$, is predicted within $\sim 10^7$ Universes with $z\ge 0.55$, corresponding to a virial mass $\ge 3\times 10^{15} M_{\odot}$, for the standard $\Lambda CDM$ (WMAP5 parameters with $2\sigma$ uncertainties).

Extensive measurements of the X-ray background (XRB) yield a reasonably reliable characterisation of its basic properties. Having resolved most of the cosmic XRB into discrete sources, the levels and spectral shapes of its main components can be used to probe both the source populations and also alternative cosmological and large-scale structure models. Recent observations of clusters seem to provide evidence that clusters formed earlier and are more abundant than predicted in the standard $\Lambda$CDM model. This motivates interest in alternative models that predict enhanced power on cluster scales. We calculate predicted levels and spectra of the superposed emission from groups and clusters of galaxies in $\Lambda$CDM and in two viable alternative non-Gaussian ($\chi^2$) and early dark energy models. The predicted levels of the contribution of clusters to the XRB in the non-Gaussian models exceed the measured level at low energies and levels of the residual XRB in the 2-8 keV band; these particular models are essentially ruled out. Our work demonstrates the diagnostic value of the integrated X-ray emission from clusters, by considering also its dependences on different metallicities, gas and temperature profiles, Galactic absorption, merger scenarios, and on a non-thermal pressure component. We also show that the XRB can be used for a upper limit for the concentration parameter value.

Mapping CMB polarization is an essential ingredient of current cosmological research. Particularly challenging is the measurement of an extremely weak B-mode polarization that can potentially yield unique insight on inflation. Achieving this objective requires very precise measurements of the secondary polarization components on both large and small angular scales. Scattering of the CMB in galaxy clusters induces several polarization effects whose measurements can probe cluster properties. Perhaps more important are levels of the statistical polarization signals from the population of clusters. Power spectra of five of these polarization components are calculated and compared with the primary polarization spectra. These spectra peak at multipoles $\ell \geq 3000$, and attain levels that are unlikely to appreciably contaminate the primordial polarization signals.

We provide a quantitative assessment of the probability distribution function of the concentration parameter of galaxy clusters. We do so by using the probability distribution function of halo formation times, calculated by means of the excursion set formalism, and a formation redshift-concentration scaling derived from results of N-body simulations. Our results suggest that the observed high concentrations of several clusters are quite unlikely in the standard Lambda CDM cosmological model, but that due to various inherent uncertainties, the statistical range of the predicted distribution may be significantly wider than commonly acknowledged. In addition, the probability distribution function of the Einstein radius of A1689 is evaluated, confirming that the observed value of ~45" +/- 5" is very improbable in the currently favoured cosmological model. If, however, a variance of ~20% in the theoretically predicted value of the virial radius is assumed, than the discrepancy is much weaker. The measurement of similarly large Einstein radii in several other clusters would pose a difficulty to the standard model. If so, earlier formation of the large scale structure would be required, in accord with predictions of some quintessence models. We have indeed verified that in a viable early dark energy model large Einstein radii are predicted in as many as a few tens of high-mass clusters.

There is some observational evidence for earlier evolution of clusters of galaxies than predicted in the standard LambdaCDM model with a Gaussian primordial density fluctuation field, and a low value for the mass variance parameter sigma_8. Particularly difficult in this model is the interpretation of possible excess CMB anisotropy on cluster scales as due to the Sunyaev-Zeldovich (S-Z) effect. We have calculated S-Z power spectra in the standard model, and in two alternative models which predict higher cluster abundance - a model with non-Gaussian PDF, and an early dark energy model. As anticipated, the levels of S-Z power in the latter two models are significantly higher than in the standard model, and in good agreement with current measurements of CMB anisotropy at high multipole values. Our results provide a sufficient basis for testing the viability of the three models by futur high quality measurements of cluster abundance and the anisotropy induced by the S-Z effect.

In the standard Lambda CDM cosmological model with a Gaussian primordial density fluctuation field, the relatively low value of the mass variance parameter (sigma_8=0.74{+0.05}{-0.06}, obtained from the WMAP 3-year data) results in a reduced likelihood that the measured level of CMB anisotropy on the scales of clusters is due to the Sunyaev-Zeldovich (S-Z) effect. To assess the feasibility of producing higher levels of S-Z power, we explore two alternative models which predict higher cluster abundance. In the first model the primordial density field has a chi^2_1 distribution, whereas in the second an early dark energy component gives rise to the desired higher cluster abundance. We carry out the necessary detailed calculations of the levels of S-Z power spectra, cluster number counts, and angular 2-point correlation function of clusters, and compare (in a self-consistent way) their predicted redshift distributions. Our results provide a sufficient basis upon which the viability of the three models may be tested by future high quality measurements.

Two little explored aspects of Compton scattering of the CMB in clusters are discussed: The statistical properties of the Sunyaev-Zeldovich (S-Z) effect in the context of a non-Gaussian density fluctuation field, and the polarization patterns in a hydrodynamcially-simulated cluster. We have calculated and compared the power spectrum and cluster number counts predicted within the framework of two density fields that yield different cluster mass functions at high redshifts. This is done for the usual Press & Schechter mass function, which is based on a Gaussian density fluctuation field, and for a mass function based on a chi^2-distributed density field. We quantify the significant differences in the respective integrated S-Z observables in these two models. S-Z polarization levels and patterns strongly depend on the non-uniform distributions of intracluster gas and on peculiar and internal velocities. We have therefore calculated the patterns of two polarization components that are produced when the CMB is doubly scattered in a simulated cluster. These are found to be very different than the patterns calculated based on spherical clusters with uniform structure and simplified gas distribution.

The main statistical properties of the Sunyaev-Zeldovich (S-Z) effect - the power spectrum, cluster number counts, and angular correlation function - are calculated and compared within the framework of two density fields which differ in their predictions of the cluster mass function at high redshifts. We do so for the usual Press and Schechter mass function, which is derived on the basis of a Gaussian density fluctuation field, and for a mass function based on a chi^2 distributed density field. These three S-Z observables are found to be very significantly dependent on the choice of the mass function. The different predictions of the Gaussian and non-Gaussian density fields are probed in detail by investigating the behaviour of the three S-Z observables in terms of cluster mass and redshift. The formation time distribution of clusters is also demonstrated to be sensitive to the underlying mass function. A semi-quantitative assessment is given of its impact on the concentration parameter and the temperature of intracluster gas.

This brief, general review the Sunyaev-Zeldovich (S-Z) effect includes a description of the calculation of the total intensity change and polarization components, a discussion of the basic properties of the S-Z power spectrum and cluster number counts, and a summary of the main observational results.

X-ray background surveys indicate the likely presence of diffuse warm gas in the Local Super Cluster (LSC), in accord with expectations from hydrodynamical simulations. We assess several other manifestations of warm LSC gas; these include anisotropy in the spatial pattern of cluster Sunyaev-Zeldovich (S-Z) counts, its impact on the CMB temperature power spectrum at the lowest multipoles, and implications on measurements of the S-Z effect in and around the Virgo cluster.

The strong dependence of the mass variance parameter, sigma_8, on the adopted cluster mass-temperature relation is explored. A recently compiled X-ray cluster catalog and various mass-temperature relations are used to derive the corresponding values of sigma_8. Calculations of the power spectrum of the CMB anisotropy induced by the Sunyaev-Zeldovich effect and cluster number counts are carried out in order to assess the need for a consistent choice of the mass-temperature scaling and the parameter sigma_8. We find that the consequences of inconsistent choice of the mass-temperature relation and sigma_8 could be quite substantial, including a considerable mis-estimation of the magnitude of the power spectrum and cluster number counts. Our results can partly explain the large scatter between published estimates of the power spectrum and number counts. We also show that the range of values of the power-law index in the scaling C_l ~ sigma_8^[6-7] deduced in previous studies is likely overestimated; we obtain a more moderate dependence, C_l ~ sigma_8^[4].

The angular power spectrum of the Sunyaev-Zeldovich (SZ) effect is calculated in the $\Lambda$CDM cosmological model with the aim of investigating its detailed dependence on the cluster population, gas morphology, and gas evolution. We calculate the power spectrum for three different mass functions, compute it within the framework of isothermal and polytropic gas distributions, and explore the effect of gas evolution on the magnitude and shape of the power spectrum. We show that it is indeed possible to explain the `excess' power measured by the CBI experiment on small angular scales as originating from the SZ effect without (arbitrary) rescaling the value of $\sigma_8$, the mass variance parameter. The need for a self-consistent choice of the basic parameters characterizing the cluster population is emphasized. In particular, we stress the need for a consistent choice of the value of $\sigma_8$ extracted from fitting theoretical models for the mass function to the observed cluster X-ray temperature function, such that it agrees with the mass-temperature relation used to evaluate the cluster Comptonization parameter. Our treatment includes the explicit spectral dependence of the thermal component of the effect, which we calculate at various frequencies. We find appreciable differences between the non-relativistic and relativistic predictions for the power spectrum even for this superposed contribution from clusters at the full range of gas temperatures.