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.     

On the relation between orbital structure and observed bar morphology

We investigate the morphological relation between the orbits of the central family of periodic orbits ( x ₁ family) and the bar itself using models of test particles moving in a barred potential. We show that different bar morphologies may have as a backbone the same set of x ₁ periodic orbits. We point out that by populating initially axisymmetric stellar discs exponentially with test particles in circular, or almost circular motion, we may end up with a response bar which reveals a shape different in crucial details from that of the individual stable x ₁ orbits. For example, a bar model in which the x ₁ orbits are pure ellipses may have a much more complicated response morphology. This depends on the particular invariant curves around x ₁, which are populated in each model.     

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

Order and chaos in spiral galaxies

Spiral galaxies contain both ordered and chaotic orbits. In normal spirals the perturbations are weak (of order 2–10%) and most orbits are ordered. The density wave theory refers mainly to linear perturbations. Nonlinear effects appear in the outer parts of the open spirals (S_b, S_c) and produce the termination of these spirals near the 4/1 resonance. On the other hand in barred spirals the perturbations are relatively strong (of order 100%). Then the outer spirals and the envelope of the bar are composed mainly of chaotic orbits, while the main body of the bar is composed of ordered orbits. The chaotic orbits of the spiral arms of strong barred galaxies are sticky, i.e. they do not escape from the galaxy for at least a Hubble time. The forms of these spirals are delineated by the unstable manifolds of the unstable periodic orbits L_1, L_2 near the ends of the bar and of other unstable periodic orbits inside and outside corotation.     

A constant dark matter halo surface density in galaxies

We confirm and extend the recent finding that the central surface density  μ_0D≡ r ₀ρ₀  of galaxy dark matter haloes, where r ₀ and  ρ₀  are the halo core radius and central density, is nearly constant and independent of galaxy luminosity. Based on the co-added rotation curves (RCs) of ∼1000 spiral galaxies, the mass models of individual dwarf irregular and spiral galaxies of late and early types with high-quality RCs, and the galaxy–galaxy weak-lensing signals from a sample of spiral and elliptical galaxies, we find that  log μ_0D= 2.15 ± 0.2  in units of  log(M_⊙ pc⁻²)  . We also show that the observed kinematics of Local Group dwarf spheroidal galaxies are consistent with this value. Our results are obtained for galactic systems spanning over 14 mag, belonging to different Hubble types and whose mass profiles have been determined by several independent methods. In the same objects, the approximate constancy of  μ_0D  is in sharp contrast to the systematical variations, by several orders of magnitude, of galaxy properties, including  ρ₀  and central stellar surface density.     

The shapes of galaxies in the Sloan Digital Sky Survey

We determine the underlying shapes of spiral and elliptical galaxies in the Sloan Digital Sky Survey Data Release 6 (SDSS DR6) from the observed distribution of projected galaxy shapes, taking into account the effects of dust extinction and reddening. We assume that the underlying shapes of spirals and ellipticals are well approximated by triaxial ellipsoids. The elliptical galaxy data are consistent with oblate spheroids, with a correlation between luminosity and ellipticity: the mean values of minor to middle axis ratios are 0.41 ± 0.03 for   M _r ≈−18  ellipticals and 0.76 ± 0.04 for   M _r ≈−22.5  ellipticals. Ellipticals show almost no dependence of axial ratio on galaxy colour, implying a negligible dust optical depth.
There is a strong variation of spiral galaxy shapes with colour indicating the presence of dust. The intrinsic shapes of spiral galaxies in the SDSS DR6 are consistent with flat discs with a mean and dispersion of thickness to diameter ratio of (21 ± 2) per cent, and a face-on ellipticity, e , of  ln( e ) =−2.33 ± 0.79  . Not including the effects of dust in the model leads to discs that are systematically rounder by up to 60 per cent. More luminous spiral galaxies tend to have thicker and rounder discs than lower luminosity spirals. Both elliptical and spiral galaxies tend to be rounder for larger galaxies.
The marginalized value of the edge-on r -band dust extinction E ₀ in spiral galaxies is   E ₀≃ 0.45  mag for galaxies of median colours, increasing to   E ₀= 1  mag for   g − r > 0.9  and   E ₀= 1.9  for the luminous and most compact galaxies, with half-light radii  <2  h ⁻¹ kpc  .     

