Three-dimensional evolution of worms and chimneys in the Galactic disc

Three-dimensional modelling of the flow of gas and plasma in a section of the Galaxy has been carried out to study the evolution and formation of Galactic chimneys and worms. It is found that clustered supernovae located on either side of the Galactic plane are sources for the formation of well-collimated chimneys, having widths of ∼     . The thick gas disc may have a role in the collimation of chimneys. Channel maps of disc gas, obtained from the simulations, show the presence of sheet-like structures running perpendicular to the Galactic plane and resembling worms. Worms are believed to result from the break-up of the shells and supershells. However, the simulations show that although some worms correlate well with the debris of broken shells/supershells, others do not. They are cold gas that has been accelerated in the disc and rise on to the thick gas disc.     

Vertical distribution of Galactic disc stars and gas constrained by a molecular cloud complex

We investigate the dynamical effects of a molecular cloud complex with a mass ∼ 10⁷ M_⊙ and a size ∼ a few 100 pc on the vertical distribution of stars and atomic hydrogen gas in a spiral galactic disc. Such massive complexes have now been observed in a number of spiral galaxies. The extended mass distribution in a complex, with an average mass density 6 times higher than the Oort limit, is shown to dominate the local gravitational field. This results in a significant redistribution or clustering of the surrounding disc components towards the mid-plane, with a resulting decrease in their vertical scaleheights.
The modified, self-consistent stellar density distribution is obtained by solving the combined Poisson equation and the force equation along the z -direction for an isothermal stellar disc on which the complex is imposed. The effect of the complex is strongest at its centre, where the stellar mid-plane density increases by a factor of 2.6 and the vertical scaleheight decreases by a factor of 3.4 compared with the undisturbed stellar disc. A surprising result is the large radial distance of ∼ 500 pc from the complex centre over which the complex influences the disc; this is due to the extended mass distribution in a complex. The complex has a comparable effect on the vertical distribution of the atomic hydrogen gas in the galactic disc. This 'pinching' or constraining effect should be detectable in the nearby spiral galaxies, as for example has been done for NGC 2403 by Sicking. Thus the gravitational field of a complex results in local corrugations of the stellar and H  i vertical scaleheights, and the galactic disc potential is highly non-uniform on scales of the intercomplex separation of ∼ 1 kpc.     

Origin and structure of the Galactic disc(s)

We examine the chemical and dynamical structure in the solar neighbourhood of a model Galaxy that is the endpoint of a simulation of the chemical evolution of the Milky Way in the presence of radial mixing of stars and gas. Although the simulation's star formation rate declines monotonically from its unique peak and no merger or tidal event ever takes place, the model replicates all known properties of a thick disc, as well as matching special features of the local stellar population such as a metal-poor extension of the thin disc that has high rotational velocity. We divide the disc by chemistry and relate this dissection to observationally more convenient kinematic selection criteria. We conclude that the observed chemistry of the Galactic disc does not provide convincing evidence for a violent origin of the thick disc, as has been widely claimed.     

Encounters between spherical galaxies – I. Systems without a dark halo

We report here on a survey of N -body simulations of encounters between spherical galaxies. Initial systems are isotropic Jaffe models. Different sets of mass ratios, impact parameters and orbital energies are studied. Both merger remnants and systems perturbed after a non-merging encounter are analysed and compared to real-life elliptical galaxies. The properties of merger remnants show a large variety. Merger remnants resulting from head-on encounters are mainly non-rotating prolate spheroids. Merger remnants from models with   J _orb≠ 0  are tri-axial or mildly oblate spheroids, supported in part by rotation. The velocity distributions are biased towards the radial direction in the prolate case and the tangential direction in the oblate case. Non-mergers are affected in various ways, depending on the orbital characteristics. We conclude that many of the global properties of real-life ellipticals can, in principle, be attributed to a merger of spherical progenitors.     

Hydrostatic equilibrium in a magnetized, warped Galactic disc

