Observations of high-redshift quasars imply the presence of supermassive black holes already at z~ 7.5. An appealing and promising pathway to their formation is the direct collapse scenario of a primordial gas in atomic-cooling haloes at z ~ 10 - 20, when the $\rm H_2$ formation is inhibited by a strong background radiation field, whose intensity exceeds a critical value, $J_{\rm crit}$. To estimate $J_{\rm crit}$, typically, studies have assumed idealized spectra, with a fixed ratio of $\rm H_{2}$ photo-dissociation rate $k_{\rm H_2}$ to the $\rm H^-$ photo-detachment rate $k_{\rm H^-}$. This assumption, however, could be too narrow in scope as the nature of the background radiation field is not known precisely. In this work we argue that the critical condition for suppressing the $\rm H_2$ cooling in the collapsing gas could be described in a more general way by a combination of $k_{\rm H_2}$ and $k_{\rm H^-}$ parameters. By performing a series of cosmological zoom-in simulations with an encompassing set of $k_{\rm H_2}$ and $k_{\rm H^-}$, we examine the gas flow by following evolution of basic parameters of the accretion flow. We test under what conditions the gas evolution is dominated by $\rm H_{2}$ and/or atomic cooling. We confirm the existence of a critical curve in the $k_{\rm H_2}-k_{\rm H^-}$ plane, and provide an analytical fit to it. This curve depends on the conditions in the direct collapse, and reveals domains where the atomic cooling dominates over the molecular cooling. Furthermore, we have considered the effect of $\rm H_{2}$ self-shielding on the critical curve, by adopting three methods for the effective column density approximation in $\rm H_{2}$. We find that the estimate of the characteristic length-scale for shielding can be improved by using $\lambda_{\rm Jeans25}$, which is 0.25 times that of the local Jeans length.

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 perform Monte-Carlo radiative transfer calculations to model the Lyman alpha properties of galaxies in high-resolution, zoom-in cosmological simulations at z ~ 6.6. The simulations include both constrained and unconstrained runs, representing respectively a highly overdense region and an average field. Different galactic wind models are used in the simulations in order to investigate the effects of these winds on the apparent Ly alpha properties of galaxies. We find that, for models including galactic winds, the Ly alpha properties of massive galaxies residing in the overdense region match well recent observations of luminous Ly alpha emitters (LAEs) at z ~ 6-7, in terms of apparent Ly alpha luminosity, Ly alpha line width and Ly alpha equivalent width distributions. Without winds, the same galaxies appear less Ly alpha bright as a result of both differences in the line profile emerging from galaxies themselves, and, in the distributions of neutral gas in the circumgalactic (CGM) and intergalactic medium (IGM). We also study the relations between apparent Ly alpha luminosity and various galaxy properties: stellar mass, star formation rate (SFR) and host halo mass. At fixed halo mass, the apparent Ly alpha luminosity of galaxies appears to depend on the large-scale environment while this is no longer true for galaxies at a given stellar mass or SFR. We provide simple linear fits to these relations that can be used for quickly constructing mock LAE samples from N-body simulations. Our results suggest that the observed luminous LAEs at z ~ 6.6 are hosted by ~10^{12} h^{-1} Mo, dark matter haloes, residing in large, overdense ionized regions.

We have modeled direct collapse of a primordial gas within dark matter halos in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disk, driving a pair of spiral shocks, subject to Kelvin-Helmholtz shear instability forming fragments; (b) the growing central core in the FLD flow quickly reaches ~10 Mo and a highly variable luminosity of 10^{38}-10^{39} erg/s, comparable to the Eddington luminosity. It experiences massive recurrent outflows driven by radiation force and thermal pressure gradients, which mix with the accretion flow and transfer the angular momentum outwards; and (c) the interplay between these processes and a massive accretion, results in photosphere at ~10 AU. We conclude that in the FLD model (1) the central object exhibits dynamically insignificant rotation and slower than adiabatic temperature rise with density; (2) does not experience fragmentation leading to star formation, thus promoting the fast track formation of a supermassive black hole (SMBH) seed; (3) inclusion of radiation force leads to outflows, resulting in the mass accumulation within the central 10^{-3} pc, which is ~100 times larger than characteristic scale of star formation. The inclusion of radiative transfer reveals complex early stages of formation and growth of the central structure in the direct collapse scenario of SMBH seed formation.

Direct collapse within dark matter (DM) halos is a promising path to form supermassive black hole (SMBH) seeds at high redshifts. The outer part of this collapse remains optically thin, and has been studied intensively using numerical simulations. However, the innermost region of the collapse is expected to become optically thick and requires us to follow the radiation field in order to understand its subsequent evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation, for isolated halos. For direct collapse within isolated DM halos, we find that (1) the photosphere forms at ~10^{-6} pc and rapidly expands outward. (2) A central core forms, with a mass of ~1 Mo, supported by thermal gas pressure gradients and rotation. (3) Growing thermal gas and radiation pressure gradients dissolve it. (4) This process is associated with a strong anisotropic outflow, and another core forms nearby and grows rapidly. (5) Typical radiation luminosity emerging from the photosphere encompassing these cores is ~5 x 10^{37}-5 x 10^{38} erg/s, of order the Eddington luminosity. (6) Two variability timescales are associated with this process: a long one, which is related to the accretion flow within the central ~10^{-4}-10^{-3} pc, and ~0.1 yr, which is related to radiation diffusion. (7) Adiabatic models have been run for comparison and their evolution differs profoundly from that of the FLD models, by forming a central geometrically-thick disk. Overall, an adiabatic equation of state is not a good approximation to the advanced stage of direct collapse, mainly because the radiation in the FLD is capable of escaping due to anisotropy in the optical depth and associated gradients.

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 employ high-resolution cosmological zoom-in simulations focusing on a high-sigma peak and an average cosmological field at $z\sim 6-12$, in order to investigate the influence of environment and baryonic feedback on galaxy evolution in the reionization epoch. Strong feedback, e.g., galactic winds, caused by elevated star formation rates (SFRs) is expected to play an important role in this evolution. We compare different outflow prescriptions: (i) constant wind velocity (CW), (ii) variable wind scaling with galaxy properties (VW), and (iii) no outflows (NW). The overdensity leads to accelerated evolution of dark matter and baryonic structures, absent in the "normal" region, and to shallow galaxy stellar mass functions at the low-mass end. Although CW shows little dependence on both environments, the more physically motivated VW model does exhibit this effect. In addition, VW can reproduce the observed specific SFR (sSFR) and the sSFR-stellar mass relation, which CW and NW fail to satisfy simultaneously. Winds also differ substantially in affecting the state of the intergalactic medium (IGM). The difference lies in volume-filling factor of hot, high-metallicity gas which is near unity for CW, while it remains confined in massive filaments for VW, and locked up in galaxies for NW. Such gas is nearly absent in the normal region. Although all wind models suffer from deficiencies, the VW model seems to be promising in correlating the outflow properties to those of host galaxies. Further constraints on the state of the IGM at high-$z$ are needed to separate different wind models.

We provide detailed comparison between the AMR code Enzo-2.4 and the SPH/N- body code GADGET-3 in the context of isolated or cosmological direct baryonic collapse within dark matter (DM) halos to form supermassive black holes. Gas flow is examined by following evolution of basic parameters of accretion flows. Both codes show an overall agreement in the general features of the collapse, however, many subtle differences exist. For isolated models, the codes increase their spatial and mass resolutions at different pace, which leads to substantially earlier collapse in SPH than in AMR cases due to higher gravitational resolution in GADGET-3. In cosmological runs, the AMR develops a slightly higher baryonic resolution than SPH during halo growth via cold accretion permeated by mergers. Still, both codes agree in the buildup of DM and baryonic structures. However, with the onset of collapse, this difference in mass and spatial resolution is amplified, so evolution of SPH models begins to lag behind. Such a delay can have effect on formation/destruction rate of H2 due to UV background, and on basic properties of host halos. Finally, isolated non-cosmological models in spinning halos, with spin parameter {\lambda} ~ 0.01 - 0.07, show delayed collapse for greater {\lambda}, but pace of this increase is faster for AMR. Within our simulation setup, GADGET-3 requires significantly larger computational resources than Enzo- 2.4 during collapse, and needs similar resources, during the pre-collapse, cosmological structure formation phase. Yet it benefits from substantially higher gravitational force and hydrodynamic resolutions, except at the end of collapse.

We investigate the ability of jets in active galactic nuclei (AGNs) to break out of the ambient gas with sufficiently large advance velocities. Using observationally estimated jet power, we analyze 28 bright elliptical galaxies in nearby galaxy clusters. Because the gas density profiles in the innermost regions of galaxies have not been resolved so far, we consider two extreme cases for temperature and density profiles. We also follow two types of evolution for the jet cocoons: being driven by the pressure inside the cocoon (Fanaroff-Riley [FR] I type), and being driven by the jet momentum (FR II type). Our main result is that regardless of the assumed form of density profiles, jets with observed powers of <~ 10^{44} erg s^{-1} are not powerful enough to evolve as FR II sources. Instead, they evolve as FR I sources and appear to be decelerated below the buoyant velocities of the cocoons when jets were propagating through the central dense regions of the host galaxies. This explains the reason why FR I sources are more frequent than FR II sources in clusters. Furthermore, we predict the sizes of X-ray cavities from the observed jet powers and compare them with the observed ones --- they are consistent within a factor of two if the FR I type evolution is realized. Finally, we find that the jets with a power >~ 10^{44} erg s^{-1} are less affected by the ambient medium, and some of them, but not all, could serve as precursors of the FR II sources.

