We have shown (Colin et al., 2019) that the acceleration of the Hubble expansion rate inferred from Type Ia supernovae (SNe Ia) is, at $3.9\sigma$ significance, a dipole approximately aligned with the CMB dipole, while its monopole component, which can be interpreted as due to a Cosmological Constant or dark energy, is consistent with zero at $1.4\sigma$. This has been challenged by Rubin & Heitlauf (2019) who assert that the dipole arises because we made an incorrect assumption about the SNe Ia light-curve parameters (viz. took them to be sample- and redshift independent), and did not allow for the motion of the Solar system (w.r.t. the 'CMB frame' in which the CMB dipole supposedly vanishes). In fact what has an even larger impact on our finding is that we reversed the inconsistent "corrections" made for the peculiar velocities of the SNe Ia host galaxies w.r.t the CMB frame, which in fact serve to bias the lever arm of the Hubble diagram towards higher inferred values of the monopole. We demonstrate that even if all such corrections are made consistently and both sample- and redshift-dependence is allowed for in the standardisation of supernova light curves, the evidence for isotropic acceleration rises to just $2.8\,\sigma$. Thus the criticism of Rubin & Heitlauf serves only to highlight that "corrections" must be made to the SNe Ia data assuming the standard $\Lambda$CDM model, in order to recover it from the data.

We study correlated fluctuations of Type~Ia supernova observables due to peculiar velocities of both the observer and the supernova host galaxies, and their impact on cosmological parameter estimation. We demonstrate using the CosmicFlows-3 dataset that at low redshifts the corrections for peculiar velocities in the JLA catalogue have been systematically underestimated. By querying a horizon-size N-body simulation we find that compared to a randomly placed Copernican observer, an observer in an environment like our local universe will see 2--5 times stronger correlations between supernovae in the JLA catalogue. Hence the covariances usually employed which assume a Copernican observer underestimate the effects of coherent motion of the supernova host galaxies. Although previous studies have suggested that this should have $<2\%$ effect on cosmological parameter estimation, we find that when peculiar velocities are treated consistently the JLA data favours significantly smaller values of matter and dark energy density than in the standard $\Lambda$CDM model. A joint fit to simultaneously determine the cosmological parameters and the bulk flow finds a bulk flow faster than 200 km~s$^{-1}$ continuing beyond 200~Mpc. This demonstrates that the local bulk flow is an essential nuisance parameter which must be included in cosmological model fitting when analysing supernova data.

Using fully GPU $N$-body simulations, we demonstrate for the first time that subhalos sink and transfer energy via dynamical friction into the centres of dwarf galaxies. This dynamical heating kicks any central massive black hole (MBH) out to tens of parsecs, especially at early epochs ($z$=1.5-3). This mechanism helps explain the observed off-center BHs in dwarf galaxies and also predicts that off-center BHs are more common in higher mass dwarf galaxies since dynamical friction becomes significantly weaker and BHs take more time to sink back towards the centres of their host galaxies. One consequence of off-center BHs during early epochs of dwarf galaxies is to quench any BH feedback.

We run state-of-the-art high resolution fully GPU N-body simulations to show that M31 harbours a dark matter core. Observational data in the inner halo of M31 provide stringent constraints on the initial conditions of our simulations. We demonstrate that an infalling satellite on a highly eccentric orbit is at the origin of the giant stellar stream and the shell-like features of M31. The infalling satellite heats up the central parts of its host and triggers a cusp-core transition in the dark matter halo of M31, generating a universal core independently of the initial conditions. Our results imply that cores in massive galaxies are a common and natural feature of cold dark matter haloes that have been initially cuspy but have accreted subhaloes on highly eccentric orbits

We performed a series of high-resolution $N$-body simulations to examine whether dark matter candidates in the form of primordial black holes (PBHs) can solve the cusp-core problem in low-mass dwarf galaxies. If some fraction of the dark matter in low-mass dwarf galaxies consists of PBHs and the rest is cold dark matter, dynamical heating of the cold dark matter by the PBHs induces a cusp-to-core transition in the total dark matter profile. The mechanism works for PBHs in the 25-100 M$_{\sun}$ mass window, consistent with the LIGO detections, but requires a lower limit on the PBH mass fraction of 1$\%$ of the total dwarf galaxy dark matter content. The cusp-to-core transition time-scale is between 1 and 8 Gyr. This time-scale is also a constant multiple of the relaxation time between cold dark matter particles and PBHs, which depends on the mass, the mass fraction and the scale radius of the initial density profile of PBHs. We conclude that dark matter cores occur naturally in halos comprised of cold dark matter and PBHs, without the need to invoke baryonic processes.

We use a fully GPU $N$-body code to demonstrate that dark matter minihalos, as a new component of globular clusters, resolve both the timing and cusp-core problems in Fornax if the five (or six) globular clusters were recently accreted ($\leq$ 3 Gyr ago) by Fornax. Under these assumptions, infall of these globular clusters does not occur and no star clusters form in the centre of Fornax in accordance with observations. We find that crossings of globular clusters that have DM minihalos near the Fornax centre induce a cusp-to-core transition of the dark matter halo and hence resolve the cusp-core problem in this dwarf galaxy. The dark matter core size depends on the frequency of globular cluster crossings. Our simulations clearly demonstrate also that between the passages, the dark matter halo can regenerate its cusp. Moreover, our models are in good agreement with constraints on the dark matter masses of globular clusters as our clusters lose a large fraction of their initial dark matter minihalos. These results provide circumstantial evidence for the universal existence of dark matter halos in globular clusters.

