Velocity and density field reconstructions of the volume of the universe within 0.05c derived from the Cosmicflows-3 catalog of galaxy distances has revealed the presence of a filamentary structure extending across ~ 0.11c. The structure, at a characteristic redshift of 12,000 km/s, has a density peak coincident with the celestial South Pole. This structure, the largest contiguous feature in the local volume and comparable to the Sloan Great Wall at half the distance, is given the name the South Pole Wall.

Tools are provided at the Extragalactic Distance Database website that provide relationships between the distances and velocities of galaxies based on smoothed versions of the velocity fields derived by the Cosmicflows program.

Cosmic Dawn II (CoDa II) is a new, fully-coupled radiation-hydrodynamics simulation of cosmic reionization and galaxy formation and their mutual impact, to redshift $z < 6$. With $4096^3$ particles and cells in a 94 Mpc box, it is large enough to model global reionization and its feedback on galaxy formation while resolving all haloes above $10^8$ M$_{\odot}$. Using the same hybrid CPU-GPU code RAMSES-CUDATON as CoDa I in Ocvirk et al. (2016), CoDa II modified and re-calibrated the subgrid star-formation algorithm, making reionization end earlier, at $z \gtrsim 6$, thereby better matching the observations of intergalactic Lyman-alpha opacity from quasar spectra and electron-scattering optical depth from cosmic microwave background fluctuations. CoDa II predicts a UV continuum luminosity function in good agreement with observations of high-z galaxies, especially at $z = 6$. As in CoDa I, reionization feedback suppresses star formation in haloes below $\sim 2 \times 10^9$ M$_{\odot}$, though suppression here is less severe, a possible consequence of modifying the star-formation algorithm. Suppression is environment-dependent, occurring earlier (later) in overdense (underdense) regions, in response to their local reionization times. Using a constrained realization of $\Lambda$CDM constructed from galaxy survey data to reproduce the large-scale structure and major objects of the present-day Local Universe, CoDa II serves to model both global and local reionization. In CoDa II, the Milky Way and M31 appear as individual islands of reionization, i.e. they were not reionized by the progenitor of the Virgo cluster, nor by nearby groups, nor by each other.

In this work an ensemble of simulated Local Group analogues is used to constrain the properties of the mass assembly history of the Milky Way (MW) and Andromeda (M31) galaxies. These objects have been obtained using the constrained simulation technique, which ensures that simulated LGs live within a large scale environment akin to the observed one. Our results are compared against a standard $\Lambda$ Cold Dark Matter ($\Lambda$CDM) series of simulations which use the same cosmological parameters. This allows us to single out the effects of the constraints on the results. We find that (a) the median constrained merging histories for M31 and MW live above the standard ones at the 1-$\sigma$ level, (b) the median formation time takes place $\approx$ 0.5 Gyr earlier than unconstrained values, while the latest major merger happens on average 1.5 Gyr earlier and (c) the probability for both LG haloes to have experienced their last major merger in the first half of the history of the Universe is $\approx$ 50% higher for the constrained pairs. These results have been estimated to be significant at the 99% confidence level by means of a Kolmogorov-Simirnov test. These results suggest that the particular environment in which the Milky Way and Andromeda form plays a role in shaping their properties, and favours earlier formation and last major merger time values in agreement with other observational and theoretical considerations.

To date at least 10 highly flattened planes of dwarf galaxies are claimed to have been discovered in the Local Universe. The origin of these planes of galaxies remains unknown. One suggestion is that they are related to the large-scale structure of the cosmic web. A recent study found that the normal of a number of these dwarf galaxy planes are very closely aligned with the eigenvector of the shear tensor corresponding to the direction of greatest collapse obtained by reconstructing the full velocity field in the linear regime. Here we extend that work by both considering an additional 5 planes beyond the 5 examined previously and by examining the alignment with respect to the quasi-linear field, a more sophisticated reconstruction, which is a better approximation on smaller (quasi-linear) scales. Our analysis recovers the previous result while not finding a significantly tight alignment with the additional 5 planes. However, the additional 5 planes normals also do not appear to be randomly oriented. We conclude that this could either be due to the normals of the new planes being poorly defined and described; the quasi-linear field at those locations being poorly constrained; or different formation mechanisms for the orientation of planes of dwarf galaxies.

Cosmicflows-3 distances and inferred peculiar velocities of galaxies have permitted the reconstruction of the structure of over and under densities within the volume extending to 0.05c. This study focuses on the under dense regions, particularly the Local Void that lies largely in the zone of obscuration and consequently has received limited attention. Major over dense structures that bound the Local Void are the Perseus-Pisces and Norma-Pavo-Indus filaments separated by 8,500 km/s. The void network of the universe is interconnected and void passages are found from the Local Void to the adjacent very large Hercules and Sculptor voids. Minor filaments course through voids. A particularly interesting example connects the Virgo and Perseus clusters, with several substantial galaxies found along the chain in the depths of the Local Void. The Local Void has a substantial dynamical effect, causing a deviant motion of the Local Group of 200-250 km/s. The combined perturbations due to repulsion from the Local Void and attraction toward the Virgo Cluster account for about 50% of the motion of the Local Group in the rest frame given by the cosmic microwave background.

A hierarchical Bayesian model is applied to the Cosmicflows-3 catalog of galaxy distances in order to derive the peculiar velocity field and distribution of matter within $z \sim 0.054$. The model assumes the $\Lambda$CDM model within the linear regime and includes the fit of the galaxy distances together with the underlying density field. By forward modeling the data, the method is able to mitigate biases inherent to peculiar velocity analyses, such as the Homogeneous Malmquist bias or the log-normal distribution of peculiar velocities. The statistical uncertainty on the recovered velocity field is about 150 km/s depending on the location, and we study systematics coming from the selection function and calibration of distance indicators. The resulting velocity field and related density fields recover the cosmography of the Local Universe which is presented in an unprecedented volume of universe 10 times larger than previously reached. This methodology open the doors to reconstruction of initial conditions for larger and more accurate constrained cosmological simulations. This work is also preparatory to larger peculiar velocity datasets coming from Wallaby, TAIPAN or LSST.

