An analytical model for the accretion of dark matter subhalos

An analytical model is developed for the mass function of cold dark matter subhalos at the time of accretion and for the distribution of their accretion times. Our model is based on the model of Zhao et al. (2009) for the median assembly histories of dark matter halos, combined with a simple log-normal distribution to describe the scatter in the main-branch mass at a given time for halos of the same final mass. Our model is simple, and can be used to predict the un-evolved subhalo mass function, the mass function of subhalos accreted at a given time, the accretion-time distribution of subhalos of a given initial mass, and the frequency of major mergers as a function of time. We test our model using high-resolution cosmological $N$-body simulations, and find that our model predictions match the simulation results remarkably well. Finally, we discuss the implications of our model for the evolution of subhalos in their hosts and for the construction of a self-consistent model to link galaxies and dark matter halos at different cosmic times.Comment: 14 pages, 10 figures (caption for figure 10 fixed). Accepted for publication in Ap

A Bayesian approach to the semi-analytic model of galaxy formation: methodology

We believe that a wide range of physical processes conspire to shape the observed galaxy population but we remain unsure of their detailed interactions. The semi-analytic model (SAM) of galaxy formation uses multi-dimensional parameterisations of the physical processes of galaxy formation and provides a tool to constrain these underlying physical interactions. Because of the high dimensionality, the parametric problem of galaxy formation may be profitably tackled with a Bayesian-inference based approach, which allows one to constrain theory with data in a statistically rigorous way. In this paper we develop a SAM in the framework of Bayesian inference. We show that, with a parallel implementation of an advanced Markov-Chain Monte-Carlo algorithm, it is now possible to rigorously sample the posterior distribution of the high-dimensional parameter space of typical SAMs. As an example, we characterise galaxy formation in the current $\Lambda$CDM cosmology using the stellar mass function of galaxies as an observational constraint. We find that the posterior probability distribution is both topologically complex and degenerate in some important model parameters, suggesting that thorough explorations of the parameter space are needed to understand the models. We also demonstrate that because of the model degeneracy, adopting a narrow prior strongly restricts the model. Therefore, the inferences based on SAMs are conditional to the model adopted. Using synthetic data to mimic systematic errors in the stellar mass function, we demonstrate that an accurate observational error model is essential to meaningful inference.Comment: revised version to match published article published in MNRA

Impact of Supernova feedback on the Tully-Fisher relation

Recent observational results found a bend in the Tully-Fisher Relation in such a way that low mass systems lay below the linear relation described by more massive galaxies. We intend to investigate the origin of the observed features in the stellar and baryonic Tully-Fisher relations and analyse the role played by galactic outflows on their determination. Cosmological hydrodynamical simulations which include Supernova feedback were performed in order to follow the dynamical evolution of galaxies. We found that Supernova feedback is a fundamental process in order to reproduce the observed trends in the stellar Tully-Fisher relation. Simulated slow rotating systems tend to have lower stellar masses than those predicted by the linear fit to the massive end of the relation, consistently with observations. This feature is not present if Supernova feedback is turned off. In the case of the baryonic Tully-Fisher relation, we also detect a weaker tendency for smaller systems to lie below the linear relation described by larger ones. This behaviour arises as a result of the more efficient action of Supernovae in the regulation of the star formation process and in the triggering of powerful galactic outflows in shallower potential wells which may heat up and/or expel part of the gas reservoir.Comment: 10 pages, 9 figures, accepted for publication in A&

Fingerprints of the hierarchical building up of the structure on the gas kinematics of galaxies

