Showing papers by "Julio F. Navarro published in 2002"

The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study

Abstract: We present a comprehensive set of convergence tests which explore the role of various numerical parameters on the equilibrium structure of a simulated dark matter halo. We report results obtained with two independent, state-of-the-art, multi-stepping, parallel N-body codes: PKDGRAV and GADGET. We find that convergent mass profiles can be obtained for suitable choices of the gravitational softening, timestep, force accuracy, initial redshift, and particle number. For softenings chosen so that particle discreteness effects are negligible, convergence in the circular velocity is obtained at radii where the following conditions are satisfied: (i) the timestep is much shorter than the local orbital timescale; (ii) accelerations do not exceed a characteristic acceleration imprinted by the gravitational softening; and (iii) enough particles are enclosed so that the collisional relaxation timescale is longer than the age of the universe. The most stringent requirement for convergence is typically that imposed on the particle number by the collisional relaxation criterion, which implies that in order to estimate accurate circular velocities at radii where the density contrast may reach $\sim 10^6$, the region must enclose of order 3000 particles (or more than a few times $10^6$ within the virial radius). Applying these criteria to a galaxy-sized $\Lambda$CDM halo, we find that the spherically-averaged density profile becomes progressively shallower from the virial radius inwards, reaching a logarithmic slope shallower than -1.2 at the innermost resolved point, $r \sim 0.005 r_{200}$, with little evidence for convergence to a power-law behaviour in the inner regions.

Simulations of Galaxy Formation in a Lambda CDM Universe II: The Fine Structure of Simulated Galactic Disks

Abstract: We present a detailed analysis of the dynamical properties of a simulated disk galaxy assembled hier- archically in theCDM cosmogony. At z = 0, two distinct dynamical components are easily identified solely on the basis of the orbital parameters of stars in the galaxy: a slowly rotating, centrally concen- trated spheroid and a disk-like component largely supported by rotation. These components are also clearly recognized in the surface brightness profile of the galaxy, which can be very well approximated by the superposition of an R 1/4 spheroid and an exponential disk. However, neither does the dynamically- identified spheroid follow de Vaucouleurs' law nor is the disk purely exponential, a result which calls for caution when estimating the importance of the disk from traditional photometric decomposition techniques. The disk may be further decomposed into a thin, dynamically cold component with stars on nearly circular orbits and a hotter, thicker component with orbital parameters transitional between the thin disk and the spheroid. Supporting evidence for the presence of distinct thick and thin disk components is found, as in the Milky Way, in the double-exponential vertical structure of the disk and in abrupt changes in the vertical velocity distribution as a function of age. The dynamical origin of these components offers intriguing clues to the assembly of spheroids and disks in the Milky Way and other spirals. The spheroid is old, and has essentially no stars younger than the time elapsed since the last major accretion event; � 8 Gyr ago for the system we consider here. The majority of thin disk stars, on the other hand, form after the merging activity is over, although a significant fraction (� 15%) of thin-disk stars are old enough to predate the last major merger event. This unexpected population of old disk stars consists mainly of the tidal debris of satellites whose orbital plane was coincident with the disk and whose orbits were circularized by dynamical friction prior to full disruption. More than half of the stars in the thick disk share this origin, part of a trend that becomes more pronounced with age: nine out of ten stars presently in the old (τ �10 Gyr) disk component were actually brought into the disk by satellites. By contrast, only one in two stars belonging to the old spheroid are tidal debris; the rest may be traced to a major merger event that dispersed the luminous progenitor at z � 1.5 and seeded the formation of the spheroid. Our results highlight the role of satellite accretion events in shaping the disk—as well as the spheroidal—component and reveal some of the clues to the assembly process of a galaxy preserved in the detailed dynamics of old stellar populations. Subject headings: cosmology, dark matter, galaxies: formation, galaxies: structure

Simulations of Galaxy Formation in a Lambda CDM Universe I: Dynamical and Photometric Properties of a Simulated Disk Galaxy

