The merging rate of cosmic structures is computed, relying on the Ansatz that they can be predicted in the initial linear density field from the coalescence of critical points with increasing smoothing scale, used here as a proxy for cosmic time. Beyond the mergers of peaks with saddle points (a proxy for halo mergers), we consider the coalescence and nucleation of all sets of critical points, including wall-saddle to filament-saddle and wall-saddle to minima (a proxy for filament and void mergers respectively), as they impact the geometry of galactic infall, and in particular filament disconnection. Analytical predictions of the one-point statistics are validated against multiscale measurements in 2D and 3D realisations of Gaussian random fields (the corresponding code being available upon request) and compared qualitatively to cosmological $N$-body simulations at early times ($z\geq 10$) and large scales ($\geq 5\, \mathrm{Mpc}/h$). The rate of filament coalescence is compared to the merger rate of haloes and the two-point clustering of these events is computed, along with their cross-correlations with critical points. These correlations are qualitatively consistent with the preservation of the connectivity of dark matter haloes, and the impact of the large scale structures on assembly bias. The destruction rate of haloes and voids as a function of mass and redshift is quantified down to $z=0$ for a $\Lambda$CDM cosmology. The one-point statistics in higher dimensions are also presented, together with consistency relations between critical point and critical event counts.

We introduce two new ways of obtaining the strength of plane-of-sky (POS) magnetic field by simultaneous use of spectroscopic Doppler-shifted lines and the information on magnetic field direction. The latter can be obtained either through polarization measurements or using the velocity gradient technique. We show the advantages that our techniques have compared to the traditional Davis-Chandrasekhar-Fermi (DCF) technique of estimating magnetic field strength from observations. The first technique that we describe in detail employs structure functions of velocity centroids and structure functions of Stokes parameters. We provide analytical expressions for obtaining magnetic field strength from observational data. We successfully test our results using synthetic observations obtained with results of MHD turbulence simulations. We measure velocity and magnetic field fluctuations at small scales using two, three and four point structure functions and compare the performance of these tools. We show that, unlike the DCF, our technique is capable of providing the detailed distribution of POS magnetic field and it can measure magnetic field strength in the presence of both velocity and magnetic field distortions arising from external shear and self-gravity. The second technique applies the velocity gradient technique to velocity channel maps in order to obtain the Alfven Mach number and uses the amplitudes of the gradients to obtain the sonic Mach number. The ratio of these two Mach numbers provides the intensity of magnetic field in the region contributing to the emission in the channel map. We test the technique and discuss obtaining the 3D distribution of POS galactic Magnetic field with it. We discuss the application of the second technique to synchrotron data.

We investigate the impact of the number of filaments connected to the nodes of the cosmic web on the physical properties of their galaxies using the Sloan Digital Sky Survey. We compare these measurements to the cosmological hydrodynamical simulations Horizon-(no)AGN and Simba. We find that more massive galaxies are more connected, in qualitative agreement with theoretical predictions and measurements in dark matter only simulation. The star formation activity and morphology of observed galaxies both display some dependence on the connectivity of the cosmic web at fixed stellar mass: less star forming and less rotation supported galaxies also tend to have higher connectivity. These results qualitatively hold both for observed and virtual galaxies, and can be understood given that the cosmic web is the main source of fuel for galaxy growth. The simulations show the same trends at fixed halo mass, suggesting that the geometry of filamentary infall impacts galaxy properties beyond the depth of the local potential well. Based on simulations, it is also found that AGN feedback is key in reversing the relationship between stellar mass and connectivity at fixed halo mass. Technically, connectivity is a practical observational proxy for past and present accretion (minor mergers or diffuse infall).

In modern cosmology, many efforts have been put to detect primordial B-mode of cosmic microwave background (CMB) polarization from the gravitational waves generated during inflation. Considering the foreground dust contamination of microwave polarization maps, it is essential to obtain a precise prediction for polarization in dust emission. In this work, we show a new method to produce synthetic maps of dust polarization in magnetized turbulent ISM from more abundant high-resolution HI data. By using Velocity Channel Gradient (VChG) technique, we are able to predict both direction and degree of dust polarization by investigating spectroscopic HI information in position-position-velocity (PPV) space. We applied our approach to The Galactic Arecibo L-band feed Array HI (GALFA-HI) data, and find a good correspondence between synthesized maps and PLANCK's polarization measurements at 353 GHz.

A recent publication by Clark et.al (2019, CX19) uses both GALFA-HI observational data and numerical simulations to address the nature of intensity fluctuations in Position-Position-Velocity (PPV) space. The study questions the validity and applicability of the statistical theory of PPV space fluctuations formulated in Lazarian & Pogosyan (2000, LP00) to HI gas and concludes that $"${\it a significant reassessment of many observational and theoretical studies of turbulence in HI}$"$. This implies that dozens of papers that used LP00 theory to explore interstellar turbulence as well as the ongoing research based on LP00 theory are in error. This situation motivates the urgency of our public response. In our Comment we explain why we believe the criticism in CX19 is based on the incorrect understanding of the LP00 theory. In particular, we illustrate that the correlation between PPV slices and dust emissions in CX19 does not properly reveal the relative importance of velocity and density fluctuations in velocity channel maps. While CX19 provides an explanation of the change of the spectral index with respect to the thickness of PPV slice based on the two-phase nature of H1 gas, we failed to see any observational support for this idea. On the contrary, we show that the observations both in two-phase HI and one phase CO show similar results. Moreover, the observed change is in good agreement with LP00 predictions and spectral indexes of velocity and density spectra that are obtained following LP00 procedures are in good agreement with the numerically confirmed expectations of compressible MHD turbulence theory. In short, we could not find any justification of the criticism of LP00 theory that is provided in CX19. On the contrary, our analysis testifies that both available observational and numerical data agree well with the predictions of LP00 theory.

Cosmic connectivity and multiplicity, i.e. the number of filaments globally or locally connected to a given cluster is a natural probe of the growth of structure and in particular of the nature of dark energy. It is also a critical ingredient driving the assembly history of galaxies as it controls mass and angular momentum accretion. The connectivity of the cosmic web is investigated here via the persistent skeleton. This tool identifies topologically the set of ridges of the cosmic landscape which allows us to investigate in details how the nodes of the cosmic web are connected together. When applied to Gaussian random fields, it is found that on average the nodes are connected to exactly $\kappa=4$ neighbours in two dimensions and ~6.1 in three dimensions. Investigating spatial dimensions up to d=11, typical departures from a cubic lattice $\kappa=2d$ are shown to scale like the power 7/4 of the dimension. These numbers strongly depend on the height of the peaks: the higher the peak the larger the connectivity. Locally, the number of filaments sticking out from a node is not equal to its number of connections. At small distance, the peak has an ellipsoidal shape so that only two filaments emerge. Further away, filaments are subdivided at bifurcation points. At the resolution of the maps, the local number of filaments is 3 on average in 2D and 4 in 3D, depending on the height of the peak. Predictions from first principles based on peak theory are shown to reproduce well the connectivity and multiplicity of Gaussian random fields and cosmological simulations. As an illustration, this connectivity is quantified in galaxy lensing convergence maps and large dark haloes catalogues. As a function of redshift and scale the mean connectivity decreases in a cosmology dependent way. As a function of halo mass it scales like 10/3 times the log of the mass.

In this work we effectively calculated the two-dimensional Minkowski functionals for cosmic microwave background (CMB) temperature maps generated by single field models of inflation with a standard kinetic term. We started with calculation of the bispectrum of initial perturbations and then calculated the two-dimensional configuration space cubic moments for temperature fluctuations. These cubic moments give rise to first order non-Gaussian correction terms to the Minkowski functionals. Thus, we developed a robust mechanism to predict the amount of non-Gaussianity generated by inflation in the CMB temperature maps using Minkowski functionals.

Recent Planck observations have revealed some of the important statistical properties of synchrotron and dust polarisation, namely, the $B$ to $E$ mode power and temperature-$E$ (TE) mode cross-correlation. In this paper, we extend our analysis in Kandel et al. (2017) that studied $B$ to $E$ mode power ratio for polarised dust emission to include TE cross-correlation and develop an analogous formalism for synchrotron signal, all using a realistic model of magnetohydrodynamical (MHD) turbulence. Our results suggest that the Planck results for both synchrotron and dust polarisation can be understood if the turbulence in the Galaxy is sufficiently sub-Alfv\'enic. We also show how $B$ to $E$ ratio as well as the TE cross-correlation can be used to study media magnetisation, compressibility, and level of density-magnetic field correlation.

The role of the cosmic web in shaping galaxy properties is investigated in the GAMA spectroscopic survey in the redshift range $0.03 \leq z \leq 0.25$. The stellar mass, $u - r$ dust corrected colour and specific star formation rate (sSFR) of galaxies are analysed as a function of their distances to the 3D cosmic web features, such as nodes, filaments and walls, as reconstructed by DisPerSE. Significant mass and type/colour gradients are found for the whole population, with more massive and/or passive galaxies being located closer to the filament and wall than their less massive and/or star-forming counterparts. Mass segregation persists among the star-forming population alone. The red fraction of galaxies increases when closing in on nodes, and on filaments regardless of the distance to nodes. Similarly, the star-forming population reddens (or lowers its sSFR) at fixed mass when closing in on filament, implying that some quenching takes place. Comparable trends are also found in the state-of-the-art hydrodynamical simulation Horizon-AGN. These results suggest that on top of stellar mass and large-scale density, the traceless component of the tides from the anisotropic large-scale environment also shapes galactic properties. An extension of excursion theory accounting for filamentary tides provides a qualitative explanation in terms of anisotropic assembly bias: at a given mass, the accretion rate varies with the orientation and distance to filaments. It also explains the absence of type/colour gradients in the data on smaller, non-linear scales.

Recent Planck observations of dust polarization in the Galaxy have revealed that the power in $E$ mode is twice that in $B$ mode. Caldwell et. al. have formulated a theoretical model in the context of magnetohydrodynamic (MHD) turbulence and found it problematic to account for this result. In particular, they concluded that there is a very narrow range of theoretical parameters that could account for the observation. This poses a problem of whether the accepted description of MHD turbulence can apply to the interstellar medium. We revisit the problem and demonstrate that MHD turbulence corresponding to the high galactic latitudes range of Alfv\'en Mach numbers, i.e. $M_A\lesssim 0.5$, can successfully explain the available results for the $E$ to $B$ mode ratio.

We study the effect of dust absorption on the recovery velocity and density spectra as well as on the anisotropies of magnetohydrodynamic turbulence using the Velocity Channel Analysis (VCA), Velocity Coordinate Spectrum (VCS) and Velocity Centroids. The dust limits volume up to an optical depth of unity. We show that in the case of the emissivity proportional to the density of emitters, the effects of random density get suppressed for strong dust absorption intensity variations arise from the velocity fluctuations only. However, for the emissivity proportional to squared density, both density and velocity fluctuations affect the observed intensities. We predict a new asymptotic regime for the spectrum of fluctuations for large scales exceeding the physical depths to unit optical depth. The spectrum gets shallower by unity in this regime. In addition, the dust absorption removes the degeneracy resulted in the universal $K^{-3}$ spectrum of intensity fluctuations of self-absorbing medium reported by Lazarian \& Pogosyan. We show that the predicted result is consistent with the available HII region emission data. We find that for sub-Alfv\'enic and trans-Alfv\'enic turbulence one can get the information about both the magnetic field direction and the fundamental Alfv\'en, fast and slow modes that constitute MHD turbulence.

Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is established using a parametrisation-independent extraction of the bias function. This average bias model is then combined with the dark matter PDF, neglecting any scatter around it: it nevertheless yields an excellent model for densities-in-cells statistics of mass tracers that is parametrised in terms of the underlying dark matter variance and three bias parameters. The procedure is validated on measurements of both the one and two point statistics of subhalo densities in the state-of-the-art Horizon Run 4 simulation showing excellent agreement for measured dark matter variance and bias parameters. Finally, it is demonstrated that this formalism allows for a joint estimation of the nonlinear dark matter variance and the bias parameters using solely the statistics of subhaloes. Having verified that galaxy counts in hydrodynamical simulations sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work provides important steps towards making theoretical predictions for density-in-cells statistics applicable to upcoming galaxy surveys like Euclid or WFIRST.

