Disk galaxies can exchange angular momentum and baryons with their host dark matter (DM) halos. These halos possess internal spin, `lambda', which is insignificant rotationally but does affect interactions between the baryonic and DM components. While statistics of prograde and retrograde spinning halos in galaxies is not available at present, the existence of such halos is important for galaxy evolution. In the previous works, we analyzed dynamical and secular evolution of stellar bars in prograde spinning halos and the DM response to the bar perturbation, and found that it is modified by the resonant interactions between the bar and the DM halo orbits. In the present work, we follow the evolution of stellar bars in retrograde halos. We find, that this evolution differs substantially from evolution in rigid unresponsive halos, discussed in the literature. First, we confirm that the bar instability is delayed progressively along the retrograde `lambda' sequence. Second, the bar evolution in the retrograde halos differs also from that in the prograde halos, in that the bars continue to grow substantially over the simulation time of 10 Gyr. The DM response is also substantially weaker compared to this response in the prograde halos. Third, using orbital spectral analysis of the DM orbital structure, we find a phenomenon we call the orbit reversal --- when retrograde DM orbits interact with the stellar bar, reverse their streaming and precession, and become prograde. This process dominates the inner halo region adjacent to the bar and allows these orbits to be trapped by the bar, thus increasing efficiency of angular momentum transfer by the Inner Lindblad Resonance. We demonstrate this reversal process explicitly in a number of examples.

We study nonlinear response of spinning dark matter (DM) halos to dynamic and secular evolution of stellar bars in the embedded galactic disks, using high-resolution numerical simulations. For a sequence of halos with the cosmological spin parameter lambda=0-0.09, and a representative angular momentum distribution, we analyze evolution of induced DM bars amplitude and quantify parameters of the response as well as trapping of DM orbits and angular momentum transfer by the main and secondary resonances. We find that (1) maximal amplitude of DM bars depends strongly on `lambda', while that of the stellar bars is indifferent to `lambda'; (2) Efficiency of resonance trapping of DM orbits by the bar increases with `lambda', and so is the mass and the volume of DM bars; (3) Contribution of resonance transfer of angular momentum to the DM halo increases with `lambda', and for larger spin, the DM halo `talks' to itself, by moving the angular momentum to larger radii --- this process is maintained by resonances; (4) Prograde and retrograde DM orbits play different roles in angular momentum transfer. The `active' part of the halo extends well beyond the bar region, up to few times the bar length in equatorial plane and away from this plane. (5) We model evolution of diskless DM halos and halos with frozen disks, and found them to be perfectly stable to any Fourier modes. Finally, further studies adopting a range of mass and specific angular momentum distributions of the DM halo will generalize the dependence of DM response on the halo spin and important implications for direct detection of DM and that of the associated stellar tracers, such as streamers.

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.

