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

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

A Unified Scaling Law in Spiral Galaxies

We investigate the origin of a unified scaling relation in spiral galaxies. Observed spiral galaxies are spread on a plane in the three-dimensionallogarithmic space of luminosity L, radius R and rotation velocity V. The plane is expressed as $L \propto (V R)^{\alpha}$ in I-passband, where $\alpha$ is a constant. On the plane, observed galaxies are distributed in an elongated region which looks like the shape of a surfboard. The well-known scaling relations, L-V (Tully-Fisher relation), V-R (also the Tully-Fisher relation) and R-L (Freeman's law), can be understood as oblique projections of the surfboard-like plane into 2-D spaces. This unified interpretation of the known scaling relations should be a clue to understand the physical origin of all the relations consistently. Furthermore, this interpretation can also explain why previous studies could not find any correlation between TF residuals and radius. In order to clarify the origin of this plane, we simulate formation and evolution of spiral galaxies with the N-body/SPH method, including cooling, star formation and stellar feedback. Initial conditions are set to isolated 14 spheres with two free parameters, such as mass and angular momentum. The CDM (h=0.5, $\Omega_0=1$) cosmology is considered as a test case. The simulations provide the following two conclusions: (a) The slope of the plane is well reproduced but the zero-point is not. This zero-point discrepancy could be solved in a low density ($\Omega_00.5) cosmology. (b) The surfboard-shaped plane can be explained by the control of galactic mass and angular momentum.Comment: Accepted for publication in ApJ Letters. 6 pages including 2 figure

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

Star Formation, Supernovae Feedback and the Angular Momentum Problem in Numerical CDM Cosmogony: Half Way There?

We present a smoothed particle hydrodynamic (SPH) simulation that reproduces a galaxy that is a moderate facsimile of those observed. The primary failing point of previous simulations of disk formation, namely excessive transport of angular momentum from gas to dark matter, is ameliorated by the inclusion of a supernova feedback algorithm that allows energy to persist in the model ISM for a period corresponding to the lifetime of stellar associations. The inclusion of feedback leads to a disk at a redshift $z=0.52$, with a specific angular momentum content within 10% of the value required to fit observations. An exponential fit to the disk baryon surface density gives a scale length within 17% of the theoretical value. Runs without feedback, with or without star formation, exhibit the drastic angular momentum transport observed elsewhere.Comment: 4 pages, 3 figures, accepted for publication in ApJ Letter