Outside galaxy clusters the competition between the inwards gravitational attraction and the outwards expansion of the Universe leads to a special radius of velocity cancellation, which is called the turn-around radius. Measurements of the turn-around radius hold promises of constraining cosmological parameters, and possibly even properties of gravity. Such a measurement is, however, complicated by the fact that the surroundings of galaxy clusters are not spherical, but instead are a complicated collection of filaments, sheets and voids. In this paper we use the results of numerically simulated universes to quantify realistic error-bars of the measurement of the turn-around radius. We find that for a $\Lambda$CDM cosmology these error-bars are typically of the order of $20\%$. We numerically simulate three different implementations of dark energy models and of a scalar dark sector interaction to address whether the turn-around radius can be used to constrain non-trivial cosmologies, and we find that only rather extreme models can be distinguished from a $\Lambda$CDM universe due to the large error-bars arising from the non-trivial cluster environments.

Dark matter dominates the properties of large cosmological structures such as galaxy clusters, and the mass profiles of the dark matter have been measured for these equilibrated structures for years using X-rays, lensing or galaxy velocities. A new method has been proposed, which should allow us to estimate a dynamical property of the dark matter, namely the velocity anisotropy. For the gas a similar velocity anisotropy is zero due to frequent collisions, however, the collisionless nature of dark matter allows it to be non-trivial. Numerical simulations have for years found non-zero and radially varying dark matter velocity anisotropies. Here we employ the method proposed by Hansen and Pifaretti (2007), and developed by Host et al. (2009) to estimate the dark matter velocity anisotropy in the bright galaxy cluster Perseus, to near 5 times the radii previously obtained. We find the dark matter velocity anisotropy to be consistent with the results of numerical simulations, however, still with large error-bars. At half the virial radius we find the velocity anisotropy to be non-zero at 1.7 standard deviations, lending support to the collisionless nature of dark matter.

The connections among galaxies, the dark matter halos where they form and the properties of the large-scale Cosmic Web still need to be completely disentangled. We use the cosmological hydrodynamical simulation TNG100 of the IllustrisTNG suite to quantify the effects played by the large-scale density field and the Cosmic Web morphology on the relation between halo mass and galaxy stellar mass. We select objects with total dynamical mass in the range $\geq 6.3\times 10^{10} h ^{-1}\, M_{\odot}$ up to a few $10^{14} h^{-1} \, M_{\odot}$ between redshift $z=4$ and redshift $z=0$. A Cosmic Web class (knot, filament, sheet, void) is assigned to each region of the volume using a density field deformation tensor-based method. We find that galaxy stellar mass strongly correlates with total dynamical mass and formation time, and more weakly with large-scale overdensity and Cosmic Web class. The latter two quantities correlate with each other, but are not entirely degenerate. Furthermore, we find that at fixed halo mass, galaxies with stellar mass lower than the median value are more likely to be found in voids and sheets, whereas galaxies with stellar mass higher than the median are more likely to be found in filaments and knots. Finally, we find that the dependence on environment is stronger for satellites than for centrals, and discuss the physical implications of these results.

We perform dark-matter-only simulations of 28 relaxed massive cluster-sized haloes for Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, to study structural differences between the models at large radii, where the impact of baryonic physics is expected to be very limited. We find that the distributions for the radial profiles of the density, ellipsoidal axis ratios, and velocity anisotropies ($\beta$) of the haloes differ considerably between the models (at the $\sim1\sigma$ level), even at $\gtrsim10\%$ of the virial radius, if the self-scattering cross section is $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$. Direct comparison with observationally inferred density profiles disfavours SIDM for $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, but in an intermediate radial range ($\sim3\%$ of the virial radius), where the impact of baryonic physics is uncertain. At this level of the cross section, we find a narrower $\beta$ distribution in SIDM, clearly skewed towards isotropic orbits, with no SIDM (90\% of CDM) haloes having $\beta>0.12$ at $7\%$ of the virial radius. We estimate that with an observational sample of $\sim30$ ($\sim10^{15}$ M$_\odot$) relaxed clusters, $\beta$ can potentially be used to put competitive constraints on SIDM, once observational uncertainties improve by a factor of a few. We study the suppression of the memory of halo assembly history in SIDM clusters. For $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, we find that this happens only in the central halo regions ($\sim1/4$ of the scale radius of the halo), and only for haloes that assembled their mass within this region earlier than a formation redshift $z_f\sim2$. Otherwise, the memory of assembly remains and is reflected in ways similar to CDM, albeit with weaker trends.

We consider the dynamics in and near galaxy clusters. Gas, dark matter and galaxies are presently falling into the clusters between approximately 1 and 5 virial radii. At very large distances, beyond 10 virial radii, all matter is following the Hubble flow, and inside the virial radius the matter particles have on average zero radial velocity. The cosmological parameters are imprinted on the infall profile of the gas, however, no method exists, which allows a measurement of it. We consider the results of two cosmological simulations (using the numerical codes RAMSES and Gadget) and find that the gas and dark matter radial velocities are very similar. We derive the relevant dynamical equations, in particular the generalized hydrostatic equilibrium equation, including both the expansion of the Universe and the cosmological background. This generalized gas equation is the main new contribution of this paper. We combine these generalized equations with the results of the numerical simulations to estimate the contribution to the measured cluster masses from the radial velocity: inside the virial radius it is negligible, and inside two virial radii the effect is below 40%, in agreement the earlier analyses for DM. We point out how the infall velocity in principle may be observable, by measuring the gas properties to distance of about two virial radii, however, this is practically not possible today.

The present day universe consists of galaxies, galaxy clusters, one-dimensional filaments and two-dimensional sheets or pancakes, all of which combine to form the cosmic web. The so called "Zeldovich pancakes", are very difficult to observe, because their overdensity is only slightly greater than the average density of the universe. Falco et al (2014) presented a method to identify Zeldovich pancakes in observational data, and these were used as a tool for estimating the mass of galaxy clusters. Here we expand and refine that observational detection method. We study two pancakes on scales of 10 Mpc, identified from spectroscopically observed galaxies near the Coma cluster, and compare with twenty numerical pancakes. We find that the observed structures have velocity dispersion about 100 km/sec, which is relatively low compared to typical groups and filaments. These velocity dispersions are consistent with those found for the numerical pancakes. We also confirm that the identified structures are in fact two-dimensional structures. Finally, we estimate the stellar to total mass of the observational pancakes to be $2 \times 10^{-4}$, within one order of magnitude, which is smaller than that of clusters of galaxies.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

The velocity distribution of dark matter near the Earth is important for an accurate analysis of the signals in terrestrial detectors. This distribution is typically extracted from numerical simulations. Here we address the possibility of deriving the velocity distribution function analytically. We derive a differential equation which is a function of radius and the radial component of the velocity. Under various assumptions this can be solved, and we compare the solution with the results from controlled numerical simulations. Our findings complement the previously derived tangential velocity distribution. We hereby demonstrate that the entire distribution function, below 0.7 v_esc, can be derived analytically for spherical and equilibrated dark matter structures.

Recent R-matrix calculations claim to produce a significant enhancement in the opacity of Fe XVII due to atomic core excitations [S. N. Nahar & A.K. Pradhan, Phys. Rev. Letters 116, 235003 (2016), arXiv:1606.02731] and assert that this enhancement is consistent with recent measurements of higher-than-predicted iron opacities [J. E. Bailey et al., Nature 517, 56 (2015)]. This comment shows that the standard opacity models which have already been directly compared with experimental data produce photon absorption cross-sections for Fe XVII that are effectively equivalent to (and in fact larger than) the new R-matrix opacities. Thus, the new R-matrix results cannot be expected to significantly impact the existing discrepancies between theory and experiment because they produce neither a "large enhancement" nor account for "missing continuum plasma opacity" relative to standard models.

We study the internal kinematics of galaxy clusters in the region beyond the sphere of virialization. Galaxies around a virialized cluster are infalling towards the cluster centre with a non-zero mean radial velocity. We develop a new formalism for describing the dynamical state of clusters, by generalizing the standard Jeans formalism with the inclusion of the peculiar infall motions of galaxies and the Hubble expansion as well as the contributions due to background cosmology. Using empirical fits to the radial profiles of density, mean radial velocity and velocity anisotropy of both a stacked cluster-mass halo and two isolated halos of a cosmological dark matter only simulation, we verify that our generalized Jeans equation correctly predicts the radial velocity dispersion out to 4 virial radii. We find that the radial velocity dispersion inferred from the standard Jeans equation is accurate up to 2 virial radii, but overestimated by \approx 20% for the stacked halo and by \approx 40% for the isolated halos, in the range \approx 2-3 virial radii. Our model depends on the logarithmic growth rate of the virial radius (function of halo mass or concentration), which we estimate in 7 different ways, and on the departure from self-similarity of the evolution of the peculiar velocity profile in virial units.

We use the exact solutions for magnetoacoustic waves in a two dimensional isothermal atmosphere with uniform inclined magnetic field to calculate the wave reflection, transmission, and conversion of slow and fast waves incident from above ($z=\infty$). This is relevant to the question of whether waves excited by flares in the solar atmosphere can penetrate the Alfv\'en/acoustic equipartition layer (which we identify as the canopy) to reach the photosphere with sufficient energy to create sunquakes. It is found that slow waves above the acoustic cutoff frequency efficiently penetrate (transmit) as acoustic (fast) waves if directed at a small attack angle to the magnetic field, with the rest converting to magnetic (slow) waves, in accord with Generalized Ray Theory. This may help explain the compact nature of seismic sources of sunquakes identified using seismic holography. The incident slow waves can also efficiently transmit at low frequency in inclined field due to the reduction in acoustic cutoff frequency (ramp effect). Incident fast (magnetic) "waves" from infinity with specified nonzero horizontal wavenumber are necessarily evanescent, but can carry energy to the equipartition level by tunnelling. It is found that this can then efficiently convert to acoustic (fast) energy that can again reach the photosphere as a travelling wave. Overall, there appear to be ample avenues for substantial compressive wave energy to penetrate the canopy and impact the photosphere.

The caustic technique for measuring mass profiles of galaxy clusters relies on the assumption of spherical symmetry. When applied to aspherical galaxy clusters, the method yields mass estimates affected by the cluster orientation. Here we employ mock redshift catalogues generated from cosmological simulations to study the effect of clusters intrinsic shape and surrounding filamentary structures on the caustic mass estimates. To this end, we develop a new method for removing perturbations from large-scale structures, modelled as the two-halo term, in a caustic analysis of stacked cluster data. We find that the cluster masses inferred from kinematical data of ~10^14 Msun clusters observed along the major axis are larger than masses from those observed along the minor axis by a factor of 1.7 within the virial radius, increasing to 1.8 within three virial radii. This discrepancy increases by 20% for the most massive clusters. In addition a smaller but still significant mass discrepancy arises when filamentary structures are present near a galaxy cluster. We find that the mean cluster mass from random sightlines is unbiased at all radii and their scatter ranges from 0.14 to 0.17 within one and three virial radii, with a 40% increase for the most massive clusters. We provide tables which estimate the caustic mass bias given observational constraints on the cluster orientation.