Hydrostatic equilibrium of the multiphase interstellar medium in the solar vicinity is reconsidered, with the regular and turbulent magnetic fields treated separately. The regular magnetic field strength required to support the gas is consistent with independent estimates, provided that energy equipartition is maintained between turbulence and random magnetic fields. Our results indicate that a mid-plane value of B ₀=4 μG for the regular magnetic field near the Sun leads to more attractive models than B ₀=2 μG . The vertical profiles of both the regular and random magnetic fields contain disc and halo components, the parameters of which we have determined. The layer at 1≲| z |≲4 kpc can be overpressured and an outflow at a speed of about 50 km s⁻¹ may occur there, presumably associated with a Galactic fountain flow, if B ₀≃2 μG .
We show that hydrostatic equilibrium in a warped disc must produce asymmetric density distributions in z , in rough agreement with H  i observations in the outer Galaxy. This asymmetry may be a useful diagnostic of the details of the warping mechanism in the Milky Way and other galaxies. We find indications that gas and magnetic field pressures are different above and below the warped midplane in the outer Galaxy, and quantify the difference in terms of turbulent velocity and/or magnetic field strength.     

Satellite kinematics – I. A new method to constrain the halo mass–luminosity relation of central galaxies

Satellite kinematics can be used to probe the masses of dark matter haloes of central galaxies. In order to measure the kinematics with sufficient signal-to-noise ratio, one uses the satellite galaxies of a large number of central galaxies stacked according to similar properties (e.g. luminosity). However, in general, the relation between the luminosity of a central galaxy and the mass of its host halo will have non-zero scatter. Consequently, this stacking results in combining the kinematics of satellite galaxies in haloes of different masses, which complicates the interpretation of the data. In this paper, we present an analytical framework to model satellite kinematics, properly accounting for this scatter and for various selection effects. We show that in the presence of scatter in the halo mass–luminosity relation, the commonly used velocity dispersion of satellite galaxies can not be used to infer a unique halo mass–luminosity relation. In particular, we demonstrate that there is a degeneracy between the mean and the scatter of the halo mass–luminosity relation. We present a new technique that can break this degeneracy, and which involves measuring the velocity dispersions using two different weighting schemes: host weighting (each central galaxy gets the same weight) and satellite weighting (each central galaxy gets a weight proportional to its number of satellites). The ratio between the velocity dispersions obtained using these two weighting schemes is sensitive to the scatter in the halo mass–luminosity relation, and can thus be used to infer a unique relation between light and mass from the kinematics of satellite galaxies.     

Morphology, photometry and kinematics of N-body bars – I. Three models with different halo central concentrations

We discuss the morphology, photometry and kinematics of the bars which have formed in three N -body simulations. These have initially the same disc and the same halo-to-disc mass ratio, but their haloes have very different central concentrations. The third model includes a bulge. The bar in the model with the centrally concentrated halo (model MH) is much stronger, longer and thinner than the bar in the model with the less centrally concentrated halo (model MD). Its shape, when viewed side-on, evolves from boxy to peanut and then to 'X'-shaped, as opposed to that of model MD, which stays boxy. The projected density profiles obtained from cuts along the bar major axis, for both the face-on and the edge-on views, show a flat part, as opposed to those of model MD which are falling rapidly. A Fourier analysis of the face-on density distribution of model MH shows very large  m=2  , 4, 6 and 8 components. Contrary to this, for model MD the components  m=6  and 8 are negligible. The velocity field of model MH shows strong deviations from axial symmetry, and in particular has wavy isovelocities near the end of the bar when viewed along the bar minor axis. When viewed edge-on, it shows cylindrical rotation, which the MD model does not. The properties of the bar of the model with a bulge and a non-centrally concentrated halo (MDB) are intermediate between those of the bars of the other two models. All three models exhibit a lot of inflow of the disc material during their evolution, so that by the end of the simulations the disc dominates over the halo in the inner parts, even for model MH, for which the halo and disc contributions were initially comparable in that region.     

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.     

Ram-pressure stripping of disc galaxies orbiting in clusters – II. Galactic wakes