Simple parameter-free analytic bias functions for the two-point correlation of densities in spheres at large separation are presented. These bias functions generalize the so-called Kaiser bias to the mildly non-linear regime for arbitrary density contrasts. The derivation is carried out in the context of large deviation statistics while relying on the spherical collapse model. A logarithmic transformation provides a saddle approximation which is valid for the whole range of densities and shown to be accurate against the 30 Gpc cube state-of-the-art Horizon Run 4 simulation. Special configurations of two concentric spheres that allow to identify peaks are employed to obtain the conditional bias and a proxy to BBKS extrema correlation functions. These analytic bias functions should be used jointly with extended perturbation theory to predict two-point clustering statistics as they capture the non-linear regime of structure formation at the percent level down to scales of about 10 Mpc/h at redshift 0. Conversely, the joint statistics also provide us with optimal dark matter two-point correlation estimates which can be applied either universally to all spheres or to a restricted set of biased (over- or underdense) pairs. Based on a simple fiducial survey, this estimator is shown to perform five times better than usual two-point function estimators. Extracting more information from correlations of different types of objects should prove essential in the context of upcoming surveys like Euclid, DESI, PFS or LSST.

We study possibility of obtaining velocity spectra by studying turbulence in an optically thick medium using velocity centroids (VCs).We find that the regime of universal, i.e. independent of underlying turbulence statistics, fluctuations discovered originally within the velocity channel analysis (VCA) carries over to the statistics of VCs. In other words, for large absorptions the VC lose their ability to reflect the spectra of turbulence. Combining our present study with the earlier studies of centroids in Esquivel & Lazarian, we conclude that centroids are applicable for studies subsonic/transsonic turbulence for the range of scales that is limited by the absorption effects. We also consider VCs based on absorption lines and define the range of their applicability. We address the problem of analytical description of spectra and anisotropies of fluctuations that are available through studies using VC. We obtain spectra and anisotropy of VC fluctuations arising from Alfv\'en, slow and fast modes that constitute the compressible MHD cascade to address the issue of anisotropy of VC statistics, and show how the VC anisotropy can be used to find the media magnetization as well as to identify and separate contributions from Alfv\'en, slow and fast modes. Our study demonstrates that VCs are complementary to the tools provided by the VCA. In order to study turbulent volume for which the resolution of single dish telescopes is insufficient, we demonstrate how the studies of anisotropy can be performed using interferometers.

We extend the velocity channel analysis (VCA), introduced by Lazarian & Pogosyan, of the intensity fluctuations in the velocity slices of position-position-velocity (PPV) spectroscopic data from Doppler broadened lines to study statistical anisotropy of the underlying velocity and density that arises in a turbulent medium from the presence of magnetic field. In particular, we study analytically how the anisotropy of the intensity correlation in the channel maps changes with the thickness of velocity channels. In agreement with the earlier VCA studies we find that the anisotropy in the thick channels reflects the anisotropy of the density field, while the relative contribution of density and velocity fluctuations to the thin velocity channels depends on the density spectral slope. We show that the anisotropies arising from Alfven, slow and fast magnetohydrodynamical modes are different, in particular, the anisotropy in PPV created by fast modes is opposite to that created by Alfven and slow modes, and this can be used to separate their contributions. We successfully compare our results with the recent numerical study of the PPV anisotropies measured with synthetic observations. We also extend our study to the medium with self-absorption as well as to the case of absorption lines. In addition, we demonstrate how the studies of anisotropy can be performed using interferometers.

In the conference presentation we have reviewed the theory of non-Gaussian geometrical measures for the 3D Cosmic Web of the matter distribution in the Universe and 2D sky data, such as Cosmic Microwave Background (CMB) maps that was developed in a series of our papers. The theory leverages symmetry of isotropic statistics such as Minkowski functionals and extrema counts to develop post- Gaussian expansion of the statistics in orthogonal polynomials of invariant descriptors of the field, its first and second derivatives. The application of the approach to 2D fields defined on a spherical sky was suggested, but never rigorously developed. In this paper we present such development treating effects of the curvature and finiteness of the spherical space $S_2$ exactly, without relying on the flat-sky approximation. We present Minkowski functionals, including Euler characteristic and extrema counts to the first non-Gaussian correction, suitable for weakly non-Gaussian fields on a sphere, of which CMB is the prime example.

Full-sky CMB maps from the 2015 Planck release allow us to detect departures from global isotropy on the largest scales. We present the first searches using CMB polarization for correlations induced by a non-trivial topology with a fundamental domain intersecting, or nearly intersecting, the last scattering surface (at comoving distance $\chi_{rec}$). We specialize to flat spaces with toroidal and slab topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology at a scale below the diameter of the last scattering surface. The limits on the radius $R_i$ of the largest sphere inscribed in the topological domain (at log-likelihood-ratio $\Delta\ln{L}>-5$ relative to a simply-connected flat Planck best-fit model) are $R_i>0.97\chi_{rec}$ for the cubic torus and $R_i>0.56\chi_{rec}$ for the slab. The limit for the cubic torus from the matched-circles search is numerically equivalent, $R_i>0.97\chi_{rec}$ (99% CL) from polarisation data alone. We also perform a Bayesian search for a Bianchi VII$_h$ geometry. In the non-physical setting where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component. However, the cosmological parameters generating this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to Planck data requires an amplitude of $-0.1\pm0.04$ compared to +1 if the model were to be correct. In the physical setting where the Bianchi parameters are fit simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models to $(\omega/H)_0<7.6\times10^{-10}$ (95% CL). [Abridged]

We consider turbulent synchrotron emitting media that also exhibits Faraday rotation and provide a statistical description of synchrotron polarization fluctuations. In particular, we consider these fluctuations as a function of the spatial separation of the direction of measurements and as a function of wavelength for the same line-of-sight. On the basis of our general analytical approach, we introduce several measures that can be used to obtain the spectral slopes and correlation scales of both the underlying magnetic turbulence responsible for emission and the spectrum of the Faraday rotation fluctuations. We show the synergetic nature of these measures and discuss how the study can be performed using sparsely sampled interferometric data. We also discuss how additional characteristics of turbulence can be obtained, including the turbulence anisotropy, the three dimensional direction of the mean magnetic field. We consider both cases when the synchrotron emission and Faraday rotation regions coincide and when they are spatially separated. Appealing to our earlier study in Lazarian and Pogosyan (2012) we explain that our new results are applicable to a wide range of spectral indexes of relativistic electrons responsible for synchrotron emission. We expect wide application of our techniques both with existing synchrotron data sets as well as with big forthcoming data sets from LOFAR and SKA.

In this work we investigate the Principal Component Analysis (PCA) sensitivity to the velocity power spectrum in high opacity regimes of the interstellar medium (ISM). For our analysis we use synthetic Position-Position-Velocity (PPV) cubes of fractional Brownian motion (fBm) and magnetohydrodynamics (MHD) simulations, post processed to include radiative transfer effects from CO. We find that PCA analysis is very different from the tools based on the traditional power spectrum of PPV data cubes. Our major finding is that PCA is also sensitive to the phase information of PPV cubes and this allows PCA to detect the changes of the underlying velocity and density spectra at high opacities, where the spectral analysis of the maps provides the universal -3 spectrum in accordance with the predictions of Lazarian \& Pogosyan (2004) theory. This makes PCA potentially a valuable tool for studies of turbulence at high opacities provided that the proper gauging of the PCA index is made. The later, however, we found to be not easy, as the PCA results change in an irregular way for data with high sonic Mach numbers. This is in contrast to synthetic Brownian noise data used for velocity and density fields that show monotonic PCA behavior. We attribute this difference to the PCA's sensitivity to Fourier phase information.

We have calculated the third order moments of scalar curvature perturbations in configuration space for different inflationary models. We developed a robust numerical technique to compute the bispectrum for different models that have some features in the inflationary potential. From the bispectrum we evaluated moments analytically in the slow roll regime while we devised a numerical mechanism to calculated these moments for non slow roll single field inflationary models with standard kinetic term that are minimally coupled to gravity. With help of these third order moments one can directly predict many non-Gaussian and geometrical measures of CBM distributions in the configuration space. Thus, we have devised a framework to calculate different third order moments and geometrical measures, e.g. Minkowski functionals or skeleton statistic, generated by different single field models of inflation.

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

The geometry of the cosmic web drives in part the spin acquisition of galaxies. This can be explained in a Lagrangian framework, by identifying the specific long-wavelength correlations within the primordial Gaussian random field which are relevant to spin acquisition. Tidal Torque Theory is revisited in the context of such anisotropic environments, biased by the presence of a filament within a wall. The point process of filament-type saddles represents it most efficiently. The constrained misalignment between the tidal and the inertia tensors in the vicinity of filament-type saddles simply explains the distribution of spin directions. This misalignment implies in particular an azimuthal orientation for the spins of more massive galaxies and a spin alignment with the filament for less massive galaxies. This prediction is found to be in qualitative agreement with measurements in Gaussian random fields and N-body simulations. It relates the transition mass to the geometry of the saddle, and accordingly predicts its measured scaling with the mass of non-linearity. Implications for galaxy formation and weak lensing are briefly discussed, as is the dual theory of spin alignments in walls.