Dynamics of the boxy elliptical galaxy NGC 1600

We use three-integral models to infer the distribution function (DF) of the boxy E3–E4 galaxy NGC 1600 from surface brightness and line-profile data on the minor and major axes. We assume axisymmetry and that the mass-to-light ratio is constant in the central ∼1 R _e. Stars in the resulting gravitational potential move mainly on regular orbits. We use an approximate third integral K from perturbation theory and write the DF as a sum of basis functions in the three integrals E , L _z and K . We then fit the projected moments of these basis functions to the kinematic observables and deprojected density, using a non-parametric algorithm. The deduced dynamical structure is radially anisotropic, with σ _θ σ _r ≈ σ _φ σ _r ≈0.7 on the major axis. Both on the minor axis and near the centre the velocity distribution is more isotropic; thus the model is flattened by equatorial radial orbits. The kinematic data are fitted without the need for a central black hole; the central mass determined previously from ground-based data therefore overestimates the actual black-hole mass. The mass-to-light ratio of the stars is M L _V =6  h ₅₀. The anisotropy structure of NGC 1600 with a radially anisotropic main body and more nearly isotropic centre is similar to that found recently in NGC 1399, 2434, 3379 and 6703, suggesting that this pattern may be common amongst massive elliptical galaxies. We discuss a possible merger origin of NGC 1600 in the light of these results.     

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.     

The dynamics of Neptune Trojan – I. The inclined orbits

The stability of Trojan type orbits around Neptune is studied. As the first part of our investigation, we present in this paper a global view of the stability of Trojans on inclined orbits. Using the frequency analysis method based on the fast Fourier transform technique, we construct high-resolution dynamical maps on the plane of initial semimajor axis a ₀ versus inclination i ₀. These maps show three most stable regions, with i ₀ in the range of  (0°, 12°), (22°, 36°)  and  (51°, 59°),  respectively, where the Trojans are most probably expected to be found. The similarity between the maps for the leading and trailing triangular Lagrange points L ₄ and L ₅ confirms the dynamical symmetry between these two points. By computing the power spectrum and the proper frequencies of the Trojan motion, we figure out the mechanisms that trigger chaos in the motion. The Kozai resonance found at high inclination varies the eccentricity and inclination of orbits, while the  ν₈  secular resonance around   i ₀∼ 44°  pumps up the eccentricity. Both mechanisms lead to eccentric orbits and encounters with Uranus that introduce strong perturbation and drive the objects away from the Trojan like orbits. This explains the clearance of Trojan at high inclination  (>60°)  and an unstable gap around  44°  on the dynamical map. An empirical theory is derived from the numerical results, with which the main secular resonances are located on the initial plane of  ( a ₀, i ₀)  . The fine structures in the dynamical maps can be explained by these secular resonances.     

The pattern speed of the bar in NGC 4596

The pattern speed is a defining parameter of any barred galaxy. A large number of model-dependent techniques have therefore been developed to derive the pattern speed. However, the only model-independent technique for measuring this quantity – the Tremaine–Weinberg method – has hitherto been applied to just one case, the SB0 galaxy NGC 936. In this paper, we apply the technique to a second system, the SBa galaxy NGC 4596. The resulting estimate for the pattern speed is Ω_p=52±13 km s⁻¹ kpc⁻¹. This result is corroborated by a spectrum obtained along the major axis of the bar in this system. The corotation radius associated with this pattern speed lies just beyond the end of the bar indicating a fast bar. Combining the bar major-axis spectra with data obtained from a Hubble Space Telescope WFPC2 image, we also find strong evidence for a nuclear disc.     

The STAGES view of red spirals and dusty red galaxies: mass-dependent quenching of star formation in cluster infall