One of the components of the cosmic web are sheets, which are commonly referred to as Zeldovich pancakes. These are structures which have only collapsed along one dimension, as opposed to filaments or galaxies and cluster, which have collapsed along two or three dimensions. These pancakes have recently received renewed interest, since they have been shown to be useful tools for an independent method to determine galaxy cluster masses. We consider sheet-like structures resulting from cosmological simulations, which were previously used to establish the cluster mass determination method, and we show through their level of equilibration, that these structures have indeed only collapsed along the one dimension. We also extract the density profiles of these pancake, which agrees acceptably well with theoretical expectations. We derive the observable velocity distribution function (VDF) analytically by generalizing the Eddington method to one dimension, and we compare with the distribution function from the numerical simulation.

Dark matter halos are built from accretion and merging. During merging some of the dark matter particles may be ejected with velocities higher than the escape velocity. We use both N-body simulations and single-particle smooth-field simulations to demonstrate that rapid changes to the mean field potential are responsible for such ejection, and in particular that dynamical friction plays no significant role in it. Studying a range of minor mergers, we find that typically between 5-15% of the particles from the smaller of the two merging structures are ejected. We also find that the ejected particles originate essentially from the small halo, and more specifically are particles in the small halo which pass later through the region in which the merging occurs.

Experimental tests are in progress to evaluate the accuracy of the modeled iron opacity at solar interior conditions, in particular to better constrain the solar abundance problem [S. Basu and H.M. Antia, Physics Reports 457, 217 (2008)]. Here we describe measurements addressing three of the key requirements for reliable opacity experiments: control of sample conditions, independent sample condition diagnostics, and verification of sample condition uniformity. The opacity samples consist of iron/magnesium layers tamped by plastic. By changing the plastic thicknesses, we have controlled the iron plasma conditions to reach i) Te=167+/-3 eV and ne=(7.1+/-1.5)e21 e/cc, ii) Te=170+/-2 eV and ne=(2.0+/-0.2)e22 e/cc, and iii) Te=196+/-6 eV and ne=(3.8+/-0.8)e22 e/cc, which were measured by magnesium tracer K-shell spectroscopy. The opacity sample non-uniformity was directly measured by a separate experiment where Al is mixed into the side of the sample facing the radiation source and Mg into the other side. The iron condition was confirmed to be uniform within their measurement uncertainties by Al and Mg K-shell spectroscopy. The conditions are suitable for testing opacity calculations needed for modeling the solar interior, other stars, and high energy density plasmas.

We report a first search for weakly interacting massive particles (WIMPs) using the background rejection capabilities of SuperCDMS. An exposure of 577 kg-days was analyzed for WIMPs with mass < 30 GeV/c2, with the signal region blinded. Eleven events were observed after unblinding. We set an upper limit on the spin-independent WIMP-nucleon cross section of 1.2e-42 cm2 at 8 GeV/c2. This result is in tension with WIMP interpretations of recent experiments and probes new parameter space for WIMP-nucleon scattering for WIMP masses < 6 GeV/c2.

