In order to attain a statistical description of the evolution of cosmic density fluctuations in agreement with results from the numerical simulations, we introduce a probability conditional formalism (CF) based on an inventory of isolated overdense regions in a density random field. This formalism is a useful tool for describing at the same time the mass function (MF) of dark haloes, their mass aggregation histories (MAHs) and merging rates (MRs). The CF focuses on virialized regions in a self-consistent way rather than in mass elements, and it offers an economical description for a variety of random fields. Within the framework of the CF, we confirm that, for a Gaussian field, it is not possible to reproduce at the same time the MF, MAH, and MR of haloes, both for a constant and moving barrier. Then, we develop an inductive method for constraining the cumulative conditional probability from a given halo MF description, and thus, using the CF, we calculate the halo MAHs and MRs. By applying this method to the MF measured in numerical simulations by Tinker et al. 2008, we find that a reasonable solution, justified by a mass conservation argument, is obtained if ones introduce a rescaling -increment by ~30% - of the virial mass used in simulations and a (slight) deviation from Gaussianity. Thus, both the MAH and MR obtained by a Monte Carlo merger tree agree now with the predictions of numerical simulations. We discuss on the necessity of rescaling the virial mass in simulations when comparing with analytical approaches on the ground of the matter not accounted as part of the halos and the halo mass limit due to numerical. Our analysis supports the presence of a diffuse dark matter component that is not taken into account in the measured halo MFs inasmuch as it is not part of the collapsed structures.

We use the Millennium I and II cosmological simulations to revisit the impact of mergers in the growth of bulges in central galaxies in the LCDM scenario. We seed galaxies within the growing CDM haloes using semi-empirical relations to assign stellar and gaseous masses, and an analytic treatment to estimate the transfer of stellar mass to the bulge of the remnant after a galaxy merger. We find that this model roughly reproduces the observed correlation between the bulge-to-total (B/T) mass ratio and stellar mass in present-day central galaxies as well as their observed demographics, although low-mass B/T<0.1 (bulgeless) galaxies might be scarce relative to the observed abundance. In our merger-driven scenario, bulges have a composite population made of (i) stars acquired from infalling satellites, (ii) stars transferred from the primary disc due to merger-induced perturbations, and (iii) newly formed stars in starbursts triggered by mergers. We find that (i) and (ii) are the main channels of mass assembly, with the first being dominant for massive galaxies, creating large bulges with different stellar populations than those of the inner discs, while the second is dominant for intermediate/low-mass galaxies creating small bulges with similar stellar populations to the inner discs. We associate the dominion of the first (second) channel to classical (pseudo) bulges, and compare the predicted fractions to observations. We remark that our treatment does not include other mechanisms of bulge growth such as intrinsic secular disc instabilities or misaligned gas accretion. We find that the evolution of the stellar and gaseous contents of the satellite as it moves towards the central galaxy is a key ingredient in setting the morphology of the remnant, and that a good match to the observed bulge demographics occurs when this evolution proceeds closely to that of the central galaxy.

An unified picture of stellar and halo mass build-up as a function of mass is presented. Inferred stellar-dark halo mass relations of galaxies, Ms-Mh, out to z=4 together with average LCDM halo mass aggregation histories (MAHs) are used for inferring average Ms growth histories, the Galaxian Hybrid Evolutionary Tracks (GHETs). The more massive the galaxy, the earlier transited in average from an active regime of Ms growth to a passive one: log(Mtran/Msun)=10.30+0.55z ("population downsizing"), where Mtran is the typical transition stellar mass. This result agrees with independent observational determinations based on the evolution of the galaxy stellar mass function decomposition into blue and red galaxies. The specific star formation rate, SSFR, predicted from the derivative of the GHET is consistent with direct measures of the SSFR for galaxies at different z's. The average GHETs of galaxies smaller than Mtran at z=0 (Ms~10^10.3 Msun) did not reach the quiescent regime, and for them, the lower the mass, the faster the later Ms growth rate ("downsizing in SSFR"). The GHETs allow to predict the transition rate in number density of active to passive population; the predicted values agree with direct estimates of growth rate in number density for the (massive) red population up to z~1. We show that LCDM-based models of disk galaxy evolution are able to reproduce the low-mass side of the Ms-Mh relation at z~0, but at higher z's disagree strongly with the GHETs: models do not reproduce the downsizing in SSFR and the high SSFR of low mass galaxies. (Abridged)

By means of the abundance matching technique, we infer the local stellar and baryonic mass-halo mass (Ms-Mh and Mb-Mh) relation for central blue and red galaxies separately in the mass range Ms~10^8.5-10^12.0 Msun. The observational inputs are the SDSS central blue and red Galaxy Stellar Mass Functions reported in Yang et al. 2009, and the measured local gas mass-Ms relations for blue and red galaxies. For the Halo Mass Function associated to central blue galaxies, the distinct LCDM one is used and set up to exclude: (i) the observed group/cluster mass function (blue galaxies are rare as centers of groups/clusters), and (ii) halos with a central major merger at resdshifts z<0.8 (dry late major mergers destroy the disks of blue galaxies). For red galaxies, we take the complement of this function to the total. The obtained mean Ms-Mh and Mb-Mh relations of central blue and red galaxies do not differ significantly from the respective relations for all central galaxies. For Mh>10^11.5 Msun, the Ms's of red galaxies tend to be higher than those of blue ones for a given Mh, the difference not being larger than 1.7. For Mh<10^11.5 Msun, this trend is inverted. For blue (red) galaxies: (a) the maximum value of fs=Ms/Mh is 0.021^{+0.016}_{-0.009} (0.034{+0.026}_{-0.015}) and it is attained atlog(Mh/Msun)~12.0 (log(Mh/Msun)~11.9); (b) fs\propto Mh (fs\propto Mh^3) at the low-mass end while at the high-mass end, fs\propto Mh^-0.4 (fs\propto Mh^-0.6). The baryon mass fractions, fb=Mb/Mh, of blue and red galaxies reach maximum values of fb=0.028^{+0.018}_{-0.011} and fb=0.034^{+0.025}_{-0.014}, respectively. For Mh<10^11.3 Msun, a much steeper dependence of fb on Mh is obtained for the red galaxies than for the blue ones. We discuss on the differences found in the fs-Mh and fb-Mh relations between blue and red galaxies in the light of of semi-empirical galaxy models.

(Abridged) By means of high-resolution cosmological simulations in the context of the LCDM scenario, the specific star formation rate (SSFR=SFR/Ms, Ms is the stellar mass)--Ms and stellar mass fraction (Fs=Ms/Mh, Mh is the halo mass)--Ms relations of low-mass galaxies (2.5< Mh/10^10 Msun <50 at redshift z=0) at different epochs are predicted. The Hydrodynamics ART code was used and some variations of the sub-grid parameters were explored. Most of simulated galaxies, specially those with the highest resolutions, have significant disk components and their structural and dynamical properties are in reasonable agreement with observations of sub-M* field galaxies. However, the SSFRs are 5-10 times smaller than the averages of several (compiled and homogenized here) observational determinations for field blue/star-forming galaxies at z<0.3 (at low masses, most of observed field galaxies are actually blue/star-forming). This inconsistency seems to remain even at z~1.5 though less drastic. The Fs of simulated galaxies increases with Mh as semi-empirical inferences show, but in absolute values the former are ~5-10 times larger than the latter at z=0; this difference increases probably to larger factors at z~1-1.5. The inconsistencies reported here imply that simulated low-mass galaxies (0.2<Ms/10^9 Msun <30 at z=0) assembled their stellar masses much earlier than observations suggest. This confirms the predictions previously found by means of LCDM-based models of disk galaxy formation and evolution for isolated low-mass galaxies (Firmani & Avila-Reese 2010), and highlight that our implementation of astrophysics into simulations and models are still lacking vital ingredients.

The emerging empirical picture of galaxy stellar mass (Ms) assembly shows that galaxy population buildup proceeds from top to down in Ms. By connecting galaxies to LCDM halos and their histories, individual (average) Ms growth tracks can be inferred. These tracks show that massive galaxies assembled their Ms the earlier the more massive the halo, and that less massive galaxies are yet actively growing in Ms, the more active the less massive is the halo. The predicted star formation rates as a function of mass and the downsizing of the typical mass that separate active galaxies from the passive ones agree with direct observational determinations. This implies that the LCDM scenario is consistent with these observations. The challenge is now to understand the baryonic physics that drives the significant and systematical shift of the stellar mass assembly of galaxies from the mass assembly of their corresponding halos (from halo upsizing to galaxy downsizing).

The observations show that less massive the galaxies are, the higher on average is their specific star formation rate (SSFR = SFR/Ms, Ms is the stellar mass). Such a trend, called the 'SSFR downsizing' (SSFR-DS) phenomenon, is seen for local and high-z (back to z~1-2) galaxy samples. We use observational data related only to disc galaxies and explore the average SSFR change with z for different masses. For Ms in the range ~10^9.5-10^10.5 Msun, the SSFR increases with (1+z) to a power that barely depends on Ms, and at all z's smaller galaxies have ever higher SSFRs. The latter strongly disagree with the LCDM hierarchical mass accretion rates. By means of self-consistent models of disc galaxy evolution inside growing LCDM halos, the effects that disc feedback-driven outflows and gas re-accretion have on the galaxy SSFR histories are explored. The parameters of the outflow and re-accretion schemes are tuned to reproduce the z~0 Mh-Ms relation inferred from observations. In the case of outflows only, the SSFR of individual model galaxies is roughly proportional to (1+z)^2.2 for all the masses with a normalization factor that depends on mass as Ms^0.1, i.e more massive galaxies have slightly larger SSFRs, contrary to the observed strong SSFR-DS trend. For the re-accretion cases, the dependence on z remains approximately the same as without re-infall, but the correlation on Ms even increases for most of the reasonable values of the model parameters. The comparison of models and observations in the SSFR-Ms plane at z~0 shows the divergent trend in SSFR as lower are the masses (upsizing vs downsizing). We explain why the models show the reported trends, and conclude that the SSFR-DS for low-mass galaxies poses a challenge for LCDM-based disc galaxy evolution models. (Abridged)