The intrinsic alignment of galaxy shapes (by means of their angular momentum) and their cross-correlation with the surrounding dark matter tidal field are investigated using the 160 000, z=1.2 synthetic galaxies extracted from the high-resolution cosmological hydrodynamical simulation Horizon-AGN. One- and two-point statistics of the spin of the stellar component are measured as a function of mass and colour. For the low-mass galaxies, this spin is locally aligned with the tidal field `filamentary' direction while, for the high-mass galaxies, it is perpendicular to both filaments and walls. The bluest galaxies of our synthetic catalog are more strongly correlated with the surrounding tidal field than the reddest galaxies, and this correlation extends up to 10 Mpc/h comoving distance. We also report a correlation of the projected ellipticities of blue, intermediate mass galaxies on a similar scale at a level of 10^(-4) which could be a concern for cosmic shear measurements. We do not report any measurable intrinsic alignments of the reddest galaxies of our sample. This work is a first step toward the use of very realistic catalog of synthetic galaxies to evaluate the contamination of weak lensing measurement by the intrinsic galactic alignments.

The kinematic analysis of dark matter and hydrodynamical simulations suggests that the vorticity in large-scale structure is mostly confined to, and predominantly aligned with their filaments, with an excess of probability of 20 per cent to have the angle between vorticity and filaments direction lower than 60 degrees relative to random orientations. The cross sections of these filaments are typically partitioned into four quadrants with opposite vorticity sign, arising from multiple flows, originating from neighbouring walls. The spins of halos embedded within these filaments are consistently aligned with this vorticity for any halo mass, with a stronger alignment for the most massive structures up to an excess of probability of 165 per cent. On large scales, adiabatic/cooling hydrodynamical simulations display the same vorticity in the gas as in the dark matter. The global geometry of the flow within the cosmic web is therefore qualitatively consistent with a spin acquisition for smaller halos induced by this large-scale coherence, as argued in Codis et al. (2012). In effect, secondary anisotropic infall (originating from the vortex-rich filament within which these lower-mass halos form) dominates the angular momentum budget of these halos. The transition mass from alignment to orthogonality is related to the size of a given multi-flow region with a given polarity. This transition may be reconciled with the standard tidal torque theory if the latter is augmented so as to account for the larger scale anisotropic environment of walls and filaments.

Filaments of the cosmic web drive spin acquisition of disc galaxies. The point process of filament-type saddle represent best this environment and can be used to revisit the Tidal Torque Theory in the context of an anisotropic peak (saddle) background split. The constrained misalignment between the tidal tensor and the Hessian of the density field generated in the vicinity of filament saddle points simply explains the corresponding transverse and longitudinal point-reflection symmetric geometry of spin distribution. It predicts in particular an azimuthal orientation of the spins of more massive galaxies and spin alignment with the filament for less massive galaxies. Its scale dependence also allows us to relate the transition mass corresponding to the alignment of dark matter halos spin relative to the direction of their neighboring filament to this geometry, and to predict accordingly it s scaling with the mass of non linearity, as was measured in simulations.

A large-scale hydrodynamical cosmological simulation, Horizon-AGN, is used to investigate the alignment between the spin of galaxies and the cosmic filaments above redshift 1.2. The analysis of more than 150 000 galaxies per time step in the redshift range 1.2<z<1.8 with morphological diversity shows that the spin of low-mass blue galaxies is preferentially aligned with their neighbouring filaments, while high-mass red galaxies tend to have a perpendicular spin. The reorientation of the spin of massive galaxies is provided by galaxy mergers, which are significant in their mass build-up. We find that the stellar mass transition from alignment to misalignment happens around 3.10^10 M_sun. Galaxies form in the vorticity-rich neighbourhood of filaments, and migrate towards the nodes of the cosmic web as they convert their orbital angular momentum into spin. The signature of this process can be traced to the properties of galaxies, as measured relative to the cosmic web. We argue that a strong source of feedback such as active galactic nuclei is mandatory to quench in situ star formation in massive galaxies and promote various morphologies. It allows mergers to play their key role by reducing post-merger gas inflows and, therefore, keeping spins misaligned with cosmic filaments.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

In the context of upcoming large-scale structure surveys such as Euclid, it is of prime importance to quantify the effect of peculiar velocities on geometric probes. Hence the formalism to compute in redshift space the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmic fields is developed. Leveraging the partial isotropy of the target statistics, the Gram-Charlier expansion of the joint probability distribution of the field and its derivatives is reformulated in terms of the corresponding anisotropic variables. In particular, the cosmic non-linear evolution of the Minkowski functionals, together with the statistics of extrema are investigated in turn for 3D catalogues and 2D slabs. The amplitude of the non-Gaussian redshift distortion correction is estimated for these geometric probes. In 3D, gravitational perturbation theory is implemented in redshift space to predict the cosmic evolution of all relevant Gram-Charlier coefficients. Applications to the estimation of the cosmic parameters sigma(z) and beta=f/b1 from upcoming surveys is discussed. Such statistics are of interest for anisotropic fields beyond cosmology.

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

Planck CMB temperature maps allow detection of large-scale departures from homogeneity and isotropy. We search for topology with a fundamental domain nearly intersecting the last scattering surface (comoving distance $\chi_r$). For most topologies studied the likelihood maximized over orientation shows some preference for multi-connected models just larger than $\chi_r$. This effect is also present in simulated realizations of isotropic maps and we interpret it as the alignment of mild anisotropic correlations with chance features in a single realization; such a feature can also exist, in milder form, when the likelihood is marginalized over orientations. Thus marginalized, the limits on the radius $R_i$ of the largest sphere inscribed in a topological domain (at log-likelihood-ratio -5) are: in a flat Universe, $R_i>0.9\chi_r$ for the cubic torus (cf. $R_i>0.9\chi_r$ at 99% CL for a matched-circles search); $R_i>0.7\chi_r$ for the chimney; $R_i>0.5\chi_r$ for the slab; in a positively curved Universe, $R_i>1.0\chi_r$ for the dodecahedron; $R_i>1.0\chi_r$ for the truncated cube; $R_i>0.9\chi_r$ for the octahedron. Similar limits apply to alternate topologies. We perform a Bayesian search for an anisotropic Bianchi VII$_h$ geometry. In a non-physical setting where the Bianchi parameters are decoupled from cosmology, Planck data favour a Bianchi component with a Bayes factor of at least 1.5 units of log-evidence: a Bianchi pattern is efficient at accounting for some large-scale anomalies in Planck data. However, the cosmological parameters are in strong disagreement with those found from CMB anisotropy data alone. In the physically motivated setting where the Bianchi parameters are fitted simultaneously with standard cosmological parameters, we find no evidence for a Bianchi VII$_h$ cosmology and constrain the vorticity of such models: $(\omega/H)_0<8\times10^{-10}$ (95% CL). [Abridged]

We study Minkowski Functionals as probes of primordial non-Gaussianity in the Cosmic Microwave Background, specifically for the estimate of the primordial `local' bi-spectrum parameter f_NL, with instrumental parameters which should be appropriate for the Planck experiment. We use a maximum likelihood approach, which we couple with various filtering methods and test thoroughly for convergence. We included the effect of inhomogeneous noise as well as astrophysical biases induced by point sources and by the contamination from the Galaxy. We find that, when Wiener filtered maps are used (rather than simply smoothed with Gaussian), the expected error on the measurement of f_NL should be as small as \Delta f_NL \simeq 10 when combining the 3 channels at 100, 143 and 217 GHz in the Planck extended mission setup. This result is fairly insensitive to the non homogeneous nature of the noise, at least for realistic hit-maps expected from Planck. We then estimate the bias induced on the measurement of f_NL by point sources in those 3 channels. With the appropriate masking of the bright sources, this bias can be reduced to a negligible level in the 100 and 143 GHz channels. It remains significant in the 217 GHz channel, but can be corrected for. The galactic foreground biases are quite important and present a complex dependence on sky coverage: making them negligible will depend strongly on the quality of the component separation methods.

We investigate the alignment of the spin of dark matter halos relative (i) to the surrounding large-scale filamentary structure, and (ii) to the tidal tensor eigenvectors using the Horizon 4pi dark matter simulation which resolves over 43 million dark matter halos at redshift zero. We detect a clear mass transition: the spin of dark matter halos above a critical mass tends to be perpendicular to the closest filament, and aligned with the intermediate axis of the tidal tensor, whereas the spin of low-mass halos is more likely to be aligned with the closest filament. Furthermore, this critical mass of 5 10^12 is redshift-dependent and scales as (1+z)^-2.5. We propose an interpretation of this signal in terms of large-scale cosmic flows. In this picture, most low-mass halos are formed through the winding of flows embedded in misaligned walls; hence they acquire a spin parallel to the axis of the resulting filaments forming at the intersection of these walls. On the other hand, more massive halos are typically the products of later mergers along such filaments, and thus they acquire a spin perpendicular to this direction when their orbital angular momentum is converted into spin. We show that this scenario is consistent with both the measured excess probabilities of alignment w.r.t. the eigen-directions of the tidal tensor, and halo merger histories. On a more qualitative level, it also seems compatible with 3D visualization of the structure of the cosmic web as traced by "smoothed" dark matter simulations or gas tracer particles. Finally, it provides extra support to the disc forming paradigm presented by Pichon et al (2011) as it extends it by characterizing the geometry of secondary infall at high redshift.

We use cosmological hydrodynamical simulations to show that a significant fraction of the gas in high redshift rare massive halos falls nearly radially to their very centre on extremely short timescales. This process results in the formation of very compact bulges with specific angular momentum a factor 5-30$smaller than the average angular momentum of the baryons in the whole halo. Such low angular momentum originates both from segregation and effective cancellation when the gas flows to the centre of the halo along well defined cold filamentary streams. These filaments penetrate deep inside the halo and connect to the bulge from multiple rapidly changing directions. Structures falling in along the filaments (satellite galaxies) or formed by gravitational instabilities triggered by the inflow (star clusters) further reduce the angular momentum of the gas in the bulge. Finally, the fraction of gas radially falling to the centre appears to increase with the mass of the halo; we argue that this is most likely due to an enhanced cancellation of angular momentum in rarer halos which are fed by more isotropically distributed cold streams. Such an increasingly efficient funnelling of low-angular momentum gas to the centre of very massive halos at high redshift may account for the rapid pace at which the most massive supermassive black holes grow to reach observed masses around $10^9$M$_\odot$ at an epoch when the Universe is barely 1 Gyr old.

We provide a theoretical description of synchrotron fluctuations arising from magnetic turbulence. We derive an expression that relates the correlation of synchrotron fluctuations for an arbitrary index of relativistic electrons with the correlations of squared fluctuations of magnetic field component perpendicular to the line of sight. The latter correlations we study assuming that the turbulence is axisymmetric. We obtain general relations valid for an arbitrary model of magnetic turbulence and analyse the relations for particular example of magnetic turbulence that is supported by numerical simulations. We predict that the synchrotron intensity fluctuations are anisotropic with larger correlation present along the direction of magnetic field. This anisotropy is dominated by the quadrupole component with the ratio between quadrupole and monopole parts being sensitive to the compressibility of underlying turbulence. Our work opens avenues for quantitative studies of magnetic turbulence in our galaxy and beyond using synchrotron emission. It also outlines the directions of how synchrotron foreground emission can be separated from cosmological signal, i.e. from CMB or highly redshifted HI emission. For the sake of completeness we also provide the expressions for the synchrotron polarization (Stocks parameters and their combinations) for the model of axisymmetric magnetic turbulence.

The formalism to compute the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmological fields is developed. Leveraging the isotropy of the target statistics, the Gram-Charlier expansion is reformulated with rotation invariant variables. This formulation allows us to track the geometrical statistics of the cosmic field to all orders. It then allows us to connect the one point statistics of the critical sets to the growth factor through perturbation theory, which predicts the redshift evolution of higher order cumulants. In particular, the cosmic non-linear evolution of the skeleton's length, together with the statistics of extrema and Euler characteristic are investigated in turn. In 2D, the corresponding differential densities are analytic as a function of the excursion set threshold and the shape parameter. In 3D, the Euler characteristics and the field isosurface area are also analytic to all orders in the expansion. Numerical integrations are performed and simple fits are provided whenever closed form expressions are not available. These statistics are compared to estimates from N-body simulations and are shown to match well the cosmic evolution up to root mean square of the density field of ~0.2. In 3D, gravitational perturbation theory is implemented to predict the cosmic evolution of all the relevant Gram-Charlier coefficients for universes with scale invariant matter distribution. The one point statistics of critical sets could be used to constrain primordial non-Gaussianities and the dark energy equation of state on upcoming cosmic surveys; this is illustrated on idealized experiments.

State-of-the-art hydrodynamical simulations show that gas inflow through the virial sphere of dark matter halos is focused (i.e. has a preferred inflow direction), consistent (i.e. its orientation is steady in time) and amplified (i.e. the amplitude of its advected specific angular momentum increases with time). We explain this to be a consequence of the dynamics of the cosmic web within the neighbourhood of the halo, which produces steady, angular momentum rich, filamentary inflow of cold gas. On large scales, the dynamics within neighbouring patches drives matter out of the surrounding voids, into walls and filaments before it finally gets accreted onto virialised dark matter halos. As these walls/filaments constitute the boundaries of asymmetric voids, they acquire a net transverse motion, which explains the angular momentum rich nature of the later infall which comes from further away. We conjecture that this large-scale driven consistency explains why cold flows are so efficient at building up high redshift thin discs from the inside out.

In the context of the geometrical analysis of weakly non Gaussian CMB maps, the 2D differential extrema counts as functions of the excursion set threshold is derived from the full moments expansion of the joint probability distribution of an isotropic random field, its gradient and invariants of the Hessian. Analytic expressions for these counts are given to second order in the non Gaussian correction, while a Monte Carlo method to compute them to arbitrary order is presented. Matching count statistics to these estimators is illustrated on fiducial non-Gaussian "Planck" data.

The full moments expansion of the joint probability distribution of an isotropic random field, its gradient and invariants of the Hessian is presented in 2 and 3D. It allows for explicit expression for the Euler characteristic in ND and computation of extrema counts as functions of the excursion set threshold and the spectral parameter, as illustrated on model examples.

We report on two quantitative, morphological estimators of the filamentary structure of the Cosmic Web, the so-called global and local skeletons. The first, based on a global study of the matter density gradient flow, allows us to study the connectivity between a density peak and its surroundings, with direct relevance to the anisotropic accretion via cold flows on galactic halos. From the second, based on a local constraint equation involving the derivatives of the field, we can derive predictions for powerful statistics, such as the differential length and the relative saddle to extrema counts of the Cosmic web as a function of density threshold (with application to percolation of structures and connectivity), as well as a theoretical framework to study their cosmic evolution through the onset of gravity-induced non-linearities.

We propose a qualitatively new mechanism for generating cosmological fluctuations from inflation. The non-equilibrium excitation of interacting scalar fields often evolves into infra-red (IR) and ultra-violet (UV) cascading, resulting in an intermediate scaling regime. We observe elements of this phenomenon in a simple model with inflaton \phi and iso-inflaton \chi fields interacting during inflation via the coupling g^2 (\phi-\phi_0)^2 \chi^2. Iso-inflaton particles are created during inflation when they become instantaneously massless at \phi=\phi_0, with occupation numbers not exceeding unity. We point out that very quickly the produced \chi particles become heavy and their multiple re-scatterings off the homogeneous condensate \phi(t) generates bremschtrahlung radiation of light inflaton IR fluctuations with high occupation numbers. The subsequent evolution of these IR fluctuations is qualitatively similar to that of the usual inflationary fluctuations, but their initial amplitude is different. The IR cascading generates a bump-shaped contribution to the cosmological curvature fluctuations, which can even dominate over the usual fluctuations for g^2>0.06. The IR cascading curvature fluctuations are significantly non-gaussian and the strength and location of the bump are model-dependent, through g^2 and \phi_0. The effect from IR cascading fluctuations is significantly larger than that from the momentary slowing-down of \phi(t). With a sequence of such bursts of particle production, the superposition of the bumps can lead to a new broad band non-gaussian component of cosmological fluctuations added to the usual fluctuations. Such a sequence of particle creation events can, but need not, lead to trapped inflation.

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