We study the evolution of the phase-space of collisionless N-body systems under repeated stirrings or perturbations, which has been shown to lead to a convergence towards a limited group of end states. This so-called attractor was previously shown to be independent of the initial system and environmental conditions. However the fundamental reason for its appearance is still unclear. It has been suggested that the origin of the attractor may be either radial infall (RI), the radial orbit instability (ROI), or energy exchange which, for instance, happens during violent relaxation. Here we examine the effects of a set of controlled perturbations, referred to as `kicks', which act in addition to the standard collisionless dynamics by allowing pre-specified instantaneous perturbations in phase-space. We first demonstrate that the attractor persists in the absence of RI and ROI by forcing the system to expand. We then consider radial velocity kicks in a rigid potential and isotropic velocity kicks, since there are no energy exchanges in these two recipes of kicks. We find that these kicks do not lead to the attractor, indicating that the energy exchange in a dynamic potential is indeed the physical mechanism responsible for the attractor.

The mass measurement of galaxy clusters is an important tool for the determination of cosmological parameters describing the matter and energy content of the Universe. However, the standard methods rely on various assumptions about the shape or the level of equilibrium of the cluster. We present a novel method of measuring cluster masses. It is complementary to most of the other methods, since it only uses kinematical information from outside the virialized cluster. Our method identifies objects, as galaxy sheets or filaments, in the cluster outer region, and infers the cluster mass by modeling how the massive cluster perturbs the motion of the structures from the Hubble flow. At the same time, this technique allows to constrain the three-dimensional orientation of the detected structures with a good accuracy. We use a cosmological numerical simulation to test the method. We then apply the method to the Coma cluster, where we find two galaxy sheets, and measure the mass of Coma to be Mvir=(9.2\pm2.4)10^{14} Msol, in good agreement with previous measurements obtained with the standard methods.

When measuring the mass profile of any given cosmological structure through internal kinematics, the distant background density is always ignored. This trick is often refereed to as the "Jeans Swindle". Without this trick a divergent term from the background density renders the mass profile undefined, however, this trick has no formal justification. We show that when one includes the expansion of the Universe in the Jeans equation, a term appears which exactly cancels the divergent term from the background. We thereby establish a formal justification for using the Jeans Swindle.

Dark matter haloes from cosmological N-body simulations typically have triaxial shapes and anisotropic velocity distributions. Recently it has been shown that the velocity anisotropy, beta, of cosmological haloes and major merger remnants depends on direction in such a way that beta is largest along the major axis and smallest along the minor axis. In this work we use a wide range of non-cosmological N-body simulations to examine halo shapes and direction-dependence of velocity anisotropy profiles. For each of our simulated haloes we define 48 cones pointing in different directions, and from the particles inside each cone we compute velocity anisotropy profiles. We find that elongated haloes can have very distinct velocity anisotropies. We group the behaviour of haloes into three different categories, that range from spherically symmetric profiles to a much more complex behaviour, where significant differences are found for beta along the major and minor axes. We encourage future studies of velocity anisotropies in haloes from cosmological simulations to calculate beta-profiles in cones, since it reveals information, which is hidden from a spherically averaged profile. Finally, we show that spherically averaged profiles often obey a linear relation between beta and the logarithmic density slope in the inner parts of haloes, but this relation is not necessarily obeyed, when properties are calculated in cones.

The Blanco Cosmology Survey (BCS) is a 60 night imaging survey of $\sim$80 deg$^2$ of the southern sky located in two fields: ($\alpha$,$\delta$)= (5 hr, $-55^{\circ}$) and (23 hr, $-55^{\circ}$). The survey was carried out between 2005 and 2008 in $griz$ bands with the Mosaic2 imager on the Blanco 4m telescope. The primary aim of the BCS survey is to provide the data required to optically confirm and measure photometric redshifts for Sunyaev-Zel'dovich effect selected galaxy clusters from the South Pole Telescope and the Atacama Cosmology Telescope. We process and calibrate the BCS data, carrying out PSF corrected model fitting photometry for all detected objects. The median 10$\sigma$ galaxy (point source) depths over the survey in $griz$ are approximately 23.3 (23.9), 23.4 (24.0), 23.0 (23.6) and 21.3 (22.1), respectively. The astrometric accuracy relative to the USNO-B survey is $\sim45$ milli-arcsec. We calibrate our absolute photometry using the stellar locus in $grizJ$ bands, and thus our absolute photometric scale derives from 2MASS which has $\sim2$% accuracy. The scatter of stars about the stellar locus indicates a systematics floor in the relative stellar photometric scatter in $griz$ that is $\sim$1.9%, $\sim$2.2%, $\sim$2.7% and$\sim$2.7%, respectively. A simple cut in the AstrOmatic star-galaxy classifier {\tt spread\_model} produces a star sample with good spatial uniformity. We use the resulting photometric catalogs to calibrate photometric redshifts for the survey and demonstrate scatter $\delta z/(1+z)=0.054$ with an outlier fraction $\eta<5$% to $z\sim1$. We highlight some selected science results to date and provide a full description of the released data products.

Dark matter haloes in cosmological N-body simulations are affected by processes such as mergers, accretion and the gravitational interaction with baryonic matter. Typically the analysis of dark matter haloes is performed in spherical or elliptical bins and the velocity distributions are often assumed to be constant within those bins. However, the velocity anisotropy, which describes differences between the radial and tangential velocity dispersion, has recently been show to have a strong dependence on direction in the triaxial halos formed in cosmological simulations. In this study we derive properties of particles in cones parallel or perpendicular to the collision axis of merger remnants. We find that the velocity anisotropy has a strong dependence on direction. The finding that the direction-dependence of the velocity anisotropy of a halo depends on the merger history, explains the existence of such trends in cosmological simulations. It also explains why a large diversity is seen in the velocity anisotropy profiles in the outer parts of high-resolution simulations of cosmological haloes.

All dark matter structures appear to follow a set of universalities, such as phase-space density or velocity anisotropy profiles, however, the origin of these universalities remains a mystery. Any equilibrated dark matter structure can be fully described by two functions, namely the radial and the tangential velocity distribution functions (VDF), and when we will understand these two then we will understand all the observed universalities. Here we demonstrate that if we know the radial VDF, then we can derive and understand the tangential VDF. This is based on simple dynamical arguments about properties of collisionless systems. We use a range of controlled numerical simulations to demonstrate the accuracy of this result. We therefore boil the question of the dark matter structural properties down to understanding the radial VDF.

Observations and numerical simulations of galaxy clusters strongly indicate that the hot intracluster x-ray emitting gas is not spherically symmetric. In many earlier studies spherical symmetry has been assumed partly because of limited data quality, however new deep observations and instrumental designs will make it possible to go beyond that assumption. Measuring the temperature and density profiles are of interest when observing the x-ray gas, however the spatial shape of the gas itself also carries very useful information. For example, it is believed that the x-ray gas shape in the inner parts of galaxy clusters is greatly affected by feedback mechanisms, cooling and rotation, and measuring this shape can therefore indirectly provide information on these mechanisms. In this paper we present a novel method to measure the three-dimensional shape of the intracluster x-ray emitting gas. We can measure the shape from the x-ray observations only, i.e. the method does not require combination with independent measurements of e.g. the cluster mass or density profile. This is possible when one uses the full spectral information contained in the observed spectra. We demonstrate the method by measuring radial dependent shapes along the line of sight for CHANDRA mock data. We find that at least 10^6 photons are required to get a 5-{\sigma} detection of shape for an x-ray gas having realistic features such as a cool core and a double powerlaw for the density profile. We illustrate how Bayes' theorem is used to find the best fitting model of the x-ray gas, an analysis that is very important in a real observational scenario where the true spatial shape is unknown. Not including a shape in the fit may propagate to a mass bias if the x-ray is used to estimate the total cluster mass. We discuss this mass bias for a class of spacial shapes.

We study the evolution of the phase-space of collisionless N-body systems under repeated stirrings or perturbations. We find convergence towards a limited solution group, in accordance with Hansen 2010, that is independent of the initial system and environmental conditions, paying particular attention to the assumed gravitational paradigm (Newtonian and MOND). We examine the effects of changes to the perturbation scheme and in doing so identify a large group of perturbations featuring radial orbit instability (ROI) which always lead to convergence. The attractor is thus found to be a robust and reproducible effect under a variety of circumstances.

When dark matter structures form and equilibrate they have to release a significant amount of energy in order to obey the virial theorem. Since dark matter is believed to be unable to radiate, this implies that some of the accreted dark matter particles must be ejected with high velocities. These ejected particles may then later hit other cosmological structures and deposit their momentum within these structures. This induces a pressure between the cosmological structures which opposes the effect of gravity and may therefore mimic a cosmological constant. We estimate the magnitude of this effect and find that it may be as large as the observed accelerated expansion. Our estimate is accurate only within a few orders of magnitude. It is therefore important to make a much more careful calculation of this redshift dependent effect, before beginning to interpret the observed accelerated expansion as a time dependent generalization of a cosmological constant.

Alfv\'en waves may be difficult to excite at the photosphere due to low ionization fraction and suffer near-total reflection at the transition region (TR). Yet they are ubiquitous in the corona and heliosphere. To overcome these difficulties, we show that they may instead be generated high in the chromosphere by conversion from reflecting fast magnetohydrodynamic waves, and that Alfv\'enic transition region reflection is greatly reduced if the fast reflection point is within a few scale heights of the TR. The influence of mode conversion on the phase of the reflected fast wave is also explored. This phase can potentially be misinterpreted as a travel speed perturbation, with implications for the practical seismic probing of active regions.

Galaxy clusters provide us with important information about the cosmology of our universe. Observations of the X-ray radiation or of the SZ effect allow us to measure the density and temperature of the hot intergalactic medium between the galaxies in a cluster, which then allow us to calculate the total mass of the galaxy cluster. However, no simple connection between the density and the temperature profiles has been identified. Here we use controlled high-resolution numerical simulations to identify a relation between the density and temperature of the gas in equilibrated galaxy clusters. We demonstrate that the temperature-density relation is a real attractor, by showing that a wide range of equilibrated structures all move towards the attractor when perturbed and subsequently allowed to relax. For structures which have undergone sufficient perturbations for this connection to hold, one can therefore extract the mass profile directly from the X-ray intensity profile.

The XMM-Newton - Blanco Cosmology Survey project (XMM-BCS) is a coordinated X-ray, optical and mid-infrared cluster survey in a field also covered by Sunyaev-Zel'dovich effect surveys by the South Pole Telescope and the Atacama Cosmology Telescope. The aim of the project is to study the cluster population in a 14 deg$^2$ field. In this work, we present a catalog of 46 X-ray selected clusters from the initial 6 deg$^2$ survey core. We describe the XMM-BCS source detection pipeline and derive physical properties of the clusters. We provide photometric redshift estimates derived from the BCS imaging data and spectroscopic redshift measurements for a low redshift subset of the clusters. We derive the cluster log N - log S relation using an approximation to the survey selection function and find it in good agreement with previous studies. We carry out an initial comparison between X-ray luminosity derived masses and masses from optical estimators from the Southern Cosmology Survey for a subset of the cluster sample. Optical masses based on cluster richness and total optical luminosity are found to be significantly higher than the X-ray values. (abridged)

The theoretical framework of cosmology is mainly defined by gravity, of which general relativity is the current model. Recent tests of general relativity within the \Lambda Cold Dark Matter (CDM) model have found a concordance between predictions and the observations of the growth rate and clustering of the cosmic web. General relativity has not hitherto been tested on cosmological scales independent of the assumptions of the \Lambda CDM model. Here we report observation of the gravitational redshift of light coming from galaxies in clusters at the 99 per cent confidence level, based upon archival data. The measurement agrees with the predictions of general relativity and its modification created to explain cosmic acceleration without the need for dark energy (f(R) theory), but is inconsistent with alternative models designed to avoid the presence of dark matter.

We report on the serendipitous discovery in the Blanco Cosmology Survey (BCS) imaging data of a z = 0.9057 galaxy that is being strongly lensed by a massive galaxy cluster at a redshift of z=0.3838. The lens (BCS J2352-5452) was discovered while examining i- and z-band images being acquired in October 2006 during a BCS observing run. Follow-up spectroscopic observations with the GMOS instrument on the Gemini South 8m telescope confirmed the lensing nature of this system. Using weak plus strong lensing, velocity dispersion, cluster richness N_200, and fitting to an NFW cluster mass density profile, we have made three independent estimates of the mass M_200 which are all very consistent with each other. The combination of the results from the three methods gives M_200 = (5.1 +- 1.3) x 10^14 M_Sun, which is fully consistent with the individual measurements. The final NFW concentration c_200 from the combined fit is c_200 = 5.4(+1.4 -1.1). We have compared our measurements of M_200 and c_200 with predictions for (a) clusters from Lambda_CDM simulations, (b) lensing selected clusters from simulations, and (c) a real sample of cluster lenses. We find that we are most compatible with the predictions for Lambda_CDM simulations for lensing clusters, and we see no evidence based on this one system for an increased concentration compared to Lambda_CDM. Finally, using the flux measured from the [OII]3727 line we have determined the star formation rate (SFR) of the source galaxy and find it to be rather modest given the assumed lens magnification.

Alfv\'en waves may be generated via mode conversion from fast magneto-acoustic waves near their reflection level in the solar atmosphere, with implications both for coronal oscillations and for active region helio-seismology. In active regions this reflection typically occurs high enough that the Alfv\'en speed $a$ greatly exceeds the sound speed $c$, well above the $a=c$ level where the fast and slow modes interact. In order to focus on the fundamental characteristics of fast/Alfv\'en conversion, stripped of unnecessary detail, it is therefore useful to freeze out the slow mode by adopting the gravitationally stratified cold MHD model $c\to0$. This provides a benchmark for fast-to-Alfv\'en mode conversion in more complex atmospheres. Assuming a uniform inclined magnetic field and an exponential Alfv\'en speed profile with density scale height $h$, the Alfv\'en conversion coefficient depends on three variables only; the dimensionless transverse-to-the-stratification wavenumber $\kappa=kh$, the magnetic field inclination from the stratification direction $\theta$, and the polarization angle $\phi$ of the wavevector relative to the plane containing the stratification and magnetic field directions. We present an extensive exploration of mode conversion in this parameter space and conclude that near-total conversion to outward-propagating Alfv\'en waves typically occurs for small $\theta$ and large $\phi$ ($80^\circ$--$90^\circ$), though it is absent entirely when $\theta$ is exactly zero (vertical field). For wavenumbers of helioseismic interest, the conversion region is broad enough to encompass the whole chromosphere.

The distribution of mass in the halos of galaxies and galaxy clusters has been probed observationally, theoretically, and in numerical simulations. Yet there is still confusion about which of several suggested parameterized models is the better representation, and whether these models are universal. We use the temperature and density profiles of the intracluster medium as measured by X-ray observations of 11 relaxed galaxy clusters to investigate mass models for the halo using a thorough Bayesian statistical analysis. We make careful comparisons between two- and three-parameter models, including the issue of a universal third parameter. We find that, of the two-parameter models, the NFW is the best representation, but we also find moderate statistical evidence that a generalized three-parameter NFW model with a freely varying inner slope is preferred, despite penalizing against the extra degree of freedom. There is a strong indication that this inner slope needs to be determined for each cluster individually, i.e. some clusters have central cores and others have steep cusps. The mass-concentration relation of our sample is in reasonable agreement with predictions based on numerical simulations.

Reducing the scatter between cluster mass and optical richness is a key goal for cluster cosmology from photometric catalogs. We consider various modifications to the red-sequence matched filter richness estimator of Rozo et al. (2009), and evaluate their impact on the scatter in X-ray luminosity at fixed richness. Most significantly, we find that deeper luminosity cuts can reduce the recovered scatter, finding that sigma_lnLX|lambda=0.63+/-0.02 for clusters with M_500c >~ 1.6e14 h_70^-1 M_sun. The corresponding scatter in mass at fixed richness is sigma_lnM|lambda ~ 0.2-0.3 depending on the richness, comparable to that for total X-ray luminosity. We find that including blue galaxies in the richness estimate increases the scatter, as does weighting galaxies by their optical luminosity. We further demonstrate that our richness estimator is very robust. Specifically, the filter employed when estimating richness can be calibrated directly from the data, without requiring a-priori calibrations of the red-sequence. We also demonstrate that the recovered richness is robust to up to 50% uncertainties in the galaxy background, as well as to the choice of photometric filter employed, so long as the filters span the 4000 A break of red-sequence galaxies. Consequently, our richness estimator can be used to compare richness estimates of different clusters, even if they do not share the same photometric data. Appendix 1 includes "easy-bake" instructions for implementing our optimal richness estimator, and we are releasing an implementation of the code that works with SDSS data, as well as an augmented maxBCG catalog with the lambda richness measured for each cluster.

We present first results of an examination of the optical properties of the galaxy populations in SZE selected galaxy clusters. Using clusters selected by the South Pole Telescope survey and deep multiband optical data from the Blanco Cosmology Survey, we measure the radial profile, the luminosity function, the blue fraction and the halo occupation number of the galaxy populations of these four clusters with redshifts ranging from 0.3 to 1. Our goal is to understand whether there are differences among the galaxy populations of these SZE selected clusters and previously studied clusters selected in the optical and the X-ray. The radial distributions of galaxies in the four systems are consistent with NFW profiles with a galaxy concentration of 3 to 6. We show that the characteristic luminosities in $griz$ bands are consistent with passively evolving populations emerging from a single burst at redshift $z=3$. The faint end power law slope of the luminosity function is found to be on average $\alpha \approx -1.2$ in griz. Halo occupation numbers (to $m^*+2$) for these systems appear to be consistent with those based on X-ray selected clusters. The blue fraction estimated to $0.36L^*$, for the three lower redshift systems, suggests an increase with redshift, although with the current sample the uncertainties are still large. Overall, this pilot study of the first four clusters provides no evidence that the galaxy populations in these systems differ significantly from those in previously studied cluster populations selected in the X-ray or the optical.

The temperature profile of hot gas in galaxies and galaxy clusters is largely determined by the depth of the total gravitational potential and thereby by the dark matter (DM) distribution. We use high-resolution hydrodynamical simulations of galaxy formation to derive a surprisingly simple relation between the gas temperature and DM properties. We show that this relation holds not just for galaxy clusters but also for equilibrated and relaxed galaxies at radii beyond the central stellar-dominated region of typically a few kpc. It is then clarified how a measurement of the temperature and density of the hot gas component can lead to an indirect measurement of the DM velocity anisotropy in galaxies. We also study the temperature relation for galaxy clusters in the presence of self-regulated, recurrent active galactic nuclei (AGN), and demonstrate that this temperature relation even holds outside the inner region of 30 kpc in clusters with an active AGN.

Cosmological simulations of dark matter structures have identified a set of universal profiles, and similar characteristics have been seen in non-cosmological simulations. It has therefore been speculated whether these profiles of collisionless systems relate to accretion and merger history, or if there is an attractor for the dark matter systems. Here we identify such a 1-dimensional attractor in the 3-dimensional space spanned by the 2 radial slopes of the density and velocity dispersion, and the velocity anisotropy. This attractor effectively removes one degree of freedom from the Jeans equation. It also allows us to speculate on a new fluid interpretation for the Jeans equation, with an effective polytropic index for the dark matter particles between 1/2 and 3/4. If this attractor solution holds for other collisionless structures, then it may hold the key to break the mass-anisotropy degeneracy, which presently prevents us from measuring the mass profiles in dwarf galaxies uniquely.

A signal from decaying dark matter (DM) can be unambiguously distinguished from spectral features of astrophysical or instrumental origin by studying its spatial distribution. We demonstrate this approach by examining the recent claim of 0912.0552 regarding the possible DM origin of the 2.5 keV line in Chandra observations of the Milky Way satellite known as Willman 1. Our conservative strategy is to adopt a relatively large dark mass for Willman 1 and relatively small dark masses for the comparison objects. We analyze archival observations by XMM-Newton of M31 and Fornax dwarf spheroidal galaxy (dSph) and Chandra observations of Sculptor dSph. By performing a conservative analysis of X-ray spectra, we show the absence of a DM decay line with parameters consistent with those of 0912.0552. For M31, the observations of the regions between 10 and 20 kpc from the center, where the uncertainties in the DM distribution are minimal, make a strong exclusion at the level above 10sigma. The minimal estimate for the amount of DM in the central 40 kpc of M31 is provided by the model of 0912.4133, assuming the stellar disk's mass to light ratio ~8 and almost constant DM density within a core of 28 kpc. Even in this case one gets an exclusion at 5.7sigma from central region of M31 whereas modeling all processed data from M31 and Fornax produces more than 14sigma exclusion. Therefore, despite possible systematic uncertainties, we exclude the possibility that the spectral feature at ~2.5 keV found in 0912.0552 is a DM decay line. We conclude, however, that the search for DM decay line, although demanding prolonged observations of well-studied dSphs, M31 outskirts and other similar objects, is rather promising, as the nature of a possible signal can be checked. An (expected) non-observation of a DM decay signal in the planned observations of Willman 1 should not discourage further dedicated observations.

Using a combination of deep 574ks Chandra data, XMM-Newton high-resolution spectra, and optical Halpha+NII images, we study the nature and spatial distribution of the multiphase plasma in M87. Our results provide direct observational evidence of `radio mode' AGN feedback in action, stripping the central galaxy of its lowest entropy gas and preventing star-formation. This low entropy gas was entrained with and uplifted by the buoyantly rising relativistic plasma, forming long "arms". These arms are likely oriented within 15-30 degrees of our line-of-sight. The mass of the uplifted gas in the arms is comparable to the gas mass in the approximately spherically symmetric 3.8 kpc core, demonstrating that the AGN has a profound effect on its immediate surroundings. The coolest X-ray emitting gas in M87 has a temperature of ~0.5 keV and is spatially coincident with Halpha+NII nebulae, forming a multiphase medium where the cooler gas phases are arranged in magnetized filaments. We place strong upper limits of 0.06 Msun/yr on the amount of plasma cooling radiatively from 0.5 keV and show that a uniform, volume-averaged heating mechanism could not be preventing the cool gas from further cooling. All of the bright Halpha filaments appear in the downstream region of the <3 Myr old shock front, at smaller radii than ~0.6'. We suggest that shocks induce shearing around the filaments, thereby promoting mixing of the cold gas with the ambient hot ICM via instabilities. By bringing hot thermal particles into contact with the cool, line-emitting gas, mixing can supply the power and ionizing particles needed to explain the observed optical spectra. Mixing of the coolest X-ray emitting plasma with the cold optical line emitting filamentary gas promotes efficient conduction between the two phases, allowing non-radiative cooling which could explain the lack of X-ray gas with temperatures under 0.5 keV.

