Dark Halo and Disk Galaxy Scaling Laws in Hierarchical Universes

We use cosmological N-body/gasdynamical simulations that include star formation and feedback to examine the proposal that scaling laws between the total luminosity, rotation speed, and angular momentum of disk galaxies reflect analogous correlations between the structural parameters of their surrounding dark matter halos. The numerical experiments follow the formation of galaxy-sized halos in two Cold Dark Matter dominated universes: the standard Omega=1 CDM scenario and the currently popular LCDM model. We find that the slope and scatter of the I-band Tully-Fisher relation are well reproduced in the simulations, although not, as proposed in recent work, as a result of the cosmological equivalence between halo mass and circular velocity: large systematic variations in the fraction of baryons that collapse to form galaxies and in the ratio between halo and disk circular velocities are observed in our numerical experiments. The Tully-Fisher slope and scatter are recovered in this model as a direct result of the dynamical response of the halo to the assembly of the luminous component of the galaxy. We conclude that models that neglect the self-gravity of the disk and its influence on the detailed structure of the halo cannot be used to derive meaningful estimates of the scatter or slope of the Tully-Fisher relation. Our models fail, however, to match the zero-point of the Tully-Fisher relation, as well as that of the relation linking disk rotation speed and angular momentum. These failures can be traced, respectively, to the excessive central concentration of dark halos formed in the Cold Dark Matter cosmogonies we explore and to the formation of galaxy disks as the final outcome of a sequence of merger events. (abridged)Comment: submitted to The Astrophysical Journa

The cosmological origin of the Tully-Fisher relation