We present an overview of the Velocity Coordinate Spectrum (VCS), a new technique for studying astrophysical turbulence that utilizes the line-of-sight statistics of Doppler-broadened spectral lines. We consider the retrieval of turbulence spectra from emission intensity observations of both high and low spatial resolution and find that the VCS allows one to study turbulence even when the emitting turbulent volume is not spatially resolved. This opens interesting prospects for using the technique for extragalactic research. VCS developed for spectral emission lines is applicable to absorption lines as well if the optical depth is used instead of intensity. VCS for absorption lines in point-source spectra benefit from effectively narrow beam and does not require dense sky coverage by sampling directions. Even strongly saturated absorption lines still carry the information about the small scale turbulence, albeit limited to the wings of a line. Combining different absorption lines one can develop tomography of the turbulence in the interstellar gas in all its complexity.

The local theory of the critical lines of 2D and 3D Gaussian fields that underline the cosmic structures is presented. In the context of cosmological matter distribution the subset of critical lines of the 3D density field serves to delineate the skeleton of the observed filamentary structure at large scales. A stiff approximation used to quantitatively describe the filamentary skeleton shows that the flux of the skeleton lines is related to the average Gaussian curvature of the 1D (2D) sections of the field, much in the same way as the density of the peaks. The distribution of the length of the critical lines with threshold is analyzed in detail, while the extended descriptors of the skeleton - its curvature and its singular points, are introduced and briefly described. Theoretical predictions are compared to measurements of the skeleton in realizations of Gaussian random fields in 2D and 3D. It is found that the stiff approximation predicts accurately the shape of the differential length, allows for analytical insight, and explicit closed form solutions. Finally, it provides a simple classification of the singular points of the critical lines: i) critical points; ii) bifurcation points; iii) slopping plateaux.

This technical paper describes a software package that was designed to produce initial conditions for large cosmological simulations in the context of the Horizon collaboration. These tools generalize E. Bertschinger's Grafic1 software to distributed parallel architectures and offer a flexible alternative to the Grafic2 software for ``zoom'' initial conditions, at the price of large cumulated cpu and memory usage. The codes have been validated up to resolutions of 4096^3 and were used to generate the initial conditions of large hydrodynamical and dark matter simulations. They also provide means to generate constrained realisations for the purpose of generating initial conditions compatible with, e.g. the local group, or the SDSS catalog.

The onset of stochasticity is measured in $\Lambda$CDM cosmological simulations using a set of classical observables. It is quantified as the local derivative of the logarithm of the dispersion of a given observable (within a set of different simulations differing weakly through their initial realization), with respect to the cosmic growth factor. In an Eulerian framework, it is shown here that chaos appears at small scales, where dynamic is non-linear, while it vanishes at larger scales, allowing the computation of a critical transition scale corresponding to ~ 3.5 Mpc/h. This picture is confirmed by Lagrangian measurements which show that the distribution of substructures within clusters is partially sensitive to initial conditions, with a critical mass upper bound scaling roughly like the perturbation's amplitude to the power 0.15. The corresponding characteristic mass, $M_{\rm crit}=2 10^{13} M_{\odot}$, is roughly of the order of the critical mass of non linearities at z=1 and accounts for the decoupling induced by the dark energy triggered acceleration. The sensitivity to detailed initial conditions spills to some of the overall physical properties of the host halo (spin and velocity dispersion tensor orientation) while other "global" properties are quite robust and show no chaos (mass, spin parameter, connexity and center of mass position). This apparent discrepancy may reflect the fact that quantities which are integrals over particles rapidly average out details of difference in orbits, while the other observables are more sensitive to the detailed environment of forming halos and reflect the non-linear scale coupling characterizing the environments of halos.

We continue our work on developing techniques for studying turbulence with spectroscopic data. We show that Doppler-broadened absorption spectral lines, in particularly, saturated absorption lines, can be used within the framework of the earlier-introduced technique termed the Velocity Coordinate spectrum (VCS). The VCS relates the statistics of fluctuations along the velocity coordinate to the statistics of turbulence, thus it does not require spatial coverage by sampling directions in the plane of the sky. We consider lines with different degree of absorption and show that for lines of optical depth less than one, our earlier treatment of the VCS developed for spectral emission lines is applicable, if the optical depth is used instead of intensity. This amounts to correlating the logarithms of absorbed intensities. For larger optical depths and saturated absorption lines, we show, that the amount of information that one can use is, inevitably, limited by noise. In practical terms, this means that only wings of the line are available for the analysis. In terms of the VCS formalism, this results in introducing an additional window, which size decreases with the increase of the optical depth. As a result, strongly saturated absorption lines carry the information only about the small scale turbulence. Nevertheless, the contrast of the fluctuations corresponding to the small scale turbulence increases with the increase of the optical depth, which provides advantages for studying turbulence combining lines with different optical depths. Combining different absorption lines one can tomography turbulence in the interstellar gas in all its complexity.

We analyze the BOOMERanG 2003 (B03) 145 GHz temperature map to constrain the amplitude of a non Gaussian, primordial contribution to CMB fluctuations. We perform a pixel space analysis restricted to a portion of the map chosen in view of high sensitivity, very low foreground contamination and tight control of systematic effects. We set up an estimator based on the three Minkowski functionals which relies on high quality simulated data, including non Gaussian CMB maps. We find good agreement with the Gaussian hypothesis and derive the first limits based on BOOMERanG data for the non linear coupling parameter f_NL as -300<f_NL<650 at 68% CL and -800<f_NL<1050 at 95% CL.

The skeleton formalism aims at extracting and quantifying the filamentary structure of the universe is generalized to 3D density fields; a numerical method for computating a local approximation of the skeleton is presented and validated here on Gaussian random fields. This method manages to trace well the filamentary structure in 3D fields such as given by numerical simulations of the dark matter distribution on large scales and is insensitive to monotonic biasing. Two of its characteristics, namely its length and differential length, are analyzed for Gaussian random fields. Its differential length per unit normalized density contrast scales like the PDF of the underlying density contrast times the total length times a quadratic Edgeworth correction involving the square of the spectral parameter. The total length scales like the inverse square smoothing length, with a scaling factor given by 0.21 (5.28+ n) where n is the power index of the underlying field. This dependency implies that the total length can be used to constrain the shape of the underlying power spectrum, hence the cosmology. Possible applications of the skeleton to galaxy formation and cosmology are discussed. As an illustration, the orientation of the spin of dark halos and the orientation of the flow near the skeleton is computed for dark matter simulations. The flow is laminar along the filaments, while spins of dark halos within 500 kpc of the skeleton are preferentially orthogonal to the direction of the flow at a level of 25%.

We discuss a new technique for studying astrophysical turbulence that utilizes the statistics of Doppler-broadened spectral lines. The technique relates the power Velocity Coordinate Spectrum (VCS), i.e. the spectrum of fluctuations measured along the velocity axis in Position-Position-Velocity (PPV) data cubes available from observations, to the underlying power spectra of the velocity/density fluctuations. Unlike the standard spatial spectra, that are function of angular wavenumber, the VCS is a function of the velocity wave number k_v ~ 1/v. We show that absorption affects the VCS to a higher degree for small k_v and obtain the criteria for disregarding the absorption effects for turbulence studies at large k_v. We consider the retrieval of turbulence spectra from observations for high and low spatial resolution observations and find that the VCS allows one to study turbulence even when the emitting turbulent volume is not spatially resolved. This opens interesting prospects for using the technique for extragalactic research. We show that, while thermal broadening interferes with the turbulence studies using the VCS, it is possible to separate thermal and non-thermal contributions. This allows a new way of determining the temperature of the interstellar gas using emission and absorption spectral lines.

We study the clustering properties of metals in the intergalactic medium (IGM) as traced by 619 CIV, 81 SiIV, N >= 10^12 cm^-2 and 316 MgII, and 82 FeII N >= 10^11.5 cm^-2 absorption components in 19 high signal-to-noise (60-100 per pixel), high resolution (R = 45000) quasar spectra. Over the redshift range probed (1.5-3.0), CIV and SiIV trace each other closely and their line-of-sight correlation functions exhibit a steep decline at large separations and a flatter profile below ~ 150 km s^-1, with a large overall bias. These features do not depend on column depth. Carrying out a detailed SPH simulation (2 X 320^3, 57 Mpc^3 comoving), we show that this behavior can not be reproduced by models in which the IGM metallicity is constant or a local function of density. However, the CIV correlation function is consistent with a model in which metals are confined within bubbles with a typical radius Rs = 2 comoving Mpc about sources of mass >= Ms = 10^12 solar masses at z=3. Our lower redshift (0.5-2) measurements of the MgII and FeII correlation functions also uncover a steep decline at large separations and a flatter profile at small separations, but the clustering is even higher, and the turn-over is shifted to ~ 75 km s^-1. Again these features do not change with column depth. We describe an analytical bubble model for these species, which come from regions that are too compact to be simulated numerically, deriving best-fit values of R_s ~ 2.4 Mpc and M_s ~ 10^12 solar masses. Equally good fits to all four species are found in a similarly biased high-redshift enrichment model in which metals are placed within 2.4 comoving Mpc of 3 x 10^9 solar mass sources at z = 7.5.

We present the cosmological parameters from the CMB intensity and polarization power spectra of the 2003 Antarctic flight of the BOOMERANG telescope. The BOOMERANG data alone constrains the parameters of the $\Lambda$CDM model remarkably well and is consistent with constraints from a multi-experiment combined CMB data set. We add LSS data from the 2dF and SDSS redshift surveys to the combined CMB data set and test several extensions to the standard model including: running of the spectral index, curvature, tensor modes, the effect of massive neutrinos, and an effective equation of state for dark energy. We also include an analysis of constraints to a model which allows a CDM isocurvature admixture.

We present the BOOMERanG-03 experiment and maps of the Stokes parameters I, Q, U of the microwave sky obtained during a 14 day balloon flight in 2003. Three regions of the southern sky were surveyed: a deep survey (~ 90 square degrees) and a shallow survey (~ 750 square degrees) at high Galactic latitudes (both centered at RA ~ 5.5 h, dec ~ -45 deg) and a survey of ~ 300 square degrees across the Galactic plane at RA ~ 9.1 h, dec ~ -47 deg. All three surveys were carried out in three wide frequency bands centered at 145, 245 and 345 GHz, with an angular resolution of ~ 10'. The 145 GHz maps of Stokes I are dominated by Cosmic Microwave Background (CMB) temperature anisotropy, which is mapped with high signal to noise ratio. The measured anisotropy pattern is consistent with the pattern measured in the same region by BOOMERanG-98 and by WMAP. The 145 GHz maps of Stokes Q and U provide a robust statistical detection of polarization of the CMB when subjected to a power spectrum analysis. This amplitude of the polarization is consistent with that of the CMB in the $\Lambda$CDM cosmological scenario. At 145 GHz, in the CMB surveys, the intensity and polarization of the astrophysical foregrounds are found to be negligible with respect to the cosmological signal. At 245 and 345 GHz we detect ISD emission correlated to the 3000 GHz IRAS/DIRBE maps, and give upper limits for any other non-CMB component. We also present intensity maps of the surveyed section of the Galactic plane. These are compared to monitors of different interstellar components, showing that a variety of emission mechanisms is present in that region.

We report measurements of the CMB polarization power spectra from the January 2003 Antarctic flight of BOOMERANG. The primary results come from six days of observation of a patch covering 0.22% of the sky centered near R.A. = 82.5 deg., Dec= -45 deg. The observations were made using four pairs of polarization sensitive bolometers operating in bands centered at 145 GHz. Using two independent analysis pipelines, we measure a non-zero <EE> signal in the range 100< l <1000 with a significance 4.8-sigma, a 2-sigma upper limit of 8.6 uK^2 for any <BB> contribution, and a 2-sigma upper limit of 7.0 uK^2 for the <EB> spectrum. Estimates of foreground intensity fluctuations and the non-detection of <BB> and <EB> signals rule out any significant contribution from galactic foregrounds. The results are consistent with a Lambda-CDM cosmology seeded by adiabatic perturbations. We note that this is the first detection of CMB polarization with bolometric detectors.