Observations reveal a `bulk flow' in the local Universe which is faster and extends to much larger scales than is expected around a typical observer in the standard $\Lambda$CDM cosmology. This is expected to result in a scale-dependent dipolar modulation of the acceleration of the expansion rate inferred from observations of objects within the bulk flow. From a maximum-likelihood analysis of the Joint Lightcurve Analysis (JLA) catalogue of Type Ia supernovae we find that the deceleration parameter, in addition to a small monopole, indeed has a much bigger dipole component aligned with the CMB dipole which falls exponentially with redshift $z$: $q_0 = q_\mathrm{m} + \vec{q}_\mathrm{d}.\hat{n}\exp(-z/S)$. The best fit to data yields $q_\mathrm{d} = -8.03$ and $S = 0.0262~(\Rightarrow d \sim 100~\mathrm{Mpc})$, rejecting isotropy ($q_\mathrm{d} = 0$) with $3.9\sigma$ statistical significance, while $q_\mathrm{m} = -0.157$ and consistent with no acceleration ($q_\mathrm{m} = 0$) at $1.4\sigma$. Thus the cosmic acceleration deduced from supernovae may be an artefact of our being non-Copernican observers, rather than evidence for a dominant component of `dark energy' in the Universe.

We re-investigate the Fornax cusp-core problem using observational results on the spatial and mass distributions of globular clusters (GCs) in order to put constraints on the dark matter profile. We model Fornax using high-resolution N-body simulations with entirely live systems, i.e. self-gravitating systems composed of stars and dark matter, which account correctly for dynamical friction and tidal effects between Fornax and the globular clusters. We test two alternative hypotheses, which are a cored and a cuspy halo for Fornax by exploring a reasonable range of initial conditions on globular clusters. For Fornax cored dark matter halo, we derive a lower limit on the core size of $r_{c}\gtrsim$ 0.5 kpc. Contrary to many previous works, we show also that for different initial conditions, a cuspy halo is not ruled out in our simulations based on observations of Fornax globular clusters.

Fornax is the most massive of the Milky Way dwarf spheroidal galaxies and has five globular clusters orbiting in a dense background of dark matter. Observational analyses suggest that globular clusters were initially much more massive and lost most of their stars to the Fornax field. We re-investigate the Fornax cusp-core problem, to clarify tensions between simulations and observations concerning the dark matter halo density profile. N-body simulations predict a centrally steep power-law density profile, while observations of the globular clusters seem to prefer that the dark matter halo density is constant at the center. For the first time, we ran pure N-body simulations with entirely live systems, i.e. self-gravitating systems composed of particles (including stars and dark matter). Only this numerical approach accounts correctly for dynamical friction and tidal effects between Fornax and the globular clusters. We show that a weak cusp ($r_{s}=$ 1.5 kpc) or a large core ($r_{c}=$ 848 pc) are not compatible with the current observed positions and masses of Fornax clusters. In contrast, a small dark matter core ($r_{c}=$ 282 pc) for Fornax naturally reproduces the cluster spatial and mass distributions over a wide range of initial globular cluster masses. We derive an upper limit of $r_{c}\lesssim$ 282 pc. This core size range favors only warm dark datter (WDM). It is also possible to obtain a compatible core size range from cold dark matter (CDM) theories, if the initial halo's central cusp is heated by gas to form a small core.

We determine the dipole in the WISE galaxy catalogue. After reducing star contamination to <0.1% by rejecting sources with high apparent motion and those close to the Galactic plane, we eliminate low redshift sources to suppress the non-kinematic, clustering dipole. We remove sources within {\pm}5{\deg} of the super-galactic plane, as well as those within 1'' of 2MRS sources at redshift z < 0.03. We enforce cuts on the source angular extent to preferentially select distant ones. As we progress along these steps, the dipole converges in direction to within 5{\deg} of the CMB dipole and its magnitude also progressively reduces but stabilises at {\sim}0.012, corresponding to a velocity >1000 km/s if it is solely of kinematic origin. However, previous studies have shown that only {\sim}70% of the velocity of the Local Group as inferred from the CMB dipole is due to sources at z < 0.03. We examine the Dark Sky simulations to quantify the prevalence of such environments and find that <2.1% of Milky Way-like observers in a {\Lambda}CDM universe should observe the bulk flow (> 240 km/s extending to z > 0.03) that we do. We construct mock catalogues in the neighbourhood of such peculiar observers in order to mimic our final galaxy selection and quantify the residual clustering dipole. After subtracting this the remaining dipole is 0.0048 {\pm} 0.0022, corresponding to a velocity of 420 {\pm} 213 km/s which is consistent with the CMB. However the sources (at z > 0.03) of such a large clustering dipole remain to be identified.

Previous studies have found our velocity in the rest frame of radio galaxies at high redshift to be substantially larger than that inferred from the CMB temperature dipole anisotropy. We construct a full sky catalogue NVSUMSS, by merging the NVSS and SUMSS catalogues and removing local sources by various means including cross-correlating with the 2MRS catalogue. We take into account both aberration and Doppler boost to deduce our velocity from the hemispherical number count asymmetry, as well as via a 3-dimensional linear estimator. Both the magnitude and direction depend on cuts made to the catalogue, e.g. on the lowest source flux, however these effects are small. With the hemispheric number count asymmetry method we obtain a velocity of 1729 $\pm$ 187 km/s i.e. about 4 times larger than that obtained from the CMB dipole, but close in direction, towards RA=149 $\pm$ 2 degree, DEC = -17 $\pm$ 12 degree. With the 3-dimensional estimator, the derived velocity is 1355 $\pm$ 174 km/s towards RA=141 $\pm$ 11 degree, DEC=-9 $\pm$ 10 degree. We assess the statistical significance of these results by constructing catalogues of random distributions and show that they are at best significant at the $2.81 \sigma$ (99.95% confidence) level.

