Dynamical stability and evolution of the discs of Sc galaxies

We examine the local stability of galactic discs against axisymmetric density perturbations with special attention to the different dynamics of the stellar and gaseous components. In particular, the discs of the Milky Way and of NGC 6946 are studied. The Milky Way is shown to be stable, whereas the inner parts of NGC 6946, a typical Sc galaxy from the Kennicutt sample, are dynamically unstable. The ensuing dynamical evolution of the composite disc is studied by numerical simulations. The evolution is so fierce that the stellar disc heats up dynamically on a short time-scale to such a high degree, which seems to contradict the morphological appearance of the galaxy. The star formation rate required to cool the disc dynamically is estimated. Even if the star formation rate in NGC 6946 is at present high enough to meet this requirement, it is argued that the discs of Sc galaxies cannot sustain such a high star formation rate for extended periods.     

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.     

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.     

The orientations of molecular clouds in the outer Galaxy: evidence for the scale of the turbulence driver?

Supernova (SN) explosions inject a considerable amount of energy into the interstellar medium (ISM) in regions with high-to-moderate star formation rates. In order to assess whether the driving of turbulence by supernovae is also important in the outer Galactic disc, where the star formation rates are lower, we study the spatial distribution of molecular cloud (MC) inclinations with respect to the Galactic plane. The latter contains important information on the nature of the mechanism of energy injection into the ISM. We analyse the spatial correlations between the position angles (PAs) of a selected sample of MCs (the largest clouds in the catalogue of the outer Galaxy published by Heyer et al). Our results show that when the PAs of the clouds are all mapped to values into the  [0°, 90°]  interval, there is a significant degree of spatial correlation between the PAs on spatial scales in the range of 100–800 pc. These scales are of the order of the sizes of individual SN shells in low-density environments such as those prevailing in the outer Galaxy and where the metallicity of the ambient gas is of the order of the solar value or smaller. These findings suggest that individual SN explosions, occurring in the outer regions of the Galaxy and in likewise spiral galaxies, albeit at lower rates, continue to play an important role in shaping the structure and dynamics of the ISM in those regions. The SN explosions we postulate here are likely associated with the existence of young stellar clusters in the far outer regions of the Galaxy and the ultraviolet emission and low levels of star formation observed with the Galaxy Evolution Explorer (GALEX) satellite in the outer regions of local galaxies.     

Chemical evolution with radial mixing

Models of the chemical evolution of our Galaxy are extended to include radial migration of stars and flow of gas through the disc. The models track the production of both iron and α-elements. A model is chosen that provides an excellent fit to the metallicity distribution of stars in the Geneva–Copenhagen survey (GCS) of the solar neighbourhood and a good fit to the local Hess diagram. The model provides a good fit to the distribution of GCS stars in the age–metallicity plane, although this plane was not used in the fitting process. Although this model's star formation rate is monotonically declining, its disc naturally splits into an α-enhanced thick disc and a normal thin disc. In particular, the model's distribution of stars in the ([O/Fe], [Fe/H]) plane resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's ridge line is entirely due to stellar migration, and there is the characteristic variation of stellar angular momentum along it that has been noted by Haywood in survey data. Radial mixing of stellar populations with high  σ _z   from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of  σ _z   with age than predicted by theory. The metallicity gradient in the interstellar medium is predicted to be steeper than in earlier models, but appears to be in good agreement with data for both our Galaxy and external galaxies. The models are inconsistent with a cut-off in the star formation rate at low gas surface densities. The absolute magnitude of the disc is given as a function of time in several photometric bands, and radial colour profiles are plotted for representative times.     

Mass models of the Milky Way

A parametrized model of the mass distribution within the Milky Way is fitted to the available observational constraints. The most important single parameter is the ratio of the scalelength R _d* of the stellar disc to R ₀. The disc and bulge dominate v _c( R ) at R ≲ R ₀ only for R _d,*/ R ₀≲0.3. Since the only knowledge we have of the halo derives from studies like the present one, we allow it to contribute to the density at all radii. When allowed this freedom, however, the halo causes changes in assumptions relating to R  ≪  R ₀ to affect profoundly the structure of the best-fitting model at R  ≫  R ₀. For example, changing the disc slightly from an exponential surface-density profile significantly changes the form of v _c( R ) at R  ≫  R ₀, where the disc makes a negligible contribution to v _c. Moreover, minor changes in the constraints can cause the halo to develop a deep hole at its centre that is not physically plausible. These problems call into question the proposition that flat rotation curves arise because galaxies have physically distinct haloes rather than outwards-increasing mass-to-light ratios.   The mass distribution of the Galaxy and the relative importance of its various components will remain very uncertain until more observational data can be used to constrain mass models. Data that constrain the Galactic force field at z ≳ R and at R  >  R ₀ are especially important.     

