The flattening and the orbital structure of early-type galaxies and collisionless N-body binary disc mergers

We use oblate axisymmetric dynamical models including dark haloes to determine the orbital structure of intermediate mass to massive early-type galaxies in the Coma galaxy cluster. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 per cent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models, we show that flattening by stellar anisotropy maximizes the entropy for a given density distribution. Collisionless disc merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of discs unless the influence of dissipational processes was significant.     

Frequency map analysis of the orbital structure in elliptical galaxies

We present an application of the frequency map analysis to an elliptical galaxy which is represented by a generalization of a double-power-law spherical mass model. The density distribution of this model varies as r ^−γ close to the centre and as r ⁻⁴ at large radii. We study the case with γ = 1, which is known as the 'weak-cusp' model and which represents well the density profile of the 'core' galaxies observed by the Hubble Space Telescope . The final objective of our work is to improve our understanding of the dynamics of elliptical galaxies in a similar way to Merritt &38; Fridman, finding the regions of stochasticity, looking for resonances that might play an important role in sustaining the triaxial morphology, and analysing the diffusion of orbits. To this end, we use the frequency map analysis of Laskar, which has been applied widely in the field of celestial mechanics but which is a relatively new technique in the area of galactic dynamics. Finally, we show some useful features of this method in understanding the global dynamical structure of the system.     

Model-independent measurements of bar pattern speeds

The pattern speed is one of the fundamental parameters that determines the structure of barred galaxies. This quantity is usually derived from indirect methods or by employing model assumptions. The number of bar pattern speeds derived using the model-independent Tremaine & Weinberg technique is still very limited. We present the results of model-independent measurements of the bar pattern speed in four galaxies ranging in Hubble type from SB0 to SBbc. Three of the four galaxies in our sample are consistent with bars being fast rotators. The lack of slow bars is consistent with previous observations and suggests that barred galaxies do not have centrally concentrated dark matter haloes. This contradicts simulations of cosmological structure formation and observations of the central mass concentration in nonbarred galaxies.     

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.     

On the global structure of self-gravitating discs for softened gravity

The effects of gravitational softening on the global structure of self-gravitating discs in centrifugal equilibrium are examined in relation to hydrodynamical/gravitational simulations. The one-parameter spline softening proposed by Hernquist & Katz is used.
It is found that if the characteristic size of a disc, r , is comparable to or less than the gravitational softening length, ε, then the cross-section of the simulated disc is significantly larger than that of a no-softening (Newtonian) disc with the same mass and angular momentum.
We demonstrate, furthermore, that if r ≲ε/2 then the scaling relation r ∝ε^3/4 holds for a given mass and specific angular momentum distribution with mass. Finally, we compare some of the theoretical results obtained in this paper and a previous one with the results of numerical Tree-SPH simulations and find qualitative agreement.     

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.     

Response of the integrals in the Tremaine–Weinberg method to multiple pattern speeds: a counter-rotating inner bar in NGC 2950?

When integrals in the standard Tremaine–Weinberg method are evaluated for the case of a realistic model of a doubly barred galaxy, their modifications introduced by the second rotating pattern are in accord with what can be derived from a simple extension of that method, based on separation of tracer's density. This extension yields a qualitative argument that discriminates between prograde and retrograde inner bars. However, the estimate of the value of inner bar's pattern speed requires further assumptions. When this extension of the Tremaine–Weinberg method is applied to the recent observation of the doubly barred galaxy NGC 2950, it indicates that the inner bar there is counter-rotating, possibly with the pattern speed of  −140 ± 50 km s⁻¹ arcsec⁻¹  . The occurrence of counter-rotating inner bars can constrain theories of galaxy formation.     

Pattern speed of the stellar bar in NGC 7479

We study the pattern speed of the bar in NGC 7479 by comparing observations with numerical simulations of gas flow in a two-dimensional gravitational potential, derived from observations. The best agreement between the observations and the modelling is achieved for the fast bar pattern speed of 27 km s⁻¹ kpc⁻¹, when the corotation radius is at 50 arcsec, i.e. 1.1 times the radial length of the bar. This result is supported by the gas and dust lane morphologies, star formation distribution, projected velocity field and overall morphology. We find that star formation is most likely to be triggered close to the large-scale shocks and dust lanes in the bar. The net gas inflow rate in the simulations at 1-kpc radius is 4–6 M⊙ yr⁻¹ at intermediate times.     

Unravelling the mystery of the M31 bar

The inclination of M31 is too close to edge-on for a bar component to be easily recognized and is not sufficiently edge-on for a boxy/peanut bulge to protrude clearly out of the equatorial plane. Nevertheless, a sufficient number of clues allow us to argue that this galaxy is barred. We use fully self-consistent N -body simulations of barred galaxies and compare them with both photometric and kinematic observational data for M31. In particular, we rely on the near-infrared photometry presented in a companion paper. We compare isodensity contours to isophotal contours and the light profile along cuts parallel to the galaxy major axis and offset towards the north, or the south, to mass profiles along similar cuts on the model. All these comparisons, as well as position–velocity diagrams for the gaseous component, give us strong arguments that M31 is barred. We compare four fiducial N -body models to the data and thus set constraints on the parameters of the M31 bar, as its strength, length and orientation. Our 'best' models, although not meant to be exact models of M31, reproduce in a very satisfactory way the main relevant observations. We present arguments that M31 has both a classical and a boxy/peanut bulge. Its pseudo-ring-like structure at roughly 50 arcmin is near the outer Lindblad resonance of the bar and could thus be an outer ring, as often observed in barred galaxies. The shape of the isophotes also argues that the vertically thin part of the M31 bar extends considerably further out than its boxy bulge, that is, that the boxy bulge is only part of the bar, thus confirming predictions from orbital structure studies and from previous N -body simulations. It seems very likely that the backbone of M31's boxy bulge is families of periodic orbits, members of the x₁-tree and bifurcating from the x₁ family at its higher order vertical resonances, such as the x1v3 or x1v4 families.     

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.