We present 3D hydrodynamical simulations of ram-pressure stripping of a disc galaxy orbiting in a galaxy cluster. In this paper, we focus on the properties of the galaxies' tails of stripped gas. The galactic wakes show a flaring width, where the flaring angle depends on the gas disc's cross-section with respect to the galaxy's direction of motion. The velocity in the wakes shows a significant turbulent component of a few     . The stripped gas is deposited in the cluster rather locally, i.e. within     from where it was stripped. We demonstrate that the most important quantity governing the tail density, length and gas mass distribution along the orbit is the galaxy's mass-loss per orbital length. This in turn depends on the ram pressure as well as the galaxy's orbital velocity.
For a sensitivity limit of     in projected gas density, we find typical tail lengths of     . Such long tails are seen even at large distances (0.5 to     ) from the cluster centre. At this sensitivity limit, the tails show little flaring, but a width similar to the gas disc's size.
Morphologically, we find good agreement with the H  i tails observed in the Virgo cluster by Chung et al. 2007 . However, the observed tails show a much smaller velocity width than predicted from the simulation. The few known X-ray and Hα tails are generally much narrower and much straighter than the tails in our simulations. Thus, additional physics like a viscous intracluster medium (ICM), the influence of cooling and tidal effects may be needed to explain the details of the observations.
We discuss the hydrodynamical drag as a heat source for the ICM but conclude that it is not likely to play an important role, especially not in stopping cooling flows.     

Structure and dynamics of galaxies with a low surface-brightness disc – I. The stellar and ionized-gas kinematics

Photometry and long-slit spectroscopy are presented for a sample of six galaxies with a low surface-brightness stellar disc and a bulge. The characterizing parameters of the bulge and disc components were derived by means of a two-dimensional photometric decomposition of the images of the sample galaxies. Their surface-brightness distribution was assumed to be the sum of the contribution of a Sérsic bulge and an exponential disc, with each component being described by elliptical and concentric isophotes of constant ellipticity and position angle. The stellar and ionized-gas kinematics were measured along the major and minor axes in half of the sample galaxies, whereas the other half was observed only along two diagonal axes. Spectra along two diagonal axes were obtained also for one of the objects with major and minor axis spectra. The kinematic measurements extend in the disc region out to a surface-brightness level  μ _R ≈ 24  mag arcsec⁻², reaching in all cases the flat part of the rotation curve. The stellar kinematics turns out to be more regular and symmetric than the ionized-gas kinematics, which often shows the presence of non-circular, off-plane and non-ordered motions. This raises the question about the reliability of the use of the ionized gas as the tracer of the circular velocity in the modelling of the mass distribution, in particular in the central regions of low surface-brightness galaxies.     

Galactic H i on the 50-au scale in the direction of three extragalactic sources observed with MERLIN

We present MERLIN observations of Galactic 21-cm H  i absorption at an angular resolution of  ∼0.1–0.2  arcsec and a velocity resolution of 0.5 km s⁻¹, in the direction of three moderately low latitude  (−8° < b < −12°)  extragalactic radio sources, 3C 111, 3C 123 and 3C 161, all of which are heavily reddened. H  i absorption is observed against resolved background emission sources up to ∼2 arcsec in extent and we distinguish details of the opacity distribution within 1–1.5 arcsec regions towards 3C 123 and 3C 161. This study is the second MERLIN investigation of small-scale structure in interstellar H  i (earlier work probed Galactic H  i in the directions of the compact sources 3C 138 and 3C 147). The 0.1-arcsec scale is intermediate between H  i absorption studies made with other fixed element interferometers with resolution of 1–10 arcsec and very long baseline interferometry studies with resolutions of 10–20 mas. At a scale of 1 arcsec (about 500 au), prominent changes in Galactic H  i opacity in excess of 1–1.5 are determined in the direction of 3C 161 with a signal-to-noise ratio of at least 10σ. Possible fluctuations in the H  i opacity at the level of about 1 are detected at the  2.5–3σ  level in the direction of 3C 123.