By means of galaxy evolutionary models, we explore the direct consequences of the LCDM cosmogony on the size evolution of galactic discs, avoiding intentionally the introduction of intermediate (uncertain) astrophysical processes. Based on the shape of the rotation curves and guided by a simplicity criterion, we adopt an average galaxy mass baryon fraction of 0.03. In order to study general behaviors, only models with the average initial conditions are analyzed. The stellar and B-band effective radii, R* and RB, of individual galaxies grow significantly with time (inside-out disc formation) with laws that are weakly dependent on mass, M*,or luminosity, LB. However, the change of R* with z at fixed M* is slow; for z<2.5, R*(M*=const) ~ (1+z)^-0.4. On the other hand, the change of RB with z at a fixed LB is strong and resembles the RB decreasing law of the individual models; roughly RB(LB=const) ~ (1+z)^-0.85 for z<0.75, and ~(1+z)^-1.1 for z>0.75. We find also that at z=0, R* ~ M*^0.38 and RB ~ LB^0.40, remaining the slopes of these relations almost the same up to z ~ 3. Our predictions are in reasonable agreement with observational inferences on the typical radius change with z of late-type galaxies more luminous (massive) than high values imposed by the selection effects. The models seem also to be consistent, within the large scatter, with the RB and LB values obtained from non complete samples of sub-L* late-type galaxies with available rest-frame photometric information at different z's. The properties and evolution of the LCDM haloes seem to be the main drivers of galaxy disc size evolution. Nevertheless, the models reveal a potential difficulty in explaining the observed steepening of the RB-LB relation with respect to the R*-M* one, an effect related to the well established color-magnitude relation.

We compare the spectral properties of 79 short and 79 long Gamma-Ray Bursts (GRBs) detected by BATSE and selected with the same limiting peak flux. Short GRBs have a low-energy spectral component harder and a peak energy slightly higher than long GRBs, but no difference is found when comparing short GRB spectra with those of the first 1-2 sec emission of long GRBs. These results confirm earlier findings for brighter GRBs. The bolometric peak flux of short GRBs correlates with their peak energy in a similar way to long bursts. Short and long GRBs populate different regions of the bolometric fluence-peak energy plane, short bursts being less energetic by a factor similar to the ratio of their durations. If short and long GRBs had similar redshift distributions, they would have similar luminosities yet different energies, which correlate with the peak energy E_peak for the population of long GRBs. We also test whether short GRBs are consistent with the E_peak-E_iso and E_peak-L_iso correlations for the available sample of short (6 events) and long (92 events) GRBs with measured redshifts and E_peak,obs: while short GRBs are inconsistent with the E_peak-E_iso correlation of long GRBs, they could follow the E_peak-L_iso correlation of long bursts. All the above indications point to short GRBs being similar to the first phases of long bursts. This suggests that a similar central engine (except for its duration) operates in GRBs of different durations.

For a sample of long GRBs with known redshift, we study the distribution of the evolutionary tracks on the rest-frame luminosity-peak energy Liso-Ep' diagram. We are interested in exploring the extension of the `Yonetoku' correlation to any phase of the prompt light curve, and in verifying how the high-signal prompt duration time, Tf, in the rest frame correlates with the residuals of such correlation (Firmani et al. 2006). For our purpose, we analyse separately two samples of time-resolved spectra corresponding to 32 GRBs with peak fluxes >1.8 phot cm^-2 s^-1 from the Swift-BAT detector, and 7 bright GRBs from the CGRO-BATSE detector previously processed by Kaneko et al. (2006). After constructing the Liso-Ep' diagram, we discuss the relevance of selection effects, finding that they could affect significantly the correlation. However, we find that these effects are much less significant in the Liso x Tf-Ep' diagram, where the intrinsic scatter reduces significantly. We apply further corrections for reducing the intrinsic scatter even more. For the sub-samples of GRBs (7 from Swift and 5 from CGRO) with measured jet break time, we analyse the effects of correcting Liso by jet collimation. We find that (i) the scatter around the correlation is reduced, and (ii) this scatter is dominated by the internal scatter of the individual evolutionary tracks. These results suggest that the time, integrated `Amati' and `Ghirlanda' correlations are consequences of the time resolved features, not of selection effects, and therefore call for a physical origin. We finally remark the relevance of looking inside the nature of the evolutionary tracks.

We explore how the slopes and scatters of the scaling relations of disk galaxies (Vm-L[-M], R-L[-M], and Vm-R) do change when moving from B to K bands and to stellar and baryonic quantities. For our compiled sample of 76 normal, non-interacting high and low surface brightness galaxies, we find some changes, which evidence evolution effects, mainly related to gas infall and star formation (SF). We also explore correlations among the (B-K) color, stellar mass fraction fs, mass M (luminosity L), and surface density (SB), as well as correlations among the residuals of the scaling relations. Some of our findings are: (i) the scale length Rb is a third parameter in the baryonic TF relation and the residuals of this relation follow a trend (slope ~-0.15) with the residuals of the Rb-Mb relation; for the stellar and K band cases, R is not anymore a third parameter and the mentioned trend disappears; (ii) among the TFRs, the B-band TFR is the most scattered; in this case, the color is a third parameter; (iii) the LSB galaxies break some observed trends, which suggest a threshold in the gas surface density Sg, below which the SF becomes independent of the gas infall rate and Sg. Our results are interpreted and discussed in the light of LCDM-based models of galaxy evolution. The models explain not only the baryonic scaling relations, but also most of the processes responsible for the observed changes in the slopes, scatters, and correlations among the residuals when changing to stellar and luminous quantities. The baryon fraction is required to be smaller than 0.05 on average. We detect some potential difficulties for the models: the observed color-M and surface density-M correlations are steeper, and the intrinsic scatter in the baryonic TFR is smaller than those predicted. [abridged]

The spectral-energy and (luminosity) correlations in long GRBs are being hotly debated to establish, first of all, their reality against possible selection effects. These are best studied in the observer planes, namely the peak energy E_peak_obs vs the fluence F or the peak flux P. In a recent paper we started to attack this problem considering all GRBs with known z and spectral properties. Here we consider instead all bursts with known E_peak_obs, irrespective of z, adding to those a sample of 100 faint BATSE bursts representative of a larger population. This allows us to construct a complete, fluence limited, sample, to study the selection/instrumental effects. We found that fainter bursts have smaller E_peak_obs than those of bright events. As a consequence, the E_peak_obs of these bursts is correlated with the fluence, though with a slope flatter than that defined by bursts with z. Selection effects, which are present, are shown not to be responsible for the existence of such a correlation. About 6% of these bursts are surely outliers of the E_peak-E_iso correlation (updated to include 83 bursts), since they are inconsistent with it for any z. E_peak_obs correlates also with P, with a slope similar to the E_peak-L_iso correlation.In this case there is only one sure outlier.The scatter of the E_peak_obs-P correlation defined by the BATSE bursts of our sample is smaller than the E_peak_obs-F correlation of the same bursts, while for the bursts with known z the E_peak-E_iso correlation is tighter than the E_peak-L_iso one. Once a very large number of bursts with E_peak_obs and z will be available, we thus expect that the E_peak-L_iso correlation will be similar to that currently found, whereas it is likely that the E_peak-E_iso correlation will become flatter and with a larger scatter.

We study the distribution of long Gamma Ray Bursts in the Ep-Eiso and in the Ep,obs-Fluence planes through an updated sample of 76 bursts, with measured redshift and spectral parameters, detected up to September 2007. We confirm the existence of a strong rest frame correlation Ep ~ Eiso^0.54+-0.01. Contrary to previous studies, no sign of evolution with redshift of the Ep-Eiso correlation (either its slope and normalisation) is found. The 76 bursts define a strong Ep,obs-Fluence correlation in the observer frame (Ep,obs ~ F^0.32+-0.05) with redshifts evenly distributed along this correlation. We study possible instrumental selection effects in the observer frame Ep,obs-Fluence plane. In particular, we concentrate on the minimum peak flux necessary to trigger a given GRB detector (trigger threshold) and the minimum fluence a burst must have to determine the value of Ep,obs (spectral analysis threshold). We find that the latter dominates in the Ep,obs-Fluence plane over the former. Our analysis shows, however, that these instrumental selection effects do not dominate for bursts detected before the launch of the Swift satellite, while the spectral analysis threshold is the dominant truncation effect of the Swift GRB sample (27 out of 76 events). This suggests that the Ep,obs-Fluence correlation defined by the pre--Swift sample could be affected by other, still not understood, selection effects. Besides we caution about the conclusions on the existence of the Ep,obs-Fluence correlation based on our Swift sample alone.

Long Gamma-Ray Bursts (GRBs) are the brightest electromagnetic explosions in the Universe, associated to the death of massive stars. As such, GRBs are potential tracers of the evolution of the cosmic massive star formation, metallicity, and Initial Mass Function. GRBs also proved to be appealing cosmological distance indicators. This opens a unique opportunity to constrain the cosmic expansion history up to redshifts 5-6. A brief review on both subjects is presented here.