(Abridged) We propose a new method to estimate the Bondi (hot gas) accretion rates, \dot M_B, onto the supermassive black holes (SMBHs) at the centres of elliptical galaxies. It can be applied even if the Bondi radius is not well-resolved in X-ray observations. This method is based on two simple assumptions: (1) hot gas outside the Bondi radius is in nearly a hydrostatic equilibrium, and (2) the gas temperature near the galaxy centre is close to the virial temperature. We apply this method to 28 bright elliptical galaxies in nearby galaxy clusters (27 of them are the brightest cluster galaxies; BCGs). We find a correlation between the Bondi accretion rates and the power of jets associated with the SMBHs over four orders of magnitude in \dot M_B. For most galaxies, the accretion rates are large enough to account for the jet powers. Our results indicate that the feedback from the active galactic nuclei (AGN) correlates with the properties of the hot gas surrounding the SMBHs. We also find that more massive SMBHs in BCGs tend to have larger specific growth rates. This may explain the hyper masses (~ 10^10 M_sun) of some of the SMBHs. Comparison between the accretion rates and the X-ray luminosities of the AGN suggests that the AGN in the BCGs are extremely radiatively inefficient compared with X-ray binaries in the Milky Way, even when their Eddington accretion ratio, \dot M_B/\dot M_Edd, exceeds 0.01. Lastly, we find a tight correlation between the Bondi accretion rates and the X-ray luminosities of cool cores. Their relation is linear and the power generated by the Bondi accretion is large enough to compensate the radiative cooling of the cool cores. Although the `classical' Bondi accretion model is a greatly oversimplified one, the correlations we find here demonstrate that the accretion onto the SMBHs reflects broadly the properties of the Bondi accretion.

We use cosmological adaptive mesh refinement (AMR) code Enzo zoom-in simulations to study the long term evolution of the collapsing gas within dark matter (DM) halos at high redshifts. This direct collapse process is a leading candidate for rapid formation of supermassive black hole (SMBH) seeds at high z. To circumvent the Courant condition at small radii, we have used the sink particle method, and focus on the evolution on scales ~0.01-10 pc. The collapse proceeds in two stages, with the secondary runaway happening within the central 10 pc, and with no detected fragmentation. The sink particles form when the collapsing gas requires additional refinement of the grid size at the highest refinement level. Their mass never exceeds ~10^3 Mo, with the sole exception of the central seed which grows dramatically to ~ 2 x 10^6 Mo in ~2 Myr, confirming the feasibility of this path to the SMBH. The time variability of angular momentum axis in the accreted gas results in the formation of two misaligned disks --- a small inner disk, and a more massive, outer disk which is inclined by ~45^o to the inner disk. The self-gravity of these disks is heavily diluted --- both disks lie within the Roche limit of the central seed. While the inner disk is geometrically thin and weakly asymmetric, the outer disk flares due to turbulent motions as a result of the massive inflow along a pair of penetrating filaments. The geometry of inflow via filaments determines the dominant and secondary Fourier modes in this disk --- these modes have a non-self-gravitational origin. We do not confirm that m=1 is a principal mode that drives the inflow in the presence of a central massive object. While the positions of the disks depend on the scale chosen to break the self-similar collapse, the overall configuration appears to be generic, and is expected to form when the central seed becomes sufficiently massive.

We analyse the environment of the supermassive black hole (SMBH) in the centre of a massive elliptical galaxy NGC 1275 in the Perseus cluster, hosting the radio source 3C 84. We focus on the young radio lobe observed inside the estimated Bondi accretion radius. We discuss the momentum balance between the jet associated with the lobe and the surrounding gas. The results are compared with the proper motion of the radio lobe obtained with the very long baseline interferometry. We find that under assumption of a high-density environment >~ 100 cm^-3), the jet power must be comparable to the Eddington luminosity --- this is clearly inconsistent with the current moderate activity of 3C 84, which indicates instead that the jet is expanding in a very low density region (<~1 cm^-3), along the rotation axis of the accretion flow. The power required for the jet to expand in the low-density environment is comparable to the past average jet power estimated from the X-ray observations. We estimate the classical Bondi accretion rate, assuming that (1) gas accretion is spherically symmetric, (2) accretion is associated with the jet environment, and (3) the medium surrounding the jet is representative of the properties of the dominant accreting gas. We find that Bondi accretion is inconsistent with the estimated jet power. This means that either accretion of the cold gas in the NGC 1275 is more efficient than that of the hot gas, or the jets are powered by the SMBH spin.

We study the early stage of the formation of seed supermassive black holes via direct collapse in dark matter (DM) halos, in the cosmological context. We perform high-resolution zoom-in simulations of such collapse at high-$z$. Using the adaptive mesh refinement code ENZO, we resolve the formation and growth of a DM halo, until its virial temperature reaches $\sim 10^4$K, atomic cooling turns on, and collapse ensues. We demonstrate that direct collapse proceeds in two stages, although they are not well separated. The first stage is triggered by the onset of atomic cooling, and leads to rapidly increasing accretion rate with radius, from $\dot M\sim 0.1\,M_\odot {\rm yr^{-1}}$ at the halo virial radius to few $M_\odot \,{\rm yr^{-1}}$, around the scale radius $R_{\rm s}\sim 30$pc of the NFW DM density profile. The second stage of the collapse commences when the gas density takes precedence over the DM density. This is associated with the gas decoupling from the DM gravitational potential. The ensuing collapse approximates that of an isothermal sphere with $\dot M ( r )\sim $const. We confirm that the gas loses its angular momentum through non-axisymmetric perturbations and gravitational torques, to overcome the centrifugal barrier. During the course of the collapse, this angular momentum transfer process happens on nearly all spatial scales, and the angular momentum vector of the gas varies with position and time. Collapsing gas also exhibits supersonic turbulent motions which suppress gas fragmentation, and are characterized by density PDF consisting of a lognormal part and a high-density power law tail.

We use high-resolution zoom-in cosmological simulations of galaxies of Romano-Diaz et al., post-processing them with a panchromatic three-dimensional radiation transfer code to obtain the galaxy UV luminosity function (LF) at z ~ 6-12. The galaxies are followed in a rare, heavily overdense region within a ~ 5-sigma density peak, which can host high-z quasars, and in an average density region, down to the stellar mass of M_star ~ 4* 10^7 Msun. We find that the overdense regions evolve at a substantially accelerated pace --- the most massive galaxy has grown to M_star ~ 8.4*10^10 Msun by z = 6.3, contains dust of M_dust~ 4.1*10^8 Msun, and is associated with a very high star formation rate, SFR ~ 745 Msun/yr.The attained SFR-M_star correlation results in the specific SFR slowly increasing with M_star. Most of the UV radiation in massive galaxies is absorbed by the dust, its escape fraction f_esc is low, increasing slowly with time. Galaxies in the average region have less dust, and agree with the observed UV LF. The LF of the overdense region is substantially higher, and contains much brighter galaxies. The massive galaxies are bright in the infrared (IR) due to the dust thermal emission, with L_IR~ 3.7*10^12 Lsun at z = 6.3, while L_IR < 10^11 Lsun for the low-mass galaxies. Therefore, ALMA can probe massive galaxies in the overdense region up to z ~ 10 with a reasonable integration time. The UV spectral properties of disky galaxies depend significantly upon the viewing angle.The stellar and dust masses of the most massive galaxy in the overdense region are comparable to those of the sub-millimetre galaxy (SMG) found by Riechers et al. at z = 6.3, while the modelled SFR and the sub-millimetre flux fall slightly below the observed one. Statistical significance of these similarities and differences will only become clear with the upcoming ALMA observations.

We have explored prevailing modes of galaxy growth for redshifts z ~ 6-14, comparing substantially overdense and normal regions of the universe, using high-resolution zoom-in cosmological simulations. Such rare overdense regions have been projected to host high-z quasars. We demonstrate that galaxies in such environments grow predominantly by a smooth accretion from cosmological filaments which dominates the mass input from major, intermediate and minor mergers. We find that by z ~6, the accumulated galaxy mass fraction from mergers falls short by a factor of 10 of the cumulative accretion mass for galaxies in the overdense regions, and by a factor of 5 in the normal environments. Moreover, the rate of the stellar mass input from mergers also lies below that of an in-situ star formation (SF) rate. The fraction of stellar masses in galaxies contributed by mergers in overdense regions is ~12%, and ~33% in the normal regions, at these redshifts. Our median SF rates for ~few X 10^9 Mo galaxies agrees well with the recently estimated rates for z ~ 7 galaxies from Spitzer's SURF-UP survey. Finally, we find that the main difference between the normal and overdense regions lies in the amplified growth of massive galaxies in massive dark matter halos. This leads to the formation of >= 10^{10} Mo galaxies due to the ~100-fold increase in mass during the above time period. Such galaxies are basically absent in the normal regions at these redshifts.

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 use numerical simulations to explore whether direct collapse can lead to the formation of SMBH seeds at high-z. We follow the evolution of gas within DM halos of 2 x 10^8 Mo and 1 kpc. We adopt cosmological density profiles and j-distributions. Our goal is to understand how the collapsing flow overcomes the centrifugal barrier and whether it is subject to fragmentation. We find that the collapse leads either to a central runaway or to off-center fragmentation. A disk-like configuration is formed inside the centrifugal barrier. For more cuspy DM distribution, the gas collapses more and experiences a bar-like perturbation and a central runaway. We have followed this inflow down to ~10^{-4} pc. The flow remains isothermal and the specific angular momentum is efficiently transferred by gravitational torques in a cascade of nested bars. This cascade supports a self-similar, disk-like collapse. In the collapsing phase, virial supersonic turbulence develops and fragmentation is damped. For larger initial DM cores the timescales become longer. In models with more organized initial rotation, a torus forms and appears to be supported by turbulent motions. The evolution depends on the competition between two timescales, corresponding to the onset of the central runaway and off-center fragmentation. For less organized rotation, the torus is greatly weakened, the central accretion timescale is shortened, and off-center fragmentation is suppressed --- triggering the central runaway even in previously `stable' models. The resulting SMBH masses lie in the range 2 x 10^4 Mo - 2 x 10^6 Mo, much higher than for Population III remnants. We argue that the above upper limit appears to be more realistic mass. Corollaries of this model include a possible correlation between SMBH and DM halo masses, and similarity between the SMBH and halo mass functions, at time of formation.