We demonstrate that growth of stellar bars in spinning dark matter halos is heavily suppressed in the secular phase of evolution, using numerical simulations of isolated galaxies. In a representative set of models, we show that for values of the cosmological spin parameter lambda > 0.03, bar growth (in strength and size) becomes increasingly quenched. Furthermore, slowdown of bar pattern speed weakens substantially with increasing `lambda,' until it ceases completely. The terminal structure of bars is affected as well, including extent and shape of their boxy/peanut bulges. The essence of this effect lies in the modified angular momentum exchange between the disk and the halo facilitated by the bar. For the first time we have demonstrated that a dark matter halo can emit and not purely absorb angular momentum. Although the halo as a whole is not found to emit, the net transfer of angular momentum from the disk to the halo is significantly reduced or completely eliminated. The paradigm shift implies that the accepted view that disks serve as sources of angular momentum and halos serve as sinks, must be revised. Halos with lambda > 0.03 are expected to form a substantial fraction, based on lognormal distribution of lambda. Dependence of secular bar evolution on halo spin, therefore, implies profound corollaries for the cosmological evolution of galactic disks.

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.

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.

We construct and evolve families of steady-state models of stellar disks embedded in live DM halos, in order to study the dynamical and secular phases of bar evolution. These models are tested against those published in the literature in order to extend them and include the gaseous component in the follow up paper. We are interested in the angular momentum (J) redistribution in the disk-halo system. We confirm the previous results and quantify for the first time the dual role that the DM halos play in the bar evolution: more centrally concentrated halos dilute the dynamical processes, such as spontaneous bar instability and vertical buckling instability, and slowdown the J transfer, while facilitating it in the secular phase. Within the corotation radius (Rcr), the disk J remains nearly constant, as long as Rcr stays within the disk -- a sign that the lost J to the outer disk and the halo is being compensated by an influx of fresh J due to the outward motion of Rcr. This is feasible as long as the bar slowdown dominates the loss of J inside Rcr. We find that in some models the bar pattern speed stalls for prolonged time periods when Rcr is located outside the disk. This phenomenon appears concurrent with the near absence of J transfer between the disk and the halo. Furthermore, we confirm that stellar bars generally display the corotation to bar size ratios in the range of ~1-1.4, but only between the times of the first buckling and Rcr leaving the disk. The corotation-to-disk size ratio emerges as an important dynamic discriminator between various stages of barred disk evolution. Finally, we analyze a number of correlations between the basic parameters of a barred disk and a halo, some already reported in the literature and some new.

We compare assembly of DM halos with and without baryons, within the context of cosmological evolution in the LCDM WMAP3 Universe (baryons+DM, BDM model, and pure DM, PDM model). In representative PDM and BDM models we find that baryons contribute decisively to the evolution of the central region, leading to an isothermal DM cusp, and to a flat DM density core -- the result of heating by dynamical friction of the substructure during a quiescent evolution epoch. This process ablates the cold gas from an embedded disk, cutting the star formation rate by ~10, and heats up the spheroidal gas and stellar components, triggering their expansion. The substructure is more resilient in the presence of baryons. The disk which formed from inside-out as gas dominated, is transformed into an intermediate Hubble type by z ~ 2 and to an early type by z ~ 0.5, based on its gas contents and spheroidal-to-disk stellar mass ratio. Only a relatively small ~20% fraction of DM particles in PDM and BDM models are bound within the radius of maximal circular velocity in the halo -- most of the DM particles perform larger radial excursions. We also find that the fraction of baryons within the halo virial radius somewhat increases during the major mergers and decreases during the minor mergers. The net effect appears to be negligible. While the substructure is being tidally-disrupted, mixing of its debris in the halo is not efficient and becomes even less so with z. The streamers formed after z ~ 1 survive largely to the present time -- an important implication for embedded disk evolution.

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.

We study the central dark matter (DM) cusp evolution in cosmological galactic halos. Models with and without baryons (baryons+DM, hereafter BDM model, and pure DM, PDM model, respectively) are advanced from identical initial conditions. The DM cusp properties are contrasted by a direct comparison of pure DM and baryonic models. We find a divergent evolution between the PDM and BDM models within the inner ~10 kpc region. The PDM model forms a R^{-1} cusp as expected, while the DM in the BDM model forms a larger isothermal cusp R^{-2} instead. The isothermal cusp is stable until z~1 when it gradually levels off. This leveling proceeds from inside out and the final density slope is shallower than -1 within the central 3 kpc (i.e., expected size of the R^{-1} cusp), tending to a flat core within ~2 kpc. This effect cannot be explained by a finite resolution of our code which produces only a 5% difference between the gravitationally softened force and the exact Newtonian force of point masses at 1 kpc from the center. Neither is it related to the energy feedback from stellar evolution or angular momentum transfer from the bar. Instead it can be associated with the action of DM+baryon subhalos heating up the cusp region via dynamical friction and forcing the DM in the cusp to flow out and to `cool' down. The process described here is not limited to low z and can be efficient at intermediate and even high z.

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.

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.

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.

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.

