We show herein that a proto-cluster of Ly$\alpha$ emitting galaxies, spectroscopically confirmed at redshift 6.5, produces a remarkable number of ionising continuum photons. We start from the Ly$\alpha$ fluxes measured in the spectra of the sources detected spectroscopically. From these fluxes we derive the ionising emissivity of continuum photons of the proto-cluster, which we compare with the ionising emissivity required to reionise the proto-cluster volume. We find that the sources in the proto-cluster are capable of ionising a large bubble, indeed larger than the volume occupied by the proto-cluster. For various calculations we have used the model AMIGA, in particular to derive the emissivity of the Lyman continuum photons required to maintain the observed volume ionised. Besides, we have assumed the ionising photons escape fraction given by AMIGA at this redshift.

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 ALMA Atacama Compact Array (ACA) 870 micron continuum maps of 28 infrared-bright SDSS quasars with Herschel/SPIRE detections at redshifts 2-4, the largest such sample ever observed with ALMA. The ACA detections are centred on the SDSS coordinates to within 1 arcsec for about 80 per cent of the sample. Larger offsets indicate that the far-infrared (FIR) emission detected by Herschel might come from a companion source. The majority of the objects (about 70 per cent) have unique ACA counterparts within the SPIRE beam down to 3-4 arcsec resolution. Only 30 per cent of the sample shows clear evidence for multiple sources with secondary counterparts contributing to the total 870 micron flux within the SPIRE beam to at least 25 per cent. We discuss the limitations of the data based on simulated pairs of point-like sources at the resolution of the ACA and present an extensive comparison of our findings with recent works on the multiplicities of sub-millimetre galaxies. We conclude that, despite the coarse resolution of the ACA, our data support the idea that, for a large fraction of FIR-bright quasars, the sub-mm emission comes from single sources. Our results suggest that, on average, optically bright quasars with strong FIR emission are not triggered by early-stage mergers but are, instead, together with their associated star formation rates, the outcome of either late-stage mergers or secular processes.

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 present charm (cosmic history agnostic reconstruction method), a novel inference algorithm that reconstructs the cosmic expansion history as encoded in the Hubble parameter $H(z)$ from SNe Ia data. The novelty of the approach lies in the usage of information field theory, a statistical field theory that is very well suited for the construction of optimal signal recovery algorithms. The charm algorithm infers non-parametrically $s(a)=\ln(\rho(a)/\rho_{\mathrm{crit}0})$, the density evolution which determines $H(z)$, without assuming an analytical form of $\rho(a)$ but only its smoothness with the scale factor $a=(1+z)^{-1}$. The inference problem of recovering the signal $s(a)$ from the data is formulated in a fully Bayesian way. In detail, we have rewritten the signal as the sum of a background cosmology and a perturbation. This allows us to determine the maximum a posteriory estimate of the signal by an iterative Wiener filter method. Applying charm to the Union2.1 supernova compilation, we have recovered a cosmic expansion history that is fully compatible with the standard $\Lambda$CDM cosmological expansion history with parameter values consistent with the results of the Planck mission.

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.

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.

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.

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.

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.

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 present an entirely analytic model for a preheated, polytropic intergalactic medium in hydrostatic equilibrium within a NFW dark halo potential in which the evolution of the halo structure between major merger events proceeds inside-out by accretion. This model is used to explain, within a standard $\Lambda$CDM cosmogony, the observed X-ray properties of nearby relaxed, non-cooling flow groups and clusters of galaxies. We find that our preferred solution to the equilibrium equations produces scaling relations in excellent agreement with observations, while simultaneously accounting for the typical structural characteristics of the distribution of the diffuse baryons. In the class of preheating models, ours stands out because it offers a unified description of the intrahalo medium for galaxy systems with total masses above $\sm 2\times 10^{13}$\msun, does not produce baryonic configurations with large isentropic cores, and reproduces faithfully the observed behavior of the gas entropy at large radii. All this is achieved with a moderate level of energy injection of about half a keV, which can be easily accommodated within the limits of the total energy released by the most commonly invoked feedback mechanisms, as well as with a polytropic index of 1.2, consistent with both many observational determinations and predictions from high-resolution gas-dynamical simulations of non-cooling flow clusters. More interestingly, our scheme offers a physical motivation for the adoption of this specific value of the polytropic index, as it is the one that best ensures the conservation after halo virialization of the balance between the total specific energies of the gas and dark matter components for the full range of masses investigated.

Relaxed dark-matter haloes are found to exhibit the same universal density profiles regardless of whether they form in hierarchical cosmologies or via spherical collapse. Likewise, the shape parameters of haloes formed hierarchically do not seem to depend on the epoch in which the last major merger took place. Both findings suggest that the density profile of haloes does not depend on their aggregation history. Yet, this possibility is apparently at odds with some correlations involving the scale radius r_s found in numerical simulations. Here we prove that the scale radius of relaxed, non-rotating, spherically symmetric haloes endowed with the universal density profile is determined exclusively by the current values of four independent, though correlated, quantities: mass, energy and their respective instantaneous accretion rates. Under this premise and taking into account the inside-out growth of haloes during the accretion phase between major mergers, we build a simple physical model for the evolution of r_s along the main branch of halo merger trees that reproduces all the empirical trends shown by this parameter in N-body simulations. This confirms the conclusion that the empirical correlations involving r_s do not actually imply the dependence of this parameter on the halo aggregation history. The present results give strong support to the explanation put forward in a recent paper by Manrique et al. (2003) for the origin of the halo universal density profile.