I review the subject of the cosmological evolution of galaxies, including different aspects of growth in disk galaxies, by focussing on the angular momentum problem, mergers, and their by-products. I discuss the alternative to merger-driven growth -- cold accretion and related issues. In the follow-up, I review possible feedback mechanisms and their role in galaxy evolution. Special attention is paid to high-redshift galaxies and their properties. In the next step, I discuss a number of processes, gas- and stellar-dynamical, within the central kiloparsec of disk galaxies, and their effect on the larger spatial scales, as well as on the formation and fuelling of the seed black holes in galactic centres at high redshifts.

Theoretical and numerical modeling of dark-matter halo assembly predicts that the most luminous galaxies at high redshift are surrounded by overdensities of fainter companions. We test this prediction with HST observations acquired by our Brightest of Reionizing Galaxies (BoRG) survey, which identified four very bright z~8 candidates as Y-dropout sources in four of the 23 non-contiguous WFC3 fields observed. We extend here the search for Y-dropouts to fainter luminosities (M_* galaxies with M_AB\sim-20), with detections at >5sigma confidence (compared to >8sigma confidence adopted earlier) identifying 17 new candidates. We demonstrate that there is a correlation between number counts of faint and bright Y-dropouts at >99.84% confidence. Field BoRG58, which contains the best bright z\sim8 candidate (M_AB=-21.3), has the most significant overdensity of faint Y-dropouts. Four new sources are located within 70arcsec (corresponding to 3.1 comoving Mpc at z=8) from the previously known brighter z\sim8 candidate. The overdensity of Y-dropouts in this field has a physical origin to high confidence (p>99.975%), independent of completeness and contamination rate of the Y-dropout selection. We modeled the overdensity by means of cosmological simulations and estimate that the principal dark matter halo has mass M_h\sim(4-7)x10^11Msun (\sim5sigma density peak) and is surrounded by several M_h\sim10^11Msun halos which could host the fainter dropouts. In this scenario, we predict that all halos will eventually merge into a M_h>2x10^14Msun galaxy cluster by z=0. Follow-up observations with ground and space based telescopes are required to secure the z\sim8 nature of the overdensity, discover new members, and measure their precise redshift.

Using a high-resolution cosmological numerical simulation, we have analyzed the evolution of galaxies at z~10 in a highly overdense region of the universe. These objects could represent the high redshift galaxies recently observed by the Hubble's WFC3, and be as well possible precursors of QSOs at z~6-7. To overcome the sampling and resolution problems in cosmological simulations, we have used the Constrained Realizations method. Our main result for z~10 shows the region of 3.5h^{-1}Mpc radius in comoving coordinates completely dominated by disk galaxies in the total mass range of >=10^9h^{-1}Mo. We have verified that the gaseous and stellar disks we identify are robust morphological features, capable of surviving the ongoing merger process at these redshifts. Below this mass range, we find a sharp decline in the disk fraction to negligible numbers. At this redshift, the disks appear to be gas-rich and the dark matter halos baryon-rich, by a factor of ~2-3 above the average fraction of baryons in the universe. The prevalence of disk galaxies in the high density peaks during the epoch of reionization is contrary to the morphology-density trend observed at low redshifts.

Observed high-redshift QSOs, at z~6, may reside in massive dark matter (DM) halos of more than 10^{12} Msun and are thus expected to be surrounded by overdense regions. In a series of 10 constrained simulations, we have tested the environment of such QSOs. Comparing the computed overdensities with respect to the unconstrained simulations of regions empty of QSOs, assuming there is no bias between the DM and baryon distributions, and invoking an observationally-constrained duty-cycle for Lyman Break Galaxies, we have obtained the galaxy count number for the QSO environment. We find that a clear discrepancy exists between the computed and observed galaxy counts in the Kim et al. (2009) samples. Our simulations predict that on average eight z~6 galaxies per QSO field should have been observed, while Kim et al. detect on average four galaxies per QSO field compared to an average of three galaxies in a control sample (GOODS fields). While we cannot rule out a small number statistics for the observed fields to high confidence, the discrepancy suggests that galaxy formation in the QSO neighborhood proceeds differently than in the field. We also find that QSO halos are the most massive of the simulated volume at z~6 but this is no longer true at z~3. This implies that QSO halos, even in the case they are the most massive ones at high redshifts, do not evolve into most massive galaxy clusters at z=0.

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.

Rapid inflows associated with early galaxy formation lead to the accumulation of self-gravitating gas in the centers of proto-galaxies. Such gas accumulations are prone to non-axisymmetric instabilities, as in the well-known Maclaurin sequence of rotating ellipsoids, which are accompanied by a catastrophic loss of angular momentum (J). Self-gravitating gas is also intuitively associated with star formation. However, recent simulations of the infall process display highly turbulent continuous flows. We propose that J-transfer, which enables the inflow, also suppresses fragmentation. Inefficient J loss by the gas leads to decay of turbulence, triggering global instabilities and renewed turbulence driving. Flow regulated in this way is stable against fragmentation, whilst staying close to the instability threshold for bar formation -- thick self-gravitating disks are prone to global instabilities before they become unstable locally. On smaller scales, the fraction of gravitationally unstable matter swept up by shocks in such a flow is a small and decreasing function of the Mach number. We conclude counterintuitively that gas able to cool down to a small fraction of its virial temperature will not fragment as it collapses. This provides a venue for supermassive black holes to form via direct infall, without the intermediary stage of forming a star cluster. Some black holes could have formed or grown in massive halos at low redshifts. Thus the fragmentation is intimately related to J redistribution within the system: it is less dependent on the molecular and metal cooling but is conditioned by the ability of the flow to develop virial, supersonic turbulence.

We review some general properties of assembling galactic dark matter (DM) halos which have a direct effect on the baryon dynamics. Specifically, we focus on the mutual dynamical feedback between baryons and DM which influence disk formation and evolution, by comparing models with and without baryons evolved from identical initial conditions. Baryons are found to be capable of modifying the DM density profiles of inner halos, leading to isothermal profiles rather than the NFW cusps. Those can be leveled off by the subsequent dynamical friction of the subhalos at lower redshifts. Furthermore, subhalos appear efficient in triggering galactic bars and ablating cold gas from the disks, and therefore quenching of star formation there.

Galactic disks can form in asymmetric potentials of the assembling dark matter (DM) halos, giving rise to the first generation of gas-rich bars. Properties of these bars differ from canonical bars analyzed so far. Moreover, rapid disk growth is associated with the influx of clumpy DM and baryons along the large-scale filaments. Subsequent interactions between this substructure and the disk can trigger generations of bars, which can explain their ubiquity in the Universe. I provide a brief summary of such bar properties and argue that they fit naturally within the broad cosmological context of a hierarchical buildup of structure in the universe.

We examine the bar instability in galactic models with an exponential disk and a cuspy dark matter (DM) halo with a Navarro-Frenk-White (NFW) cosmological density profile. We construct equilibrium models from a 3-integral composite distribution function that are subject to the bar instability. We generate a sequence of models with a range of mass resolution from 1.8K to 18M particles in the disk and 10K to 100M particles in the halo along with a multi-mass model with an effective resolution of ~10^10 particles. We describe how mass resolution affects the bar instability, including its linear growth phase, the buckling instability, pattern speed decay through the resonant transfer of angular momentum to the DM halo, and the possible destruction of the halo cusp. Our higher resolution simulations show a converging spectrum of discrete resonance interactions between the bar and DM halo orbits. As the pattern speed decays, orbital resonances sweep through most of the DM halo phase space and widely distribute angular momentum among the halo particles. The halo does not develop a flat density core and preserves the cusp, except in the region dominated by gravitational softening. The formation of the bar increases the central stellar density and the DM is compressed adiabatically increasing the halo central density by 1.7X. Overall, the evolution of the bar displays a convergent behavior for halo particle numbers between 1M and 10M particles, when comparing bar growth, pattern speed evolution, the DM halo density profile and a nonlinear analysis of the orbital resonances. Higher resolution simulations clearly illustrate the importance of discrete resonances in transporting the angular momentum from the bar to the halo.

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 examine the bar instability in models with an exponential disk and a cuspy NFW-like dark matter (DM) halo inspired by cosmological simulations. Bar evolution is studied as a function of numerical resolution in a sequence of models spanning 10K to 100M DM particles - including a multi-mass model with an effective resolution of 10G. The goal is to find convergence in dynamical behaviour. We characterize the bar growth, the buckling instability, pattern speed decay through resonant transfer of angular momentum, and possible destruction of the DM halo cusp. Overall, most characteristics converge in behaviour in detail for halos containing more than 10M particles. Notably, the formation of the bar does not destroy the density cusp in this case. These higher resolution simulations clearly illustrate the importance of discrete resonances in transporting angular momentum from the bar to the halo.