We investigate the effect of an assembly history on the evolution of dark matter (DM) halos of 10^{12} Msun/h using Constrained Realizations of random Gaussian fields. Five different realizations of a DM halo with distinct merging histories were constructed and evolved. Our main results are: A halo evolves via a sequence of quiescent phases of a slow mass accretion intermitted by violent episodes of major mergers. In the quiescent phases, the density is well fitted by an NFW profile, the inner scale radius Rs and the mass enclosed within it remain constant, and the virial radius (Rvir) grows linearly with the expansion parameter "a". Within each quiescent phase the concentration parameter ("c") scales as "a", and the mass accretion history (Mvir) is well described by the Tasitsiomi etal. fitting formula. In the violent phases the halos are not in a virial equilibrium and both Rs and Rvir grow discontinuously. The violent episodes drive the halos from one NFW dynamical equilibrium to another. The final structure of a halo, including "c", depends on the degree of violence of the major mergers and on their number. Next, we find a distinct difference between the behavior of various NFW parameters taken as averages over an ensemble of halos and those of individual halos. Moreover, the simple scaling relations c--Mvir do not apply to the entire evolution of individual halos, and so is the common notion that late forming halos are less concentrated than early forming ones. The entire evolution of the halo cannot be fitted by single analytical expressions.

The formation and structure of dark matter halos is studied by constrained simulations. A series of experiments of the formation of a 10^12 Msun/h halo is designed to study the dependence of the density profile on its merging history. We find that the halo growth consist of several quiescent phases intermitted by violent events, with the density well approximated by the NFW profile during the former phases. We find that (1) the NFW scale radius R_s stays constant during the quiescent phase and grows abruptly during the violent one. In contrast, the virial radius grows linearly during the quiescent and abruptly during the violent phases. (2) The central density stays unchanged during the quiescent phase while dropping abruptly during the violent phase, and it does not reflect the formation time of the halo. (3) The clear separation of the evolution of an individual halo into quiescent and violent phases implies that its entire evolution cannot be fitted by simple scaling relations.

Evolution of stellar bars in disk galaxies is accompanied by dynamical instabilities and secular changes. Following the vertical buckling instability, the bars are known to weaken dramatically and develop a pronounced boxy/peanut shape when observed edge-on. Using high-resolution N-body simulations of stellar disks embedded in live axisymmetric dark matter halos, we have investigated the long-term changes in the bar morphology, specifically the evolution of the bar size, its vertical structure and exchange of angular momentum. We find that following the initial buckling, the bar resumes its growth from deep inside the corotation radius and follows the Ultra-Harmonic resonance thereafter. We also find that this secular bar growth triggers a spectacular secondary vertical buckling instability which leads to the appearance of characteristic boxy/peanut/X-shaped bulges. The secular bar growth is crucial for the recurrent buckling to develop. Furthermore, the secondary buckling is milder, persists over a substantial period of time, ~3 Gyr, and can have observational counterparts. Overall, the stellar bars show recurrent behavior in their properties and evolve by increasing their linear and vertical extents, both dynamically and secularly. We also demonstrate explicitly that the prolonged growth of the bar is mediated by continuous angular momentum transfer from the disk to the surrounding halo, and that this angular momentum redistribution is resonant in nature -- a large number of lower resonances trap disk and halo particles and this trapping is robust, in a broad agreement with the earlier results in the literature.

The formation and structure of dark matter (DM) halos is studied by means of constrained realizations of Gaussian fields using N-body simulations. A series of experiments of the formation of a 10^{12} Msun halo is designed to study the dependence of the density profile on its merging history. We confirm that the halo growth consists of violent and quiescent phases, with the density well approximated by the Navarro-Frenk-White (NFW) profile during the latter phases. We find that (1) the NFW scale radius R_s stays constant during the quiescent phase and grows abruptly during the violent one. In contrast, the virial radius grows linearly during the quiescent and abruptly during the violent phases. (2) The central density stays unchanged during the quiescent phase while dropping abruptly during the violent phase. (3) The value of \rs reflects the violent merging history of the halo, and depends on the number of violent events and their fractional magnitudes, independent of the time and order of these events. It does not reflect the formation time of the halo. (4) The fractional change in R_s is a nonlinear function of the fractional absorbed kinetic energy within R_s in a violent event.