We consider the interaction between a relativistic fireball and material assumed to be still located just outside the progenitor star. Only a small fraction of the expected mass is sufficient to efficiently decelerate the fireball, leading to dissipation of most of its kinetic energy. Since the scattering optical depths are still large at distances comparable to the progenitor radius, the dissipated energy is trapped in the system, accelerating it to relativistic velocities. The process resembles the birth of another fireball at radii R~1e11 cm, not far from the transparency radius, and with a starting bulk Lorentz factors Gamma_c~10. As seen in the observer frame, this "re--generated" fireball appears collimated within an angle theta_j=1/Gamma_c. If the central engine works intermittently, the funnel can, at least partially, refill and the process can repeat itself. We discuss how this idea can help solving some open issues of the more conventional internal shock scenario for interpreting the Gamma-Ray Burst properties.

We study the spectral and energetics properties of 47 long-duration gamma-ray bursts (GRBs) with known redshift, all of them detected by the Swift satellite. Due to the narrow energy range (15-150 keV) of the Swift-BAT detector, the spectral fitting is reliable only for fitting models with 2 or 3 parameters. As high uncertainty and correlation among the errors is expected, a careful analysis of the errors is necessary. We fit both the power law (PL, 2 parameters) and cut--off power law (CPL, 3 parameters) models to the time-integrated spectra of the 47 bursts, and present the corresponding parameters, their uncertainties, and the correlations among the uncertainties. The CPL model is reliable only for 29 bursts for which we estimate the nuf_nu peak energy Epk. For these GRBs, we calculate the energy fluence and the rest- frame isotropic-equivalent radiated energy, Eiso, as well as the propagated uncertainties and correlations among them. We explore the distribution of our homogeneous sample of GRBs on the rest-frame diagram E'pk vs Eiso. We confirm a significant correlation between these two quantities (the "Amati" relation) and we verify that, within the uncertainty limits, no outliers are present. We also fit the spectra to a Band model with the high energy power law index frozen to -2.3, obtaining a rather good agreement with the "Amati" relation of non-Swift GRBs.

We study 7 Gamma Ray Bursts (GRBs), detected both by the BATSE instrument, on-board the Compton Gamma Ray Observatory, and by the Wide Field Camera (WFC), on-board BeppoSAX. These bursts have measured spectroscopic redshifts and are a sizeable fraction of the bursts defining the correlation between the peak energy E_peak (i.e. the peak of the vFv spectrum) and the total prompt isotropic energy E_iso (the so called "Amati" relation). Recent theoretical interpretations of this correlation assume that black-body emission dominates the time resolved spectra of GRBs, even if, in the time integrated spectrum, its presence may be hidden by the change of its temperature and by the dilution of a possible non-thermal power law component. We perform a time resolved spectral analysis, and show that the sum of a power-law and a black-body gives acceptable fits to the time dependent spectra within the BATSE energy range, but overpredicts the flux in the WFC X-ray range. Moreover, a fit with a cutoff power-law plus a black-body is consistent with the WFC data, but the black-body component contributes a negligible fraction of the total flux. On the contrary, we find that fitting the spectra with a Band model or a simple cutoff power-law model yields an X-ray flux and spectral slope which well matches the WFC spectra.

In their recent paper, Campana et al. (2007) found that 5 bursts, among those detected by Swift, are outliers with respect to the E_peak-E_gamma ("Ghirlanda") correlation. We instead argue that they are not.

The early X-ray light curve of Gamma Ray Bursts (GRBs) is complex, and shows a typical steep-flat-steep behaviour. The time T_a at which the flat (plateau) part ends may bear some important physical information, especially if it plays the same role of the so called jet break time t_jet. To this aim, stimulated by the recent analysis of Willingale et al., we have assembled a sample of GRBs of known redshifts, spectral parameters of the prompt emission, and T_a. By using T_a as a jet angle indicator, and then estimating the collimation corrected prompt energetics, we find a correlation between the latter quantity and the peak energy of the prompt emission. However, this correlation has a large dispersion, similar to the dispersion of the Amati correlation and it is not parallel to the Ghirlanda correlation. Furthermore, we show that the correlation itself results mainly from the dependence of the jet opening angle on the isotropic prompt energy, with the time T_a playing no role, contrary to what we find for the jet break time t_jet. We also find that for the bursts in our sample T_a weakly correlates with E_iso of the prompt emission, but that this correlation disappears when considering all bursts of known redshift and T_a. There is no correlation between T_a and the isotropic energy of the plateau phase.

Recently, a tight correlation among three quantities that characterize the prompt emission of long Gamma-Ray Bursts (GRBs) with known redshift z, was discovered (Firmani et al. 2006). We use this correlation to construct the Hubble diagram (HD) with a sample of 19 GRBs in the broad range of z=0.17-4.5, and carry out a full statistical analysis to constrain cosmological parameters (CPs). To optimally solve the problem of circularity, a Bayesian approach is applied. The main result is that the concordance LambdaCDM cosmology is fully consistent with the GRB data at the level of several tests. If we assume the Lambda cosmology, then we find Om_M=0.31^{+0.09}_{-0.08} and Om_Lambda=0.80^{+0.20}_{-0.30}$ (1sigma); the flat-geometry case is within 1sigma. Assuming flatness, we find Om_M=0.29^{+0.08}_{-0.06}, and fixing Om_M=0.28, we obtain a dark energy equation of state parameter w=-1.07^{+0.25}_{-0.38}, i.e. the ambdaCDM model (w=-1) is within 1sigma. Given the low number of usable GRBs we cannot yet constrain well the possible evolution of w=w(z). However, the case w(z)=-1 (LambdaCDM) is consistent at the 68.3% CL with GRBs. It is shown also how a broad range of z's in the used sample improves the determination of CPs from the HD, which is the case of GRBs as distance indicators.

The flat decay phase in the first 1e2-1e4 seconds of the X-ray light curve of Gamma Ray Bursts (GRBs) has not yet found a convincing explanation. The fact that the optical and X-ray lightcurves are often different, with breaks at different times, makes contrived any explanation based on the same origin for both the X-ray and optical fluxes. We here assume that the central engine can be active for a long time, producing shells of decreasing bulk Lorentz factors Gamma. We also assume that the internal dissipation of these late shells produces a continuous and smooth emission (power-law in time), usually dominant in X-rays and sometimes in the optical. When Gamma of the late shells is larger than 1/theta_j, where theta_j is the jet opening angle, we see only a portion of the emitting surface. Eventually, Gamma becomes smaller than 1/theta_j, and the entire emitting surface is visible. Thus there is a break in the light curve when Gamma=1/theta_j, which we associate to the time at which the plateau ends. After the steeply decaying phase which follows the early prompt, we see the sum of two emission components: the "late-prompt" emission (due to late internal dissipation), and the "real afterglow" emission (due to external shocks). A variety of different optical and X-ray light curves are then possible, explaining why the X-ray and the optical light curves often do not track each other (but sometimes do), and often they do not have simultaneous breaks.

We test the spectral-energy correlation including the new bursts detected (mostly) by Swift with firm measurements of their redshifts and peak energy. The problem of identifying the jet breaks is discussed in the complex and multibreak/flaring X-ray light curves observed by Swift. We use the optical data as the most reliable source for the identification of the jet break, since the X-ray flux may be produced by a mechanism different from the external shocks between the fireball and the circumburst medium, which are responsible for the optical afterglow. We show that the presence of an underlying SN event in XRF 050416A requires a break to occur in the afterglow optical light curve at around the expected jet break time. The possible presence of a jet break in the optical light curve of GRB 050401 is also discussed. We point out that, for measuring the jet break, it is mandatory that the optical light curve extends after the epoch where the jet break is expected. The interpretation of the early optical breaks in GRB 050922C and GRB 060206 as jet breaks is controversial because they might instead correspond to the flat-to-steep decay transition common in the early X-ray light curves. All the 16 bursts coming from Swift are consistent with the E_p-E_gamma and E_p-E_iso-t_jet correlation. No outlier is found to date. Moreover, the small dispersion of these correlations, confirmed also by the Swift bursts, strengthens the case of using GRBs as standard candles.

Ghirlanda et al. (2004) derived the collimation-corrected energy E_gamma for a sample of 15 bursts under the assumption of a homogeneous circumburst density profile. They found a correlation (the so-called Ghirlanda correlation) between E_gamma and the rest frame peak energy of the nuF_nu prompt spectrum (E_p). Nava et al. (2006) showed that, assuming a circumburst density distribuited with a r^-2 wind profile, the Ghirlanda correlation remains tight and becomes linear. This implies that: i) it remains linear also in the comoving frame, no matter the distribution of bulk Lorentz factors, ii) it entails that different bursts have the same number of relevant photons. We have updated these findings including recently detected bursts (21 in total), stressing the two important implications.

We studied the behaviour of the optical afterglow lightcurves of a sample of 24 Gamma--Ray Bursts (GRBs) with known redshift and published estimates of the optical extinction in the source frame, detected before the SWIFT satellite launch. We found an unexpected clustering of the optical luminosities at 12 hours in the source frame. The distribution of the optical luminosities is narrower than the distribution of X-ray luminosities at the same time. Few (3) bursts stand apart from the main optical distribution, being fainter by a factor of about 15. We also analysed the optical luminosities of the SWIFT burst with known redshift finding that the luminosity distribution is similar to the pre SWIFT GRBs one, even if they have a different mean redshift. These results can suggest the existence of a family of intrinsically optically under--luminous dark GRBs.