Ring-shaped morphologies of nuclear star-forming regions within the central 40-200 pc of disk galaxies have been barely resolved so far in three composite Sy2 nuclei, the Sy2 Circinus galaxy and in three non-AGN galaxies. Such morphologies resemble those of the standard 1 kpc-size nuclear rings that lie in the inner Lindblad resonance regions of disk galaxies and, if they have a similar origin, represent recent radial gas inflows tantalisingly close to the central supermassive black holes. We aim to identify the population of such ultra-compact nuclear rings (UCNRs) and study their properties in relation to those of the host galaxies. From archival Hubble Space Telescope UV and Halpha images and from dust structure maps of the circumnuclear regions in nearby galaxies, we analyse the morphology of the star formation and dust, specifically searching for ring structures on the smallest observable scales. In a sample of 38 galaxies studied, we have detected a total of four new UCNRs, 30-130 pc in radius, in three different galaxies. Including our confirmation of a previous UCNR detection, this yields a UCNR fraction of roughly 10%, although our sample is neither complete nor unbiased. For the first time we resolve UCNRs in two LINERs. Overall the UCNR phenomenon appears widespread and limited neither to late-type galaxies nor exclusively to AGN hosts.

Dark matter (DM) halos formed in CDM cosmologies seem to be characterized by a power law phase-space density profile. The density of the DM halos is often fitted by the NFW profile but a better fit is provided by the Sersic fitting formula. These relations are empirically derived from cosmological simulations of structure formation but have not yet been explained on a first principle basis. Here we solve the Jeans equation under the assumption of a spherical DM halo in dynamical equilibrium, that obeys a power law phase space density and either the NFW-like or the Sersic density profile. We then calculate the velocity anisotropy, beta(r), analytically. Our main result is that for the NFW-like profile the beta - gamma relation is not a linear one (where gamma is the logarithmic derivative of the density rho[r]). The shape of beta(r) depends mostly on the ratio of the gravitational to kinetic energy within the NFW scale radius R_s. For the Sersic profile a linear beta - gamma relation is recovered, and in particular for the Sersic index of n = 6.0 case the linear fit of Hansen & Moore is reproduced. Our main result is that the phase-space density power law, the Sersic density form and the linear beta - gamma dependence constitute a consistent set of relations which obey the spherical Jeans equation and as such provide the framework for the dynamical modeling of DM halos.

This paper has been removed by arXiv admin because it was an erroneous duplicate of astro-ph/9411031.

Some aspects of disk-halo interactions for models of in and out of equilibrium disk galaxies are reviewed. Specifically, we focus on disk-halo resonant interaction without and in the presence of a gas component. Another issue is the disk growth within an assembling triaxial dark matter halo. We argue that while the triaxiality is the result of the merger process and the radial orbit instability, it is the developing chaos that damps the first generation of bars and washes out the halo prolateness. This chaos is triggered by the gravitational quadrupole interaction(s) in the system and supported by a number of other processes which are characteristic of baryons.

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 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.

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 review our recent work on the formation and evolution of disks with in triaxial dark matter (DM) halos by means of numerical simulations, including star formation and feedback from stellar evolution. The growing disks are strongly in fluenced by shapes of DM halos and modify them in turn. Disk parameters are in a broad agreement with those in the local universe. Gas-rich stellar bars grow in tandem with the disk and facilitate the angular momentum redistribution in the system and radial gas inflow. Nested bars appear to form as a by-product. Interactions between various non-axisymmetric components -- bars, disks and halos lead to decay of bars or washing out of ellipticity in the inner halo.

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.

Unified schemes of active galactic nuclei (AGN) require an obscuring dusty torus around the central engine. The compact sizes (only a few pc) determined in recent high-resolution observations require that the obscuring matter be clumpy and located inside the region where the black-hole gravity dominates over the galactic bulge. This location is in line with the scenario depicting the torus as the region of the clumpy wind coming off the accretion disk in which the clouds are dusty and optically thick. We study here the outflow scenario within the framework of hydromagnetic disk winds, incorporating the cloud properties determined from detailed modeling of the IR emission from clumpy tori. We find that torus clouds were likely detected in recent water maser observations of NGC 3079. In the wind scenario, the AGN main dynamic channel for release of accreted mass seems to be switching at low luminosities from torus outflow to radio jets. The torus disappears when the bolometric luminosity decreases below about \E{42} erg/sec because the accretion onto the central black hole can no longer sustain the required cloud outflow rate. This disappearance seems to have been observed in both LINERs and radio galaxies. With further luminosity decrease, suppression of cloud outflow spreads radially inward from the disk's dusty, molecular region into its atomic, ionized zone, resulting in disappearance of the broad emission line region at lower luminosities, yet to be determined.

We study the growth of galactic disks in live triaxial DM halos. The halos have been assembled through constrained realizations method and evolved from the linear regime using cosmological simulations. The `seed' disks have been inserted at redshift z=3 and increased in mass tenfold over various time periods, ~1-3 Gyr, with the halo responding quasi-adiabatically to this process. We follow the dynamical and secular evolution of the disk-halo system and analyze changes in the most important parameters, like 3-D DM shapes, stellar and DM radial density profiles, stellar bar development, etc. We find that a growing disk is responsible for washing out the halo prolateness and for diluting its flatness over a period of time comparable to the disk growth. Moreover, we find that a disk which contributes more to the overall rotation curve in the system is also more efficient in axisymmetrizing the halo, without accelerating the halo figure rotation. The observational corollary is that the maximal disks probably reside in nearly axisymmetric halos, while disks whose rotation is dominated by the halo at all radii are expected to reside in more prolate halos. The halo shape is sensitive to the final disk mass, but is independent of how the seed disk is introduced into the system. We also expect that the massive disks are subject to a bar instability, while light disks have this instability damped by the halo triaxiality. Implications to the cosmological evolution of disks embedded in asymmetric halos are discussed and so are the corollaries for the observed fraction of stellar bars. Finally, the halo responds to the stellar bar by developing a gravitational wake -- a `ghost' bar of its own which is almost in-phase with that in the disk.

Context: The morphology of massive star formation in the central regions of galaxies is an important tracer of the dynamical processes that govern the evolution of disk, bulge, and nuclear activity. Aims: We present optical imaging of the central regions of a sample of 73 spiral galaxies in the H alpha line and in optical broad bands, and derive information on the morphology of massive star formation. Methods: We obtained images with the William Herschel Telescope, mostly at a spatial resolution of below one second of arc. For most galaxies, no H alpha imaging is available in the literature. We outline the observing and data reduction procedures, list basic properties, and present the I-band and continuum-subtracted H alpha images. We classify the morphology of the nuclear and circumnuclear H alpha emission and explore trends with host galaxy parameters. Results: We confirm that late-type galaxies have a patchy circumnuclear appearance in H alpha, and that nuclear rings occur primarily in spiral types Sa-Sbc. We identify a number of previously unknown nuclear rings, and confirm that nuclear rings are predominantly hosted by barred galaxies. Conclusions: Other than in stimulating nuclear rings, bars do not influence the relative strength of the nuclear H alpha peak, nor the circumnuclear H alpha morphology. Even though our selection criteria led to an over-abundance of galaxies with close massive companions, we do not find any significant influence of the presence or absence of a close companion on the relative strength of the nuclear H alpha peak, nor on the H alpha morphology around the nucleus.

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.

We analyze the evolution of stellar bars in galactic disks in mildly triaxial flat-core and cuspy CDM halos. We use tailored simulations of rigid and live halos which include the feedback from disk/bar onto the halo in order to test the work by El-Zant & Shlosman (2002). The latter used the Liapunov exponents to analyze the fate of bars in analytical asymmetric halos. We find: (1) The bar growth is similar in all rigid axisymmetric and triaxial halos. (2) Bars in live models vertically buckle and form a pseudobulge with a boxy/peanut shape. (3) In live axisymmetric halos, the bar strength varies little during the secular evolution. The bar pattern speed anticorrelates with the halo core size. The bar strength is larger for smaller disk-to-halo mass ratios within disk radii, the bar size correlates with the halo core sizes, and the bar pattern speeds -- with the halo central mass concentration. Bars embedded in live triaxial halos have a starkly different fate: they dissolve on ~1.5-5 Gyr due to the onset of chaos over continuous zones, leaving behind a weak oval distortion. The onset of chaos is related to the halo triaxiality, the fast rotating bar and the halo cuspiness. Before the bar dissolves, the region outside it develops strong spiral structures, especially in the live triaxial halos. (4) More angular momentum is absorbed by the triaxial halos as compared to the axisymmetric models and its exchange is mediated by resonances. (5) Cuspy halos are more susceptible than flat-core halos to having their prolateness washed out by the bar. We analyze these results in terms of the stability of trajectories and development of chaos. We set constraints on the triaxiality of DM halos by comparing our predictions to recent observations of bars out to z~1.

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.