We present a measurement of the temperature-polarization angular cross power spectrum, <TE>, of the Cosmic Microwave Background. The result is based on $\sim 200$ hours of data from 8 polarization sensitive bolometers operating at 145 GHz during the 2003 flight of BOOMERANG. We detect a significant <TE> correlation in the $\ell$-range between 50 and 950 with a statistical significance > 3.5 sigma. Contamination by polarized foreground emission and systematic effects are negligible in comparison with statistical uncertainty. The spectrum is consistent with previous detections and with the "concordance model" that assumes adiabatic initial conditions. This is the first measurement of <TE> using bolometric detectors.

We report on observations of the Cosmic Microwave Background (CMB) obtained during the January 2003 flight of Boomerang . These results are derived from 195 hours of observation with four 145 GHz Polarization Sensitive Bolometer (PSB) pairs, identical in design to the four 143 GHz Planck HFI polarized pixels. The data include 75 hours of observations distributed over 1.84% of the sky with an additional 120 hours concentrated on the central portion of the field, itself representing 0.22% of the full sky. From these data we derive an estimate of the angular power spectrum of temperature fluctuations of the CMB in 24 bands over the multipole range (50 < l < 1500). A series of features, consistent with those expected from acoustic oscillations in the primordial photon-baryon fluid, are clearly evident in the power spectrum, as is the exponential damping of power on scales smaller than the photon mean free path at the epoch of last scattering (l > 900). As a consistency check, the collaboration has performed two fully independent analyses of the time ordered data, which are found to be in excellent agreement.

In this presentation we review the link between the statistics of intensity fluctuations in spectral line data cubes with underlying statistical properties of turbulence in the interstellar medium. Both the formalism of Velocity Channel Analysis for optically thin lines and its extension to the lines with self-absorption is described. We demonstrate that by observing optically thin lines from cold gas in sufficiently narrow (thin) velocity channels one may recover the scaling of the stochastic velocities from turbulent cascade, in particular, Kolmogorov velocities give $K^{-2.7}$ contribution to the intensity power spectrum. Synthetically increasing the channel thickness separates out the underlying density inhomogeneities of the gas. Effects of self absorption, on the other hand, retain the velocity signature even for integrated lines. As a result, intensity fluctuations tend to show universal but featureless scaling of the power $\propto K^{-3}$ over the range of scales.

We discuss the nature of the possible high-l excess in the Cosmic Microwave Background (CMB) anisotropy power spectrum observed by the Cosmic Background Imager (CBI). We probe the angular structure of the excess in the CBI deep fields and investigate whether it could be due to the scattering of CMB photons by hot electrons within clusters, the Sunyaev-Zeldovich (SZ) effect. We estimate the density fluctuation parameters for amplitude, sigma_8, and shape, Gamma, from CMB primary anisotropy data and other cosmological data. We use the results of two separate hydrodynamical codes for Lambda-CDM cosmologies, consistent with the allowed sigma_8 and Gamma values, to quantify the expected contribution from the SZ effect to the bandpowers of the CBI experiment and pass simulated SZ effect maps through our CBI analysis pipeline. The result is very sensitive to the value of sigma_8, and is roughly consistent with the observed power if sigma_8 ~ 1. We conclude that the CBI anomaly could be a result of the SZ effect for the class of Lambda-CDM concordance models if sigma_8 is in the upper range of values allowed by current CMB and Large Scale Structure (LSS) data.

Polarization observations of the cosmic microwave background with the Cosmic Background Imager from September 2002 to May 2004 provide a significant detection of the E-mode polarization and reveal an angular power spectrum of polarized emission showing peaks and valleys that are shifted in phase by half a cycle relative to those of the total intensity spectrum. This key agreement between the phase of the observed polarization spectrum and that predicted based on the total intensity spectrum provides support for the standard model of cosmology, in which dark matter and dark energy are the dominant constituents, the geometry is close to flat, and primordial density fluctuations are predominantly adiabatic with a matter power spectrum commensurate with inflationary cosmological models.

We review the currrent cosmic parameter determinations of relevance to inflation using the WMAP-1year, Boomerang, CBI, Acbar and other CMB data. The basic steps in the pipelines which determine the bandpowers from the raw data from which these estimations are made are summarized. We forecast how the precision is likely to improve with more years of WMAP in combination with future ground-based experiments and with Planck. We address whether the current data indicates strong breaking from uniform acceleration through the relatively small region of the inflaton potential that the CMB probes, manifest in the much-discussed running spectral index or in even more radical braking/breaking scenarios. Although some weak ``anomalies'' appear in the current data, the statistical case is not there. However increased precision, at the high multipole end and with polarization measurements, will significantly curtail current freedom.

Quantitative description of the statistics of intensity fluctuations within spectral line data cubes introduced in our earlier work is extended to the absorbing media. A possibility of extracting 3D velocity and density statistics from both integrated line intensity as well as from the individual channel maps is analyzed. We find that absorption enables the velocity effects to be seen even if the spectral line is integrated over frequencies. This regime that is frequently employed in observations is characterized by a non-trivial relation between the spectral index of velocities and the spectral index of intensity fluctuations. For instance when density is dominated by fluctuations at large scales, i.e. when correlations scale as r^{-\gamma}, \gamma<0, the intensity fluctuations exhibit a universal spectrum of fluctuations ~K^{-3} over a range of scales. When small scale fluctuations of density contain most of the energy, i.e. when correlations scale as r^{-\gamma}, \gamma>0, the resulting spectrum of the integrated lines depends on the scaling of the underlying density and scales as K^{-3+\gamma}. We show that if we take the spectral line slices that are sufficiently thin we recover our earlier results for thin slice data without absorption. As the result we extend the Velocity Channel Analysis (VCA) technique to optically thick lines enabling studies of turbulence in molecular clouds. In addition, the developed mathematical machinery enables a quantitative approach to solving other problems that involved statistical description of turbulence within emitting and absorbing gas.

Two years of microwave background observations with the Cosmic Background Imager (CBI) have been combined to give a sensitive, high resolution angular power spectrum over the range 400 < l < 3500. This power spectrum has been referenced to a more accurate overall calibration derived from WMAP. The data cover 90 deg^2 including three pointings targeted for deep observations. The uncertainty on the l > 2000 power previously seen with the CBI is reduced. Under the assumption that any signal in excess of the primary anisotropy is due to a secondary Sunyaev-Zeldovich anisotropy in distant galaxy clusters we use CBI, ACBAR, and BIMA data to place a constraint on the present-day rms mass fluctuation sigma_8. We present the results of a cosmological parameter analysis on the l < 2000 primary anisotropy data which show significant improvements in the parameters as compared to WMAP alone, and we explore the role of the small-scale cosmic microwave background data in breaking parameter degeneracies.

In the 1998-99 flight, BOOMERanG has produced maps of $\sim 4 %$ of the sky at high Galactic latitudes, at frequencies of 90, 150, 240 and 410 GHz, with resolution $\simgt 10'$. The faint structure of the Cosmic Microwave Background at horizon and sub-horizon scales is evident in these maps. These maps compare well to the maps recently obtained at lower frequencies by the WMAP experiment. Here we compare the amplitude and morphology of the structures observed in the two sets of maps. We also outline the polarization sensitive version of BOOMERanG, which was flown early this year to measure the linear polarization of the microwave sky at 150, 240 and 350 GHz.