Gamma Ray Bursts (GRBs) are among the most powerful sources in the Universe: they emit up to 10^54 erg in the hard X-ray band in few tens of seconds. The cosmological origin of GRBs has been confirmed by several spectroscopic measurements of their redshifts, distributed in the range 0.1-6.3. These two properties make GRBs very appealing to investigate the far Universe. The energetics implied by the observed fluences and redshifts span at least four orders of magnitudes. Therefore, at first sight, GRBs are all but standard candles. But there are correlations among some observed quantities which allow us to know the total energy or the peak luminosity emitted by a specific burst with a great accuracy. Through these correlations, GRBs become "known" candles to constrain the cosmological parameters. One of these correlation is between the rest frame peak spectral energy E_peak and the total energy emitted in gamma--rays E_gamma, properly corrected for the collimation factor. Another correlation, discovered very recently, relates the total GRB luminosity L_iso, its peak spectral energy E_peak and a characteristic timescale T_0.45, related to the variability of the prompt emission. It is based only on prompt emission properties, it is completely phenomenological, model independent and assumption--free. The constraints found through these correlations on the Omega_M and Omega_Lambda parameters are consistent with the concordance model. The present limited sample of bursts and the lack of low redshift events, necessary to calibrate these correlations, makes the cosmological constraints obtained with GRBs still large compared to those obtained with other cosmological probes (e.g. SNIa or CMB). However, the newly born field of GRB--cosmology is very promising for the future.

Tight constraints on cosmological parameters can be obtained with standard candles spanning a range of redshifts as large as possible. We propose to treat SN Ia and long Gamma-Ray Bursts (GRBs) as a single class of candles. Taking advantage of the recent release of the Supernova Legacy Survey and {\it the recent finding of a tight correlation among the energetics and other prompt gamma-ray emission properties of GRBs}, we are able to standardize the luminosities/energetics of both classes of objects. In this way we can jointly use GRB and SNIa as cosmological probes to constrain Omega_m and Omega_L and the Dark Energy equation of state parameters through the same Bayesian method that we have, so far, applied to GRBs alone. Despite the large disparity in number (115 SNIa versus 19 GRBs) we show that the constraints on Omega_m and Omega_L are greatly improved by the inclusion of GRBs. More importantly, the result of the combined sample is in excellent agreement with the Lambda-CDM concordance cosmological model and does not require an evolving equation of state for the Dark Energy.

GRB 980425 and GRB 031203 are apparently two outliers with respect to the correlation between the isotropic equivalent energy E_iso emitted in the prompt radiation phase and the peak frequency E_peak of the spectrum in a vF(v) representation (the so-called Amati relation). We discuss if these two bursts are really different from the others or if their location in the E_iso-E_peak plane is the result of other effects, such as viewing them off-axis, or through a scattering screen, or a misinterpretation of their spectral properties. The latter case seems particularly interesting after GRB 060218, that, unlike GRB 031203 and GRB 980425, had a prompt emission detected both in hard and soft X-rays which lasted ~2800 seconds. This allowed to determine its E_peak and total emitted energy. Although it shares with GRB 031203 the total energetics, it is not an outlier with respect to the Amati correlation. We then investigate if a hard-to-soft spectral evolution in GRB 031203 and GRB 980425, consistent with all the observed properties, can give rise to a time integrated spectrum with an E_peak consistent with the Amati relation.

We report the discovery of a correlation among three prompt emission properties of GRBs. These are the isotropic peak luminosity L_iso, the peak energy of the time-integrated prompt emission spectrum E_pk, and the ``high signal" timescale T_0.45, previously used to characterize the variability behavior of bursts. In the rest frame of the source the found correlation reads L_iso\propto E_pk^1.62 T_0.45^-0.49. We find other strong correlations, but at the cost of increasing the number of variables, involving the variability and the isotropic energy of the prompt emission. With respect to the previous tight correlations found in GRBs the newly found correlation does not require any information from the afterglow phase of the bursts, nor any model-dependent assumption. In the popular scenario in which we are receiving beamed radiation originating in a fireball pointing at us, the found correlation preserves its form in the comoving frame. This helps to explain the small scatter of the correlation, and underlines the role of the local brightness (i.e. the brightness of the visible fraction of the fireball surface). This correlation has been found for 19 objects, and it is hard to establish if any selection bias affects it. Its connection with the prompt local brightness is promising, but a solid physical understanding is still to be found. Despite all that, we find that some properties of the correlation, which we discuss, support its true existence, and this has important implications for the GRB physics. Furthermore, it is possible to use such correlation as an accurate redshift estimator, and its tightness will allow us to use it as a tool to constrain the cosmological parameters (abridged)

As soon as it was realized that long GRBs lie at cosmological distances, attempts have been made to use them as cosmological probes. Besides their use as lighthouses, a task that presents mainly the technological challenge of a rapid deep high resolution follow-up, researchers attempted to find the Holy Grail: a way to create a standard candle from GRB observables. We discuss here the attempts and the discovery of the Ghirlanda correlation, to date the best method to standardize the GRB candle. Together with discussing the promises of this method, we will underline the open issues, the required calibrations and how to understand them and keep them under control. Even though GRB cosmology is a field in its infancy, ongoing work and studies will clarify soon if and how GRBs will be able to keep up to the promises.

We studied the optical afterglows of the 24 pre-Swift Gamma-Ray Bursts with known spectroscopic redshift and published estimates of the optical extinction in the source frame. We find an unexpected clustering of the optical afterglow luminosities measured 12 hours (source frame time) after the trigger. For 21 out of 24 bursts, the distribution of the optical luminosities is narrower than the distribution of the X-ray luminosities, and even narrower than the distribution of the ratio between the monochromatic optical luminosities and the total isotropic emitted prompt energy. Three bursts stand apart from the distribution of the other sources, being underluminous by a factor ~15. We compare this result with the somewhat analogous result concerning the luminosity of the X-ray afterglows studied by Gendre & Boer. For all our GRBs we construct the optical to X-ray spectral energy distribution. For all but a minority of them, the optical and the X-ray emissions are consistent with being produced by the same radiation process. We discuss our results in the framework of the "standard" external shock synchrotron model. Finally, we consider the behavior of the first GRBs of known redshifts detected by Swift. We find that these Swift GRBs entirely confirm our findings.

The use of Gamma Ray Bursts as ``standard candles'' has been made possible by the recent discovery of a very tight correlation between their rest frame intrinsic properties. This correlation relates the GRB prompt emission peak spectral energy E_peak to the energy E_gamma corrected for the collimation angle theta_jet of these sources. The possibility to use GRBs to constrain the cosmological parameters and to study the nature of Dark Energy are very promising.

We present the constraints on the cosmological parameters obtained with the $E_{\rm peak}$--$E_{\gamma}$ correlation found with the most recent sample of 19 GRBs with spectroscopically measured redshift and well determined prompt emission spectral and afterglow parameters. We compare our results obtained in the two possible uniform jet scenarios, i.e. assuming a homogeneous density profile (HM) or a wind density profile (WM) for the circumburst medium. Better constraints on $\Omega_{M}$ and $\Omega_{\Lambda}$ are obtained with the (tighter) $E_{\rm peak}$--$E_{\gamma}$ correlation derived in the wind density scenario. We explore the improvements to the constraints of the cosmological parameters that could be reached with a large sample, $\sim$ 150 GRBs, in the future. We study the possibility to calibrate the slope of these correlations. Our optimization analysis suggests that $\sim 12$ GRBs with redshift $z\in(0.9,1.1)$ can be used to calibrate the $E_{\rm peak}$--$E_{\gamma}$ with a precision better than 1%. The same precision is expected for the same number of bursts with $z\in(0.45,0.75)$. This result suggests that we do not necessarily need a large sample of low z GRBs for calibrating the slope of these correlations.

Recently, Liang & Zhang (2005) found a tight correlation involving only observable quantities, namely the isotropic emitted energy $E_{\gamma,iso}$, the energy of the peak of the prompt spectrum $E^\prime_{p}$, and the jet break time $t^\prime_{j}$. This phenomenological correlation can have a first explanation in the framework of jetted fireballs, whose semiaperture angle $\theta_{j}$ is measured by the jet break time $t^\prime_{j}$. By correcting $E_{\gamma, iso}$ for the angle $\theta_{j}$ one obtains the so called Ghirlanda correlation linking the collimation corrected energy $E_\gamma$ and $E^\prime_{p}$. There are two ways to derive $\theta_{j}$ from $t^\prime_{j}$ in the standard scenario, corresponding to an homogeneous or to a wind-like circumburst medium. We show that the Ghirlanda correlation with a wind-like medium is as tight as (if not tighter) than the Ghirlanda correlation found in the case of an homogeneous medium. There are hence two Ghirlanda correlations, both entirely consistent with the phenomenological Liang & Zhang relation. We consider the difference between the observed correlations and the ones one would see in the comoving frame (i.e. moving with the same bulk Lorentz factor of the fireball). Since both $E_{p}$ and $E_\gamma$ transform in the same way, the wind-like Ghirlanda relation, which is linear, remains linear also in the comoving frame, no matter the distribution of bulk Lorentz factors. Instead, in the homogeneous density case, one is forced to assume the existence of a strict relation between the bulk Lorentz factor and the total energy, which in turn put constraints on the radiation mechanisms of the prompt emission. The wind-like Ghirlanda correlation, being linear, corresponds to different bursts having the same number of photons.

High-resolution LCDM cosmological N-body simulations are used to study the properties of galaxy-size dark halos in different environments (cluster, void, and "field"). Halos in clusters and their surroundings have a median spin parameter ~1.3 times lower, and tend to be more spherical and to have less aligned internal angular momentum than halos in voids and the field. For halos in clusters the concentration parameters decrease on average with mass with a slope of ~0.1; for halos in voids these concentrations do not change with mass. For masses <5 10^11 M_sh^-1, halos in clusters are on average ~30-40% more concentrated and have ~2 times higher central densities than halos in voids. When comparing only parent halos, the differences are less pronounced but they are still significant. The Vmax-and Vrms-mass relations are shallower and more scattered for halos in clusters than in voids, and for a given Vmax or Vrms, the mass is smaller at z=1 than at z=0 in all the environments. At z=1, the differences in the halo properties with environment almost dissapear, suggesting this that the differences were stablished mainly after z~1. The halos in clusters undergo more dramatic changes than those in the field or the voids. The differences with environment are owing to (i) the dependence of halo formation time on environment, and (ii) local effects as tidal stripping and the tumultuos histories that halos suffer in high-density regions. We calculate seminumerical models of disk galaxy evolution in halos with the properties found for the different environments. For a given disk mass, the galaxy disks have higher surface density, larger Vd,max and secular bulge-to-disk ratio, lower gas fraction, and are redder as one goes from cluster to void environments, in rough agreement with observations. (abridged)