High-velocity stars from decay of small stellar systems

In this paper we present numerical results on the decay of small stellar systems under different initial conditions (multiplicity 3 ≤  N  ≤ 10, and various mass spectra, initial velocities and initial configurations). The numerical treatment uses the CHAIN1 code (Mikkola &38; Aarseth). Particular attention is paid to the distribution of high-velocity escapers: we define these as stars with velocity above 30 km s⁻¹. These numerical experiments confirm that small N -body systems are dynamically unstable and produce cascades of escapers in the process of their decay. It is shown that the fraction of stars that escape from small dense stellar systems with an escape velocity greater than 30 km s⁻¹ is ∼1 per cent for all systems treated here. This relatively small fraction must be considered in relation to the rate of star formation in the Galaxy in small groups: this could explain some moderately high-velocity stars observed in the Galactic disc and possibly some young stars with relatively high metallicity in the thick disc.     

Effect of angular momentum distribution on gravitational loss-cone instability in stellar clusters around a massive black hole

We study the small perturbations in spherical and thin disc stellar clusters surrounding a massive black hole. Because of the black hole, stars with sufficiently low angular momentum escape from the system through the loss cone. We show that the stability properties of spherical clusters crucially depend on whether the distribution of stars is monotonic or non-monotonic in angular momentum. It turns out that only non-monotonic distributions can be unstable. At the same time, instability in disc clusters is possible for both types of distribution.     

The local density of matter mapped by Hipparcos

We determine the velocity distribution and space density of a volume-complete sample of A and F stars, using parallaxes and proper motions from the Hipparcos satellite. We use these data to solve for the gravitational potential vertically in the local Galactic disc, by comparing the Hipparcos measured space density with predictions from various disc models. We derive an estimate of the local dynamical mass density of 0.102±0.010 pc⁻³, which may be compared with an estimate of 0.095 M_⊙ pc⁻³ in visible disc matter. Our estimate is found to be in reasonable agreement with other estimates by Crézé et al. and Pham, also based on Hipparcos data. We conclude that there is no compelling evidence for significant amounts of dark matter in the disc.     

Revisiting the relations: Galactic thin disc age–velocity dispersion relation

The velocity dispersion of stars in the solar neighbourhood thin disc increases with time after star formation. Nordström et al. performed the most recent observations to constrain the age–velocity dispersion relation. They fitted the age–velocity dispersion relations of each Galactic cardinal direction space velocity component, U (towards the Galactic Centre), V (in the direction of Galactic rotation) and W (towards the North Galactic Pole), with power laws and interpreted these as evidence for continuous heating of the disc in all directions throughout its lifetime. We revisit these relations with their data and use the results of Famaey et al. to show that structure in the local velocity distribution function distorts the in-plane ( U and V ) velocity distributions away from Gaussian so that a dispersion is not an adequate parametrization of their functions. The age–σ _W relation can however be constrained because the sample is well phase-mixed vertically. We do not find any local signature of the stellar warp in the Galactic disc. Vertical disc heating does not saturate at an early stage. Our new result is that a power law is not required by the data: disc heating models that saturate after ∼4.5 Gyr are equally consistent with observations.     

Simulations of minor mergers – II. The phase-space structure of thick discs

We analyse the phase-space structure of simulated thick discs that are the result of a 5:1 mass-ratio merger between a disc galaxy and a satellite. Our main goal is to establish what would be the imprints of a merger origin for the Galactic thick disc. We find that the spatial distribution predicted for thick-disc stars is asymmetric, seemingly in agreement with recent observations of the Milky Way thick disc. Near the Sun, the accreted stars are expected to rotate more slowly, to have broad velocity distributions and to occupy preferentially the wings of the line-of-sight velocity distributions. The majority of the stars in our model thick discs have low eccentricity orbits (in clear reference to the pre-existing heated disc) which give rise to a characteristic (sinusoidal) pattern for their line-of-sight velocities as a function of galactic longitude. The z -component of the angular momentum of thick-disc stars provides a clear discriminant between stars from the pre-existing disc and those from the satellite, particularly at large radii. These results are robust against the particular choices of initial conditions made in our simulations.     

