BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Pierluigi Monaco

Pierluigi Monaco appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GAastro-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

40works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2026arXiv

Accelerating cosmological simulations on GPUs: a step towards sustainability and green-awareness

The increasing complexity and scale of cosmological N-body simulations, driven by astronomical surveys like Euclid, call for a paradigm shift towards more sustainable and energy-efficient high-performance computing (HPC). The rising energy consumption of supercomputing facilities poses a significant environmental and financial challenge. In this work, we build upon a recently developed GPU implementation of pinocchio, a widely-used tool for the fast generation of dark matter (DM) halo catalogues, to investigate energy consumption. Using a different resource configuration, we confirmed the time-to-solution behavior observed in a companion study, and we use these runs to compare time-to-solution with energy-to-solution. By profiling the code on various HPC platforms with a newly developed implementation of the Power Measurement Toolkit (PMT), we demonstrate an 8x reduction in energy-to-solution and 8x speed-up in time-to-solution compared to the CPU-only version. Taken together, these gains translate into an overall efficiency improvement of up to 64x. Our results show that the GPU-accelerated pinocchio not only achieves substantial speed-up, making the generation of large-scale mock catalogues more tractable, but also significantly reduces the energy footprint of the simulations. This work represents an step towards ``green-aware" scientific computing in cosmology, proving that performance and sustainability can be simultaneously achieved.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dust evolution in zoom-in cosmological simulations of galaxy formation

We present cosmological zoom-in hydro-dynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain size distribution, with the MUPPI star formation and feedback sub-resolution model. Our dust evolution model follows carbon and silicate dust separately. To distinguish differences induced by the chaotic behaviour of simulations from those genuinely due to different simulation set-up, we run each model six times, after introducing tiny perturbations in the initial conditions. With this method, we discuss the role of various dust-related physical processes and the effect of a few possible approximations adopted in the literature. Metal depletion and dust cooling affect the evolution of the system, causing substantial variations in its stellar, gas and dust content. We discuss possible effects on the Spectral Energy Distribution of the significant variations of the size distribution and chemical composition of grains, as predicted by our simulations during the evolution of the galaxy. We compare dust surface density, dust-to-gas ratio and small-to-big grain mass ratio as a function of galaxy radius and gas metallicity predicted by our fiducial run with recent observational estimates for three disc galaxies of different masses. The general agreement is good, in particular taking into account that we have not adjusted our model for this purpose.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Testing the Reliability of Fast Methods for Weak Lensing Simulations: WL-MOKA on PINOCCHIO

The generation of simulated convergence maps is of key importance in fully exploiting weak lensing by Large Scale Structure (LSS) from which cosmological parameters can be derived. In this paper we present an extension of the PINOCCHIO code which produces catalogues of dark matter haloes so that it is capable of simulating weak lensing by LSS. Like WL-MOKA, the method starts with a random realisation of cosmological initial conditions, creates a halo catalogue and projects it onto the past-light-cone, and paints in haloes assuming parametric models for the mass density distribution within them. Large scale modes that are not accounted for by the haloes are constructed using linear theory. We discuss the systematic errors affecting the convergence power spectra when Lagrangian Perturbation Theory at increasing order is used to displace the haloes within PINOCCHIO, and how they depend on the grid resolution. Our approximate method is shown to be very fast when compared to full ray-tracing simulations from an N-Body run and able to recover the weak lensing signal, at different redshifts, with a few percent accuracy. It also allows for quickly constructing weak lensing covariance matrices, complementing PINOCCHIO's ability of generating the cluster mass function and galaxy clustering covariances and thus paving the way for calculating cross covariances between the different probes. This work advances these approximate methods as tools for simulating and analysing surveys data for cosmological purposes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Accuracy of Weak Lensing Simulations

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Rise of Active Galactic Nuclei in the GAlaxy Evolution and Assembly semi-analytic model