[Abridged] Very recently, relations between the peak energy of Gamma-Ray Burst burst spectra, the isotropic-equivalent energy of the burst, and the radiated energy of the burst have been found. In a way that is exactly analogous to the way in which the relation between the peak luminosity and the rate of decline of the light curve of Type Ia supernovae can be used to make Type Ia supernovae excellent standard candles for cosmology, so too, the relations between Epk, Eiso, and Egamma point toward a methodology for using GRBs as excellent standard candles for cosmology. In addition, GRBs occur over the broad redshift range from z=0.1 to at least z=4.5, and both they and their afterglows are easily detectable out to z > 8. Thus GRBs show great promise as cosmological ``yardsticks'' to measure the rate of expansion of the universe over time, and therefore the properties of dark energy (i.e., Omega_M, Omega_Lambda, w_0, and w_a). We describe a concept for a possible MIDEX-class mission dedicated to using GRBs to constrain the properties of dark energy that would obtain these quantities for > 800 bursts in the redshift range 0.1 >~ z <~ 10 during a 2-year mission. This burst sample would enable both Omega_M and w_0 to be determined to +/- 0.07 and +/- 0.06 (68% CL), respectively, and w_a to be significantly constrained.

We probe the existence of the E_peak-E_iso correlation in long GRBs using a sample of 442 BATSE bursts with known E_peak and with redshift estimated through the lag-luminosity correlation. This sample confirms that the rest frame peak energy is correlated with the isotropic equivalent energy. The distribution of the scatter of the points around the best fitting line is similar to that obtained with the 27 bursts with spectroscopic redshifts. We interpret the scatter in the E_peak-E_iso plane as due to the opening angle distribution of GRB jets. By assuming that the collimation corrected energy correlates with E_peak we can derive the observed distribution of the jet opening angles, which turns out to be log-normal with a peak value of ~6.5 degrees.

The isotropic luminosity function (LF) and formation rate history (FRH) of long GRBs is by the first time constrained by using jointly both the observed GRB peak-flux and redshift distributions. Our results support an evolving LF and a FRH that keeps increasing after z=2. We discuss some interesting implications related to these results.

Apparently, Gamma-Ray Bursts (GRBs) are all but standard candles. Their emission is collimated into a cone and the received flux depends on the cone aperture angle. Fortunately we can derive the aperture angle through an achromatic steepening of the lightcurve of the afterglow, and thus we can measure the "true" energetics of the prompt emission. Ghirlanda et al. (2004) found that this collimation-corrected energy correlates tightly with thefrequency at which most of the radiation of the prompt is emitted. Through this correlation we can infer the burst energy accurately enough for a cosmological use. Using the best known 15 GRBs we find very encouraging results that emphasize the cosmological GRB role. Probing the universe with high accuracy up to high redshifts, GRBs establish a new insight on the cosmic expanding acceleration history and accomplish the role of "missing link" between the Cosmic Microwave Background and type Ia supernovae, motivating the most optimistic hopes for what can be obtained from the bursts detected by SWIFT.

The recently discovered correlation between the rest frame GRB peak spectral energy $E_{\rm peak}$ and the collimation corrected energy $E_\gamma$ in long GRBs is potentially very important, yet awaits confirmation from an independent sample. It may help to shed light on the radiation mechanism of the prompt GRB phase and on the way -- and in which form -- the energy is released from the central engine. We here present some additional evidence for the correlation (two new bursts) and re-derive the best-fit parameters. The tightness of the correlation is confirmed (sigma=0.1 dex). We show that this correlation allows us, for the first time, to use GRBs as cosmological probes to constrain the expansion history of the universe.

After many years of speculation, recent observations have confirmed the association of gamma-ray bursts with core-collapse supernova explosions from massive stars. This association carries with it important consequences. The burst relativistic jet has to propagate through the cold dense stellar material before it reaches the transparency radius and the burst photons are produced. This propagation is likely to affect the initial properties of the jet, shaping it and changing its energy composition. The variability injected at the base of the jet is also likely to be erased by the jet-star interaction. Despite this, GRBs seem to have remarkably predictable properties once the radiative phase sets in, as emphasized by the recent discovery of several tight correlation between spectral, geometric and energetic properties of the jet. In this contribution we discuss the jet interaction with the star, emphasizing its time-dependent properties and the resulting energy distribution. We finally emphasize the surprising predictability of jet and radiation properties outside the star and underline its implication for standardizing the GRB candle.

We derive the peak luminosity - peak energy (L_iso - E_peak) correlation using 22 long Gamma-Ray Bursts (GRBs) with firm redshift measurements. We find that its slope is similar to the correlation between the time integrated isotropic emitted energy E_iso and E_peak (Amati et al. 2002). For the 15 GRBs in our sample with estimated jet opening angle we compute the collimation corrected peak luminosity L_gamma, and find that it correlates with E_peak. This has, however, a scatter larger than the correlation between E_peak and E_gamma (the time integrated emitted energy, corrected for collimation; Ghirlanda et al. 2004), which we ascribe to the fact that the opening angle is estimated through the global energetics. We have then selected a large sample of 442 GRBs with pseudo--redshifts, derived through the lag-luminosity relation, to test the existence of the L_iso-E_peak correlation. With this sample we also explore the possibility of a correlation between time resolved quantities, namely L_iso,p and the peak energy at the peak of emission E_peak,p.

We present a new method aimed to handle long Gamma-Ray Burst (GRBs) as cosmic rulers. The recent discovery of a tight correlation between the collimation corrected GRB energy and the peak of the gamma-ray spectrum has opened the possibility to use GRBs as a new category of standard candles. Unfortunately, because of the lack of low-z GRBs, up to now this correlation is obtained from high-z GRBs with the consequence that it depends on the cosmological parameters we pretend to constrain. Hopefully this circularity problem will be solved when, in a few years, the low-z GRB sample will be increased enough. In the meanwhile we present here a new Bayesian method that eases the aforesaid circularity problem, and allows to introduce new constrains on the cosmological (Om,OL) diagram as well as to explore the universe kinematics up to z~3. The method we propose offers the further advantage to make handy the problem of the (Om,OL) loitering line singularity which inevitably appears when standard candles with z>2 are used. The combination of GRB with SN Ia data makes the popular LambdaCDM cosmology more consistent with the Hubble diagram at a 68% confidence level. For a flat cosmology we find Om=0.28\pm0.03 for the combined GRB+SN Ia data set. Correspondingly, the transition redshift between cosmic deceleration and acceleration is z_T=0.73\pm0.09, slightly larger than the value found by considering SNe Ia alone. We briefly discuss our results also in terms of non--LambdaCDM dark energy models.

The best measure of the Universe should be done using a standard "ruler" at any redshift. Type Ia Supernovae (SN Ia) probe the universe up to z$\sim$1.5, while the Cosmic Microwave Background (CMB) primary anisotropies concern basically $z\sim$1000. Apparently, Gamma--Ray Bursts (GRBs) are all but standard candles. However, their emission is collimated and the collimation--corrected energy correlates tightly with the frequency at which most of the radiation of the prompt is emitted, as found by Ghirlanda et al. (2004). Through this correlation we can infer the burst energy accurately enough to probe the intermediate redshift ($z<10$) Universe. Using the best known 15 GRBs we find very encouraging results that emphasize the cosmological GRB role. A combined fit with SN Ia yields $\Omega_{\rm M}=0.37\pm0.10$ and $\Omega_{\Lambda}=0.87\pm 0.23$. Assuming in addition a flat Universe, the parameters are constrained to be $\Omega_{\rm M}=0.29\pm0.04$ and $\Omega_{\Lambda}=0.71\pm 0.05$. GRBs accomplish the role of "missing link" between SN Ia and CMB primary anisotropies. They can provide a new insight on the cosmic effects of dark energy, complementary to the one supplied by CMB secondary anisotropies through the Integrated Sachs Wolfe effect. The unexpected Standard Candle cosmological role of GRBs motivates the most optimistic hopes for what can be obtained when the GRB-dedicated satellite, Swift, will be launched.

We constrain the isotropic luminosity function (LF) and formation rate of long gamma-ray bursts (GRBs) by fitting models 'jointly' to both the observed differential peak flux and redshift distributions. We find evidence supporting an evolving LF, where the luminosity scales as (1+z)^delta with an optimal delta of 1.0\pm 0.2. For a single power-law LF, the best slope is ~ -1.57 with an upper luminosity of 10^53.3 erg/s, while the best slopes for a double power-law LF are approximately -1.6 and -2.6 with a break luminosity of 10^52.7 erg/s. Our finding implies a jet model intermediate between the universal structured epsilon(theta) proportional to theta^-2 model and the quasi-universal Gaussian structured model. For the uniform jet model our result is compatible with an angle distribution between 2 and 15 grades. Our best constrained GRB formation rate histories increase from z=0 to z=2 by a factor of ~30 and then continue increasing slightly. We connect these histories to that of the cosmic star formation history, and compare with observational inferences up to z~6. GRBs could be tracing the cosmic rates of both the normal and obscured star formation regimes. We estimate a current GRB event rate in the Milky Way of ~5 10^-5 yr^-1, and compare it with the birthrate of massive close WR+BH binaries with orbital periods of hours. The agreement is rather good suggesting that these systems could be the progenitors of the long GRBs.

