What Makes the Family of Barred Disc Galaxies So Rich: Damping Stellar Bars in Spinning Haloes

We model and analyse the secular evolution of stellar bars in spinning dark matter (DM) haloes with the cosmological spin lambda ~ 0 -- 0.09. Using high-resolution stellar and DM numerical simulations, we focus on angular momentum exchange between stellar discs and DM haloes of various axisymmetric shapes --- spherical, oblate and prolate. We find that stellar bars experience a diverse evolution which is guided by the ability of parent haloes to absorb angular momentum lost by the disc through the action of gravitational torques, resonant and non-resonant. We confirm the previous claim that dynamical bar instability is accelerated via resonant angular momentum transfer to the halo. Our main findings relate to the long-term, secular evolution of disc-halo systems: with an increasing lambda, bars experience less growth and dissolve after they pass through the vertical buckling instability. Specifically, with an increasing halo spin, (1) The vertical buckling instability in stellar bars colludes with inability of the inner halo to absorb angular momentum --- this emerges as the main factor weakening or destroying bars in spinning haloes; (2) Bars lose progressively less angular momentum, and their pattern speeds level off; (3) Bars are smaller, and for lambda >= 0.06 cease their growth completely following buckling; (4) Bars in lambda > 0.03 haloes have ratio of corotation-to-bar radii, R_CR / R_b > 2, and represent so-called slow bars which do not show offset dust lanes. We provide a quantitative analysis of angular momentum transfer in disc-halo systems, and explain the reasons for absence of growth in fast spinning haloes and its observational corollaries. We conclude that stellar bar evolution is substantially more complex than anticipated, and bars are not as resilient as has been considered so far.Comment: 15 pages., 11 figures, MNRAS, in pres

Gas Feedback on Stellar Bar Evolution

We analyze evolution of live disk-halo systems in the presence of various gas fractions, f_gas less than 8% in the disk. We addressed the issue of angular momentum (J) transfer from the gas to the bar and its effect on the bar evolution. We find that the weakening of the bar, reported in the literature, is not related to the J-exchange with the gas, but is caused by the vertical buckling instability in the gas-poor disks and by a steep heating of a stellar velocity dispersion by the central mass concentration (CMC) in the gas-rich disks. The gas has a profound effect on the onset of the buckling -- larger f_gas brings it forth due to the more massive CMCs. The former process leads to the well-known formation of the peanut-shaped bulges, while the latter results in the formation of progressively more elliptical bulges, for larger f_gas. The subsequent (secular) evolution of the bar differs -- the gas-poor models exhibit a growing bar while gas-rich models show a declining bar whose vertical swelling is driven by a secular resonance heating. The border line between the gas-poor and -rich models lies at f_gas ~ 3% in our models, but is model-dependent and will be affected by additional processes, like star formation and feedback from stellar evolution. The overall effect of the gas on the evolution of the bar is not in a direct J transfer to the stars, but in the loss of J by the gas and its influx to the center that increases the CMC. The more massive CMC damps the vertical buckling instability and depopulates orbits responsible for the appearance of peanut-shaped bulges. The action of resonant and non-resonant processes in gas-poor and gas-rich disks leads to a converging evolution in the vertical extent of the bar and its stellar dispersion velocities, and to a diverging evolution in the bulge properties.Comment: 12 pages, 12 figures, accepted for publication by the Astrophysical Journal. Minor corrections following the referee repor

Evolution of Phase-Space Density in Dark Matter Halos

The evolution of the phase-space density profile in dark matter (DM) halos is investigated by means of constrained simulations, designed to control the merging history of a given DM halo. Halos evolve through a series of quiescent phases of a slow accretion intermitted by violent events of major mergers. In the quiescent phases the density of the halo closely follows the NFW profile and the phase-space density profile, Q(r), is given by the Taylor & Navarro power law, r^{-beta}, where beta ~ 1.9 and stays remarkably stable over the Hubble time. Expressing the phase-space density by the NFW parameters, Q(r)=Qs (r/Rs)^{-beta}, the evolution of Q is determined by Qs. We have found that the effective mass surface density within Rs, Sigma_s = rhos Rs, remains constant throughout the evolution of a given DM halo along the main branch of its merging tree. This invariance entails that Qs ~ Rs^{-5/2} and Q(r) ~ Sigma_s^{-1/2} Rs^{-5/2} (r/ Rs)^{-beta}. It follows that the phase-space density remains constant, in the sense of Qs=const., in the quiescent phases and it decreases as Rs^{-5/2} in the violent ones. The physical origin of the NFW density profile and the phase-space density power law is still unknown. Yet, the numerical experiments show that halos recover these relations after the violent phases. The major mergers drive Rs to increase and Qs to decrease discontinuously while keeping Qs Rs^{5/2} = const. The virial equilibrium in the quiescent phases implies that a DM halos evolves along a sequence of NFW profiles with constant energy per unit volume (i.e., pressure) within Rs.Comment: 7 pages, 5 figures, accepted by the Astrophysical Journal. Revised, 2 figures adde