We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas accretion onto supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted onto the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disc instabilities in star forming discs. We compare our predictions with the luminosity function (LF) observed for Active Galactic Nuclei (AGN) and we confirm that a non-instantaneous accretion timescale (either in the form of a low-angular momentum reservoir or as an assumed light curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies or from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies ($> 10^{11} M_\odot$) without changing their overall assembly history. These winds also affect the SFR of lower mass galaxies, resulting in a too large fraction of passive galaxies at $< 10^{10} M_\odot$. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than $10^8 M_\odot$ are already in a self-regulated state as inferred from observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Toward a robust inference method for the galaxy bispectrum: likelihood function and model selection

The forthcoming generation of galaxy redshift surveys will sample the large-scale structure of the Universe over unprecedented volumes with high-density tracers. This advancement will make robust measurements of three-point clustering statistics possible. In preparation for this improvement, we investigate how several methodological choices can influence inferences based on the bispectrum about galaxy bias and shot noise. We first measure the real-space bispectrum of dark-matter haloes extracted from 298 N-body simulations covering a volume of approximately $1000 h^{-3} \mathrm{Gpc}^3$. We then fit a series of theoretical models based on tree-level perturbation theory to the numerical data. To achieve this, we estimate the covariance matrix of the measurement errors by using 10,000 mock catalogues generated with the Pinocchio code. We study how the model constraints are influenced by the binning strategy for the bispectrum configurations and by the form of the likelihood function. We also use Bayesian model-selection techniques to single out the optimal theoretical description of our data. We find that a three-parameter bias model combined with Poissonian shot noise is necessary to model the halo bispectrum up to scales of $k_\mathrm{max}\lesssim 0.08 h \mathrm{Mpc}^{-1}$, although fitting formulae that relate the bias parameters can be helpful to reduce the freedom of the model without compromising accuracy. Our data clearly disfavour local Eulerian and local Lagrangian bias models and do not require corrections to Poissonian shot noise. We anticipate that model-selection diagnostics will be particularly useful to extend the analysis to smaller scales as, in this case, the number of model parameters will grow significantly large.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fast numerical method to generate halo catalogs in modified gravity (part I): second-order Lagrangian Perturbation Theory

We present and test a new numerical method to determine second-order Lagrangian displacement fields in the context of modified gravity (MG) theories. We start from the extension of Lagrangian Perturbation Theory to a class of MG models that can be described by a parametrized Poisson equation, with the introduction of a scale-dependent function. We exploit fast Fourier transforms to compute the full source term of the differential equation for the second-order Lagrangian displacement field. We compare its mean to the source term computed for specific configurations for which a k-dependent solution can be found numerically. We choose the configuration that best matches the full source term, thus obtaining an approximate factorization of the second-order displacement field as the space term valid for standard gravity times a k-dependent, second-order growth factor $D_2(k,t)$. This approximation is used to compute second order displacements for particles. The method is tested against N-body simulations run with standard and $f(R)$ gravity: we rely on the results of a friends-of-friends code run on the N-body snapshots to assign particles to halos, then compute the halo power spectrum. We find very consistent results for the two gravity theories: second-order LPT (2LPT) allows to recover the halo power spectrum of N-body simulations within $\sim 10\%$ precision to $k\sim 0.2-0.4\ h\ {\rm Mpc}^{-1}$, as well as halo positions, with an error that is a fraction of the inter-particle distance. We show that, when considering the same level of non-linearity in the density field, the performance of 2LPT with MG is the same (within $1\%$) as the one obtained for the standard $Λ$CDM model with General Relativity. When implemented in a computer code, this formulation of 2LPT can quickly generate dark matter distributions with $f(R)$ gravity, and can easily be extended to other MG theories.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Impact of AGN feedback on galaxies and their multiphase ISM across cosmic time

We present simulations of galaxy formation, based on the GADGET-3 code, in which a sub-resolution model for star formation and stellar feedback is interfaced with a new model for AGN feedback. Our sub-resolution model describes a multiphase ISM, accounting for hot and cold gas within the same resolution element: we exploit this feature to investigate the impact of coupling AGN feedback energy to the different phases of the ISM over cosmic time. Our fiducial model considers that AGN feedback energy coupling is driven by the covering factors of the hot and cold phases. We perform a suite of cosmological hydrodynamical simulations of disc galaxies ($M_{\rm halo, \, DM} \simeq 2 \cdot 10^{12}$ M$_{\odot}$, at $z=0$), to investigate: $(i)$ the effect of different ways of coupling AGN feedback energy to the multiphase ISM; $(ii)$ the impact of different prescriptions for gas accretion (i.e. only cold gas, both cold and hot gas, with the additional possibility of limiting gas accretion from cold gas with high angular momentum); $(iii)$ how different models of gas accretion and coupling of AGN feedback energy affect the coevolution of supermassive BHs and their host galaxy. We find that at least a share of the AGN feedback energy has to couple with the diffuse gas, in order to avoid an excessive growth of the BH mass. When the BH only accretes cold gas, it experiences a growth that is faster than in the case in which both cold and hot gas are accreted. If the accretion of cold gas with high angular momentum is reduced, the BH mass growth is delayed, the BH mass at $z=0$ is reduced by up to an order of magnitude, and the BH is prevented from accreting below $z \lesssim 2$, when the galaxy disc forms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Approximate methods for the generation of dark matter halo catalogs in the age of precision cosmology

Precision cosmology has recently triggered new attention on the topic of approximate methods for the clustering of matter on large scales, whose foundations date back to the period from late '60s to early '90s. Indeed, although the prospect of reaching sub-percent accuracy in the measurement of clustering poses a challenge even to full N-body simulations, an accurate estimation of the covariance matrix of clustering statistics, not to mention the sampling of parameter space, requires usage of a large number (hundreds in the most favourable cases) of simulated (mock) galaxy catalogs. Combination of few N-body simulations with a large number of realizations performed with approximate methods gives the most promising approach to solve these problems with a reasonable amount of resources. In this paper I review this topic, starting from the foundations of the methods, then going through the pioneering efforts of the '90s, and finally presenting the latest extensions and a few codes that are now being used in present-generation surveys and thoroughly tested to assess their performance in the context of future surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Kinetic AGN Feedback Effects on Cluster Cool Cores Simulated using SPH

We implement novel numerical models of AGN feedback in the SPH code GADGET-3, where the energy from a supermassive black hole (BH) is coupled to the surrounding gas in the kinetic form. Gas particles lying inside a bi-conical volume around the BH are imparted a one-time velocity (10,000 km/s) increment. We perform hydrodynamical simulations of isolated cluster (total mass 10^14 /h M_sun), which is initially evolved to form a dense cool core, having central T<10^6 K. A BH resides at the cluster center, and ejects energy. The feedback-driven fast wind undergoes shock with the slower-moving gas, which causes the imparted kinetic energy to be thermalized. Bipolar bubble-like outflows form propagating radially outward to a distance of a few 100 kpc. The radial profiles of median gas properties are influenced by BH feedback in the inner regions (r<20-50 kpc). BH kinetic feedback, with a large value of the feedback efficiency, depletes the inner cool gas and reduces the hot gas content, such that the initial cool core of the cluster is heated up within a time 1.9 Gyr, whereby the core median temperature rises to above 10^7 K, and the central entropy flattens. Our implementation of BH thermal feedback (using the same efficiency as kinetic), within the star-formation model, cannot do this heating, where the cool core remains. The inclusion of cold gas accretion in the simulations produces naturally a duty cycle of the AGN with a periodicity of 100 Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Spectral Slope and Escape Fraction of Bright Quasars at $z \sim 3.8$: the Contribution to the Cosmic UV Background

We use a sample of 1669 QSOs ($r<20.15$, $3.6<z<4.0$) from the BOSS survey to study the intrinsic shape of their continuum and the Lyman continuum photon escape fraction (f$_{esc}$), estimated as the ratio between the observed flux and the expected intrinsic flux (corrected for the intergalactic medium absorption) in the wavelength range 865-885 Å rest-frame. Modelling the intrinsic QSO continuum shape with a power-law, $F_λ\proptoλ^{-γ}$, we find a median $γ=1.30$ (with a dispersion of $0.38$, no dependence on the redshift and a mild intrinsic luminosity dependence) and a mean f$_{esc}=0.75$ (independent of the QSO luminosity and/or redshift). The f$_{esc}$ distribution shows a peak around zero and a long tail of higher values, with a resulting dispersion of $0.7$. If we assume for the QSO continuum a double power-law shape (also compatible with the data) with a break located at $λ_{\rm br}=1000$ Å and a softening $Δγ=0.72 $ at wavelengths shorter than $λ_{\rm br}$, the mean f$_{esc}$ rises to $=0.82$. Combining our $γ$ and f$_{esc}$ estimates with the observed evolution of the AGN luminosity function (LF) we compute the AGN contribution to the UV ionizing background (UVB) as a function of redshift. AGN brighter than one tenth of the characteristic luminosity of the LF are able to produce most of it up $z\sim 3$, if the present sample is representative of their properties. At higher redshifts a contribution of the galaxy population is required. Assuming an escape fraction of Lyman continuum photons from galaxies between $5.5$ and $7.6\%$, independent of the galaxy luminosity and/or redshift, a remarkably good fit to the observational UVB data up to $z\sim 6$ is obtained. At lower redshift the extrapolation of our empirical estimate agrees well with recent UVB observations, dispelling the so-called Photon Underproduction Crisis.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interpreting the possible break in the Black Hole - Bulge mass relation

Recent inspections of local available data suggest that the almost linear relation between the stellar mass of spheroids ($M_{\rm sph}$) and the mass of the super massive Black Holes (BHs) residing at their centres, shows a break below $M_{\rm sph} \sim 10^{10}\ {\rm M}_\odot$, with a steeper, about quadratic relation at smaller masses. We investigate the physical mechanisms responsible for the change in slope of this relation, by comparing data with the results of the semi-analytic model of galaxy formation MORGANA, which already predicted such a break in its original formulation. We find that the change of slope is mostly induced by effective stellar feedback in star-forming bulges. The shape of the relation is instead quite insensitive to other physical mechanisms connected to BH accretion such as disc instabilities, galaxy mergers, Active Galactic Nucleus (AGN) feedback, or even the exact modelling of accretion onto the BH, direct or through a reservoir of low angular momentum gas. Our results support a scenario where most stars form in the disc component of galaxies and are carried to bulges through mergers and disc instabilities, while accretion onto BHs is connected to star formation in the spheroidal component. Therefore, a model of stellar feedback that produces stronger outflows in star-forming bulges than in discs will naturally produce a break in the scaling relation. Our results point to a form of co-evolution especially at lower masses, below the putative break, mainly driven by stellar feedback rather than AGN feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy Cosmology: Comparison of Galaxy Formation Models

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A semi-analytic model comparison: testing cooling models against hydrodynamical simulations

We compare predictions of cooled masses and cooling rates from three stripped-down Semi-Analytic Models (SAMs) of galaxy formation with the results of N-body+SPH simulations with gas particle mass of 3.9x10^6 Msun, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ~2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z=0 to within ~20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved halos is subtracted out, SAMs tend to under-predict by ~30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses SAMs tend to over-predict cooling rates, though the numerical result may be affected by the use of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z~2, cold gas accumulates very quickly in the simulated halos. This accumulation is reproduced by SAMs with varying degrees of accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Environmental dependence of bulge-dominated galaxy sizes in hierarchical models of galaxy formation. Comparison with the local Universe

We compare state-of-the-art semi-analytic models of galaxy formation as well as advanced sub-halo abundance matching models with a large sample of early-type galaxies from SDSS at z < 0.3. We focus our attention on the dependence of median sizes of central galaxies on host halo mass. The data do not show any difference in the structural properties of early-type galaxies with environment, at fixed stellar mass. All hierarchical models considered in this work instead tend to predict a moderate to strong environmental dependence, with the median size increasing by a factor of about 1.5-3 when moving from low to high mass host haloes. At face value the discrepancy with the data is highly significant, especially at the cluster scale, for haloes above log Mhalo > 14. The convolution with (correlated) observational errors reduces some of the tension. Despite the observational uncertainties, the data tend to disfavour hierarchical models characterized by a relevant contribution of disc instabilities to the formation of spheroids, strong gas dissipation in (major) mergers, short dynamical friction timescales, and very short quenching timescales in infalling satellites. We also discuss a variety of additional related issues, such as the slope and scatter in the local size-stellar mass relation, the fraction of gas in local early-type galaxies, and the general predictions on satellite galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Galactic outflow and diffuse gas properties at z>=1 using different baryonic feedback models

We measure and quantify properties of galactic outflows and diffuse gas at $z \geq 1$ in cosmological hydrodynamical simulations. Our novel sub-resolution model, MUPPI, implements supernova feedback using fully local gas properties, where the wind velocity and mass loading are not given as input. We find the following trends at $z = 2$ by analysing central galaxies having a stellar mass higher than $10^{9} M_{\odot}$. The outflow velocity and mass outflow rate ($\dot{M}_{\rm out}$) exhibit positive correlations with galaxy mass and with the star formation rate (SFR). However, most of the relations present a large scatter. The outflow mass loading factor ($η$) is between $0.2 - 10$. The comparison Effective model generates a constant outflow velocity, and a negative correlation of $η$ with halo mass. The number fraction of galaxies where outflow is detected decreases at lower redshifts, but remains more than $80 \%$ over $z = 1 - 5$. High SF activity at $z \sim 2 - 4$ drives strong outflows, causing the positive and steep correlations of velocity and $\dot{M}_{\rm out}$ with SFR. The outflow velocity correlation with SFR becomes flatter at $z = 1$, and $η$ displays a negative correlation with halo mass in massive galaxies. Our study demonstrates that both the MUPPI and Effective models produce significant outflows at $\sim 1 / 10$ of the virial radius; at the same time shows that the properties of outflows generated can be different from the input speed and mass loading in the Effective model. Our MUPPI model, using local properties of gas in the sub-resolution recipe, is able to develop galactic outflows whose properties correlate with global galaxy properties, and consistent with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Properties of barred spiral disks in hydrodynamical cosmological simulations

We present a quantification of the properties of bars in two N-body+SPH cosmological simulations of spiral galaxies, named GA and AqC. The initial conditions were obtained using the zoom-in technique and represent two dark matter (DM) halos of $2-3\times10^{12}\ {\rm M}_\odot$, available at two different resolutions. The resulting galaxies are presented in the companion paper of Murante et al. (2014). We find that the GA galaxy has a bar of length $8.8$ kpc, present at the two resolution levels even though with a slightly different strength. Classical bar signatures (e.g. pattern of streaming motions, high $m=2$ Fourier mode with roughly constant phase) are consistently found at both resolutions. Though a close encounter with a merging satellite at $z\sim0.6$ (mass ratio $1:50$) causes a strong, transient spiral pattern and some heating of the disk, we find that bar instability is due to secular process, caused by a low Toomre parameter $Q\lesssim1$ due to accumulation of mass in the disk. The AqC galaxy has a slightly different history: it suffers a similar tidal disturbance due to a merging satellite at $z\sim0.5$ but with a mass ratio of $1:32$, that triggers a bar in the high-resolution simulation, while at low resolution the merging is found to take place at a later time, so that both secular evolution and merging are plausible triggers for bar instability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating realistic disk galaxies with a novel sub-resolution ISM model

We present results of cosmological simulations of disk galaxies carried out with the GADGET-3 TreePM+SPH code, where star formation and stellar feedback are described using our MUlti Phase Particle Integrator (MUPPI) model. This description is based on simple multi-phase model of the interstellar medium at unresolved scales, where mass and energy flows among the components are explicitly followed by solving a system of ordinary differential equations. Thermal energy from SNe is injected into the local hot phase, so as to avoid that it is promptly radiated away. A kinetic feedback prescription generates the massive outflows needed to avoid the over-production of stars. We use two sets of zoomed-in initial conditions of isolated cosmological halos with masses (2-3) * 10^{12} Msun, both available at several resolution levels. In all cases we obtain spiral galaxies with small bulge-over-total stellar mass ratios (B/T \approx 0.2), extended stellar and gas disks, flat rotation curves and realistic values of stellar masses. Gas profiles are relatively flat, molecular gas is found to dominate at the centre of galaxies, with star formation rates following the observed Schmidt-Kennicutt relation. Stars kinematically belonging to the bulge form early, while disk stars show a clear inside-out formation pattern and mostly form after redshift z=2. However, the baryon conversion efficiencies in our simulations differ from the relation given by Moster et al. (2010) at a 3 sigma level, thus indicating that our stellar disks are still too massive for the Dark Matter halo in which they reside. Results are found to be remarkably stable against resolution. This further demonstrates the feasibility of carrying out simulations producing a realistic population of galaxies within representative cosmological volumes, at a relatively modest resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The definition of environment and its relation to the quenching of galaxies at z=1-2 in a hierarchical Universe

A well calibrated method to describe the environment of galaxies at all redshifts is essential for the study of structure formation. Such a calibration should include well understood correlations with halo mass, and the possibility to identify galaxies which dominate their potential well (centrals), and their satellites. Focusing on z = 1 and 2 we propose a method of environmental calibration which can be applied to the next generation of low to medium resolution spectroscopic surveys. Using an up-to-date semi-analytic model of galaxy formation, we measure the local density of galaxies in fixed apertures on different scales. There is a clear correlation of density with halo mass for satellite galaxies, while a significant population of low mass centrals is found at high densities in the neighbourhood of massive haloes. In this case the density simply traces the mass of the most massive halo within the aperture. To identify central and satellite galaxies, we apply an observationally motivated stellar mass rank method which is both highly pure and complete, especially in the more massive haloes where such a division is most meaningful. Finally we examine a test case for the recovery of environmental trends: the passive fraction of galaxies and its dependence on stellar and halo mass for centrals and satellites. With careful calibration, observationally defined quantities do a good job of recovering known trends in the model. This result stands even with reduced redshift accuracy, provided the sample is deep enough to preserve a wide dynamic range of density.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Research Note on the Implementation of Star Formation and Stellar Feedback in Semi-Analytic Models

We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, disc instabilities or AGN feedback. We run these versions on identical samples of dark matter (DM) haloes extracted from high-resolution N-body simulations in order to perform both statistical analysis and object-by-object comparisons. We compare our results with previous work based on stripped-down versions of the same SAMs including only gas cooling, and show that all feedback models provide coherent modifications in the distribution of baryons between the various gas phases. In particular, we find that the predicted hot gas fractions are considerably increased by up to a factor of three, while the corresponding cold gas fractions are correspondingly decreased, and a significant amount of mass is ejected from the DM halo. Nonetheless, we also find relevant differences in the predicted properties of model galaxies among the three SAMs: these deviations are more relevant at mass scales comparable to that of our own Galaxy, and are reduced at larger masses, confirming the varying impact of stellar feedback at different mass scales. We also check the effect of enhanced star formation events (i.e. starbursts modes), defined in connection with galaxy mergers. We find that, in general, these episodes have a limited impact in the overall star formation histories of model galaxies, even in massive DM halos where merger-driven star formation has often been considered very important.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bias deconstructed: Unravelling the scale dependence of halo bias using real space measurements

We explore the scale dependence of halo bias using real space cross-correlation measurements in N-body simulations and in Pinocchio, an algorithm based on Lagrangian Perturbation Theory. Recent work has shown how to interpret such real space measurements in terms of k-dependent bias in Fourier space, and how to remove the k-dependence to reconstruct the k-independent peak-background split halo bias parameters. We compare our reconstruction of the linear bias, which requires no free parameters, with previous estimates from N-body simulations which were obtained directly in Fourier space at large scales, and find very good agreement. Our reconstruction of the quadratic bias is similarly parameter-free, although in this case there are no previous Fourier space measurements to compare with. Our analysis of N-body simulations explicitly tests the predictions of the excursion set peaks (ESP) formalism of Paranjape et al. (2013) for the scale dependence of bias; we find that the ESP predictions accurately describe our measurements. In addition, our measurements in Pinocchio serve as a useful, successful consistency check between Pinocchio and N-body simulations that is not accessible to traditional measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Galactic Winds in Cosmological Simulations of the Circumgalactic Medium

(Abridged) We explore new observationally-constrained sub-resolution models of galactic outflows and investigate their impact on the circumgalactic medium (CGM) over redshifts z = 2 - 4. We perform cosmological hydrodynamic simulations, including star formation, chemical enrichment, and four cases of SNe-driven outflows: no wind (NW), an energy-driven constant velocity wind (CW), a radially varying wind (RVWa) where the outflow velocity has a positive correlation with galactocentric distance (r), and a RVW with additional dependence on halo mass (RVWb). Overall, we find that the outflows expel metal-enriched gas away from galaxies, significantly quench star formation, and enrich the CGM. At z = 2, the radial profiles of gas properties around galaxy centers are most sensitive to the choice of the wind model for halo masses (10^9 - 10^11) M_sun. We infer that the RVWb model is similar to the NW case, except that it substantially enriches the CGM: the carbon metallicity (Z_C) is 10 times higher in RVWb than in NW at r > R_200; and the warm gas of 10^4 - 10^5 K and delta < 100 is enriched to 50 times higher than in NW. We also find that the impact of models CW and RVWa are similar, with the following differences. RVWa causes a greater suppression of star formation rate at z < 5, and has a higher fraction of low-density (delta < 10), warm-hot (10^4 - 10^6 K) gas than in CW. Outflows in CW produce a higher and earlier enrichment of some IGM phases than RVWa. By visual inspection, we note that the RVWa model shows more pronounced bipolar outflows and galactic disks. We present fitting formulae for [Z_C-delta] and [Z_C-r], also for the abundance of CIV as a function of r. We predict observational diagnostics to distinguish between different outflow scenarios: Z_C of the CGM gas at r = (30 - 300) kpc/h comoving, and CIV fraction of the inner gas at r < (4 - 5) kpc/h comoving.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hierarchical Origins of Observed Galaxy Morphology

Galaxies grow primarily via accretion-driven star formation in discs and merger-driven growth of bulges. These processes are implicit in semi-analytical models of galaxy formation, with bulge growth in particular relating directly to the hierarchical build-up of halos and their galaxies. In this paper, we consider several implementations of two semi-analytical models. Focusing on implementations in which bulges are formed during mergers only, we examine the fractions of elliptical galaxies and both passive and star-forming disk galaxies as functions of stellar and halo mass, for central and satellite systems. This is compared to an observational cross-matched SDSS+RC3 z ~ 0 sample of galaxies with accurate visual morphological classifications and M_{stellar} > 10^10.5 M_{sol}. The models qualitatively reproduce the observed increase of elliptical fraction with stellar mass, and with halo mass for central galaxies, supporting the idea that observed ellipticals form during major mergers. However, the overall elliptical fraction produced by the models is much too high compared with the z ~ 0 data. Since the "passive" -- i.e. non-star-forming -- fractions are approximately reproduced, and since the fraction which are star-forming disc galaxies is also reproduced, the problem is that the models overproduce ellipticals at the expense of passive S0 and spiral galaxies. Bulge-growth implementations (tuned to reproduce simulations) which allow the survival of residual discs in major mergers still destroy too much of the disc. Increasing the lifetime of satellites, or allowing significant disc regrowth around merger remnants, merely increases the fraction of star-forming disc galaxies. Instead, it seems necessary to reduce the mass ratios of merging galaxies, so that most mergers produce modest bulge growth in disc-galaxy remnants instead of ellipticals. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A warm mode of gas accretion on forming galaxies

We present results from high--resolution cosmological hydrodynamical simulations of a Milky--Way-sized halo, aimed at studying the effect of feedback on the nature of gas accretion. Simulations include a model of inter-stellar medium and star formation, in which SN explosions provide effective thermal feedback. We distinguish between gas accretion onto the halo, which occurs when gas particles cross the halo virial radius, and gas accretion onto the central galaxy, which takes place when gas particles cross the inner one-tenth of the virial radius. Gas particles can be accreted through three different channels, depending on the maximum temperature value, $T_{\rm max}$, reached during the particles' past evolution: a cold channel for $T_{\rm max}<2.5 \times 10^5$ K, a hot one for $T>10^6$K, and a warm one for intermediate values of $T_{\rm max}$. We find that the warm channel is at least as important as the cold one for gas accretion onto the central galaxy. This result is at variance with previous findings that the cold mode dominates gas accretion at high redshift. We ascribe this difference to the different supernova feedback scheme implemented in our simulations. While results presented so far in the literature are based on uneffective SN thermal feedback schemes and/or the presence of a kinetic feedback, our simulations include only effective thermal feedback. We argue that observational detections of a warm accretion mode in the high--redshift circum-galactic medium would provide useful constraints on the nature of the feedback that regulates star formation in galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Schmidt-Kennicutt relations in SPH simulations of disc galaxies with effective thermal feedback from supernovae

We study several versions of the Schmidt-Kennicutt (SK) relation obtained for isolated spiral galaxies in TreeSPH simulations run with the GADGET3 code including the novel MUlti-Phase Particle Integrator (MUPPI) algorithm for star formation and stellar feedback. [...] The standard SK relation between surface densities of cold (neutral+molecular) gas and star formation rate of simulated galaxies shows a steepening at low gas surface densities, starting from a knee whose position depends on disc gas fraction: for more gas-rich discs the steepening takes place at higher surface densities. Because gas fraction and metallicity are typically related, this environmental dependence mimics the predictions of models where the formation of H2 is modulated by metallicity. The cold gas surface density at which HI and molecular gas surface densities equate can range from ~10 up to 34 Msun/pc^2. As expected, the SK relation obtained using molecular gas shows much smaller variations among simulations. We find that disc pressure is not well represented by the classical external pressure of a disc in vertical hydrostatic equilibrium. Instead is well fit by the expression P_fit = Sigma_cold sigma_cold kappa / 6, where the three quantities on the right-hand side are cold gas surface density, vertical velocity dispersion and epicyclic frequency. When the "dynamical" SK relation, i.e. the relation that uses gas surface density divided by orbital time, is considered, we find that all of our simulations stay on the same relation. We interpret this as a manifestation of the equilibrium between energy injection and dissipation in stationary galaxy discs, when energetic feedback is effective and pressure is represented by the expression given above. These findings further support the idea that a realistic model of the structure of galaxy discs should take into account energy injection by SNe. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The effect of thermally pulsating asymptotic giant branch stars on the evolution of the rest-frame near-infrared galaxy luminosity function

We address the fundamental question of matching the rest-frame K-band luminosity function (LF) of galaxies over the Hubble time using semi-analytic models, after modification of the stellar population modelling. We include the Maraston evolutionary synthesis models, that feature a higher contribution by the Thermally Pulsating - Asymptotic Giant Branch (TP-AGB) stellar phase, into three different semi-analytic models, namely the De Lucia and Blaizot version of the Munich model, MORGANA and the Menci model. We leave all other input physics and parameters unchanged. We find that the modification of the stellar population emission can solve the mismatch between models and the observed rest-frame K-band luminosity from the brightest galaxies derived from UKIDSS data at high redshift. For all explored semi-analytic models this holds at the redshifts - between 2 and 3 - where the discrepancy was recently pointed out. The reason for the success is that at these cosmic epochs the model galaxies have the right age (~1 Gyr) to contain a well-developed TP-AGB phase which makes them redder without the need of changing their mass or age. At the same time, the known overestimation of the faint end is enhanced in the K-band when including the TP-AGB contribution. At lower redshifts (z<2) some of the explored models deviate from the data. This is due to too short merging timescales and inefficient 'radio-mode' AGN feedback. Our results show that a strong evolution in mass predicted by hierarchical models is compatible with no evolution on the bright-end of the K-band LF from z=3 to the local universe. This means that, at high redshifts and contrary to what is commonly accepted, K-band emission is not necessarily a good tracer of galaxy mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Convergence of Galaxy Properties with Merger Tree Temporal Resolution

Dark matter halo merger trees are now routinely extracted from cosmological simulations of structure formation. These trees are frequently used as inputs to semi-analytic models of galaxy formation to provide the backbone within which galaxy formation takes place. By necessity, these merger trees are constructed from a finite set of discrete "snapshots" of the N-body simulation and so have a limited temporal resolution. To date, there has been little consideration of how this temporal resolution affects the properties of galaxies formed within these trees. In particular, the question of how many snapshots are needed to achieve convergence in galaxy properties has not be answered. Therefore, we study the convergence in the stellar and total baryonic masses of galaxies, distribution of merger times, stellar mass functions and star formation rates in the Galacticus model of galaxy formation as a function of the number of "snapshot" times used to represent dark matter halo merger trees. When utilizing snapshots between z=20 and z=0, we find that at least 128 snapshots are required to achieve convergence to within 5% for galaxy masses. This convergence is obtained for mean quantities averaged over large samples of galaxies - significant variance for individual galaxies remains even when using very large numbers of snapshots. We find only weak dependence of the rate of convergence on the distribution of snapshots in time - snapshots spaced uniformly in the expansion factor, uniformly in the logarithm of expansion factor or uniformly in the logarithm of critical overdensity for collapse work equally well in almost all cases. We provide input parameters to Galacticus which allow this type of convergence study to be tuned to other simulations and to be carried out for other galaxy properties.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The other side of Bulge Formation in a Lambda-CDM cosmology: Bulgeless Galaxies in the Local Universe

We study the physical properties, formation histories, and environment of galaxies without a significant "classical" spheroidal component, as predicted by semi-analytical models of galaxy formation and evolution. This work is complementary to the analysis presented in De Lucia et al., (2011), where we focus on the relative contribution of various physical mechanisms responsible for bulge assembly in a Lambda-CDM cosmology. We find that the fraction of bulgeless galaxies is a strong decreasing function of stellar mass: they represent a negligible fraction of the galaxy population with M* > 10^12 Msun, but dominate at M* < 10^10 Msun. We find a clear dichotomy in this galaxy population, between central galaxies of low-mass dark matter haloes, and satellite galaxies in massive groups/clusters. We show that bulgeless galaxies are relatively young systems, that assemble most of their mass at low-redshift, but they can also host very old stellar populations. Since galaxy-galaxy mergers are assumed to lead to the formation of a spheroidal component, in our models these galaxies form preferentially in low-mass haloes that host a small number of satellites galaxies. We show that the adopted modelling for galaxy mergers represents a key ingredient in determining the actual number of bulgeless galaxies. Our results show that these galaxies are not a rare population in theoretical models: at z~0, galaxies with no classical bulge (but often including galaxies with the equivalent of pseudo-bulges) account for up to 14% of the galaxies with 10^11 < M*/Msun < 10^12.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Times, environments and channels of bulge formation in a LambdaCDM cosmology

We analyze predictions from two independently developed galaxy formation models to study the mechanisms, environments, and characteristic times of bulge formation in a LambdaCDM cosmogony. For each model, we test different prescriptions for bulge formation in order to quantify the relative importance of different channels. Our results show that the strong correlation between galaxy and halo mass for central galaxies, and the richer merger history of more massive systems naturally give rise to a strong correlation between galaxy mass and morphology, and between halo mass and morphological type of central galaxies. Large fractions of the bulge mass are acquired through major and minor mergers, but disc instability plays an important role, particularly for intermediate mass galaxies. We find that the modelling of disc instability events, as well as of the galaxy merger times, can affect significantly the timing of bulge formation, and the relative importance of different channels. Bulge dominated galaxies acquire their morphology through major mergers, but this can be modified by cooling of gas from the surrounding hot halo. We find that disc regrowth is a non negligible component of the evolution of bulge dominated galaxies, particularly for low to intermediate masses, and at high redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A decline and fall in the future of Italian Astronomy?

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A semi-analytic model comparison - gas cooling and galaxy mergers

We use stripped-down versions of three semi-analytic galaxy formation models to study the influence of different assumptions about gas cooling and galaxy mergers. By running the three models on identical sets of merger trees extracted from high-resolution cosmological N-body simulations, we are able to perform both statistical analyses and halo-by-halo comparisons. Our study demonstrates that there is a good statistical agreement between the three models used here, when operating on the same merger trees, reflecting a general agreement in the underlying framework for semi-analytic models. We also show, however, that various assumptions that are commonly adopted to treat gas cooling and galaxy mergers can lead to significantly different results, at least in some regimes. In particular, we find that the different models adopted for gas cooling lead to similar results for mass scales comparable to that of our own Galaxy. Significant differences, however, arise at larger mass scales. These are largely (but not entirely) due to different treatments of the `rapid cooling' regime, and different assumptions about the hot gas distribution. At this mass regime, the predicted cooling rates can differ up to about one order of magnitude, with important implications on the relative weight that these models give to AGN feedback in order to counter-act excessive gas condensation in relatively massive haloes at low redshift. Different assumptions in the modelling of galaxy mergers can also result in significant differences in the timings of mergers, with important consequences for the formation and evolution of massive galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A sub-resolution multiphase interstellar medium model of star formation and SNe energy feedback

We present a new multi-phase sub-resolution model for star formation and feedback in SPH numerical simulations of galaxy formation. Our model, called MUPPI (MUlti-Phase Particle Integrator), describes each gas particle as a multi-phase system, with cold and hot gas phases, coexisting in pressure equilibrium, and a stellar component. Cooling of the hot tenuous gas phase feeds the cold gas phase. Stars are formed out of molecular gas with a given efficiency, which scales with the dynamical time of the cold phase. Our prescription for star formation is not based on imposing the Schmidt-Kennicutt relation, which is instead naturally produced by MUPPI. Energy from supernova explosions is deposited partly into the hot phase of the gas particles, and partly to that of neighboring particles. Mass and energy flows among the different phases of each particle are described by a set of ordinary differential equations which we explicitly integrate for each gas particle, instead of relying on equilibrium solutions. This system of equations also includes the response of the multi-phase structure to energy changes associated to the thermodynamics of the gas. We apply our model to two isolated disk galaxy simulations and two spherical cooling flows. MUPPI is able to reproduce the Schmidt-Kennicutt relation for disc galaxies. It also reproduces the basic properties of the inter-stellar medium in disc galaxies, the surface densities of cold and molecular gas, of stars and of star formation rate, the vertical velocity dispersion of cold clouds and the flows connected to the galactic fountains. Quite remarkably, MUPPI also provides efficient stellar feedback without the need to include a scheme of kinetic energy feedback. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The active and passive populations of Extremely Red Objects

[abridged] The properties of galaxies with the reddest observed R-K colors (Extremely Red Objects, EROs), including their apparent division into passive and obscured active objects with roughly similar number densities, are a known challenge for models of galaxy formation. We produce mock catalogues generated by interfacing the predictions of the semi-analytical MORGANA model for the evolution of galaxies in a Lambda-CDM cosmology with the spectro-photometric + radiative transfer code GRASIL and Infrared (IR) template library to show that the model correctly reproduces number counts, redshift distributions and active fractions of R-K>5 sources. We test the robustness of our results against different dust attenuations and, most importantly, against the inclusion of TP-AGB stars in Simple Stellar Populations used to generate galaxy spectra, and find that the inclusion of TP-AGBs has a relevant effect, in that it allows to increase by a large factor the number of very red active objects at all color cuts. We find that though the most passive and the most obscured active galaxies have a higher probability of being selected as EROs, many EROs have intermediate properties and the population does not show bimodality in specific star formation rate (SSFR). We predict that deep observations in the Far-IR, from 100 to 500 micron, are the most efficient way to constrain the SSFR of these objects; we give predictions for future Herschel observations. Finally, we test whether a simple evolutionary sequence for the formation of z=0 massive galaxies, going through a sub-mm-bright phase and then a ERO phase, are typical in this galaxy formation model. We find that this sequence holds for ~25 per cent of z=0 massive galaxies, while the model typically shows a more complex connection between sub-mm, ERO and massive galaxies. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Faint Lyman-Break galaxies as a crucial test for galaxy formation models

It has recently been shown that galaxy formation models within the LambdaCDM cosmology predict that, compared to the observed population, small galaxies (with stellar masses < 10^{11} M_sun) form too early, are too passive since z ~ 3 and host too old stellar populations at z=0. We then expect an overproduction of small galaxies at z > 4 that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the MORGANA galaxy formation model and GRASIL spectro-photometric + radiative transfer code to generate mock catalogues of deep fields observed with HST-ACS. We add observational noise and the effect of Lyman-alpha emission, and perform color-color selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i-dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of z_{850} ~ 27 V-dropouts at z ~ 5. We quantify the properties of these "excess" galaxies and discuss the implications: these galaxies are hosted in dark matter halos with circular velocities in excess of 100 km s^{-1}, and their suppression may require a deep re-thinking of stellar feedback processes taking place in galaxy formation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Luminosity function and radial distribution of Milky Way Satellites in a LCDM Universe

We study the luminosity function and the radial distribution of satellite galaxies within Milky Way sized haloes as predicted in Cold Dark Matter based models of galaxy formation, making use of numerical N-body techniques as well as three different semi-analytic model (SAMs) galaxy formation codes. We extract merger trees from very high-resolution dissipationless simulations of four Galaxy-sized DM haloes, and use these as common input for the semi-analytic models. We present a detailed comparison of our predictions with the observational data recently obtained on the Milky Way satellite luminosity function (LF). We find that semi-analytic models with rather standard astrophysical ingredients are able to reproduce the observed luminosity function over six orders of magnitude in luminosity, down to magnitudes as faint as M_V=-2. We also perform a comparison with the actual observed number of satellites as a function of luminosity, by applying the selection criteria of the SDSS survey to our simulations instead of correcting the observations for incompleteness. Using this approach we again find good agreement for both the luminosity and radial distributions of MW satellites. We investigate which physical processes in our models are responsible for shaping the predicted satellite LF, and find that tidal destruction, suppression of gas infall by a photo-ionizing background, and supernova feedback all make important contributions. We conclude that the number and luminosity of Milky Way satellites can be naturally accounted for within the (Lambda) Cold Dark Matter paradigm, and this should no longer be considered a problem.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Many Manifestations of Downsizing: Hierarchical Galaxy Formation Models confront Observations

[abridged] It has been widely claimed that several lines of observational evidence point towards a "downsizing" (DS) of the process of galaxy formation over cosmic time. This behavior is sometimes termed "anti-hierarchical", and contrasted with the "bottom-up" assembly of the dark matter structures in Cold Dark Matter models. In this paper we address three different kinds of observational evidence that have been described as DS: the stellar mass assembly, star formation rate and the ages of the stellar populations in local galaxies. We compare a broad compilation of available data-sets with the predictions of three different semi-analytic models of galaxy formation within the Lambda-CDM framework. In the data, we see only weak evidence at best of DS in stellar mass and in star formation rate. We find that, when observational errors on stellar mass and SFR are taken into account, the models acceptably reproduce the evolution of massive galaxies, over the entire redshift range that we consider. However, lower mass galaxies are formed too early in the models and are too passive at late times. Thus, the models do not correctly reproduce the DS trend in stellar mass or the archaeological DS, while they qualitatively reproduce the mass-dependent evolution of the SFR. We demonstrate that these discrepancies are not solely due to a poor treatment of satellite galaxies but are mainly connected to the excessively efficient formation of central galaxies in high-redshift haloes with circular velocities ~100-200 km/s. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

A Lagrangian Dynamical Theory for the Mass Function of Cosmic Structures: I Dynamics

A new theory for determining the mass function of cosmic structures is presented. It relies on a realistic treatment of collapse dynamics. Gravitational collapse is analyzed in the Lagrangian perturbative framework. Lagrangian perturbations provide an approximation of truncated type, i.e. small-scale structure is filtered out. The collapse time is suitably defined as the instant at which orbit crossing takes place. The convergence of the Lagrangian series in predicting the collapse time of a homogeneous ellipsoid is demonstrated; it is also shown that third-order calculations are necessary in predicting collapse. Then, the Lagrangian prediction, with a correction for quasi-spherical perturbations, can be used to determine the collapse time of a homogeneous ellipsoid in a fast and precise way. Furthermore, ellipsoidal collapse can be considered as a particular truncation of the Lagrangian series. Gaussian fields with scale-free power spectra are then considered. The Lagrangian series for the collapse time is found to converge when the collapse time is not large. In this case, ellipsoidal collapse gives a fast and accurate approximation of the collapse time; spherical collapse is found to poorly reproduce the collapse time, even in a statistical sense. Analytical fits of the distribution functions of the inverse collapse times, as predicted by the ellipsoid model and by third-order Lagrangian theory, are given. These will be necessary for a determination of the mass function, which will be given in paper II.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

A Lagrangian Dynamical Theory for the Mass Function of Cosmic Structures: II Statistics

The statistical tools needed to obtain a mass function from realistic collapse time estimates are presented. Collapse dynamics has been dealt with in paper I of this series by means of the powerful Lagrangian perturbation theory and the simple ellipsoidal collapse model. The basic quantity considered here is the inverse collapse time F; it is a non-linear functional of the initial potential, with a non-Gaussian distribution. In the case of sharp k-space smoothing, it is demonstrated that the fraction of collapsed mass can be determined by extending to the F process the diffusion formalism introduced by Bond et al. (1991). The problem is then reduced to a random walk with a moving absorbing barrier, and numerically solved; an accurate analytical fit, valid for small and moderate resolutions, is found. For Gaussian smoothing, the F trajectories are strongly correlated in resolution. In this case, an approximation proposed by Peacock & Heavens (1990) can be used to determine the mass functions. Gaussian smoothing is preferred, as it optimizes the performances of dynamical predictions and stabilizes the F trajectories. The relation between resolution and mass is treated at a heuristic level, and the consequences of this approximation are discussed. The resulting mass functions, compared to the classical Press & Schechter (1974) one, are shifted toward large masses (confirming the findings of Monaco 1995), and tend to give more intermediate-mass objects at the expense of small-mass objects. However, the small-mass part of the mass function, which depends on uncertain dynamics and is likely to be affected by uncertainties in the resolution--mass relation, is not considered a robust prediction of this theory.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

THE LOCAL GALAXY DENSITY AND THE BULGE-TO-DISK RATIO OF DISK GALAXIES

Relying on samples of disk galaxies for which a detailed photometric bulge/disk decomposition has been provided in the literature, we examine the dependence of the bulge-to-disk luminosity ratio (B/D) on the blue absolute luminosity and on the environmental density. In our statistical analysis of various B/D data sets we pay particular attention to disentangling the role played by the galaxy morphology--galaxy density relation. Besides, we focus our attention on nearby ($z<0.01$) galaxies, for which we can provide a three-dimensional characterization of the local galaxy density. We find that the observed tendency of galaxies to have greater B/D with increasing galaxy density simply reflects the average decline of B/D towards later morphological types together with the morphology--density relation. This relation tends to give rise also to a greater proportion of bright bulges in denser regions, because the decrease of B/D towards later types is mostly due to a dimming of the bulge rather than to a brightening of the disk. But when we remove the effect induced by the morphology--density relation, we detect no clear evidence of a dependence of B/D on galaxy density. Furthermore, B/D turns out to be substantially unrelated to the blue absolute magnitude of the galaxy. We briefly discuss to what extent our results (partially) disagree with previous claims.

0 citations0 reviewsNo code linkedOpen access