The epoch corresponding to a redshift of z $\sim 6.5$ is close to full re-ionisation of the Universe, and early enough to provide an intriguing environment to observe the early stage of large-scale structure formation. It is also en epoch that can be used to verify the abundance of a large population of low luminosity star-forming galaxies, that are deemed responsible for cosmic re-ionisation. Here, we present the results of follow-up multi-object spectroscopy using OSIRIS at Gran Telescopio Canarias (GTC) of 16 Ly$\alpha$ emitter (LAE) candidates discovered in the Subaru/XMM Newton Deep Survey. We have securely confirmed 10 LAEs with sufficient signal-to-noise ratio of the Ly$\alpha$ emission line. The inferred star formation rates of the confirmed LAEs are on the low side, within the range 0.9-4.7 M$_{\odot}$ yr$^{-1}$. However, they show relatively high Ly$\alpha$ rest frame equivalent widths. Finally we have shown that the mechanical energy released by the star formation episodes in these galaxies is enough to create holes in the neutral hydrogen medium such that Lyman continuum photons can escape to the intergalactic medium, thus contributing to the re-ionisation of the Universe.

We present the results of our search for the faint galaxies near the end of the Reionisation Epoch. This has been done using very deep OSIRIS images obtained at the Gran Telescopio Canarias (GTC). Our observations focus around two close, massive Lyman Alpha Emitters (LAEs) at redshift 6.5, discovered in the SXDS field within a large-scale overdense region (Ouchi et al. 2010). The total GTC observing time in three medium band filters (F883w35, F913w25 and F941w33) is over 34 hours covering $7.0\times8.5$ arcmin$^2$ (or $\sim30,000$ Mpc$^3$ at $z=6.5$). In addition to the two spectroscopically confirmed LAEs in the field, we have identified 45 other LAE candidates. The preliminary luminosity function derived from our observations, assuming a spectroscopic confirmation success rate of $\frac{2}{3}$ as in previous surveys, suggests this area is about 2 times denser than the general field galaxy population at $z=6.5$. If confirmed spectroscopically, our results will imply the discovery of one of the earliest protoclusters in the universe, which will evolve to resemble the most massive galaxy clusters today.

We combine observational data on a dozen independent cosmic properties at high-$z$ with the information on reionization drawn from the spectra of distant luminous sources and the cosmic microwave background (CMB) to constrain the interconnected evolution of galaxies and the intergalactic medium since the dark ages. The only acceptable solutions are concentrated in two narrow sets. In one of them reionization proceeds in two phases: a first one driven by Population III stars, completed at $z\sim 10$, and after a short recombination period a second one driven by normal galaxies, completed at $z\sim 6$. In the other set both kinds of sources work in parallel until full reionization at $z\sim 6$. The best solution with double reionization gives excellent fits to all the observed cosmic histories, but the CMB optical depth is 3-$\sigma$ larger than the recent estimate from the Planck data. Alternatively, the best solution with single reionization gives less good fits to the observed star formation rate density and cold gas mass density histories, but the CMB optical depth is consistent with that estimate. We make several predictions, testable with future observations, that should discriminate between the two reionization scenarios. As a byproduct our models provide a natural explanation to some characteristic features of the cosmic properties at high-$z$, as well as to the origin of globular clusters.

The modeling of galaxy formation and reionization, two central issues of modern cosmology, relies on the accurate follow-up of the intergalactic medium (IGM). Unfortunately, owing to the complex nature of this medium, the differential equations governing its ionization state and temperature are only approximate. In this paper, we improve these master equations. We derive new expression for the distinct composite inhomogeneous IGM phases, including all relevant ionizing/recombining and cooling/heating mechanisms, taking into account inflows/outflows into/from halos, and using more accurate recombination coefficients. Furthermore, to better compute the source functions in the equations we provide an analytic procedure for calculating the halo mass function in ionized environments, accounting for the bias due to the ionization state of their environment. Such an improved treatment of IGM evolution is part of a complete realistic model of galaxy formation presented elsewhere.