Induced Nested Galactic Bars Inside Assembling Dark Matter Halos

We investigate the formation and evolution of nested bar systems in disk galaxies in a cosmological setting by following the development of an isolated dark matter (DM) and baryon density perturbation. The disks form within the assembling triaxial DM halos and the feedback from the stellar evolution is accounted for in terms of supernovae and OB stellar winds. Focusing on a representative model, we show the formation of an oval disk and of a first generation of nested bars with characteristic sub-kpc and a few kpc sizes. The system evolves through successive dynamical couplings and decouplings, forcing the gas inwards and settles in a state of resonant coupling. The inflow rate can support a broad range of activity within the central kpc, from quasar- to Seyfert-types, supplemented by a vigorous star formation as a by-product. The initial bar formation is triggered in response to the tidal torques from the triaxial DM halo, which acts as a finite perturbation. This first generation of bars does not survive for more than 4--5 Gyr: by that time the secondary bar has totally dissolved, while the primary one has very substantially weakened, reduced to a fat oval. This evolution is largely due to chaos introduced by the interaction of the multiple non-axisymmetric components.Comment: 4 pages, 4 figures, 1 mpeg animation. To be published by the Astrophysical Journal Letters. The animation can be found at http://www.pa.uky.edu/~shlosman/research/galdyn/movies.html Replaced with an updated version (small text corrections

Disk Evolution and Bar Triggering Driven by Interactions with Dark Matter Substructure

We study formation and evolution of bar-disk systems in fully self-consistent cosmological simulations of galaxy formation in the LCDM WMAP3 Universe. In a representative model we find that the first generation of bars form in response to the asymmetric dark matter (DM) distribution (i.e., DM filament) and quickly decay. Subsequent bar generations form and are destroyed during the major merger epoch permeated by interactions with a DM substructure (subhalos). A long-lived bar is triggered by a tide from a subhalo and survives for ~10 Gyr. The evolution of this bar is followed during the subsequent numerous minor mergers and interactions with the substructure. Together with intrinsic factors, these interactions largely determine the stellar bar evolution. The bar strength and its pattern speed anticorrelate, except during interactions and when the secondary (nuclear) bar is present. For about 5 Gyr bar pattern speed increases substantially despite the loss of angular momentum to stars and cuspy DM halo. We analyze the evolution of stellar populations in the bar-disk and relate them to the underlying dynamics. While the bar is made mainly of an intermediate age, ~5-6 Gyr, disk stars at z=0, a secondary nuclear bar which surfaces at z~0.1 is made of younger, ~1-3 Gyr stars.Comment: 5 pages, 5 figures, accepted for publication in ApJ Letter

Structure Formation Inside Triaxial Dark Matter Halos: Galactic Disks, Bulges and Bars

