The coalescence of invariant manifolds and the spiral structure of barred galaxies

In a previous paper (Voglis et al., Paper I), we demonstrated that, in a rotating galaxy with a strong bar, the unstable asymptotic manifolds of the short-period family of unstable periodic orbits around the Lagrangian points L ₁ or L ₂ create correlations among the apocentric positions of many chaotic orbits, thus supporting a spiral structure beyond the bar. In this paper, we present evidence that the unstable manifolds of all the families of unstable periodic orbits near and beyond corotation contribute to the same phenomenon. Our results refer to a N -body simulation, a number of drawbacks of which, as well as the reasons why these do not significantly affect the main results, are discussed. We explain the dynamical importance of the invariant manifolds as due to the fact that they produce a phenomenon of 'stickiness' slowing down the rate of chaotic escape in an otherwise non-compact region of the phase space. We find a stickiness time of the order of 100 dynamical periods, which is sufficient to support a long-living spiral structure. Manifolds of different families become important at different ranges of values of the Jacobi constant. The projections of the manifolds of all the different families in the configuration space produce a pattern due to the 'coalescence' of the invariant manifolds. This follows closely the maxima of the observed   m = 2  component near and beyond corotation. Thus, the manifolds support both the outer edge of the bar and the spiral arms.     

Gas-driven evolution of stellar orbits in barred galaxies

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 disc (stars+gas) and a live halo. The results are compared with a pure stellar model, subject to nearly identical initial conditions, which are chosen so 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 disc. 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 disc 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 centre. 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 centre. Milder instability results in a smaller bulge, the basic parameters of which are in agreement with observations. We also find that the overall evolution in the model with a gas component is accelerated because of the larger central mass concentration and the resulting decrease in the characteristic dynamical time.     

Model-based pattern speed estimates for 38 barred galaxies

We have modelled 38 barred galaxies by using near-infrared and optical data from the Ohio State University Bright Spiral Galaxy Survey. We constructed the gravitational potentials of the galaxies from H -band photometry, assuming a constant mass-to-light ratio. The halo component we choose corresponds to the so-called universal rotation curve. In each case, we used the response of gaseous and stellar particle disc to rigidly rotating potential to determine the pattern speed.
We find that the pattern speed of the bar depends roughly on the morphological type. The average value of corotation resonance radius to bar radius,     , increases from 1.15 ± 0.25 in types SB0/a–SBab to 1.44 ± 0.29 in SBb and 1.82 ± 0.63 in SBbc–SBc. Within the error estimates for the pattern speed and bar radius, all galaxies of type SBab or earlier have a fast bar     , whereas the bars in later type galaxies include both fast and slow rotators. Of 16 later type galaxies with a nominal value of     , there are five cases, where the fast-rotating bar is ruled out by the adopted error estimates.
We also study the correlation between the parameter     and other galactic properties. The clearest correlation is with the bar size: the slowest bars are also the shortest bars when compared to the galaxy size. A weaker correlation is seen with bar strength in a sense that slow bars tend to be weaker. These correlations leave room for a possibility that the determined pattern speed in many galaxies corresponds to actually that of the spiral, which rotates more slowly than the bar. No clear correlation is seen with either the galaxy luminosity or the colour.     

Numerical simulations of interacting gas-rich barred galaxies: vertical impact of small companions

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

A bar in the inner halo of barred galaxies – I. Structure and kinematics of a representative model

N -body simulations argue that the inner haloes of barred galaxies should not be spherical, nor even axisymmetric, but triaxial. The departure from sphericity is the strongest near the centre and decreases outwards; typical axial ratios for the innermost parts are of the order of 0.8. The halo shape is prolate-like in the inner parts up to a certain radius and then turns to oblate-like. I call this inner halo structure the 'halo bar' and analyse here in depth its structure and kinematics in a representative model. It is always considerably shorter than the disc bar. It lags the disc bar by only a few degrees at all radii and the difference between the two bar phases increases with distance from the centre. The two bars turn with roughly the same pattern speed. This means that the halo bar is a slow bar, since its corotation radius is much larger than its length. The bisymmetric component in the halo continues well outside the halo bar in the form of an open spiral, trailing behind the disc bar. The inner parts of the halo display some mean rotation in the same sense as the disc rotation. This is more important for particles nearer to the equatorial plane and decreases with increasing distance from it, but is always much smaller than the disc rotation.     

Invariant manifolds, phase correlations of chaotic orbits and the spiral structure of galaxies

In the presence of a strong   m = 2  component in a rotating galaxy, the phase-space structure near corotation is shaped to a large extent by the invariant manifolds of the short-period family of unstable periodic orbits terminating at L ₁ or L ₂. The main effect of these manifolds is to create robust phase correlations among a number of chaotic orbits large enough to support a spiral density wave outside corotation. The phenomenon is described theoretically by soliton-like solutions of a Sine–Gordon equation. Numerical examples are given in an N -body simulation of a barred spiral galaxy. In these examples, we demonstrate how the projection of unstable manifolds in configuration space reproduces essentially the entire observed bar–spiral pattern.     

Orbital structure in N-body models of barred-spiral galaxies

We study the orbital structure in a series of self-consistent N -body configurations simulating rotating barred galaxies with spiral and ring structures. We perform frequency analysis in order to measure the angular and the radial frequencies of the orbits at two different time snapshots during the evolution of each N -body system. The analysis is done separately for the regular and the chaotic orbits. We thereby identify the various types of orbits, determine the shape and percentages of the orbits supporting the bar and the ring/spiral structures, and study how the latter quantities change during the secular evolution of each system. Although the frequency maps of the chaotic orbits are scattered, we can still identify concentrations around resonances. We give the distributions of frequencies of the most important populations of orbits. We explore the phase-space structure of each system using projections of the 4D surfaces of section. These are obtained via the numerical integration not only of the orbits of test particles, but also of the real N -body particles. We thus identify which domains of the phase space are preferred and which are avoided by the real particles. The chaotic orbits are found to play a major role in supporting the shape of the outer envelope of the bar as well as the rings and the spiral arms formed outside corotation.     

Rotation curves of luminous spiral galaxies

We have investigated the stellar light distribution and the rotation curves of high‐luminosity spiral galaxies in the local Universe. The sample contains 30 high‐quality extended Hα and H I rotation curves. The stellar disk scale‐length of these objects was measured or taken from the literature. We find that in the outermost parts of the stellar disks of these massive objects, the rotation curves agree with the Universal Rotation Curve (Salucci et al. 2007), however a few rotation curves of the sample show a divergence (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application of the global modal approach to the spiral galaxies

We have tested the applicability of the global modal approach in the density wave theory of spiral structure for a sample of spiral galaxies with measured axisymmetric background properties. We report here the results of the simulations for four galaxies: NGC 488, NGC 628, NGC 1566, and NGC 3938. Using the observed radial distributions for the stellar velocity dispersions and the rotation velocities we have constructed the equilibrium models for the galactic disks in each galaxy and implemented two kinds of stability analyses - the linear global analysis and 2D-nonlinear simulations. In general, the global modal approach is able to reproduce the observed properties of the spiral arms in the galactic disks. The growth of spirals in the galactic disks can be physically understood in terms of amplification by over-reflection at the corotation resonance. Our results support the global modal approach as a theoretical explanation of spiral structure in galaxies. This revised version was published online in August 2006 with corrections to the Cover Date.