We study the regeneration of stellar bars triggered by a tidal interaction, using numerical simulations of either purely stellar or stellar+gas disc galaxies. We find that interactions which are sufficiently strong to regenerate the bar in the purely stellar models do not lead to a regeneration in the dissipative models, owing to the induced gas inflow in those models. In models in which the bar can be regenerated, we find a tight correlation between the strength and the pattern speed of the induced bar. This relation can be explained by a significant radial redistribution of angular momentum in the disc due to the interaction, similar to the processes and correlations found for isolated barred spirals. We furthermore show that the regenerated bars show the same dynamical properties as their isolated counterparts.

We investigate the dynamical effects of an interaction between an initially barred galaxy and a small spherical companion using an N-body/SPH algorithm. In the models described here the small companion passes through the disc of the larger galaxy near-perpendicular to its plane. The impact positions and times are varied with respect to the phase of the bar and the dynamical evolution of the disc. The interactions produce expanding ring structures, offset bars, spokes,and other asymmetries in the stars and gas. These characteristic signatures of the interaction are present in the disc for about 1 Gyr. We find that in some cases it is possible to destroy the bar while keeping the disc structure. In general the central impacts cause larger damage to the bar and the disc than the peripheral ones. The interaction tends to accelerate the transition from a strongly barred galaxy to a weakly or non-barred galaxy. The final disc morphology is determined more by the impact position relative to the bar rather than the impact time.

Under certain conditions, sub-kpc nuclear bars form inside large-scale stellar bars of disk galaxies. These secondary bars spend a fraction of their lifetime in a dynamically-decoupled state, tumbling in the gravitational field of the outer bars. We analyze the flow pattern in such nested bar systems and find that secondary bars differ fundamentally from their large-scale counterparts. In particular the gas flow across the bar-bar interface in these systems can be more chaotic or more regular in nature, and, contrary to predictions, has no difficulty in penetrating the secondary bars. The outer parts of both short and long nuclear bars (with respect to their corotation) appear to be depopulated of gas, while deep inside them the flow exhibits low Mach numbers and follows ovally-shaped orbits with little dissipation. We find that gas-dominated and star-dominated nuclear bars avoid the bar-bar interface, making both types of bars short relative to their corotation. Furthermore, our earlier work has shown that dynamically-coupled secondary bars exhibit a similarly relaxed low-dissipation flow as well. Therefore, no large-scale shocks form in the nuclear bars, and consequently, no offset dust lanes are expected there. We find that offset dust lanes cannot be used in the search for secondary (nuclear) bars.

We find that nuclear rings in barred galaxies can be subject to a new type of non-self-gravitational dynamical instability. The instability leads to the formation of gaseous molecular bars with pattern speeds which are substantially slower than speeds of the primary stellar bars. This spectacular decoupling of nuclear bars from the underlying gravitational potential is triggered but is not driven by the gas viscosity. We find that low-viscosity systems can spend a substantial period of time in a fully decoupled state, with the nuclear bar slowly tumbling in the gravitational field of the primary bar. Higher viscosity systems form nuclear bars which librate about the primary bar. The shape of a nuclear bar, i.e., its eccentricity, correlates strongly with the angle between the bars. We also find that such decoupling, partial or full, most probably will be associated with bursts of star formation and with gas inflow across the inner (Lindblad) resonance zone towards smaller radii.

Two continuum subtracted H-alpha HST frames of M101 are used to determine the positions, angular sizes and absolute fluxes of 237 HII regions using a semi-automated technique. From these we have constructed the luminosity and diameter distribution functions. We repeat this process on the images after artificially reducing the linear resolution to that typically obtained with ground based imaging. We find substantial differences in the luminosity function and diameter distribution. The measured internal properties, such as central surface brightness and radial gradient are dominated by the PSF at linear resolutions less than roughly 40 pc FWHM. From the ground such resolutions are currently only obtainable for the nearest galaxies. We find evidence for two regimes of clustering of the HII regions and diffuse emission, suggesting two different regimes of star formation in late type spiral galaxies. This property in combination with the blending that occurs at ground based resolutions might be responsible for the break in the luminosity function which is usually found at log L=38.6 erg/s. From the high resolution HST data we find a luminosity function slope of alpha=-1.74+/-0.08. We also observe a flattening at luminosities log L<36.7 erg/s. For the diameter distribution we find a characteristic scale of D0=29.2 pc from an exponential fit. However, a scale free power law with index beta=-2.84+/-0.16 provides a better fit to the data.