Abbreviated: We highlight the remarkable evolution in the CMB power spectrum over the past few years, and in the cosmological parameters for minimal inflation models derived from it. Grand unified spectra (GUS) show pre-WMAP optimal bandpowers are in good agreement with each other and with the one-year WMAP results, which now dominate the L < 600 bands. GUS are used to determine calibrations, peak/dip locations and heights, and damping parameters. These CMB experiments significantly increased the case for accelerated expansion in the early universe (the inflationary paradigm) and at the current epoch (dark energy dominance) when they were combined with `prior' probabilities on the parameters. A minimal inflation parameter set is applied in the same way to the evolving data. Grid-based and and Monte Carlo Markov Chain methods are shown to give similar values, highly stable over time and for different prior choices, with the increasing precision best characterized by decreasing errors on uncorrelated parameter eigenmodes. After marginalizing over the other cosmic and experimental variables for a weak+LSS prior, the pre-WMAP data of Jan03 cf. the post-WMAP data of Mar03 give Omega_{tot} =1.03^{+0.05}_{-0.04} cf. 1.02^{+0.04}_{-0.03}. Adding the flat prior, n_s =0.95^{+0.07}_{-0.04} cf. 0.97^{+0.02}_{-0.02}, with < 2\sigma evidence for a log variation of n_s. The densities have concordance values. The dark energy pressure-to-density ratio is not well constrained by our weak+LSS prior, but adding SN1 gives w_Q < -0.7. We find \sigma_8 = 0.89^{+0.06}_{-0.07} cf. 0.86^{+0.04}_{-0.04}, implying a sizable SZ effect; the high L power suggest \sigma_8 \sim 0.94^{+0.08}_{-0.16} is needed to be SZ-compatible.

BOOMERANG is a balloon-borne telescope designed for long duration (LDB) flights around Antarctica. The second LDB Flight of BOOMERANG took place in January 2003. The primary goal of this flight was to measure the polarization of the CMB. The receiver uses polarization sensitive bolometers at 145 GHz. Polarizing grids provide polarization sensitivity at 245 and 345 GHz. We describe the BOOMERANG telescope noting changes made for 2003 LDB flight, and discuss some of the issues involved in the measurement of polarization with bolometers. Lastly, we report on the 2003 flight and provide an estimate of the expected results.

We review the first science results from the Arcminute Cosmology Bolometer Array Receiver (ACBAR); a multifrequency millimeter-wave receiver optimized for observations of the Cosmic Microwave Background (CMB) and the Sunyaev-Zel'dovich (SZ) effect in clusters of galaxies. ACBAR was installed on the 2 m Viper telescope at the South Pole in January 2001 and the results presented here incorporate data through July 2002. We present the power spectrum of the CMB at 150 GHz over the range \ell = 150 - 3000 measured by ACBAR as well as estimates for the values of the cosmological parameters within the context of Lambda-CDM models. We find that the inclusion of Omega_Lambda greatly improves the fit to the power spectrum. We also observe a slight excess of small-scale anisotropy at 150 GHz; if interpreted as power from the SZ effect of unresolved clusters, the measured signal is consistent with CBI and BIMA within the context of the SZ power spectrum models tested.

We report on the cosmological parameters derived from observations with the Cosmic Background Imager (CBI), covering 40 square degrees and the multipole range 300 < l < 3500. The angular scales probed by the CBI correspond to structures which cover the mass range from 10^14 to 10^17 M_sun, and the observations reveal, for the first time, the seeds that gave rise to clusters of galaxies. These unique, high-resolution observations also show damping in the power spectrum to l ~ 2000, which we interpret as due to the finite width of the photon-baryon decoupling region and the viscosity operating at decoupling. Because the observations extend to much higher l the CBI results provide information complementary to that probed by the Boomerang, DASI, Maxima, and VSA experiments. As the observations are pushed to higher multipoles no anomalies relative to standard models appear, and extremely good consistency is found between the cosmological parameters derived for the CBI observations over the range 610 < l < 2000 and observations at lower l [abridged].

We report measurements of anisotropy in the cosmic microwave background radiation over the multipole range l ~ 200 - 3500 with the Cosmic Background Imager based on deep observations of three fields. These results confirm the drop in power with increasing l first reported in earlier measurements with this instrument, and extend the observations of this decline in power out to l \~ 2000. The decline in power is consistent with the predicted damping of primary anisotropies. At larger multipoles, l = 2000 - 3500, the power is 3.1 sigma greater than standard models for intrinsic microwave background anisotropy in this multipole range, and 3.5 sigma greater than zero. This excess power is not consistent with expected levels of residual radio source contamination but, for sigma_8 >~ 1, is consistent with predicted levels due to a secondary Sunyaev-Zeldovich anisotropy. Further observations are necessary to confirm the level of this excess and, if confirmed, determine its origin.

Using the Cosmic Background Imager, a 13-element interferometer array operating in the 26-36 GHz frequency band, we have observed 40 sq deg of sky in three pairs of fields, each ~ 145 x 165 arcmin, using overlapping pointings (mosaicing). We present images and power spectra of the cosmic microwave background radiation in these mosaic fields. We remove ground radiation and other low-level contaminating signals by differencing matched observations of the fields in each pair. The primary foreground contamination is due to point sources (radio galaxies and quasars). We have subtracted the strongest sources from the data using higher-resolution measurements, and we have projected out the response to other sources of known position in the power-spectrum analysis. The images show features on scales ~ 6 - 15 arcmin, corresponding to masses ~ (5 - 80)*10^{14} Msun at the surface of last scattering, which are likely to be the seeds of clusters of galaxies. The power spectrum estimates have a resolution Delta-l = 200 and are consistent with earlier results in the multipole range l <~ 1000. The power spectrum is detected with high signal-to-noise ratio in the range 300 <~ l <~ 1700. For 1700 <~ l <~ 3000 the observations are consistent with the results from more sensitive CBI deep-field observations. The results agree with the extrapolation of cosmological models fitted to observations at lower l, and show the predicted drop at high l (the "damping tail").

In a previous work Lazarian and Pogosyan suggested a technique to extract velocity and density statistics, of interstellar turbulence, by means of analysing statistics of spectral line data cubes. In this paper we test that technique, by studying the effect of correlation between velocity and density fields, providing a systematic analysis of the uncertainties arising from the numerics, and exploring the effect of a linear shear. We make use of both compressible MHD simulations and synthetic data to emulate spectroscopic observations and test the technique. With the same synthetic spectroscopic data, we also studied anisotropies of the two point statistics and related those anisotropies with the magnetic field direction. This presents a new technique for magnetic field studies. The results show that the velocity and density spectral indices measured are consistent with the analytical predictions. We identified the dominant source of error with the limited number of data points along a given line of sight. We decrease this type of noise by increasing the number of points and by introducing Gaussian smoothing. We argue that in real observations the number of emitting elements is essentially infinite and that source of noise vanishes.

We present millimeter-wave observations of three extra-galactic and six Galactic sources in the Southern sky. Observations were made at 90, 150, 240 and 400 GHz with resolutions of 18, 10, 14 and 12 arcmin respectively during the 1998 Antarctic long duration balloon flight of BOOMERANG. Observations were also made with the SEST telescope, at 90 and 150 GHz with resolutions of 57 and 35 arcsec respectively. These observations can be used for calibrations of Cosmic Microwave Background (CMB) experiments as well as an understanding of the physical processes of the sources.

Full suite of the present day Cosmic Microwave background (CMB) data, when combined with weak prior information on the Hubble constant and the age of the Universe, or the Large-Scale structure, provides strong indication for a non-zero density of the vacuum-like dark energy in our universe. This result independently supports the conclusions from Supernovae Ia (SN1a) data. When the model parameter space is extended to allow for the range of the equation of state parameter w_Q for the dynamical field Q which gives rise to dark energy, the CMB data is found to give a weak upper bound w_Q < -0.4 at 95% CL, however combined with SN1a data it strongly favours w_Q < -0.8, consistent with Lambda-term like behaviour.

We report an investigation of cosmological parameters based on the measurements of anisotropy in the cosmic microwave background radiation (CMB) made by ACBAR. We use the ACBAR data in concert with other recent CMB measurements to derive Bayesian estimates of parameters in inflation-motivated adiabatic cold dark matter models. We apply a series of additional cosmological constraints on the shape and amplitude of the density power spectrum, the Hubble parameter and from supernovae to further refine our parameter estimates. Previous estimates of parameters are confirmed, with sensitive measurements of the power spectrum now ranging from \ell \sim 3 to 2800. Comparing individual best model fits, we find that the addition of \Omega_\Lambda as a parameter dramatically improves the fits. We also use the high-\ell data of ACBAR, along with similar data from CBI and BIMA, to investigate potential secondary anisotropies from the Sunyaev-Zeldovich effect. We show that the results from the three experiments are consistent under this interpretation, and use the data, combined and individually, to estimate \sigma_8 from the Sunyaev-Zeldovich component.

We report the most complete analysis to date of observations of the Cosmic Microwave Background (CMB) obtained during the 1998 flight of BOOMERANG. We use two quite different methods to determine the angular power spectrum of the CMB in 20 bands centered at l = 50 to 1000, applying them to 50% more data than has previously been analyzed. The power spectra produced by the two methods are in good agreement with each other, and constitute the most sensitive measurements to date over the range 300 < l < 1000. The increased precision of the power spectrum yields more precise determinations of several cosmological parameters than previous analyses of BOOMERANG data. The results continue to support an inflationary paradigm for the origin of the universe, being well fit by a 13.5 Gyr old, flat universe composed of approximately 5% baryonic matter, 30% cold dark matter, and 65% dark energy, with a scale invariant initial density perturbations.

We review the recently published results from the CBI's first season of observations. Angular power spectra of the CMB were obtained from deep integrations of 3 single fields covering a total of 3 deg^2 and 3 shallower surveys of overlapping (mosaiced) fields covering a total of 40 deg^2. The observations show a damping of the anisotropies at high-l as expected from the standard scenarios of recombination. We present parameter estimates obtained from the data and discuss the significance of an excess at l>2000 observed in the deep fields.

Boomerang, Maxima, DASI, CBI and VSA significantly increase the case for accelerated expansion in the early universe (the inflationary paradigm) and at the current epoch (dark energy dominance), especially when combined with data on high redshift supernovae (SN1) and large scale structure (LSS). There are ``7 pillars of Inflation'' that can be shown with the CMB probe, and at least 5, and possibly 6, of these have already been demonstrated in the CMB data: (1) a large scale gravitational potential; (2) acoustic peaks/dips; (3) damping due to shear viscosity; (4) a Gaussian (maximally random) distribution; (5) secondary anisotropies; (6) polarization. A 7th pillar, anisotropies induced by gravity wave quantum noise, could be too small. A minimal inflation parameter set, \omega_b,\omega_{cdm}, \Omega_{tot}, \Omega_Q,w_Q,n_s,\tau_C, \sigma_8}, is used to illustrate the power of the current data. We find the CMB+LSS+SN1 data give \Omega_{tot} =1.00^{+.07}_{-.03}, consistent with (non-baroque) inflation theory. Restricting to \Omega_{tot}=1, we find a nearly scale invariant spectrum, n_s =0.97^{+.08}_{-.05}. The CDM density, \Omega_{cdm}{\rm h}^2 =.12^{+.01}_{-.01}, and baryon density, \Omega_b {\rm h}^2 = >.022^{+.003}_{-.002}, are in the expected range. (The Big Bang nucleosynthesis estimate is 0.019\pm 0.002.) Substantial dark (unclustered) energy is inferred, \Omega_Q \approx 0.68 \pm 0.05, and CMB+LSS \Omega_Q values are compatible with the independent SN1 estimates. The dark energy equation of state, crudely parameterized by a quintessence-field pressure-to-density ratio w_Q, is not well determined by CMB+LSS (w_Q < -0.4 at 95% CL), but when combined with SN1 the resulting w_Q < -0.7 limit is quite consistent with the w_Q=-1 cosmological constant case.

We derive cosmological parameters from the CBI measurements of the Cosmic Microwave Background (CMB) angular power spectrum. Our results provide an independent confirmation of the standard $\Omega_{\rm tot}=1$ $\Lambda$CDM model within the adiabatic, inflationary paradigm. Above $\ell=2000$ the observations show evidence of power in excess of that expected in the standard models. We use hydrodynamical simulations to show how Sunyaev-Zeldovich Effect (SZE) may account for the excess power for models with fluctuation amplitude $\sigma_8\sim 1$ which is in the high end of the range allowed by the primary CMB observations.

We describe an algorithm for the extraction of the angular power spectrum of an intensity field, such as the cosmic microwave background (CMB), from interferometer data. This new method, based on the gridding of interferometer visibilities in the aperture plane followed by a maximum likelihood solution for bandpowers, is much faster than direct likelihood analysis of the visibilities, and deals with foreground radio sources, multiple pointings, and differencing. The gridded aperture-plane estimators are also used to construct Wiener-filtered images using the signal and noise covariance matrices used in the likelihood analysis. Results are shown for simulated data. The method has been used to determine the power spectrum of the cosmic microwave background from observations with the Cosmic Background Imager, and the results are given in companion papers.

We search the BOOMERanG maps of the anisotropy of the Cosmic Microwave Background (CMB) for deviations from gaussianity. In this paper we focus on analysis techniques in pixel-space, and compute skewness, kurtosis and Minkowski functionals for the BOOMERanG maps and for gaussian simulations of the CMB sky. We do not find any significant deviation from gaussianity in the high galactic latitude section of the 150 GHz map. We do find deviations from gaussianity at lower latitudes and at 410 GHz, and we ascribe them to Galactic dust contamination. Using non-gaussian simulations of instrumental systematic effects, of foregrounds, and of sample non-gaussian cosmological models, we set upper limits to the non-gaussian component of the temperature field in the BOOMERanG maps. For fluctuations distributed as a 1 DOF $\chi^2$ mixed to the main gaussian component our upper limits are in the few % range.

This paper presents a measurement of the angular power spectrum of the Cosmic Microwave Background from l=75 to l=1025 (~10' to 5 degrees) from a combined analysis of four 150 GHz channels in the BOOMERANG experiment. The spectrum contains multiple peaks and minima, as predicted by standard adiabatic-inflationary models in which the primordial plasma undergoes acoustic oscillations. These results significantly constrain the values of Omega_tot, Omega_b h^2, Omega_c h^2 and n_s.

We describe the BOOMERanG experiment and its main result, i.e. the measurement of the large scale curvature of the Universe. BOOMERanG is a balloon-borne microwave telescope with sensitive cryogenic detectors. BOOMERanG has measured the angular distribution of the Cosmic Microwave Background on $\sim 3%$ of the sky, with a resolution of $\sim 10$ arcmin and a sensitivity of $\sim 20 \mu K$ per pixel. The resulting image is dominated by hot and cold spots with rms fluctuations $\sim 80 \mu K$ and typical size of $\sim 1^o$. The detailed angular power spectrum of the image features three peaks and two dips at $\ell = (213^{+10}_{-13}), (541^{+20}_{-32}), (845^{+12}_{-25})$ and $\ell = (416^{+22}_{-12}), (750^{+20}_{-750})$, respectively. Such very characteristic spectrum can be explained assuming that the detected structures are the result of acoustic oscillations in the primeval plasma. In this framework, the measured pattern constrains the density parameter $\Omega$ to be $0.85 < \Omega < 1.1$ (95% confidence interval). Other cosmological parameters, like the spectral index of initial density fluctuations, the density parameter for baryons, dark matter and dark energy, are detected or constrained by the BOOMERanG measurements and by other recent CMB anisotropy experiments. When combined with other cosmological observations, these results depict a new, consistent, cosmological scenario.

HI data cubes are sources of unique information on interstellar turbulence. Doppler shifts due to supersonic motions contain information on turbulent velocity field which is otherwise difficult to obtain. However, the problem of separation of velocity and density fluctuations within HI data cubes is far from being trivial. Analytical description of the emissivity statistics of channel maps (velocity slices) in Lazarian & Pogosyan (2000) showed that the relative contribution of the density and velocity fluctuations depends on the thickness of the velocity slice. In particular, power-law assymptotics of the emissivity fluctuations change when the dispersion of the velocity at the scale under study becomes of the order of the velocity slice thickness (integrated width of the channel map). These results are the foundations of the Velocity-Channel Analysis (VCA) technique which allows to determine velocity and density statistics using 21-cm data cubes. The VCA has been successfully tested using data cubes obtained via compressible magnetohydrodynamic simulations and applied to Galactic and Magellanic Clouds data. As a tool it has become much more sophisticated recently when effects of absorption were accounted for. The systematic studies of vast 21-cm data sets to correlate the variations in the turbulence statistics with the astrophysical activity is on the agenda. This should allow to determine the interstellar energy injection mechanisms. Going beyond the VCA, we discuss other tools, namely, genus and anisotropy analysis. The first characterises the topology of HI, while the second provides magnetic field directions. We show a few applications of these new tools to HI data and MHD simulations.

Three peaks and two dips have been detected in the power spectrum of the cosmic microwave background from the BOOMERANG experiment, at $\ell \sim 210, 540, 840$ and $\ell \sim 420, 750$, respectively. Using model-independent analyses, we find that all five features are statistically significant and we measure their location and amplitude. These are consistent with the adiabatic inflationary model. We also calculate the mean and variance of the peak and dip locations and amplitudes in a large 7-dimensional parameter space of such models, which gives good agreement with the model-independent estimates, and forecast where the next few peaks and dips should be found if the basic paradigm is correct. We test the robustness of our results by comparing Bayesian marginalization techniques on this space with likelihood maximization techniques applied to a second 7-dimensional cosmological parameter space, using an independent computational pipeline, and find excellent agreement: $\Omega_{\rm tot} = 1.02^{+0.06}_{-0.05}$ {\it vs.} $1.04 \pm 0.05$, $\Omega_b h^2 = 0.022^{+0.004}_{-0.003}$ {\it vs.} $0.019^{+0.005}_{-0.004}$, and $n_s = 0.96^{+0.10}_{-0.09}$ {\it vs.} $0.90 \pm 0.08$. The deviation in primordial spectral index $n_s$ is a consequence of the strong correlation with the optical depth.

Recent results from BOOMERANG-98 and MAXIMA-1, taken together with COBE-DMR, provide consistent and high signal-to-noise measurements of the CMB power spectrum at spherical harmonic multipole bands over $2<\ell\lta800$. Analysis of the combined data yields 68% (95%) confidence limits on the total density, $\Omega_{\rm {tot}}\simeq 1.11 \pm 0.07 (^{+0.13}_{-0.12})$, the baryon density, $\Omega_b h^2\simeq 0.032^{+0.005}_{-0.004} (^{+0.009}_{-0.008})$, and the scalar spectral tilt, $n_s\simeq1.01^{+0.09}_{-0.07} (^{+0.17}_{-0.14})$. These data are consistent with inflationary initial conditions for structure formation. Taken together with other cosmological observations, they imply the existence of both non-baryonic dark matter and dark energy in the universe.

In this paper we test the results of a recent analytical study by Lazarian and Pogosyan, on the statistics of emissivity in velocity channel maps, in the case of realistic density and velocity fields obtained from numerical simulations of MHD turbulence in the interstellar medium (ISM). To compensate for the lack of well-developed inertial ranges in the simulations due to the limited resolution, we apply a procedure for modifying the spectral slopes of the fields while preserving the spatial structures. We find that the density and velocity are moderately correlated in space and prove that the analytical results by Lazarian and Pogosyan hold in the case when these fields obey the fluid conservation equations. Our results imply that the spectra of velocity and density can be safely recovered from the position-position-velocity (PPV) data cubes available through observations, and confirm that the relative contributions of the velocity and density fluctuations to those of the emissivity depend on the velocity resolution used and on the steepness of the density spectral index. Furthermore, this paper supports previous reports that an interpretation of the features in the PPV data cubes as simple density enhancements (i.e., ``clouds'') can be often erroneous, as we observe that changes in the velocity statistics substantially modify the statistics of emissivity within the velocity data cubes.

We describe here an iterative method for jointly estimating the noise power spectrum from a CMB experiment's time-ordered data, together with the maximum-likelihood map. We test the robustness of this method on simulated Boomerang datasets with realistic noise.

We report the first results from the BOOMERanG experiment, which mapped at 90, 150, 240 and 410 GHz a wide (3%) region of the microwave sky with minimal local contamination. From the data of the best 150 GHz detector we find evidence for a well defined peak in the power spectrum of temperature fluctuations of the Cosmic Microwave Background, localized at $\ell = 197 \pm 6$, with an amplitude of $(68 \pm 8) \mu K_{CMB}$. The location, width and amplitude of the peak is suggestive of acoustic oscillations in the primeval plasma. In the framework of inflationary adiabatic cosmological models the measured spectrum allows a Bayesian estimate of the curvature of the Universe and of other cosmological parameters. With reasonable priors we find $\Omega = (1.07 \pm 0.06)$ and $n_s = (1.00 \pm 0.08)$ (68%C.L.) in excellent agreement with the expectations from the simplest inflationary theories. We also discuss the limits on the density of baryons, of cold dark matter and on the cosmological constant.

BOOMERanG has recently resolved structures on the last scattering surface at redshift $\sim$ 1100 with high signal to noise ratio. We review the technical advances which made this possible, and we focus on the current results for maps and power spectra, with special attention to the determination of the total mass-energy density in the Universe and of other cosmological parameters.

We show how estimates of parameters characterizing inflation-based theories of structure formation localized over the past year when large scale structure (LSS) information from galaxy and cluster surveys was combined with the rapidly developing cosmic microwave background (CMB) data, especially from the recent Boomerang and Maxima balloon experiments. All current CMB data plus a relatively weak prior probability on the Hubble constant, age and LSS points to little mean curvature (Omega_{tot} = 1.08\pm 0.06) and nearly scale invariant initial fluctuations (n_s =1.03\pm 0.08), both predictions of (non-baroque) inflation theory. We emphasize the role that degeneracy among parameters in the L_{pk} = 212\pm 7 position of the (first acoustic) peak plays in defining the $\Omega_{tot}$ range upon marginalization over other variables. Though the CDM density is in the expected range (\Omega_{cdm}h^2=0.17\pm 0.02), the baryon density Omega_bh^2=0.030\pm 0.005 is somewhat above the independent 0.019\pm 0.002 nucleosynthesis estimate. CMB+LSS gives independent evidence for dark energy (Omega_\Lambda=0.66\pm 0.06) at the same level as from supernova (SN1) observations, with a phenomenological quintessence equation of state limited by SN1+CMB+LSS to w_Q<-0.7 cf. the w_Q=-1 cosmological constant case.

We describe the implications of cosmic microwave background (CMB) observations and galaxy and cluster surveys of large scale structure (LSS) for theories of cosmic structure formation, especially emphasizing the recent Boomerang and Maxima CMB balloon experiments. The inflation-based cosmic structure formation paradigm we have been operating with for two decades has never been in better shape. Here we primarily focus on a simplified inflation parameter set, {omega_b,omega_{cdm},Omega_{tot}, Omega_\Lambda,n_s,\tau_C, \sigma_8}. Combining all of the current CMB+LSS data points to the remarkable conclusion that the local Hubble patch we can access has little mean curvature (Omega_{tot}=1.08\pm 0.06) and the initial fluctuations were nearly scale invariant (n_s=1.03\pm 0.08), both predictions of (non-baroque) inflation theory. The baryon density is found to be slightly larger than that preferred by independent Big Bang Nucleosynthesis estimates (omega_b=0.030\pm 0.005 cf. 0.019\pm 0.002). The CDM density is in the expected range (omega_{cdm}=0.17 \pm 0.02). Even stranger is the CMB+LSS evidence that the density of the universe is dominated by unclustered energy akin to the cosmological constant (Omega_\Lambda=0.66\pm 0.06), at the same level as that inferred from high redshift supernova observations. We also sketch the CMB+LSS implications for massive neutrinos.

The past, present and future of cosmic microwave background (CMB) anisotropy research is discussed, with emphasis on the Boomerang and Maxima balloon experiments. These data are combined with large scale structure (LSS) information and high redshift supernova (SN1) observations to explore the inflation-based cosmic structure formation paradigm. Here we primarily focus on a simplified inflation parameter set, {omega_b,omega_{cdm},Omega_{tot}, Omega_Q,w_Q, n_s,tau_C, sigma_8}. After marginalizing over the other cosmic and experimental variables, we find the current CMB+LSS+SN1 data gives Omega_{tot}=1.04\pm 0.05, consistent with (non-baroque) inflation theory. Restricting to Omega_{tot}=1, we find a nearly scale invariant spectrum, n_s =1.03 \pm 0.07. The CDM density, omega_{cdm}=0.17\pm 0.02, is in the expected range, but the baryon density, omega_b=0.030\pm 0.004, is slightly larger than the current nucleosynthesis estimate. Substantial dark energy is inferred, Omega_Q\approx 0.68\pm 0.05, and CMB+LSS Omega_Q values are compatible with the independent SN1 estimates. The dark energy equation of state, parameterized by a quintessence-field pressure-to-density ratio w_Q, is not well determined by CMB+LSS (w_Q<-0.3 at 95%CL), but when combined with SN1 the resulting w_Q<-0.7 limit is quite consistent with the w_Q=-1 cosmological constant case. Though forecasts of statistical errors on parameters for current and future experiments are rosy, rooting out systematic errors will define the true progress.

We present a numerical study of topological descriptors of initially Gaussian and scale-free density perturbations evolving via gravitational instability in an expanding universe. We carefully evaluate and avoid numerical contamination in making accurate measurements on simulated fields on a grid in a finite box. Independent of extent of non linearity, the measured Euler number of the excursion set at the percolation threshold, $\delta_c$, is positive and nearly equal to the number of isolated components, suggesting that these structures are trees. Our study of critical point counts reconciles the clumpy appearance of the density field at $\delta_c$ with measured filamentary local curvature. In the Gaussian limit, we measure $|\delta_c|> \sigma$ in contrast to widely held belief that $|\delta_c| \sim \sigma$, where $\sigma^2$ is the variance of the density field.

We propose and implement a fast, universally applicable method for extracting the angular power spectrum C_l from CMB temperature maps by first estimating the correlation function \xi(\theta). Our procedure recovers the C_l's using N^2 (but potentially N logN), operations, where N is the number of pixels. This is in contrast with standard maximum likelihood techniques which require N^3 operations. Our method makes no special assumptions about the map, unlike present fast techniques which rely on symmetries of the underlying noise matrix, sky coverage, scanning strategy, and geometry. This enables for the first time the analysis of megapixel maps without symmetries. The key element of our technique is the accurate multipole decomposition of \xi(\theta). The C_l error bars and cross-correlations are found by a Monte-Carlo approach. We applied our technique to a large number of simulated maps with Boomerang sky coverage in 81000 pixels. We used a diagonal noise matrix, with approximately the same amplitude as Boomerang. These studies demonstrate that our technique provides an unbiased estimator of the C_l's. Even though our method is approximate, the error bars obtained are nearly optimal, and converged only after few tens of Monte-Carlo realizations. Our method is directly applicable for the non-diagonal noise matrix. This, and other generalizations, such as minimum variance weighting schemes, polarization, and higher order statistics are also discussed.

The anisotropy of the cosmic microwave background radiation contains information about the contents and history of the universe. We report new limits on cosmological parameters derived from the angular power spectrum measured in the first Antarctic flight of the BOOMERANG experiment. Within the framework of inflation-motivated adiabatic cold dark matter models, and using only weakly restrictive prior probabilites on the age of the universe and the Hubble expansion parameter $h$, we find that the curvature is consistent with flat and that the primordial fluctuation spectrum is consistent with scale invariant, in agreement with the basic inflation paradigm. We find that the data prefer a baryon density $\Omega_b h^2$ above, though similar to, the estimates from light element abundances and big bang nucleosynthesis. When combined with large scale structure observations, the BOOMERANG data provide clear detections of both dark matter and dark energy contributions to the total energy density $\Omega_{\rm {tot}}$, independent of data from high redshift supernovae.

The blackbody radiation left over from the Big Bang has been transformed by the expansion of the Universe into the nearly isotropic 2.73K Cosmic Microwave Background. Tiny inhomogeneities in the early Universe left their imprint on the microwave background in the form of small anisotropies in its temperature. These anisotropies contain information about basic cosmological parameters, particularly the total energy density and curvature of the universe. Here we report the first images of resolved structure in the microwave background anisotropies over a significant part of the sky. Maps at four frequencies clearly distinguish the microwave background from foreground emission. We compute the angular power spectrum of the microwave background, and find a peak at Legendre multipole $\ell_{peak}=(197 \pm 6)$, with an amplitude $DT_{200}=(69 \pm 8)\mu K$. This is consistent with that expected for cold dark matter models in a flat (euclidean) Universe, as favoured by standard inflationary scenarios.

The distribution of atomic hydrogen in the Galactic plane is usually mapped using the Doppler shift of 21cm emission line. We calculate the emission spectrum in velocity slices of data (channel maps) and derive its dependence on the statistics of HI velocity and density fields. We find that (a) if the density spectrum is steep, i.e. n<-3, the short-wave asymptotics of the emissivity spectrum is dominated by velocity fluctuations; (b) the velocity fluctuations make the emission spectra shallower, provided that the data slices are sufficiently thin. In other words, turbulent velocity creates small scale structure that can erroneously be identified as clouds. The contribution of fluctuations in warm HI is suppressed relative to cold component when velocity channels used are narrower than warm HI thermal velocity and small angular scale fluctuations are measured. We calculate how emission spectra vary with the change of velocity slice thickness and show that the observational 21cm data is consistent with the explanation that intensity fluctuations within individual channel maps are generated by a turbulent velocity field. As the thickness of velocity slices increases density fluctuations get to dominate the emissivity. This allows to disentangle velocity and density statistics. Application of our technique to the Galactic and SMC data reveals spectra of density and velocity with the power law index close to -11/3. This is a Kolmogorov index, but the explanation of the spectrum appealing to the Kolmogorov-type cascade faces substantial difficulties. We generalize our treatment for the case of a statistical study of turbulence inside individual clouds. The mathematical machinery developed is applicable to other emission lines.

CMB anisotropy measurements have brought the issue of global topology of the universe from the realm of theoretical possibility to within the grasp of observations. The global topology of the universe modifies the correlation properties of cosmic fields. In particular, strong correlations are predicted in CMB anisotropy patterns on the largest observable scales if the size of the Universe is comparable to the distance to the CMB last scattering surface. We describe in detail our completely general scheme using a regularized method of images for calculating such correlation functions in models with nontrivial topology, and apply it to the computationally challenging compact hyperbolic spaces. Our procedure directly sums over images within a specified radius, ideally many times the diameter of the space, effectively treats more distant images in a continuous approximation, and uses Cesaro resummation to further sharpen the results. At all levels of approximation the symmetries of the space are preserved in the correlation function. This new technique eliminates the need for the difficult task of spatial eigenmode decomposition on these spaces. Although the eigenspectrum can be obtained by this method if desired, at a given level of approximation the correlation functions are more accurately determined. We use the 3-torus example to demonstrate that the method works very well. We apply it to power spectrum as well as correlation function evaluations in a number of compact hyperbolic (CH) spaces. Application to the computation of CMB anisotropy correlations on CH spaces, and the observational constraints following from them, are given in a companion paper.

We calculate the CMB anisotropy in compact hyperbolic universe models using the regularized method of images described in paper-I, including the 'line-of-sight `integrated Sachs-Wolfe' effect, as well as the last-scattering surface terms. We calculate the Bayesian probabilities for a selection of models by confronting our theoretical pixel-pixel temperature correlation functions with the COBE-DMR data. Our results demonstrate that strong constraints on compactness arise: if the universe is small compared to the `horizon' size, correlations appear in the maps that are irreconcilable with the observations. This conclusion is qualitatively insensitive to the matter content of the universe, in particular, the presence of a cosmological constant. If the universe is of comparable size to the 'horizon', the likelihood function is very dependent upon orientation of the manifold wrt the sky. While most orientations may be strongly ruled out, it sometimes happens that for a specific orientation the predicted correlation patterns are preferred over those for the conventional infinite models. The full Bayesian analysis we use is the most complete statistical test that can be done on the COBE maps, taking into account all possible signals and their variances in the theoretical skies, in particular the high degree of anisotropic correlation that can exist. We show that standard visual measures for comparing theoretical predictions with the data such as the isotropized power spectrum $C_\ell$ are not so useful in small compact spaces because of enhanced cosmic variance associated with the breakdown of statistical isotropy.