We present redshifts and optical richness properties of 21 galaxy clusters uniformly selected by their Sunyaev-Zel'dovich signature. These clusters, plus an additional, unconfirmed candidate, were detected in a 178 square-degree area surveyed by the South Pole Telescope in 2008. Using griz imaging from the Blanco Cosmology Survey and from pointed Magellan telescope observations, as well as spectroscopy using Magellan facilities, we confirm the existence of clustered red-sequence galaxies, report red-sequence photometric redshifts, present spectroscopic redshifts for a subsample, and derive R_200 radii and M_200 masses from optical richness. The clusters span redshifts from 0.15 to greater than 1, with a median redshift of 0.74; three clusters are estimated to be at z > 1. Redshifts inferred from mean red-sequence colors exhibit 2% RMS scatter in sigma_z/(1+z) with respect to the spectroscopic subsample for z < 1. We show that M_200 cluster masses derived from optical richness correlate with masses derived from South Pole Telescope data and agree with previously derived scaling relations to within the uncertainties. Optical and infrared imaging is an efficient means of cluster identification and redshift estimation in large Sunyaev-Zel'dovich surveys, and exploiting the same data for richness measurements, as we have done, will be useful for constraining cluster masses and radii for large samples in cosmological analysis.

We propose an X-ray mission called Xenia to search for decaying superweakly interacting Dark Matter particles (super-WIMP) with a mass in the keV range. The mission and its observation plan are capable of providing a major break through in our understanding of the nature of Dark Matter (DM). It will confirm, or reject, predictions of a number of particle physics models by increasing the sensitivity of the search for decaying DM by about two orders of magnitude through a wide-field imaging X-ray spectrometer in combination with a dedicated observation program. The proposed mission will provide unique limits on the mixing angle and mass of neutral leptons, right handed partners of neutrinos, which are important Dark Matter candidates. The existence of these particles is strongly motivated by observed neutrino flavor oscillations and the problem of baryon asymmetry of the Universe. In super-WIMP models, the details of the formation of the cosmic web are different from those of LambdaCDM. The proposed mission will, in addition to the search for decaying Dark Matter, provide crucial insight into the nature of DM by studying the structure of the "cosmic web". This will be done by searching for missing baryons in emission, and by using gamma-ray bursts as backlight to observe the warm-hot intergalactic media in absorption.

We demonstrate that all properties of the hot X-ray emitting gas in galaxy clusters are completely determined by the underlying dark matter (DM) structure. Apart from the standard conditions of spherical symmetry and hydrostatic equilibrium for the gas, our proof is based on the Jeans equation for the DM and two simple relations which have recently emerged from numerical simulations: the equality of the gas and DM temperatures, and the almost linear relation between the DM velocity anisotropy profile and its density slope. For DM distributions described by the NFW or the Sersic profiles, the resulting gas density profile, the gas-to-total-mass ratio profile, and the entropy profile are all in good agreement with X-ray observations. All these profiles are derived using zero free parameters. Our result allows us to predict the X-ray luminosity profile of a cluster in terms of its DM content alone. As a consequence, a new strategy becomes available to constrain the DM morphology in galaxy clusters from X-ray observations. Our results can also be used as a practical tool for creating initial conditions for realistic cosmological structures to be used in numerical simulations.

Context. The sterile neutrino is an excellent dark matter candidate, which can be searched for in a wide range of astrophysical sites. It has previously been shown that the optimal search strategy is to consider dwarf galaxies belonging to the Milky Way. Aims. We search for line emission from decaying dark matter. Methods. We analyse publicly available Chandra X-ray observations of the dwarf galaxy Draco. Results. The Draco and blank sky (background) spectra are nearly identical in shape, which allows us to conclude: i) Dwarf spheroidals are ideal for studying dark matter X-ray emission since the baryonic noise is impressively low, ii) there is very little room for line emission, which leads to constraints in the mass-mixing angle parameter space of the sterile neutrino. We compare the standard flux derivation method to a very conservative rebinning approach. The resulting constraints are strongly dependent on the chosen method.

Cosmological simulations of dark matter structures have shown that the equilibrated dark matter structures have a fairly small angular momentum. It appears from these N-body simulations that the radial profile of the angular momentum has an almost universal behavior, even if the different dark matter structures have experienced very different formation and merger histories. We suggest a perturbed Jeans equation, which includes a rotational term. This is done under a reasonable assumed form of the change in the distribution function. By conjecturing that the (new) subdominant rotation term must be proportional to the (old) dominant mass term, we find a clear connection, which is in rather good agreement with the results of recent high resolution simulations. We also present a new connection between the radial profiles of the angular momentum and the velocity anisotropy, which is also in fair agreement with numerical findings. Finally we show how the spin parameter $\lambda$ increases as a function of radius.

We use the abundance and weak lensing mass measurements of the SDSS maxBCG cluster catalog to simultaneously constrain cosmology and the richness--mass relation of the clusters. Assuming a flat \LambdaCDM cosmology, we find \sigma_8(\Omega_m/0.25)^{0.41} = 0.832\pm 0.033 after marginalization over all systematics. In common with previous studies, our error budget is dominated by systematic uncertainties, the primary two being the absolute mass scale of the weak lensing masses of the maxBCG clusters, and uncertainty in the scatter of the richness--mass relation. Our constraints are fully consistent with the WMAP five-year data, and in a joint analysis we find \sigma_8=0.807\pm 0.020 and \Omega_m=0.265\pm 0.016, an improvement of nearly a factor of two relative to WMAP5 alone. Our results are also in excellent agreement with and comparable in precision to the latest cosmological constraints from X-ray cluster abundances. The remarkable consistency among these results demonstrates that cluster abundance constraints are not only tight but also robust, and highlight the power of optically-selected cluster samples to produce precision constraints on cosmological parameters.

Generalized Ray Theory (GRT) provides a simple description of MHD mode transmission and conversion between magnetoacoustic fast and slow waves and is directly applicable to solar active regions. Here it is tested in a simple two-dimensional, isothermal, gravitationally-stratified model with inclined magnetic field using previously published exact solutions and found to perform very well.

The density profile of simulated dark matter structures is fairly well-established, and several explanations for its characteristics have been put forward. In contrast, the radial variation of the velocity anisotropy has still not been explained. We suggest a very simple origin, based on the shapes of the velocity distributions functions, which are shown to differ between the radial and tangential directions. This allows us to derive a radial variation of the anisotropy profile which is in good agreement with both simulations and observations. One of the consequences of this suggestion is that the velocity anisotropy is entirely determined once the density profile is known. We demonstrate how this explains the origin of the \gamma-\beta relation, which is the connection between the slope of the density profile and the velocity anisotropy. These findings provide us with a powerful tool, which allows us to close the Jeans equations.

The radial density profile of dark matter structures has been observed to have an almost universal behaviour in numerical simulations, however, the physical reason for this behaviour remains unclear. It has previously been shown that if the pseudo phase-space density, rho/sigma_d^epsilon, is a beautifully simple power-law in radius, with the "golden values" epsilon=3 and d=r (i.e., the phase-space density is only dependent on the radial component of the velocity dispersion), then one can analytically derive the radial variation of the mass profile, dispersion profile etc. That would imply, if correct, that we just have to explain why rho/sigma^3_r ~r^{-alpha}, and then we would understand everything about equilibrated DM structures. Here we use a set of simulated galaxies and clusters of galaxies to demonstrate that there are no such golden values, but that each structure instead has its own set of values. Considering the same structure at different redshifts shows no evolution of the phase-space parameters towards fixed points. There is also no clear connection between the halo virialized mass and these parameters. This implies that we still do not understand the origin of the profiles of dark matter structures.