The existence of critical points for the peculiar velocity field is a natural feature of the correlated vector field. These points appear at the junctions of velocity domains with different orientations of their averaged velocity vectors. Since peculiar velocities are the important cause of the scatter in the Hubble expansion rate, we propose that a more precise determination of the Hubble constant can be made by restricting analysis to a subsample of observational data containing only the zones around the critical points of the peculiar velocity field, associated with voids and saddle points. On large-scales the critical points, where the first derivative of the gravitational potential vanishes, can easily be identified using the density field and classified by the behavior of the Hessian of the gravitational potential. We use high-resolution N-body simulations to show that these regions are stable in time and hence are excellent tracers of the initial conditions. Furthermore, we show that the variance of the Hubble flow can be substantially minimized by restricting observations to the subsample of such regions of vanishing velocity instead of aiming at increasing the statistics by averaging indiscriminately using the full data sets, as is the common approach.

We propose a new probe to test the nature of gravity at various redshifts through large-scale cosmological observations. We use our void catalog, extracted from the Sloan Digital Sky Survey (SDSS, DR10), to trace the distribution of matter along the lines of sight to SNe Ia that are selected from the Union 2 catalog. We study the relation between SNe Ia luminosities and convergence and also the peculiar velocities of the sources. We show that the effects, on SNe Ia luminosities, of convergence and of peculiar velocities predicted by the theory of general relativity and theories of modified gravities are different and hence provide a new probe of gravity at various redshifts. We show that the present sparse large-scale data does not allow us to determine any statistically- significant deviation from the theory of general relativity but future more comprehensive surveys should provide us with means for such an exploration.

We propose that the accretion of a dwarf spheroidal galaxy provides a common origin for the giant southern stream and the warp of M31. We run about 40 full N-body simulations with live M31, infalling galaxies with varying masses and density profiles, and cosmologically-plausible initial orbital parameters. Excellent agreement with a full range of observational data is obtained for a model in which a dark-matter-rich dwarf spheroidal, whose trajectory lies on the thin plane of corotating satellites of M31, is accreted from its turnaround radius of about 200 kpcs into M31 at approximately 3 Gyrs ago. The satellite is disrupted as it orbits in the potential well of the galaxy and forms the giant stream and in return heats and warps the disk of M31. We show that our cosmologically-motivated model is favoured by the kinematic data over the phenomenological models in which the satellite starts its infall from a close distance of M31. Our model predicts that the remnant of the disrupted satellite resides in the region of the North-Eastern shelf of M31. The results here suggest that the surviving satellites of M31 that orbit on the same thin plane, as the disrupted satellite once did, could have all been accreted from an intergalactic filament.

We present catalogues of voids for the SDSS DR7 redshift survey and for Millennium I simulation mock data. We aim to compare the observations with simulations based on a $\Lambda$CDM model and a semi-analytic galaxy formation model. We use the void statistics as a test for these models. We assembled a mock catalogue that closely resembles the SDSS DR7 catalogue and carried out a parallel statistical analysis of the observed and simulated catalogue. We find that in the observation and the simulation, voids tend to be equally spherical. The total volume occupied by the voids and their total number are slightly larger in the simulation than in the observation. We find that large voids are less abundant in the simulation and the total luminosity of the galaxies contained in a void with a given radius is higher on average than observed by SDSS DR7 survey. We expect these discrepancies to be even more significant in reality than found here since the present value of $\sigma_8$ given by WMAP7 is lower than the value of 0.9 used in the Millennium I simulation. The reason why the simulation fails to produce enough large and dark voids might be the failure of certain semi-analytic galaxy formation models to reduce the small-scale power of $\Lambda$CDM and to produce sufficient power on large scales.

We generate the peculiar velocity field for the 2MASS Redshift Survey (2MRS) catalog using an orbit-reconstruction algorithm. The reconstructed velocities of individual objects in 2MRS are well-correlated with the peculiar velocities obtained from high-precision observed distances within 3,000 km/s. We estimate the mean matter density to be 0.31 +/- 0.05 by comparing observed to reconstructed velocities in this volume. The reconstructed motion of the Local Group in the rest frame established by distances within 3,000 km/s agrees with the observed motion and is generated by fluctuations within this volume, in agreement with observations. Then, we reconstruct the velocity field of 2MRS in successively larger radii, to study the problem of convergence towards the CMB dipole. We find that less than half of the amplitude of the CMB dipole is generated within a volume enclosing the Hydra-Centaurus-Norma supercluster at around 40 Mpc/h. Although most of the amplitude of the CMB dipole seems to be recovered by 120 Mpc/h, the direction does not agree and hence we observe no convergence up to this scale. We develop a statistical model which allows us to estimate cosmological para meters from the reconstructed growth of convergence of the velocity of the Local Group towards the CMB dipole motion. For scales up to 60 Mpc/h, assuming a Local Group velocity of 627 km/s, we estimate Omega_m h^2 = 0.11 +/- 0.06 and sigma_8=0.9 +/- 0.4, in agreement with WMAP5 measurements at the 1-sigma level. However, for scales up to 100 Mpc/h, we obtain Omega_m h^2 = 0.08 +/- 0.03 and sigma_8=1.0 +/- 0.4, which agrees at the 1 to 2-sigma level with WMAP5 results. (abridged)

Dark matter substructure around nearby galaxies provides an interesting opportunity for confusion-free indirect detection of dark matter. We calculate the boost over a smooth background distribution of dark matter for gamma-ray emission from dark matter self-annihilations in tidal structure in M31, assuming a cross-section inversely proportional to the relative velocities of the dark matter particles as proposed by the Sommerfeld effect. The low velocity of the material in the structure results in a significant increase in gamma-ray emission compared to both the background halo and the predicted emission for a velocity-independent cross section. We also calculate the expected signal for Fermi, for reasonable choices of the dark matter parameters. We find that for a cross section proportional to the inverse-square of the relative velocity, the enhancement to the annihilation rate is sufficient to test the velocity dependence of the cross section by spatial correlation with the stellar component of the stream, given sufficient detector sensitivity.