The central regions of disk galaxies are hosts to supermassive black holes whose masses show a tight correlation with the properties of surrounding stellar bulges. While the exact origin of this dependency is not clear, it can be related to the very basic properties of dark matter halos and the associated gas and stellar dynamics in the central kpc of host galaxies. In this review we discuss some of the recent developments in modeling the wide spectrum of dynamical processes which can be affiliated with the above phenomena, such as the structure of molecular tori in AGN, structure formation in triaxial halos, and dissipative and non-dissipative dynamics in nested bar systems, with a particular emphasis on decoupling of gaseous nuclear bars. We also briefly touch on the subject of fueling the nuclear starbursts and AGN.

A substantial fraction of barred galaxies host additional nuclear bars which tumble with pattern speeds exceeding those of the large-scale (primary) stellar bars. We have investigated the mechanism of formation and dynamical decoupling in such nested bars which include gaseous (secondary) nuclear bars within the full size galactic disks, hosting a double inner Lindblad resonance. Becoming increasingly massive and self-gravitating, the nuclear bars lose internal (circulation) angular momentum to the primary bars and increase their strength. Developing chaos within these bars triggers a rapid gas collapse -- bar contraction. During this time period, the secondary bar pattern speed Omega_s~a^{-1}, where "a" stands for the bar size. As a result, Omega_s increases dramatically until a new equilibrium is reached (if at all), while the gas specific angular momentum decreases -- demonstrating the dynamical decoupling of nested bars. Viscosity, and therefore the gas presence, appears to be a necessary condition for the prograde decoupling of nested bars. This process maintains an inflow rate of ~1 M_o/yr over ~10^8 yrs across the central 200 pc and has important implications for fueling the nuclear starbursts and AGN.

One third of present-day spirals host optically visible strong bars that drive their dynamical evolution. However, the fundamental question of how bars evolve over cosmological times has yet to be addressed, and even the frequency of bars at intermediate redshifts remains controversial. We investigate the frequency of bars out to z~1.0 drawing on a sample of 1590 galaxies from the GEMS survey, which provides morphologies from HST ACS two-color images, and highly accurate redshifts from the COMBO-17 survey. We identify spiral galaxies using the Sersic index, concentration parameter, and rest-frame color. We characterize bars and disks by fitting ellipses to F606W and F850LP images, taking advantage of the two bands to minimize bandpass shifting. We exclude highly inclined (i>60 deg) galaxies to ensure reliable morphological classifications, and apply completeness cuts of M_v <= -19.3 and -20.6. More than 40% of the bars that we detect have semi major axes a<0.5" and would be easily missed in earlier surveys without the small PSF of ACS. The bars that we can reliably detect are fairly strong (with ellipticities e>=0.4) and have a in the range ~1.2-13 kpc. We find that the optical fraction of such strong bars remains at ~(30% +- 6%) from the present-day out to look-back times of 2-6 Gyr (z~0.2-0.7) and 6-8 Gyr (z~0.7-1.0); it certainly shows no sign of a drastic decline at z>0.7. Our findings of a large and similar bar fraction at these three epochs favor scenarios in which cold gravitationally unstable disks are already in place by z~1, and where on average bars have a long lifetime (well above 2 Gyr). The distributions of structural bar properties in the two slices are, however, not statistically identical and therefore allow for the possibility that the bar strengths and sizes may evolve over time.

Bars drive the dynamical evolution of disk galaxies by redistributing mass and angular momentum, and they are ubiquitous in present-day spirals. Early studies of the Hubble Deep Field reported a dramatic decline in the rest-frame optical bar fraction f_opt to below 5% at redshifts z>0.7, implying that disks at these epochs are fundamentally different from present-day spirals. The GEMS bar project, based on ~8300 galaxies with HST-based morphologies and accurate redshifts over the range 0.2-1.1, aims at constraining the evolution and impact of bars over the last 8 Gyr. We present early results indicating that f_opt remains nearly constant at ~30% over the range z=0.2-1.1,corresponding to lookback times of ~2.5-8 Gyr. The bars detected at z>0.6 are primarily strong with ellipticities of 0.4-0.8. Remarkably, the bar fraction and range of bar sizes observed at z>0.6 appear to be comparable to the values measured in the local Universe for bars of corresponding strengths. Implications for bar evolution models are discussed.

Young stellar bars in disk galaxies experience a vertical buckling instability which terminates their growth and thickens them, resulting in a characteristic peanut/boxy shape when viewed edge on. Using N-body simulations of galactic disks embedded in live halos, we have analyzed the bar structure throughout this instability and found that the outer third of the bar dissolves completely while the inner part (within the vertical inner Lindblad resonance) becomes less oval. The bar acquires the frequently observed peanut/boxy-shaped isophotes. We also find that the bar buckling is responsible for a mass injection above the plane, which is subsequently trapped by specific 3-D families of periodic orbits of particular shapes explaining the observed isophotes, in line with previous work. Using a 3-D orbit analysis and surfaces of sections, we infer that the outer part of the bar is dissolved by a rapidly widening stochastic region around its corotation radius -- a process related to the bar growth. This leads to a dramatic decrease in the bar size, decrease in the overall bar strength and a mild increase in its pattern speed, but is not expected to lead to a complete bar dissolution. The buckling instability appears primarily responsible for shortening the secular diffusion timescale to a dynamical one when building the boxy isophotes. The sufficiently long timescale of described evolution, ~1 Gyr, can affect the observed bar fraction in local universe and at higher redshifts, both through reduced bar strength and the absence of dust offset lanes in the bar.

Sand, Treu, & Ellis (2002) have measured the central density profile of cluster MS2137-23 with gravitational lensing and velocity dispersion and removed the stellar contribution with a reasonable M/L. The resulting dark matter distribution within r<50 kpc was fitted by a density cusp of r^{-beta} with beta=0.35. This stands in an apparent contradiction to the CDM prediction of beta~1, and the disagreement worsens if adiabatic compression of the dark matter by the infalling baryons is considered. Following El-Zant, Shlosman & Hoffman (2001), we argue that dynamical friction acting on galaxies moving within the dark matter background counters the effect of adiabatic compression by transfering the orbital energy of galaxies to the dark matter, thus heating up and softening the central density cusp. Using N-body simulations of massive solid clumps moving in clusters we show that indeed the inner dark matter distribution flattens (with beta approx 0.35 for a cluster like MS2137-23) when the galaxies spiral inward. We find as a robust result that while the dark matter distribution becomes core-like, the overall mass distribution preserves its cuspy nature, in agreement with X-ray and lensing observations of clusters.

A substantial fraction of disk galaxies is double-barred. We analyze the dynamical stability of such nested bar systems by means of Liapunov exponents,by fixing a generic model and varying the inner (secondary) bar mass. We show that there exists a critical mass below which the secondary bar cannot sustain its own orbital structure, and above which it progressively destroys the outer (primary) bar-supporting orbits. In this critical state, a large fraction of the trajectories (regular and chaotic) are aligned with either bar, suggesting the plausibility of long-lived dynamical states when secondary-to-primary bar mass ratio is of the order of a few percent. Qualitatively similar results are obtained by varying the size of the secondary bar, within certain limits, while keeping its mass constant. In both cases, an important role appears to be played by chaotic trajectories which are trapped around (especially) the primary bar for long periods of time.

We discuss different aspects of nested bar dynamics and its effect on the gas flow and fueling of Active Galactic Nuclei. Specifically we focus on the dynamical decoupling between the primary and secondary bars and the gas flow across the bar-bar interface. We analyze the nuclear gaseous bar formation when gas gravity can be neglected or when it dominates. Finally, we discuss the possible effect of flat core, triaxial, dark halos on the formation of galactic bulges and supermassive black holes (SBHs) and argue in favor of SBH-bulge-halo correlation.

The masses of supermassive black holes (SBHs) show correlations with bulge properties in disk and elliptical galaxies. We study the formation of galactic structure within flat-core triaxial haloes and show that these correlations can be understood within the framework of a baryonic component modifying the orbital structure in the underlying potential. In particular, we find that terminal properties of bulges and their central SBHs are constrained by the destruction of box orbits in the harmonic cores of dark haloes and the emergence of progressively less eccentric loop orbits there. SBH masses, M_SBH, should exhibit a tighter correlation with bulge velocity dispersions, sigma_B, than with bulge masses, M_B, in accord with observations, if there is a significant scatter in the M_H-sigma_H relation for the halo. In the context of this model the observed M_SBH-sigma_B relation implies that haloes should exhibit a Faber-Jackson type relationship between their masses and velocity dispersions. The most important prediction of our model is that halo properties determine the bulge and SBH parameters. The model also has important implications for galactic morphology and the process of disk formation.