We use high-resolution cosmological simulations that include the effects of gasdynamics and star formation to investigate the origin of the Tully-Fisher relation in the standard Cold Dark Matter cosmogony. Luminosities are computed for each model galaxy using their full star formation histories and the latest spectrophotometric models. We find that at z=0 the stellar mass of model galaxies is proportional to the total baryonic mass within the virial radius of their surrounding halos. Circular velocity then correlates tightly with the total luminosity of the galaxy, reflecting the equivalence between mass and circular velocity of systems identified in a cosmological context. The slope of the relation steepens slightly from the red to the blue bandpasses, and is in fairly good agreement with observations. Its scatter is small, decreasing from \~0.45 mag in the U-band to ~0.34 mag in the K-band. The particular cosmological model we explore here seems unable to account for the zero-point of the correlation. Model galaxies are too faint at z=0 (by about two magnitudes) if the circular velocity at the edge of the luminous galaxy is used as an estimator of the rotation speed. The Tully-Fisher relation is brighter in the past, by about ~0.7 magnitudes in the B-band at z=1, at odds with recent observations of z~1 galaxies. We conclude that the slope and tightness of the Tully-Fisher relation can be naturally explained in hierarchical models but that its normalization and evolution depend strongly on the star formation algorithm chosen and on the cosmological parameters that determine the universal baryon fraction and the time of assembly of galaxies of different mass.Comment: 5 pages, 4 figures included, submitted to ApJ (Letters

On the origin of the Tully-Fisher relation

We discuss the origin of the Tully-Fisher (TF) relation using the $N$-body/SPH method, which includes cooling, star formation and stellar feedback of energy, mass and metals. We consider initially rotating overdense spheres, and trace formation processes of disk galaxies from $z=25$ to $z=0$ in the Cold Dark Matter (CDM) cosmology. To clarify the origin of the TF relation, we simulate formation of 14 galaxies with different masses and spin parameters, and compute observable values, such as the total magnitude and the line-width. We find that the simulated galaxies reproduce the slope and scatter of the TF relation: the slope is originated in the difference of total galactic masses, and the scatter is produced by the difference of initial spin parameters. As well as the TF relation, observed features of spiral galaxies, such as the exponential light-profile and the flat rotation curve, are reproduced in our simulations, which were assumed {\it a priori} in past semi-analytical approaches.Comment: 11 pages, including 6 figures, submitted to Ap

Simulations of galaxy formation in a Λ cold dark matter universe : I : dynamical and photometric properties of a simulated disk galaxy.

We present a detailed analysis of the dynamical and photometric properties of a disk galaxy simulated in the cold dark matter (CDM) cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z1 that lead to the formation of two distinct dynamical components: a spheroid of mostly old stars and a rotationally supported disk of younger stars. The surface brightness profile is very well approximated by the superposition of an R1/4 spheroid and an exponential disk. Each photometric component contributes a similar fraction of the total luminosity of the system, although less than a quarter of the stars form after the last merger episode at z1. In the optical bands the surface brightness profile is remarkably similar to that of Sab galaxy UGC 615, but the simulated galaxy rotates significantly faster and has a declining rotation curve dominated by the spheroid near the center. The decline in circular velocity is at odds with observation and results from the high concentration of the dark matter and baryonic components, as well as from the relatively high mass-to-light ratio of the stars in the simulation. The simulated galaxy lies 1 mag off the I-band Tully-Fisher relation of late-type spirals but seems to be in reasonable agreement with Tully-Fisher data on S0 galaxies. In agreement with previous simulation work, the angular momentum of the luminous component is an order of magnitude lower than that of late-type spirals of similar rotation speed. This again reflects the dominance of the slowly rotating, dense spheroidal component, to which most discrepancies with observation may be traced. On its own, the disk component has properties rather similar to those of late-type spirals: its luminosity, its exponential scale length, and its colors are all comparable to those of galaxy disks of similar rotation speed. This suggests that a different form of feedback than adopted here is required to inhibit the efficient collapse and cooling of gas at high redshift that leads to the formation of the spheroid. Reconciling, without fine-tuning, the properties of disk galaxies with the early collapse and high merging rates characteristic of hierarchical scenarios such as CDM remains a challenging, yet so far elusive, proposition

Stars beyond Galaxies: The Origin of Extended Luminous Halos around Galaxies

(Abridged) We use numerical simulations to investigate the origin and structure of the luminous halos that surround isolated galaxies. These stellar structures extend out to several hundred kpc away from a galaxy, and consist of stars shed by merging subunits during the many accretion events that characterize the hierarchical assembly of galaxies. Such origin suggests that outer luminous halos are ubiquitous and that they should appear as an excess of light over extrapolations of the galaxy's inner profile beyond its traditional luminous radius. The mass profile of the accreted stellar component is well approximated by a model where the logarithmic slope steepens monotonically with radius; from -3 at the luminous edge of the galaxy to -4 or steeper near the virial radius of the system. Such spatial distribution is consistent with that of Galactic and M31 globular clusters, suggesting that many of the globulars were brought in by accretion events, in a manner akin to the classic Searle-Zinn scenario. The outer stellar spheroid is supported by a velocity dispersion tensor with a substantial and radially increasing radial anisotropy. These properties distinguish the stellar halo from the dark matter component, which is more isotropic in velocity space, as well as from some tracers of the outer spheroid such as satellite galaxies. Most stars in the outer halo formed in progenitors that have since merged with the central galaxy; very few stars in the halo are contributed by satellites that survive as self-bound entities at the present. These features are in reasonable agreement with recent observations of the outer halo of the MW, of M31, and of other isolated spirals, and suggest that all of these systems underwent an early period of active merging, as envisioned in hierarchical models of galaxy formation.Comment: Submitted to MNRAS, 13 pages, 12 figure

The Effects of a Photoionizing UV Background on the Formation of Disk Galaxies

We use high resolution N-body/gasdynamical simulations to investigate the effects of a photoionizing UV background on the assembly of disk galaxies in hierarchically clustering universes. We focus on the mass and rotational properties of gas that can cool to form centrifugally supported disks in dark matter halos of different mass. Photoheating can significantly reduce the amount of gas that can cool in galactic halos. Depending on the strength of the UV background field, the amount of cooled gas can be reduced by up to $50\%$ in systems with circular speeds in the range $80$-$200$ \kms. The magnitude of the effect, however, is not enough to solve the ``overcooling'' problem that plagues hierarchical models of galaxy formation if the UV background is chosen to be consistent with estimates based on recent observations of QSO absorption systems. Photoionization has little effect on the collapse of gas at high redshift and affects preferentially gas that is accreted at late times. Since disks form inside-out, accreting higher angular momentum gas at later times, disks formed in the presence of a UV background have spins that are even smaller than those formed in simulations that do not include the effects of photoionization. This exacerbates the angular momentum problem that afflicts hierarchical models of disk formation. We conclude that photoionization cannot provide the heating mechanism required to reconcile hierarchically clustering models with observations. Energy feedback and enrichment processes from the formation and evolution of stars must therefore be indispensable ingredients for any successful model of the formation of disk galaxies.Comment: 36 pages, w/ embedded figures, submitted to ApJ. Also available at http://penedes.as.arizona.edu/~jfn/preprints/dskform.ps.g

The Angular Momentum Distribution of Gas and Dark Matter in Galactic Halos

(Abridged) We report results of a series of non radiative N-body/SPH simulations in a LCDM cosmology. We find that the spin of the baryonic component is on average larger than that of the dark matter (DM) component and we find this effect to be more pronounced at lower redshifts. A significant fraction f of gas has negative angular momentum and this fraction is found to increase with redshift. We describe a toy model in which the tangential velocities of particles are smeared by Gaussian random motions. This model is successful in explaining some of the angular momentum properties. We compare and contrast various techniques to determine the angular momentum distributions (AMDs). We show that broadening of velocity dispersions is unsuitable for making comparisons between gas and DM. We smooth the angular momentum of the particles over a fixed number of neighbors. We find that an analytical function based on gamma distribution can be used to describe a wide variety of profiles, with just one parameter \alpha. The distribution of the shape parameter $\alpha$ for both gas and DM follows roughly a log-normal distribution. The mean and standard deviation of log(\alpha) for gas is -0.04 and 0.11 respectively. About 90-95% of halos have \alpha<1.3, while exponential disks in NFW halos would require 1.3<\alpha<1.6. This implies that a typical halo in simulations has an excess of low angular momentum material as compared to that of observed exponential disks, a result which is consistent with the findings of earlier works. \alpha for gas is correlated with that of DM but they have a significant scatter =1.09 \pm 0.2. \alpha_Gas is also biased towards slightly higher values compared to \alpha_DM.Comment: 19 pages, 32 figures (replaced to correct a typo in the authors field in the above line, paper unchanged