Minimizing the scatter between cluster mass and accessible observables is an important goal for cluster cosmology. In this work, we introduce a new matched filter richness estimator, and test its performance using the maxBCG cluster catalog. Our new estimator significantly reduces the variance in the L_X-richness relation, from \sigma_{\ln L_X}^2=(0.86\pm0.02)^2 to \sigma_{\ln L_X}^2=(0.69\pm0.02)^2. Relative to the maxBCG richness estimate, it also removes the strong redshift dependence of the richness scaling relations, and is significantly more robust to photometric and redshift errors. These improvements are largely due to our more sophisticated treatment of galaxy color data. We also demonstrate the scatter in the L_X-richness relation depends on the aperture used to estimate cluster richness, and introduce a novel approach for optimizing said aperture which can be easily generalized to other mass tracers.

We measure the logarithmic scatter in mass at fixed richness for clusters in the maxBCG cluster catalog, an optically selected cluster sample drawn from SDSS imaging data. Our measurement is achieved by demanding consistency between available weak lensing and X-ray measurements of the maxBCG clusters, and the X-ray luminosity--mass relation inferred from the 400d X-ray cluster survey, a flux limited X-ray cluster survey. We find \sigma_{\ln M|N_{200}}=0.45^{+0.20}_{-0.18} (95% CL) at N_{200} ~ 40, where N_{200} is the number of red sequence galaxies in a cluster. As a byproduct of our analysis, we also obtain a constraint on the correlation coefficient between \ln Lx and \ln M at fixed richness, which is best expressed as a lower limit, r_{L,M|N} >= 0.85 (95% CL). This is the first observational constraint placed on a correlation coefficient involving two different cluster mass tracers. We use our results to produce a state of the art estimate of the halo mass function at z=0.23 -- the median redshift of the maxBCG cluster sample -- and find that it is consistent with the WMAP5 cosmology. Both the mass function data and its covariance matrix are presented.

The internal dynamics of a dark matter structure may have the remarkable property that the local temperature in the structure depends on direction. This is parametrized by the velocity anisotropy beta which must be zero for relaxed collisional structures, but has been shown to be non-zero in numerical simulations of dark matter structures. Here we present a method to infer the radial profile of the velocity anisotropy of the dark matter halo in a galaxy cluster from X-ray observables of the intracluster gas. This non-parametric method is based on a universal relation between the dark matter temperature and the gas temperature which is confirmed through numerical simulations. We apply this method to observational data and we find that beta is significantly different from zero at intermediate radii. Thus we find a strong indication that dark matter is effectively collisionless on the dynamical time-scale of clusters, which implies an upper limit on the self-interaction cross-section per unit mass sigma/m < 1 cm2/g. Our results may provide an independent way to determine the stellar mass density in the central regions of a relaxed cluster, as well as a test of whether a cluster is in fact relaxed.

We present measurements of the excess mass-to-light ratio measured aroundMaxBCG galaxy clusters observed in the SDSS. This red sequence cluster sample includes objects from small groups with masses ranging from ~5x10^{12} to ~10^{15} M_{sun}/h. Using cross-correlation weak lensing, we measure the excess mass density profile above the universal mean \Delta \rho(r) = \rho(r) - \bar{\rho} for clusters in bins of richness and optical luminosity. We also measure the excess luminosity density \Delta l(r) = l(r) - \bar{l} measured in the z=0.25 i-band. For both mass and light, we de-project the profiles to produce 3D mass and light profiles over scales from 25 kpc/ to 22 Mpc/h. From these profiles we calculate the cumulative excess mass M(r) and excess light L(r) as a function of separation from the BCG. On small scales, where \rho(r) >> \bar{\rho}, the integrated mass-to-light profile may be interpreted as the cluster mass-to-light ratio. We find the M/L_{200}, the mass-to-light ratio within r_{200}, scales with cluster mass as a power law with index 0.33+/-0.02. On large scales, where \rho(r) ~ \bar{\rho}, the M/L approaches an asymptotic value independent of cluster richness. For small groups, the mean M/L_{200} is much smaller than the asymptotic value, while for large clusters it is consistent with the asymptotic value. This asymptotic value should be proportional to the mean mass-to-light ratio of the universe <M/L>. We find <M/L>/b^2_{ml} = 362+/-54 h (statistical). There is additional uncertainty in the overall calibration at the ~10% level. The parameter b_{ml} is primarily a function of the bias of the L <~ L_* galaxies used as light tracers, and should be of order unity. Multiplying by the luminosity density in the same bandpass we find \Omega_m/b^2_{ml} = 0.02+/-0.03, independent of the Hubble parameter.

We consider direct dark matter detection rates and investigate the difference between a standard Maxwell-Boltzmann velocity distribution and a "realistic" distribution like the ones extracted from numerical N-body simulations. Sizable differences are observed when such results are compared to the standard Maxwell-Boltzmann distribution. For a light target both the total rate and the annual modulation are reduced by ~25%. For a heavy target the total rate is virtually unchanged, whereas the annual modulation is modified by up to 50%, depending on the WIMP mass and detector energy threshold. We also consider the effect of a possible velocity anisotropy, and the effect is found to be largest for a light target For the realistic velocity distribution the anisotropy may reduce the annual modulation, in contrast to the Maxwell-Boltzmann case.

The Universe contains approximately 6 times more dark matter than normal baryonic matter, and a directly observed fundamental difference between dark matter and baryons would both be significant for our understanding of dark matter structures and provide us with information about the basic characteristics of the dark matter particle. We discuss one distinctive feature of dark matter structures in equilibrium, namely the property that a local dark matter temperature may depend on direction. This is in stark contrast to baryonic gases. We used X-ray observations of two nearby, relaxed galaxy clusters, under the assumptions of hydrostatic equilibrium and identical dark matter and gas temperatures in the outer cluster region, to measure this dark matter temperature anisotropy beta_dm, with non-parametric Monte Carlo methods. We find that beta_dm is greater than the value predicted for baryonic gases, beta_gas=0, at more than 3 sigma confidence. The observed value of the temperature anisotropy is in fair agreement with the results of cosmological N-body simulations and shows that the equilibration of the dark matter particles is not governed by local point-like interactions in contrast to baryonic gases.

Imaging data from the Sloan Digital Sky Survey are used to characterize the population of galaxies in groups and clusters detected with the MaxBCG algorithm. We investigate the dependence of Brightest Cluster Galaxy (BCG) luminosity, and the distributions of satellite galaxy luminosity and satellite color, on cluster properties over the redshift range 0.1 < z < 0.3. The size of the dataset allows us to make measurements in many bins of cluster richness, radius and redshift. We find that, within r_200 of clusters with mass above 3e13 h-1 M_sun, the luminosity function of both red and blue satellites is only weakly dependent on richness. We further find that the shape of the satellite luminosity function does not depend on cluster-centric distance for magnitudes brighter than ^{0.25}M_i - 5log(h) < -19. However, the mix of faint red and blue galaxies changes dramatically. The satellite red fraction is dependent on cluster-centric distance, galaxy luminosity and cluster mass, and also increases by ~5% between redshifts 0.28 and 0.2, independent of richness. We find that BCG luminosity is tightly correlated with cluster richness, scaling as L_{BCG} ~ M_{200}^{0.3}, and has a Gaussian distribution at fixed richness, with sigma_{log L} ~ 0.17 for massive clusters. The ratios of BCG luminosity to total cluster luminosity and characteristic satellite luminosity scale strongly with cluster richness: in richer systems, BCGs contribute a smaller fraction of the total light, but are brighter compared to typical satellites. This study demonstrates the power of cross-correlation techniques for measuring galaxy populations in purely photometric data.

We show that the two basic assumptions of the model recently proposed by Manrique and coworkers for the universal density profile of cold dark matter (CDM) halos, namely that these objects grow inside out in periods of smooth accretion and that their mass profile and its radial derivatives are all continuous functions, are both well understood in terms of the very nature of CDM. Those two assumptions allow one to derive the typical density profile of halos of a given mass from the accretion rate characteristic of the particular cosmology. This profile was shown by Manrique and coworkers to recover the results of numerical simulations. In the present paper, we investigate its behavior beyond the ranges covered by present-day N-body simulations. We find that the central asymptotic logarithmic slope depends crucially on the shape of the power spectrum of density perturbations: it is equal to a constant negative value for power-law spectra and has central cores for the standard CDM power spectrum. The predicted density profile in the CDM case is well fitted by the 3D S\'ersic profile over at least 10 decades in halo mass. The values of the S\'ersic parameters depend on the mass of the structure considered. A practical procedure is provided that allows one to infer the typical values of the best NFW or S\'ersic fitting law parameters for halos of any mass and redshift in any given standard CDM cosmology.

Numerous ongoing experiments aim at detecting WIMP dark matter particles from the galactic halo directly through WIMP-nucleon interactions. Once such a detection is established a confirmation of the galactic origin of the signal is needed. This requires a direction-sensitive detector. We show that such a detector can measure the velocity anisotropy beta of the galactic halo. Cosmological N-body simulations predict the dark matter anisotropy to be nonzero, beta~0.2. Baryonic matter has beta=0 and therefore a detection of a nonzero beta would be strong proof of the fundamental difference between dark and baryonic matter. We estimate the sensitivity for various detector configurations using Monte Carlo methods and we show that the strongest signal is found in the relatively few high recoil energy events. Measuring beta to the precision of ~0.03 will require detecting more than 10^4 WIMP events with nuclear recoil energies greater than 100 keV for a WIMP mass of 100 GeV and a 32S target. This number corresponds to ~10^6 events at all energies. We discuss variations with respect to input parameters and we show that our method is robust to the presence of backgrounds and discuss the possible improved sensitivity for an energy-sensitive detector.

We consider the possibility of constraining the lifetime of radiatively decaying dark matter particles in clusters of galaxies inspired from generic axions of the Kaluza-Klein type. Such axions have been invoked as a possible explanation for the coronal X-ray emission from the Sun. These axions, or similar particles, can be produced inside stars and some of them remain confined by the deep gravitational potential of clusters of galaxies. Specifically, we consider regions within merging galaxy clusters (Abell 520 and the "Bullet Cluster"), where gravitational lensing observations have identified massive, but baryon poor, structures. From an analysis of X-ray observations of these mass concentrations, and the expected photon spectrum of decaying solar KK-axions, we derive lower limits to the lifetime of such axions of tau approx. 10^23 sec. However, if KK-axions constitute less than a few percent of the dark matter mass, this lifetime constraint is similar to that derived from solar KK-axions.

It has recently been shown that relaxed spherically symmetric dark matter halos develop from the inside out, by permanently adapting their inner structure to the boundary conditions imposed by the current accretion rate. Such a growth allows one to infer the typical density profiles of halos. Here we follow the same approach to infer the typical spherically averaged profiles of the main structural and kinematic properties of triaxial, anisotropic, rotating halos. Specifically, we derive their density, spatial velocity dispersion, phase-space density, anisotropy and specific angular momentum profiles. The results obtained are in agreement with available data on these profiles from \nbody simulations.