We calculate by means of the Press-Schechter formalism the density profile developed by dark-matter halos during accretion, i.e., the continuous aggregation of small clumps. We find that the shape of the predicted profile is similar to that shown by halos in high-resolution cosmological simulations. Furthermore, the mass-concentration relation is correctly reproduced at any redshift in all the hierarchical cosmologies analyzed, except for very large halo masses. The role of major mergers, which can cause the rearrangement of the halo structure through violent relaxation, is also investigated. We show that, as a result of the boundary conditions imposed by the matter continuously infalling into the halo during the violent relaxation process, the shape of the density profile emerging from major mergers is essentially identical to the shape the halo would have developed through pure accretion. This result explains why, according to high-resolution cosmological simulations, relaxed halos of a given mass have the same density profile regardless of whether they have had a recent merger or not, and why both spherical infall and hierarchical assembly lead to very similar density profiles. Finally, we demonstrate that the density profile of relaxed halos is not affected either by the capture of clumps of intermediate mass.

We construct a physically motivated analytical model for the quasar luminosity function, based on the joint star formation and feeding of massive black holes suggested by the observed correlation between the black hole mass and the stellar mass of the hosting spheroids. The parallel growth of massive black holes and host galaxies is assumed to be triggered by major mergers of haloes. The halo major merger rate is computed in the frame of the extended Press-Schechter model. The evolution of black holes on cosmological timescales is achieved by the integration of the governing set of differential equations, established from a few reasonable assumptions that account for the distinct (Eddington-limited or supply-limited) accretion regimes. Finally, the typical lightcurves of the reactivated quasars are obtained under the assumption that, in such accretion episodes, the fall of matter onto the black hole is achieved in a self-regulated stationary way. The predicted quasar luminosity function is compared to the luminosity functions of the 2dF QSO sample and other, higher redshift data. We find good agreement in all cases, except for z<1 where the basic assumption of our model is likely to break down.

The presence on the Virgo cluster outskirts of spiral galaxies with gas deficiencies as strong as those of the inner galaxies stripped by the intracluster medium has led us to explore the possibility that some of these peripheral objects are not newcomers. A dynamical model for the collapse and rebound of spherical shells under the point mass and radial flow approximations has been developed to account for the amplitude of the motions in the Virgo I cluster (VIC) region. According to our analysis, it is not unfeasible that galaxies far from the cluster, including those in a gas-deficient group well to its background, went through its core a few Gyr ago. The implications would be: (1) that the majority of the HI-deficient spirals in the VIC region might have been deprived of their neutral hydrogen by interactions with the hot intracluster medium; and (2) that objects spending a long time outside the cluster cores might keep the gas deficient status without altering their morphology.

High-resolution N-body simulations of hierarchical cosmologies have shown that the density and velocity dispersion profiles of dark-matter haloes display well-definite universal forms whose origin remains unknown. In the present paper, we calculate the internal structure of haloes expected to arise in any such cosmologies by simply taking into account that halo growth proceeds through an alternate sequence of discrete major mergers and long periods of gentle accretion. Major mergers cause the violent relaxation of the system subject to the boundary conditions imposed by accreting layers beginning to fall in at that moment. Accretion makes the system develop inside-out from the previous seed according to the spherical infall model. The predicted structure is in very good agreement with the results of numerical simulations, particularly for moderate and low mass haloes. We find strong indications that the slight departure observed in more massive systems is not due to a poorer theoretical prediction, but to the more marked effects of the limited resolution used in the simulations on the empirical profiles. This may have important consequences on the reported universality of halo structure.

The growth of Super Massive Black Holes and the parallel development of activity in galactic nuclei are implemented in an analytic code of hierarchical clustering. The evolution of the luminosity function of quasars and AGN will be computed with special attention paid to the connection between quasars and Seyfert galaxies. One of the major interests of the model is the parallel study of quasar formation and evolution and the History of Star Formation.

We derive the atomic hydrogen content for a total of 1900 spirals in the fields of eighteen nearby clusters. By comparing the HI deficiency distributions of the galaxies inside of and outside of one Abell radius of each cluster, we find that two thirds of the clusters in our sample show a dearth of neutral gas in their interiors. Possible connections between the gaseous deficiency and the characteristics of both the underlying galaxies and their environment are investigated in order to gain insight into the mechanisms responsible for HI depletion. While we do not find a statistically significant variation of the fraction of HI-deficient spirals in a cluster with its global properties, a number of correlations emerge that argue in favor of the interplay between spiral disks and their environment. In the clusters in which neutral gas deficiency is pronounced, we see clear indications that the degree of HI depletion is related to the morphology of the galaxies and not to their optical size: early-type and, probably, dwarf spirals are more easily emptied of gas than the intermediate Sbc--Sc types. The wealth of 21-cm data collected for the Virgo region has made it possible to study the 2D pattern of HI deficiency in that cluster. We also find evidence that gas-poor spirals in HI-deficient clusters move on orbits more radial than those of the gas-rich objects. The implications of all these results on models of how galaxies interact with their environment are reviewed. Hydrodynamic effects appear as the most plausible cause of HI removal. (abridged)