We present an Analytic Model of Intergalactic-medium and GAlaxy evolution since the dark ages. AMIGA is in the spirit of the popular semi-analytic models of galaxy formation, although it does not use halo merger trees but interpolates halo properties in grids that are progressively built. This strategy is less memory-demanding and allows one to start the modeling at redshifts high enough and halo masses low enough to have trivial boundary conditions. The number of free parameters is minimized by making the causal connection between physical processes usually treated as independent from each other, which leads to more reliable predictions. But the strongest points of AMIGA are: i) the inclusion of molecular cooling and metal-poor, population III (Pop III) stars, with the most dramatic feedback, and ii) the accurate follow-up of the temperature and volume filling factor of neutral, singly, and doubly ionized regions, taking into account the distinct halo mass functions in those environments. We find the following general results. Massive Pop III stars determine the IGM metallicity and temperature, and the growth of spheroids and disks is self-regulated by that of massive black holes developed from the remnants of those stars. Yet, the properties of normal galaxies and active galactic nuclei appear to be quite insensitive to Pop III star properties owing to the much higher yield of ordinary stars compared to Pop III stars and the dramatic growth of MBHs when normal galaxies begin to develop, which cause the memory loss of the initial conditions.

We establish a one-to-one correspondence between virialized haloes and their seeds, namely peaks with a given density contrast at appropriate Gaussian-filtering radii, in the initial Gaussian random density field. This fixes a rigorous formalism for the analytic derivation of halo properties from the linear power spectrum of density perturbations in any hierarchical cosmology. The typical spherically averaged density profile and mass function of haloes so obtained match those found in numerical simulations.

Comparing the excursion set and CUSP formalisms for the derivation of the halo mass function, we investigate the role of the mass definition in the properties of the multiplicity function of cold dark matter (CDM) haloes. We show that the density profile for haloes formed from triaxial peaks that undergo ellipsoidal collapse and virialisation is such that the ratio between the mean inner density and the outer local density is essentially independent of mass. This causes that, for suited values of the spherical overdensity $\Delta$ and the linking length $b$, SO and FoF masses are essentially equivalent to each other and the respective multiplicity functions are essentially the same. The overdensity for haloes having undergone ellipsoidal collapse is the same as if they had formed according to the spherical top-hat model, which leads to a value of $b$ corresponding to the usual virial overdensity, $\Delta_{vir}$, equal to $\sim$ 0.2. The multiplicity function resulting from such mass definitions, expressed as a function of the top-hat height for spherical collapse, is very approximately universal in all CDM cosmologies. The reason for this is that, for such mass definitions, the top-hat density contrast for ellipsoidal collapse and virialisation is close to a universal value, equal to $\sim$ 0.9 times the usual top-hat density contrast for spherical collapse.

In a recent paper, Salvador-Sol\'e et al. (2012) have derived the typical inner structure of dark matter haloes from that of peaks in the initial random Gaussian density field, determined by the power-spectrum of density perturbations characterising the hierarchical cosmology under consideration. In the present paper, we extend this formalism to the typical kinematics and triaxial shape of haloes. Specifically, we establish the link between such halo properties and the power-spectrum of density perturbations through the typical shape of peaks. The trends of the predicted typical halo shape, pseudo phase-space density and anisotropy profiles are in good agreement with the results of numerical simulations. Our model sheds light on the origin of the power-law-like pseudo phase-space density profile for virialised haloes.

We derive the density profile for collisionless dissipationless dark matter haloes in hierarchical cosmologies making use of the Secondary Infall (SI) model. The novelties are: i) we deal with triaxial virialised objects; ii) their seeds in the linear regime are peaks endowed with {\it unconvolved} spherically averaged density profiles according to the peak formalism; iii) the initial peculiar velocities are taken into account; and iv) accreting haloes are assumed to develop from the inside out, keeping the instantaneous inner system unaltered. The validity of this latter assumption is accurately checked by comparing analytical predictions on such a growth with the results of numerical simulation. We show that the spherically averaged density profile of virialised objects can be inferred with no need to specify their shape. The {\it typical} spherically averaged halo density profile is inferred, down to arbitrarily small radii, from the power-spectrum of density perturbations. The predicted profile in the $\Lambda$CDM cosmology is approximately described by an Einasto profile, meaning that it does not have a cusp but rather a core, where the inner slope slowly converges to zero. Down to one hundredth the total radius, the profile has the right NFW and Einasto forms, being close to the latter down to a radius of about four orders of magnitude less. The inner consistency of the model implies that the density profiles of haloes harbour no information on their past aggregation history. This would explain why major mergers do not alter the typical density profile of virialised objects formed by SI and do not invalidate the peak formalism based on such a formation.