N-body simulations sample their initial conditions on an initial particle distribution, which for cosmological simulations is usually a glass or grid, whilst a Poisson distribution is used for galaxy models, spherical collapse etc. These pre-initial conditions have inherent correlations, noise due to discreteness and preferential alignments, whilst the glass distribution is poorly defined and computationally expensive to construct. We present a novel particle distribution which can be useful as a pre-initial condition for N-body simulations, using a simple construction based on a ``quaquaversal'' tiling of space. This distribution has little preferred orientation (i.e. is statistically isotropic), has a rapidly vanishing large scale power-spectrum (P(k) ~ k^4), and is trivial to create. It should be particularly useful for warm dark matter and cold collapse simulations.

Gravitational lensing observations of galaxy clusters have identified dark matter ``blobs'' with remarkably low baryonic content. We use such a system to probe the particle nature of dark matter with X-ray observations. We also study high resolution X-ray grating spectra of a cluster of galaxies. From these grating spectra we improve the conservative constraints on a particular dark matter candidate, the sterile neutrino, by more than one order of magnitude. Based on these conservative constraints obtained from Cosmic X-ray observations alone, the low mass (m_s < 10keV) and low mixing angle (sin^2(2\theta) 10^{-6}) sterile neutrino is still a viable dark matter candidate.

Extended gas haloes around galaxies are a ubiquitous prediction of galaxy formation scenarios. However, the density profiles of this hot halo gas is virtually unknown, although various profiles have been suggested on theoretical grounds. In order to quantitatively address the gas profile, we compare galaxies from direct cosmological simulations with analytical solutions of the underlying gas equations. We find remarkable agreement between simulations and theoretical predictions. We present an expression for this gas profile with a non-trivial dependence on the total mass profile. This expression is useful when setting up equilibrium galaxy models for numerical experiments.

We consider the possibility of constraining decaying dark matter by looking out through the Milky Way halo. Specifically we use Chandra blank sky observations to constrain the parameter space of sterile neutrinos. We find that a broad band in parameter space is still open, leaving the sterile neutrino as an excellent dark matter candidate.

One can solve the Jeans equation analytically for equilibrated dark matter structures, once given two pieces of input from numerical simulations. These inputs are 1) a connection between phase-space density and radius, and 2) a connection between velocity anisotropy and density slope, the \alpha-\beta relation. The first (phase-space density v.s. radius) has already been analysed through several different simulations, however the second (\alpha-\beta relation) has not been quantified yet. We perform a large set of numerical experiments in order to quantify the slope and zero-point of the \alpha-\beta relation. We find strong indication that the relation is indeed an attractor. When combined with the assumption of phase-space being a power-law in radius, this allows us to conclude that equilibrated dark matter structures indeed have zero central velocity anisotropy \beta_0 = 0, central density slope of \alpha_0 = -0.8, and outer anisotropy of \beta_\infty = 0.5.

One can solve the Jeans equation analytically for equilibrated dark matter structures, once given two pieces of input from numerical simulations. These inputs are 1) a connection between phase-space density and radius, and 2) a connection between velocity anisotropy and density slope, the \alpha-\beta relation. The first (phase-space density v.s. radius) has been analysed through several different simulations, however the second (\alpha-\beta relation) has not been quantified yet. We perform a large set of numerical experiments in order to quantify the slope and zero-point of the \alpha-\beta relation. When combined with the assumption of phase-space being a power-law in radius this allows us to conclude that equilibrated dark matter structures indeed have zero central velocity anisotropy, central density slope of \alpha_0 = -0.8, and outer anisotropy of approximately \beta_\infinity = 0.5.

We identify a universal relation between the radial density slope \alpha (r) and the velocity anisotropy \beta (r) for equilibrated structures. This relation holds for a variety of systems, including disk galaxy mergers, spherical collapses, cold dark matter (CDM) halos both with and without cooling. We argue that the shape of the relation is reasonable from fundamental principles when the dark matter or stars are assumed to obey Tsallis statistics, and in that case we fit the \alpha - \beta relation with just one free parameter. One can use this result to close the Jeans equations, for example to construct mass models of elliptical galaxies from observational data or to tune dark matter direct detection experiments. We also predict the asymptotic central slope and anisotropy of CDM halos to be approximately -1 and 0.

Several general trends have been identified for equilibrated, self-gravitating collisionless systems, such as density or anisotropy profiles. These are integrated quantities which naturally depend on the underlying velocity distribution function (VDF) of the system. We study this VDF through a set of numerical simulations, which allow us to extract both the radial and the tangential VDF. We find that the shape of the VDF is universal, in the sense that it depends only on two things namely the dispersion (radial or tangential) and the local slope of the density. Both the radial and the tangential VDF's are universal for a collection of simulations, including controlled collisions with very different initial conditions, radial infall simulation, and structures formed in cosmological simulations.

Imaging data from the Sloan Digital Sky Survey is used to measure the empirical size-richness relation for a large sample of galaxy clusters. Using population subtraction methods, we determine the radius at which the cluster galaxy number density is 200/Omega_m times the mean galaxy density, without assuming a model for the radial distribution of galaxies in clusters. If these galaxies are unbiased on Mpc scales, this galaxy-density-based R_200 reflects the characteristic radii of clusters. We measure the scaling of this characteristic radius with richness over an order of magnitude in cluster richness, from rich clusters to poor groups. We use this information to examine the radial profiles of galaxies in clusters as a function of cluster richness, finding that the concentration of the galaxy distribution decreases with richness and is systematically lower than the concentrations measured for dark matter profiles in N-body simulations. Using these scaled radii, we investigate the behavior of the cluster luminosity function, and find that it is well matched by a Schechter function for galaxies brighter than M_r = -18 only after the central galaxy has been removed. We find that the luminosity function varies with richness and with distance from the cluster center, underscoring the importance of using an aperture that scales with cluster mass to compare physically equivalent regions of these different systems. We note that the lowest richness systems in our catalog have properties consistent with those expected of the earliest-forming halos; our cluster-finding algorithm, in addition to reliably finding clusters, may be efficient at finding fossil groups.

Neutrinos may possibly violate the spin-statistics theorem, and hence obey Bose statistics or mixed statistics despite having spin half. We find the generalized equilibrium distribution function of neutrinos which depends on a single fermi-bose parameter, \kappa, and interpolates continuously between the bosonic and fermionic distributions when \kappa changes from -1 to +1. We consider modification of the Big Bang Nucleosynthesis (BBN) in the presence of bosonic or partly bosonic neutrinos. For pure bosonic neutrinos the abundances change (in comparison with the usual Fermi-Dirac case) by -3.2% for 4He (which is equivalent to a decrease of the effective number of neutrinos by \Delta N_\nu = - 0.6), +2.6% for 2H and -7% for 7Li. These changes provide a better fit to the BBN data. Future BBN studies will be able to constrain the fermi-bose parameter to \kappa > 0.5, if no deviation from fermionic nature of neutrinos is found. We also evaluate the sensitivity of future CMB and LSS observations to the fermi-bose parameter.

(Abriged version) The Sunyaev-Zel'dovich (SZ) effect of galaxy clusters is a tool to measure three quantities: Compton parameter, electron temperature, and cluster peculiar velocity. However, a major problem is non-removed contamination by astrophysical sources that emit in the SZ frequencies. This includes interstellar dust emission, infra-red (IR) galaxies, and radio sources in addition to primary Cosmic Microwave Background (CMB) anisotropies. The three former contaminations induce systematic shifts in the three SZ parameters. In this study, we carefully estimated, both for a large beam experiment (namely Planck Surveyor) and a small beam experiment (ACT-like), the systematic errors that result if a fraction of the expected levels of emission from dust, IR galaxies, and radio sources remains non-removed. ...

We propose a novel component to the understanding of the temperature structure of galaxy clusters which does not rely on any heating or cooling mechanism. The new ingredient is the use of non-extensive thermo-statistics which is based on the natural generalization of entropy for systems with long-range interactions. Such interactions include gravity and attraction or repulsion due to charges. We explain that there is growing theoretical indications for the need of this generalization for large cosmological structures. The observed pseudo temperature is generally different from the true thermodynamic temperature, and we clarify the connection between the two. We explain that this distinction is most important in the central part of the cluster where the density profile is most shallow. We show that the observed pseudo temperature may differ up to a factor 2/5 from the true thermodynamic temperature, either larger or smaller. In general the M-T and L-T relations will be affected, and the central DM slope derived through hydrostatic equilibrium may be either more shallow or steeper. We show how the true temperature can be extracted correctly either from the spectrum or from the shape of the Doppler broadening of spectral lines.

We present an analytical and numerical study of the velocity distribution function of self gravitating collisionless particles, which include dark matter and star clusters. We show that the velocity distribution derived through the Eddington's formula is identical to the analytical one derived directly from the generalized entropy of non-extensive statistical mechanics. This implies that self gravitating collisionless structures are to be described by non-extensive thermo-statistics. We identify a connection between the density slope of dark matter structures, \gamma, from \rho ~ r^{-\gamma}, and the entropic index, q, from the generalized entropy, S_q. Our numerical result confirms the analytical findings of earlier studies and clarifies which is the correct connection between the density slope and the entropic index. We use this result to conclude that from a fundamental statistical mechanics point of view the central density slope of self gravitating collisionless dark matter structures is not constrained, and even cored dark matter structures are allowed with \gamma = 0. We find that the outer density slope is bounded by \gamma= 10/3.

In the near future the Planck satellite will gather impressive information about the anisotropies of the cosmic microwave background and about the galaxy clusters which perturb that signal. We will here review the ability of Planck to extract information about these galaxy clusters through the Sunyaev-Zeldovich imprint. We will conclude that Planck will provide a catalogue of galaxy clusters, which will be very useful for future targeted observations. We will explain why Planck will not be very good in extracting detailed information about individual clusters, except for the dominating Compton parameter, $y$, which will be measured to a few percent for individual clusters. In this last point I am being rather conservative, but that will leave space for pleasant surprises when the Planck data will be analysed.

We make a simple analytical study of radial profiles of dark matter structures, with special attention to the question of the central radial density profile. We let our theoretical assumptions be guided by results from numerical simulations, and show that at any radius where both the radial density profile, rho, and the phase-space-like density profile, rho/sigma^epsilon, are exact power laws, the only allowed density slopes in agreement with the spherical symmetric and isotropic Jeans equation are in the range 1< beta <3, where beta = - dln(rho)/dln(r). We also allow for a radial variation of these power laws, as well as anisotropy, and show how this allows for more shallow central slopes.

Galaxy clusters are not iso-thermal, and the radial temperature dependence will affect the cluster parameters derived through the observation of the Sunyaev-Zeldovich (SZ) effect. We show that the derived peculiar velocity will be systematically shifted by 10-20%. For future all-sky surveys one cannot rely on the observationally expensive X-ray observations to remove this systematic error, but one should instead reach for sufficient angular resolution to perform a deprojection in the SZ spectra. The Compton weighted electron temperature is accurately derived through SZ observations.

Populations of relativistic electrons in clusters of galaxies should imprint a characteristic spectral distortion signature on the Cosmic Microwave Background photons - the relativistic Sunyaev-Zeldovich (SZ) effect. We investigate how sensitive PLANCK and successor experiments will be in detecting or constraining the total optical depth tau of relativistic electrons. We find an expected sensitivity of PLANCK of sigma_tau ~ 2 10^-5 for an individual cluster, which is too insensitive to detect even optimistic scenarios, which predict tau ~ 10^-6 to 10^-7. However, by stacking the PLANCK-signal of all SZ detectable clusters the relativistic SZ effect can be statistically probed to a sensitivity level of sigma_tau ~ 10^-7. PLANCK successor experiments will be able to probe even conservative scenarios of relativistic electron populations in clusters of galaxies.