Through observations and modeling, we demonstrate how the recently discovered large-scale bar in NGC 5248 generates spiral structure which extends from 10 kpc down to 100 pc, fuels star formation on progressively smaller scales, and drives disk evolution. Deep inside the bar, two massive molecular spirals cover nearly 180 degrees in azimuth, show streaming motions of 20--40 km/s and feed a starburst ring of super star clusters at 375 pc. They also connect to two narrow K-band spirals which delineate the UV-bright star clusters in the ring. The data suggest that the K-band spirals are young, and the starburst has been triggered by a bar-driven spiral density wave (SDW). The latter may even have propagated closer to the center where a second Halpha ring and a dust spiral are found. The molecular and HST data support a scenario where stellar winds and supernovae efficiently clear out gas from dense star-forming regions on timescales less than a few Myrs. We have investigated the properties of massive CO spirals within the framework of bar-driven SDWs, incorporating the effect of gas self-gravity. We find good agreement between the model predictions and the observed morphology, kinematics, and pitch angle of the spirals. This combination of observations and modeling provides the best evidence to date for a strong dynamical coupling between the nuclear region and the surrounding disk. It also confirms that a low central mass concentration, which may be common in late-type galaxies, is particularly favorable to the propagation of a bar-driven gaseous SDW deep into the central region of the galaxy, whereas a large central mass concentration favors other processes, such as the formation and decoupling of nuclear bars.

For decades, the grand-design SAB spiral galaxy NGC 5248 has been postulated to host a short bar of semi-major axis 22" (1.6 kpc). From dynamical and morphological arguments, we argue, however, that its spiral structure is being driven by a large-scale bar whose corotation radius lies at around 115" (8.6 kpc).Our estimate is based partially on a deep R-band image which reveals that the feature previously thought to be an inclined disk is in fact an extended stellar bar. The bar is embedded within a fainter outer disk visible out to a radius of 230" (17.2 kpc). The bar has a deprojected ellipticity of 0.44 and a semi-major axis of 95" (7.1 kpc). The classical grand-design spirals of NGC 5248, prominent in B, R, and K light, lie on the leading edge of the large-scale stellar bar and are accompanied by concave dust lanes out to at least 70". The offset between the dust and young stars is consistent with our understanding of gas flows in barred galaxies, where shocks along the leading edges of a moderately strong bar compress the gas to form massive young stars. While in many strongly barred galaxies, optical spiral arms are prominent outside the bar but not within it, NGC 5248 illustrates how intense star formation along a moderately strong bar can lead to conspicuous open spiral arms within the bar itself. NGC 5248 also provides a clear example of how a large-scale stellar bar embedded within a faint outer optical disk can be misidentified as an inclined disk when imaging studies lack the sensitivity to detect the actual outer disk. We discuss the implications for the estimated bar fraction at higher redshifts.

We investigate the stability of trajectories in barred galaxies with mildly triaxial halos by means of Liapunov exponents. This method is perfectly suitable for time-dependent 3D potentials where surfaces of sections and other simple diagnostics are not applicable. We find that when halos are centrally-concentrated most trajectories starting near the plane containing the bar become chaotic. Moreover, the shape of many of the remaining regular trajectories do not match the bar density distribution, being too round. Therefore, time-independent self-consistent solutions are highly unlikely to be found. When the non-rotating non-axisymmetric perturbation in the potential reaches 10%, almost all trajectories integrated are chaotic and have large Liapunov exponents. No regular trajectories aligned with the bar have been found. Hence, if the evolution of the density figure is directly related to the characteristic timescale of orbital instability, bar dissolution would take place on a timescale of few dynamical times. The slowly rotating non-axisymmetric contribution to the potential required for the onset of widespread chaotic behavior is remarkably small. Systems consisting of centrally-concentrated axisymmetric halos and stellar bars thus appear to be structurally unstable, and small (1%) deviations from perfect axisymmetry should result in a bar dissolution on a timescale significantly smaller than the Hubble time. Since halos found in CDM simulations of large scale structure are both centrally-concentrated and triaxial it is unlikely that stellar bars embedded in such halos would form and survive unless the halos are modified during the formation of the baryonic component.

We discuss different aspects of nested bar systems, both observational and theoretical. Such systems consist of a large-scale primary bar leading to the formation of a sub-kpc size secondary bar, whose pattern speed differs substantially from that of the main bar. Specifically, we focus on the origin and gravitational decoupling of nested bars, and their characteristic gas flows on scales of about 100 pc -- 10 kpc, with and without gas gravity. We find that the gas response in nested bars differs profoundly from that in single bars, and that no offset dust lanes form in the secondary bars. We also discuss briefly the importance of nested bar systems for redshifts corresponding to galaxy formation epoch.

The central kpc of barred galaxies exhibits a wealth of morphological information on different components with clear dynamical consequences. These include nuclear rings, spirals, bars, and more. We argue that this morphology is driven by large-scale stellar bars and analyze its consequences for gas dynamics and the distribution of star-forming regions. Specifically, we focus on gas flows in nested bar systems and study their origin, as well as the gravitational decoupling of gaseous nuclear bars with and without gas self-gravity. We find that the gas response in nested bars differs profoundly from that in single bars, and that no offset dust lanes form in the nuclear bars.

We describe selected highlights of the international conference ``The Central Kiloparsec of Starbursts and AGN: The La Palma Connection'', held in May 2001 in Los Cancajos on the Spanish island of La Palma.

We analyze the observed properties of nested and single stellar bar systems in disk galaxies. The 112 galaxies in our sample comprise the largest matched Seyfert vs. non-Seyfert sample of nearby galaxies with complete near-infrared or optical imaging sensitive to lengthscales ranging from tens of pc to tens of kpc. We find that a significant fraction of the sample galaxies, 17% +/- 4%, has more than one bar, and that 28% +/- 5% of barred galaxies have nested bars. The bar fractions appear to be stable according to reasonable changes in our adopted bar criteria. For these nested bars, we detect a clear division in length between the large-scale (primary) bars and small-scale (secondary) bars, both in absolute and normalized (to the size of the galaxy) length. We argue that this bimodal distribution can be understood within the framework of the inner Lindblad resonances (ILRs), which are located where the gravitational potential of the innermost galaxy switches effectively from 3D to 2D. While primary bars are found to correlate with the host galaxy sizes, no such correlation is observed for the secondary bars. Moreover, we find that secondary bars differ morphologically from single bars. Overall, our matched Seyfert and non-Seyfert samples show a statistically significant excess of bars among the Seyfert galaxies at practically all lengthscales. We confirm our previous result that Seyfert galaxies always show a preponderance of "thick" bars compared to the bars in non-Seyfert galaxies. Finally, no correlation is observed between the presence of a bar and that of companion galaxies, even relatively bright ones.

Line-driven winds from stars and accretion disks are accelerated by scattering in numerous line transitions. The wind is believed to adopt a unique critical solution, out of the infinite variety of shallow and steep solutions. We study the inherent dynamics of the transition towards the critical wind. A new runaway wind mechanism is analyzed in terms of radiative-acoustic (Abbott) waves which are responsible for shaping the wind velocity law and fixing the mass loss. Three different flow types result, depending on the location of perturbations. First, if the shallow solution is perturbed sufficiently far downstream, a single critical point forms in the flow, which is a barrier for Abbott waves, and the solution tends to the critical one. Second, if the shallow solution is perturbed upstream from this critical point, mass overloading results, and the critical point is shifted inwards. This wind exhibits a broad, stationary region of decelerating flow and its velocity law has kinks. Third, for perturbations even further upstream, the overloaded wind becomes time-dependent, and develops shocks and dense shells.

Observations and theory suggest that line driven winds from hot stars and luminous accretion disks adopt a unique, critical solution which corresponds to maximum mass loss rate. We analyze the numerical stability of the infinite family of shallow wind solutions, which resemble solar wind breezes, and their transition to the critical wind. Shallow solutions are sub-critical with respect to radiative (or Abbott) waves. These waves can propagate upstream through shallow winds at high speeds. If the waves are not accounted for in the Courant time step, numerical runaway results. The outer boundary condition is equally important for wind stability. Assuming pure outflow conditions, as is done in the literature, triggers runaway of shallow winds to the critical solution or to accretion flow.

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.

Coupling the dark matter (DM) and baryonic components in galaxies via dynamical friction (by allowing the latter to be clumpy) leads to energy input into the DM. The resulting expanding tendency in the DM component overcomes the competing effect of adiabatic contraction, resulting from the shrinking gas component, and leads to a nearly constant density core. The baryonic mass inflow also leads to the coupling of the structure at the kpc scale and below and that of the structure of the DM at the 10 kpc scale and above, with rate of mass inflow being dependent on the initial halo central concentration. The latter also determines the final halo core radius.