We investigate the properties of optically passive spirals and dusty red galaxies in the A901/2 cluster complex at redshift ∼0.17 using rest-frame near-ultraviolet–optical spectral energy distributions, 24-μm infrared data and Hubble Space Telescope morphologies from the STAGES data set. The cluster sample is based on COMBO-17 redshifts with an rms precision of  σ _cz ≈ 2000 km s⁻¹  . We find that 'dusty red galaxies' and 'optically passive spirals' in A901/2 are largely the same phenomenon, and that they form stars at a substantial rate, which is only four times lower than that in blue spirals at fixed mass. This star formation is more obscured than in blue galaxies and its optical signatures are weak. They appear predominantly in the stellar mass range of  log  M _*/M_⊙=[10, 11]  where they constitute over half of the star-forming galaxies in the cluster; they are thus a vital ingredient for understanding the overall picture of star formation quenching in clusters. We find that the mean specific star formation rate (SFR) of star-forming galaxies in the cluster is clearly lower than in the field, in contrast to the specific SFR properties of blue galaxies alone, which appear similar in cluster and field. Such a rich red spiral population is best explained if quenching is a slow process and morphological transformation is delayed even more. At  log  M _*/M_⊙ < 10  , such galaxies are rare, suggesting that their quenching is fast and accompanied by morphological change. We note that edge-on spirals play a minor role; despite being dust reddened they form only a small fraction of spirals independent of environment.     

Contribution to the search for binaries among Am stars – VIII. New spectroscopic orbits of eight systems and statistical study of a sample of 91 Am stars

Red clump giant (RCG) stars can be used as distance indicators to trace the mass distribution of the Galactic bar. We use RCG stars from 44 bulge fields from the OGLE-II microlensing collaboration data base to constrain analytic triaxial models for the Galactic bar. We find the bar major-axis is oriented at an angle of 24°–27° to the Sun–Galactic Centre line-of-sight. The ratio of semimajor and semiminor bar axis scalelengths in the Galactic plane   x ₀, y ₀  , and vertical bar scalelength z ₀, is   x ₀ :  y ₀ :  z ₀= 10 : 3.5 : 2.6  , suggesting a slightly more prolate bar structure than the working model of Gerhard which gives the scalelength ratios as   x ₀ :  y ₀ :  z ₀= 10 : 4 : 3  .     

The pattern speed of the OH/IR stars in the Milky Way

We show how the continuity equation can be used to determine pattern speeds in the Milky Way Galaxy (MWG). This method, first discussed by Tremaine & Weinberg in the context of external galaxies, requires projected positions, ( l , b ), and line-of-sight velocities for a spatially complete sample of relaxed tracers. If the local standard of rest (LSR) has a zero velocity in the radial direction ( u _LSR), then the quantity that is measured is  Δ V ≡Ω_p R ₀- V _LSR  , where Ω_p is the pattern speed of the non-axisymmetric feature, R ₀ is the distance of the Sun from the Galactic centre and V _LSR is the tangential motion of the LSR, including the circular velocity. We use simple models to assess the reliability of the method for measuring a single, constant pattern speed of either a bar or spiral in the inner MWG. We then apply the method to the OH/IR stars in the ATCA/VLA OH 1612-MHz survey of Sevenster et al., finding  Δ V =252±41 km s^-1,  if   u _LSR=0  . Assuming further that   R ₀=8 kpc  and   V _LSR=220 km s^-1,  this gives  Ω_p=59±5 km s^-1 kpc^-1  with a possible systematic error of perhaps 10 km s⁻¹ kpc⁻¹. The non-axisymmetric feature for which we measure this pattern speed must be in the disc of the MWG.     

Polarization of synchrotron emission from relativistic reconfinement shocks