We present high angular and velocity resolution two-dimensional kinematic observations in the spectral lines of H alpha and CO J=1-0 of the circumnuclear starburst region in the barred spiral galaxy M100, and compare them with kinematics derived from our previously published numerical modeling. The H alpha data, fully sampled and at sub-arcsecond resolution, show a rotation curve that is rapidly rising in the central ~140 pc, and stays roughly constant, at the main disk value, further out. Non-circular motions are studied from the H alpha and CO data by detailed consideration of the velocity fields, residual velocity fields after subtraction of the rotation curve, and sets of position-velocity diagrams. These motions are interpreted as the kinematic signatures of gas streaming along the inner part of the bar, and of density wave streaming motions across a two-armed mini-spiral. Comparison with a two-dimensional velocity field and rotation curve derived from our 1995 dynamical model shows good qualitative and quantitative agreement for the circular and non-circular kinematic components. Both morphology and kinematics of this region require the presence of a double inner Lindblad Resonance in order to explain the observed twisting of the near-infrared isophotes and the gas velocity field. These are compatible with the presence of a global density wave driven by the moderately strong stellar bar in this galaxy. We review recent observational and modeling results on the circumnuclear region in M100, and discuss the implications for bar structure and gas dynamics in the core of M100 and other disk galaxies.

Using an N-body+SPH code we have performed numerical simulations to investigate the dynamical effects of an interaction between an initially barred galaxy and a small spherical companion. In the models described here the small companion passes through the disk of the larger galaxy perpendicular to its plane. The impact positions and times are varied with respect to the evolutionary phase of the bar and disc. The interactions produce expanding ring structures, offset bars, spokes, and other asymmetries in the stars and gas. They also affect the strength and pattern speed of the bar.

The evolution in shape of an isolated density enhancement in the early universe is studied through numerical simulations. The formation scenarios of a cold dissipationless collapse and that of a slow accumulation of gas in a dark matter halo are examined. In the later case the conversion of gas to stars and the accompanying energy feedback from stellar winds and SNII are taken into account. It is found that with the more realistic initial conditions the radial orbit instability (ROI) is able to operate at higher values of the virial coefficient than previously believed. In the cold dissipationless collapse the system becomes fully axisymmetric on a timescale that is less than one collapse time. A prolate or oblate figure is obtained for low and high rotations, respectively. The inclusion of gas and star formation leads to halos that are more triaxial, while the stellar systems that form in them are oblate, independent of the initial rotation.

We carry out a detailed orbit analysis of gravitational potentials selected at different times from an evolving self-consistent model galaxy consisting of a two-component disk (stars+gas) and a live halo. The results are compared with a pure stellar model, subject to nearly identical initial conditions, which are chosen as to make the models develop a large scale stellar bar. The bars are also subject to hose-pipe (buckling) instability which modifies the vertical structure of the disk. The diverging morphological evolution of both models is explained in terms of gas radial inflow, the resulting change in the gravitational potential at smaller radii, and the subsequent modification of the main families of orbits, both in and out of the disk plane. We find that dynamical instabilities become milder in the presence of the gas component, and that the stability of planar and 3D stellar orbits is strongly affected by the related changes in the potential -- both are destabilized with the gas accumulation at the center. This is reflected in the overall lower amplitude of the bar mode and in the substantial weakening of the bar, which appears to be a gradual process. The vertical buckling of the bar is much less pronounced and the characteristic peanut shape of the galactic bulge almost disappears when there is a substantial gas inflow towards the center. Milder instability results in a smaller bulge whose basic parameters are in agreement with observations. We also find that the overall evolution in the model with a gas component is accelerated due to the larger central mass concentration and resulting decrease in the characteristic dynamical time.