Motivated by the ~7 % dipole anisotropy in the distribution of ultra-high energy cosmic-rays (UHECRs) above 8 EeV, we explore the anisotropy induced by the large scale structure, using constrained simulations of the local Universe and taking into account the effect of magnetic fields. The value of the intergalactic magnetic field (IGMF) is critical as it determines the UHECR cosmic horizon. We calculate the UHECR sky maps for different values of the IGMF variance and show the effect of the UHECR horizon on the observed anisotropy. The footprint of the local ($\lesssim350$ Mpc) Universe on the UHECR background, a small angular scale enhancement in the Northern Hemisphere, is seen. At 11.5 EeV (the median value of the energy bin at which the dipole has been reported), the LSS-induced dipole amplitude is $A_1\sim10\%$, for IGMF in the range [0.3-3] nG for protons, helium and nitrogen, compatible with the rms value derived from the cosmic power spectrum. However at these energies the UHECRs are also influenced by the Galactic Magnetic Field (GMF) and we discuss its effect on the LSS-induced anisotropy.

Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting different regions at different times. To establish the causal connection between reionization and suppression, we must take this local variation into account. We analyze the results of CoDa (`Cosmic Dawn') I, the first fully-coupled radiation-hydrodynamical simulation of reionization and galaxy formation in the Local Universe, in a volume large enough to model reionization globally but with enough resolving power to follow all atomic-cooling galactic halos in that volume. For every halo identified at a given time, we find the redshift at which the surrounding IGM reionized, along with its instantaneous star formation rate (`SFR') and baryonic gas-to-dark matter ratio ($M_\text{gas}/M_\text{DM}$). The average SFR per halo with $M < 10^9 \text{ M}_\odot$ was steady in regions not yet reionized, but declined sharply following local reionization. For $M > 10^{10} \text{ M}_\odot$, this SFR continued through local reionization, increasing with time, instead. For $10^9 < M < 10^{10} \text{ M}_\odot$, the SFR generally increased modestly through reionization, followed by a modest decline. In general, halo SFRs were higher for regions that reionized earlier. A similar pattern was found for $M_\text{gas}/M_\text{DM}$, which declined sharply following local reionization for $M < 10^9 \text{ M}_\odot$. Local reionization time correlates with local matter overdensity, which determines the local rates of structure formation and ionizing photon consumption. The earliest patches to develop structure and reionize ultimately produced more stars than they needed to finish and maintain their own reionization, exporting their `surplus' starlight to help reionize regions that developed structure later.

The local Universe provides a unique opportunity for testing cosmology and theories of structure formation. To facilitate this opportunity we present a new method for the reconstruction of the quasi-linear matter density and velocity fields from galaxy peculiar velocities and apply it to the Cosmicflows-2 data. The method consists of constructing an ensemble of cosmological simulations, constrained by the standard cosmological model and the observational data. The quasi-linear density field is the geometric mean and variance of the fully non-linear density fields of the simulations. The main nearby clusters (Virgo, Centaurus, Coma), superclusters (Shapley, Perseus-Pisces) and voids (Dipole Repeller) are robustly reconstructed. Galaxies are born biased with respect to the underlying dark matter distribution. Using our quasi-linear framework we demonstrate that the luminosity-weighted density field derived from the 2M++ redshift compilations is non-linearly biased with respect to the matter density field. The bias diminishes in the linear regime.

Today's galaxies experienced cosmic reionization at different times in different locations. For the first time, reionization ($50\%$ ionized) redshifts, $z_R$, at the location of their progenitors are derived from new, fully-coupled radiation-hydrodynamics simulation of galaxy formation and reionization at $z > 6$, matched to N-body simulation to z = 0. Constrained initial conditions were chosen to form the well-known structures of the local universe, including the Local Group and Virgo, in a (91 Mpc)$^3$ volume large enough to model both global and local reionization. Reionization simulation CoDa I-AMR, by CPU-GPU code EMMA, used (2048)$^3$ particles and (2048)$^3$ initial cells, adaptively-refined, while N-body simulation CoDa I-DM2048, by Gadget2, used (2048)$^3$ particles, to find reionization times for all galaxies at z = 0 with masses $M(z=0)\ge 10^8 M_\odot$. Galaxies with $M(z=0) \gtrsim 10^{11} M_\odot$ reionized earlier than the universe as a whole, by up to $\sim$ 500 Myrs, with significant scatter. For Milky-Way-like galaxies, $z_R$ ranged from 8 to 15. Galaxies with $M(z=0) \lesssim 10^{11} M_\odot$ typically reionized as late or later than globally-averaged $50\%$ reionization at $\langle z_R\rangle =7.8$, in neighborhoods where reionization was completed by external radiation. The spread of reionization times within galaxies was sometimes as large as the galaxy-to-galaxy scatter. The Milky Way and M31 reionized earlier than global reionization but later than typical for their mass, neither dominated by external radiation. Their most massive progenitors at $z>6$ had $z_R$ = 9.8 (MW) and 11 (M31), while their total masses had $z_R$ = 8.2 (both).

The fully nonlinear gravitationally induced trajectories of a nearly complete set of galaxies, groups, and clusters in the Local Supercluster are constructed in a Numerical Action Method (NAM) model constrained by data from the CosmicFlows survey and various distance indicators. We add the gravity field due to inhomogeneities external to the sample sphere by making use of larger scale peculiar flow measurements. Assignments of total masses were made to find the best overall set of mutual attractions, as determined by a goodness criterion based on present day radial velocities, individually for the Virgo Cluster, M31, and the Milky Way (MW), and via a mass-to-light ratio relationship for other masses. The low median chi-square found indicates the model fits the present day velocity flow well, but a slightly high mean chi-square may indicate that some masses underwent complex orbits. The best fit, when setting the value of $H_0$ to the CosmicFlows value of 75 km/s/Mpc and the WMAP value for $\Omega_m=0.244$ consistent with that $H_0$, occurs with the following parameters: $\Omega_{orphan} = 0.077\pm0.016$, $M/L_K $= $40 \pm 2$ $L_{10}^{0.15}$ $M_\odot/L_\odot$ ($L_{10}$ is K-band luminosity in units of $10^{10} L_\odot$), a Virgo mass of $6.3\pm0.8 \times 10^{14} M_\odot$ ($M/L_K = 113\pm15 M_\odot/L_\odot$), and a mass for the MW plus M31 of $5.15\pm0.35 \times 10^{12} M_\odot$. The best constant mass-to-light ratio is $M/L_K = 58\pm3 M_\odot/L_\odot$. The Virgocentric turnaround radius is $7.3\pm0.3$ Mpc. We explain several interesting trends in peculiar motions for various regions now that we can construct the 3D orbital histories.

We revisit the estimation of the combined mass of the Milky Way and Andromeda (M31), which dominate the mass of the Local Group. We make use of an ensemble of 30,190 halo pairs from the Small MultiDark simulation, assuming a $\Lambda$CDM (Cosmological Constant with Cold Dark Matter) cosmology, to investigate the relationship between the bound mass and parameters characterising the orbit of the binary and their local environment with the aid of machine learning methods (artificial neural networks, ANN). Results from the ANN are most successful when information about the velocity shear is provided, which demonstrates the flexibility of machine learning to model physical phenomena and readily incorporate new information as it becomes available. The resulting estimate for the Local Group mass, when shear information is included, is $4.9 \times 10^{12} M_\odot$, with an error of $\pm0.8 \times 10^{12} M_\odot$ from the 68% uncertainty in observables, and a 68% confidence interval of $^{+1.3}_{-1.4} \times 10^{12}M_\odot$ from the intrinsic scatter from the differences between the model and simulation masses. We also consider a recently reported large transverse velocity of M31 relative to the Milky Way, and produce an alternative mass estimate of $3.6\pm0.3\pm1.4 \times 10^{12}M_\odot$. Although different methods predict similar values for the most likely mass of the LG, application of ANN compared to the Timing Argument reduces the scatter in the log mass by over half when tested on samples from the simulation.

The influence of the large scale structure on host halos may be studied by examining the angular infall pattern of subhalos. In particular, since warm and cold dark matter cosmologies predict different abundances and internal properties for halos at the low mass end of the mass function, it is interesting to examine if there are differences in how these low mass halos are accreted. The accretion events are defined as the moment a halo becomes a substructure, namely when it crosses its host's virial radius. We quantify the cosmic web at each point by the shear tensor and examine where, with respect to its eigenvectors, such accretion events occur in cold ($\Lambda$CDM) and warm (1keV sterile neutrino WDM) dark matter cosmological models. We find that the CDM and WDM subhalos are preferentially accreted along the principal axis of the shear tensor corresponding to the direction of weakest collapse. The beaming strength is modulated by the host and subhalo masses and by the redshift at which the accretion event occurs. Although strongest for the most massive hosts and subhalos at high redshift, the preferential infall is found to be always aligned with the axis of weakest collapse, thus we say that it has universal nature. We compare the strength of beaming in the WDM cosmology with the one found in the $\Lambda$CDM scenario. While the main findings remain the same, the accretion in the WDM model for the most massive host halos appears more beamed than in $\Lambda$CDM cosmology across all the redshifts.

The three-dimensional gravitational velocity field within z~0.1 has been modeled with the Wiener filter methodology applied to the Cosmicflows-3 compilation of galaxy distances. The dominant features are a basin of attraction and two basins of repulsion. The major basin of attraction is an extension of the Shapley concentration of galaxies. One basin of repulsion, the Dipole Repeller, is located near the anti-apex of the cosmic microwave background dipole. The other basin of repulsion is in the proximate direction toward the 'Cold Spot' irregularity in the cosmic microwave background. It has been speculated that a vast void might contribute to the amplitude of the Cold Spot from the integrated Sachs-Wolfe effect.

The network of filaments with embedded clusters surrounding voids seen in maps derived from redshift surveys and reproduced in simulations has been referred to as the cosmic web. A complementary description is provided by considering the shear in the velocity field of galaxies. The eigenvalues of the shear provide information on whether a region is collapsing in three dimensions, the condition for a knot, expanding in three-dimensions, the condition for a void, or in the intermediate condition of a filament or sheet. The structures that are quantitatively defined by the eigenvalues can be approximated by iso-contours that provide a visual representation of the cosmic velocity (V) web. The current application is based on radial peculiar velocities from the Cosmicflows-2 collection of distances. The three-dimensional velocity field is constructed using the Wiener filter methodology in the linear approximation. Eigenvalues of the velocity shear are calculated at each point on a grid. Here, knots and filaments are visualized across a local domain of diameter ~0.1c.

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web -- depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper we bring twelve of these methods together and apply them to the same data set in order to understand how they compare. In general these cosmic web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore one would not {\it a priori} expect agreement between different techniques however, many of these methods do converge on the identification of specific features. In this paper we study the agreements and disparities of the different methods. For example, each method finds that knots inhabit higher density regions than filaments, etc. and that voids have the lowest densities. For a given web environment, we find substantial overlap in the density range assigned by each web classification scheme. We also compare classifications on a halo-by-halo basis; for example, we find that 9 of 12 methods classify around a third of group-mass haloes (i.e. $M_{\rm halo}\sim10^{13.5}h^{-1}M_{\odot}$) as being in filaments. Lastly, so that any future cosmic web classification scheme can be compared to the 12 methods used here, we have made all the data used in this paper public.

Observations of the Milky Way (MW), M31, and their vicinity, known as the Local Group (LG), can provide clues about the sources of reionization. We present a suite of radiative transfer simulations based on initial conditions provided by the Constrained Local UniversE Simulations (CLUES) project that are designed to recreate the Local Universe, including a realistic MW-M31 pair and a nearby Virgo. Our box size (91 Mpc) is large enough to incorporate the relevant sources of ionizing photons for the LG. We employ a range of source models, mimicking the potential effects of radiative feedback for dark matter haloes between $10^{8}-10^{9}$ M$_{\odot}$. Although the LG mostly reionizes in an inside-out fashion, the final 40 per cent of its ionization shows some outside influence. For the LG satellites, we find no evidence that their redshift of reionization is related to the present-day mass of the satellite or the distance from the central galaxy. We find that less than 20 per cent of present-day satellites for MW and M31 have undergone any star formation prior to the end of global reionization. Approximately five per cent of these satellites could be classified as fossils, meaning the majority of star formation occurred at these early times. The more massive satellites have more cumulative star formation prior to the end of global reionization, but the scatter is significant, especially at the low-mass end. Present-day mass and distance from the central galaxy are poor predictors for the presence of ancient stellar populations in satellite galaxies.

The Wiener Filter (WF) technique enables the reconstruction of density and velocity fields from observed radial peculiar velocities. This paper aims at identifying the optimal design of peculiar velocity surveys within the WF framework. The prime goal is to test the dependence of the quality of the reconstruction on the distribution and nature of data points. Mock datasets, extending to 250 Mpc/h, are drawn from a constrained simulation that mimics the local Universe to produce realistic mock catalogs. Reconstructed fields obtained with these mocks are compared to the reference simulation. Comparisons, including residual distributions, cell-to-cell and bulk velocities, imply that the presence of field data points is essential to properly measure the flows. The fields reconstructed from mocks that consist only of galaxy cluster data points exhibit poor quality bulk velocities. In addition, the quality of the reconstruction depends strongly on the grouping of individual data points into single points to suppress virial motions in high density regions. Conversely, the presence of a Zone of Avoidance hardly affects the reconstruction. For a given number of data points, a uniform sample does not score any better than a sample with decreasing number of data points with the distance. The best reconstructions are obtained with a grouped survey containing field galaxies: Assuming no error, they differ from the simulated field by less than 100 km/s up to the extreme edge of the catalogs or up to a distance of three times the mean distance of data points for non-uniform catalogs. The overall conclusions hold when errors are added.

In the standard (LCDM) model of cosmology the universe has emerged out of an early homogeneous and isotropic phase. Structure formation is associated with the growth of density irregularities and peculiar velocities. Our Local Group is moving with respect to the cosmic microwave background (CMB) with a velocity 631+/-20 km s-1 and participates in a bulk flow that extends out to distances of at least 20,000 km s-1. Since the discovery of the CMB dipole, the implicit assumption was that excesses in the abundance of galaxies induce the Local Group motion. Yet, underdense regions push as much as overdensities attract but they are deficient of light and consequently difficult to chart. It was suggested a decade ago that an underdensity in the northern hemisphere roughly 15,000 km s-1 away is a significant actor in the local flow. Here we report on kinematic evidence for such an underdensity. We map the large scale 3D velocity field using a Wiener filter reconstruction from the Cosmicflows-2 dataset of peculiar velocities, and identify the attractors and repellers that dominate the local dynamics. We show here that the local flow is dominated by a single attractor -associated with the Shapley Concentration- and a single previously unidentified repeller. Multipole expansion of the local flow provides further support for the existence and role played by the attractor and repeller. The bulk flow (i.e. dipole moment) is closely (anti)aligned with the repeller at a distance of 16,000+/-4,500 km s-1. The expansion eigenvector of the shear tensor (quadrupole moment) is closely aligned with the Shapley Attractor out to 7,000 km s-1. The close alignment of the local bulk flow with the repeller provides further support for its dominant role in shaping the local flow. This Dipole Repeller is predicted to be associated with a void in the distribution of galaxies.

Cosmography, the study and making of maps of the universe or cosmos, is a field where visual representation benefits from modern three-dimensional visualization techniques and media. At the extragalactic distance scales, visualization is contributing in understanding the complex structure of the local universe, in terms of spatial distribution and flows of galaxies and dark matter. In this paper, we report advances in the field of extragalactic cosmography obtained using the SDvision visualization software in the context of the Cosmicflows Project. Here, multiple visualization techniques are applied to a variety of data products: catalogs of galaxy positions and galaxy peculiar velocities, reconstructed velocity field, density field, gravitational potential field, velocity shear tensor viewed in terms of its eigenvalues and eigenvectors, envelope surfaces enclosing basins of attraction. These visualizations, implemented as high-resolution images, videos, and interactive viewers, have contributed to a number of studies: the cosmography of the local part of the universe, the nature of the Great Attractor, the discovery of the boundaries of our home supercluster of galaxies Laniakea, the mapping of the cosmic web, the study of attractors and repellers.

The mass of the Local Group (LG) is a crucial parameter for galaxy formation theories. However, its observational determination is challenging - its mass budget is dominated by dark matter which cannot be directly observed. To meet this end the posterior distributions of the LG and its massive constituents have been constructed by means of constrained and random cosmological simulations. Two priors are assumed - the LCDM model that is used to set up the simulations and a LG model,which encodes the observational knowledge of the LG and is used to select LG-like objects from the simulations. The constrained simulations are designed to reproduce the local cosmography as it is imprinted onto the Cosmicflows-2 database of velocities. Several prescriptions are used to define the LG model, focusing in particular on different recent estimates of the tangential velocity of M31. It is found that (a) different $v_{tan}$ choices affect the peak mass values up to a factor of 2, and change mass ratios of $M_{M31}$ to $M_{MW}$ by up to 20%; (b) constrained simulations yield more sharply peaked posterior distributions compared with the random ones; (c) LG mass estimates are found to be smaller than those found using the timing argument; (d) preferred MW masses lie in the range of $(0.6 - 0.8)\times10^{12} M_{\odot}$ whereas (e) $M_{M31}$ is found to vary between $(1.0 - 2.0)\times10^{12} M_{\odot}$, with a strong dependence on the $v_{tan}$ values used.

Cosmic reionization by starlight from early galaxies affected their evolution, thereby impacting reionization, itself. Star formation suppression, for example, may explain the observed underabundance of Local Group dwarfs relative to N-body predictions for Cold Dark Matter. Reionization modelling requires simulating volumes large enough [~(100Mpc)^3] to sample reionization "patchiness", while resolving millions of galaxy sources above ~10^8 Msun , combining gravitational and gas dynamics with radiative transfer. Modelling the Local Group requires initial cosmological density fluctuations pre-selected to form the well-known structures of the local universe today. Cosmic Dawn ("CoDa") is the first such fully-coupled, radiation-hydrodynamics simulation of reionization of the local universe. Our new hybrid CPU-GPU code, RAMSES-CUDATON, performs hundreds of radiative transfer and ionization rate-solver timesteps on the GPUs for each hydro-gravity timestep on the CPUs. CoDa simulated (91Mpc)^3 with 4096^3 particles and cells, to redshift 4.23, on ORNL supercomputer Titan, utilizing 8192 cores and 8192 GPUs. Global reionization ended slightly later than observed. However, a simple temporal rescaling which brings the evolution of ionized fraction into agreement with observations also reconciles ionizing flux density, cosmic star formation history, CMB electron scattering optical depth and galaxy UV luminosity function with their observed values. Photoionization heating suppressed the star formation of haloes below ~2 x 10^9 Msun , decreasing the abun- dance of faint galaxies around MAB_1600 = [-10,-12]. For most of reionization, star formation was dominated by haloes between 10^10 - 10^11 Msun , so low-mass halo suppression was not reflected by a distinct feature in the global star formation history. (Abridged)

An open question in cosmology and the theory of structure formation is to what extent does environment affect the properties of galaxies and haloes. The present paper aims at shedding light on this problem. The paper focuses on the analysis of a dark matter only simulation and it addresses the issue of how the environment affects the abundance of haloes, which are are assigned four attributes: their virial mass, an ambient density calculated with an aperture that scales with $R_{vir}$ ($\Delta_M$), a fixed-aperture ($\Delta_R$) ambient density, and a cosmic web classification (i.e. voids, sheets, filaments, and knots, as defined by the V--web algorithm). $\Delta_M$ is the mean density around a halo evaluated within a sphere of a radius of $5$\rvir, where \rvir\ is the virial radius. $\Delta_R$ is the density field Gaussian smoothed with $R=4h^{-1}$Mpc, evaluated at the center of the halo. The main result of the paper is that the difference between haloes in different web elements stems from the difference in their mass functions, and does not depend on their adaptive-aperture ambient density. A dependence on the fixed-aperture ambient density is induced by the cross correlation between the mass of a halo and its fixed-aperture ambient density.

While the Virgo cluster is the nearest galaxy cluster and therefore the best observed one, little is known about its formation history. In this paper, a set of cosmological simulations that resemble the Local Universe is used to shed the first light on this mystery. The initial conditions for these simulations are constrained with galaxy peculiar velocities of the second catalog of the Cosmicflows project using algorithms developed within the Constrained Local UniversE Simulation project. Boxes of 500 Mpc/h on a side are set to run a series of dark matter only constrained simulations. In each simulation, a unique dark matter halo can be reliably identified as Virgo's counterpart. The properties of these Virgo halos are in agreement at a 10-20% level with the global properties of the observed Virgo cluster. Their zero-velocity masses agree at one-sigma with the observational mass estimate. In all the simulations, the matter falls onto the Virgo objects along a preferential direction that corresponds to the observational filament and the slowest direction of collapse. A study of the mass accretion history of the Virgo candidates reveals the most likely formation history of the Virgo cluster, namely a quiet accretion over the last 7 Gigayears.

The goodness-of-fit (GoF) of the Cosmicflows-2 (CF2) database of peculiar velocities with the LCDM standard model of cosmology is presented. Standard application of the Chi^2 statistics of the full database, of its 4,838 data points, is hampered by the small scale nonlinear dynamics which is not accounted for by the (linear regime) velocity power spectrum. The bulk velocity constitutes a highly compressed representation of the data which filters out the small scales non-linear modes. Hence the statistics of the bulk flow provides an efficient tool for assessing the GoF of the data given a model. The particular approach introduced here is to use the (spherical top-hat window) bulk velocity extracted from the Wiener filter reconstruction of the 3D velocity field as a linear low pass filtered highly compressed representation of the CF2 data. An ensemble 2250 random linear realizations of the WMAP/LCDM model has been used to calculate the bulk velocity auto-covariance matrix. We find that the CF2 data is consistent with the WMAP/LCDM model to better than the 2 sigma confidence limits. This provides a further validation that the CF2 database is consistent with the standard model of cosmology.

Determining the precise value of the tangential component of the velocity of M31 is a non trivial astrophysical issue, that relies on complicated modeling. This has recently lead to con- flicting estimates, obtained by several groups that used different methodologies and assump- tions. This letter addresses the issue by computing a Bayesian posterior distribution function of this quantity, in order to measure the compatibility of those estimates with LambdaCDM. This is achieved using an ensemble of local group (LG) look-alikes collected from a set of Con- strained Simulations (CSs) of the local Universe, and a standard unconstrained LambdaCDM. The latter allows us to build a control sample of LG-like pairs and to single out the influence of the environment in our results. We find that neither estimate is at odds with LambdaCDM; how- ever, whereas CSs favour higher values of vtan , the reverse is true for estimates based on LG samples gathered from unconstrained simulations, overlooking the environmental element

Near field cosmology is practiced by studying the Local Group (LG) and its neighbourhood. The present paper describes a framework for simulating the near field on the computer. Assuming the LCDM model as a prior and applying the Bayesian tools of the Wiener filter (WF) and constrained realizations of Gaussian fields to the Cosmicflows-2 (CF2) survey of peculiar velocities, constrained simulations of our cosmic environment are performed. The aim of these simulations is to reproduce the LG and its local environment. Our main result is that the LG is likely a robust outcome of the LCDM scenario when subjected to the constraint derived from CF2 data, emerging in an environment akin to the observed one. Three levels of criteria are used to define the simulated LGs. At the base level, pairs of halos must obey specific isolation, mass and separation criteria. At the second level the orbital angular momentum and energy are constrained and on the third one the phase of the orbit is constrained. Out of the 300 constrained simulations 146 LGs obey the first set of criteria, 51 the second and 6 the third. The robustness of our LG factory enables the construction of a large ensemble of simulated LGs. Suitable candidates for high resolution hydrodynamical simulations of the LG can be drawn from this ensemble, which can be used to perform comprehensive studies of the formation of the LG

Tidal torque theory suggests that galaxies gain angular momentum in the linear stage of structure formation. Such a theory predicts alignments between the spin of haloes and tidal shear field. However, non-linear evolution and angular momentum acquisition may alter this prediction significantly. In this paper, we use a reconstruction of the cosmic shear field from observed peculiar velocities combined with spin axes extracted from galaxies within $115\, \mathrm{Mpc} $ ($\sim8000 \, {\mathrm {km}}{\mathrm s}^{-1}$) from 2MRS catalog, to test whether or not galaxies appear aligned with principal axes of shear field. Although linear reconstructions of the tidal field have looked at similar issues, this is the first such study to examine galaxy alignments with velocity-shear field. Ellipticals in the 2MRS sample, show a statistically significant alignment with two of the principal axes of the shear field. In general, elliptical galaxies have their short axis aligned with the axis of greatest compression and perpendicular to the axis of slowest compression. Spiral galaxies show no signal. Such an alignment is significantly strengthened when considering only those galaxies that are used in velocity field reconstruction. When examining such a subsample, a weak alignment with the axis of greatest compression emerges for spiral galaxies as well. This result indicates that although velocity field reconstructions still rely on fairly noisy and sparse data, the underlying alignment with shear field is strong enough to be visible even when small numbers of galaxies are considered - especially if those galaxies are used as constraints in the reconstruction.

We use a cosmological simulation of the formation of the Local Group to explore the origin of age and metallicity gradients in dwarf spheroidal galaxies. We find that a number of simulated dwarfs form "outside-in", with an old, metal-poor population that surrounds a younger, more concentrated metal-rich component, reminiscent of dwarf spheroidals like Sculptor or Sextans. We focus on a few examples where stars form in two populations distinct in age in order to elucidate the origin of these gradients. The spatial distributions of the two components reflect their diverse origin; the old stellar component is assembled through mergers, but the young population forms largely in situ. The older component results from a first episode of star formation that begins early but is quickly shut off by the combined effects of stellar feedback and reionization. The younger component forms when a late accretion event adds gas and reignites star formation. The effect of mergers is to disperse the old stellar population, increasing their radius and decreasing their central density relative to the young population. We argue that dwarf-dwarf mergers offer a plausible scenario for the formation of systems with multiple distinct populations and, more generally, for the origin of age and metallicity gradients in dwarf spheroidals.

This paper combines observational datasets and cosmological simulations to generate realistic numerical replicas of the nearby Universe. These latter are excellent laboratories for studies of the non-linear process of structure formation in our neighborhood. With measurements of radial peculiar velocities in the Local Universe (cosmicflows-2) and a newly developed technique, we produce Constrained Local UniversE Simulations (CLUES). To assess the quality of these constrained simulations, we compare them with random simulations as well as with local observations. The cosmic variance, defined as the mean one-sigma scatter of cell-to-cell comparison between two fields, is significantly smaller for the constrained simulations than for the random simulations. Within the inner part of the box where most of the constraints are, the scatter is smaller by a factor 2 to 3 on a 5 Mpc/h scale with respect to that found for random simulations. This one-sigma scatter obtained when comparing the simulated and the observation-reconstructed velocity fields is only 104 +/- 4 km/s i.e. the linear theory threshold. These two results demonstrate that these simulations are in agreement with each other and with the observations of our neighborhood. For the first time, simulations constrained with observational radial peculiar velocities resemble the Local Universe up to a distance of 150 Mpc/h on a scale of a few tens of megaparsecs. When focusing on the inner part of the box, the resemblance with our cosmic neighborhood extends to a few megaparsecs (< 5 Mpc/h). The simulations provide a proper Large Scale environment for studies of the formation of nearby objects.

Superclusters of galaxies can be defined kinematically from local evaluations of the velocity shear tensor. The location where the smallest eigenvalue of the shear is positive and maximal defines the center of a basin of attraction. Velocity and density fields are reconstructed with Wiener Filter techniques. Local velocities due to the density field in a restricted region can be separated from external tidal flows, permitting the identification of boundaries separating inward flows toward a basin of attraction and outward flows. This methodology was used to define the Laniakea Supercluster that includes the Milky Way. Large adjacent structures include Perseus-Pisces, Coma, Hercules, and Shapley but current kinematic data are insufficient to capture their full domains. However there is a small region trapped between Laniakea, Perseus-Pisces, and Coma that is close enough to be reliably characterized and that satisfies the kinematic definition of a supercluster. Because of its shape, it is given the name the Arrowhead Supercluster. This entity does not contain any major clusters. A characteristic dimension is ~25 Mpc and the contained mass is only ~10^15 Msun.

Recent observational studies have demonstrated that the majority of satellite galaxies tend to orbit their hosts on highly flattened, vast, possibly co-rotating planes. Two nearly parallel planes of satellites have been confirmed around the M31 galaxy and around the Centaurus A galaxy, while the Milky Way also sports a plane of satellites. It has been argued that such an alignment of satellites on vast planes is unexpected in the standard ({\Lambda}CDM) model of cosmology if not even in contradiction to its generic predictions. Guided by {\Lambda}CDM numerical simulations, which suggest that satellites are channeled towards hosts along the axis of the slowest collapse as dictated by the ambient velocity shear tensor, we re-examine the planes of local satellites systems within the framework of the local shear tensor derived from the Cosmicflows-2 dataset. The analysis reveals that the Local Group and Centaurus A reside in a filament stretched by the Virgo cluster and compressed by the expansion of the Local Void. Four out of five thin planes of satellite galaxies are indeed closely aligned with the axis of compression induced by the Local Void. Being the less massive system, the moderate misalignment of the Milky Way's satellite plane can likely be ascribed to its greater susceptibility to tidal torques, as suggested by numerical simulations. The alignment of satellite systems in the local universe with the ambient shear field is thus in general agreement with predictions of the {\Lambda}CDM model.

The cosmic web that characterizes the large-scale structure of the Universe can be quantified by a variety of methods. For example, large redshift surveys can be used in combination with point process algorithms to extract long curvilinear filaments in the galaxy distribution. Alternatively, given a full 3D reconstruction of the velocity field, kinematic techniques can be used to decompose the web into voids, sheets, filaments and knots. In this paper we look at how two such algorithms - the Bisous model and the velocity shear web - compare with each other in the local Universe (within 100 Mpc), finding good agreement. This is both remarkable and comforting, given that the two methods are radically different in ideology and applied to completely independent and different data sets. Unsurprisingly, the methods are in better agreement when applied to unbiased and complete data sets, like cosmological simulations, than when applied to observational samples. We conclude that more observational data is needed to improve on these methods, but that both methods are most likely properly tracing the underlying distribution of matter in the Universe.

Full sky surveys of peculiar velocity are arguably the best way to map the large scale structure out to distances of a few times 100 Mpc/h. Using the largest and most accurate ever catalog of galaxy peculiar velocities "Cosmicflows-2", the large scale structure has been reconstructed by means of the Wiener filter and constrained realizations assuming as a Bayesian prior model the LCDM model with the WMAP inferred cosmological parameters. The present paper focuses on studying the bulk flow of the local flow field, defined as the mean velocity of top-hat spheres with radii ranging out to R=500 Mpc/h. The estimated large scale structures, in general, and the bulk flow, in particular, are determined by the tension between the observational data and the assumed prior model. A prerequisite for such an analysis is the requirement that the estimated bulk flow is consistent with the prior model. Such a consistency is found here. At R=50(150) Mpc/h the estimated bulk velocity is 250+/-21 (239+/-38) km/s. The corresponding cosmic variance at these radii is 126(60)km/s, which implies that these estimated bulk flows are dominated by the data and not by the assumed prior model. The estimated bulk velocity is dominated by the data out to R~200 Mpc/h, where the cosmic variance on the individual Supergalactic Cartesian components (of the r.m.s. values) exceeds the variance of the Constrained Realizations by at least a factor of 2. The supergalactic SGX and SGY components of the CMB dipole velocity are recovered by the Wiener filter velocity field down to a very few km/s. The SGZ component of the estimated velocity, the one that is most affected by the Zone of Avoidance, is off by 126 km/s (an almost 2 sigma discrepancy).

We search for vast planes of satellites (VPoS) in a high resolution simulation of the Local Group performed by the CLUES project, which improves significantly the resolution of former similar studies. We use a simple method for detecting planar configurations of satellites, and validate it on the known plane of M31. We implement a range of prescriptions for modelling the satellite populations, roughly reproducing the variety of recipes used in the literature, and investigate the occurence and properties of planar structures in these populations. The structure of the simulated satellite systems is strongly non-random and contains planes of satellites, predominantly co-rotating, with, in some cases, sizes comparable to the plane observed in M31 by Ibata et al.. However the latter is slightly richer in satellites, slightly thinner and has stronger co-rotation, which makes it stand out as overall more exceptional than the simulated planes, when compared to a random population. Although the simulated planes we find are generally dominated by one real structure, forming its backbone, they are also partly fortuitous and are thus not kinematically coherent structures as a whole. Provided that the simulated and observed planes of satellites are indeed of the same nature, our results suggest that the VPoS of M31 is not a coherent disc and that one third to one half of its satellites must have large proper motions perpendicular to the plane.

We seek to understand the relationship between galaxy properties and their local environment, which calls for a proper formulation of the notion of environment. We analyse the GIMIC suite of cosmological hydrodynamical simulations within the framework of the cosmic web as formulated by Hoffman et al., focusing on properties of simulated DM haloes and luminous galaxies with respect to voids, sheets, filaments and knots - the four elements of the cosmic web. We find that the mass functions of haloes depend on environment, which drives other environmental dependence of galaxy formation. The web shapes the halo mass function, and through the strong dependence of the galaxy properties on the mass of their host haloes, it also shapes the galaxy-(web) environment dependence.

We use high resolution simulations of the formation of the local group post-processed by a radiative transfer code for UV photons, to investigate the reionization of the satellite populations of an isolated Milky Way-M31 galaxy pair in a variety of scenarios. We use an improved version of ATON which includes a simple recipe for radiative feedback. In our baseline models, reionization is initiated by low mass, radiatively regulated haloes at high redshift, until more massive haloes appear, which then dominate and complete the reionization process. We investigate the relation between reionization history and present-day positions of the satellite population. We find that the average reionization redshift (zr) of satellites is higher near galaxy centers (MW and M31). This is due to the inside-out reionization patterns imprinted by massive haloes within the progenitor during the EoR, which end up forming the center of the galaxy. Thanks to incomplete dynamical mixing during galaxy assembly, these early patterns survive down to present day, resulting in a a clear radial gradient in the average satellites reionization redshift, up to the virial radius of MW and M31 and beyond. In the lowest emissivity scenario, the outer satellites are reionized about 180 Myr later than the inner satellites. This delay decreases with increasing source model emissivity, or in the case of external reionization by Virgo or M31, because reionization happens faster overall, and becomes spatially quasi-uniform at the highest emissivity.

Galaxies congregate in clusters and along filaments, and are missing from large regions referred to as voids. These structures are seen in maps derived from spectroscopic surveys that reveal networks of structure that are interconnected with no clear boundaries. Extended regions with a high concentration of galaxies are called 'superclusters', although this term is not precise. There is, however, another way to analyse the structure. If the distance to each galaxy from Earth is directly measured, then the peculiar velocity can be derived from the subtraction of the mean cosmic expansion, the product of distance times the Hubble constant, from observed velocity. The peculiar velocity is the line-of-sight departure from the cosmic expansion and arises from gravitational perturbations; a map of peculiar velocities can be translated into a map of the distribution of matter. Here we report a map of structure made using a catalogue of peculiar velocities. We find locations where peculiar velocity flows diverge, as water does at watershed divides, and we trace the surface of divergent points that surrounds us. Within the volume enclosed by this surface, the motions of galaxies are inward after removal of the mean cosmic expansion and long range flows. We define a supercluster to be the volume within such a surface, and so we are defining the extent of our home supercluster, which we call Laniakea. A key component of this paper is an accompanying movie that can be viewed (also in 3D) and downloaded at http://irfu.cea.fr/laniakea or http://vimeo.com/pomarede/laniakea

We examine the angular infall pattern of subhaloes onto host haloes in the context of the large-scale structure. We find that this infall pattern is essentially driven by the shear tensor of the ambient velocity field. Dark matter subhaloes are found to be preferentially accreted along the principal axis of the shear tensor which corresponds to the direction of weakest collapse. We examine the dependence of this preferential infall on subhalo mass, host halo mass and redshift. Although strongest for the most massive hosts and the most massive subhaloes at high redshift, the preferential infall of subhaloes is effectively universal in the sense that its always aligned with the axis of weakest collapse of the velocity shear tensor. It is the same shear tensor that dictates the structure of the cosmic web and hence the shear field emerges as the key factor that governs the local anisotropic pattern of structure formation. Since the small (sub-Mpc) scale is strongly correlated with the mid-range (~10 Mpc) scale - a scale accessible by current surveys of peculiar velocities - it follows that findings presented here open a new window into the relation between the observed large scale structure unveiled by current surveys of peculiar velocities and the preferential infall direction of the Local Group. This may shed light on the unexpected alignments of dwarf galaxies seen in the Local Group

The evolution of the large scale distribution of matter in the universe is often characterized by the density field. Here we take a complimentary approach and characterize it using the cosmic velocity field, specifically the deformation of the velocity field. The deformation tensor is decomposed into its symmetric component (known as the "shear tensor") and its anti-symmetric part (the "vorticity"). Using a high resolution cosmological simulation we examine the relative orientations of the shear and the vorticity as a function of spatial scale and redshift. The shear is found to be remarkable stable to the choice of scale, while the vorticity is found to quickly decay with increasing spatial scale or redshift. The vorticity emerges out of the linear regime randomly oriented with respect to the shear eigenvectors. Non-linear evolution drives the vorticity to lie within the plane defined by the eigenvector of the fastest collapse. Within that plane the vorticity first gets aligned with the middle eigenvector and then it moves to be preferentially aligned with the third eigenvector, of slowest collapse. Finally a scale of "non-linearity" to be used when calculating properties of the non-linear deformation tensor at different redshifts is suggested.

We explore the impact of cosmic reionization on nearby isolated dwarf galaxies using a compilation of SFHs estimated from deep HST data and a cosmological hydrodynamical simulation of the Local Group. The nearby dwarfs show a wide diversity of star formation histories; from ancient systems that have largely completed their star formation $\sim 10$ Gyr ago to young dwarfs that have formed the majority of their stars in the past $\sim 5$ Gyr to two-component systems characterized by the overlap of comparable numbers of old and young stars. Taken as an ensemble, star formation in nearby dwarfs dips to lower-than-average rates at intermediate times ($4<t$/Gyr $<8$), a feature that we trace in the simulation to the effects of cosmic reionization. Reionization heats the gas and drives it out of the shallow potential wells of low mass halos, affecting especially those below a sharp mass threshold that corresponds to a virial temperature of $\sim 2 \times 10^4 $ $\mathrm{K}$ at $z_{\rm reion}$. The loss of baryons leads to a sharp decline in the star forming activity of early-collapsing systems, which, compounded by feedback from early star formation, empties halos of gas and leaves behind systems where a single old stellar component prevails. In halos below the threshold at $z_{\rm reion}$, reionization heating leads to a delay in the onset of star formation that lasts until the halo grows massive enough to allow some of the remaining gas to cool and form stars. Young stellar components therefore dominate in dwarfs whose halos assemble late and thus form few stars before reionization. Two-component systems may be traced to late mergers of individual examples of the two aforementioned cases. The relative dearth of intermediate-age stars in nearby dwarfs might thus be the clearest signature yet identified of the imprint of cosmic reionization on the star formation history of dwarf galaxies.

We contend that a single power law halo mass distribution is appropriate for direct matching to the stellar masses of observed Local Group dwarf galaxies, allowing the determination of the slope of the stellar mass-halo mass relation for low mass galaxies. Errors in halo masses are well defined as the Poisson noise of simulated local group realisations, which we determine using constrained local universe simulations (CLUES). For the stellar mass range 10$^7$<M*<10$^8$M$_\odot$, for which we likely have a complete census of observed galaxies, we find that the stellar mass-halo mass relation follows a power law with slope of 3.1, significantly steeper than most values in the literature. The steep relation between stellar and halo masses indicates that Local Group dwarf galaxies are hosted by dark matter halos with a small range of mass. Our methodology is robust down to the stellar mass to which the census of observed Local Group galaxies is complete, but the significant uncertainty in the currently measured slope of the stellar-to halo mass relation will decrease dramatically if the Local Group completeness limit was $10^{6.5}$M$\odot$ or below, highlighting the importance of pushing such limit to lower masses and larger volumes.

We review how dark matter is distributed in our local neighbourhood from an observational and theoretical perspective. We will start by describing first the dark matter halo of our own galaxy and in the Local Group. Then we proceed to describe the dark matter distribution in the more extended area known as the Local Universe. Depending on the nature of dark matter, numerical simulations predict different abundances of substructures in Local Group galaxies, in the number of void regions and in the abundance of low rotational velocity galaxies in the Local Universe. By comparing these predictions with the most recent observations, strong constrains on the physical properties of the dark matter particles can be derived. We devote particular attention to the results from the Constrained Local UniversE Simulations (CLUES) project, a special set of simulations whose initial conditions are constrained by observational data from the Local Universe. The resulting simulations are designed to reproduce the observed structures in the nearby universe. The CLUES provides a numerical laboratory for simulating the Local Group of galaxies and exploring the physics of galaxy formation in an environment designed to follow the observed Local Universe. It has come of age as the numerical analogue of Near-Field Cosmology.

Peculiar velocities, obtained from direct distance measurements, are data of choice to achieve constrained simulations of the Local Universe reliable down to a scale of a few Megaparsecs. Unlike redshift surveys, peculiar velocities are direct tracers of the underlying gravitational field as they trace both baryonic and dark matter. This paper presents the first attempt to use solely observational peculiar velocities to constrain cosmological simulations of the nearby universe. In order to set up Initial Conditions, a Reverse Zel'dovich Approximation (RZA) is used to displace constraints from their positions at z=0 to their precursors' locations at higher redshifts. An additional new feature replaces original observed radial peculiar velocity vectors by their full 3D reconstructions provided by the Wiener-Filter (WF) estimator. Subsequently, the Constrained Realization of Gaussian fields technique (CR) is applied to build various realizations of the Initial Conditions. The WF/RZA/CR method is first tested on realistic mock catalogs built from a reference simulation similar to the Local Universe. These mocks include errors on peculiar velocities, on data-point positions and a large continuous zone devoid of data in order to mimic galactic extinction. Large scale structures are recovered with a typical accuracy of 5 Megaparsecs/h in position, the best realizations reaching a 2-3 Mpc/h precision, the limit imposed by the RZA linear theory. Then, the method is applied to the first observational radial peculiar velocity catalog of the project Cosmicflows. This paper is a proof of concept that the WF/RZA/CR method can be applied to observational peculiar velocities to successfully build constrained Initial Conditions.

We present the results of a set of numerical simulations aimed at studying reionization at galactic scale. We use a high resolution simulation of the formation of the Milky Way-M31 system to simulate the reionization of the local group. The reionization calculation was performed with the post-processing radiative transfer code ATON and the underlying cosmological simulation was performed as part of the CLUES project. We vary the source models to bracket the range of source properties used in the literature. We investigate the structure and propagation of the galatic ionization fronts by a visual examination of our reionization maps. Within the progenitors we find that reionization is patchy, and proceeds locally inside out. The process becomes patchier with decreasing source photon output. It is generally dominated by one major HII region and 1-4 additional isolated smaller bubbles, which eventually overlap. Higher emissivity results in faster and earlier local reionization. In all models, the reionization of the Milky Way and M31 are similar in duration, i.e. between 203 Myr and 22 Myr depending on the source model, placing their zreion between 8.4 and 13.7. In all models except the most extreme, the MW and M31 progenitors reionize internally, ignoring each other, despite being relatively close to each other even during the epoch of reionization. Only in the case of strong supernova feedback suppressing star formation in haloes less massive than 10^9 M_sun, and using our highest emissivity, we find that the MW is reionized by M31.

The large-scale structure of the Universe is characterised by a web-like structure made of voids, sheets, filaments, and knots. The structure of this so-called cosmic web is dictated by the local velocity shear tensor. In particular, the local direction of a filament should be strongly aligned with e3, the eigenvector associated with the smallest eigenvalue of the tensor. That conjecture is tested here on the basis of a cosmological simulation. The cosmic web delineated by the halo distribution is probed by a marked point process with interactions (the Bisous model), detecting filaments directly from the halo distribution (P-web). The detected P-web filaments are found to be strongly aligned with the local e3: the alignment is within 30 degree for 80% of the elements. This indicates that large-scale filaments defined purely from the distribution of haloes carry more than just morphological information, although the Bisous model does not make any prior assumption on the underlying shear tensor. The P-web filaments are also compared to the structure revealed from the velocity shear tensor itself (V-web). In the densest regions, the P- and V-web filaments overlap well (90%), whereas in lower density regions, the P-web filaments preferentially mark sheets in the V-web.

We revise the mass estimate of the Local Group (LG) when Dark Energy (in the form of the Cosmological Constant) is incorporated into the Timing Argument (TA) mass estimator for the Local Group (LG). Assuming the age of the Universe and the Cosmological Constant according to the recent values from the Planck CMB experiment, we find the mass of the LG to be M_TAL = (4.73 +- 1.03) x 10^{12} M_sun, which is 13% higher than the classical TA mass estimate. This partly explains the discrepancy between earlier results from LCDM simulations and the classical TA. When a similar analysis is performed on 16 LG-like galaxy pairs from the CLUES simulations, we find that the scatter in the ratio of the virial to the TA estimated mass is given by M_vir/M_TAL = 1.04 +-0.16. Applying it to the LG mass estimation we find a calibrated M_vir = (4.92 +- 1.08 (obs) +- 0.79 (sys)) x 10^{12} M_sun.

The large scale structure of the universe is a complex web of clusters, filaments, and voids. Its properties are informed by galaxy redshift surveys and measurements of peculiar velocities. Wiener Filter reconstructions recover three-dimensional velocity and total density fields. The richness of the elements of our neighborhood are revealed with sophisticated visualization tools. A key component of this paper is an accompanying movie that can be viewed and downloaded at : http://irfu.cea.fr/cosmography or http://vimeo.com/pomarede/cosmography .The ability to translate and zoom helps the viewer follow structures in three dimensions and grasp the relationships between features on different scales while retaining a sense of orientation. The ability to dissolve between scenes provides a technique for comparing different information, for example, the observed distribution of galaxies, smoothed representations of the distribution accounting for selection effects, observed peculiar velocities, smoothed and modeled representations of those velocities, and inferred underlying density fields. The agreement between the large scale structure seen in redshift surveys and that inferred from reconstructions based on the radial peculiar velocities of galaxies strongly supports the standard model of cosmology where structure forms from gravitational instabilities and galaxies form at the bottom of potential wells.

The differences between cold (CDM) and warm (WDM) dark matter in the formation of a group of galaxies is examined by running two identical simulations where in the WDM case the initial power spectrum has been altered to mimic a 1keV dark matter particle. The CDM initial conditions were constrained to reproduce at z = 0 the correct local environment within which a "Local Group" (LG) of galaxies may form. Two significant differences between the two simulations are found. While in the CDM case a group of galaxies that resembles the real LG forms, the WDM run fails to reproduce a viable LG, instead forming a diffuse group which is still expanding at z = 0. This is surprising since, due to the suppression of small scale power in its power spectrum, WDM is naively expected to only affect the collapse of small haloes and not necessarily the dynamics on a scale of a group of galaxies. Furthermore the concentration of baryons in halo center's is greater in CDM than in WDM and the properties of the disks differ.

Recent observations constrained the tangential velocity of M31 with respect to the Milky Way (MW) to be v_tan<34.4 km/s and the radial velocity to be in the range v_rad=-109+/- 4.4 km/s (van der Marel et al. 2012). In this study we use a large volume high resolution N-body cosmological simulation (Bolshoi) together with three constrained simulations to statistically study this kinematics in the context of the LCDM. The comparison of the ensembles of simulated pairs with the observed LG at the 1-sigma level in the uncertainties has been done with respect to the radial and tangential velocities, the reduced orbital energy (e_tot), angular momentum (l_orb) and the dimensionless spin parameter, lambda. Our main results are: (i) the preferred radial and tangential velocities for pairs in LCDM are v_rad=-80+/-20 km/s, v_tan=50+/-10 km/s, (ii) pairs around that region are 3 to 13 times more common than pairs within the observational values, (iii) 15%to 24% of LG-like pairs in LCDM have energy and angular momentum consistent with observations while (iv) 9% to 13% of pairs in the same sample show similar values in the inferred dimensionless spin parameter. It follows that within current observational uncertainties the quasi-conserved quantities that characterize the orbit of the LG, i.e. e_tot, r_orb and lambda, do not challenge the standard LCDM model, but the model is in tension with regard to the actual values of the radial and tangential velocities. This might hint to a problem of the LCDM model to reproduce the observed LG.

We use dark matter only and full hydrodynamical Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group to study the density profile of subhaloes of the simulated Milky Way and Andromeda galaxies. We show that the Einasto model provides the best description of the subhaloes' density profile, as opposed to the more commonly used NFW profile or any generalisation of it. We further find that the Einasto shape parameter \nEin\ is strongly correlated with the total subhalo mass, pointing towards the notion of a non-universality of the subhaloes' density profile. We observe that the effect of mass loss due to tidal stripping, in both the dark matter only and the hydrodynamical run, is the reduction of the shape parameter \nEin\ between the infall and the present time. Assuming now that the dSphs of our Galaxy follow the Einasto profile and using the maximum and minimum values of \nEin\ from our hydrodynamical simulation as a gauge, we can improve the observational constraints on the \Rmax-\Vmax\ pairs obtained for the brightest satellite galaxies of the Milky Way. When considering only the subhaloes with $-13.2\lesssim M_V\lesssim-8.8$, i.e. the range of luminosity of the classical dwarfs, we find that all our simulated objects are consistent with the observed dSphs if their haloes follow the Einasto model with $1.6\lesssim n_{\rm E} \lesssim5.3$. The numerically motivated Einasto profile for the observed dSphs will alleviate the recently presented "massive failures" problem.

The Reverse Zeldovich Approximation (RZA) is a reconstruction method which allows to estimate the cosmic displacement field from galaxy peculiar velocity data and to constrain initial conditions for cosmological simulations of the Local Universe. In this paper, we investigate the effect of different observational errors on the reconstruction quality of this method. For this, we build a set of mock catalogues from a cosmological simulation, varying different error sources like the galaxy distance measurement error (0 - 20%), the sparseness of the data points, and the maximum catalogue radius (3000 - 6000 km/s). We perform the RZA reconstruction of the initial conditions on these mock catalogues and compare with the actual initial conditions of the simulation. We also investigate the impact of the fact that only the radial part of the peculiar velocity is observationally accessible. We find that the sparseness of a dataset has the highest detrimental effect on RZA reconstruction quality. Observational distance errors also have a significant influence, but it is possible to compensate this relatively well with Wiener Filter reconstruction. We also investigate the effect of different object selection criteria and find that distance catalogues distributed randomly and homogeneously across the sky (such as spiral galaxies selected for the Tully-Fisher method) allow for a higher reconstruction quality than if when data is preferentially drawn from massive objects or dense environments (such as elliptical galaxies). We find that the error of estimating the initial conditions with RZA is always dominated by the inherent non-linearity of data observed at z=0 rather than by the combined effect of the observational errors. Even an extremely sparse dataset with high observational errors still leads to a good reconstruction of the initial conditions on a scale of about 5 Mpc/h.

We propose a new method to recover the cosmological initial conditions of the presently observed galaxy distribution, which can serve to run constrained simulations of the Local Universe. Our method, the Reverse Zeldovich Approximation (RZA), can be applied to radial galaxy peculiar velocity data and extends the previously used Constrained Realizations (CR) method by adding a Lagrangian reconstruction step. The RZA method consists of applying the Zeldovich approximation in reverse to galaxy peculiar velocities to estimate the cosmic displacement field and the initial linear matter distribution from which the present-day Local Universe evolved.We test our method with a mock survey taken from a cosmological simulation. We show that the halo peculiar velocities at z = 0 are close to the linear prediction of the Zeldovich approximation, if a grouping is applied to the data to remove virial motions. We find that the addition of RZA to the CR method significantly improves the reconstruction of the initial conditions. The RZA is able to recover the correct initial positions of the velocity tracers with a median error of only 1.36 Mpc/h in our test simulation. For realistic sparse and noisy data, this median increases to 5 Mpc/h. This is a significant improvement over the previous approach of neglecting the displacement field, which introduces errors on a scale of 10 Mpc/h or even higher. Applying the RZA method to the upcoming high-quality observational peculiar velocity catalogues will generate much more precise constrained simulations of the Local Universe.

In previous works we proposed the Reverse Zeldovich Approximation (RZA) method, which can be used to estimate the cosmological initial conditions underlying the galaxy distribution in the Local Universe using peculiar velocity data. In this paper, we apply the technique to run constrained cosmological simulations from the RZA-reconstructed initial conditions, designed to reproduce the large-scale structure of the Local Universe. We test the method with mock peculiar velocity catalogues extracted from a reference simulation. We first reconstruct the initial conditions of this reference simulation using the mock data, and then run the reconstructed initial conditions forward in time until z=0. We compare the resulting constrained simulations with the original simulation at z=0 to test the accuracy of this method. We also compare them with constrained simulations run from the mock data without the addition of RZA, i.e. using only the currently established constrained realizations (CR) method. Our re-simulations are able to correctly recover the evolution of the large-scale structure underlying the data. The results show that the addition of RZA to the CR method significantly improves both the reconstruction of the initial conditions and the accuracy of the obtained constrained resimulations. Haloes from the original simulation are recovered in the re-simulations with an average accuracy of about 2 Mpc/h on their position and a factor of 2 in mass, down to haloes with a mass of approx 10^14 M_odot/h. In comparison, without RZA the re-simulations recover only the most massive haloes with masses of about 5x10^14 M_odot/h and higher, and with a systematic shift on their position of about approx 10 Mpc/h due to the cosmic displacement field. We show that with the additional Lagrangian reconstruction step introduced by the RZA, this shift can be removed.

In the standard model of cosmology, structure emerges out of non-rotational flow and the angular momentum of collapsing halos is induced by tidal torques. The growth of halo angular momentum in the linear and quasi-linear phases is associated with a shear, curl-free, flow and it is well described within the linear framework of tidal torque theory (TTT). However, TTT is rendered irrelevant as haloes approach turn around and virialization. At that stage the flow field around halos has non-zero vorticity. Using a cosmological simulation, we have examined the importance of the curl of the velocity field (vorticity) in determining halo spin, finding a strong alignment between the two. We have also examined the alignment of vorticity with the principle axes of the shear tensor, finding that it tends to be perpendicular to the axis along which material is collapsing fastest (e1). This behavior is independent of halo masses and cosmic web environment. Our results agree with previous findings on the tendency of halo spin to be perpendicular to e1, and of the spin of (simulated) halos and (observed) galaxies to be aligned with the large-scale structure. Our results imply that angular momentum growth proceeds in two distinct phases. In the first phase angular momentum emerges out of a shear, curl-free, potential flow, as described by TTT. In the second phase, in which haloes approach virialization, the angular momentum emerges out of a vortical flow and halo spin becomes strongly aligned with the vorticity of the ambient flow field.

We present a cosmography analysis of the Local Universe based on the recently released Two-Micron All-Sky Redshift Survey (2MRS). Our method is based on a Bayesian Networks Machine Learning algorithm (the Kigen-code) which self-consistently samples the initial density fluctuations compatible with the observed galaxy distribution and a structure formation model given by second order Lagrangian perturbation theory (2LPT). From the initial conditions we obtain an ensemble of reconstructed density and peculiar velocity fields which characterize the local cosmic structure with high accuracy unveiling nonlinear structures like filaments and voids in detail. Coherent redshift space distortions are consistently corrected within 2LPT. From the ensemble of cross-correlations between the reconstructions and the galaxy field and the variance of the recovered density fields we find that our method is extremely accurate up to k ~ 1 h Mpc^-1 and still yields reliable results down to scales of about 3-4 h^-1 Mpc. The motion of the local group we obtain within ~ 80 h^-1 Mpc (v_LG=522+-86 km s^-1, l_LG=291^o +- 16^o, b_LG=34^o+-8^o) is in good agreement with measurements derived from the CMB and from direct observations of peculiar motions and is consistent with the predictions of LambdaCDM.

We use a cosmological simulation of the formation of the Local Group of Galaxies to identify a mechanism that enables the removal of baryons from low-mass halos without appealing to feedback or reionization. As the Local Group forms, matter bound to it develops a network of filaments and pancakes. This moving web of gas and dark matter drifts and sweeps a large volume, overtaking many halos in the process. The dark matter content of these halos is unaffected but their gas can be efficiently removed by ram-pressure. The loss of gas is especially pronounced in low-mass halos due to their lower binding energy and has a dramatic effect on the star formation history of affected systems. This "cosmic web stripping" may help to explain the scarcity of dwarf galaxies compared with the numerous low-mass halos expected in \Lambda CDM and the large diversity of star formation histories and morphologies characteristic of faint galaxies. Although our results are based on a single high-resolution simulation, it is likely that the hydrodynamical interaction of dwarf galaxies with the cosmic web is a crucial ingredient so far missing from galaxy formation models.

The alignment of DM halos and the surrounding large scale structure (LSS) is examined in the context of the cosmic web. Halo spin, shape and the orbital angular momentum of subhaloes is investigated relative to the LSS using the eigenvectors of the velocity shear tensor evaluated on a grid with a scale of 1 Mpc/h, deep within the non-linear regime. Knots, filaments, sheets and voids are associated with regions that are collapsing along 3, 2, 1 or 0 principal directions simultaneously. Each halo is tagged with a web classification (i.e. knot halo, filament halo, etc) according to the nature of the collapse at the halo's position. The full distribution of shear eigenvalues is found to be substantially different from that tagged to haloes, indicating that the observed velocity shear is significantly biased. We find that larger mass haloes live in regions where the shear is more isotropic, namely the expansion or collapse is more spherical. A correlation is found between the halo's shape and the eigenvectors of the shear tensor, with the longest (shortest) axis of the halo's shape being aligned with the slowest (fastest) collapsing eigenvector. This correlation is web independent, suggesting that the velocity shear is a fundamental tracer of the halo alignment. A similar result is found for the alignment of halo spin with the cosmic web. It has been shown that high mass haloes exhibit a spin flip with respect to the LSS: we find the mass at which this spin flip occurs is web dependent and not universal as suggested previously. Although weaker than haloes, subhalo orbits too exhibit an alignment with the LSS, providing a possible insight into the highly correlated co-rotation of the Milky Way's satellite system. The present study suggests that the velocity shear tensor constitutes the natural framework for studying the directional properties of the non-linear LSS and of halos and galaxies.

A new approach for the classification of the cosmic web is presented. In extension of the previous work of Hahn et al. (2007) and Forero-Romero et al. (2009) the new algorithm is based on the analysis of the velocity shear tensor rather than the gravitational tidal tensor. The procedure consists of the construction of the the shear tensor at each (grid) point in space and the evaluation of its three eigenvectors. A given point is classified to be either a void, sheet, filament or a knot according to the number of eigenvalues above a certain threshold, 0, 1, 2, or 3 respectively. The threshold is treated as a free parameter that defines the web. The algorithm has been applied to a dark matter only, high resolution simulation of a box of side-length 64$h^{-1}$Mpc and N = $1024^3$ particles with the framework of the WMAP5/LCDM model. The resulting velocity based cosmic web resolves structures down to <0.1$h^{-1}$Mpc scales, as opposed to the ~1$h^{-1}$Mpc scale of the tidal based web. The under-dense regions are made of extended voids bisected by planar sheets, whose density is also below the mean. The over-dense regions are vastly dominated by the linear filaments and knots. The resolution achieved by the velocity based cosmic web provides a platform for studying the formation of halos and galaxies within the framework of the cosmic web.

In this work we investigate the nonlinear and nonlocal relation between cosmological density and peculiar velocity fields. Our goal is to provide an algorithm for the recon- struction of the nonlinear velocity field from the fully nonlinear density. We find that including the gravitational tidal field tensor using second order Lagrangian perturba- tion theory (2LPT) based upon an estimate of the linear component of the nonlinear density field significantly improves the estimate of the cosmic flow in comparison to linear theory not only in the low density, but also and more dramatically in the high density regions. In particular we test two estimates of the linear component: the log- normal model and the iterative Lagrangian linearisation. The present approach relies on a rigorous higher order Lagrangian perturbation theory analysis which incorpo- rates a nonlocal relation. It does not require additional fitting from simulations being in this sense parameter free, it is independent of statistical-geometrical optimisation and it is straightforward and efficient to compute. The method is demonstrated to yield an unbiased estimator of the velocity field on scales ~> 5 Mpc/h with closely Gaussian distributed errors. Moreover, the statistics of the divergence of the peculiar velocity field is extremely well recovered showing a good agreement with the true one from N-body simulations. The typical errors of about 10 km/s (1 sigma confidence intervals) are reduced by more than 80% with respect to linear theory in the scale range between 5 and 10 Mpc/h in high density regions ({\delta} > 2). We also find that iterative Lagrangian linearisation is significantly superior in the low density regime with respect to the lognormal model.

We use the recently proposed scale-free mass estimators to determine the masses of the Milky Way (MW) and Andromeda (M31) galaxy in a dark matter only Constrained Local UniversE Simulation (CLUES). While these mass estimators work rather well for isolated spherical host systems, we examine here their applicability to a simulated binary system with a unique satellite population similar to the observed satellites of MW and M31. We confirm that the scale-free estimators work also very well in our simulated Local Group galaxies with the right number of satellites which follow the observed radial distribution. In the isotropic case and under the assumption that the satellites are tracking the total gravitating mass, the power-law index of the radial satellite distribution $N(<r)\propto r^{3-\gamma}$ is directly related to the host's mass profile $M(<r)\propto r^{1-\alpha}$ as $\alpha=\gamma-2$. The use of this relation for any given $\gamma$ leads to highly accurate mass estimations which is a crucial point for observer, since they do not know a priori the mass profile of the MW and M31 haloes. We discuss possible bias in the mass estimators and conclude that the scale-free mass estimators can be satisfactorily applied to the real MW and M31 system.

We present all-sky simulated Fermi maps of gamma-rays from dark matter decay and annihilation in the Local Universe. The dark matter distribution is obtained from a constrained cosmological simulation of the neighboring large-scale structure provided by the CLUES project. The dark matter fields of density and density squared are then taken as an input for the Fermi observation simulation tool to predict the gamma-ray photon counts that Fermi would detect in 5 years of all-sky survey for given dark matter models. Signal-to-noise sky maps have also been obtained by adopting the current Galactic and isotropic diffuse background models released by the Fermi collaboration. We point out the possibility for Fermi to detect a dark matter gamma-ray signal in local extragalactic structures. In particular, we conclude here that Fermi observations of nearby clusters (e.g. Virgo and Coma) and filaments are expected to give stronger constraints on decaying dark matter compared to previous studies. As an example, we find a significant signal-to-noise ratio in dark matter models with a decay rate fitting the positron excess as measured by PAMELA. This is the first time that dark matter filaments are shown to be promising targets for indirect detection of dark matter. On the other hand, the prospects for detectability of annihilating dark matter in local extragalactic structures are less optimistic even with extreme cross-sections. We make the dark matter density and density squared maps available online at http://www.clues-project.org/articles/darkmattermaps.html

We use a 64$h^{-1}$Mpc dark matter (DM) only cosmological simulation to examine the large scale orientation of haloes and substructures with respect the cosmic web. A web classification scheme based on the velocity shear tensor is used to assign to each halo in the simulation a web type: knot, filament, sheet or void. Using $\sim10^6$ haloes that span ~3 orders of magnitude in mass the orientation of the halo's spin and the orbital angular momentum of subhaloes with respect to the eigenvectors of the shear tensor is examined. We find that the orbital angular momentum of subhaloes tends to align with the intermediate eigenvector of the velocity shear tensor for all haloes in knots, filaments and sheets. This result indicates that the kinematics of substructures located deep within the virialized regions of a halo is determined by its infall which in turn is determined by the large scale velocity shear, a surprising result given the virilaized nature of haloes. The non-random nature of subhalo accretion is thus imprinted on the angular momentum measured at z = 0. We also find that haloes' spin axis is aligned with the third eigenvector of the velocity shear tensor in filaments and sheets: the halo spin axis points along filaments and lies in the plane of cosmic sheets.

This paper presents an analysis of the local peculiar velocity field based on the Wiener Filter reconstruction method. We used our currently available catalog of distance measurements containing 1,797 galaxies within 3000 km/s: Cosmicflows-1. The Wiener Filter method is used to recover the full 3D peculiar velocity field from the observed map of radial velocities and to recover the underlying linear density field. The velocity field within a data zone of 3000 km/s is decomposed into a local component that is generated within the data zone and a tidal one that is generated by the mass distribution outside that zone. The tidal component is characterized by a coherent flow toward the Norma-Hydra-Centaurus (Great Attractor) region while the local component is dominated by a flow toward the Virgo Cluster and away from the Local Void. A detailed analysis shows that the local flow is predominantly governed by the Local Void and the Virgo Cluster plays a lesser role. The analysis procedure was tested against a mock catalog. It is demonstrated that the Wiener Filter accurately recovers the input velocity field of the mock catalog on the scale of the extraction of distances and reasonably recovers the velocity field on significantly larger scales. The Bayesian Wiener Filter reconstruction is carried out within the ?CDM WMAP5 framework. The Wiener Filter reconstruction draws particular attention to the importance of voids in proximity to our neighborhood. The prominent structure of the Local Supercluster is wrapped in a horseshoe collar of under density with the Local Void as a major component.

Using Constrained Local UniversE Simulations (CLUES) of the formation of the Local Group in a cosmological context we investigate the recently highlighted problem that the majority of the most massive dark subhaloes of the Milky Way are too dense to host any of its bright satellites. In particular, we examine the influence of baryonic processes and find that they leave a twofold effect on the relation between the peak of the rotation curve and its position (Vmax and Rmax). Satellites with a large baryon fraction experience adiabatic contraction thus decreasing Rmax while leaving Vmax more or less unchanged. Subhaloes with smaller baryon fractions undergo a decrease in Vmax possibly due to outflows of material. Furthermore, the situation of finding subhaloes in simulations that lie outside the confidence interval for possible hosts of the bright MW dwarf spheroidals, appears to be far more prominent in cosmologies with a high sigma8 normalisation and depends on the mass of the host. We conclude that the problem cannot be simply solved by including baryonic processes and hence demands further investigations.

The outer haloes of the Milky Way and Andromeda galaxies contain as much important information on their assembly and formation history as the properties of the discs resident in their centres. In this paper we have used the Constrained Local UniversE Simulation project to disentangle the stellar and DM component of three galaxies that resemble the MW, M31 and M33 using both DM only and DM + gas-dynamical simulations. Stars that are accreted in substructures and then stripped follow a completely different radial distribution than the stripped DM: the stellar halo is much more centrally concentrated than DM. In order to understand how the same physical process can lead to different z=0 radial profiles, we examined the potential at accretion of each stripped particle. We found that star particles sit at systematically higher potentials than DM, making them harder to strip. We then searched for a threshold in the potential of accreted particles \phi_[th], above which DM particles behave as star particles. We found such a threshold at >16 \phi_{subhalo}, where \phi_{subhalo} is the potential at a subhaloes edge at the time of accretion. Thus a rule as simple as selecting particles according to their potential at accretion is able to reproduce the effect that the complicated physics of star formation has on the stellar distribution. This niversal result reproduces the stellar halo to an accuracy of within ~2%. Studies which make use of DM particles as a proxy for stars will undoubtedly miscalculate their proper radial distribution and structure unless particles are selected according to their potential at accretion. Furthermore, we have examined the time it takes to strip a given star or DM particle after accretion. We find that, owing to their higher binding energies, stars take longer to be stripped than DM. The stripped DM halo is thus considerably older than the stripped stellar halo.

Using a dark matter only Constrained Local UniversE Simulation (CLUES) we examine the existence of subhaloes that change their affiliation from one of the two prominent hosts in the Local Group (i.e. the Milky Way and the Andromeda galaxy) to the other, and call these objects "renegade subhaloes". In light of recent claims that the two Magellanic Clouds (MCs) may have originated from another region (or even the outskirts) of the Local Group or that they have been spawned by a major merger in the past of the Andromeda galaxy, we investigate the nature of such events. However, we cannot confirm that renegade subhaloes enter as deep into the potential well of their present host nor that they share the most simplest properties with the MCs, namely mass and relative velocity. Our simulation rather suggests that these renegade subhaloes appear to be flying past one host before being pulled into the other. A merger is not required to trigger such an event, it is rather the distinct environment of our simulated Local Group facilitating such behavior. Since just a small fraction of the full z=0 subhalo population are renegades, our study indicates that it will be intrinsically difficult to distinguish them despite clear differences in their velocity, radial distribution, shape and spin parameter distributions.

We make detailed theoretical predictions for the assembly properties of the Local Group (LG) in the standard LambdaCDM cosmological model. We use three cosmological N-body dark matter simulations from the CLUES project, which are designed to reproduce the main dynamical features of the matter distribution down to the scale of a few Mpc around the LG. Additionally, we use the results of an unconstrained simulation with a sixty times larger volume to calibrate the influence of cosmic variance. We characterize the Mass Aggregation History (MAH) for each halo by three characteristic times, the formation, assembly and last major merger times. A major merger is defined by a minimal mass ratio of 10:1. We find that the three LGs share a similar MAH with formation and last major merger epochs placed on average \approx 10 - 12 Gyr ago. Between 12% and 17% of the halos in the mass range 5 x 10^11 Msol/h < M_h < 5 x 10^12 Msol/h have a similar MAH. In a set of pairs of halos within the same mass range, a fraction of 1% to 3% share similar formation properties as both halos in the simulated LG. An unsolved question posed by our results is the dynamical origin of the MAH of the LGs. The isolation criteria commonly used to define LG-like halos in unconstrained simulations do not narrow down the halo population into a set with quiet MAHs, nor does a further constraint to reside in a low density environment. The quiet MAH of the LGs provides a favorable environment for the formation of disk galaxies like the Milky Way and M31. The timing for the beginning of the last major merger in the Milky Way dark matter halo matches with the gas rich merger origin for the thick component in the galactic disk. Our results support the view that the specific large and mid scale environment around the Local Group play a critical role in shaping its MAH and hence its baryonic structure at present.

Observed high-redshift QSOs, at z~6, may reside in massive dark matter (DM) halos of more than 10^{12} Msun and are thus expected to be surrounded by overdense regions. In a series of 10 constrained simulations, we have tested the environment of such QSOs. Comparing the computed overdensities with respect to the unconstrained simulations of regions empty of QSOs, assuming there is no bias between the DM and baryon distributions, and invoking an observationally-constrained duty-cycle for Lyman Break Galaxies, we have obtained the galaxy count number for the QSO environment. We find that a clear discrepancy exists between the computed and observed galaxy counts in the Kim et al. (2009) samples. Our simulations predict that on average eight z~6 galaxies per QSO field should have been observed, while Kim et al. detect on average four galaxies per QSO field compared to an average of three galaxies in a control sample (GOODS fields). While we cannot rule out a small number statistics for the observed fields to high confidence, the discrepancy suggests that galaxy formation in the QSO neighborhood proceeds differently than in the field. We also find that QSO halos are the most massive of the simulated volume at z~6 but this is no longer true at z~3. This implies that QSO halos, even in the case they are the most massive ones at high redshifts, do not evolve into most massive galaxy clusters at z=0.

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

We study the differences and similarities in the luminosities of bound, infalling and the so-called backsplash (Gill et al. 2005) galaxies of the Milky Way and M31 using a hydrodynamical simulation performed within the Constrained Local UniversE Simulation (CLUES) project. The simulation models the formation of the Local Group within a self-consistent cosmological framework. We find that even though backsplash galaxies passed through the virial radius of their host halo and hence may have lost a (significant) fraction of their mass, their stellar populations are hardly affected. This leaves us with comparable luminosity functions for infalling and backsplash galaxies and hence little hope to decipher their past (and different) formation and evolutionary histories by luminosity measurements alone. Nevertheless, due to the tidal stripping of dark matter we find that the mass-to-light ratios have changed when comparing the various populations against each other: they are highest for the infalling galaxies and lowest for the bound satellites with the backsplash galaxies in-between.

Using a high resolution DM simulation of the Local Group (LG), conducted within the framework of the Constrained Local UniversE Simulation (CLUES) project, we investigate the nature of how satellites of the MW and M31 are accreted. Satellites of these galaxies are accreted anisotropically, entering the virial radius of their hosts from specific "spots" with respect to the large scale structure. Furthermore, material which is tidally stripped from these satellites is also, at z=0, distributed anisotropically and is characterized by an ellipsoidal sub-volume embedded in the halo. The angular pattern created by the locus of satellite infall points and the projected stripped DM is investigated within a coordinate system determined by the location of the LG companion and the simulated Virgo cluster across concentric shells ranging from 0.1 to 5 r_vir. Remarkably, the principal axis of the ellipsoidal sub-volume shows a coherent alignment extending from well within the halo to a few r_vir. A spherical harmonics transform applied to the angular distributions confirms the visual impression: namely, the angular distributions of both the satellites entry points and stripped DM is dominated by the l=2 quadrupole term, whose major principal axis is aligned across the shells considered. It follows that the structure of the outer (r>0.5 r_vir) regions of the main halos is closely related to the cosmic web within which it is embedded. Given the hypothesis that a large fraction of the stellar halo of the Milky Way has been accreted from satellites, the present results can be applied to the stellar halo of the MW and M31. We predict that the remnants of tidally stripped satellites should be embedded in streams of material composed of dark matter and stars. The present results can therefore shed light on the existence of satellites embedded within larger streams of matter, such as the Segue 2 satellite.

We compare the substructure evolution in pure dark matter (DM) halos with those in the presence of baryons (PDM and BDM). The prime halos have been analyzed by Romano-Diaz et al (2009). Models have been evolved from identical initial conditions using Constrained Realizations, including star formation and feedback. A comprehensive catalog of subhalos has been compiled and properties of subhalos analyzed in the mass range of 10^8 Mo - 10^11 Mo. We find that subhalo mass functions are consistent with a single power law, M_sbh^{alpha}, but detect a nonnegligible shift between these functions, alpha -0.86 for the PDM, and -0.98 for the BDM. Overall, alpha const. in time with variations of +-15%. Second, we find that the radial mass distribution of subhalos can be approximated by a power law, R^{gamma} with a steepening around the radius of a maximal circular velocity, Rvmax, in the prime halos. Gamma ~-1.5 for the PDM and -1 for the BDM, inside Rvmax, and is steeper outside. We detect little spatial bias between the subhalo populations and the DM of the main halos. The subhalo population exhibits much less triaxiality with baryons, in tandem with the prime halo. Finally, we find that, counter-intuitively, the BDM population is depleted at a faster rate than the PDM one within the central 30kpc of the prime. Although the baryons provide a substantial glue to the subhalos, the main halos exhibit the same trend. This assures a more efficient tidal disruption of the BDM subhalos. This effect can be reversed for a more efficient feedback from stellar evolution and supermassive black holes, which will expel baryons from the center and decrease the concentration of the prime halo. We compare our results with via Lactea and Aquarius simulations and other published results.

We present the first detailed structure formation and radiative transfer simulations of the reionization history of our cosmic neighbourhood. To this end, we follow the formation of the Local Group of galaxies and nearby clusters by means of constrained simulations, which use the available observational constraints to construct a representation of those structures which reproduces their actual positions and properties at the present time. We find that the reionization history of the Local Group is strongly dependent on the assumed photon production efficiencies of the ionizing sources, which are still poorly constrained. If sources are relatively efficient, i.e. the process is 'photon-rich', the Local Group is primarily ionized externally by the nearby clusters. Alternatively, if the sources are inefficient, i.e. reionization is 'photon-poor' the Local Group evolves largely isolated and reionizes itself. The mode of reionization, external vs. internal, has important implications for the evolution of our neighbourhood, in terms of e.g. its satellite galaxy populations and primordial stellar populations. This therefore provides an important avenue for understanding the young universe by detailed studies of our nearby structures.

The local universe is the best known part of our universe. Within the CLUES project (http://clues-project.org - Constrained Local UniversE Simulations) we perform numerical simulations of the evolution of the local universe. For these simulations we construct initial conditions based on observational data of the galaxy distribution in the local universe. Here we review the technique of these constrained simulations. In the second part we summarize our predictions of a possible Warm Dark Matter cosmology for the observed local distribution of galaxies and the local spectrum of mini-voids as well as a study of the satellite dynamics in a simulated Local Group.

We use two simulations performed within the Constrained Local UniversE Simulation (CLUES) project to study both the shape and radial alignment of (the dark matter component of) subhaloes; one of the simulations is a dark matter only model while the other run includes all the relevant gas physics and star formation recipes. We find that the involvement of gas physics does not have a statistically significant effect on either property -- at least not for the most massive subhaloes considered in this study. However, we observe in both simulations including and excluding gasdynamics a (pronounced) evolution of the dark matter shapes of subhaloes as well as of the radial alignment signal since infall time. Further, this evolution is different when positioned in the central and outer regions of the host halo today; while subhaloes tend to become more aspherical in the central 50% of their host's virial radius, the radial alignment weakens in the central regime while strengthening in the outer parts. We confirm that this is due to tidal torquing and the fact that subhaloes at pericentre move too fast for the alignment signal to respond.

We use a simulation performed within the Constrained Local UniversE Simulation (CLUES) project to study a realistic Local Group-like object. We employ this group as a numerical laboratory for studying the evolution of the population of its subhaloes from the point of view of the effects it may have on the origin of different types of dwarf galaxies. We focus on the processes of tidal stripping of the satellites, their interaction, merging and grouping before infall. The tidal stripping manifests itself in the transition between the phase of mass accretion and mass loss seen in most subhaloes, which occurs at the moment of infall on to the host halo, and the change of the shape of their mass function with redshift. Although the satellites often form groups, they are loosely bound within them and do not interact with each other. The infall of a large group could however explain the observed peculiar distribution of the Local Group satellites, but only if it occurred recently. Mergers between prospective subhaloes are significant only during an early stage of evolution, i.e. more than 7 Gyr ago, when they are still outside the host haloes. Such events could thus contribute to the formation of more distant early type Milky Way companions. Once the subhaloes enter the host halo the mergers become very rare.

We examine the properties of satellites found in high resolution simulations of the local group. We use constrained simulations designed to reproduce the main dynamical features that characterize the local neighborhood, i.e. within tens of Mpc around the Local Group (LG). Specifically, a LG-like object is found located within the 'correct' dynamical environment and consisting of three main objects which are associated with the Milky Way, M31 and M33. By running two simulations of this LG from identical initial conditions - one with and one without baryons modeled hydrodynamically - we can quantify the effect of gas physics on the $z=0$ population of subhaloes in an environment similar to our own. We find that above a certain mass cut, $M_{\rm sub} > 2\times10^{8}h^{-1} M_{\odot}$ subhaloes in hydrodynamic simulations are more radially concentrated than those in simulations with out gas. This is caused by the collapse of baryons into stars that typically sit in the central regions of subhaloes, making them denser. The increased central density of such a subhalo, results in less mass loss due to tidal stripping than the same subhalo simulated with only dark matter. The increased mass in hydrodynamic subhaloes with respect to dark matter ones, causes dynamical friction to be more effective, dragging the subhalo towards the centre of the host. This results in these subhaloes being effectively more radially concentrated then their dark matter counterparts.

We compare assembly of DM halos with and without baryons, within the context of cosmological evolution in the LCDM WMAP3 Universe (baryons+DM, BDM model, and pure DM, PDM model). In representative PDM and BDM models we find that baryons contribute decisively to the evolution of the central region, leading to an isothermal DM cusp, and to a flat DM density core -- the result of heating by dynamical friction of the substructure during a quiescent evolution epoch. This process ablates the cold gas from an embedded disk, cutting the star formation rate by ~10, and heats up the spheroidal gas and stellar components, triggering their expansion. The substructure is more resilient in the presence of baryons. The disk which formed from inside-out as gas dominated, is transformed into an intermediate Hubble type by z ~ 2 and to an early type by z ~ 0.5, based on its gas contents and spheroidal-to-disk stellar mass ratio. Only a relatively small ~20% fraction of DM particles in PDM and BDM models are bound within the radius of maximal circular velocity in the halo -- most of the DM particles perform larger radial excursions. We also find that the fraction of baryons within the halo virial radius somewhat increases during the major mergers and decreases during the minor mergers. The net effect appears to be negligible. While the substructure is being tidally-disrupted, mixing of its debris in the halo is not efficient and becomes even less so with z. The streamers formed after z ~ 1 survive largely to the present time -- an important implication for embedded disk evolution.

Using high-resolution simulations within the Cold and Warm Dark Matter models we study the evolution of small scale structure in the Local Volume, a sphere of 8 Mpc radius around the Local Group. We compare the observed spectrum of mini-voids in the Local Volume with the spectrum of mini-voids determined from the simulations. We show that the \LWDM model can easily explain both the observed spectrum of mini-voids and the presence of low-mass galaxies observed in the Local Volume, provided that all haloes with circular velocities greater than 20 km/s host galaxies. On the contrary within the LCDM model the distribution of the simulated mini-voids reflects the observed one if haloes with maximal circular velocities larger than 35 km/s host galaxies. This assumption is in contradiction with observations of galaxies with circular velocities as low as 20 km/s in our Local Universe. A potential problem of the LWDM model could be the late formation of the haloes in which the gas can be efficiently photo-evaporated. Thus star formation is suppressed and low-mass haloes might not host any galaxy at all.

Using constrained simulations of the local Universe for generic cold dark matter and for 1keV warm dark matter, we investigate the difference in the abundance of dark matter halos in the local environment. We find that the mass function within 20 Mpc/h of the Local Group is ~2 times larger than the universal mass function in the 10^9-10^13 M_odot/h mass range. Imposing the field of view of the on-going HI blind survey ALFALFA in our simulations, we predict that the velocity function in the Virgo-direction region exceeds the universal velocity function by a factor of 3. Furthermore, employing a scheme to translate the halo velocity function into a galaxy velocity function, we compare the simulation results with a sample of galaxies from the early catalog release of ALFALFA. We find that our simulations are able to reproduce the velocity function in the 80-300 km/s velocity range, having a value ~10 times larger than the universal velocity function in the Virgo-direction region. In the low velocity regime, 35-80 km/s, the warm dark matter simulation reproduces the observed flattening of the velocity function. On the contrary, the simulation with cold dark matter predicts a steep rise in the velocity function towards lower velocities; for V_max=35 km/s, it forecasts ~10 times more sources than the ones observed. If confirmed by the complete ALFALFA survey, our results indicate a potential problem for the cold dark matter paradigm or for the conventional assumptions about energetic feedback in dwarf galaxies.

Using a suite of N-body simulations in different Cold Dark Matter (CDM) scenarios, with cosmological constant (\LCDM) and without (OCDM, SCDM), we study the Hubble flow (\sigh) in Local Volumes (LV) around Local Group (LG) like objects found in these simulations, and compare the numerical results with the most recent observations. We show that \LCDM and OCDM models exhibit the same behavior of \sigh. Hence, we demonstrate that the observed coldness of the Hubble flow is not likely to be a manifestation of the dark energy, contrary to previous claims. The coldness does not constitute a problem by itself but it poses a problem to the standard \LCDM model only if the mean density within the Local Volume is greater than twice the mean matter cosmic density. The lack of blueshifted galaxies in the LV, outside of the LG can be considered as another manifestation of the coldness of the flow. Finally, we show that the main dynamical parameter that affects the coldness of the flow is the relative isolation of the LG, and the absence of nearby Milky Way like objects within a distance of about $3\mpc$.

A dynamical classification of the cosmic web is proposed. The large scale environment is classified into four web types: voids, sheets, filaments and knots. The classification is based on the evaluation of the deformation tensor, i.e. the Hessian of the gravitational potential, on a grid. The classification is based on counting the number of eigenvalues above a certain threshold, lambda_th at each grid point, where the case of zero, one, two or three such eigenvalues corresponds to void, sheet, filament or a knot grid point. The collection of neighboring grid points, friends-of-friends, of the same web attribute constitutes voids, sheets, filaments and knots as web objects. A simple dynamical consideration suggests that lambda_th should be approximately unity, upon an appropriate scaling of the deformation tensor. The algorithm has been applied and tested against a suite of (dark matter only) cosmological N-body simulations. In particular, the dependence of the volume and mass filling fractions on lambda_th and on the resolution has been calculated for the four web types. Also, the percolation properties of voids and filaments have been studied. Our main findings are: (a) Already at lambda_th = 0.1 the resulting web classification reproduces the visual impression of the cosmic web. (b) Between 0.2 < lambda_th < 0.4, a system of percolated voids coexists with a net of interconected filaments. This suggests a reasonable choice for lambda_th as the parameter that defines the cosmic web. (c) The dynamical nature of the suggested classification provides a robust framework for incorporating environmental information into galaxy formation models, and in particular the semi-analytical ones.

We study formation and evolution of bar-disk systems in fully self-consistent cosmological simulations of galaxy formation in the LCDM WMAP3 Universe. In a representative model we find that the first generation of bars form in response to the asymmetric dark matter (DM) distribution (i.e., DM filament) and quickly decay. Subsequent bar generations form and are destroyed during the major merger epoch permeated by interactions with a DM substructure (subhalos). A long-lived bar is triggered by a tide from a subhalo and survives for ~10 Gyr. The evolution of this bar is followed during the subsequent numerous minor mergers and interactions with the substructure. Together with intrinsic factors, these interactions largely determine the stellar bar evolution. The bar strength and its pattern speed anticorrelate, except during interactions and when the secondary (nuclear) bar is present. For about 5 Gyr bar pattern speed increases substantially despite the loss of angular momentum to stars and cuspy DM halo. We analyze the evolution of stellar populations in the bar-disk and relate them to the underlying dynamics. While the bar is made mainly of an intermediate age, ~5-6 Gyr, disk stars at z=0, a secondary nuclear bar which surfaces at z~0.1 is made of younger, ~1-3 Gyr stars.

We use a set of high-resolution simulations of scale-free Einstein-de Sitter cosmologies to investigate the logarithmic slope of the phase-space density profile $Q(r) = \rho(r)/\sigma^3(r)$ of dark matter (DM) haloes. The initial conditions for the simulations are determined by a power law power spectrum of the form $P(k) \propto k^n$. We compute the Q(r) profiles using the radial, tangential and full velocity dispersion, and the velocity anisotropy parameter, $\beta(r)$. We express Q(r) as a single power-law $Q(r) \propto r^\alpha$ and derive a median slope $\alpha$ in each simulation and for each definition of Q. Our main findings are: 1. The various Q(r) profiles follow a power law to a good approximation. 2. The slopes depend on the concentration parameter c of the DM haloes, where for $c \gtrsim 10$ the slopes steepen with rising concentration and for $c \lesssim 10$ the trend flattens and even turns around. 3. The asymptotic value of $\beta$ as $r\to R_{\mathrm{vir}}$ increases with the value of c. 4. In accordance with Zait et al. 2007 $\alpha_{\mathrm{rad}}$ becomes more negative as the asymptotic value of $\beta$ at the virial radius increases. 5. This introduces a weak dependence of the $Q(r)$ slopes on the slope of the power spectrum.

We study the central dark matter (DM) cusp evolution in cosmological galactic halos. Models with and without baryons (baryons+DM, hereafter BDM model, and pure DM, PDM model, respectively) are advanced from identical initial conditions. The DM cusp properties are contrasted by a direct comparison of pure DM and baryonic models. We find a divergent evolution between the PDM and BDM models within the inner ~10 kpc region. The PDM model forms a R^{-1} cusp as expected, while the DM in the BDM model forms a larger isothermal cusp R^{-2} instead. The isothermal cusp is stable until z~1 when it gradually levels off. This leveling proceeds from inside out and the final density slope is shallower than -1 within the central 3 kpc (i.e., expected size of the R^{-1} cusp), tending to a flat core within ~2 kpc. This effect cannot be explained by a finite resolution of our code which produces only a 5% difference between the gravitationally softened force and the exact Newtonian force of point masses at 1 kpc from the center. Neither is it related to the energy feedback from stellar evolution or angular momentum transfer from the bar. Instead it can be associated with the action of DM+baryon subhalos heating up the cusp region via dynamical friction and forcing the DM in the cusp to flow out and to `cool' down. The process described here is not limited to low z and can be efficient at intermediate and even high z.

We study the distribution function (DF) of dark matter particles in haloes of mass range 10^{14}--10^{15}\Msun. In the numerical part of this work we measure the DF for a sample of relaxed haloes formed in the simulation of a standard \LambdaCDM model. The DF is expressed as a function of energy E and the absolute value of the angular momentum L, a form suitable for comparison with theoretical models. By proper scaling we obtain the results that do not depend on the virial mass of the haloes. We demonstrate that the DF can be separated into energy and angular momentum components and propose a phenomenological model of the DF in the form f_{E}(E)[1+L^{2}/(2L_{0}^{2})]^{-\beta_{\infty}+\beta_{0}}L^{-2\beta_{0}}. This formulation involves three parameters describing the anisotropy profile in terms of its asymptotic values (\beta_{0} and \beta_{\infty}) and the scale of transition between them (L_{0}). The energy part f_{E}(E) is obtained via inversion of the integral for spatial density. We provide a straightforward numerical scheme for this procedure as well as a simple analytical approximation for a typical halo formed in the simulation. The DF model is extensively compared with the simulations: using the model parameters obtained from fitting the anisotropy profile, we recover the DF from the simulation as well as the profiles of the dispersion and kurtosis of radial and tangential velocities. Finally, we show that our DF model reproduces the power-law behaviour of phase space density Q=\rho(r)/\sigma^{3}(r).

We explore the dependence of the radial alignment of subhalos on the mass of the host halo they orbit in. As the effect is seen on a broad range of scales including massive clusters as well as galactic systems it only appears natural to explore this phenomenon by means of cosmological simulations covering the same range in masses. We have 25 well resolved host dark matter halos at our disposal ranging from 10^15 Msun/h down to 10^12 Msun/h each consisting of order of a couple of million particles within the virial radius. We observe that subhalos tend to be more spherical than isolated objects. Both the distributions of sphericity and triaxiality of subhalos are Gaussian distributed with peak values of s approx. 0.80 and T approx. 0.56, irrespective of host mass. Interestingly we note that the radial alignment is independent of host halo mass and the distribution of \cos\theta (i.e. the angle between the major axis E_a of each subhalo and the radius vector of the subhalo in the reference frame of the host) is well fitted by a simple power law P(\cos\theta) proportional to \cos^4\theta with the same fitting parameters for all host halos.

To study the local Hubble flow, we have run constrained dark matter (DM) simulations of the Local Group (LG) in the concordance LCDM and OCDM cosmologies, with identical cosmological parameters apart from the Lambda term. The simulations were performed within a computational box of 64 h^{-1}Mpc centred on the LG. The initial conditions were constrained by the observed peculiar velocities of galaxies and positions of X-ray nearby clusters of galaxies. The simulations faithfully reproduce the nearby large scale structure, and in particular the Local Supercluster and the Virgo cluster. LG-like objects have been selected from the DM halos so as to closely resemble the dynamical properties of the LG. Both the LCDM and OCDM simulations show very similar local Hubble flow around the LG-like objects. It follows that, contrary to recent statements, the dark energy (DE) does not manifest itself in the local dynamics.

Dark matter (DM) halos formed in CDM cosmologies seem to be characterized by a power law phase-space density profile. The density of the DM halos is often fitted by the NFW profile but a better fit is provided by the Sersic fitting formula. These relations are empirically derived from cosmological simulations of structure formation but have not yet been explained on a first principle basis. Here we solve the Jeans equation under the assumption of a spherical DM halo in dynamical equilibrium, that obeys a power law phase space density and either the NFW-like or the Sersic density profile. We then calculate the velocity anisotropy, beta(r), analytically. Our main result is that for the NFW-like profile the beta - gamma relation is not a linear one (where gamma is the logarithmic derivative of the density rho[r]). The shape of beta(r) depends mostly on the ratio of the gravitational to kinetic energy within the NFW scale radius R_s. For the Sersic profile a linear beta - gamma relation is recovered, and in particular for the Sersic index of n = 6.0 case the linear fit of Hansen & Moore is reproduced. Our main result is that the phase-space density power law, the Sersic density form and the linear beta - gamma dependence constitute a consistent set of relations which obey the spherical Jeans equation and as such provide the framework for the dynamical modeling of DM halos.

We study the distinct effects of Dark Matter and Dark Energy on the future evolution of nearby large scale structures using constrained N-body simulations. We contrast a model of Cold Dark Matter and a Cosmological Constant (LCDM) with an Open CDM (OCDM) model with the same matter density Omega_m =0.3 and the same Hubble constant h=0.7. Already by the time the scale factor increased by a factor of 6 (29 Gyr from now in LCDM; 78 Gyr from now in OCDM) the comoving position of the Local Group is frozen. Well before that epoch the two most massive members of the Local Group, the Milky Way and Andromeda, will merge. However, as the expansion rates of the scale factor in the two models are different, the Local Group will be receding in physical coordinates from Virgo exponentially in a LCDM model and at a roughly constant velocity in an OCDM model. More generally, in comoving coordinates the future large scale structure will look like a sharpened image of the present structure: the skeleton of the cosmic web will remain the same, but clusters will be more `isolated' and the filaments will become thinner. This implies that the long-term fate of large scale structure as seen in comoving coordinates is determined primarily by the matter density. We conclude that although the LCDM model is accelerating at present due to its Dark Energy component while the OCDM model is non accelerating, their large scale structure in the future will look very similar in comoving coordinates.

The evolution of the phase-space density profile in dark matter (DM) halos is investigated by means of constrained simulations, designed to control the merging history of a given DM halo. Halos evolve through a series of quiescent phases of a slow accretion intermitted by violent events of major mergers. In the quiescent phases the density of the halo closely follows the NFW profile and the phase-space density profile, Q(r), is given by the Taylor & Navarro power law, r^{-beta}, where beta ~ 1.9 and stays remarkably stable over the Hubble time. Expressing the phase-space density by the NFW parameters, Q(r)=Qs (r/Rs)^{-beta}, the evolution of Q is determined by Qs. We have found that the effective mass surface density within Rs, Sigma_s = rhos Rs, remains constant throughout the evolution of a given DM halo along the main branch of its merging tree. This invariance entails that Qs ~ Rs^{-5/2} and Q(r) ~ Sigma_s^{-1/2} Rs^{-5/2} (r/ Rs)^{-beta}. It follows that the phase-space density remains constant, in the sense of Qs=const., in the quiescent phases and it decreases as Rs^{-5/2} in the violent ones. The physical origin of the NFW density profile and the phase-space density power law is still unknown. Yet, the numerical experiments show that halos recover these relations after the violent phases. The major mergers drive Rs to increase and Qs to decrease discontinuously while keeping Qs Rs^{5/2} = const. The virial equilibrium in the quiescent phases implies that a DM halos evolves along a sequence of NFW profiles with constant energy per unit volume (i.e., pressure) within Rs.

On the basis of a large scale 'adiabatic', namely non-radiative and non-dissipative, cosmological smooth particle hydrodynamic simulation we compare the entropy profiles of the gas and the dark matter (DM) in galaxy clusters. The quantity K_g = T_g \rho_g^{-2/3} provides a measure for the entropy of the intra-cluster gas. By analogy with the thermodynamic variables of the gas the velocity dispersion of the DM is associated with a formal temperature and thereby K_DM = \sigma_DM^2 \rho_DM^{-2/3} is defined. This DM entropy is related to the DM phase space density by K_DM \propto Q_DM^{-2/3}. In accord with other studies the DM phase space density follows a power law behaviour, Q_DM \propto r^{-1.82}, which corresponds to K_DM \propto r^{1.21}. The simulated intra-cluster gas has a flat entropy core within (0.8 \pm 0.4) R_s, where R_s is the NFW scale radius. The outer profile follows the DM behaviour, K_g \propto r^{1.21}, in close agreement with X-ray observations. Upon scaling the DM and gas densities by their mean cosmological values we find that outside the entropy core a constant ratio of K_g / K_{DM} = 0.71 \pm 0.18 prevails. By extending the definition of the gas temperature to include also the bulk kinetic energy the ratio of the DM and gas extended entropy is found to be unity for r > 0.8 R_s. The constant ratio of the gas thermal entropy to that of the DM implies that observations of the intra-cluster gas can provide an almost direct probe of the DM.

It has been recently claimed that there is no correlation between the distribution of galaxies and their peculiar velocities within the Local Volume (LV), namely a sphere of R=7/h Mpc around the Local Group (LG). It has been then stated that this implies that either locally dark matter is not distributed in the same way as luminous matter, or peculiar velocities are not due to fluctuations in mass. To test that statement a set of constrained N-body cosmological simulations, designed to reproduce the main observed large scale structure, have been analyzed. The simulations were performed within the flat-Lambda, open and flat matter only CDM cosmogonies. Two unconstrained simulations of the flat-Lambda and open CDM models were performed for comparison. LG-like objects have been selected so as to mimic the real LG environment. The local gravitational field due to all halos found within each LV is compared with the exact gravitational field induced by all matter in the simulation. We conclude that there is no correlation between the exact and the local gravitational field obtained by pairwise newtonian forces between halos. Moreover, the local gravitational field is uncorrelated with the peculiar velocities of halos. The exact gravitational field has a linear correlation with peculiar velocities but the proportionality constant relating the velocity with gravitational field falls below the prediction of the linear theory. Upon considering all matter inside the LVs, the exact and local gravitational accelerations show a much better correlation, but with a considerable scatter independent on the cosmological models. The main conclusion is that the lack of correlation between the local gravitation and the peculiar velocity fields around LG-like objects is naturally expected in the CDM cosmologies.

We investigate the effect of an assembly history on the evolution of dark matter (DM) halos of 10^{12} Msun/h using Constrained Realizations of random Gaussian fields. Five different realizations of a DM halo with distinct merging histories were constructed and evolved. Our main results are: A halo evolves via a sequence of quiescent phases of a slow mass accretion intermitted by violent episodes of major mergers. In the quiescent phases, the density is well fitted by an NFW profile, the inner scale radius Rs and the mass enclosed within it remain constant, and the virial radius (Rvir) grows linearly with the expansion parameter "a". Within each quiescent phase the concentration parameter ("c") scales as "a", and the mass accretion history (Mvir) is well described by the Tasitsiomi etal. fitting formula. In the violent phases the halos are not in a virial equilibrium and both Rs and Rvir grow discontinuously. The violent episodes drive the halos from one NFW dynamical equilibrium to another. The final structure of a halo, including "c", depends on the degree of violence of the major mergers and on their number. Next, we find a distinct difference between the behavior of various NFW parameters taken as averages over an ensemble of halos and those of individual halos. Moreover, the simple scaling relations c--Mvir do not apply to the entire evolution of individual halos, and so is the common notion that late forming halos are less concentrated than early forming ones. The entire evolution of the halo cannot be fitted by single analytical expressions.

We test the Secondary Infall Model (SIM) by direct comparison with the results of N-body simulations. Eight cluster-size and six galactic-size dark matter haloes have been selected at $z=0$ and re-simulated with high resolution. Based on their density profiles at the initial redshift, we compute their evolution by the SIM, assuming a simple prescription for the angular momentum. A comparison of the density profiles obtained by the SIM and the numerical experiments at $z=5$, 1 and 0 shows that, for most of the haloes at most epochs, the SIM reproduces the simulated mater distribution with a typical fractional deviation of less than 40 per cent over more than six order of magnitudes in the density. It is also found that, within the SIM framework, most of the diversity in the shape of the density profiles at $z=0$ arises from the scatter in the primordial initial conditions rather than the scatter in the angular momentum distribution. A crude optimization shows that a similar degree of agreement is obtained for galactic and cluster haloes, but the former seem to require slightly higher amounts of angular momentum. Our main conclusion is that the SIM provides a viable dynamical model for predicting the structure and evolution of the density profile of dark matter haloes.

We extend to large scales a method proposed in previous work that reconstructs non-parametrically the primordial power spectrum from cosmic microwave background data at high resolution. The improvement is necessary to account for the non-gaussianity of the Wilkinson Microwave Anisotropy Probe (WMAP) likelihood due primarily to cosmic variance. We assume the concordance LambdaCDM cosmology, utilise a smoothing prior and perform Monte Carlo simulations around an initial power spectrum that is scale-free and with spectral index n_s=0.97, very close to the concordance spectrum. The horizon scale for the model we are considering corresponds to the wavenumber k_h=4.52X10^{-4} Mpc^{-1}. We find some evidence for the presence of features and we quantify the probabilities of exceeding the observed deviations in WMAP data with respect to the fiducial models. We detect the following marginal departures from a scale-free (spectral index n_s=0.97) initial spectrum: a cut-off at 0.0001<k<0.001 Mpc^{-1} at 79.5% (92%), a dip at 0.001<k<0.003 Mpc^{-1} at 87.2% (98%) and a bump at 0.003<k<0.004 Mpc^{-1} at 90.3% (55.5%) confidence level. These frequentist confidence levels are calculated by integrating over the distribution of the Monte Carlo reconstructions built around the fiducial models. The frequentist analysis finds the low k cutoff of the estimated power spectrum to be about 2.5 sigma away from the n_s=0.97 model, while in the Bayesian analysis the model is about 1.5 sigma away from the estimated spectrum. (The sigma's are different for the two different methods.)

A large scale SPH+N-body simulation (GADGET) of the concordance LCDM universe is used to investigate orientation and angular momentum of galaxy clusters at z=0 in connection with their recent accretion histories. The basic cluster sample comprises the 3000 most massive friends-of-friends halos found in the 500 Mpc/h simulation box. Two disjoint sub-samples are constructed, using the mass ratio of the two most massive progenitors at z=0.5 m_2 / m_1 (m_1 < m_2), namely a recent major merger sample and a steady accretion mode sample. The mass of clusters in the merger sample is on average ~43% larger than the mass of the two progenitors (m_1 + m_2), whereas in the steady accretion mode sample a smaller increase of ~25% is found. The separation vector connecting the two most massive progenitor halos at z=0.5 is strongly correlated with the orientation of the cluster at z=0. The angular momentum of the clusters in the recent major merger sample tends to be parallel to orbital angular momentum of the two progenitors, whereas the angular momentum of the steady accretion mode sample is mainly determined by the angular momentum of the most massive progenitor. The long range correlations for the major and the minor principal axes of cluster pairs extend to distances of ~100 Mpc/h. Weak angular momentum correlations are found for distances < 20 Mpc/h. Within these ranges the major axes tend to be aligned with the connecting line of the cluster pairs whereas minor axes and angular momenta tend to be perpendicular to this line. A separate analysis of the two sub-samples reveals that the long range correlations are independent of the mass accretion mode. Thus orientation and angular momentum of galaxy clusters is mainly determined by the accretion along the filaments independently of the particular accretion mode.

The formation and structure of dark matter halos is studied by constrained simulations. A series of experiments of the formation of a 10^12 Msun/h halo is designed to study the dependence of the density profile on its merging history. We find that the halo growth consist of several quiescent phases intermitted by violent events, with the density well approximated by the NFW profile during the former phases. We find that (1) the NFW scale radius R_s stays constant during the quiescent phase and grows abruptly during the violent one. In contrast, the virial radius grows linearly during the quiescent and abruptly during the violent phases. (2) The central density stays unchanged during the quiescent phase while dropping abruptly during the violent phase, and it does not reflect the formation time of the halo. (3) The clear separation of the evolution of an individual halo into quiescent and violent phases implies that its entire evolution cannot be fitted by simple scaling relations.

The formation and structure of dark matter (DM) halos is studied by means of constrained realizations of Gaussian fields using N-body simulations. A series of experiments of the formation of a 10^{12} Msun halo is designed to study the dependence of the density profile on its merging history. We confirm that the halo growth consists of violent and quiescent phases, with the density well approximated by the Navarro-Frenk-White (NFW) profile during the latter phases. We find that (1) the NFW scale radius R_s stays constant during the quiescent phase and grows abruptly during the violent one. In contrast, the virial radius grows linearly during the quiescent and abruptly during the violent phases. (2) The central density stays unchanged during the quiescent phase while dropping abruptly during the violent phase. (3) The value of \rs reflects the violent merging history of the halo, and depends on the number of violent events and their fractional magnitudes, independent of the time and order of these events. It does not reflect the formation time of the halo. (4) The fractional change in R_s is a nonlinear function of the fractional absorbed kinetic energy within R_s in a violent event.

Sand, Treu, & Ellis (2002) have measured the central density profile of cluster MS2137-23 with gravitational lensing and velocity dispersion and removed the stellar contribution with a reasonable M/L. The resulting dark matter distribution within r<50 kpc was fitted by a density cusp of r^{-beta} with beta=0.35. This stands in an apparent contradiction to the CDM prediction of beta~1, and the disagreement worsens if adiabatic compression of the dark matter by the infalling baryons is considered. Following El-Zant, Shlosman & Hoffman (2001), we argue that dynamical friction acting on galaxies moving within the dark matter background counters the effect of adiabatic compression by transfering the orbital energy of galaxies to the dark matter, thus heating up and softening the central density cusp. Using N-body simulations of massive solid clumps moving in clusters we show that indeed the inner dark matter distribution flattens (with beta approx 0.35 for a cluster like MS2137-23) when the galaxies spiral inward. We find as a robust result that while the dark matter distribution becomes core-like, the overall mass distribution preserves its cuspy nature, in agreement with X-ray and lensing observations of clusters.

We use $N$-body simulations to measure mass functions in flat cosmological models with quintessence characterized by constant $w$ with $w=-1$, -2/3 and -1/2. The results are compared to the predictions of the formula proposed by Jenkins et al. at different redshifts, in terms of FOF masses as well as Abell masses appropriate for direct comparison to observations. The formula reproduces quite well the mass functions of simulated haloes in models with quintessence. We use the cluster mass function data at a number of redshifts from Carlberg et al. to constrain $\Omega_0$, $\sigma_8$ and $w$. The best fit is obtained in the limit $w \to 0$, but none of the values of $w$ in the considered range $-1 \le w < 0$ can actually be excluded. However, the adopted value of $w$ affects significantly the constraints in the $\Omega_0-\sigma_8$ plane. Taking into account the dependence on $w$ we find $\Omega_0=0.32 \pm 0.15$ and $\sigma_8=0.85_{-0.12}^{+0.38}$ (68% c.l.). Since less negative $w$ push the confidence regions toward higher $\Omega_0$ and lower $\sigma_8$ we conclude that relaxing the assumption of $w=-1$ typically made in such comparisons may resolve the discrepancy between recent cluster mass function results (yielding rather low $\Omega_0$ and high $\sigma_8$) and most other estimates. The fact that high $w$ values are preferred may however also point towards some unknown systematics in the data or the model with constant $w$ being inadequate.

Using high resolution N-body simulations we address the problem of emptiness of giant 20 Mpc/h diameter voids found in the distribution of bright galaxies. Are the voids filled by dwarf galaxies? Do cosmological models predict too many small dark matter haloes inside the voids? Can the problems of cosmological models on small scales be addressed by studying the abundance of dwarf galaxies inside voids? We find that voids in the distribution of 10^12 Msun/h haloes (expected galactic magnitudes ~ M_*) are almost the same as the voids in 10^11 Msun/h haloes. Yet, much smaller haloes with masses 10^9 Msun/h and circular velocities v_circ about 20 km/s readily fill the voids: there should be almost 1000 of these haloes in a 20 Mpc/h void. A typical void of diameter 20 Mpc/h contains about 50 haloes with v_circ > 50 km/s. The haloes are arranged in a pattern, which looks like a miniature Universe: it has the same structural elements as the large-scale structure of the galactic distribution of the Universe. There are filaments and voids; larger haloes are at the intersections of filaments. The only difference is that all masses are four orders of magnitude smaller. There is severe (anti)bias in the distribution of haloes, which depends on halo mass and on the distance from the centre of the void. Large haloes are more antibiased and have a tendency to form close to void boundaries. The mass function of haloes in voids is different from the ``normal'' mass function. It is much steeper for high masses resulting in very few M33-type galaxies (v_circ about 100 km/s). We present an analytical approximation for the mass function of haloes in voids.

The theory of structure formation predicts that galaxies form within extended massive dark matter halos built from smaller pieces that collided and merged, resulting in the hierarchy of galaxies, groups, and clusters observed today. Here we present constrained numerical simulations designed to match the observed local universe as well as possible and to study the formation, evolution and present day properties of such dark matter halos in different environments. Simulations have been done with different versions of the Adaptive Refinement Tree (ART) N-body code, some of them including eulerian gasdynamics.

We generalize the spherical collapse model for the formation of bound objects to apply in a Universe with arbitrary positive cosmological constant. We calculate the critical condition for collapse of an overdense region and give exact values of the characteristic densities and redshifts of its evolution.

We apply a new algorithm, called the Unbiased Minimal Variance (hereafter UMV) estimator, to reconstruct the cosmic density and peculiar velocity fields in our local universe from the SEcat catalog of peculiar velocities comprising both early (ENEAR) and late type (SFI) galaxies. The reconstructed fields are compared with those predicted from the IRAS PSCz galaxy redshift survey to constrain the value of beta = Omega_m^{0.6}/b, where Omega_m and b are the mass density and the bias parameters. The comparison of the density and velocity fields is carried out within the same methodological framework, and leads, for the first time, to consistent values of beta, yielding beta = 0.57_{-0.13}^{+0.11} and beta = 0.51 +/- {0.06}, respectively. We find that the distribution of the density and velocity residuals, relative to their respective errors, is consistent with a Gaussian distribution with sigma approximately 1, indicating that the density field predicted from the PSCz is an acceptable fit to that deduced from the peculiar velocities of the SEcat galaxies.

We evaluate the success of linear tidal-torque theory (TTT) in predicting galactic-halo spin using a cosmological N-body simulation with thousands of well-resolved haloes. The proto-haloes are identified by tracing today's haloes back to the initial conditions. The TTT predictions for the proto-haloes match, on average, the spin amplitudes of today's virialized haloes if linear growth is assumed until ~t_0/3, or 55-70 per cent of the halo effective turn-around time. This makes it a useful qualitative tool for understanding certain average properties of galaxies, such as total spin and angular-momentum distribution within haloes, but with a random scatter of the order of the signal itself. Non-linear changes in spin direction cause a mean error of ~50 degrees in the TTT prediction at t_0, such that the linear spatial correlations of spins on scales >~ 1 Mpc/h are significantly weakened by non-linear effects. This questions the usefulness of TTT for predicting intrinsic alignments in the context of gravitational lensing. We find that the standard approximations made in TTT, including a second order expansion of the Zel'dovich potential and a smoothing of the tidal field, provide close-to-optimal results.

We investigate the cross-talk between the two key components of tidal-torque theory, the inertia (I) and shear (T) tensors, using a cosmological N-body simulation with thousands of well-resolved haloes. We find that the principal axes of I and T are strongly aligned, even though I characterizes the proto-halo locally while T is determined by the large-scale structure. Thus, the resultant galactic spin, which plays a key role in galaxy formation, is only a residual due to ~10 per cent deviations from perfect alignment of T and I. The T-I correlation induces a weak tendency for the proto-halo spin to be perpendicular to the major axes of T and I, but this correlation is erased by non-linear effects at late times, making the observed spins poor indicators of the initial shear field. However, the T-I correlation implies that the shear tensor can be used for identifying the positions and boundaries of proto-haloes in cosmological initial conditions -- a missing piece in galaxy formation theory. The typical configuration is of a prolate proto-halo lying perpendicular to a large-scale high-density ridge, with the surrounding voids inducing compression along the major and intermediate inertia axes of the proto-halo. This leads to a transient sub-halo filament along the large-scale ridge, whose sub-clumps then flow along the filament and merge into the final halo. The centres of proto-haloes tend to lie in ~1 sigma over-density regions, but their association with linear density maxima smoothed on galactic scales is vague: only ~40 per cent of the proto-haloes contain peaks within them. Several other characteristics distinguish proto-haloes from density peaks, e.g., they tend to compress along two principal axes while many peaks compress along three axes.

A general method of deprojecting two-dimensional images to reconstruct the three dimensional structure of the projected object --specifically X-ray, Sunyaev-Zel'dovich (SZ) and gravitational lensing maps of rich clusters of galaxies -- assuming axial symmetry (Zaroubi et. al. 98), is considered. Here we test the applicability of the method for realistic, numerically simulated galaxy clusters, viewed from three orthogonal projections at four redshift outputs. We demonstrate that the assumption of axial symmetry is a good approximation for the 3D structure in this ensemble of galaxy clusters. Applying the method, we demonstrate that a unique determination of the cluster inclination angle is possible from comparison between the SZ and X-ray images and, independently, between SZ and surface density maps. Moreover, the results from these comparisons are found to be consistent with each other and with the full 3D structure inclination angle determination. The radial dark matter and gas density profiles as calculated from the actual and reconstructed 3D distributions show a very good agreement. The method is also shown to provide a direct determination of the baryon fraction in clusters, independent of the cluster inclination angle.

We present results of gasdynamics+N-body constrained cosmological simulations of the Local Supercluster region (LSC; about 30/h Mpc around the Virgo cluster), which closely mimic the real Universe within 100 Mpc by imposing constraints from Mark III catalog of galaxy peculiar velocities. The simulations are used to study the properties and possible observational signatures of intergalactic medium in the LSC region. We find that, in agreement with previous unconstrained simulations, about 30% of the gas in this region is in the warm/hot phase at T~10^5-10^7 K, and about 40% in the diffuse phase at T<10^5 K in low-density regions. The best prospects for detection of the warm/hot intergalactic medium of the LSC located in filaments and in the vicinity of virialized regions of groups and clusters are through absorption in resonant lines of OVII and OVIII in soft X-rays and in the OVI doublet in UV. If intergalactic gas in filaments (rho/<rho>~ 1-10) is enriched to typical metallicities of >0.05, the column densities of OVI, OVII, and OVIII along a random line of sight near the North Galactic Pole, especially near the supergalactic plane, have a significant probability to be in the range detectable by the current (FUSE,XMM) and future (Constellation-X) instruments. (Abridged)

Coupling the dark matter (DM) and baryonic components in galaxies via dynamical friction (by allowing the latter to be clumpy) leads to energy input into the DM. The resulting expanding tendency in the DM component overcomes the competing effect of adiabatic contraction, resulting from the shrinking gas component, and leads to a nearly constant density core. The baryonic mass inflow also leads to the coupling of the structure at the kpc scale and below and that of the structure of the DM at the 10 kpc scale and above, with rate of mass inflow being dependent on the initial halo central concentration. The latter also determines the final halo core radius.

We present cosmological simulations which closely mimic the real Universe within 100Mpc of the Local Group. The simulations, called Constrained Simulations, reproduce the large-scale density field with major nearby structures, including the Local Group, the Coma and Virgo clusters, the Great Attractor, the Perseus-Pices, and the Local Supercluster, in approximately correct locations. The MARK III survey of peculiar velocities of the observed structures inside 80Mpc/h sphere is used to constrain the initial conditions. Fourier modes on scales larger then 5Mpc/h are dominated by the constraints, while small scale waves are random. The main aim of this paper is the structure of the Local Supercluster (LSC; 30Mpc/h around the Virgo cluster) and the Local Group environment. We find that at the current epoch most of the mass (7.5e14Msun/h) of the LSC is located in a filament roughly centered on the Virgo cluster and extending over 40Mpc/h. The simulated Local Group (LG) is located in an adjacent smaller filament, which is not a part of the main body of the LSC, and has a peculiar velocity of 250kms toward the Virgo cluster. The peculiar velocity field in the LSC region is complicated and is drastically different from the field assumed in the Virgocentric infall models. The peculiar velocity flow in the vicinity of the LG in the simulation is ``cold'': the peculiar line-of-sight velocity dispersion within 7Mpc/h of the LG is less than 60km/s, comparable to the observed velocity dispersion of nearby galaxies.

N-body simulations and analytical calculations of the gravitational collapse in an expanding universe predict that halos should form with a diverging inner density profile, the cusp. There are some observational indications that the dark matter distribution in galaxies might be characterized by a finite core. This `core catastrophe' has prompted a search for alternatives to the CDM cosmogony. It is shown here that the discrepancy between theory and observations can be very naturally resolved within the standard CDM model, provided that gas is not initially smoothly distributed in the dark matter halo, but rather is concentrated in clumps of mass $\geq 0.01 %$ the total mass of the system. Dynamical friction acting on these lumps moving in the background of the dark matter particles, dissipates the clumps orbital energy and deposits it in the dark matter. Using Monte-Carlo simulations, it is shown that the dynamical friction provides a strong enough drag, and that with realistic baryonic mass fractions, the available orbital energy of the clumps is sufficient to heat the halo and turn the primordial cusp into a finite, non-diverging core --- overcoming the competing effect of adiabatic contraction due to gravitational influence of the shrinking baryonic component. Depending on the initial conditions, the total density distribution may either become more or less centrally concentrated. Possible consequences of the proposed mechanism for other problems in the CDM model and for the formation and early evolution of the baryonic component of galaxies are also briefly discussed.

We present a method for decomposing the cosmological velocity field in a given volume into its divergent component due to the density fluctuations inside the volume, and its tidal component due to the matter distribution outside the volume. The input consists of the density and velocity fields that are reconstructed either by POTENT or by Wiener Filter from a survey of peculiar velocities. The tidal field is further decomposed into a bulk velocity and a shear field. The method is applied here to the Mark III data within a sphere of radius 60 Mpc/h about the Local Group, and to the SFI data for comparison. We find that the tidal field contributes about half of the Local-Group velocity with respect to the CMB, with the tidal bulk velocity pointing to within ~ 30 degrees of the CMB dipole. The eigenvector with the largest eigenvalue of the shear tensor is aligned with the tidal bulk velocity to within ~ 40 degrees. The tidal field thus indicates the important dynamical role of a super attractor of mass (2-5) x 10^17 M_sun/h Omega^0.4 at ~ 150 Mpc/h, coinciding with the Shapley Concentration. There is also a hint for the dynamical role of two big voids in the Supergalactic Plane. The results are consistent for the two data sets and the two methods of reconstruction.

We generalize the spherical collapse model for the formation of dark matter halos to apply in a universe with arbitrary positive cosmological constant. We calculate the critical condition for collapse of an overdense region and give exact values of the characteristic densities and redshifts of its evolution.

The uncovering of the large-scale structure that is hidden in the Zone of Avoidance is essential to our quest for the understanding of the dynamics of the nearby universe. Dedicated and sophisticated observations of the Zone of Avoidance have to be matched by equally sophisticated theoretical tools for the reconstruction of the large-scale structure from sparse, noisy and very incomplete data. A general approach for the estimation of cosmological parameters and the reconstruction of the density and velocity fields is presented here, based on maximum likelihood analysis, Wiener filtering and constrained realizations. Radial velocity surveys and the IRAS 1.2 Jy redshift survey are used to reconstruct the Zone of Avoidance.

Using the formalism of the spherical infall model the structure of collapsed and virialized dark halos is calculated for a variety of scale-free initial conditions. In spite of the scale-free cosmological nature of the problem, the collapse of individual objects is not self-similar. Unlike most of previous calculations the dynamics used here relies only on adiabatic invariants and not on self-similarity. The paper focuses on the structure of the innermost part of the collapsed halos and addresses the problem of central density cusps. The slopes of density profiles at 1% of virial radius are calculated for a variety of cosmological models and are found to vary with the mass of the halos and power spectrum of the initial conditions. The inner slopes range between r^-2.3 and r^-2 with the limiting case of r^-2 reached for the largest masses. The steep cusps found here correspond to the limiting case where all particles move on radial orbits. The introduction of angular momentum will make the density profile shallower. We expect this to resolve the discrepancy found between the calculated profiles and the ones found in high resolution N-body simulations, where the exponent ranges from -0.5 to -1.5. The robust prediction here is that collisionless gravitational collapse in an expanding universe is expected to form density cups and not halos with a core structure.

We estimate the mass density fluctuations power spectrum (PS) on large scales by applying a maximum likelihood technique to the peculiarvelocity data of the recently completed redshift-distance survey of early-type galaxies (ENEAR). The general results are in agreement with the high amplitude power spectra found from similar analysis of other independent all-sky catalogs of peculiar velocity data such as MARK III and SFI. For Lambda & Open CDM COBE normalized PS models, the best-fit parameters are confined by a contour approximately defined by Omega h^{1.3}=0.377+-0.08 and Omega h^{0.88}=0.517+-0.083, respectively. Gamma-shape models, free of COBE normalization, resultsin the weak constraint of $\Gamma \geq 0.17$ and in the rather stringent constraint of sigma_8 Omega^{0.6}=1.0+-0.25. All quoted uncertainties refer to 3-sigma confidence-level. The calculated PS is used as a prior for Wiener reconstruction of the density field at different resolutions and the three-dimensional velocity field within a volume of radius ~80 Mpc/h. All major structures in the nearby universe are recovered and are well matched to those predicted from all-sky redshift surveys.

Using eigenmode expansion of the Mark-3 and SFI surveys of cosmological radial velocities a goodness-of-fit analysis is applied on a mode-by-mode basis. This differential analysis complements theBayesian maximum likelihood analysis that finds the most probable model given the data. Analyzing the surveys with their corresponding most likely models from the CMB-like family of models, as well as with the currently popular Lambda-CDM model, reveals a systematic inconsistency of the data with these `best' models. There is a systematic trend of the cumulative chi^2 to increase with the mode number (where the modes are sorted by decreasing order of the eigenvalues). This corresponds to a decrease of the chi^2 with the variance associated with a mode, and hence with its effective scale. It follows that the differential analysis finds that on small (large) scales the global analysis of all the modes `puts' less (more) power than actually required by the data. This observed trend might indicate one of the followings: a. The theoretical model (i.e. power spectrum) or the error model (or both) have an excess of power on large scales; b. Velocity bias; c. The velocity data suffers from still uncorrected systematic errors.

Radial velocity surveys are examined in terms of eigenmode analysis within the framework of CDM-like family of models. Rich surveys such as MARK III and SFI, which consist of more than $10^{3}$ radial velocities, are found to have a few tens of modes that are not noise dominated. Poor surveys, which have only a few tens of radial velocities, are noise dominated across the eigenmode spectrum. In particular, the bulk velocity of such surveys has been found to be dominated by the more noisy modes. The MARK III and SFI are well fitted by a tilted flat CDM model found by a maximum likelihood analysis and a $\chi^{2}$ statistics . However, a mode-by-mode inspection shows that a substantial fraction of the modes lie outside the $90% $ confidence level. This implies that although globally the CDM-like family of models seems to be consistent with radial velocity surveys, in detail it does not. This might indicate a need for a revised power spectrum or for some non-trivial biasing scheme.

Unlike the substantial coherent bulk motion in which our local patch of the Cosmos is participating, the amplitude of the random motions around this large scale flow seems to be surprisingly low. Attempts to invoke global explanations to account for this coldness of the local cosmic velocity field have not yet been succesfull. Here we propose a different view on this cosmic dilemma, stressing the repercussions of our cosmic neighbourhood embodying a rather uncharacteristic region of the Cosmos. Suspended between two huge mass concentrations, the Great Attractor region and the Perseus-Pisces chain, we find ourselves in a region of relatively low density yet with a very strong tidal shear. By means of constrained realizations of our local Universe, based on Wiener-filtered reconstructions inferred from the Mark III catalogue of galaxy peculiar velocities, we show that indeed this configuration may induce locally cold regions. Hence, the coldness of the local flow may be a cosmic variance effect.

Several studies have indicated that the local cosmic velocity field is rather cold, in particular in the regions outside the massive, virialized clusters of galaxies. If our local cosmic environment is taken to be a representative volume of the Universe, the repercussion of this finding is that either we live in a low-$\Omega$ Universe and/or that the galaxy distribution is a biased reflection of the underlying mass distribution. Otherwise, the pronounced nature of the observed galaxy distribution would be irreconcilable with the relatively quiet flow of the galaxies. Here we propose a different view on this cosmic dilemma, stressing the fact that our cosmic neighbourhood embodies a region of rather particular dynamical properties, and henceforth we are apt to infer flawed conclusions with respect to the global Universe. Suspended between two huge mass concentrations, the Great Attractor region and the Perseus-Pisces chain, we find ourselves in a region of relatively low density yet with a very strong tidal shear. This tidal field induces a local velocity field with a significant large-scale bulk flow but a low small-scale velocity dispersion. By means of constrained realizations of our local Universe, consisting of Wiener-filtered reconstructions inferred from the Mark III catalogue of galaxy peculiar velocities in combination with appropriate spectrally determined fluctuations, we study the implications for our local velocity field. We find that we live near a local peak in the distribution of the cosmic Mach number, $|v_{bulk}|/\sigma_v$, and that our local cosmic niche is located in the tail of the Mach number distribution function.

We present an alternative, Bayesian method for large-scale reconstruction from observed peculiar velocity data. The method stresses a rigorous treatment of the random errors and it allows extrapolation into poorly sampled regions in real space or in k-space. A likelihood analysis is used to determine the fluctuation power spectrum, followed by a Wiener Filter (WF) analysis to obtain the minimum-variance mean fields of velocity and mass density. Constrained Realizations (CR) are then used to sample the statistical scatter about the WF mean field. The WF/CR method is applied as a demonstration to the Mark III data with 1200 km/s, 900 km/s, and 500 km/s resolutions. The main reconstructed structures are consistent with those extracted by the POTENT method. A comparison with the structures in the distribution of IRAS 1.2Jy galaxies yields a general agreement. The reconstructed velocity field is decomposed into its divergent and tidal components relative to a cube of +/-8000 km/s centered on the Local Group. The divergent component is very similar to the velocity field predicted from the distribution of IRAS galaxies. The tidal component is dominated by a bulk flow of 194 +/- 32 km/s towards the general direction of the Shapley concentration, and it also indicates a significant quadrupole.

We observed the large and faint planetary nebula (PN) around IN Com in H-alpha and [NII] light with a coronagraphic CCD device on the Wise Observatory reflector blocking the light from the central star. Our goal was to provide a second image of the object with which to confirm the features seen in the only published photograph from the paper reporting the discovery of this object. The nebula is extremely faint, but a combination of images totalling ~one and a half hours of exposure shows it fairly well. A novel image processing algorithm has been applied to the noisy image in order to reveal faint extended details of the images. The algorithm is based on a non-linear self-adaptive filter applied to the wavelet transform of the image. The nebula is not round or elliptical, but shows a two-lobed and possibly three-lobed morphology, as well as a peculiar hole-like feature East of the central star. There is definite East-West and slightly less definite North-South asymmetry.

We consider a general method of deprojecting 2D images to reconstruct the 3D structure of the projected object, assuming axial symmetry. The method consists of the application of the Fourier Slice Theorem to the general case where the axis of symmetry is not necessarily perpendicular to the line of sight, and is based on an extrapolation of the image Fourier transform into the so-called cone of ignorance. The method is specifically designed for the deprojection of X-ray, Sunyaev-Zeldovich (SZ) and gravitational lensing maps of rich clusters of galaxies. For known values of the Hubble constant, H0, and inclination angle, the quality of the projection depends on how exact is the extrapolation in the cone of ignorance. In the case where the axis of symmetry is perpendicular to the line of sight and the image is noise-free, the deprojection is exact. Given an assumed value of H0, the inclination angle can be found by matching the deprojected structure out of two different images of a given cluster, e.g., SZ and X-ray maps. However, this solution is degenerate with respect to its dependence on the assumed H0, and a third independent image of the given cluster is needed to determine H0 as well. The application of the deprojection algorithm to upcoming SZ, X-ray and weak lensing projected mass images of clusters will serve to determine the structure of rich clusters, the value of H0, and place constraints on the physics of the intra-cluster gas and its relation to the total mass distribution.

The linear algorithm of the Wiener filter and constrained realizations (CRs) of Gaussian random fields is extended here to perform non-linear CRs. The procedure consists of: (1) Using low resolution data to constrain a high resolution realization of the underlying field, as if the linear theory is valid; (2) Taking the linear CR backwards in time, by the linear theory, to set initial conditions for N-body simulations; (3) Forwarding the field in time by an N-body code. An intermediate step is introduced to `linearize' the low resolution data. The non-linear CR can be applied to any observational data set that is quasi-linearly related to the underlying field. Here it is applied to the IRAS 1.2Jy catalog using 846 data points within a sphere of 6000 km/s, to reconstruct the full non-linear large scale structure of our `local' universe. The method is tested against mock IRAS surveys, taken from random non-linear realizations. A detaile analysis of the reconstructed non-linear structure is presented.

The power spectrum (PS) of mass density fluctuations, independent of `biasing', is estimated from the Mark III catalog of peculiar velocities using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the probability of the model given the data. The method has been tested using detailed mock catalogs. It has been applied to generalized CDM models with and without COBE normalization. The robust result for all the models is a relatively high PS, with $P(k) \Omega^{1.2} = (4.8 \pm 1.5) \times 10^3 (Mpc/h)^3$ at $k=0.1 h/Mpc$. An extrapolation to smaller scales using the different CDM models yields $\sigma_8 \Omega^{0.6} = 0.88 \pm 0.15$. The peak is weakly constrained to the range $0.02 \leq k \leq 0.06 h/Mpc$. These results are consistent with a direct computation of the PS (Kolatt & Dekel 1996). When compared to galaxy-density surveys, the implied values for $\beta$ ($\equiv \Omega^{0.6}/b$) are of order unity to within 25%. The parameters of the COBE-normalized, flat CDM model are confined by a 90% likelihood contour of the sort $\Omega h_{50}^\mu n^\nu = 0.8 \pm 0.2$, where $\mu = 1.3$ and $\nu = 3.4, 2.0$ for models with and without tensor fluctuations respectively. For open CDM the powers are $\mu = 0.95$ and $\nu = 1.4$ (no tensor fluctuations). A $\Gamma$-shape model free of COBE normalization yields only a weak constraint: $\Gamma = 0.4 \pm 0.2$.

The mass density field as extracted from peculiar velocities in our cosmological neighborhood is mapped back in time to the CMB in two ways. First, the density power spectrum ($P_k$) is translated into a temperature angular power spectrum of sub-degree resolution ($C_l$) and compared to observations. Second, the local density field is translated into a temperature map in a patch on the last-scattering surface of a distant observer. A likelihood analysis of the Mark III peculiar velocity data have constrained the range of parameters for $P_k$ within the family of COBE-normalized CDM models (Zaroubi et al 1996), favoring a slight tilt in the initial spectrum, $n<1$. The corresponding range of $C_l$'s is plotted against current observations, indicating that the CMB data can tighten the constraints further: only models with ``small'' tilt ($n\sim 0.9$) and ``high'' baryonic content ($\Omega_b \sim 0.1$) could survive the two data sets simultaneously. The local mass density field that has been recovered from the velocities via a Wiener method is convolved with a Boltzmann calculation to recover $10'$ resolution temperature maps as viewed from different directions. The extent of the CMB patch and the amplitude of fluctuations depend on the choice of cosmological parameters, e.g., the local $100\hmpc$ sphere corresponds to $90'$ to $30'$ at the CMB for $\Omega$ between 1 and 0 respectively. The phases of the temperature map are correlated with those of the density field, contrary to the contribution of the Sachs-Wolfe effect alone. This correlation suggests the possibility of an inverse reconstruction of the underlying density field from CMB data with interesting theoretical implications.

The power spectrum (PS) of {\it mass} density fluctuations, independent of ``biasing", is estimated from the Mark3 Catalog of Peculiar Velocities of galaxies using Bayesian statistics. A parametric model is assumed for the PS, and the free parameters are determined by maximizing the posterior probability of the model given the data. The method has been tested using detailed mock catalogs. We use generalized CDM models with and without COBE normalization. The robust result for all the models tested is $P(k) \Omega^{1.2} = (4.1 \pm 0.7)x10^3 (Mpc h^{-1})^3$ at $k=0.1 h Mpc^{-1}$, with the peak constrained to the range $0.03 < k < 0.06 h Mpc^{-1}$. It is consistent with a direct computation of the PS (Kolatt \& Dekel 1996). When compared to galaxy-density surveys, the implied values for $\beta$ ($\equiv \Omega^{0.6}/b$) are of order unity to within 25\%. A $\Gamma$-shape model, free of COBE normalization, is constrained by the velocity data to $\Gamma=0.5\pm 0.15$ and $\sigma_8 \Omega^{0.6}=0.85\pm 0.1$. Within the family of COBE-normalized CDM models, the best tilted model ($\Omega=1$, $n \approx 0.84 h_{50}^{-0.65}$) and the best open model ($n=1$, $\Omega \approx 0.75 h_{50}^{-1.0}$) are more likely than the best $\Lambda$ model ($n=1$, $\Lambda=1-\Omega$, $\Omega\approx 0.70 h_{50}^{-1.2}$). The most likely CDM model with $\Omega\leq 1$ is found to be of $\Omega=1$, $h\approx 0.75$, and a {\it tilted} spectrum of $n= 0.8\pm 0.02$ with tensor fluctuations.The tightest constraint for the tilted-$\Lambda$ models is of the sort $\Omega h_{50}^{1.2} n^\nu = 0.7 \pm 0.08$, with $\nu=3.8$ and $1.85$ with and without tensor fluctuations respectively.

We apply a Wiener filter to the two-year COBE DMR data. The resulting sky map has significantly reduced noise levels compared to the raw data: the most prominent hot and cold spots are significant at the 4-sigma level. Furthermore, the entire covariance matrix of the errors in the filtered sky map is known, and it is therefore possible to make constrained realizations of the microwave sky with the correct a posteriori probability distribution. The filtered DMR sky map is used to make predictions for the Tenerife experiment. Two prominent features are predicted in a region of the sky not yet analyzed by the Tenerife group. The presence of these features is a robust prediction of the standard cosmological paradigm; if these features are not observed, some of our fundamental assumptions must be incorrect.

The linear algorithm of the the full non-linear large scale structure of Gaussian random fields is extended here to to perform non-linear CRs. The procedure consists of: (1) Using linear CR of low resolution data to construct a high resolution realization of the underlying field, as if the linear theory is valid; (2) Taking the linear CR backwards in time, by the linear theory, to set initial conditions for an N-body simulation; (3) Forwarding the field in time by an N-body code. An intermediate step might be introduced to `linearize' the low resolution data. The non-linear CR can be applied to any observational data set that is linearly related to the underlying field. Here it is applied to the IRAS 1.2Jy catalog using 843 data points within a sphere of $6000 Km/s$, to reconstruct the full non-linear large scale structure of our `local' universe.

The applicability of the highly idealized secondary infall model to `realistic' initial conditions is investigated. The collapse of proto-halos seeded by $3\sigma$ density perturbations to an Einstein--de Sitter universe is studied here for a variety of scale-free power spectra with spectral indices ranging from $n=1$ to $-2$. Initial conditions are set by the constrained realization algorithm and the dynamical evolution is calculated both analytically and numerically. The analytic calculation is based on the simple secondary infall model where spherical symmetry is assumed. A full numerical simulation is performed by a Tree N-body code where no symmetry is assumed. A hybrid calculation has been performed by using a monopole term code, where no symmetry is imposed on the particles but the force is approximated by the monopole term only. The main purpose of using such code is to suppress off-center mergers. In all cases studied here the rotation curves calculated by the two numerical codes are in agreement over most of the mass of the halos, excluding the very inner region, and these are compared with the analytically calculated ones. The main result obtained here, reinforces the foundings of many N-body experements, is that the collapse proceeds 'gently' and not {\it via} violent relaxation. There is a strong correlation of the final energy of individual particles with the initial one. In particular we find a preservation of the ranking of particles according to their binding energy. In cases where the analytic model predicts non-increasing rotation curves its predictions are confirmed by the simulations. Otherwise, sensitive dependence on initial conditions is found and the analytic model fails completely.

The formalism of Wiener filtering is developed here for the purpose of reconstructing the large scale structure of the universe from noisy, sparse and incomplete data. The method is based on a linear minimum variance solution, given data and an assumed \prior model which specifies the covariance matrix of the field to be reconstructed. While earlier applications of the Wiener filter have focused on estimation, namely suppressing the noise in the measured quantities, we extend the method here to perform both prediction and dynamical reconstruction. The Wiener filter is used to predict the values of unmeasured quantities, such as the density field in un-sampled regions of space, or to deconvolve blurred data. The method is developed, within the context of linear gravitational instability theory, to perform dynamical reconstruction of one field which is dynamically related to some other observed field. This is the case, for example, in the reconstruction of the real space galaxy distribution from its redshift distribution When the field to be reconstructed is a Gaussian random field, such as the primordial perturbation field predicted by the canonical model of cosmology, the Wiener filter can be pushed to its fullest potential. In such a case the Wiener estimator coincides with the Bayesian estimator designed to maximize the {\it posterior} probability. The Wiener filter can be also derived by assuming a quadratic regularization function, in analogy with the `Maximum Entropy' method. The mean field obtained by the minimal variance solution can be supplemented with constrained realizations of the Gaussian field to

The clustering in redshift space is studied here to first order within the framework of gravitational instability. The distortion introduced by the peculiar velocities of galaxies results in anisotropy in the galaxy distribution and mode-mode coupling when analyzed in Fourier space. An exact linear calculation of the full covariance matrix in both the real and Fourier space is presented here. The explicit dependence on $\Omeg$ and the biasing parameter is calculated and its potential use as a probe of these parameters is analyzed. It is shown that Kaiser's formalism can be applied only to a data set that subtends a small solid angle on the sky, and therefore cannot be used in the case of all sky surveys. The covariance matrix in the real space is calculated explicitly for {\it CDM} model, where the behavior along and perpendicular to the line of sight is shown.

We present a new method for recovering the cosmological density, velocity, and potential fields from all-sky redshift catalogues. The method is based on an expansion of the fields in orthogonal radial (Bessel) and angular (spherical harmonic) functions. In this coordinate system, peculiar velocities introduce a coupling of the radial harmonics describing the density field in redshift space but leave the angular modes unaffected. In the harmonic transform space, this radial coupling is described by a distortion matrix which can be computed analytically within the context of linear theory; the redshift space harmonics can then be converted to their real space values by inversion of this matrix. Statistical noise is mitigated by regularizing the matrix inversion with a Wiener filter. The method yields a minimum variance estimate of the density field in real space. In this coordinate system, the minimum variance harmonics of the peculiar velocity and potential fields are related to those of the density field by simple linear transformations. Tests of the method with simulations of a CDM universe and comparison with previously proposed methods demonstrate it to be a very promising new reconstruction method for the local density and velocity field. A first application to the 1.2 Jy IRAS redshift survey is presented.

We derive an optimal linear filter to suppress the noise from the COBE DMR sky maps for a given power spectrum. We then apply the filter to the first-year DMR data, after removing pixels within $20^\circ$ of the Galactic plane from the data. The filtered data have uncertainties 12 times smaller than the noise level of the raw data. We use the formalism of constrained realizations of Gaussian random fields to assess the uncertainty in the filtered sky maps. In addition to improving the signal-to-noise ratio of the map as a whole, these techniques allow us to recover some information about the CMB anisotropy in the missing Galactic plane region. From these maps we are able to determine which hot and cold spots in the data are statistically significant, and which may have been produced by noise. In addition, the filtered maps can be used for comparison with other experiments on similar angular scales.

The analysis of whole-sky galaxy surveys commonly suffers from the problems of shot-noise and incomplete sky coverage (e.g. at the Zone of Avoidance). The orthogonal set of spherical harmonics is utilized here to expand the observed galaxy distribution. We show that in the framework of Bayesian statistics and Gaussian random fields the $4 \pi$ harmonics can be recovered and the shot-noise can be removed, giving the most probable picture of the underlying density field. The correction factor from observed to reconstructed harmonics turns out to be the well-known Wiener filter (the ratio of signal to signal+noise), which is also derived by requiring minimum variance. We apply the method to the projected 1.2 Jy IRAS survey. The reconstruction confirms the connectivity of the Supergalactic Plane across the Galactic Plane (at Galactic longitude $l \sim 135^o$ and $l \sim 315^o$) and the Puppis cluster behind the Galactic Plane ($ l \sim 240^o$). The method can be extended to 3-D in both real and redshift space, and applied to other cosmic phenomena such as the COBE Microwave Background maps.

The large scale structure of the universe is reconstructed by means of the Wiener filter and constrained realizations of Gaussian fields. The density field constructed from optically selected galaxies (Hudson \cite{Hudson}) has been sampled out to a distance of $80\hmpc$. This data set is used for the reconstruction of the underlying primordial perturbation field. A preliminary reconstruction analysis from the IRAS $1.9Jy$ density field is presented here for the sake of comparison. The perturbation field within a cube of $160\hmpc$ on the side, centered on the Local Group, and with Gaussian smoothing of length $10\hmpc$ is reconstructed, assuming the standard CDM model. In particular, the density field at the Zone of Avoidance is studied. A preliminary maximum likelihood analysis performed on Hudson's data shows that it is consistent with a tilted CDM with $n\sim -0.5\pm 0.3$ and $\sigma_8\sim 1.4\pm 0.15$