We consider for the first time the ability of present-day cosmic microwave background (CMB) anisotropies data to determine the primordial helium mass fraction, Y_p. We find that CMB data alone gives the confidence interval 0.160 < Y_p < 0.501 (at 68% c.l.). We analyse the impact on the baryon abundance as measured by CMB and discuss the implications for big bang nucleosynthesis. We identify and discuss correlations between the helium mass fraction and both the redshift of reionization and the spectral index. We forecast the precision of future CMB observations, and find that Planck alone will measure Y_p with error-bars of 5%. We point out that the uncertainty in the determination of the helium fraction will have to be taken into account in order to correctly estimate the baryon density from Planck-quality CMB data.

Galaxies and galaxy clusters are observed to have a rather non-trivial radial behaviour. The observations show that the radial profiles change from one power-law profile near the centre to another power-law profile in the outer region. We present a simple explanation for this complex behaviour by finding the analytical solutions to the governing hydrodynamic equations. We see that the origin of this complexity is the collisional nature of the baryonic plasma, possibly related to a turbulence-enhanced viscosity.

We present a user-friendly tool for the analysis of data from Sunyaev-Zeldovich effect observations. The tool is based on the stochastic method of simulated annealing, and allows the extraction of the central values and error-bars of the 3 SZ parameters, Comptonization parameter, y, peculiar velocity, v_p, and electron temperature, T_e. The f77-code SASZ will allow any number of observing frequencies and spectral band shapes. As an example we consider the SZ parameters for the COMA cluster.

A remarkable result of the Wilkinson Microwave Anisotropy Probe (WMAP) observations is that the universe was significantly reionized at large redshifts. The standard explanation is that massive stars formed early and reionized the universe around redshift z=17. Here we explore an alternative possibility, in which the universe was reionized in two steps. An early boost of reionization is provided by a decaying sterile neutrino, whose decay products, relativistic electrons, result in partial ionization of the smooth gas. We demonstrate that a neutrino with a mass of m_nu ~ 200 MeV and a decay time of t ~ 4 * 10^{15} s can account for the electron scattering optical depth tau=0.16 measured by WMAP without violating existing astrophysical limits on the cosmic microwave and gamma ray backgrounds. Reionization is then completed by subsequent star formation at lower redshifts. This scenario alleviates constraints on structure formation models with reduced small-scale power, such as those with a running or tilted scalar index, or warm dark matter models.

We study the origin of cores and density profiles of gaseous baryonic structures in cosmology. By treating the baryons as a viscous gas, we find that both spheres and disks are possible solutions. We find analytically that the density profiles have inner and outer solutions, which in general are different. For disks we identify a central core, with density profile rho_d = constant, and the outer profile rho_d ~ r^{-3}. For spherical structures we find the profile rho_s ~ r^{-6}. In the presence of a dominating central black hole we find the inner profile rho ~ r^{-3/2}. When the mass is dominated by a dark matter component then the baryonic density profile will depend on the dark matter profile, and we point out how one can use this connection to infer the DM profile directly by observing the baryonic density profile.

We present here results on the composite luminosity functions of galaxies in the clusters of galaxies selected from the SDSS Cut and Enhance cluster catalog (CE; Goto et al. 2001). We construct the composite luminosity function in the five SDSS bands, u,g,r,i and z, using 204 CE clusters ranging from z=0.02 to z=0.25. We use photometric redshifts to construct composite luminosity functions. The robustness of the weighting scheme was tested using Monte Carlo simulation. We find the slope of composite LFs become flatter toward redder color band. Comparing with field LFs of SDSS (Blanton et al. 2001), cluster LFs have brighter characteristic magnitude and flatter slopes in g, r, i and z bands. These results are consistent with the hypothesis that the cluster LF has two distinct underlying populations i.e. the bright end of the LF is dominated by bright early types that follow a gaussian--like luminosity distribution, while the faint--end of the cluster LF is a steep power--law like function dominated by star--forming (bluer) galaxies. We also study the composite luminosity functions for early--type and late--type galaxies using profile fits, a concentration parameter and u-r color to classify galaxy morphology. The strong dependence of LF on galaxy morphology is found.The faint end slope of the LF is always flatter for early--type galaxies than late--type regardless of passband and methodology. These results are consistent with the hypothesis that the cluster region are dominated by bright elliptical galaxies.

We investigate the evolution of the fractions of late type cluster galaxies as a function of redshift, using one of the largest, most uniform cluster samples available. The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created using a single automated cluster finding algorithm on uniform data from a single telescope, with accurate CCD photometry, thus, minimizing selection biases. We use four independent methods to analyze the evolution of the late type galaxy fraction. Specifically, we select late type galaxies based on: restframe g-r color, u-r color, galaxy profile fitting and concentration index. The first criterion corresponds to the one used in the classical Butcher-Oemler analyses. The last three criteria are more sensitive to the morphological type of the galaxies. In all four cases, we find an increase in the fraction of late type galaxies with increasing redshift, significant at the 99.9% level. The results confirm that cluster galaxies do change colors with redshift (the Butcher-Oemler effect) and, in addition, they change their morphology to later-type toward higher redshift -- indicating a morphological equivalent of the Butcher-Oemler effect. We also find a tendency of richer clusters to have lower fractions of late type galaxies. The trend is consistent with a ram pressure stripping model, where richer clusters have more effective ram pressure due to their higher temperature.

We present a catalog of 799 clusters of galaxies in the redshift range z_est = 0.05 - 0.3 selected from ~400 deg^2 of early SDSS commissioning data along the celestial equator. The catalog is based on merging two independent selection methods -- a color-magnitude red-sequence maxBCG technique (B), and a Hybrid Matched-Filter method (H). The BH catalog includes clusters with richness \Lambda >= 40 (Matched-Filter) and N_gal >= 13 (maxBCG), corresponding to typical velocity dispersion of \sigma_v >~ 400 km s^{-1} and mass (within 0.6 h^{-1) Mpc radius) >~ 5*10^{13} h^{-1} M_sun. This threshold is below Abell richness class 0 clusters. The average space density of these clusters is 2*10^{-5} h^3 Mpc^{-3}. All NORAS X-ray clusters and 53 of the 58 Abell clusters in the survey region are detected in the catalog; the 5 additional Abell clusters are detected below the BH catalog cuts. The cluster richness function is determined and found to exhibit a steeply decreasing cluster abundance with increasing richness. We derive observational scaling relations between cluster richness and observed cluster luminosity and cluster velocity dispersion; these scaling relations provide important physical calibrations for the clusters. The catalog can be used for studies of individual clusters, for comparisons with other sources such as X-ray clusters and AGNs, and, with proper correction for the relevant selection functions, also for statistical analyses of clusters.

The Sunyaev-Zel'dovich (SZ) effect of galaxy clusters is characterized by three parameters: Compton parameter, electron temperature and cluster peculiar velocity. In the present study we consider the problem of extracting these parameters using multi-frequency SZ observations only. We show that there exists a parameter degeneracy which can be broken with an appropriate choice of frequencies. As a result we discuss the optimal choice of observing frequencies from a theoretical point of view. Finally, we analyze the systematic errors (of the order micro K) on the SZ measurement introduced by finite bandwidths, and suggest a possible method of reducing these errors.

The mass function of clusters of galaxies is determined from 400 deg^2 of early commissioning imaging data of the Sloan Digital Sky Survey; ~300 clusters in the redshift range z = 0.1 - 0.2 are used. Clusters are selected using two independent selection methods: a Matched Filter and a red-sequence color magnitude technique. The two methods yield consistent results. The cluster mass function is compared with large-scale cosmological simulations. We find a best-fit cluster normalization relation of sigma_8*omega_m^0.6 = 0.33 +- 0.03 (for 0.1 ~< omega_m ~< 0.4), or equivalently sigma_8 = (0.16/omega_m)^0.6. The amplitude of this relation is significantly lower than the previous canonical value, implying that either omega_m is lower than previously expected (omega_m = 0.16 if sigma_8 = 1) or sigma_8 is lower than expected (sigma_8 = 0.7 if omega_m = 0.3). The best-fit mass function parameters are omega_m = 0.19 (+0.08,-0.07) and sigma_8 = 0.9 (+0.3,-0.2). High values of omega_m (>= 0.4) and low sigma_8 (=< 0.6) are excluded at >~ 2 sigma.

