No Arabic abstract
We study the behaviors of galactic disks in triaxial halos both numerically and analytically to see if warps can be excited and sustained in triaxial potentials. We consider the following two scenarios: 1) galactic disks that are initially tilted relative to the equatorial plane of the halo (for a pedagogical purpose), and 2) tilted infall of dark matter relative to the equatorial plane of the disk and the halo. With numerical simulations of 100,000 disk particles in a fixed halo potential, we find that in triaxial halos, warps can be excited and sustained just as in spherical or axisymmetric halos but they show some oscillatory behaviors and even can be transformed to a polar-ring system if the halo has a prolate-like triaxiality. The non-axisymmetric component of the halo causes the disk to nutate, and the differential nutation between the inner and outer parts of the disk generally makes the magnitude of the warp slightly diminish and fluctuate. We also find that warps are relatively weaker in oblate and oblate-like triaxial halos, and since these halos are the halo configurations of disk galaxies inferred by cosmological simulations, our results are consistent with the fact that most of the observed warps are quite weak. We derive approximate formulae for the torques exerted on the disk by the triaxial halo and the dark matter torus, and with these formulae we successfully describe the behaviors of the disks in our simulations. The techniques used in deriving these formulae could be applied for realistic halos with more complex structures.
Using spectral methods, we analyse the orbital structure of prolate/triaxial dark matter (DM) halos in N-body simulations to understand the processes that drive the evolution of shapes of DM halos and elliptical galaxies in which central masses are grown. A longstanding issue is whether the change in the shapes of DM halos is the result of chaotic scattering of box orbits, or whether they change shape adiabatically in response to the evolving galactic potential. We use orbital frequencies to classify orbits, to quantify orbital shapes, and to identify resonant orbits and chaotic orbits. The frequency-based method overcomes the limitations of Lyapunov exponents which are sensitive to numerical discreteness effects. Regardless of the distribution of the baryonic component, the shape of a DM halo changes primarily due to changes in the shapes of individual orbits within a given family. Orbits with small pericentric radii are more likely to change both their orbital type and shape than orbits with large pericentric radii. Whether the evolution is regular (and reversible) or chaotic (and irreversible), depends primarily on the radial distribution of the baryonic component. The growth of an extended baryonic component of any shape results in a regular rather than chaotic change in orbital populations. In contrast the growth of a massive and compact central component results in chaotic scattering of a significant fraction of both box and long-axis tube orbits. The growth of a disk causes a significant fraction of halo particles to become trapped by major global orbital resonances. Despite the fact that shape of a DM halo is always quite oblate following the growth of a central baryonic component, a significant fraction of its orbit population has characteristics of its triaxial or prolate progenitor (ABRIDGED).
We analyze warps in the nearby edge-on spiral galaxies observed in the {em Spitzer/IRAC} 4.5 micron band. In our sample of 24 galaxies we find evidence of warp in 14 galaxies. We estimate the observed onset radii for the warps in a subsample of 10 galaxies. The dark matter distribution in each of these galaxies are calculated using the mass distribution derived from the observed light distribution and the observed rotation curves. The theoretical predictions of the onset radii for the warps are then derived by applying a self-consistent linear response theory to the obtained mass models for 6 galaxies with rotation curves in the literature. By comparing the observed onset radii to the theoretical ones, we find that discs with constant thickness can not explain the observations; moderately flaring discs are needed. The required flaring is consistent with the observations. Our analysis shows that the onset of warp is not symmetric in our sample of galaxies. We define a new quantity called the onset-asymmetry index and study its dependence on galaxy properties. The onset asymmetries in warps tend to be larger in galaxies with smaller disc scale lengths. We also define and quantify the global asymmetry in the stellar light distribution, that we call the edge-on asymmetry in edge-on galaxies. It is shown that in most cases the onset asymmetry in warp is actually anti-correlated with the measured edge-on asymmetry in our sample of edge-on galaxies and this could plausibly indicate that the surrounding dark matter distribution is asymmetric.
The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MWs stellar halo claim a detection of an isotropic or tangential dip at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of dips in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halos mass when using spherical Jeans equation modeling.
We use N-body hydrodynamical simulations to study the structure of disks in triaxial potentials resembling CDM halos. Our analysis focuses on the accuracy of the dark mass distribution inferred from rotation curves derived from simulated long-slit spectra. We consider a massless disk embedded in a halo with axis ratios of 0.5:0.6:1.0 and with its rotation axis aligned with the minor axis of the halo. Closed orbits for the gaseous particles deviate from coplanar circular symmetry, resulting in a variety of long-slit rotation curve shapes, depending on the orientation of the disk relative to the line of sight. Rotation curves may thus differ significantly from the spherically-averaged circular velocity profile of the dark matter halo. Solid-body rotation curves--typically interpreted as a signature of a constant density core in the dark matter distribution--are obtained about 25% of the time for random orientations although the dark matter follows the cuspy density profile proposed by Navarro, Frenk & White (NFW). We conclude that the discrepancies reported between the shape of the rotation curve of low surface brightness galaxies and the structure of CDM halos may be resolved once the complex effects of halo triaxiality on the dynamics of the gas component is properly taken into account.
We have analyzed high resolution N-body simulations of dark matter halos, focusing specifically on the evolution of angular momentum. We find that not only is individual particle angular momentum not conserved, but the angular momentum of radial shells also varies over the age of the Universe by up to factors of a few. We find that torques from external structure are the most likely cause for this distribution shift. Since the model of adiabatic contraction that is often applied to model the effects of galaxy evolution on the dark-matter density profile in a halo assumes angular momentum conservation, this variation implies that there is a fundamental limit on the possible accuracy of the adiabatic contraction model in modeling the response of DM halos to the growth of galaxies.