Recent observational and theoretical works have suggested that the Tully-Fisher Relation might be generalised to include dispersion-dominated systems by combining the rotation and dispersion velocity in the definition of the kinematical indicator. Mergers and interactions have been pointed out as responsible of driving turbulent and disordered gas kinematics, which could generate Tully-Fisher Relation outliers. We intend to investigate the gas kinematics of galaxies by using a simulated sample which includes both, gas disc-dominated and spheroid-dominated systems. Cosmological hydrodynamical simulations which include a multiphase model and physically-motivated Supernova feedback were performed in order to follow the evolution of galaxies as they are assembled. Both the baryonic and stellar Tully-Fisher relations for gas disc-dominated systems are tight while, as more dispersion-dominated systems are included, the scatter increases. We found a clear correlation between $\sigma / V_{\rm rot}$ and morphology, with dispersion-dominated systems exhibiting the larger values ($> 0.7$). Mergers and interactions can affect the rotation curves directly or indirectly inducing a scatter in the Tully-Fisher Relation larger than the simulated evolution since $z \sim 3$. Kinematical indicators which combine rotation velocity and dispersion velocity can reduce the scatter in the baryonic and the stellar mass-velocity relations. Our findings also show that the lowest scatter in both relations is obtained if the velocity indicators are measured at the maximum of the rotation curve. Moreover, the rotation velocity estimated at the maximum of the gas rotation curve is found to be the best proxy for the potential well regardless of morphology.Comment: 16 pages, 10 figures, accepted for publication in A&

On the Angular Correlation Function of SZ Clusters : Extracting cosmological information from a 2D catalog

We discuss the angular correlation function of Sunyaev-Zel'dovich (SZ)-detected galaxy clusters as a cosmological probe. As a projection of the real-space cluster correlation function, the angular function samples the underlying SZ catalog redshift distribution. It offers a way to study cosmology and cluster evolution directly with the two-dimensional catalog, even before extensive follow-up observations, thereby facilitating the immediate scientific return from SZ surveys. As a simple illustration of the information content of the angular function, we examine its dependence on the parameter pair Om_m, sigma_8 in flat cosmologies. We discuss sources of modeling uncertainty and consider application to the future Planck SZ catalog, showing how these two parameters and the normalization of the SZ flux-mass relation can be simultaneously found when the local X-ray cluster abundance constraint is included.Comment: 11 pages, 5 figures. A&A, 410, 767; corrected typo, published versio

The Impact of Cooling and Feedback on Disc Galaxies

We present detailed, analytical models for the formation of disc galaxies to investigate the impact that cooling and feedback have on their structural properties. In particular, we investigate which observables extracted directly from the models are best suited as virial mass estimators, and to what extent they allow the recovery of the model input parameters regarding the feedback and cooling efficiencies. Contrary to naive expectations, the luminosities and circular velocities of disc galaxies are extremely poor indicators of total virial mass. Instead, we show that the product of disc scale length and rotation velocity squared yields a much more robust estimate. We show that feedback can cause a narrow correlation between galaxy mass fraction and halo spin parameter, similar to that found recently by van den Bosch, Burkert and Swaters from an analysis of dwarf galaxy rotation curves. Finally we investigate the impact that cooling and feedback have on the colors, metallicities, star formation histories and Tully-Fisher relation of disc galaxies.Comment: 20 pages, 12 figures. To be published in MNRA

The baryonic Tully-Fisher relation and galactic outflows

Most of the baryons in the Universe are not in the form of stars and cold gas in galaxies. Galactic outflows driven by supernovae/stellar winds are the leading mechanism for explaining this fact. The scaling relation between galaxy mass and outer rotation velocity (also known as the baryonic Tully-Fisher relation, BTF) has recently been used as evidence against this viewpoint. We use a LCDM based semi-analytic disk galaxy formation model to investigate these claims. In our model, galaxies with less efficient star formation and higher gas fractions are more efficient at ejecting gas from galaxies. This is due to the fact that galaxies with less efficient star formation and higher gas fractions tend to live in dark matter haloes with lower circular velocities, from which less energy is required to escape the potential well. In our model the intrinsic scatter in the BTF is 0.15 dex, and mostly reflects scatter in dark halo concentration. The observed scatter, equal to 0.24 dex, is dominated by measurement errors. The best estimate for the intrinsic scatter is that it is less than 0.15 dex, and thus our LCDM based model (which does not include all possible sources of scatter) is only just consistent with this. In our model, gas rich galaxies, at fixed virial velocity (V_vir), with lower stellar masses have lower baryonic masses. This is consistent with the expectation that galaxies with lower stellar masses have had less energy available to drive an outflow. However, when the outer rotation velocity (V_flat) is used the correlation has the opposite sign, with a slope in agreement with observations. This is due to scatter in the relation between V_flat and V_vir. In summary, contrary to some previous claims, we show that basic features of the BTF are consistent with a LCDM based model in which the low efficiency of galaxy formation is determined by galactic outflows.Comment: 7 pages, 4 figures, accepted to MNRA