We consider the flow-equations for the 3 slow-roll parameters n_S (scalar spectral index), r (tensor to scalar ratio), and dn_S/dlnk (running of the spectral index). We show that the combination of these flow-equations with the observational bounds from cosmic microwave background and large scale structure allows one to put a lower bound on the fourth derivative of the inflationary potential, M_P^4(V''''/V) > -0.02.

We study the effect of imperfect subtraction of the Sunyaev-Zel'dovich effect (SZE) using a robust and non-parametric method to estimate the SZE residual in the Planck channels. We include relativistic corrections to the SZE, and present a simple fitting formula for the SZE temperature dependence for the Planck channels. We show how the relativistic corrections constitute a serious problem for the estimation of the kinematic SZE component from Planck data, since the key channel to estimate the kinematic component of the SZE, at 217 GHz, will be contaminated by a non-negligible thermal SZE component. The imperfect subtraction of the SZE will have an effect on both the Planck cluster catalogue and the recovered CMB map. In the cluster catalogue, the relativistic corrections are not a major worry for the estimation of the total cluster flux of the thermal SZE component, however, they must be included in the SZE simulation when calculating the selection function and completeness level. The power spectrum of the residual at 353 GHz, where the intensity of the thermal SZE is maximum, does not contribute significantly to the power spectrum of the CMB. We calculate the non-Gaussian signal due to the SZE residual in the 353 GHz CMB map using a simple Gaussianity estimator, and this estimator detects a 4.25-sigma non-Gaussian signal at small scales, which could be mistaken for a primordial non-Gaussian signature. The other channels do not show any significant departure from Gaussianity with our estimator.

We explore some of the consequences of Dark Matter-photon interactions on structure formation, focusing on the evolution of cosmological perturbations and performing both an analytical and a numerical study. We compute the cosmic microwave background anisotropies and matter power spectrum in this class of models. We find, as the main result, that when Dark Matter and photons are coupled, Dark Matter perturbations can experience a new damping regime in addition to the usual collisional Silk damping effect. Such Dark Matter particles (having quite large photon interactions) behave like Cold Dark Matter or Warm Dark Matter as far as the cosmic microwave background anisotropies or matter power spectrum are concerned, respectively. These Dark Matter-photon interactions leave specific imprints at sufficiently small scales on both of these two spectra, which may allow to put new constraints on the acceptable photon-Dark Matter interactions. Under the conservative assumption that the abundance of 10^12 M_sol galaxies is correctly given by Cold Dark Matter, and without any knowledge of the abundance of smaller objects, we obtain the limit on the ratio of the Dark Matter-photon cross section to the Dark Matter mass sigma_{gamma-DM} / m_DM < 10^-6 sigma_Thomson / 100 GeV \sim 6 * 10^-33 cm^2 GeV^-1 .

Direct observations provide constraints on the first two derivatives of the inflaton potential in slow roll models. We discuss how present day observations, combined with the flow equations in slow roll parameter space, provide a non-trivial constraint on the third derivative of the inflaton potential. We find a lower bound on the third derivative of the inflaton potential V'''/V > -0.2. We also show that unless the third derivative of the inflaton potential is unreasonably large, then one predicts the tensor to scalar ratio, r, to be bounded from below r > 3 * 10^{-6}.

We discuss a new method of finding the cluster temperatures which is independent of distance and therefore very useful for distant clusters. The hot gas of electrons in clusters of galaxies scatters and distorts the cosmic microwave background radiation in a well determined way. This Sunyaev-Zel'dovich (SZ) effect is a useful tool for extracting information about clusters such as their peculiar radial velocity and optical depth. Here we show how the temperature of the cluster can be inferred from the SZ effect, in principle without use of X-ray data. We use recent millimetre observation of Abell 2163 to determine for the first time a cluster temperature using SZ observations only. The result T_e = 26^+34_-19 keV at 68% confidence level (at 95% c.l. we find T>1.5 keV) is in reasonable agreement with the X-ray results, T_e =12.4^+2.8_-1.9 keV.

We study three-flavor neutrino oscillations in the early universe in the presence of neutrino chemical potentials. We take into account all effects from the background medium, i.e. collisional damping, the refractive effects from charged leptons, and in particular neutrino self-interactions that synchronize the neutrino oscillations. We find that effective flavor equilibrium between all active neutrino species is established well before the big-bang nucleosynthesis (BBN) epoch if the neutrino oscillation parameters are in the range indicated by the atmospheric neutrino data and by the large mixing angle (LMA) MSW solution of the solar neutrino problem. For the other solutions of the solar neutrino problem, partial flavor equilibrium may be achieved if the angle theta_13 is close to the experimental limit tan^2(theta_13)<0.065. In the LMA case, the BBN limit on the nu_e degeneracy parameter, |xi_nu|< 0.07, now applies to all flavors. Therefore, a putative extra cosmic radiation contribution from degenerate neutrinos is limited to such low values that it is neither observable in the large-scale structure of the universe nor in the anisotropies of the cosmic microwave background radiation. Existing limits and possible future measurements, for example in KATRIN, of the absolute neutrino mass scale will provide unambiguous information on the cosmic neutrino mass density, essentially free of the uncertainty of the neutrino chemical. potentials.

Sterile neutrinos may be one of the best Warm Dark Matter candidates we have today. Both lower and upper bounds on the mass of the sterile neutrino come from astronomical observations. We show that the proper inclusion of the neutrino momentum distribution reduces the allowed region to be 2.6 keV < m < 5 keV for the simplest models. A search for a spectral line with E=m/2 is thus more interesting than ever before.

We combine complementary datasets to constrain dark energy. Using standard Big Bang Nucleosynthesis and the observed abundances of primordial nuclides to put constraints on $\Omega_Q$ at temperatures near $T \sim 1MeV$, we find the strong constraint $\Omega_Q({MeV}) < 0.045$ at $2\sigma$ c.l.. Under the assumption of flatness, using results from Cosmic Microwave Background (CMB) anisotropy measurements, high redshift supernovae (SN-Ia) observations and data from local cluster abundances we put a new constraint on the equation of state parameter $w_Q < -0.85$ at 68% c.l..

The presence of dark matter in the solar neighbourhood can be tested in the framework of stellar evolution theory by using the new results of helioseismology. If weakly interacting massive particles accumulate in the center of the Sun, they can provide an additional mechanism for transferring energy from the solar core. The presence of these particles produces a change in the local luminosity of the Sun of the order of 0.1%, which is now within the reach of seismic solar experiments. We find that typical WIMPs with a mass of 60 GeV, annihilation cross-section of 10^{-32} cm^3/sec, and scalar scattering interaction between 10^{-36} cm^2 and 10^{-41} cm^2, the lower end of this range being that suggested by the DAMA direct detection experiment, modify the radial profile of the square of the sound speed in the solar core by \~ 0.1 %. This level of change is strongly constrained by the most recent helioseismological results. Significantly more massive WIMPs have a much smaller effect on the solar core temperature profile. However future helioseismological experiments have the potential exploring most, if not all, of the WIMP parameter space accessible to current and planned direct detection experiments.

Recent estimates of cosmological parameters derived from Cosmic Microwave Background (CMB) anisotropies are based on the assumption that we know the precise amount of energy density in relativistic particles in the universe, $\omega_{rel}$, at all times. There are, however, many possible mechanisms that can undermine this assumption. In this paper we investigate the effect that removing this assumption has on the determination of the various cosmological parameters. We obtain fairly general bounds on the redshift of equality, $z_{eq}= \omega_{m}/\omega_{rel}=3100_{-400}^{+600}$. We show that $\omega_{rel}$ is nearly degenerate with the amount of energy in matter, $\omega_m$, and that its inclusion in CMB parameter estimation also affects the present constraints on other parameters such as the curvature or the scalar spectral index of primordial fluctuations. This degeneracy has the effect of limiting the precision of parameter estimation from the MAP satellite, but it can be broken by measurements on smaller scales such as those provided by the Planck satellite mission.

We constrain the spectrum of primordial curvature perturbations by using recent Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) data. Specifically, we consider CMB data from the COBE, Boomerang and Maxima experiments, the real space galaxy power spectrum from the IRAS PSCz survey, and the linear matter power spectrum inferred from Ly-alpha forest spectra. We study the case of single field slow roll inflationary models, and we extract bounds on the scalar spectral index, n, the tensor to scalar ratio, r, and the running of the scalar spectral index, dlnk, for various combinations of the observational data. We find that CMB data, when combined with data from Lyman-alpha forest, place strong constraints on the inflationary parameters. Specifically, we obtain n \approx 0.9, r < 0.3 and dlnk \approx 0, indicating that single field hybrid models are ruled out.

This paper has been withdrawn by the authors. This paper started with a wrong term describing elastic scattering. The correct form of the kinetic equations for the density matrix are derived e.g. in A. Dolgov, Yad. Fiz. 33 (1981) 1309; K. Enqvist, K. Kainulainen and J. Maalampi, Nucl. Phys. B 349 (1991) 754; G. Sigl and G. Raffelt, Nucl. Phys. B 406 (1993) 423. When using the correct form of the kinetic equations, the result is in good agreement with old results, A. Dolgov and S. Hansen, hep-ph/0009083.

We perform a likelihood analysis of the recent results on the anisotropy of Cosmic Microwave Background Radiation from the BOOMERanG and DASI experiments to show that they single out an effective number of neutrinos in good agreement with standard Big Bang Nucleosynthesis. We also consider degenerate Big Bang Nucleosynthesis to provide new bounds on effective relativistic degrees of freedom $N_\nu$ and, in particular, on neutrino chemical potential $\xi_\alpha$. When including Supernova Ia data we find, at $2\sigma$, $N_\nu \leq 7$ and $-0.01 \leq \xi_e \leq 0.22$, $|\xi_{\mu,\tau}|\leq 2.6$.

In the past years 'quintessence' models have been considered which can produce the accelerated expansion in the universe suggested by recent astronomical observations. One of the key differences between quintessence and a cosmological constant is that the energy density in quintessence, $\Omega_\phi$, could be a significant fraction of the overall energy even in the early universe, while the cosmological constant will be dynamically relevant only at late times. We use standard Big Bang Nucleosynthesis and the observed abundances of primordial nuclides to put constraints on $\Omega_\phi$ at temperatures near $T \sim 1MeV$. We point out that current experimental data does not support the presence of such a field, providing the strong constraint $\Omega_\phi(MeV) < 0.045$ at $2\sigma$ C.L. and strengthening previous results. We also consider the effect a scaling field has on CMB anisotropies using the recent data from Boomerang and DASI, providing the CMB constraint $\Omega_\phi \le 0.39$ at $2\sigma$ during the radiation dominated epoch.

We show that massive sterile neutrinos mixed with the ordinary ones may be produced in the early universe in the right amount to be natural warm dark matter particles. Their mass should be below 40 keV and the corresponding mixing angles sin^2 2\theta > 10^{-11} for mixing with \nu_\mu or \nu_\tau, while mixing with \nu_e is slightly stronger bounded with mass less than 30 keV.

The spectral distortion of the cosmic background radiation produced by the inverse Compton scattering on hot electrons in clusters of galaxies (thermal Sunyaev-Zeldovich effect) is calculated for arbitrary optical depth and electron temperature. The distortion is found by a numerical solution of the exact Boltzmann equation for the photon distribution function. In the limit of small optical depth and low electron temperature our results confirm the previous analyses. In the opposite limits, our method is the only one that permits to make accurate calculations.

We constrain the primordial fluctuation spectrum P(k) by using the new data on the Cosmic Microwave Background (CMB) from the Boomerang and MAXIMA experiments. Our study is based on slow-roll inflationary models, and we consider the possibility of a running spectral index. Specifically, we expand the power spectrum P(k) to second order in ln(k), thus allowing the power spectrum to ``bend'' in k-space. We show that allowing the power spectrum to bend erases the ability of the present data to measure the tensor to scalar perturbation ratio. Moreover, if the primordial baryon density \Omega_b h^2 is as low as found from Big Bang nucleosynthesis (BBN), the data favor a negative bending of the power spectrum, corresponding to a bump-like feature in the power spectrum around a scale of k=0.004 Mpc^-1.

We review current astrophysical bounds on MeV sterile neutrinos, and then we discuss why a sterile keV neutrino is a natural warm dark matter candidate.

Cosmological and astrophysical effects of heavy (10 - 200 MeV) sterile Dirac neutrinos, mixed with the active ones, are considered. The bounds on mass and mixing angle from both supernovae and big-bang nucleosynthesis are presented.

The low second acoustic peak in the recent Boomerang data may indicate a cosmological baryon density which is larger than allowed by standard big bang nucleosynthesis. We show that the decay of the tau-neutrino: nu_tau -> nu_e + phi, where nu_e is the electron neutrino and phi is a scalar, essentially can assure agreement between BBN calculations and light element observations for a large baryon density.

Constraints on the lifetime of the heavy sterile neutrino, that was proposed as a possible interpretation of the KARMEN anomaly, are derived from primordial nucleosynthesis and SN 1987A. Together with the recent experimental bounds on the nu_s lifetime, SN 1987A completely excludes this interpretation. Nucleosynthesis arguments permit a narrow window for the lifetime in the interval 0.1-0.2 sec. If nu_s possesses an anomalous interaction with nucleons, the SN bounds may not apply, while the nucleosynthesis ones would remain valid.

Distortion of the equilibrium spectra of cosmic neutrinos due to interaction with hotter electrons and positrons in the primeval cosmic plasma is considered. The set of integro-differential kinetic equations for neutrinos is accurately numerically solved. The relative corrections to neutrino energy densities are approximately 0.9% for $\nu_e$ and 0.4% for $\nu_\mu$ and $\nu_\tau$. This effect results in $1.4 \cdot 10^{-4}$ increase in the primordial $^4 He$ abundance.