Fossil groups in the Millennium Simulation Evolution of the brightest galaxies

Aims. We create a catalogue of simulated fossil groups and study their properties, in particular the merging histories of their first-ranked galaxies. We compare the simulated fossil group properties with those of both simulated non-fossil and observed fossil groups. Methods. Using simulations and a mock galaxy catalogue, we searched for massive (>5 x 10(13) h(-1) M-circle dot) fossil groups in the Millennium Simulation Galaxy Catalogue. In addition, we attempted to identify observed fossil groups in the Sloan Digital Sky Survey Data Release 6 using identical selection criteria. Results. Our predictions on the basis of the simulation data are: (a) fossil groups comprise about 5.5% of the total population of groups/clusters with masses larger than 5 x 10(13) h(-1) M-circle dot. This fraction is consistent with the fraction of fossil groups identified in the SDSS, after all observational biases have been taken into account; (b) about 88% of the dominant central objects in fossil groups are elliptical galaxies that have a median R-band absolute magnitude of similar to-23.5-5 log h, which is typical of the observed fossil groups known in the literature; (c) first-ranked galaxies of systems with M > 5 x 10(13) h(-1) M-circle dot, regardless of whether they are either fossil or non-fossil, are mainly formed by gas-poor mergers; (d) although fossil groups, in general, assembled most of their virial masses at higher redshifts in comparison with non-fossil groups, first-ranked galaxies in fossil groups merged later, i.e. at lower redshifts, compared with their non-fossil-group counterparts. Conclusions. We therefore expect to observe a number of luminous galaxies in the centres of fossil groups that show signs of a recent major merger.

The influence of halo assembly on galaxies and galaxy groups