An Investigation of Sloan Digital Sky Survey Imaging Data and Multi-Band Scaling Relations of Spiral Galaxies (with Dynamical Information)

We have compiled a sample of 3041 spiral galaxies with multi-band gri imaging from the Sloan Digital Sky Survey (SDSS) Data Release 7 and available galaxy rotational velocities derived from HI line widths. We compare the data products provided through the SDSS imaging pipeline with our own photometry of the SDSS images, and use the velocities (V) as an independent metric to determine ideal galaxy sizes (R) and luminosities (L). Our radial and luminosity parameters improve upon the SDSS DR7 Petrosian radii and luminosities through the use of isophotal fits to the galaxy images. This improvement is gauged via VL and RV relations whose respective scatters are reduced by ~8% and ~30% compared to similar relations built with SDSS parameters. The tightest VRL relations are obtained with the i-band radius, R235i, measured at 23.5 mag/arcsec^-2, and the luminosity L235i, measured within R235i. Our VRL scaling relations compare well, both in scatter and slope, with similar studies (such comparisons however depend sensitively on the nature and size of the compared samples). The typical slopes, b, and observed scatters, sigma, of the i-band VL, RL and RV relations are bVL=0.27+/-0.01, bRL=0.41+/-0.01, bRV=1.52+/-0.07, and sigmaVL=0.074, sigmaRL=0.071, sigmaRV=0.154 dex. Similar results for the SDSS g and r bands are also provided. Smaller scatters may be achieved for more pruned samples. We also compute scaling relations in terms of the baryonic mass (stars + gas), Mbar, ranging from 10^8.7 Msol to 10^11.6 Msol. Our baryonic velocity-mass (VM) relation has slope 0.29+/-0.01 and a measured scatter sigma_meas = 0.076 dex. While the observed VL and VM relations have comparable scatter, the stellar and baryonic VM relations may be intrinsically tighter, and thus potentially more fundamental, than other VL relations of spiral galaxies.Comment: Submitted to MNRAS, comments welcom

Cluster abundances and S-Z power spectra: effects of non-Gaussianity and early dark energy

In the standard Lambda CDM cosmological model with a Gaussian primordial density fluctuation field, the relatively low value of the mass variance parameter (sigma_8=0.74{+0.05}{-0.06}, obtained from the WMAP 3-year data) results in a reduced likelihood that the measured level of CMB anisotropy on the scales of clusters is due to the Sunyaev-Zeldovich (S-Z) effect. To assess the feasibility of producing higher levels of S-Z power, we explore two alternative models which predict higher cluster abundance. In the first model the primordial density field has a chi^2_1 distribution, whereas in the second an early dark energy component gives rise to the desired higher cluster abundance. We carry out the necessary detailed calculations of the levels of S-Z power spectra, cluster number counts, and angular 2-point correlation function of clusters, and compare (in a self-consistent way) their predicted redshift distributions. Our results provide a sufficient basis upon which the viability of the three models may be tested by future high quality measurements.Comment: 12 pages, 5 figures, accepted for publication in MNRA

Clues on the origin of galactic angular momentum from looking at galaxy pairs

We search for correlations between the spin in pairs of spiral galaxies, to study if the angular momentum gain for each galaxy was the result of tidal torques imprint by the same tidal field. To perform our study we made use of a sample of galaxy pairs identified using the Sloan Digital Sky Survey. We find a weak, but statistically significant correlation between the spin magnitude of neighbouring galaxies, but no clear alignment between their orientation. We show that events such as interactions with close neighbours play an important role in the value of the spin for the final configuration, as we find these interactions tend to reduce the value of the $\lambda$ spin parameter of late-type galaxies considerably, with dependence on the morphology of the neighbour. This implies that the original tidal field for each pair could have been similar, but the redistribution of angular momentum at later stages of evolution is important.Comment: 10 pages, 4 figures, 1 table. Replaced to match the version accepted for publication in MNRA