No Arabic abstract
We have constructed realistic, self-consistent models of triaxial elliptical galaxies embedded in triaxial dark matter halos. Self-consistent solutions by means of the standard orbital superposition technique introduced by Schwarzschild were found in each of the three cases studied. Chaotic orbits were found to be important in all of the models, and their presence was shown to imply a possible slow evolution of the shapes of the halos. The equilibrium velocity distribution is reproduced by a Lorentzian function better than by a Gaussian. Our results demonstrate for the first time that triaxial dark matter halos can co-exist with triaxial galaxies.
Galactic disks in triaxial dark matter halos become deformed by the elliptical potential in the plane of the disk in such a way as to counteract the halo ellipticity. We develop a technique to calculate the equilibrium configuration of such a disk in the combined disk-halo potential, which is based on the method of Jog (2000) but accounts for the radial variation in both the halo potential and the disk ellipticity. This crucial ingredient results in qualitatively different behavior of the disk: the disk circularizes the potential at small radii, even for a reasonably low disk mass. This effect has important implications for proposals to reconcile cuspy halo density profiles with low surface brightness galaxy rotation curves using halo triaxiality. The disk ellipticities in our models are consistent with observational estimates based on two-dimensional velocity fields and isophotal axis ratios.
Wave dark matter ($psi$DM), which satisfies the Schrodinger-Poisson equation, has recently attracted substantial attention as a possible dark matter candidate. Numerical simulations have in the past provided a powerful tool to explore this new territory of possibility. Despite their successes to reveal several key features of $psi$DM, further progress in simulations is limited, in that cosmological simulations so far can only address formation of halos below $sim 2times 10^{11} M_odot$ and substantially more massive halos have become computationally very challenging to obtain. For this reason, the present work adopts a different approach in assessing massive halos by constructing wave-halo solutions directly from the wave distribution function. This approach bears certain similarity with the analytical construction of particle-halo (cold dark matter model). Instead of many collisionless particles, one deals with one single wave that has many non-interacting eigenstates. The key ingredient in the wave-halo construction is the distribution function of the wave power, and we use several halos produced by structure formation simulations as templates to determine the wave distribution function. Among different models, we find the fermionic King model presents the best fits and we use it for our wave-halo construction. We have devised an iteration method for constructing the nonlinear halo, and demonstrate its stability by three-dimensional simulations. A Milky-Way-sized halo has also been constructed, and the inner halo is found flatter than the NFW profile. These wave-halos have small-scale interferences both in space and time producing time-dependent granules. While the spatial scale of granules varies little, the correlation time is found to increase with radius by one order of magnitude across the halo.
We construct self-consistent dynamical models for disk galaxies with triaxial, cuspy halos. We begin with an equilibrium, axisymmetric, disk-bulge-halo system and apply an artificial acceleration to the halo particles. By design, this acceleration conserves energy and thereby preserving the systems differential energy distribution even as its phase space distribution function is altered. The halo becomes triaxial but its spherically-averaged density profile remains largely unchanged. The final system is in equilibrium, to a very good approximation, so long as the halos shape changes adiabatically. The disk and bulge are ``live while the halo is being deformed; they respond to the changing gravitational potential but also influence the deformation of the halo. We test the hypothesis that halo triaxiality can explain the rotation curves of low surface brightness galaxies by modelling the galaxy F568-3.
We use the usual method of Schwarzschild to construct self-consistent solutions for the triaxial de Zeeuw & Carollo (1996) models with central density cusps. ZC96 models are triaxial generalisations of spherical $gamma$-models of Dehnen whose densities vary as $r^{-gamma}$ near the center and $r^{-4}$ at large radii and hence, possess a central density core for $gamma=0$ and cusps for $gamma > 0$. We consider four triaxial models from ZC96, two prolate triaxials: $(p, q) = (0.65, 0.60)$ with $gamma = 1.0$ and 1.5, and two oblate triaxials: $(p, q) = (0.95, 0.60)$ with $gamma = 1.0$ and 1.5. We compute 4500 orbits in each model for time periods of $10^{5} T_{D}$. We find that a large fraction of the orbits in each model are stochastic by means of their nonzero Liapunov exponents. The stochastic orbits in each model can sustain regular shapes for $sim 10^{3} T_{D}$ or longer, which suggests that they diffuse slowly through their allowed phase-space. Except for the oblate triaxial models with $gamma =1.0$, our attempts to construct self-consistent solutions employing only the regular orbits fail for the remaining three models. However, the self-consistent solutions are found to exist for all models when the stochastic and regular orbits are treated in the same way because the mixing-time, $sim10^{4} T_{D}$, is shorter than the integration time, $10^{5} T_{D}$. Moreover, the ``fully-mixed solutions can also be constructed for all models when the stochastic orbits are fully mixed at 15 lowest energy shells. Thus, we conclude that the self-consistent solutions exist for our selected prolate and oblate triaxial models with $gamma = 1.0$ and 1.5.
This talk provides a progress report on an extended collaboration which has aimed to address two basic questions, namely: Should one expect to see cuspy, triaxial galaxies in nature? And can one construct realistic cuspy, triaxial equilibrium models that are robust? Three technical results are described: (1) Unperturbed chaotic orbits in cuspy triaxial potentials can be extraordinarily sticky, much more so than orbits in many other three-dimensional potentials. (2) Even very weak perturbations can be important by drastically reducing, albeit not completely eliminating, this stickiness. (3) A simple toy model facilitates a simple understanding of why black holes and cusps can serve as an effective source of chaos. These results suggest that, when constructing models of galaxies using Schwarzschilds method or any analogue thereof, astronomers would be well advised to use orbital building blocks that have been perturbed by `noise or other weak irregularities, since such building blocks are likely to be more nearly time-independent than orbits evolved in the absence of all perturbations.