We study the polarization properties of relativistic reconfinement shocks with chaotic magnetic fields. Using our hydrodynamical model of their structure, we calculate synthetic polarization maps, longitudinal polarization profiles and discuss the spatially averaged polarization degree as a function of jet half-opening angle  Θ _j   , jet Lorentz factor  Γ _j   and observer inclination angle to the jet axis  θ_obs  . We find that for  θ_obs≲Θ _j   the wave electric vectors are parallel in the vicinity of the structure ends and perpendicular in between, while for  θ_obs > Θ _j   the polarization can only be perpendicular. The spatially averaged polarization degree does not exceed 30 per cent. Parallel average polarization, with polarization degrees lower than 10 per cent, has been found for  θ_obs < Θ _j   under the condition  Γ _j Θ _j > 1  . As earlier works predicted the parallel polarization from relativistic conical shocks, we explain our results by discussing conical shocks with divergent upstream flow.     

The radial Tully–Fisher relation for spiral galaxies – I

We find a new Tully–Fisher-like relation for spiral galaxies holding at different galactocentric radii. This radial Tully–Fisher relation allows us to investigate the distribution of matter in the optical regions of spiral galaxies. This relation, applied to three different samples of rotation curves of spiral galaxies, directly proves that: (i) the rotation velocity of spirals is a good measure of their gravitational potential and both the rotation curve's amplitudes and profiles are well predicted by galaxy luminosity, (ii) the existence of a dark component, less concentrated than the luminous one, and (iii) a scaling law, according to which, inside the disc optical size:   M _dark/ M _lum= 0.5( L _B /10¹¹ L _{B ⊙})^−0.7  .     

X-ray observations of the galaxy cluster PKS 0745−191: to the virial radius, and beyond

We measure X-ray emission from the outskirts of the cluster of galaxies PKS 0745−191 with Suzaku , determining radial profiles of density, temperature, entropy, gas fraction and mass. These measurements extend beyond the virial radius for the first time, providing new information about cluster assembly and the diffuse intracluster medium out to  ∼1.5  r ₂₀₀( r ₂₀₀≃ 1.7 Mpc ≃ 15 arcmin  ). The temperature is found to decrease by roughly 70 per cent from 0.3 to  1 r ₂₀₀  . We also see a flattening of the entropy profile near the virial radius and consider the implications this has for the assumption of hydrostatic equilibrium when deriving mass estimates. We place these observations in the context of simulations and analytical models to develop a better understanding of non-gravitational physics in the outskirts of the cluster.     

Gas motions in the plane of the spiral galaxy NGC 3631

The velocity field of the nearly face-on galaxy NGC 3631, derived from observations in the H α line and H  i radio line, is analysed to study perturbations related to the spiral structure of the galaxy. We confirm our previous conclusion that the line-of-sight velocity field gives evidence of the wave nature of the observed two-armed spiral structure. Fourier analysis of the observed velocity field is used to determine the location of corotation of the spiral structure of this galaxy, and the radius of corotation R _c is found to be about 42 arcsec, or 3.2 kpc. The vector velocity field of the gas in the plane of the disc is restored, and, taking into account that we previously investigated vertical motions, we now have a full three-dimensional gaseous velocity field of the galaxy. We show clear evidence of the existence of two anticyclonic and four cyclonic vortices near corotation in a frame of reference rotating with the spiral pattern. The centres of the anticyclones lie between the observed spiral arms. The cyclones lie close to the observed spirals, but their centres are shifted from the maxima in brightness.     

Survival and disruption of subsystems during a cold collapse

Cold collapse of a cluster composed of small identical clumps, each of which is in virial equilibrium, is considered. Since the clumps have no relative motion with respect to each other initially, the cluster collapses under its own gravity. At the first collapse of the cluster, most of the clumps are destroyed, but some survive. In order to find the condition for the clumps to survive, we made a systematic study in two-parameter space: the number of the clumps N _c and the size of the clump r _v . We obtained the condition N _c ≫ 1 and n _k  ≥ 1, where n _k is related to r _v and the initial radius of the cluster R _ini through the relation R _ini/ r _v  = 2 N ^{( n k +5)/6}_c. A simple analytical argument supports the numerical result. This n _k corresponds to the index of the power spectrum of the density fluctuation in the cosmological hierarchical clustering, and thus our result may suggest that in the systems smaller than 2/Ω h ²)Mpc, the first violent collapse is strong enough to sweep away all the substructures that exist before the collapse.