In this paper, we study the variations of groups (galaxy properties according to the assembly history in Sloan Digital Sky Survey Data Release 6 (SDSS-DR6) selected groups. Using mock SDSS group catalogues, we find two suitable indicators of group formation time: (i) the isolation of the group, defined as the distance to the nearest neighbour ill terms of its virial radius and 00 the concentration. measured as the groups inner density calculated using the fifth nearest bright galaxy to the groups centre. Groups Within narrow ranges of Mass ill the mock catalO-Lie show increasing Ifl-OLIP alle With isolation and concentration. However, in the observational data the stellar age, as indicated by the spectral type, only shows a correlation with concentration. We study groups of similar mass and different assembly history. finding important differences ill their galaxy population. Particularly, ill high-mass SDSS groups. the number of members. mass-to-light ratios, red galaxy fractions and the magnitude difference between the brightest and second-brightest group galaxies, show different trends as a function of isolation and concentration, even when it is expected that the latter two quantities correlate with group age. Conversely. low-mass SDSS groups appear to be less sensitive to their assembly history. The correlations detected in the SDSS are not consistent with the trends measured in the mock catalogues. However, discrepancies can he explained in terms of the disagreement found in the a-e-isolation trends, suggesting that the model might be overestimating the effects of, environment, We discuss how the modelling of the cold gas ill satellite galaxies could be responsible for this problem. These results call be Used to improve our Understanding of the evolution of galaxies ill high-density environments.

The dynamical state of galaxy groups and their luminosity content

We analyse the dependence of the luminosity function (LF) of galaxies in groups on group dynamical state. We use the Gaussianity of the velocity distribution of galaxy members as a measurement of the dynamical equilibrium of groups identified in the Sloan Digital Sky Survey Data Release 7 by Zandivarez & Martinez. We apply the Anderson-Darling goodness-of-fit test to distinguish between groups according to whether they have Gaussian or non-Gaussian velocity distributions, i.e. whether they are relaxed or not. For these two subsamples, we compute the (0.1)r-band LF as a function of group virial mass and group total luminosity. For massive groups, , we find statistically significant differences between the LF of the two subsamples: the LFs of groups that have Gaussian velocity distributions have a brighter characteristic absolute magnitude (similar to 0.3 mag) and a steeper faint-end slope (similar to 0.25). We detect a similar effect when comparing the LF of bright [M-0.1r(group) - 5log(h) < -23.5] Gaussian and non-Gaussian groups. Our results indicate that, for massive/luminous groups, the dynamical state of the system is directly related to the luminosity of its galaxy members.

Modeling the cosmological co-evolution of supermassive black holes and galaxies

In this Thesis, we investigate the cosmological co-evolution of supermassive black holes (BHs), Active Galactic Nuclei (AGN) and their hosting dark matter (DM) halos and galaxies, within the standard CDM scenario. We analyze both analytic, semi-analytic and hybrid techniques and use the most recent observational data available to constrain the assumptions underlying our models. First, we focus on very simple analytic models where the assembly of BHs is directly related to the merger history of DM haloes. For this purpose, we implement the two original analytic models of Wyithe & Loeb 2002 and Wyithe & Loeb 2003, compare their predictions to the AGN luminosity function and clustering data, and discuss possible modifications to the models that improve the match to the observation. Then we study more sophisticated semi-analytic models in which however the baryonic physics is neglected as well. Finally we improve the hybrid simulation of De Lucia & Blaizot 2007, adding new semi-analytical prescriptions to describe the BH mass accretion rate during each merger event and its conversion into radiation, and compare the derived BH scaling relations, fundamental plane and mass function, and the AGN luminosity function with observations. All our results support the following scenario: • The cosmological co-evolution of BHs, AGN and galaxies can be well described within the CDM model. • At redshifts z & 1, the evolution history of DM halo fully determines the overall properties of the BH and AGN populations. The AGN emission is triggered mainly by DM halo major mergers and, on average, AGN shine at their Eddington luminosity. • At redshifts z . 1, BH growth decouples from halo growth. Galaxy major mergers cannot constitute the only trigger to accretion episodes in this phase. • When a static hot halo has formed around a galaxy, a fraction of the hot gas continuously accretes onto the central BH, causing a low-energy “radio” activity at the galactic centre, which prevents significant gas cooling and thus limiting the mass of the central galaxies and quenching the star formation at late time. • The cold gas fraction accreted by BHs at high redshifts seems to be larger than at low redshifts.

X-ray properties and evolution of high-redshift AGN, and the gas content of host galaxies

In this thesis, I have investigated the evolution of the high-redshift (z > 3) AGN population by collecting data from some of the major Chandra and XMM-Newton surveys. The final sample (141 sources) is one of the largest selected at z> 3 in the X- rays and it is characterised by a very high redshift completeness (98%). I derived the spectral slopes and obscurations through a spectral anaysis and I assessed the high-z evolution by deriving the luminosity function and the number counts of the sample. The best representation of the AGN evolution is a pure density evolution (PDE) model: the AGN space density is found to decrease by a factor of 10 from z=3 to z=5. I also found that about 50% of AGN are obscured by large column densities (logNH > 23). By comparing these data with those in the Local Universe, I found a positive evolution of the obscured AGN fraction with redshift, especially for luminous (logLx > 44) AGN. I also studied the gas content of z < 1 AGN-hosting galaxies and compared it with that of inactive galaxies. For the first time, I applied to AGN a method to derive the gas mass previously used for inactive galaxies only. AGN are found to live preferentially in gas-rich galaxies. This result on the one hand can help us in understanding the AGN triggering mechanisms, on the other hand explains why AGN are preferentially hosted by star-forming galaxies.

Tracing AGN accretion and star formation in Herschel galaxies

In this Thesis, we study the physical properties and the cosmic evolution of AGN and their host galaxies since z∼3. Our analysis exploits samples of star forming galaxies detected with Herschel at far-IR wavelengths (from 70 up to 500 micron) in different extragalactic surveys, such as COSMOS and the deep GOODS (South and North) fields. The broad-band ancillary data available in COSMOS and the GOODS fields, allows us to implement Herschel and Spitzer photometry with multi-wavelength ancillary data. We perform a multicomponent SED-fitting decomposition to decouple the emission due to star formation from that due to AGN accretion, and to estimate both host-galaxy parameters (such as stellar mass, M* and star formation rate, SFR), and nuclear intrinsic bolometric luminosities. We use the individual estimates of AGN bolometric luminosity obtained through SED-fitting decomposition to reconstruct the redshit evolution of the AGN bolometric luminosity function since z∼3. The resulting trends are used to estimate the overall AGN accretion rate density at different cosmic epochs and to trace the first ever estimate of the AGN accretion history from an IR survey. Later on, we focus our study on the connection between AGN accretion and integrated galaxy properties. We analyse the relationships of AGN accretion with galaxy properties in the SFR-M* plane and at different cosmic epochs. Finally, we infer what is the parameter that best correlates with AGN accretion, comparing our results with previous studies and discussing their physical implications in the context of current scenarios of AGN/galaxy evolution.

Fossil groups in the Millennium simulation Their environment and its evolution

Context. Fossil systems are defined to be X- ray bright galaxy groups ( or clusters) with a two- magnitude difference between their two brightest galaxies within half the projected virial radius, and represent an interesting extreme of the population of galaxy agglomerations. However, the physical conditions and processes leading to their formation are still poorly constrained. Aims. We compare the outskirts of fossil systems with that of normal groups to understand whether environmental conditions play a significant role in their formation. We study the groups of galaxies in both, numerical simulations and observations. Methods. We use a variety of statistical tools including the spatial cross- correlation function and the local density parameter Delta(5) to probe differences in the density and structure of the environments of "" normal"" and "" fossil"" systems in the Millennium simulation. Results. We find that the number density of galaxies surrounding fossil systems evolves from greater than that observed around normal systems at z = 0.69, to lower than the normal systems by z = 0. Both fossil and normal systems exhibit an increment in their otherwise radially declining local density measure (Delta(5)) at distances of order 2.5 r(vir) from the system centre. We show that this increment is more noticeable for fossil systems than normal systems and demonstrate that this difference is linked to the earlier formation epoch of fossil groups. Despite the importance of the assembly time, we show that the environment is different for fossil and non- fossil systems with similar masses and formation times along their evolution. We also confirm that the physical characteristics identified in the Millennium simulation can also be detected in SDSS observations. Conclusions. Our results confirm the commonly held belief that fossil systems assembled earlier than normal systems but also show that the surroundings of fossil groups could be responsible for the formation of their large magnitude gap.

Large change in the predicted number of small halos due to a small amplitude oscillating inflaton potential

A smooth inflaton potential is generally assumed when calculating the primordial power spectrum, implicitly assuming that a very small oscillation in the inflaton potential creates a negligible change in the predicted halo mass function. We show that this is not true. We find that a small oscillating perturbation in the inflaton potential in the slow-roll regime can alter significantly the predicted number of small halos. A class of models derived from supergravity theories gives rise to inflaton potentials with a large number of steps and many trans-Planckian effects may generate oscillations in the primordial power spectrum. The potentials we study are the simple quadratic (chaotic inflation) potential with superimposed small oscillations for small field values. Without leaving the slow-roll regime, we find that for a wide choice of parameters, the predicted number of halos change appreciably. For the oscillations beginning in the 10(7)-10(8) M(circle dot) range, for example, we find that only a 5% change in the amplitude of the chaotic potential causes a 50% suppression of the number of halos for masses between 10(7)-10(8) M(circle dot) and an increase in the number of halos for masses <10(6) M(circle dot) by factors similar to 15-50. We suggest that this might be a solution to the problem of the lack of observed dwarf galaxies in the range 10(7)-10(8) M(circle dot). This might also be a solution to the reionization problem where a very large number of Population III stars in low mass halos are required.

REVISITING THE FOSSIL GROUP CANDIDATES UGC 842 AND NGC 6034

We present a new insight on NGC 6034 and UGC 842, two groups of galaxies previously reported in the literature as being fossil groups. The study is based on optical photometry and spectroscopy obtained with the CTIO Blanco telescope and Sloan Digital Sky Survey archival data. We find that NGC 6034 is embedded in a large structure, dominated by three rich clusters and other small groups. Its first and next four ranked galaxies have magnitude differences in the r band and projected distances which violate the optical criteria to classify it as a fossil group. We confirm that the UGC 842 group is a fossil group, but with about half the velocity dispersion that is reported in previous works. The velocity distribution of its galaxies reveals the existence of two structures in its line of sight, one with sigma(nu) similar to 223 km s(-1) and another with sigma(nu) similar to 235 km s(-1), with a difference in velocity of similar to 820 km s(-1). The main structure is dominated by passive galaxies, while these represent similar to 60% of the second structure. The X-ray temperature for the intragroup medium of a group with such a velocity dispersion is expected to be kT similar to 0.5-1 keV, against the observed value of kT similar to 1.9 keV reported in the literature. This result makes UGC 842 a special case among fossil groups because (1) it represents more likely the interaction between two small groups, which warms the intragroup medium and/or (2) it could constitute evidence that member galaxies lost energy in the process of spiraling toward the group center, and decreased the velocity dispersion of the system. As far as we know, UGC 842 is the first low-mass fossil group studied in detail.

The 1000 brightest HIPASS galaxies: Newly cataloged galaxies

The H I Parkes All-Sky Survey (HIPASS) is a blind 21 cm survey for extragalactic neutral hydrogen, covering the whole southern sky. The HIPASS Bright Galaxy Catalog (BGC) is a subset of HIPASS and contains the 1000 H I brightest (peak flux density) galaxies. Here we present the 138 HIPASS BGC galaxies that had no redshift measured prior to the Parkes multibeam H I surveys. Of the 138 galaxies, 87 are newly cataloged. Newly cataloged is defined as having no optical ( or infrared) counterpart in the NASA/IPAC Extragalactic Database. Using the Digitized Sky Survey, we identify optical counterparts for almost half of the newly cataloged galaxies, which are typically of irregular or Magellanic morphological type. Several H I sources appear to be associated with compact groups or pairs of galaxies rather than an individual galaxy. The majority ( 57) of the newly cataloged galaxies lie within 10degrees of the Galactic plane and are missing from optical surveys as a result of confusion with stars or dust extinction. This sample also includes newly cataloged galaxies first discovered by Henning et al. in the H I shallow survey of the zone of avoidance. The other 30 newly cataloged galaxies escaped detection because of their low surface brightness or optical compactness. Only one of these, HIPASS J0546-68, has no obvious optical counterpart, as it is obscured by the Large Magellanic Cloud. We find that the newly cataloged galaxies with -b->10degrees are generally lower in H I mass and narrower in velocity width compared with the total HIPASS BGC. In contrast, newly cataloged galaxies behind the Milky Way are found to be statistically similar to the entire HIPASS BGC. In addition to these galaxies, the HIPASS BGC contains four previously unknown H I clouds.

Gaia photometry for white dwarfs

Context. White dwarfs can be used to study the structure and evolution of the Galaxy by analysing their luminosity function and initial mass function. Among them, the very cool white dwarfs provide the information for the early ages of each population. Because white dwarfs are intrinsically faint only the nearby (~ 20 pc) sample is reasonably complete. The Gaia space mission will drastically increase the sample of known white dwarfs through its 5-6 years survey of the whole sky up to magnitude V = 20-25. Aims. We provide a characterisation of Gaia photometry for white dwarfs to better prepare for the analysis of the scientific output of the mission. Transformations between some of the most common photometric systems and Gaia passbands are derived. We also give estimates of the number of white dwarfs of the different galactic populations that will be observed. Methods. Using synthetic spectral energy distributions and the most recent Gaia transmission curves, we computed colours of three different types of white dwarfs (pure hydrogen, pure helium, and mixed composition with H/He = 0.1). With these colours we derived transformations to other common photometric systems (Johnson-Cousins, Sloan Digital Sky Survey, and 2MASS). We also present numbers of white dwarfs predicted to be observed by Gaia. Results. We provide relationships and colourcolour diagrams among different photometric systems to allow the prediction and/or study of the Gaia white dwarf colours. We also include estimates of the number of sources expected in every galactic population and with a maximum parallax error. Gaia will increase the sample of known white dwarfs tenfold to about 200 000. Gaia will be able to observe thousands of very cool white dwarfs for the first time, which will greatly improve our understanding of these stars and early phases of star formation in our Galaxy.

The Besançon Galaxy Model renewed. I. Constraints on the Galactic thin disc evolution from Tycho data

Context. The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allow to test hypotheses on the star formation history, star evolution, as well as chemical and dynamical evolution of the Galaxy. Aims. The new version of the Besanc¸on Galaxy Model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the Initial Mass Function (IMF) and Star Formation Rate (SFR) of the Galactic disc. Methods. We present a new strategy for the generation of thin disc stars which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. We keep in this new scheme the local dynamical self-consistency as in Bienayme et al (1987). We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. Results. For the first time, a full sky comparison is performed between BGM and data. This strategy allows to constrain the IMF slope at high masses which is found to be close to 3.0, excluding a shallower slope such as Salpeter"s one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 M pc−3 implying that there is no compelling evidence for significant amount of dark matter in the disc. While the model is fitted to Tycho2 data, a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. Conclusions. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

Gaia photometry for white dwarfs

Context. White dwarfs can be used to study the structure and evolution of the Galaxy by analysing their luminosity function and initial mass function. Among them, the very cool white dwarfs provide the information for the early ages of each population. Because white dwarfs are intrinsically faint only the nearby (~ 20 pc) sample is reasonably complete. The Gaia space mission will drastically increase the sample of known white dwarfs through its 5-6 years survey of the whole sky up to magnitude V = 20-25. Aims. We provide a characterisation of Gaia photometry for white dwarfs to better prepare for the analysis of the scientific output of the mission. Transformations between some of the most common photometric systems and Gaia passbands are derived. We also give estimates of the number of white dwarfs of the different galactic populations that will be observed. Methods. Using synthetic spectral energy distributions and the most recent Gaia transmission curves, we computed colours of three different types of white dwarfs (pure hydrogen, pure helium, and mixed composition with H/He = 0.1). With these colours we derived transformations to other common photometric systems (Johnson-Cousins, Sloan Digital Sky Survey, and 2MASS). We also present numbers of white dwarfs predicted to be observed by Gaia. Results. We provide relationships and colour-colour diagrams among different photometric systems to allow the prediction and/or study of the Gaia white dwarf colours. We also include estimates of the number of sources expected in every galactic population and with a maximum parallax error. Gaia will increase the sample of known white dwarfs tenfold to about 200 000. Gaia will be able to observe thousands of very cool white dwarfs for the first time, which will greatly improve our understanding of these stars and early phases of star formation in our Galaxy.

The Besançon Galaxy Model renewed. I. Constraints on the Galactic thin disc evolution from Tycho data

Context. The understanding of Galaxy evolution can be facilitated by the use of population synthesis models, which allow to test hypotheses on the star formation history, star evolution, as well as chemical and dynamical evolution of the Galaxy. Aims. The new version of the Besanc¸on Galaxy Model (hereafter BGM) aims to provide a more flexible and powerful tool to investigate the Initial Mass Function (IMF) and Star Formation Rate (SFR) of the Galactic disc. Methods. We present a new strategy for the generation of thin disc stars which assumes the IMF, SFR and evolutionary tracks as free parameters. We have updated most of the ingredients for the star count production and, for the first time, binary stars are generated in a consistent way. We keep in this new scheme the local dynamical self-consistency as in Bienayme et al (1987). We then compare simulations from the new model with Tycho-2 data and the local luminosity function, as a first test to verify and constrain the new ingredients. The effects of changing thirteen different ingredients of the model are systematically studied. Results. For the first time, a full sky comparison is performed between BGM and data. This strategy allows to constrain the IMF slope at high masses which is found to be close to 3.0, excluding a shallower slope such as Salpeter"s one. The SFR is found decreasing whatever IMF is assumed. The model is compatible with a local dark matter density of 0.011 M pc−3 implying that there is no compelling evidence for significant amount of dark matter in the disc. While the model is fitted to Tycho2 data, a magnitude limited sample with V<11, we check that it is still consistent with fainter stars. Conclusions. The new model constitutes a new basis for further comparisons with large scale surveys and is being prepared to become a powerful tool for the analysis of the Gaia mission data.

The VIMOS-VLT Deep Survey: the evolution of type-1 AGN

Quasars and AGN play an important role in many aspects of the modern cosmology. Of particular interest is the issue of the interplay between AGN activity and formation and evolution of galaxies and structures. Studies on nearby galaxies revealed that most (and possibly all) galaxy nuclei contain a super-massive black hole (SMBH) and that between a third and half of them are showing some evidence of activity (Kormendy and Richstone, 1995). The discovery of a tight relation between black holes mass and velocity dispersion of their host galaxy suggests that the evolution of the growth of SMBH and their host galaxy are linked together. In this context, studying the evolution of AGN, through the luminosity function (LF), is fundamental to constrain the theories of galaxy and SMBH formation and evolution. Recently, many theories have been developed to describe physical processes possibly responsible of a common formation scenario for galaxies and their central black hole (Volonteri et al., 2003; Springel et al., 2005a; Vittorini et al., 2005; Hopkins et al., 2006a) and an increasing number of observations in different bands are focused on collecting larger and larger quasar samples. Many issues remain however not yet fully understood. In the context of the VVDS (VIMOS-VLT Deep Survey), we collected and studied an unbiased sample of spectroscopically selected faint type-1 AGN with a unique and straightforward selection function. Indeed, the VVDS is a large, purely magnitude limited spectroscopic survey of faint objects, free of any morphological and/or color preselection. We studied the statistical properties of this sample and its evolution up to redshift z 4. Because of the contamination of the AGN light by their host galaxies at the faint magnitudes explored by our sample, we observed that a significant fraction of AGN in our sample would be missed by the UV excess and morphological criteria usually adopted for the pre-selection of optical QSO candidates. If not properly taken into account, this failure in selecting particular sub-classes of AGN could, in principle, affect some of the conclusions drawn from samples of AGN based on these selection criteria. The absence of any pre-selection in the VVDS leads us to have a very complete sample of AGN, including also objects with unusual colors and continuum shape. The VVDS AGN sample shows in fact redder colors than those expected by comparing it, for example, with the color track derived from the SDSS composite spectrum. In particular, the faintest objects have on average redder colors than the brightest ones. This can be attributed to both a large fraction of dust-reddened objects and a significant contamination from the host galaxy. We have tested these possibilities by examining the global spectral energy distribution of each object using, in addition to the U, B, V, R and I-band magnitudes, also the UV-Galex and the IR-Spitzer bands, and fitting it with a combination of AGN and galaxy emission, allowing also for the possibility of extinction of the AGN flux. We found that for 44% of our objects the contamination from the host galaxy is not negligible and this fraction decreases to 21% if we restrict the analysis to a bright subsample (M1450 <-22.15). Our estimated integral surface density at IAB < 24.0 is 500 AGN per square degree, which represents the highest surface density of a spectroscopically confirmed sample of optically selected AGN. We derived the luminosity function in B-band for 1.0 < z < 3.6 using the 1/Vmax estimator. Our data, more than one magnitude fainter than previous optical surveys, allow us to constrain the faint part of the luminosity function up to high redshift. A comparison of our data with the 2dF sample at low redshift (1 < z < 2.1) shows that the VDDS data can not be well fitted with the pure luminosity evolution (PLE) models derived by previous optically selected samples. Qualitatively, this appears to be due to the fact that our data suggest the presence of an excess of faint objects at low redshift (1.0 < z < 1.5) with respect to these models. By combining our faint VVDS sample with the large sample of bright AGN extracted from the SDSS DR3 (Richards et al., 2006b) and testing a number of different evolutionary models, we find that the model which better represents the combined luminosity functions, over a wide range of redshift and luminosity, is a luminosity dependent density evolution (LDDE) model, similar to those derived from the major Xsurveys. Such a parameterization allows the redshift of the AGN density peak to change as a function of luminosity, thus fitting the excess of faint AGN that we find at 1.0 < z < 1.5. On the basis of this model we find, for the first time from the analysis of optically selected samples, that the peak of the AGN space density shifts significantly towards lower redshift going to lower luminosity objects. The position of this peak moves from z 2.0 for MB <-26.0 to z 0.65 for -22< MB <-20. This result, already found in a number of X-ray selected samples of AGN, is consistent with a scenario of “AGN cosmic downsizing”, in which the density of more luminous AGN, possibly associated to more massive black holes, peaks earlier in the history of the Universe (i.e. at higher redshift), than that of low luminosity ones, which reaches its maximum later (i.e. at lower redshift). This behavior has since long been claimed to be present in elliptical galaxies and it is not easy to reproduce it in the hierarchical cosmogonic scenario, where more massive Dark Matter Halos (DMH) form on average later by merging of less massive halos.