Dual AGNs are natural byproducts of hierarchical mergers of galaxies in the LambdaCDM cosmogony. Recent observations have shown that only a small fraction (~ 0.1%-1%) of AGNs at redshift z<~ 0.3 are dual with kpc-scale separations, which is rather low compared to the high merger rate of galaxies. Here we construct a phenomenological model to estimate the number density of dual AGNs and its evolution according to the observationally-estimated major merger rates of galaxies and various scaling relations on the properties of galaxies and their central massive black holes. We show that our model reproduces the observed frequency and separation distribution of dual AGNs provided that significant nuclear activities are triggered only in gas-rich progenitor galaxies with central massive black holes and only when the nuclei of these galaxies are roughly within the half-light radii of their companion galaxies. Under these constraints, the observed low dual AGN frequency is consistent with the relatively high merger rate of galaxies and supports the hypothesis that major mergers lead to AGN/QSO activities. We also predict that the number of kpc-scale dual AGNs decreases with increasing redshift and only about 0.02%--0.06% of AGNs are dual AGNs with double-peaked narrow line features at redshifts of z 0.5-1.2. Future observations of high-redshift dual AGNs would provide a solid test for this prediction.

The question of the transition to global isotropy from our anisotropic local Universe is studied using the Union 2 catalogue of Type Ia supernovae (SNe Ia). We construct a "residual" statistic sensitive to systematic shifts in their brightness in different directions and use this to search in different redshift bins for a preferred direction on the sky in which the SNe Ia are brighter or fainter relative to the 'standard' LCDM cosmology. At low redshift (z<0.05) we find that an isotropic model such as LCDM is barely consistent with the SNe Ia data at 2-3 sigma. A complementary maximum likelihood analysis of peculiar velocities confirms this finding -- there is a bulk flow of around 260 km/sec at z \sim 0.06, which disagrees with LCDM at 1-2 sigma. Since the Shapley concentration is believed to be largely responsible for this bulk flow, we make a detailed study of the infall region: the SNe Ia falling away from the Local Group towards Shapley are indeed significantly dimmer than those falling towards us and on to Shapley. Convergence to the CMB rest frame must occur well beyond Shapley (z>0.06) so the low redshift bulk flow can systematically bias any reconstruction of the expansion history of the Universe. At high redshifts z>0.15 the agreement between the SNe Ia data and the isotropic LCDM model does improve, however, the sparseness and low quality of the data means that LCDM cannot be singled out as the preferred cosmological model.

We carry out fully 3-dimensional simulations of evolution from self-similar, spherically symmetric linear perturbations of a Cold Dark Matter dominated Einstein-de Sitter universe. As a result of the radial orbit instability, the haloes which grow from such initial conditions are triaxial with major-to-minor axis ratios of order 3:1. They nevertheless grow approximately self-similarly in time. In all cases they have power-law density profiles and near-constant velocity anisotropy in their inner regions. Both the power-law index and the value of the velocity anisotropy depend on the similarity index of the initial conditions, the former as expected from simple scaling arguments. Halo structure is thus not "universal" but remembers the initial conditions. On larger scales the density and anisotropy profiles show two characteristic scales, corresponding to particles at first pericentre and at first apocentre after infall. They are well approximated by the NFW model only for one value of the similarity index. In contrast, at all radii within the outer caustic the pseudo phase-space density can be fit by a single power law with an index which depends only very weakly on the similarity index of the initial conditions. This behaviour is very similar to that found for haloes formed from LCDM initial conditions and so can be considered approximately universal.

We develop a theoretical framework for describing the hierarchical structure of the phase space of cold dark matter haloes, due to gravitationally bound substructures. Because it includes the full hierarchy of the cold dark matter initial conditions and is hence complementary to the halo model, the stable clustering hypothesis is applied for the first time here to the small-scale phase space structure. As an application, we show that the particle dark matter annihilation signal could be up to two orders of magnitude larger than that of the smooth halo within the Galactic virial radius. The local boost is inversely proportional to the smooth halo density, and thus is O(1) within the solar radius, which could translate into interesting signatures for dark matter direct detection experiments: The temporal correlation of dark matter detection can change by a factor of 2 in the span of 10 years, while there will be significant correlations in the velocity space of dark matter particles. This can introduce O(1) uncertainty in the direction of local dark matter wind, which was believed to be a benchmark of directional dark matter searches or the annual modulation signal.