Using the model for (bottom-up) hierarchical halo growth recently developed by Salvador-Sol\'e et al. (2012), we derive the typical spherically averaged density profile for haloes with several relevant masses in the concordant warm dark matter (\Lambda WDM) cosmology with non-thermal sterile neutrinos of two different masses. The predicted density profiles become flat at small radii, as expected from the effects of the spectrum cutoff. The core cannot be resolved, however, because the non-null particle velocity yields the fragmentation of minimum mass protohaloes in small nodes, which invalidates the model at the corresponding radii.

In two previous papers (Salvador-Sol\'e 2012a,b), it was shown that: i) the typical structural and kinematic properties of haloes in (bottom-up) hierarchical cosmologies endowed with random Gaussian density perturbations of dissipationless collisionless dark matter emerge naturally from the typical properties of peaks in the primordial density field and ii) halo statistics are well described by the peak formalism. In the present paper, we use these results to model halo substructure. Specifically, making use of the peak formalism and the fact that accreting haloes evolve from the inside-out, we derive the subhalo mass abundance and number density profile per infinitesimal mass for subhaloes of different masses, as a function of the subhalo maximum circular velocity or mass, before and after the tidal truncation of subhaloes by the host potential well. The subhalo properties obtained by assuming that subhaloes are mainly made of diffuse particles are in very good agreement with those found in current high-resolution N-body simulations. We also predict the subhalo properties in the opposite extreme case, likely better suited for the real universe in CDM cosmologies, that haloes are made of subhaloes within subhaloes at all scales.

This paper has been withdrawn owing a re-arrangement of two previously submitted papers. The new version of the theoretical work on the triaxial shape of dark matter haloes can be found at the ArXiv astro-ph list (CO) as article 1104.2905

We perform a simple test showing that accreting halos grow inside-out and have density profiles that are statistically indistinguishable from those of halos having undergone major mergers. The test also reveals that the ability for the NFW analytical expression to fit the density profile of halos spoils towards large masses at any redshift. On view of these results, we analyze the M-c relations on extreme mass and redshift domains predicted by usual toy models based on numerical data and one physical model based on the two assumptions above confirmed by the test.

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

In a series of recent papers, a new formalism has been developed that explains the inner structure of dark matter halos as collisionless, dissipationless systems assembled through mergers and accretion at the typical cosmological rate. Nearby ellipticals are also collisionless, dissipationless systems assembling their mass through mergers, but contrarily to the former structures they do not continuously accrete external matter because they are shielded by their host halos. Here we explore the idea that the infall of their own matter ejected within the halo on the occasion of a violent merger can play a role similar to external accretion in halos. The predicted stellar mass density profile fits the observed one, and the empirical total mass density profile is also recovered.

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

We consider the dynamical evolution of small-scale substructure in clusters within two extreme alternate scenarios for their possible origin: 1) the accretion of groups (or small clusters) on quasi-radial orbits, and 2) the merger of clusters of similar masses, followed by the decoupling of their dense cores. Using simple analytical arguments and checking with numerical computations, we show that objects are destroyed by the tidal field of the global cluster potential if their mean density is small compared to the mean cluster density within the radius of closest approach of the group or detached core. Accreted groups and small clusters are thus tidally disrupted in one cluster crossing. Since the cores of clusters are much denser than groups, they are considerably more robust to tides, but the least massive are destroyed or severely stripped by tides, while the others are brought to the cluster center by dynamical friction (and subsequently merge) in less than one orbit. The longest lived substructures are detached cores, roughly ten times less massive than the cluster, starting in near-circular orbits beyond $1 \, h^{-1} \, \rm Mpc$ from the cluster center.