Mapping the substructure in the Galactic halo with the next generation of astrometric satellites

We run numerical simulations of the disruption of satellite galaxies in a Galactic potential to build up the entire stellar halo, in order to investigate what the next generation of astrometric satellites will reveal by observing the halo of the Milky Way. We generate artificial DIVA , FAME and GAIA halo catalogues, in which we look for the signatures left by the accreted satellites. We develop a method based on the standard Friends-of-Friends algorithm applied to the space of integrals of motion. We find this simple method can recover about 50 per cent of the different accretion events, when the observational uncertainties expected for GAIA are taken into account, even when the exact form of the Galactic potential is unknown. The recovery rate for DIVA and FAME is much smaller, but these missions, like GAIA , should be able to test the hierarchical formation paradigm on our Galaxy by measuring the amount of halo substructure in the form of nearby kinematically cold streams with, for example, a two-point correlation function in velocity space.     

ECHA J0843.3–7905: Discovery of an ''old'' classical T Tauri star in the η Chamaeleontis cluster

Cold, dense clouds of gas have been proposed to explain the dark matter in Galactic haloes, and have also been invoked in the Galactic disc as an explanation for the excess faint submillimetre sources detected by SCUBA. Even if their dust-to-gas ratio is only a small percentage of that in conventional gas clouds, these dense systems would be opaque to visible radiation. We examine the possibility that the data sets of microlensing experiments searching for massive compact halo objects can also be used to search for occultation signatures by such clouds. We compute the rate and time-scale distribution of stellar transits by clouds in the Galactic disc and halo. We find that, for cloud parameters typically advocated by theoretical models, thousands of transit events should already exist within microlensing survey data sets. We examine the seasonal modulation in the rate caused by the Earth's orbital motion and find it provides an excellent probe of whether detected clouds are of disc or halo origin.     

Gravitational loss-cone instability in stellar systems with retrograde orbit precession

We study spherical and disc clusters in a near-Keplerian potential of galactic centres or massive black holes. In such a potential orbit precession is commonly retrograde, that is, the direction of the orbit precession is opposite to the orbital motion. It is assumed that stellar systems consist of nearly-radial orbits. We show that if there is a loss-cone at low angular momentum (e.g. due to consumption of stars by a black hole), an instability similar to loss-cone instability in plasma may occur. The gravitational loss-cone instability is expected to enhance black hole feeding rates. For spherical systems, the instability is possible for the number of spherical harmonics   l ≥ 3  . If there is some amount of counter-rotating stars in flattened systems, they generally exhibit the instability independent of azimuthal number m . The results are compared with those obtained recently by Tremaine for distribution functions monotonically increasing with angular momentum.
The analysis is based on simple characteristic equations describing small perturbations in a disc or a sphere of stellar orbits highly elongated in radius. These characteristic equations are derived from the linearized Vlasov equations (combining the collisionless Boltzmann kinetic equation and the Poisson equation), using the action-angle variables. We use two techniques for analysing the characteristic equations: the first one is based on preliminary finding of neutral modes, and the second one employs a counterpart of the plasma Penrose–Nyquist criterion for disc and spherical gravitational systems.     

Stellar populations in the Canis Major overdensity

We performed a photometric multicolour survey of the core of the Canis Major overdensity at     , reaching   V ∼ 22  and covering  0.3 × 1.0  arcmin². The main aim is to unravel the complex mixture of stellar populations toward this Galactic direction, where in the recent past important signatures of an accretion event have been claimed to be detected. While our previous investigations were based on disjointed pointings aimed at revealing the large-scale structure of the third Galactic Quadrant, we now focus on a complete coverage of a smaller field centred on the Canis Major overdensity. A large wavelength baseline, in the UBVRI bands, allows us to build up a suite of colour–colour and colour–magnitude diagrams, providing a much better diagnostic tool to disentangle the stellar populations of the region. In fact, the simple use of one colour–magnitude diagram, widely employed in all the previous studies defending the existence of the Canis Major galaxy, does not allow one to separate the effects of the different parameters (reddening, age, metallicity and distance) involved in the interpretation of data, forcing to rely on heavy modelling. In agreement with our previous studies, in the same general region of the Milky Way, we recognize a young stellar population compatible with the expected structure and extension of the Local (Orion) and Outer (Norma–Cygnus) spiral arms in the Third Galactic Quadrant. Moreover, we interpret the conspicuous intermediate-age metal-poor population as belonging to the Galactic thick disc, distorted by the effect of strong disc warping at this latitude, and to the Galactic halo.     