We review the main theoretical aspects of the structure formation paradigm which impinge upon wide angle surveys: the early universe generation of gravitational metric fluctuations from quantum noise in scalar inflaton fields; the well understood and computed linear regime of CMB anisotropy and large scale structure (LSS) generation; the weakly nonlinear regime, where higher order perturbation theory works well, and where the cosmic web picture operates, describing an interconnected LSS of clusters bridged by filaments, with membranes as the intrafilament webbing. Current CMB+LSS data favour the simplest inflation-based $\Lambda$CDM models, with a primordial spectral index within about 5% of scale invariant and $\Omega_\Lambda \approx 2/3$, similar to that inferred from SNIa observations, and with open CDM models strongly disfavoured. The attack on the nonlinear regime with a variety of N-body and gas codes is described, as are the excursion set and peak-patch semianalytic approaches to object collapse. The ingredients are mixed together in an illustrative gasdynamical simulation of dense supercluster formation.

Simple analytic arguments are used to understand the predominantly filamentary web in the large-scale distribution of galaxies. Numerical simulations of superclusters are performed to assess the feasibility of directly mapping the intracluster webbing, in particular with weak gravitational lensing using Wide Field Cameras.

The measurements of CMB anisotropy have opened up a window for probing the global topology of the universe on length scales comparable to and beyond the Hubble radius. For compact topologies, the two main effects on the CMB are: (1) the breaking of statistical isotropy in characteristic patterns determined by the photon geodesic structure of the manifold and (2) an infrared cutoff in the power spectrum of perturbations imposed by the finite spatial extent. We present a completely general scheme using the regularized method of images for calculating CMB anisotropy in models with nontrivial topology, and apply it to the computationally challenging compact hyperbolic topologies. This new technique eliminates the need for the difficult task of spatial eigenmode decomposition on these spaces. We estimate a Bayesian probability for a selection of models by confronting the theoretical pixel-pixel temperature correlation function with the COBE-DMR data. Our results demonstrate that strong constraints on compactness arise: if the universe is small compared to the `horizon' size, correlations appear in the maps that are irreconcilable with the observations. If the universe is of comparable size, the likelihood function is very dependent upon orientation of the manifold wrt the sky. While most orientations may be strongly ruled out, it sometimes happens that for a specific orientation the predicted correlation patterns are preferred over the conventional infinite models.

The measurements of CMB anisotropy have opened up a window for probing the global topology of the universe on length scales comparable to and beyond the Hubble radius. We have developed a new method for calculating the CMB anisotropy in models with nontrivial topology and apply it to open universe models with compact spatial topology. We conduct a Bayesian probability analysis for a selection of models which confronts the theoretical pixel-pixel temperature correlation function with the COBE-DMR data. Our results demonstrate that strong constraints on compactness arise: if the universe is small compared to the `horizon' size, correlations appear in the maps that are irreconcilable with the observations.