We simulate the growth of isolated dark matter haloes from self-similar and spherically symmetric initial conditions. Our N-body code integrates the geodesic deviation equation in order to track the streams and caustics associated with individual simulation particles. The radial orbit instability causes our haloes to develop major-to-minor axis ratios approaching 10 to 1 in their inner regions. They grow similarly in time and have similar density profiles to the spherical similarity solution, but their detailed structure is very different. The higher dimensionality of the orbits causes their stream and caustic densities to drop much more rapidly than in the similarity solution. This results in a corresponding increase in the number of streams at each point. At 1% of the turnaround radius (corresponding roughly to the Sun's position in the Milky Way) we find of order 10^6 streams in our simulations, as compared to 10^2 in the similarity solution. The number of caustics in the inner halo increases by a factor of several, because a typical orbit has six turning points rather than one, but caustic densities drop by a much larger factor. This reduces the caustic contribution to the annihilation radiation. For the region between 1% and 50% of the turnaround radius, this is 4% of the total in our simulated haloes, as compared to 6.5% in the similarity solution. Caustics contribute much less at smaller radii. These numbers assume a 100 GeV c^-2 neutralino with present-day velocity dispersion 0.03 cm s^-1, but reducing the dispersion by ten orders of magnitude only doubles the caustic luminosity. We conclude that caustics will be unobservable in the inner parts of haloes. Only the outermost caustic might potentially be detectable.

Most of the mass content of dark matter haloes is expected to be in the form of tidal debris. The density of debris is not constant, but rather can grow due to formation of caustics at the apocenters and pericenters of the orbit, or decay as a result of phase mixing. In the phase space, the debris assemble in a hierarchy which is truncated by the primordial temperature of dark matter. Understanding this phase structure can be of significant importance for the interpretation of many astrophysical observations and in particular dark matter detection experiments. With this purpose in mind, we develop a general theoretical framework to describe the hierarchical structure of the phase space of cold dark matter haloes. We do not make any assumption of spherical symmetry and/or smooth and continuous accretion. Instead, working with correlation functions in the action-angle space, we can fully account for the hierarchical structure (predicting a two-point correlation function ~ (\Delta J)^{-1.6} in the action space), as well as the primordial discreteness of the phase space. As an application, we estimate the boost to the dark matter annihilation signal due to the structure of the phase space within virial radius: the boost due to the hierarchical tidal debris is of order unity, whereas the primordial discreteness of the phase structure can boost the total annihilation signal by up to an order of magnitude. The latter is dominated by the regions beyond 20% of the virial radius, and is largest for the recently formed haloes with the least degree of phase mixing.

Gravitational collapse of dark matter, merger of dark matter haloes and tidal disruption of satellites are among processes which lead to the formation of fine and dense dark matter shells, also known as dark matter caustics. The putative weakly interacting species which may form the dark matter are expected to strongly annihilate in these dense regions of the Milky Way halo and generate in particular antiprotons and positrons. We derive the flux of these rare antimatter particles at the Earth and show that it depends significantly on the cut-off radius of the dark matter distribution at the galactic centre. Boost factors of ~30 are found with respect to a smooth NFW profile for high-energy antiprotons and low-energy positrons if this cut-off radius is taken to be 300 pc -- a somewhat extreme value though. This yields a detectable antiproton signal around hundreds of Gev in models where the annihilation cross section today is enhanced by non--perturbative effects as in the generic case of a heavy Wino. However, dark matter caustics cannot provide a better explanation for the HEAT excess reported above ~10 GeV than a smooth NFW or isothermal cored distribution.

Anisotropic emission of gravitational waves during the merger or formation of black holes can lead to the ejection of these black holes from their host galaxies. A recoiled black hole which moves on an almost radial bound orbit outside the virial radius of its central galaxy, in the cold dark matter background, reaches its apapsis in a finite time. The low value of dark matter velocity dispersion at high redshifts and also the black hole velocity near the apapsis passage yield a high-density wake around these black holes. Gamma-ray emission can result from the enhancement of dark matter annihilation in these wakes. The diffuse high-energy gamma-ray background from the ensemble of such black holes in the Hubble volume is also evaluated.

A new scenario for the emission of high-energy gamma-rays from dark matter annihilation around massive black holes is presented. A black hole can leave its parent halo, by means of gravitational radiation recoil, in a merger event or in the asymmetric collapse of its progenitor star. A recoiled black hole which moves on an almost-radial orbit outside the virial radius of its central halo, in the cold dark matter background, reaches its apapsis in a finite time. Near or at the apapsis passage, a high-density wake extending over a large radius of influence, forms around the black hole. It is shown that significant gamma-ray emission can result from the enhancement of neutralino annihilation in these wakes. At its apapsis passage, a black hole is shown to produce a flash of high-energy gamma-rays whose duration is determined by the mass of the black hole and the redshift at which it is ejected. The ensemble of such black holes in the Hubble volume is shown to produce a diffuse high-energy gamma-ray background whose magnitude is compared to the diffuse emission from dark matter haloes alone.

Motion of a continuous fluid can be decomposed into an "incompressible" rearrangement, which preserves the volume of each infinitesimal fluid element, and a gradient map that transfers fluid elements in a way unaffected by any pressure or elasticity (the polar decomposition of Y. Brenier). The Euler equation describes a system whose kinematics is dominated by the incompressible rearrangement. The opposite limit, in which the incompressible component is negligible, corresponds to the Zel'dovich approximation, a model of motion of self-gravitating fluid in cosmology. We present a method of approximate reconstruction of the large-scale proper motions of matter in the Universe from the present-day mass density field. The method is based on recovering the corresponding gradient transfer map. We discuss its algorithmics, tests of the method against mock cosmological catalogues, and its application to observational data, which result in tight constraints on the mean mass density Omega_m and age of the Universe.

Lagrangian reconstruction of large-scale peculiar velocity fields can be strongly affected by observational biases. We develop a thorough analysis of these systematic effects by relying on specially selected mock catalogues. For the purpose of this paper, we use the MAK reconstruction method, although any other Lagrangian reconstruction method should be sensitive to the same problems. We extensively study the uncertainty in the mass-to-light assignment due to luminosity incompleteness, and the poorly-determined relation between mass and luminosity. The impact of redshift distortion corrections is analyzed in the context of MAK and we check the importance of edge and finite-volume effects on the reconstructed velocities. Using three mock catalogues with different average densities, we also study the effect of cosmic variance. In particular, one of them presents the same global features as found in observational catalogues that extend to 80 Mpc/h scales. We give recipes, checked using the aforementioned mock catalogues, to handle these particular observational effects, after having introduced them into the mock catalogues so as to quantitatively mimic the most densely sampled currently available galaxy catalogue of the nearby universe. Once biases have been taken care of, the typical resulting error in reconstructed velocities is typically about a quarter of the overall velocity dispersion, and without significant bias. We finally model our reconstruction errors to propose an improved Bayesian approach to measure Omega_m in an unbiased way by comparing the reconstructed velocities to the measured ones in distance space, even though they may be plagued by large errors. We show that, in the context of observational data, a nearly unbiased estimator of Omega_m may be built using MAK reconstruction.

The peculiar velocity reconstruction methods allow one to have a deeper insight into the distribution of dark matter: both to measure mean matter density and to obtain the primordial density fluctuations. We present here the Monge-Ampere-Kantorovitch method applied to mock catalogues mimicking in both redshift and distance catalogues. After having discussed the results obtained for a class of biases that may be corrected for, we focus on the systematics coming from the unknown distribution of unobserved mass and from the statistical relationship between mass and luminosity. We then show how to use these systematics to put constraints on the dark matter distribution. Finally a preliminary application to an extended version (c z < 3000 km/s) of the Neighbour Galaxy Catalogue is presented. We recover the peculiar velocities in our neighbourhood and present a preliminary measurement of the local Omega_M.

Cold dark matter haloes are populated by caustics, which are yet to be resolved in N-body simulations or observed in the Universe. Secondary infall model provides a paradigm for the study of caustics in "typical" haloes assuming that they have had no major mergers and have grown only by smooth accretion. This is a particular characteristic of the smallest dark matter haloes of about 10^{-5} Mo, which although "atypical" contain no substructures and could have survived until now with no major mergers. Thus using this model as the first guidline, we evaluate the neutralino self-annihilation flux for these haloes. Our results show that caustics could leave a distinct sawteeth signature on the differential and cumulative fluxes coming from the outer regions of these haloes. The total annihilation signal from the regions away from the centre can be boosted by about forty percents.

Caustics are high-density structures that form frequently in collisionless media. Under self-gravity, cold dark matter flows focus onto caustics which are yet to be resolved in numerical simulations and or observed in the real world. If detected, caustics would provide a strong evidence for dark matter and would rule out alternative models such as those with modified dynamics. Here, we demonstrate how they might be observed in the weak lensing data. We evaluate the shear distortion and show that its radial profile is marked by a characteristic sawtooth pattern due to the caustics in dark matter haloes that form by selfsimilar accretion. We discuss the observational complications, mainly due to the poor knowledge of the virial radii of the haloes and demonstrate that a superposition of about 300 cluster-size haloes would give a signal-to-noise ratio which is sufficiently large for the detection of caustics with ground-based observations. This number is reduced to 100 for space-based observations. These bounds can be easily achieved by the ongoing wide field optical surveys such as CFHTLS and the future space-based projects SNAP and DUNE which have to be accompanied by a X-ray follow-up of the selected clusters for a precise determination of their virial radii.

Peculiar velocities thoughout the region of the local supercluster are reconstructed by two different orbit-retracing methods. The requirement of the optimal correlation between the radial components of reconstructed velocities and the observed peculiar velocities derived from our extensive new catalog of distances puts stringent constraints on the values of the cosmological parameters. Our constraints intersect those from studies of microwave background fluctuations and statistical properties of galaxy clustering: the ensemble of constraints are consistent with Omega_m=0.22\pm 0.02. While motions throughout the Local Supercluster provide a measure of the mean ratio of mass to light, there can be large local fluctuations. Our reconstruction of the infall velocities in the immediate vicinity of the Virgo Cluster shows that there is a mass-to-light anomaly of a factor of 3 to 6 between groups in the general field environment and the heavily populated Virgo Cluster.

Caustics are formally singular structures, with infinite density, that form in collisionless media. The non-negligible velocity dispersion of dark matter particles renders their density finite. We evaluate the maximum density of the caustics within the framework of secondary infall model of formation of dark matter haloes. The result is then used to demonstrate that caustics can be probed by properly stacking the weak-lensing signal of about 600 haloes. CFHTLS accompanied by X-ray observations and the space-based experiments like SNAP or DUNE can provide us with the required statistics and hence a way of distinguishing between the viable dark matter particle candidates. The extension of our results to more realistic models including the effects of mergers of haloes is briefly outlined.

Cold dark matter haloes are populated by high-density structures with sharply-peaked profiles known as caustics which have not yet been resolved by 3-dimensional numerical simulations. Here, we derive semi-analytic expressions for the density profiles near caustics in haloes which form by self-similar accretions of dark matter with infinitesimal velocity dispersion. A simple rescaling shows that similarly to the case of absolutely cold medium these profiles are universal: they are valid for all caustics and irrespective of physical parameters of the halo. We derive the maximum density of the caustics and show that it depends on the velocity dispersion and the caustic location. We show that both the absolute and relative thickness of the caustic monotonically decrease toward the center of the halo while the maximum density grows. This indicates that the radial component of the thermal velocities decreases in the inner streams, i.e. the collisionless medium cools down in the radial direction descending to the center of the halo. Finally, we demonstrate that there can be a significant contribution to the emission measure from dark matter particle annihilation in the caustics.

The Monge-Ampere-Kantorovich (MAK) reconstruction is tested against cosmological N-body simulations. Using only the present mass distribution sampled with particles, and the assumption of homogeneity of the primordial distribution, MAK recovers for each particle the non-linear displacement field between its present position and its Lagrangian position on a primordial uniform grid. To test the method, we examine a standard LCDM N-body simulation with Gaussian initial conditions and 6 models with non-Gaussian initial conditions: a chi-squared model, a model with primordial voids and four weakly non-Gaussian models. Our extensive analyses of the Gaussian simulation show that the level of accuracy of the reconstruction of the nonlinear displacement field achieved by MAK is unprecedented, at scales as small as about 3 Mpc. In particular, it captures in a nontrivial way the nonlinear contribution from gravitational instability, well beyond the Zel'dovich approximation. This is also confirmed by our analyses of the non-Gaussian samples. Applying the spherical collapse model to the probability distribution function of the divergence of the displacement field, we also show that from a well-reconstructed displacement field, such as that given by MAK, it is possible to accurately disentangle dynamical contributions induced by gravitational clustering from possible initial non-Gaussianities, allowing one to efficiently test the non-Gaussian nature of the primordial fluctuations. In addition, a simple application of MAK using the Zel'dovich approximation allows one to also recover accurately the present-day peculiar velocity field on scales of about 8 Mpc.

A reconstruction method for recovering the initial conditions of the Universe starting from the present galaxy distribution is presented which guarantees uniqueness of solutions. We show how our method can be used to obtain the peculiar velocities of a large number of galaxies, hence trace galaxies orbits back in time and obtain the entire past dynamical history of the Universe above scales where multi-streaming has not occurred. When tested against a 128^3 LambdaCDM simulation in a box of 200 Mpc/h length, we obtain 60% exact reconstruction on scales above 6 Mpc/h. We apply our method to a real galaxy redshift catalogue, the updated NBG (Nearby Galaxies), containing 1483 galaxies, groups and clusters in a radius of 30 Mpc/h, and reconstruct the peculiar velocity fields in the local neighbourhood. Our reconstructed distances are well-matched to the observed values outside the collapsed regions if Omega_m(t)= 0.20 exp(-0.26 (t-13)) where t is the age of the Universe in Gyrs.

The observed infall of galaxies into the Virgo Cluster puts strong constraints on the mass of the cluster. A non-parametric fully non-linear description of the infall can be made with orbit reconstructions based on Numerical Action Methods. The mass of the cluster is determined to be $1.2 \times 10^{15} M_{\odot}$. The mass-to-light ratio for the cluster is found to be seven times higher than the mean ratio found across the region within V=3000 km/s.

Possible indirect detection of neutralino, through its gamma-ray annihilation product, by the forthcoming GLAST satellite from our galactic halo, M31, M87 and the dwarf galaxies Draco and Sagittarius is studied. Gamma-ray fluxes are evaluated for the two representative energy thresholds, 0.1 GeV and 1.0 GeV, at which the spatial resolution of GLAST varies considerably. Apart from dwarfs which are described either by a modified Plummer profile or by a tidally-truncated King profiles, fluxes are compared for halos with central cusps and cores. It is demonstrated that substructures, irrespective of their profiles, enhance the gamma-ray emission only marginally. The expected gamma-ray intensity above 1 GeV at high galactic latitudes is consistent with the residual emission derived from EGRET data if the density profile has a central core and the neutralino mass is less than 50 GeV, whereas for a central cusp only a substantial enhancement would explain the observations. From M31, the flux can be detected above 0.1 GeV and 1.0 GeV by GLAST only if the neutralino mass is below 300 GeV and if the density profile has a central cusp, case in which a significant boost in the gamma-ray emission is produced by the central black hole. For Sagittarius, the flux above 0.1 GeV is detectable by GLAST provided the neutralino mass is below 50 GeV. From M87 and Draco the fluxes are always below the sensitivity limit of GLAST.

We present new results on the angular momentum evolution of dark matter halos. Halos, from N-body simulations, are classified according to their mass growth histories into two categories: the accretion category contains halos whose mass has varied continuously and smoothly, while the merger category contains halos which have undergone sudden and significant mass variations (greater than 1/3 of their initial mass per event). We observe for the growth of angular momentum that it continues well into the nonlinear regime, that it is correlated with mass growth and that it occurs in two phases: a fast followed by a slow phase. The spin parameter decreases with time for the accretion catalog whereas slightly increases for the merger catalog, but does not change when two catalogs are considered together, and furthermore it does not depend on halo mass or cosmological model. We also develop a simple model for the formation of a disk galaxy similar to Milky Way, and conclude that our own halo must have captured satellites in order to acquire the required angular momentum and to achieve most of the disk around z ~ 1.6. Angular momentum distribution indicates that at z ~ 1.6 only 22% of the halos have angular momentum of magnitude comparable to that of disk galaxies in the mass range of $(0.1-5)*10^{11}$ of solar mass, clearly insufficient to explain the present observed abundance of these objects.

We show that the deterministic past history of the Universe can be uniquely reconstructed from the knowledge of the present mass density field, the latter being inferred from the 3D distribution of luminous matter, assumed to be tracing the distribution of dark matter up to a known bias. Reconstruction ceases to be unique below those scales -- a few Mpc -- where multi-streaming becomes significant. Above 6 Mpc/h we propose and implement an effective Monge-Ampere-Kantorovich method of unique reconstruction. At such scales the Zel'dovich approximation is well satisfied and reconstruction becomes an instance of optimal mass transportation, a problem which goes back to Monge (1781). After discretization into N point masses one obtains an assignment problem that can be handled by effective algorithms with not more than cubic time complexity in N and reasonable CPU time requirements. Testing against N-body cosmological simulations gives over 60% of exactly reconstructed points. We apply several interrelated tools from optimization theory that were not used in cosmological reconstruction before, such as the Monge-Ampere equation, its relation to the mass transportation problem, the Kantorovich duality and the auction algorithm for optimal assignment. Self-contained discussion of relevant notions and techniques is provided.

A method for the reconstruction of the primordial density fluctuation field is presented. Various previous approaches to this problem rendered {\it non-unique} solutions. Here, it is demonstrated that the initial positions of dark matter fluid elements, under the hypothesis that their displacement is the gradient of a convex potential, can be reconstructed uniquely. In our approach, the cosmological reconstruction problem is reformulated as an assignment problem in optimisation theory. When tested against numerical simulations, our scheme yields excellent reconstruction on scales larger than a few megaparsecs.

A review of the present observational and theoretical status of elliptical galaxies is presented with the aim to clarify whether the monolithic or the hierarchical, is a more viable scenario for the origin of these structures. We describe the dynamical structure of elliptical galaxies using photometric and spectroscopic data, in particular 3D observations from integral field spectrographs, with emphasis on properties such as brightness distribution, velocity profiles, central structures like "cuspy" profiles and "cores", as well as central supermassive black holes. We also report on the main relations between these quantities like the fundamental plane and colour-luminosity diagram. We present observational evidences for the presence of dark matter in the elliptical galaxies and examine the theories of galaxy formation within the framework of a cold dark matter cosmological model. We discuss the formation of large-scale structure, Press-Schechter theory and universal density profile of dark matter halos. Subsequently, gas dynamics, star formation, feedback, angular momentum, morphology and the epoch of galaxy formation are studied and comparison between disks and spheroids are made. Valuable insights into the formation epoch of elliptical galaxies are provided by deriving the mean metallicities and ages at different redshifts through the study of different population synthesis models (single-burst and evolutionary models) in particular the magnesium to iron ratio.

Reconstructing the density fluctuations in the early Universe that evolved into the distribution of galaxies we see today is a challenge of modern cosmology [ref.]. An accurate reconstruction would allow us to test cosmological models by simulating the evolution starting from the reconstructed state and comparing it to the observations. Several reconstruction techniques have been proposed [8 refs.], but they all suffer from lack of uniqueness because the velocities of galaxies are usually not known. Here we show that reconstruction can be reduced to a well-determined problem of optimisation, and present a specific algorithm that provides excellent agreement when tested against data from N-body simulations. By applying our algorithm to the new redshift surveys now under way [ref.], we will be able to recover reliably the properties of the primeval fluctuation field of the local Universe and to determine accurately the peculiar velocities (deviations from the Hubble expansion) and the true positions of many more galaxies than is feasible by any other method. A version of the paper with higher-quality figures is available at http://www.obs-nice.fr/etc7/nature.pdf

The minimal energy variations of a directed polymer with tilted columnar disorder in two dimensions are shown numerically to obey a multiscaling at short distances which crosses over to global simple scaling at large distances. The scenario is analogous to that of structure functions in bifractal Burgers' turbulence. Some scaling properties are predicted from extreme value statistics. The multiscaling disappears for zero tilt.

A {\it stochastic adhesion} model is introduced, with the purpose of describing the formation and evolution of mildly nonlinear structures, such as sheets and filaments, in the intergalactic medium (IGM), after hydrogen reionization. The model is based on replacing the overall force acting on the baryon fluid -- as it results from the composition of local gravity, pressure gradients and Hubble drag -- by a mock external force, self-consistently calculated from first-order perturbation theory. A small kinematic viscosity term prevents shell-crossing on small scales (which arises because of the approximate treatment of pressure gradients). The emerging scheme is an extension of the well-known adhesion approximation for the dark matter dynamics, from which it only differs by the presence of a small-scale `random' force, characterizing the IGM. Our algorithm is the ideal tool to obtain the skeleton of the IGM distribution, which is responsible for the structure observed in the low-column density Ly$\alpha$ forest in the absorption spectra of distant quasars.

This paper seeks to check the validity of the "apparent fractal conjecture" (Ribeiro 2001ab: gr-qc/9909093, astro-ph/0104181), which states that the observed power-law behaviour for the average density of large-scale distribution of galaxies arises when some observational quantities, selected by their relevance in average density profile determination, are calculated along the past light cone. Implementing these conditions in the proposed set of observational relations profoundly changes the behaviour of many observables in the standard cosmological models. In particular, the average density becomes observationally inhomogeneous, even in the spatially homogeneous spacetime of standard cosmology, change which was already analysed by Ribeiro (1992b, 1993, 1994, 1995: astro-ph/9910145) for a non-perturbed model. Here we derive observational relations in a perturbed Einstein-de Sitter cosmology by means of the perturbation scheme proposed by Abdalla and Mohayaee (1999: astro-ph/9810146), where the scale factor is expanded in power series to yield perturbative terms. The differential equations derived in this perturbative context, and other observables necessary in our analysis, are solved numerically. The results show that our perturbed Einstein-de Sitter cosmology can be approximately described by a decaying power-law like average density profile, meaning that the dust distribution of this cosmology has a scaling behaviour compatible with the power-law profile of the density-distance correlation observed in the galaxy catalogues. These results show that, in the context of this work, the apparent fractal conjecture is correct.

We use de Vaucouleurs' power-law density-distance relation, to study a hierarchical perturbation of the Friedmann universe. We solve the Einstein equation and obtain the density contrast and the amplification factor for the perturbation. It is shown that, scale-invariant inhomogeneities decay in Einstein-de-Sitter universe. On the contrary, in an open universe, the inhomogeneities grow. For low values of $\Omega$, amplification peaks sharply and the fluctuations can grow by up to a factor of $10^{13}$ from the recombination to the present time. Our analysis of the closed universe further confirms that, unlike the common belief that perturbations grow faster with increasing $\Omega$, scale-invariant perturbations amplify with decreasing $\Omega$. This result is consistent with the very low density expected for a hierarchical universe.

The evolution of inhomogeneities in a spherical collapse model is studied by expanding the Einstein equation in powers of inverse radial parameter. In the linear regime, the density contrast is obtained for flat, closed and open universes. In addition to the usual modes, an infinite number of new growing modes are contained in the solutions for pressureless open and closed universes. In the nonlinear regime, we obtain the leading growing modes in closed forms for a flat universe and also, in the limits of small and large times, for an open universe.

We derive the time evolution of the density contrast to all orders of perturbation theory, by solving the Einstein equation for scale-invariant fluctuations. These fluctuations are represented by an infinite series in inverse powers of the radial parameter. In addition to the standard growing modes, we find infinitely many more new growing modes for open and closed universes.

We write a non-relativistic Lagrangian for a hierarchical universe. The equations of motion are solved numerically and the evolution of the fractal dimension is obtained for different initial conditions. We show that our model is homogeneous at the time of the last scattering, but evolves into a self-similar universe with a remarkably constant fractal dimension. We also show that the Hubble law is implied by this model and make an estimate for the age of the universe.

It is shown that density fluctuations obey a scaling law in an open Friedmann universe. In a flat universe, the fluctuations are not scale-invariant. We compute the growth rate of adiabatic scale-invariant density fluctuations in flat, open and inflationary universes. We find that, given a sufficiently long time, the density perturbations decay away in the Einstein-de-Sitter universe. On the contrary, the rapid growth of the density instabilities makes an open universe inhomogeneous in a time scale comparable to the age of our universe. We also find that the fluctuations grow exponentially in a flat inflationary universe.