We study the structure and substructure of halos obtained in N-body simulations for a Lambda Cold Dark Matter (LCDM) cosmology with non-Gaussian initial conditions (NGICs). The initial statistics are lognormal in the gravitational potential field with positive (LNp) and negative (LNn) skewness; the sign of the skewness is conserved by the density field, and the power spectrum is the same for all the simulations. Our aim is not to test a given non-Gaussian statistics, but to explore the generic effect of positive- and negative-skew statistics on halo properties. From our low-resolution simulations, we find that LNp (LNn) halos are systematically more (less) concentrated than their Gaussian counterparts. This result is confirmed by our Milky Way- and cluster-sized halos resimulated with high-resolution. In addition, they show inner density profiles that depend on the statistics: the innermost slopes of LNp (LNn) halos are steeper (shallower) than those obtained from the corresponding Gaussian halos. A subhalo population embedded in LNp halos is more susceptible to destruction than its counterpart inside Gaussian halos. On the other hand, subhalos in LNn halos tend to survive longer than subhalos in Gaussian halos. The spin parameter probability distribution of LNp (LNn) halos is skewed to smaller (larger) values with respect to the Gaussian case. Our results show how the statistics of the primordial density field can influence some halo properties, opening this the possibility to constrain, although indirectly, the primordial statistics at small scale.

High-resolution numerical simulations were performed to study the structure and substructure of Milky Way- and cluster-sized halos in a LCDM cosmology with self-interacting (SI) dark particles, where the particle cross section, +AFw-sig, is assumed constant or inversely proportional to the relative velocity. We conclude that the cuspy halo problem at galaxy scales of the LCDM cosmogony can be solved in the latter case. In this case, the inner density of galaxy-cluster halos is only slightly affected by the SI, resulting in agreement with observational constraints. At the same time, the subhalo population in galaxy and cluster SI-CDM halos remains roughly similar to that seen on collisionless CDM halos.

The study of formation and evolution of galaxies is reviewed, making emphasis on the physical factors which are important to understand the origin of the galaxy Hubble sequence. We concentrate on predictions of the hierarchical Cold Dark Matter (CDM) scenario and the confrontation with observations. The mass assembling of the CDM halos, the baryonic processes within them, and the evolution of disks and spheroids are described. The successes and shortcomings are discussed. Disk evolution seems to be a quiescent and extended process driven by the cosmological initial conditions, while spheroids are formed probably in violent events, where several astrophysical processes are competing.

For a compiled sample of disk galaxies with available photometry (B and K bands), velocity line-widths and HI integral fluxes, several parameters which trace the luminous, baryonic and dark matter contents were inferred. We investigated how these parameters do vary with different galaxy properties, and confronted the results with predictions of galaxy evolutionary models in the context of the LCDM cosmogony. The ratio of disk-to-total maximum circular velocity, vd/vt, depends mainly on the central disk surface density Sig_d (or surface brightness, SB), increasing roughly as Sig_d^0.15. While a fraction of high SB galaxies have a vd/vt ratio corresponding to the maximum disk solution, the low SB are completely dark matter dominated. The trend is similar for the models, although they have slightly smaller vd/vt ratios than observations. An analysis of residuals of the vd/vt-Sig_d relation shows that vd/vt tends to decrease as the galaxy is redder, more luminous (massive), and of earlier type. The models allow us to explain the physics of these results, which imply a connexion between halo structure and luminous properties. The estimated dynamical-to-baryon mass and dynamical mass-to-light ratios at a given radius, for observations and models, decrease with the SB and do not correlate with the galaxy scale, contrary to what has been reported in previous works, based on the analysis of rotation curve shapes. We discuss this difference and state the importance to solve the controversy on whether the dark and luminous contents in disk galaxies depend on SB or luminosity. The broad agreement between models and observations favors the LCDM scenario. However, the excess of dark matter inside the optical region of disk galaxies remains as the main difficulty. (abridged)

The kinetic energy dissipation rate in the turbulent ISM of disk galaxies is a key ingredient in galaxy evolution models since it determines the effectiveness of large-scale star formation (SF) feedback. Using magneto-hydro-dynamic simulations, we find that the ISM dissipates efficiently the turbulent kinetic energy injected by sources of stellar nature. Thus, the SF process may be self-regulated by an energy balance only at the level of the disk ISM. The use of the self-regulation SF mechanism in galaxy evolutionary models, where disks form inside growing Cold Dark Matter halos, allows to predict the SF history of disk galaxies, including the Milky Way and the solar neighborhood, as well as the contribution of the whole population of disk galaxies to the cosmic SF history. The results are encouraging.

We study the structure of Milky Way (MW)- and cluster-sized halos in a Lambda Cold Dark Matter (CDM) cosmology with self-interacting (SI) dark particles. The cross section per unit of particle mass has the form sigma = sig_0(1/v_100)^alpha, where sig_0 is a constant in units of cm^2/gr and v_100 is the relative velocity in units of 100 km/s. Different values for sigma with alpha= 0 or 1 were used. For small values of sigma = const. (sig_0<0.5), the core density of the halos at z=0 is typically higher at a given mass for lower values of sig_0 or, at a given sig_0, for lower masses. For values of sig_0 as high as 3.0, the halos may undergo the gravothermal catastrophe before z=0. When alpha = 1, the core density of cluster- and MW-sized halos is similar. Using sigma = 0.5-1.0x(1/v_100), our predictions agree with the central densities and the core scaling laws of halos both inferred from the observations of dwarf and LSB galaxies and clusters of galaxies. The cumulative Vmax-functions of subhalos in MW-sized halos with (sig_0,alpha) = (0.1,0.0), (0.5,0.0) and (0.5,1.0) agree roughly with observations (luminous satellites) for Vmax > 30 km/s, while at Vmax = 20 km/s the functions are a factor 5-8 higher, similar to the CDM predictions. The halos with SI have slightly more specific angular momentum at a given mass shell and are rounder than their CDM counterparts. We conclude that the introduction of SI particles with sigma \propto 1/v_100 may remedy the cuspy core problem of the CDM cosmogony, while the subhalo population number remains similar to that of the CDM halos.

We propose a new numerical procedure to simulate a single dark halo of any size and mass in a hierarchical framework coupling the extended Press-Schechter formalism (EPSF) to N-body simulations. The procedure consists of assigning cosmological initial conditions to the particles of a single halo with a EPSF technique and following only the dynamical evolution using a serial N-body code. The computational box is fixed with a side of $0.5 h^{-1}$ Mpc. This allows to simulate galaxy cluster halos using appropriate scaling relations, to ensure savings in computing time and code speed. The code can describe the properties of halos composed of collisionless or collisional dark matter. For collisionless Cold Dark Matter (CDM) particles the NFW profile is reproduced for galactic halos as well as galaxy cluster halos. Using this numerical technique we study some characteristics of halos assumed to be isolated or placed in a cosmological context in presence of weak self-interacting dark matter: the soft core formation and the core collapse. The self-interacting dark matter cross section per unit mass is assumed to be inversely proportional to the particle collision velocity: $\sigma/m_{x} \propto 1/v$.

We present high resolution Halpha rotation curves of 4 late-type dwarf galaxies and 2 low surface brightness galaxies (LSB) for which accurate HI rotation curves are available from the literature. Observations are carried out at Telescopio Nazionale Galileo (TNG). For LSB F583-1 an innovative dispersing element was used, the Volume Phase Holographic (VPH) with a dispersion of about 0.35 A/pxl. We find good agreement between the Halpha data and the HI observations and conclude that the HI data for these galaxies suffer very little from beam smearing. We show that the optical rotation curves of these dark matter dominated galaxies are best fitted by the Burkert profile. In the centers of galaxies, where the N-body simulations predict cuspy cores and fast rising rotation curves, our data seem to be in better agreement with the presence of soft cores.

We discuss the shape and decomposition of rotation curves (RCs) of galaxies formed within growing cold dark matter halos. The outer RC shape correlates mainly with the surface brightness (SB), the luminous mass fraction, fd, and the bulge fraction. In order the shapes of RC depend significantly on luminosity, fd should be a strong function of mass (feedback). For the preferred values of fd (\lesssim 0.03), the high SB models can be maximum disks only when the halos have a shallow core. The low SB models are sub-maximum disks. The residuals of the baryonic Tully-Fisher (TF) and disk mass-radius relations show a clear anti-correlation among them, but when one passes to the TF and luminosity-radius relations, the anti-correlation almost disappears. Therefore, the observed lack of correlation among the residuals of the last two relations should not be interpreted as an evidence of sub-maximal disks.

We find that the density distribution of high and low surface brightness (SB) galaxy disks formed within LCDM halos under the assumption of detailed angular momentum conservation is in rough agreement with observations. The luminous-to-dark matter ratio of the model galaxies increases with SB. The lowest SB models tend to be minimum disks, but the high SB models hardly attain the maximum disk solution. With the introduction of shallow cores in the halos, high SB models become maximum disks. The shallow cores also help to improve the inner SB profiles of bulge-less low SB models and the zero-point of the Tully-Fisher relation. The models predict well this relation and its scatter, as well as the small correlation among the residuals of this and the luminosity-radius relation, in spite of the dependence of the rotation curve shapes on SB.

We present high spatial resolution H\alpha rotation curves of late-type dwarf and LSB galaxies. From our analysis we find good agreement between our H\alpha data and the HI observations taken from the literature, concluding that the HI rotation curves for these galaxies suffer very little from beam smearing. A preliminary analysis of our data rules out the CDM model in the inner regions of these galaxies.