High resolution cosmological N-body simulations show that the density profiles of dark matter halos in hierarchical cosmogonies are universal, with low mass halos typically denser than more massive ones. This mass-density correlation is interpreted as reflecting the earlier formation of less massive objects. We investigate this hypothesis in the light of formation times defined as the epoch at which halos experience their last major merger. We find that the characteristic density and the scale radius of halos are essentially proportional, respectively, to the critical density of the universe and the virial radius at the time of their formation. These two relations are consistent with the following simple evolutionary picture. Violent relaxation caused by major mergers rearrange the structure of halos leading to a universal dimensionless density profile. Between major mergers, halos gradually grow through the accretion of surrounding layers by keeping the central part steady and only expanding their virial radius as the critical density of the universe diminishes.

High-resolution N-body simulations of hierarchical clustering in a wide variety of cosmogonies show that the density profiles of dark matter halos are universal, with low mass halos being denser than their more massive counterparts. This mass-density correlation is interpreted as reflecting the earlier typical formation time of less massive objects. We investigate this hypothesis in the light of formation times defined as the epoch at which halos experience their last major merger. Such halo formation times are calculated by means of a modification of the extended Press & Schechter formalism which includes a phenomenological frontier, Delta_m, between tiny and notable relative mass captures leading to the distinction between merger and accretion. For Delta_m=0.6, we confirm that the characteristic density of halos is essentially proportional to the mean density of the universe at their time of formation. Yet, proportionality with respect to the critical density yields slightly better results for open universes. In addition, we find that the scale radius of halos is also essentially proportional to their virial radius at the time of formation. We show that these two relations are consistent with the following simple scenario. Violent relaxation caused by mergers rearranges the structure of halos leading to the same density profile with universal values of the dimensionless characteristic density and scale radius. Between mergers, halos grow gradually through the accretion of surrounding layers by keeping their central parts steady and expanding their virial radius as the critical density of the universe diminishes.

In two previous papers a semi-analytical model was presented for the hierarchical clustering of halos via gravitational instability from peaks in a random Gaussian field of density fluctuations. This model is better founded than the extended Press-Schechter model, which is known to agree with numerical simulations and to make similar predictions. The specific merger rate, however, shows a significant departure at intermediate captured masses. The origin of this was suspected as being the rather crude approximation used for the density of nested peaks. Here, we seek to verify this suspicion by implementing a more accurate expression for the latter quantity which accounts for the correlation among peaks. We confirm that the inclusion of the peak-peak correlation improves the specific merger rate, while the good behavior of the remaining quantities is preserved.

With the present paper we conclude the presentation of a semianalytical model of hierarchical clustering of bound virialized objects formed by gravitational instability from a random Gaussian field of density fluctuations. In paper I, we introduced the basic tool, the so-called ``confluent system'' formalism, able to follow the evolution of bound virialized objects in the peak model. This was applied to derive the mass function of objects. In the present paper, we calculate other important quantities characterizing the growth history of objects. This model is compared with a similar one obtained by Lacey \& Cole (1993) following the Press \& Schechter (1974) approach. The interest of the new modeling presented here is twofold: 1) it is formally better justified as far as peaks are more reasonable seeds of objects than the indetermined regions in the Press \& Schechter prescription, and 2) it distinguishes between merger and accretion enabling us to unambiguously define the formation and destruction of objects and to estimate growth rates and characteristic times not available in the previous approach.

This is the first paper of a series of two devoted to develop a practical method to describe the growth history of bound virialized objects in the gravitational instability scenario without resorting to $N$-body simulations. Here we present the basic tool of this method, ``the confluent system formalism'', which allows us to follow the filtering evolution of peaks in a random Gaussian field of density fluctuations. This is applied to derive the theoretical mass function of objects within the peak model framework. Along the process followed for the derivation of this function, we prove that the Gaussian window is the only one consistent with the peak model ansatz. We also give a well justified derivation of the density of peaks with density contrast upcrossing a given threshold in infinitesimal ranges of scale and correct this scale function for the cloud-in-cloud effect. Finally, we characterize the form of the mass vs. scale and the critical overdensity vs. collapse time relations which are physically consistent with the peak model in an Einstein-de Sitter universe with density field endowed with different power spectra. The result is a fully justified semianalytical mass function which is close to the Press \& Schechter (1974) one giving good fits to $N$-body simulations. But the interest of the confluent system formalism is not merely formal. It allows us to distinguish between accretion and merger events, which is essential for the detailed modelling of the clustering process experienced by objects.