Measurement of the epicycle frequency in the Galactic disc and initial velocities of open clusters

A new method to measure the epicycle frequency κ in the Galactic disc is presented. We make use of the large data base on open clusters completed by our group to derive the observed velocity vector (amplitude and direction) of the clusters in the Galactic plane. In the epicycle approximation, this velocity is equal to the circular velocity given by the rotation curve, plus a residual or perturbation velocity, of which the direction rotates as a function of time with the frequency κ. Due to the non-random direction of the perturbation velocity at the birth time of the clusters, a plot of the present-day direction angle of this velocity as a function of the age of the clusters reveals systematic trends from which the epicycle frequency can be obtained. Our analysis considers that the Galactic potential is mainly axis-symmetric, or in other words, that the effect of the spiral arms on the Galactic orbits is small; in this sense, our results do not depend on any specific model of the spiral structure. The values of κ that we obtain provide constraints on the rotation velocity of the disc; in particular, V ₀ is found to be  230 ± 15 km s⁻¹  even if the short scale  ( R ₀= 7.5 kpc)  of the Galaxy is adopted. The measured κ at the solar radius is  43 ± 5 km s⁻¹ kpc ⁻¹  . The distribution of initial velocities of open clusters is discussed.     

Resolving the kinematic distance ambiguity of southern massive young stellar object candidates

We investigate the use of H  i data to resolve the near/far ambiguity in kinematic distances of massive young stellar object (MYSO) candidates. Kinematic distances were obtained from ¹³CO 1-0 (and N₂H⁺) spectral-line observations with the Mopra Telescope towards 94 candidates selected from the Red MSX Source (RMS) survey in the fourth Galactic quadrant  (282° < l < 350°)  . H  i data from the Southern Galactic Plane Survey (SGPS) were used in conjunction with the H  i self-absorption (SA) technique to determine the near or far distance. We resolved the kinematic distance ambiguity to 70 per cent of the sources. We can also simultaneously solve for any multiple line-of-sight component sources. We discuss the advantages and disadvantages of this technique in comparison with other methods, and also perform confidence checks on the reliability of using the H  i SA technique.
We examined the projected location of these sources in both the Galactic plane and longitude–velocity ( l – v ) diagrams to ascertain any recognizable spiral arm pattern. Although no obvious spiral pattern was found when compared to that proposed by Cordes and Lazio, far distance sources tended to lie on or near spiral arm loci. Near distance sources, however, had peculiar velocity uncertainties consistent with the separation between the spiral arms themselves. The l – v plot shows a more consistent picture, with tangent points of the spiral arm loci easily seen.
We conclude that using the H  i SA technique to determine kinematic distance ambiguities is a quick and reliable method in most cases, with an 80 per cent success rate in determining the correct designation of whether an object is at the near or far distance.     

The asymmetric structure of the Galactic halo

Using the stellar photometry catalogue based on the latest data release (DR4) of the Sloan Digital Sky Survey (SDSS), a study of the Galactic structure using star counts is carried out for selected areas of the sky. The sample areas are selected along a circle at a Galactic latitude of +60°, and 10 strips of high Galactic latitude along different longitudes. Direct statistics of the data show that the surface densities of ℓ from 180° to 360° are systematically higher than those of ℓ from 0° to 180°, defining a region of overdensity (in the direction of Virgo) and another one of underdensity (in the direction of Ursa Major) with respect to an axisymmetric model. It is shown by comparing the results from star counts in the ( g − r ) colour that the density deviations are due to an asymmetry of the stellar density in the halo. Theoretical models for the surface density profile are built and star counts are performed using a triaxial halo of which the parameters are constrained by observational data. Two possible reasons for the asymmetric structure are discussed.