We investigate the dynamical response of stellar orbits in a rotating barred galaxy potential to the perturbation by a nuclear gaseous ring. The change in 3D periodic orbit families is examined as the gas accumulates between the inner Lindblad resonances. It is found that the phase space allowable to the x2 family of orbits is substantially increased and a vertical instability strip appears with the growing mass of the ring. A significant distortion of the x1 orbits is observed in the vicinity of the ring, which leads to the intersection between orbits with different values of the Jacobi integral. We also examine the dependence of the orbital response to the eccentricity and alignment of the ring with the bar. Misalignment between an oval ring and a bar can leave observational footprints in the form of twisted near-infrared isophotes in the vicinity of the ring. It is inferred that a massive nuclear ring acts to weaken and dissolve the stellar bar exterior to the ring, whereas only weakly affecting the orbits interior to the inner Lindblad resonances. Consequences for gas evolution in the circumnuclear regions of barred galaxies are discussed as well.

High-resolution NIR and optical images are used to constrain a dynamical model of the circumnuclear star forming (SF) region in the barred galaxy M100 (=NGC 4321). Subarcsecond resolution allowed us to distinguish important morphological details which are easily misinterpreted when using images at lower resolution. Small leading arms observed in our K-band image of the nuclear region are reproduced in the gas flow in our model, and lead us to believe that part of the K light comes from young stars, which trace the gas flow.

We investigate the dynamical response of stellar orbits in a rotating barred galaxy potential to the perturbation by a nuclear gaseous ring. The change in 3D periodic orbit families is examined as the gas accumulates near the inner Lindblad resonance. It is found that the x2/x3 loop extends to higher Jacobi energy and a vertical instability strip forms in each family. These strips are connected by a symmetric/anti-symmetric pair of 2:2:1 3D orbital families. A significant distortion of the x1 orbits is observed in the vicinity of the ring, which leads to the intersection between orbits over a large range of the Jacobi integral. We also find that a moderately elliptical ring oblique to the stellar bar produces significant phase shifts in the x1 orbital response.

We present new high-resolution observations of the center of the late-type spiral M100 (NGC 4321) supplemented by 3D numerical modeling of stellar and gas dynamics, including star formation (SF). NIR imaging has revealed a stellar bar, previously inferred from optical and 21 cm observations, and an ovally-shaped ring-like structure in the plane of the disk. The K isophotes become progressively elongated and skewed to the position angle of the bar (outside and inside the `ring') forming an inner bar-like region. The galaxy exhibits a circumnuclear starburst in the inner part of the K `ring'. Two maxima of the K emission have been observed to lie symmetrically with respect to the nucleus and equidistant from it slightly leading the stellar bar. We interpret the twists in the K isophotes as being indicative of the presence of a double inner Lindblad resonance (ILR) and test this hypothesis by modeling the gas flow in a self-consistent gas + stars disk embedded in a halo, with an overall NGC4321-like mass distribution. We have reproduced the basic morphology of the region (the bar, the large scale trailing shocks, two symmetric K peaks corresponding to gas compression maxima which lie at the caustic formed by the interaction of a pair of trailing and leading shocks in the vicinity of the inner ILR, both peaks being sites of SF, and two additional zones of SF corresponding to the gas compression maxima, referred usually as `twin peaks').

New optical and NIR K-band images of the inner 3 kpc region of the nearby Virgo spiral M100 (NGC 4321) display remarkable morphological changes with wavelength. While in the optical the light is dominated by a circumnuclear zone of enhanced star formation, the features in the 2.2\mum image correspond to a newly discovered kpc-size stellar bar, and a pair of leading arms situated inside an ovally-shaped region. Only 5% of the K flux is emitted in antisymmetric structures. This indicates that the morphology seen in the NIR is dominated by a global density wave. Making a first-order correction for the presence of localized dust extinction in K using the I-K image, we find that the observed leading arm morphology is slightly hidden by dust. Possible mechanisms responsible for the optical and NIR morphology are discussed, and tests are proposed to discriminate between them. Conclusions are supported with an evolutionary stellar population model reproducing the optical and NIR colors in a number of star forming zones. We argue that this morphology is compatible with the presence of a pair of inner Lindblad resonances in the region, and show this explicitly in an accompanying paper by detailed numerical modeling. The observed phenomena may provide insight into physical processes leading to central activity in galaxies.