N-body simulations and analytical calculations of the gravitational collapse in an expanding universe predict that halos should form with a diverging inner density profile, the cusp. There are some observational indications that the dark matter distribution in galaxies might be characterized by a finite core. This `core catastrophe' has prompted a search for alternatives to the CDM cosmogony. It is shown here that the discrepancy between theory and observations can be very naturally resolved within the standard CDM model, provided that gas is not initially smoothly distributed in the dark matter halo, but rather is concentrated in clumps of mass $\geq 0.01 %$ the total mass of the system. Dynamical friction acting on these lumps moving in the background of the dark matter particles, dissipates the clumps orbital energy and deposits it in the dark matter. Using Monte-Carlo simulations, it is shown that the dynamical friction provides a strong enough drag, and that with realistic baryonic mass fractions, the available orbital energy of the clumps is sufficient to heat the halo and turn the primordial cusp into a finite, non-diverging core --- overcoming the competing effect of adiabatic contraction due to gravitational influence of the shrinking baryonic component. Depending on the initial conditions, the total density distribution may either become more or less centrally concentrated. Possible consequences of the proposed mechanism for other problems in the CDM model and for the formation and early evolution of the baryonic component of galaxies are also briefly discussed.

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.

Using all available major samples of Seyfert galaxies and their corresponding control samples of closely matched non-active galaxies, we find that the bar ellipticities (or axial ratios) in Seyfert galaxies are systematically different from those in non-active galaxies. Overall, there is a deficiency of bars with large ellipticities (i.e., `fat' or `weak' bars) in Seyferts, compared to non-active galaxies. Accompanied with a large dispersion due to small number statistics, this effect is strictly speaking at the `2sigma' level. To obtain this result, the active galaxy samples of near-infrared surface photometry were matched to those of normal galaxies in type, host galaxy ellipticity, absolute magnitude, and, to some extent, in redshift. We discuss possible theoretical explanations of this phenomenon within the framework of galactic evolution, and, in particular, of radial gas redistribution in barred galaxies. Our conclusions provide further evidence that Seyfert hosts differ systematically from their non-active counterparts on scales of a few kpc.

Line-driven winds from hot stars and accretion disks are thought to adopt a unique, critical solution which corresponds to maximum mass loss rate and a particular velocity law. We show that in the presence of negative velocity gradients, radiative-acoustic (Abbott) waves can drive shallow wind solutions towards larger velocities and mass loss rates. Perturbations introduced downstream from the wind critical point lead to convergence towards the critical solution. By contrast, low-lying perturbations cause evolution towards a mass-overloaded solution, developing a broad deceleration region in the wind. Such a wind differs fundamentally from the critical solution. For sufficiently deep-seated perturbations, overloaded solutions become time-dependent and develop shocks and shells.

We present a deep ASCA observation of a Broad Absorption Line Quasar (BALQSO) PG0946+301. The source was clearly detected in one of the gas imaging spectrometers, but not in any other detector. If BALQSOs have intrinsic X-ray spectra similar to normal radio-quiet quasars, our observations imply that there is Thomson thick X-ray absorption (N_H >~10^{24} cm^{-2}) toward PG0946+301. This is the largest column density estimated so far toward a BALQSO. The absorber must be at least partially ionized and may be responsible for attenuation in the optical and UV. If the Thomson optical depth toward BALQSOs is close to one, as inferred here, then spectroscopy in hard X-rays with large telescopes like XMM would be feasible.

A hydromagnetic (MHD) wind from a clumpy molecular accretion disk is invoked to explain observations of warm absorbing (WA) gas in UVX from Sy galaxies. This paper focuses on two issues: (1) compatibility of kinematics and dynamics of MHD wind with the observed properties of WAs; and (2) relationship between the UVX absorptions. We provide an in-depth comparison between the MHD model and the Sy 1 galaxy NGC 5548, which at high spectral resolution exhibits a number of discrete UV absorption components. We find that: (1) the total column densities of Ovii, Oviii and H, are reproduced by constraining the UV ion column densities of Civ and Nv in each component to lie within a factor of 2 of their observed values and optimizing over the possible sets of component ionization states and Civ column densities; (2) the WA exists in the outer part of the wind and is not a continuation of the flow in the BLR; and (3) the WA extends in radial and polar directions and is ionization-stratified. X-ray absorption is found to be heavily biased towards smaller r, and UV absorption originates at larger distances from the central continuum source. We show that the discrete absorption components along the line-of-sight are intrinsically clumpy. Density differences between kinematic components result in a range of ionization and recombination timescales. We further test the applicability of the MHD wind to WAs in general, by constructing a quasi-continuous flow model, and extending it to arbitrary aspect angles. We estimate the fraction of Sy 1s having detectable WAs with larger Ovii column density than Oviii, and the range of total H column densities. We also find that the ratio of Ovii to Oviii optical depths can serve as a new diagnostic of AGN aspect angle.

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.

We analyze formation of grand-design two-arm spiral structure in the nuclear regions of disk galaxies. Such morphology has been recently detected in a number of objects using high-resolution near-infrared observations. Motivated by the observed (1) continuity between the nuclear and kpc-scale spiral structures, and by (2) low arm-interarm contrast, we apply the density wave theory to explain the basic properties of the spiral nuclear morphology. In particular, we address the mechanism for the formation, maintenance and the detailed shape of nuclear spirals. We find, that the latter depends mostly on the shape of the underlying gravitational potential and the sound speed in the gas. Detection of nuclear spiral arms provides diagnostics of mass distribution within the central kpc of disk galaxies. Our results are supported by 2D numerical simulations of gas response to the background gravitational potential of a barred stellar disk. We investigate the parameter space allowed for the formation of nuclear spirals using a new method for constructing a gravitational potential in a barred galaxy, where positions of resonances are prescribed.

We present a detailed study of the bar fraction in the CfA sample of Seyfert galaxies, and in a carefully selected control sample of non-active galaxies, to investigate the relation between the presence of bars and of nuclear activity. To avoid the problems related to bar classification in the RC3, e.g., subjectivity, low resolution and contamination by dust, we have developed an objective bar classification method, which we conservatively apply to our new sub-arcsecond resolution near-infrared imaging data set (Peletier et al. 1999). We are able to use stringent criteria based on radial profiles of ellipticity and major axis position angle to determine the presence of a bar and its axial ratio. Concentrating on non-interacting galaxies in our sample for which morphological information can be obtained, we find that Seyfert hosts are barred more often (79% +/- 7.5%) than the non-active galaxies in our control sample (59% +/- 9%), a result which is at the 2.5 sigma significance level. The fraction of non-axisymmetric hosts becomes even larger when interacting galaxies are taken into account. We discuss the implications of this result for the fueling of central activity by large-scale bars. This paper improves on previous work by means of imaging at higher spatial resolution and by the use of a set of stringent criteria for bar presence, and confirms that the use of NIR is superior to optical imaging for detection of bars in disk galaxies.

We present new high-resolution near-infrared observations in the J,H and K bands, obtained to study the properties of Seyfert host galaxies. The dataset consists of images in the three bands of practically the entire CfA sample of Seyfert galaxies, and K-band images of a control sample of non-active, `normal', galaxies, matched to the Seyfert sample in the distribution of type and inclination. The spatial resolution and sampling of the new images is a factor 2 better than previously published K-band data. In this paper, we present the data in the form of profiles of surface brightness and color, ellipticity and major axis position angle, as well as greyscale maps of surface brightness in H or K and both J-H and H-K colors. We compare our surface brightness and color profiles with the literature, and find good agreement. Our data are discussed in detail in three subsequent publications, where we analyze the morphologies of Seyfert and normal hosts, quantify the strength of nonaxisymmetric features in disks and their relationship to nuclear activity, address the question of bar fraction in Seyferts and normal galaxies, and analyze the color information in the framework of emission mechanisms in Seyfert 1s and 2s, and in non-active galaxies.

We analyze the dynamics of 2D stationary line-driven winds from accretion disks in cataclysmic variables (CVs), by generalizing the Castor, Abbott and Klein theory. In paper 1, we have solved the wind Euler equation, derived its two eigenvalues, and addressed the solution topology and wind geometry. Here, we focus on mass loss and velocity laws. We find that disk winds, even in luminous novalike variables, have low optical depth, even in the strongest driving lines. This suggests that thick-to-thin transitions in these lines occur. For disks with a realistic radial temperature, the mass loss is dominated by gas emanating from the inner decade in r. The total mass loss rate associated with a luminosity 10 Lsun is 10^{-12} Msun/yr, or 10^{-4} of the mass accretion rate. This is one order of magnitude below the lower limit obtained from P Cygni lines, when the ionizing flux shortwards of the Lyman edge is supressed. The difficulties with such small mass loss rates in CVs are principal, and confirm our previous work. We conjecture that this issue may be resolved by detailed nonLTE calculations of the line force within the context of CV disk winds, and/or better accounting for the disk energy distribution and wind ionization structure. We find that the wind velocity profile is well approximated by the empirical law used in kinematical modeling. The acceleration length scale is given by the footpoint radius of the wind streamline in the disk. This suggests an upper limit of 10 Rwd to the acceleration scale, which is smaller by factors of a few as compared to values derived from line fitting.