This work is focused on the preliminary results of the observations of H\alpha rotation curves for some of the late type dwarf and LSB galaxies carried out at the TNG telescope. In light of the observational data and of the N-body simulations we have performed recently, we discuss some of the implications on the nature of the DM particles and the formation of the dark haloes.

We use a cosmological galactic evolutionary approach to model the Milky Way. A detailed treatment of the mass aggregation and dynamical history of the growing dark halo is included, together with a self consistent physical treatment for the star formation processes within the growing galactic disc. This allows us to calculate the temporal evolution of star and gas surface densities at all galactic radii, in particular, the star formation history (SFH) at the solar radius. A large range of cosmological mass aggregation histories (MAHs) is capable of producing a galaxy with the present day properties of the Milky Way. The resulting SFHs for the solar neighbourhood bracket the available observational data for this feature, the most probable MAH yielding the optimal comparison with these observations. We also find that the rotation curve for our Galaxy implies the presence of a constant density core in its dark matter halo.

(Abridged) Using high-resolution cosmological N-body simulations, we study how the density profiles of dark matter halos are affected by the filtering of the density power spectrum below a given scale length and by the introduction of a thermal velocity dispersion. In the warm dark matter (WDM) scenario, both the free-streaming scale, R_f, and the velocity dispersion, v_w, are determined by the mass m_w of the WDM particle. We find that v_w is too small to affect the density profiles of WDM halos. Down to the resolution attained in our simulations, there is not any significant difference in the density profiles and concentrations of halos obtained in simulations with and without the inclusion of v_w. The density profiles of halos with masses down to ~0.01 the filtering mass M_f can be described by the NFW shape; significant soft cores are not formed. Nevertheless, the concentrations of these halos are lower than those of the CDM counterparts and are approximately independent of mass. The lower concentrations of WDM halos with respect to their CDM counterparts can be accounted for their late formation epoch. Overall, our results point to a series of advantages of a WDM model over the CDM one. In addition to solving the substructure problem, a WDM model with R_f~0.16 Mpc (m_w~0.75 kev; flat cosmology with Omega_L=h=0.7) also predicts concentrations, a Tully-Fisher relation, and formation epochs for small halos which seems to be in better agreement with observations, relative to CDM predictions.

Using a cosmological galactic evolutionary approach to model the Milky Way, we calculate the star formation history (SFH) of the solar neighborhood. The good agreement we obtain with the observational inferences suggests that our physical model describes accurately the long term/large spatial trends of the local and global Milky Way SFH. In this model, star formation is triggered by disk gravitational instabilities and self-regulated by an energy balance in the ISM. The drivers of the SFH are the cosmological gas infall rate and the gas surface density determined by the primordial spin parameter. A LambdaCDM cosmology was used throughout.

We investigate the density profiles and evolution of weakly self-interacting cold dark matter haloes using a numerical code based on the collisional Boltzmann equation. This approach is alternative to N-body techniques in following the dynamical evolution of haloes in the cosmological context and taking into account particle self-interaction. The physical case with a cross section inversely proportional to the relative velocity of the colliding particles is modelled with an unprecedented resolution, spanning five orders of magnitude on the radius for each halo. The modelled haloes cover a mass range from dwarf galaxies to galaxy clusters. We find that for \sigma v_{100} \approx 10^{-24} cm^2/GeV, where \sigma is the cross section per unit mass and v_{100} is the collision velocity in units of 100 km/s, soft cores in good agreement with observations on galactic as well as on galaxy cluster scales are obtained. Remarkably, the observed nearly invariance of the halo central density with mass is reproduced.

We highlight some results from disk galaxy evolution models conceived within a cosmological context. When disk mergers and strong disk-halo feedback are omitted, several properties and correlations of disk galaxies seem to be related to initial conditions given by the CDM model, for example, the intensive galaxy properties, the disk Hubble sequence, and the Tully-Fisher relation in the infrared bands.

One of the predictions of the standard CDM is that dark haloes have centrally divergent density profiles. An extensive body of rotation curve observations of dwarf and low surface brightness galaxies shows the dark haloes of those systems to be characterized by soft constant density central cores. Several physical processes have been proposed to produce soft cores in dark haloes, each one with different scaling properties. With the aim of discriminating among them we have examined the rotation curves of dark matter dominated dwarf and low surface brightness galaxies and the inner mass profiles of two clusters of galaxies lacking a central cD galaxy and with evidence of soft cores in the centre. The core radii and central densities of these haloes scale in a well defined manner with the depth of their potential wells, as measured through the maximum circular velocity. As a result of our analysis we identify self-interacting CDM as a viable solution to the core problem, where a non-singular isothermal core is formed in the halo center surrounded by a Navarro, Frenk, & White profile in the outer parts. We show that this particular physical situation predicts core radii in agreement with observations. Furthermore, using the observed scalings, we derive an expression for the minimum cross section (\sigma) which has an explicit dependence with the halo dispersion velocity (v). If m_x is the mass of the dark matter particle: \sigma/m_x ~4 10^-25 (v/100 km s^-1)^-1 cm^2/Gev.

Using observational data for a sample of dark matter dominated galaxies and two cluster of galaxies, we have found that the central halo density does not depend on its mass, and the core radius is roughly proportional to the maximum rotation velocity. A good agreement with these scaling laws is obtained for CDM halos whose dense inner parts were expanded by gravothermal instabilities if the particles efficiently self-interact only in these parts. We find that the particle cross-section is inversely proportional to the velocity dispersion.

We study the evolution of disk galaxies using galaxy evolutionary models with initial and boundary conditions linked to the hierarchical formation scenario. Disks galaxies are modeled locally within growing cold dark matter halos and including a physical model for star formation. We focus our attention on predictions of the star formation history, size and surface brightness evolution and the evolution of the H- and B-band Tully-Fisher relations. Comparisons with available observational data are presented.

Within the framework of the cold dark matter (CDM) cosmogony, a central cusp in the density profiles of virialized dark haloes is predicted. This prediction disagrees with the soft inner halo mass distribution inferred from observations of dwarf and low surface brightness galaxies, and some clusters of galaxies. By analysing data for some of these objects, we find that the halo central density is nearly independent of the mass from galactic to galaxy cluster scales with an average value of around 0.02 M_sun/pc^3. We show that soft cores can be produced in the CDM haloes by introducing a lower cut-off in the power spectra of fluctuations and assuming high orbital thermal energies during halo formation. However, the scale invariance of the halo central density is not reproduced in these cases. The introduction of self-interaction in the CDM particles offers the most attractive alternative to the core problem. We propose gravothermal expansion as a possible mechanism to produce soft cores in the CDM haloes with self-interacting particles. A global thermodynamical equilibrium can explain the central density scale invariance. We find a minimum cross section capable of establishing isothermal cores in agreement with the observed shallow cores. If \sigma and m_x are the cross section and mass of the dark matter particle, and v is the halo velocity dispersion, then (\sigma/m_x) \~ 4 10^{-25} (100 km s^{-1}/v) cm^2/GeV.

We propose and show that the dependence of the internal face-on dust absorption upon B-band luminosity reported by Wang & Heckman (1996) for a sample of late-type galaxies may explain the empirical color-magnitude and color Tully-Fisher relations of disk galaxies. Thus, in order to explain these relations, it is not necessary to evoke star formation and gas infall efficiencies dependent on the mass of the galaxy system. After applying the Wang & Heckman's luminosity-dependent extinction to hierarchical inside-out galaxy formation and evolution models where the star formation and gas infall efficiencies do not significantly depend on mass, we succesfully predict the observed Tully-Fisher relations in the H and B bands as well as the color-magnitude relation.

We present results of extensive model calculations of disk galaxy evolution within an hierarchical inside-out formation scenario. We first compare properties of the dark halos identified in a cosmological N-body simulation with predictions of a seminumerical method based on an extended collapse model and find a good agreement. We also study the properties of the halos in dependence on their environment. We then describe detailed modelling of the formation and evolution of disks within the growing isolated cold dark matter halos and predictions for the main properties, correlations and evolutionary features of normal disk galaxies. The possible shortcomings of the scenario are discussed.

We used galaxy evolutionary models in a hierarchical inside-out formation scenario to study the origin of the main properties and correlations of disk galaxies. We found that most of these properties and correlations are the result of three (cosmological) initial factors and their dispersions: the virial mass, the halo mass aggregation history (MAH), and the angular momentum given through the spin parameter \lambda. The MAH determines mainly the halo structure and the color indexes while \lambda determines mainly the surface brightness and the bulge-to-disk ratio. We calculated star formation (SF) using a gravitational instability criterion and a self-regulation mechanism in the turbulent ISM. The efficiency of SF in this model is almost independent from the mass. We show that the luminosity-dependent dust absorption empirically determined by Wang & Heckman explains the observed color-magnitude and color Tully-Fisher (TF) relations without the necessity of introducing a mass-dependent SF efficiency. The disks in centrifugal equilibrium form within growing CDM halos with a gas accretion rate proportional to the MAH. The disks present exponential surface density and brightness profiles, negative radial color index gradients, and nearly flat rotation curves. We also calculated the secular formation of a bulge due to gravitational instabilities in the stellar disk. The intensive properties of our models agree with the observational data and the trends of the Hubble sequence are reproduced. The predicted infrared TF and luminosity-radius relations also agree with observations. The main shortcomings of our inside-out hierarchical models are the excessive radial color gradients and the dark halo dominion in the rotation curve decompositions.