Abstract: We analyze the properties of a disk galaxy simulated with unprecedented numerical resolution in the Lambda CDM cosmogony. The galaxy is assembled through a number of high-redshift mergers followed by a period of quiescent accretion after z~1 which lead to the formation of two distinct dynamical components: a spheroid and a disk. The surface brightness profile is very well approximated by the superposition of an R^{1/4} spheroid and an exponential disk. The surface brightness profile is remarkably similar to that of Sab galaxy UGC615, 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. The angular momentum of the luminous component is an order of magnitude lower than that of late-type spirals of similar rotation speed. This reflects the dominance of the slowly-rotating, dense spheroidal component. The disk component, on the other hand, has properties rather similar to those of late-type spirals. This suggests that a different form of feedback than adopted in this simulation is required to inhibit the efficient collapse and cooling of gas at high redshift that leads to the formation of the spheroid. Reconciling the properties of disk galaxies with the early collapse and high merging rates characteristic of hierarchical scenarios such as Lambda CDM remains a challenging, yet so far elusive, proposition.

The Structural Evolution of Substructure

Abstract: We investigate the evolution of substructure in cold dark matter halos using N-body simulations of tidal stripping of substructure halos (subhalos) within a static host potential. We find that halos modeled following the Navarro, Frenk & White (NFW) mass profile lose mass continuously due to tides from the massive host, leading to the total disruption of satellite halos with small tidal radii. The structure of stripped NFW halos depends mainly on the fraction of mass lost, and can be expressed in terms of a simple correction to the original NFW profile. We apply these results to substructure in the Milky Way, and conclude that the dark matter halos surrounding its dwarf spheroidal (dSph) satellites have circular velocity curves that peak well beyond the luminous radius at velocities significantly higher than expected from the stellar velocity dispersion. Our modeling suggests that the true tidal radii of dSphs lie well beyond the putative tidal cutoff observed in the surface brightness profile, suggesting that the latter are not really tidal in origin but rather features in the light profile of limited dynamical relevance. For Draco, in particular, our modeling implies that its tidal radius is much larger than derived by Irwin & Hatzidimitriou (1995), lending support to the interpretation of recent Sloan survey data by Odenkirchen et al. (2001). Similarly, our model suggests that Carina's halo has a peak circular velocity of ~55 km/s, which may help explain how this small galaxy has managed to retain enough gas to undergo several bursts of star formation. Our results imply a close correspondence between the most massive subhalos expected in a CDM universe and the known satellites of the Milky Way, and suggest that only subhalos with peak circular velocities below 35 km/s lack readily detectable luminous counterparts.

The Hierarchical Origin of Galaxy Morphologies

Abstract: We report first results from a series of N-body/gasdynamical simulations designed to study the origin of galaxy morphologies in a cold dark matter-dominated universe. The simulations include star formation and feedback and have numerical resolution sufficiently high to allow for a direct investigation of the morphology of simulated galaxies. We find, in agreement with previous theoretical work, that the presence of the main morphological components of galaxies--disks, spheroids, bars--is regulated by the mode of gas accretion and intimately linked to discrete accretion events. In the case we present, disks arise from the smooth deposition of cooled gas at the center of dark halos, spheroids result from the stirring of preexisting disks during mergers, and bars are triggered by tides generated by satellites. This demonstrates that morphology is a transient phenomenon within the lifetime of a galaxy and that the Hubble sequence reflects the varied accretion histories of galaxies in hierarchical formation scenarios. In particular, we demonstrate directly that disk/bulge systems can be built and rebuilt by the smooth accretion of gas onto the remnant of a major merger and that the present-day remnants of late dissipative mergers between disks are spheroidal stellar systems with structure resembling that of field ellipticals. The perplexing variety of galaxy morphologies is thus highly suggestive of--and may actually even demand--a universe where structures have evolved hierarchically.