We analyze the dynamics of 2-D stationary, line-driven winds from accretion disks in cataclysmic variable stars. The driving force is that of line radiation pressure, in the formalism developed by Castor, Abbott & Klein for O stars. Our main assumption is that wind helical streamlines lie on straight cones. We find that the Euler equation for the disk wind has two eigenvalues, the mass loss rate and the flow tilt angle with the disk. Both are calculated self-consistently. The wind is characterized by two distinct regions, an outer wind launched beyond four white dwarf radii from the rotation axis, and an inner wind launched within this radius. The inner wind is very steep, up to 80 degrees with the disk plane, while the outer wind has a typical tilt of 60 degrees. In both cases the ray dispersion is small. We, therefore, confirm the bi-conical geometry of disk winds as suggested by observations and kinematical modeling. The wind collimation angle appears to be robust and depends only on the disk temperature stratification. The flow critical points lie high above the disk for the inner wind, but close to the disk photosphere for the outer wind. Comparison with existing kinematical and dynamical models is provided. Mass loss rates from the disk as well as wind velocity laws are discussed in a subsequent paper.

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 analyze the results of long-term observations of broad-line region (BLR) in the Sy1 galaxy NGC 5548 and provide a critical comparison with the predictions of a hydromagnetically-driven outflow model of Emmering, Blandford and Shlosman. This model was used to generate a time series of CIV line profiles that have responded to a time varying continuum. The model includes cloud emission anisotropy, cloud obscuration, a CLOUDY-generated emissivity function and a narrow-line component used to generate the line profiles, and is driven with continuum input based on the monitoring campaigns of NGC 5548. The line strengths, profiles and lags are compared with the observations. It reproduces the basic features of CIV line variability in this AGN without varying the model parameters. The best fit model provides the effective size, the dominant geometry, the emissivity distribution and the 3D velocity field of the CIV BLR and constrains the mass of the central BH. The inner part of the wind appears to be responsible for the anisotropically emitted CIV line, while its outer part remains dusty and molecular, thus having similar spectral characteristics to a molecular torus, although its dynamics is fundamentally different. In addition, our model predicts a differential response across the CIV line profile, producing a red-side-first response in the mid-wings followed by the blue mid-wing and by the line core. Based on the comparison of data and model cross-correlation functions and 1D and 2D transfer functions, we find that the rotating outflow model is compatible with observations of the BLR in NGC 5548.

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.

We investigate gas dynamics in the presence of a double inner Lindblad resonance within a barred disk galaxy. Using an example of a prominent spiral, M100, we reproduce the basic central morphology, including four dominant regions of star formation corresponding to the compression maxima in the gas. These active star forming sites delineate an inner boundary (so-called nuclear ring) of a rather broad oval detected in the near infrared. We find that inclusion of self-gravitational effects in the gas is necessary in order to understand its behavior in the vicinity of the resonances and its subsequent evolution. The self-gravity of the gas is also crucial to estimate the effect of a massive nuclear ring on periodic orbits in the stellar bar.

We study the eclipsing nova-like variable V347 Pup by matching its UV emission line profiles in and out of eclipse to synthetic lines using a 3D kinematic and radiation transfer model. Our results support the accretion disk origin of winds in non-magnetic CVs as opposite to the WD origin. Our main point concerns the importance of rotation for the UV emission line shapes in such systems. In particular, we show that the narrowing of the UV emission lines in V347 Pup during eclipse can be easily explained by the eclipse of the innermost part of the wind by the secondary and the resulting reduction in the contribution of rotational broadening to the width of the lines. During the eclipse, the residual line flux is very sensitive to the maximal temperature of disk radiation. Good fits for reasonable mass-loss rates have been obtained for maximum disk temperatures of 50,000 degrees. This constraint was imposed either by leveling off the inner disk temperature profiles, in agreement with recent observations of some nova-like objects, or by assuming that the accretion disk does not extend to the surface of the white dwarf, in which case V347 up would be an intermediate polar. In anticipation of high-speed spectrophotometry of CVs by the HST, we provide numerical model of a time-resolved eclipse of V347 Pup or similar such system to be verified by future observations.

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').

We use a free-energy minimization approach to describe the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Our approach is fully nonlinear and stems from the Landau-Ginzburg theory of phase transitions. Here we examine higher than 2nd-harmonic disturbances applied to Maclaurin spheroids, the corresponding bifurcating sequences, and their relation to nonlinear fission processes. The triangle and ammonite sequences bifurcate from the two 3rd-harmonic neutral points on the Maclaurin sequence while the square and one-ring sequences bifurcate from two of the three known 4th harmonic neutral points. In the other three cases, secular instability does not set in at the corresponding bifurcation points because the sequences stand and terminate at higher energies relative to the Maclaurin sequence. There is no known bifurcating sequence at the point of 3rd-harmonic dynamical instability. Our nonlinear approach easily identifies resonances between the Maclaurin sequence and various multi-fluid-body sequences that cannot be detected by linear stability analyses. Resonances appear as first-order phase transitions at points where the energies of the two sequences are nearly equal but the lower energy state belongs to one of the multi-fluid-body sequences.

In Papers I and II, we have used a free-energy minimization approach that stems from the Landau-Ginzburg theory of phase transitions to describe in simple and clear physical terms the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Here we investigate the secular and dynamical 3rd-harmonic instabilities that appear first on the Jacobi sequence of incompressible zero- vorticity ellipsoids. Poincare found a bifurcation point on the Jacobi sequence where a 3rd-harmonic mode becomes neutral. A sequence of pear-shaped equilibria branches off at this point but stands at higher energies. Therefore, the Jacobi ellipsoids remain secularly and dynamically stable. Cartan found that dynamical 3rd-harmonic instability also sets in at the Jacobi-pear bifurcation. We find that Cartan's instability leads to differentially rotating objects and not to uniformly rotating pear-shaped equilibria. We demonstrate that the pear-shaped sequence exists at higher energies and at higher rotation relative to the Jacobi sequence. The Jacobi ellipsoid undergoes fission on a secular time scale and a short-period binary is produced. The classical fission hypothesis of binary-star formation of Poincare and Darwin is thus feasible.

We use a free-energy minimization approach to describe the secular and dynamical instabilities as well as the bifurcations along equilibrium sequences of rotating, self-gravitating fluid systems. Our approach is fully nonlinear and stems from the Ginzburg-Landau theory of phase transitions. In this paper, we examine fourth-harmonic axisymmetric disturbances in Maclaurin spheroids and fourth-harmonic nonaxisymmetric disturbances in Jacobi ellipsoids. These two cases are very similar in the framework of phase transitions. Irrespective of whether a nonlinear first-order phase transition occurs between the critical point and the higher turning point or an apparent second-order phase transition occurs beyond the higher turning point, the result is fission (i.e. ``spontaneous breaking'' of the topology) of the original object on a secular time scale: the Maclaurin spheroid becomes a uniformly rotating axisymmetric torus and the Jacobi ellipsoid becomes a binary. The presence of viscosity is crucial since angular momentum needs to be redistributed for uniform rotation to be maintained. The phase transitions of the dynamical systems are briefly discussed in relation to previous numerical simulations of the formation and evolution of protostellar systems.

We have previously introduced the parameter `alpha' as an indicator of stability to m=2 nonaxisymmetric modes in rotating, self-gravitating, axisymmetric, gaseous and stellar systems. This parameter can be written as a function of the total rotational kinetic energy, the total gravitational potential energy, and as a function of the topology/connectedness and the geometric shape of a system. Here we extend the stability criterion to nonaxisymmetric equilibrium systems, such as ellipsoids and elliptical disks and cylinders. We test the validity of this extension by considering predictions for previously published, gaseous and stellar, nonaxisymmetric models. The above formulation and critical values account accurately for the stability properties of m=2 modes in gaseous Riemann S-type ellipsoids (including the Jacobi and Dedekind ellipsoids) and elliptical Riemann disks as well as in stellar elliptical Freeman disks and cylinders: all these systems are dynamically stable except for the stellar elliptical Freeman disks that exhibit a relatively small region of m=2 dynamical instability.

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.

We analyze previous results on the stability of uniformly and differentially rotating, self-gravitating, gaseous and stellar, axisymmetric systems to derive a new stability criterion for the appearance of toroidal, m=2 Intermediate (I) and bar modes. In the process, we demonstrate that the bar modes in stellar systems and the m=2 I-modes in gaseous systems have many common physical characteristics and only one substantial difference: because of the anisotropy of the stress tensor, dynamical instability sets in at lower rotation in stellar systems. This difference is reflected also in the new stability criterion. The new stability parameter "alpha" is formulated first for uniformly rotating systems and is based on the angular momentum content rather than on the energy content of a system. For stability of stellar systems "alpha" = 0.254-0.258, while it is in the range of 0.341-0.354 for gaseous systems. For uniform rotation, one can write "alpha" as a function of the total (rotational) kinetic and gravitational energies, and of a parameter which is characteristic of the topology/connectedness and the geometric shape of a system. "Alpha" can be extended to and calculated for a variety of differentially rotating, gaseous and stellar, axisymmetric disk and spheroidal models whose equilibrium structures and stability characteristics are known. We also estimate "alpha" for gaseous toroidal models and for stellar disk systems embedded in an inert or responsive "halo". We find that the new stability criterion holds equally well for all these previously published axisymmetric models.

A free-energy minimization approach is used to address the secular & dynamical instabilities & the bifurcations along sequences of rotating, self-gravitating fluid and stellar systems. Our approach stems from the Landau-Ginzburg theory of phase transitions. We focus on the Maclaurin sequence of oblate spheroids & on the effects of second-harmonic disturbances. Second-order phase transitions appear on the Maclaurin sequence also at the points of dynamical instability & of bifurcation of the Dedekind sequence. Distinguishing characteristic of each second-order phase transition is the (non)conservation of an integral of motion (e.g. circulation) which determines the appearance of the transition. Circulation is not conserved in stellar systems because the stress-tensor gradient terms that appear in the Jeans equations of motion include viscosity-like off-diagonal terms of the same order of magnitude as the conventional pressure gradient terms. This explains why the Jacobi bifurcation is a point of dynamical instability in stellar systems but only a point of secular instability in viscous fluids. The second-order phase transitions are discussed in relation to the dynamical instability of stellar systems, the lambda-transition of liquid He-4, the second-order phase transition in superconductivity & the mechanism of spontaneous symmetry breaking.