We investigate the formation and evolution of galactic disks immersed in assembling live DM halos. Disk/halo components have been evolved from the cosmological initial conditions and represent the collapse of an isolated density perturbation. The baryons include gas (which participates in star formation [SF]) and stars. The feedback from the stellar energy release onto the ISM has been implemented. We find that (1) The growing triaxial halo figure tumbling is insignificant and the angular momentum (J) is channeled into the internal circulation; (2) Density response of the disk is out of phase with the DM, thus diluting the inner halo flatness and washing out its prolateness; (3) The total J is neathly conserved, even in models accounting for feedback; (4) The specific J for the DM is nearly constant, while that for baryons is decreasing; (5) Early stage of disk formation resembles the cat's cradle -- a small amorphous disk fueled via radial string patterns; (6) The initially puffed up gas component in the disk thins when the SF rate drops below ~5 Mo/yr; (7) About 40%-60% of the baryons remain outside the SF region; (8) Rotation curves appear to be flat and account for the observed disk/halo contributions; (9) A range of bulge-dominated to bulgeless disks was obtained; Lower density threshold for SF leads to a smaller, thicker disk; Gravitational softening in the gas has a substantial effect on various aspects of galaxy evolution and mimics a number of intrinsic processes within the ISM; (10) The models are characterized by an extensive bar-forming activity; (11) Nuclear bars, dynamically coupled and decoupled form in response to the gas inflow along the primary bars.Comment: 18 pages, 16 figures, accepted by the Astrophysical Journal. Minor revisions. The high-resolution figures can be found at http://www.pa.uky.edu/~shlosman/research/galdyn/figs07a

Dissecting Galaxy Formation: II. Comparing Substructure in Pure Dark Matter and Baryonic Models

We compare the substructure evolution in pure dark matter (DM) halos with those in the presence of baryons (PDM and BDM). The prime halos have been analyzed by Romano-Diaz et al (2009). Models have been evolved from identical initial conditions using Constrained Realizations, including star formation and feedback. A comprehensive catalog of subhalos has been compiled and properties of subhalos analyzed in the mass range of 10^8 Mo - 10^11 Mo. We find that subhalo mass functions are consistent with a single power law, M_sbh^{alpha}, but detect a nonnegligible shift between these functions, alpha -0.86 for the PDM, and -0.98 for the BDM. Overall, alpha const. in time with variations of +-15%. Second, we find that the radial mass distribution of subhalos can be approximated by a power law, R^{gamma} with a steepening around the radius of a maximal circular velocity, Rvmax, in the prime halos. Gamma ~-1.5 for the PDM and -1 for the BDM, inside Rvmax, and is steeper outside. We detect little spatial bias between the subhalo populations and the DM of the main halos. The subhalo population exhibits much less triaxiality with baryons, in tandem with the prime halo. Finally, we find that, counter-intuitively, the BDM population is depleted at a faster rate than the PDM one within the central 30kpc of the prime. Although the baryons provide a substantial glue to the subhalos, the main halos exhibit the same trend. This assures a more efficient tidal disruption of the BDM subhalos. This effect can be reversed for a more efficient feedback from stellar evolution and supermassive black holes, which will expel baryons from the center and decrease the concentration of the prime halo. We compare our results with via Lactea and Aquarius simulations and other published results.Comment: 12 pages, 9 figures, to be published by the Astrophysical Journa

Dark Matter Halos and Evolution of Bars in Disk Galaxies: Varying Gas Fraction and Gas Spatial Resolution

We conduct numerical experiments by evolving gaseous/stellar disks embedded in live dark matter halos aiming at quantifying the effect of gas spatial resolution and gas content on the bar evolution. Model sequences have been constructed using different resolution, and gas fraction has been varied along each sequence within fgas=0%-50%, keeping the disk and halo properties unchanged. We find that the spatial resolution becomes important with an increase in `fgas'. For the higher resolution model sequences, we observe a bimodal behavior in the bar evolution with respect to the gas fraction, especially during the secular phase of this evolution. The switch from the gas-poor to gas-rich behavior is abrupt and depends on the resolution used. The diverging evolution has been observed in nearly all basic parameters characterizing bars, such as the bar strength, central mass concentration, vertical buckling amplitude, size, etc. We find that the presence of the gas component severely limits the bar growth and affects its pattern speed evolution. Gas-poor models display rapidly decelerating bars, while gas-rich models exhibit bars with constant or even slowly accelerating tumbling. The gas-rich models have bar corotation (CR) radii within the disk at all times, in constrast with gas-poor and purely stellar disks. The CR-to-bar size ratio is less than 2 for gas rich-models. We have confirmed that the disk angular momentum within the CR remains unchanged in the gas-poor models, as long as the CR stays within the disk, but experiences a sharp drop before leveling off in the gas-rich models. Finally, we discuss a number of observed correlations between various parameters of simulated bars, e.g., bar sizes and gas fractions, bar strength and buckling amplitude, bar strength and its size, etc.Comment: 13 pages, 13 figures, to be published by the Astrophysical Journal; minor changes following the referee repor