The (large angle) COBE DMR data can be used to probe the global topology of our universe on scales comparable to and just beyond the present ``horizon''. For compact topologies, the two main effects on the CMB are: [1] the breaking of statistical isotropy in characteristic patterns determined by the photon geodesic structure of the manifold and [2] an infrared cutoff in the power spectrum of perturbations imposed by the finite spatial extent. To make a detailed confrontation of these effects with the COBE maps requires the computation of the pixel-pixel temperature correlation function for each topology and for each orientation of it relative to the sky. We present a general technique using the method of images for doing this in compact hyperbolic (CH) topologies which does not require spatial eigenmode decomposition. We demonstrate that strong constraints on compactness follow from [2] and that these limits can be improved by exploiting the details of the geodesic structure for each individual topology ([1]), as we show for the flat 3-torus and selected CH models.

This paper presents a statistical explanation of filament formation in the galactic atomic hydrogen. Recently developed technique allows to determine the 3D spectrum of random HI density. We claim that even in the absence of dynamical factors the Gaussian field corresponding to the measured values of the spectrum of HI density should exhibit filamentary structure, the existence of which has long been claimed.

Observations indicate galaxies are distributed in a filament-dominated web-like structure. Numerical experiments at high and low redshift of viable structure formation theories also show filament-dominance. We present a simple quantitative explanation of why this is so, showing that the final-state web is actually present in embryonic form in the overdensity pattern of the initial fluctuations, with nonlinear dynamics just sharpening the image. The web is largely defined by the position and primordial tidal fields of rare events in the medium, with the strongest filaments between nearby clusters whose tidal tensors are nearly aligned. Applications of the cosmic web theory to observations include probing cluster-cluster bridges by weak gravitational lensing, X-rays, and the Sunyaev-Zeldovich effect and probing high redshift galaxy-galaxy bridges by low column density Lyman alpha absorption lines.

We discuss the structure of clusters in a class of flat cosmological models with the fraction of mass \Omega_{CDM} ~0.8 in cold dark matter, and the rest in hot dark matter in the form of massive neutrinos. We consider such Cold+Hot Dark Matter (CHDM) models with one, two or three massive neutrinos, with total mass ~4.6eV. Neutrinos of such low mass cannot constitute halos of galaxies and groups, but only of clusters of galaxies. The limit on the density of neutrinos in the central parts of galaxy clusters is estimated from the phase space den- sity constraints. The ratio of the neutrino density to that of CDM through the cluster is found analytically. It appears that the density of neutrinos is suppressed within the Abell radius. However, neutrinos contribute ~20% of the mass density to the cluster halo. Our numerical simulations match analytical results. The simulations indicate that the cluster halo dark matter density profile has the power-law slope ~-2.5 which is close to that in the model with cosmological constant. We also found that in the CHDM models the velocity dispersion is almost constant across the cluster. This is quite different from the model with cosmological constant or the open model where the velocity dispersion falls in the cluster outskirts. We discuss X-ray emission and weak gravitational lensing by clusters in the model. We input the found spherically symmetrical fit to the CHDM mass density profile and the X-ray surface brightness for the cluster A2256 into simple equation of hydrostatic equilibrium of the hot gas. X-ray temperature derived this way departures from both the data and actual prediction of the model, which give almost constant temperature. We found also that the problem of high baryonic fraction in clusters is not resolved in the CHDM models.

Mixed cold-hot dark matter cosmological models (CHDM) with $\Omega_{tot}=1$, approximately flat initial spectrum of adiabatic perturbations and 1, 2 or 3 types of massive neutrinos are compared and tested using recent observational data. The models with 2 or 3 neutrino types of equal mass permit as the best fit larger values of both the Hubble constant ($H_0\le 60$ for 2 types, $H_0\le 65$ for 3 types) and the total $\Omega_{\nu}$ (up to 0.3 for 3 types) than the model with 1 massive type. Also, they have less problems with abundances of early compact objects including $Ly-\alpha$ clouds.

The mixed cold-hot dark matter cosmological model (CHDM) with $\Omega_{tot}=1$ and a falling power-law initial spectrum of Gaussian adiabatic perturbations ($n>1$) is tested using recent obserbational data. It is shown that its fit to the data becomes worse with the growth of $n-1$, and may be considered as unreasonable for $n>1.1$ for all possible values of the Hubble constant. Thus, the CHDM model with a falling initial spectrum is worse than the same model with the approximately flat $(|n-1|<0.1)$ spectrum. On the other hand, the CHDM model provides a rather good fit to the data if $n$ lies in the range $(0.9-1.0)$, the Hubble constant $H_0 < 60$ km/s/Mpc ($H_0 < 55$ for $n=1$) and the neutrino energy density $\Omega_{\nu}< 0.25$. So, the CHDM model provides the best possibility for the realization of the simplest variants of the inflationary scenario having the effective slope $n\approx (0.95-0.97)$ between galaxy and horizon scales, including a modest contribution of primordial gravitational wave background to large-angle $\Delta T/T$ fluctuations of the cosmic microwave background (resulting in the increase of their total {\it rms} amplitude by $(5-10)%$) expected in some variants. A classification of cosmological models according to the number of fundamental parameters used to fit observational data is presented, too.

We study the development of gravitational instability in the strongly non-linear regime. For this purpose we use a number of statistical indicators such as filamentary statistics, spectrum of overdense/underdense regions and the void probability function, each of which probes a particular aspect of gravitational clustering. We use these statistical indicators to discriminate between different approximations to gravitational instability which we test against N-body simulations. The approximations which we test are, the truncated Zel'dovich approximation (TZ), the adhesion model (AM), and the frozen flow (FF) and linear potential (LP) approximations. Of these we find that FF and LP break down relatively early, soon after the non-linear length scale exceeds $R_*$ -- the mean distance between peaks of the gravitational potential. The reason for this break down is easy to understand, particles in FF are constrained to follow the streamlines of the initial velocity field. Shell crossing is absent in this case and structure gradually freezes as particles begin to collect near minima of the gravitational potential. In LP particles follow the lines of force of the primordial potential, oscillating about its minima at late times when the non-linear length scale $k_{\rm NL}^{-1}\simeq R_*$. Unlike FF and LP the adhesion model (and to some extent TZ) continues to give accurate results even at late times when $k_{\rm NL}^{-1} \ge R_*$. This is because both AM and TZ use