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 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

Accuracy of Mesh Based Cosmological Hydrocodes: Tests and Corrections

We perform a variety of tests to determine the numerical resolution of the cosmological TVD eulerian code developed by Ryu et al (1993). Tests include 512^3 and 256^3 simulations of a Pk=k^{-1} spectrum to check for self-similarity and comparison of results with those from higher resolution SPH and grid-based calculations (Frenk et al 1998). We conclude that in regions where density gradients are not produced by shocks the code degrades resolution with a Gaussian smoothing (radius) length of 1.7 cells. At shock caused gradients (for which the code was designed) the smoothing length is 1.1 cells. Finally, for \beta model fit clusters, we can approximately correct numerical resolution by the transformation R^2_{core}\to R^2_{core}-(C\Delta l)^2, where \Delta l is the cell size and C=1.1-1.7. When we use these corrections on our previously published computations for the SCDM and \Lambda CDM models we find luminosity weighted, zero redshift, X-ray cluster core radii of (210\pm 86, 280\pm 67)h^{-1}kpc, respectively, which are marginally consistent with observed (Jones & Forman 1992) values of 50-200h^{-1}kpc. Using the corrected core radii, the COBE normalized SCDM model predicts the number of bright L_x>10^{43}erg/s clusters too high by a factor of \sim 20 and the \Lambda CDM model is consistent with observations.Comment: ApJ in press (1999

Solar comparison spectra, 1.0-2.5 mu, from altitudes 1.5-12.5 km

Solar and telluric infrared spectra from altitudes between 1.5 and 12.5 k

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 in Isolated Disk Galaxies. I. Models and Characteristics of Nonlinear Gravitational Collapse

We model gravitational collapse leading to star formation in a wide range of isolated disk galaxies using a three-dimensional, smoothed particle hydrodynamics code. The model galaxies include a dark matter halo and a disk of stars and isothermal gas. Absorbing sink particles are used to directly measure the mass of gravitationally collapsing gas. They reach masses characteristic of stellar clusters. In this paper, we describe our galaxy models and numerical methods, followed by an investigation of the gravitational instability in these galaxies. Gravitational collapse forms star clusters with correlated positions and ages, as observed, for example, in the Large Magellanic Cloud. Gravitational instability alone acting in unperturbed galaxies appears sufficient to produce flocculent spiral arms, though not more organized patterns. Unstable galaxies show collapse in thin layers in the galactic plane; associated dust will form thin dust lanes in those galaxies, in agreement with observations. (abridged)Comment: 49 pages, 22 figures, to appear in ApJ (July, 2005), version with high quality color images can be fond in http://research.amnh.org/~yuexing/astro-ph/0501022.pd

Physical Bias of Galaxies From Large-Scale Hydrodynamic Simulations

We analyze a new large-scale ($100h^{-1}$Mpc) numerical hydrodynamic simulation of the popular $\Lambda$CDM cosmological model, including in our treatment dark matter, gas and star-formation, on the basis of standard physical processes. The method, applied with a numerical resolution of $<200h^{-1}$kpc (which is still quite coarse for following individual galaxies, especially in dense regions), attempts to estimate where and when galaxies form. We then compare the smoothed galaxy distribution with the smoothed mass distribution to determine the "bias" defined as $b\equiv (\delta M/M)_{gal}/(\delta M/M)_{total}$ on scales large compared with the code numerical resolution (on the basis of resolution tests given in the appendix of this paper). We find that (holding all variables constant except the quoted one) bias increases with decreasing scale, with increasing galactic age or metallicity and with increasing redshift of observations. At the $8h^{-1}$Mpc fiducial comoving scale bias (for bright regions) is 1.35 at $z=0$ reaching to 3.6 at $z=3$, both numbers being consistent with extant observations. We also find that $(10-20)h^{-1}$Mpc voids in the distribution of luminous objects are as observed (i.e., observed voids are not an argument against CDM-like models) and finally that the younger systems should show a colder Hubble flow than do the early type galaxies (a testable proposition). Surprisingly, little evolution is found in the amplitude of the smoothed galaxy-galaxy correlation function (as a function of {\it comoving} separation). Testing this prediction vs observations will allow a comparison between this work and that of Kauffmann et al which is based on a different physical modelingmethod.Comment: in press, ApJ, 26 latex pages plus 7 fig

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

Does the Superior Colliculus Control Perceptual Sensitivity or Choice Bias during Attention? Evidence from a Multialternative Decision Framework

Distinct networks in the forebrain and the midbrain coordinate to control spatial attention. The critical involvement of the superior colliculus (SC)—the central structure in the midbrain network—in visuospatial attention has been shown by four seminal, published studies in monkeys (Macaca mulatta) performing multialternative tasks. However, due to the lack of a mechanistic framework for interpreting behavioral data in such tasks, the nature of the SC's contribution to attention remains unclear. Here we present and validate a novel decision framework for analyzing behavioral data in multialternative attention tasks. We apply this framework to re-examine the behavioral evidence from these published studies. Our model is a multidimensional extension to signal detection theory that distinguishes between two major classes of attentional mechanisms: those that alter the quality of sensory information or “sensitivity,” and those that alter the selective gating of sensory information or “choice bias.” Model-based simulations and model-based analyses of data from these published studies revealed a converging pattern of results that indicated that choice-bias changes, rather than sensitivity changes, were the primary outcome of SC manipulation. Our results suggest that the SC contributes to attentional performance predominantly by generating a spatial choice bias for stimuli at a selected location, and that this bias operates downstream of forebrain mechanisms that enhance sensitivity. The findings lead to a testable mechanistic framework of how the midbrain and forebrain networks interact to control spatial attention