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.     

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.     

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.     

Vertical dynamics of disc galaxies in modified Newtonian dynamics

We investigate the possibility of discriminating between modified Newtonian dynamics (MOND) and Newtonian gravity with dark matter, by studying the vertical dynamics of disc galaxies. We consider models with the same circular velocity in the equatorial plane (purely baryonic discs in MOND and the same discs in Newtonian gravity embedded in spherical dark matter haloes), and we construct their intrinsic and projected kinematical fields by solving the Jeans equations under the assumption of a two-integral distribution function. We find that the vertical velocity dispersion of deep MOND discs can be much larger than in the equivalent spherical Newtonian models. However, in the more realistic case of high surface density discs, this effect is significantly reduced, casting doubt on the possibility of discriminating between MOND and Newtonian gravity with dark matter by using current observations.     

Stellar and gas dynamics of late-type barred-spiral galaxies: NGC 3359, a test case

We study the dynamics of a model for the late-type barred-spiral galaxy NGC 3359 by using both observational and numerical techniques. The results of our modelling are compared with photometric and kinematical data. The potential used is estimated directly from observations of the galaxy. It describes with a single potential function, a barred-spiral system with an extended spiral structure. Thus, the study of the dynamics in this potential has an interest by itself. We apply orbital theory and response models for the study of the stellar component, and smoothed particle hydrodynamics for modelling the gas. In particular, we examine the pattern speed of the system and the orbital character (chaotic or ordered) of the spiral arms. We conclude that the spiral pattern rotates slowly, in the sense that its corotation is close to or even beyond the end of the arms. Although a single, slow pattern speed could, under certain assumptions, characterize the whole disc, the comparison with the observational data indicates that probably the bar and the spirals have different angular velocities. In our two pattern speeds model, the best fit is obtained with a bar ending close to its 4:1 resonance and a more slowly rotating spiral. Assuming an 11 Mpc distance to the galaxy, a match of our models with the observed data indicates a pattern speed of about  39 km s⁻¹ kpc⁻¹  for the bar and about  15 km s⁻¹ kpc⁻¹  for the spiral. We do not find any indication for a chaotic character of the arms in this barred-spiral system. The flow in the region of the spirals can best be described as a regular 'precessing-ellipses flow'.     

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.     

Rotation and anisotropy of galaxies revisited

The use of the tensor virial theorem (TVT) as a diagnostic of anisotropic velocity distributions in galaxies is revisited. The TVT provides a rigorous global link between velocity anisotropy, rotation and shape, but the quantities appearing in it are not easily estimated observationally. Traditionally, use has been made of a centrally averaged velocity dispersion and the peak rotation velocity. Although this procedure cannot be rigorously justified, tests on model galaxies show that it works surprisingly well. With the advent of integral-field spectroscopy it is now possible to establish a rigorous connection between the TVT and observations. The TVT is reformulated in terms of sky-averages, and the new formulation is tested on model galaxies.     

The dynamics of S0 galaxies and their Tully–Fisher relation

This paper investigates the detailed dynamical properties of a relatively homogeneous sample of disc-dominated S0 galaxies, with a view to understanding their formation, evolution and structure. By using high signal-to-noise ratio long-slit spectra of edge-on systems, we have been able to reconstruct the complete line-of-sight velocity distributions of stars along the major axes of the galaxies. From these data, we have derived both model distribution functions (the phase density of their stars) and the approximate form of their gravitational potentials.
The derived distribution functions are all consistent with these galaxies being simple disc systems, with no evidence for a complex formation history. Essentially no correlation is found between the characteristic mass scalelengths and the photometric scalelengths in these galaxies, suggesting that they are dark-matter dominated even in their inner parts. Similarly, no correlation is found between the mass scalelengths and asymptotic rotation speed, implying a wide range of dark matter halo properties.
By comparing their asymptotic rotation speeds with their absolute magnitudes, we find that these S0 galaxies are systematically offset from the Tully–Fisher relation for later-type galaxies. The offset in luminosity is what one would expect if star formation had been suddenly switched off a few Gyr ago, consistent with a simple picture in which these S0s were created from ordinary later-type spirals which were stripped of their star-forming interstellar medium when they encountered a dense cluster environment.     

Determining the character of motion in quiet and active galaxies with a satellite companion

A galaxy model with a satellite companion is used to study the character of motion for stars moving in the x ‐y plane. It is observed that a large part of the phase plane is covered by chaotic orbits. The percentage of chaotic orbits increases when the galaxy has a dense nucleus of massM_n. The presence of the dense nucleus also increases the stellar velocities near the center of the galaxy. For small values of the distance R between the two bodies, low energy stars display a chaotic region near the centre of the galaxy, when the dense nucleus is present, while for larger values of R the motion in active galaxies is regular for low energy stars. Our results suggest that in galaxies with a satellite companion, the chaotic character of motion is not only a result of galactic interaction but also a result caused by the dense nucleus. Theoretical arguments are used to support the numerical outcomes. We follow the evolution of the galaxy, as mass is transported adiabatically from the disk to the nucleus. Our numerical results are in satisfactory agreement with observational data from M51‐type binary galaxies (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)