We predict the structure and dynamics of disc galaxies using galaxy evolution models within a hierarchical formation scenario The halo mass aggregation histories, for a Lambda CDM model, were generated and used to calculate the virialization of dark matter (DM) haloes. A diversity of halo density profiles were obtained, the most typical one being close to the NFW profile. We modeled the formation of discs in centrifugal equilibrium within the evolving DM haloes using gas accretion rates proportional to the halo mass aggregation rates, and assuming detailed angular momentum conservation. We calculated the gravitational interactions between halo and disc, and the hydrodynamics, star formation, and evolution of the galaxy discs. We found that the slope and zero-point of the infrared Tully-Fisher relations (TFR) may be explained as a direct consequence of the cosmological initial conditions. This relation is almost independent of the assumed disc mass fraction. The rms scatter of the TFR originates mainly from the scatter in the DM halo structure and, to a minor extension, from the dispersion of the primordial spin parameter. The scatter obtained does not disagree with the observational estimates. Our models allow us to understand why the residuals of the TFR do not correlate significantly with disc size or surface brightness (SB), and why low and high SB galaxies have the same TFR. The correlations between gas fraction and SB, and between scale length and V_max agree with those observed. The discs present nearly exponential SB distributions. The shape of the rotation curves changes with the SB and is nearly flat for most cases. The rotation curve decompositions show a dominance of DM down to very small radii. The introduction of shallow cores in the DM halo attenuates this difficulty.

Using galaxy evolutionary models in a hierarchical formation scenario, we predict the structure, dynamics and evolution of disk galaxies in a LCDM universe. Our models include star formation and hydrodynamics of the ISM. We find that the Tully-Fisher relation (TFR) in the I and H bands is an imprint of the mass-velocity relation of the cosmological dark halos. The scatter of the TFR originates mainly from the scatter in the dark halo structure and, to a minor extension, from the dispersion of the primordial spin parameter lambda. Our models allow us to explain why low and high surface brightness galaxies have the same TFR. The disk gas fractions predicted agree with the observations. The disks formed within the growing halos have nearly exponential surface brightness and flat rotation curves. Towards high redshifts, the zero-point of the TFR in the H band increases while in the B-band it slightly decreases.

(Abridged) We study the outer density profiles of dark matter haloes predicted by a generalized secondary infall model and observed in a N-body cosmological simulation of a \Lambda CDM model. We find substantial systematic variations in shapes and concentrations of the halo profiles as well as a strong correlation of the profiles with the environment. In the N-body simulation, the average outer slope of the density profiles, \beta (\rho\propto r^{-\beta}), of isolated haloes is \approx 2.9; 68% of these haloes have values of \beta between 2.5 and 3.8. Haloes in dense environments of clusters are more concentrated and exhibit a broad distribution of \beta with values larger than for isolated haloes . Contrary to what one may expect, the haloes contained within groups and galaxy systems are less concentrated and have flatter outer density profiles than the isolated haloes. The concentration decreases with M_h, but its scatter for a given mass is substantial. The mass and circular velocity of the haloes are strongly correlated: M_h \propto V_m^{\alpha} with \alpha ~ 3.3 (isolated) and ~3.5 (haloes in clusters). For M_h=10^12M_sun the rms deviations from these relations are \Delta logM_h=0.12 and 0.18, respectively. Approximately 30% of the haloes are contained within larger haloes or have massive companions (larger than ~0.3 the mass of the current halo) within 3 virial radii. The remaining 70% of the haloes are isolated objects. The distribution of \beta as well as the concentration-mass and M_h-V_m relations for the isolated haloes agree very well with the predictions of our seminumerical approach which is based on a generalization of the secondary infall model and on the extended Press-Schechter formalism.

Making use of a seminumerical method we develop a scenario of disk galaxy formation and evolution in the framework of inflationary cold dark matter (CDM) cosmologies. Within the virializing dark matter halos, disks in centrifugal equilibrium are built-up and their galactic evolution is followed through an approach which considers the gravitational interactions among the galaxy components, the turbulence and energy balance of the ISM, the star formation (SF) process due to disk gravitational instabilities, the stellar evolution and the secular formation of a bulge. We find that the main properties and correlations of disk galaxies are determined by the mass, the hierarchical mass aggregation history and the primordial angular momentum. The models follow the same trends across the Hubble sequence than the observed galaxies. The predicted TF relation is in good agreement with the observations except for the standart CDM. While the slope of this relation remains almost constant up to intermediate redshifts, its zero-point decreases in the H-band and slightly increases in the B-band. A maximum in the SF rate for most of the models is attained at $z\sim 1.5-2.5$.

We study the evolution of disk galaxies within the frame of the cold dark matter (CDM) cosmologies. The hydrodynamics of a centrifugally supported gaseous disk and the growth of a stellar disk are calculated in detail taking into account the energy balance of the ISM and the gravitational instabilities that concern gas and stars. The halo density profile is derived from the primordial cosmological conditions and its gravitational contraction produced by the disk is included. Several features of the spiral galaxies at different redshifts are predicted, and the main factors which influence on these features are found. A strong evidence is provided that the Tully-Fisher relation is an imprint of the primordial cosmological conditions.

We study the formation of hot spheroidal systems within the frame of a scenario where galaxy formation and evolution is related to the gentle mass aggregation history and primordial angular momentum of protogalaxies, both defined by the cosmological initial conditions. We explore two cases: (1) the hot spheroidal system forms from the dynamical instabilities of the stellar disks, and (2) the spheroidal systems are formed during the dissipative collapse of the gas before falling to the disk in centrifugal equilibrium. In the former case a good agreement with observations for late type galaxies is found. In the second case, contrary to recent claims, we find that the tidal stability criterion is not easily reached. The gas that dissipatively collapses within the dark matter halos should be very clumpy, and the clumps very dense, in order to avoid the tidal destruction of the star formation unities.

We use a semianalytical approach and a CDM cosmological model to study the gravitational collapse and virialization, the structure, as well as the global and statistical properties of isolated dark matter galactic halos which emerge from primordial Gaussian fluctuations. From the statistical properties of the density fluctuation field the possible mass aggregation histories (MAHs) are generated, and these histories are used as the initial conditions of the gravitational collapse. To calculate the structure of the virialized systems we have generalized the secondary infall model. Although a range of halo structures is obtained, the average density profiles agree with the profile derived as a fitting formula to results of N-body cosmological simulations. The comparison of the density profiles with the observational data is disscused, and some possible solutions to the disagreement found in the inner regions are discussed. The results of our approach, after considering the gravitational dragging of the baryon matter that forms a central disk in centrifugal equilibrium, show that the Tully-Fisher relation and its scatter can be explained through the initial cosmological conditions. The sigma8=1 SCDM model produces galaxies with high velocities when compared to observations, but when the SCDM power spectrum is normalized to sigma8=0.57 an excellent agreement with the observable TF relation is found, suggesting that this relation is the natural extension to galactic scales of the observed galaxy distribution power spectrum. The theoretical TF scatter is close to the measured one. The slope of the TF relation is practically invariant with respect to the spin parameter lambda.

We present a scenario for the formation and evolution of disk galaxies within the framework of an inflationary CDM universe, and we compare the results with observations ranking from the present-day up to z~1. Galactic disks are built-up inside-out by gas infall with accretion rates driven by the cosmological mass aggregation history (MAH). We generate the MAHs for a Gaussian density fluctuation field, and we calculate the gravitational collapse and virialization of the dark halos. Assuming detailed angular momentum conservation, disks in centrifugal equilibrium are built-up within them. The disk galactic evolution is followed through a physically self-consistent approach. The main disk galaxy properties and their correlations are determined by the mass, the MAH, and the spin parameter. The models give exponential disk surface brightness (SB) profiles with realistic parameters (including LSB galaxies), nearly flat rotation curves, and negative gradients in the B-V radial profile. The main trends across the Hubble sequence are reproduced. A maximum in the star formation rate for most of the models is attained at z~ 1.5-2.5 (this rate is approximately 2.5-4.0 times larger than the present one). The scale radii and the bulge-to-total ratios decrease with z, while the central SB increase. The B-band TF relation remains almost the same at different redshifts. Our scenario reveals that the cosmological initial conditions are able to determine the main properties of disk galaxies and predict evolutionary features for the present-day dominant galaxy population that are in agreement with very recent deep field observational studies.

The galactic evolutionary theory is now sufficiently mature to provide information about the galaxy formation phase. From evolutionary models we establish a link between the present features of late type galaxies and the protogalactic density fluctuations. This link is used to estimate the fluctuation power spectrum amplitude at galactic scales. The agreement with the power spectrum derived from galaxy distribution and extrapolated to galactic scales, is satisfactory. Meanwhile, the comparison with respect to the standard cold dark matter power spectrum normalized to COBE satellite meassurements reveals only a marginal agreement. Among the most interesting results, the Tully-Fisher relation and the constant central brightness pointed out by Freeman (1970) appear to be the natural consequences of the initial cosmological conditions.

The spherical gravitational collapse and virialization of arbitrary density fluctuations in an expanding universe is studied. In the context of the standard cosmological model and the peak $ansatz$, disk galaxies are supposed to be the product of "extended" gravitational collapses of regions surrounding local-maxima of the density fluctuation field. The substructure over the main profile (merging) is considered as a second-order phenomenon. The range and distribution of possible mass growth histories or accretion regimes for a given present-day mass are estimated and used as the initial condition for galaxy formation. Considering the effect of dark halo contraction due to disk formation, the final rotation curves are calculated. If the structural properties of dark halos are important in establishing the Hubble sequence's properties of visible galaxies, then this sequence is defined not only by the mass, but among other possible parameters, by the accretion regime.

The formation and evolution of disk galaxies in the cosmological context is studied. We consider the observable properties of disk galaxies and treat the disk formation and galactic evolutionary processes in a self-consistent fashion. We find the matter accretion regime to be the dominant ingredient in establishing the Hubble sequence. The accretion regime is a phenomenon directly related to the statistical properties of the primordial density fluctuations from which disk galaxies emerged. In this paper we assumed a Gaussian density field and a standard CDM model.