The Radial Velocity Experiment (RAVE) is a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A). The 6th and final data release (DR6 or FDR) is based on 518387 observations of 451783 unique stars. RAVE observations were taken between 12 April 2003 and 4 April 2013. Here we present the genesis, setup and data reduction of RAVE as well as wavelength-calibrated and flux-normalized spectra and error spectra for all observations in RAVE DR6. Furthermore, we present derived spectral classification and radial velocities for the RAVE targets, complemented by cross matches with Gaia DR2 and other relevant catalogs. A comparison between internal error estimates, variances derived from stars with more than one observing epoch and a comparison with radial velocities of Gaia DR2 reveals consistently that 68% of the objects have a velocity accuracy better than 1.4 km/s, while 95% of the objects have radial velocities better than 4.0 km/s. Stellar atmospheric parameters, abundances and distances are presented in subsequent publication. The data can be accessed via the RAVE Web (http://rave-survey.org) or the Vizier database.

We present part 2 of the 6th and final Data Release (DR6 or FDR) of the Radial Velocity Experiment (RAVE), a magnitude-limited (9<I<12) spectroscopic survey of Galactic stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra (R~7500) cover the Ca-triplet region (8410-8795A) and span the complete time frame from the start of RAVE observations on 12 April 2003 to their completion on 4 April 2013. In the second of two publications, we present the data products derived from 518387 observations of 451783 unique stars using a suite of advanced reduction pipelines focussing on stellar atmospheric parameters, in particular purely spectroscopically derived stellar atmospheric parameters (Teff, log(g), and the overall metallicity), enhanced stellar atmospheric parameters inferred via a Bayesian pipeline using Gaia DR2 astrometric priors, and asteroseismically calibrated stellar atmospheric parameters for giant stars based on asteroseismic observations for 699 K2 stars. In addition, we provide abundances of the elements Fe, Al, and Ni, as well as an overall [alpha/Fe] ratio obtained using a new pipeline based on the GAUGUIN optimization method that is able to deal with variable signal-to-noise ratios. The RAVE DR6 catalogs are cross matched with relevant astrometric and photometric catalogs, and are complemented by orbital parameters and effective temperatures based on the infrared flux method. The data can be accessed via the RAVE Web site (http://rave-survey.org) or the Vizier database.

The conventional approach to orbit trapping at Lindblad resonances via a pendulum equation fails when the parent of the trapped orbits is too circular. The problem is explained and resolved in the context of the Torus Mapper and a realistic Galaxy model. Tori are computed for orbits trapped at both the inner and outer Lindblad resonances of our Galaxy. At the outer Lindblad resonance, orbits are quasiperiodic and can be accurately fitted by torus mapping. At the inner Lindblad resonance, orbits are significantly chaotic although far from ergodic, and each orbit explores a small range of tori obtained by torus mapping.

Torus mapping yields constants of motion for stars trapped at a resonance. Each such constant of motion yields a system of contours in velocity space at the Sun and neighbouring points. If Jeans' theorem applied to resonantly trapped orbits, the density of stars in velocity space would be equal at all intersections of any two contours. A quantitative measure of the violation of this principal is defined and used to assess various pattern speeds for a model of the bar recently fitted to observations of interstellar gas. Trapping at corotation of a bar with pattern speed in the range 33-36 /Gyr is favoured and trapping at the outer Lindblad resonance is disfavoured. As one moves around the Sun the structure of velocity space varies quite rapidly, both as regards the observed star density and the zones of trapped orbits. The data seem consistent with trapping at corotation.

Dynamical models will be key to exploitation of the incoming flood data for our Galaxy. Modelling techniques are reviewed with an emphasis on f(J) modelling.

We address an issue that arises when self-consistently flattened dynamical stellar systems are constructed by adopting a distribution function (DF) that depends on the action integrals. The velocity distribution at points on the symmetry axis is controlled by the the dependence of the DF on just one action, while at points off the symmetry axis two actions are involved. Consequently, the physical requirement that the velocity distribution evolves continuously in the neighbourhood of the symmetry axis restricts the functional forms of acceptable DFs. An algorithm for conforming to this restriction is presented and used to construct a variety of flattened models.

10,000 simulations of 1000-particle realisations of the same cluster are computed by direct force summation. Over three crossing times the original Poisson noise is amplified more than tenfold by self-gravity. The cluster's fundamental dipole mode is strongly excited by Poisson noise, and this mode makes a major contribution to driving diffusion of stars in energy. The diffusive flow through action space is computed for the simulations and compared with the predictions of both local-scattering theory and the Balescu-Lenard (BL) equation. The predictions of local-scattering theory are qualitatively wrong because the latter neglects self-gravity. These results imply that local-scattering theory is of little value. Future work on cluster evolution should employ either N-body simulation or the BL equation. However, significant code development will be required to make use of the BL equation practicable.

The principal results of the classic analysis of the shearing sheet and swing amplification by Julian & Toomre (1966) are re-derived in a more accessible way and then used to gain a better quantitative understanding of the dynamics of stellar discs. The axisymmetric limit of the shearing sheet is derived and used to re-derive Kalnajs' 1965 dispersion relation and Toomre's 1964 stability criterion for axisymmetric disturbances. Using the shearing sheet to revisit Toomre's important 1969 paper on the group velocity implied by Lin-Shu-Kalnajs dispersion relation, we discover that two rather than one wavepackets emerges inside corotation: one each side of the inner Lindblad resonance. Although LSK dispersion relation provides useful interpretations of both wavepackets, the shearing sheet highlights the limitations of the LSK approach to disc dynamics. Disturbances by no means avoid an annulus around corotation, as the LSK dispersion relation implies. While disturbances of the shearing sheet have a limited life in real space, they live on much longer in velocity space, which Gaia allows us to probe extensively. C++ code is provided to facilitate applications of winding spiral waves.

A new family of self-consistent DF-based models of stellar systems is explored. The stellar component of the models is described by a distribution function (DF) depending on the action integrals, previously used to model the Fornax dwarf spheroidal galaxy (dSph). The stellar component may cohabit with either a dark halo, also described by a DF, or with a massive central black hole. In all cases we solve for the model's self-consistent potential. Focussing on spherically symmetric models, we show how the stellar observables vary with the anisotropy prescribed by the DF, with the dominance and nature of the dark halo, and with the mass of the black hole. We show that precise fits to the observed surface brightness profiles of four globular clusters can be obtained for a wide range of prescribed velocity anisotropies. We also obtain precise fits to the observed projected densities of four dSphs. Finally, we present a three-component model of the Scupltor dSph with distinct DFs for the red and blue horizontal branch stars and the dark matter halo.

We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zV_z phase plane, for both V_phi and V_R, and recover the remarkable "phase spiral" discovered by Gaia. We identify the anticipated quadrupole signature in zV_z of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the alpha-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zV_z phase plane lasts >2 Gyr but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited < 0.5 Gyr ago by an object like Sgr with total mass 3 x 10^10 Msun (stripped down from 5 x 10^10 Msun when it first entered the halo) passing through the plane.

Fitting equilibrium dynamical models to observational data is an essential step in understanding the structure of the gaseous hot haloes that surround our own and other galaxies. However, the two main categories of models that are used in the literature are poorly suited for this task: (i) simple barotropic models are analytic and can therefore be adjusted to match the observations, but are clearly unrealistic because the rotational velocity $v_\phi(R,z)$ does not depend on the distance $z$ from the galactic plane, while (ii) models obtained as a result of cosmological galaxy formation simulations are more realistic, but are impractical to fit to observations due to high computational cost. Here we bridge this gap by presenting a general method to construct axisymmetric baroclinic equilibrium models of rotating galactic coronae in arbitrary external potentials. We consider in particular a family of models whose equipressure surfaces in the $(R,z)$ plane are ellipses of varying axis ratio. These models are defined by two one-dimensional functions, the axial ratio of pressure $q_{\rm axis}(z)$ and the value of the pressure $P_{\rm axis}(z)$ along the galaxy's symmetry axis. These models can have a rotation speed $v_\phi(R,z)$ that realistically decreases as one moves away from the galactic plane, and can reproduce the angular momentum distribution found in cosmological simulations. The models are computationally cheap to construct and can thus be used in fitting algorithms. We provide a python code that given $q_{\rm axis}(z)$, $P_{\rm axis}(z)$ and $\Phi(R,z)$ returns $\rho(R,z)$, $T(R,z)$, $P(R,z)$, $v_\phi(R,z)$. We show a few examples of these models using the Milky Way as a case study.

The classical theory of cluster relaxation is unsatisfactory because it involves the Coulomb logarithm. The Balescu-Lenard (BL) equation provides a rigorous alternative that has no ill-defined parameter. Moreover, the BL equation, unlike classical theory, includes the cluster's self-gravity. A heuristic argument is given that indicates that relaxation does not occur predominantly through two-particle scattering and is enhanced by self-gravity. The BL equation is adapted to a spherical system and used to estimate the flux through the action space of isochrone clusters with different velocity anisotropies. A range of fairly different secular behaviours is found depending on the fraction of radial orbits. Classical theory is also used to compute the corresponding classical fluxes. The BL and classical fluxes are very different because (a) the classical theory materially under-estimates the impact of large-scale collectively amplified fluctuations and (b) only the leading terms in an infinite sum for the BL flux are computed. A complete theory of cluster relaxation likely requires that the sum in the BL equation be decomposed into a sum over a finite number of small wavenumbers complemented by an integral over large wavenumbers analogous to classical theory.

A simple model is presented of the formation of the spiral the (z,v_z) phase plane of solar-neighbourhood stars that was recently discovered in Gaia data. The key is that the frequency Omega_z at which stars oscillate vertically depends on angular momentum about the z axis in addition to the amplitude of the star's vertical oscillations. Spirals should form in both <v_phi> and <v_R> whenever a massive substructure, such as the Sgr dwarf galaxy, passes through the Galactic plane. The model yields similar spirals to those observed in both <v_phi> and <v_R>. The primary driver is the component of the tidal force that lies in the plane. We investigate the longevity of the spirals and the mass of the substructure, but the approximations inherent in the model make quantitative results unreliable. The work relies heavily on a self-consistent, multi-component model of our Galaxy produced by the AGAMA package for f(J) modelling.

We present new dynamical models of dwarf spheroidal galaxies (dSphs) in which both the stellar component and the dark halo are described by analytic distribution functions that depend on the action integrals. In their most general form these distribution functions can represent axisymmetric and possibly rotating stellar systems. Here, as a first application, we model the Fornax dSph, limiting ourselves, for simplicity, to the non rotating, spherical case. The models are compared with state-of-the-art spectroscopic and photometric observations of Fornax, exploiting the knowledge of the line-of-sight velocity distribution of the models and accounting for the foreground contamination from the Milky Way. The model that best fits the structural and kinematic properties of Fornax has a cored dark halo, with core size $r_{\rm c}\simeq1.03$ kpc. The dark-to-luminous mass ratio is $(M_{\rm dm}/M_{\star})|_{R_{\rm eff}}\simeq9.6$ within the effective radius $R_{\rm eff} \simeq 0.62\,$kpc and $(M_{\rm dm}/M_{\star})|_{3 {\rm kpc}} \simeq 144$ within 3 kpc. The stellar velocity distribution is isotropic almost over the full radial range covered by the spectroscopic data and slightly radially anisotropic in the outskirts of the stellar distribution. The dark-matter annihilation $J$-factor and decay $D$-factor are, respectively, $\log_{10}(J$ $[$GeV$^2$ cm$^{-5}])\simeq18.34$ and $\log_{10}(D$ $[$GeV cm$^{-2}])\simeq18.55$, for integration angle $\theta = 0.5^{\circ}$. This cored halo model of Fornax is preferred, with high statistical significance, to both models with a Navarro, Frenk and White dark halo and simple mass-follows-light models.

We combine parallaxes from the first Gaia data release with the spectrophotometric distance estimation framework for stars in the fifth RAVE survey data release. The combined distance estimates are more accurate than either determination in isolation - uncertainties are on average two times smaller than for RAVE-only distances (three times smaller for dwarfs), and 1.4 times smaller than TGAS parallax uncertainties (two times smaller for giants). We are also able to compare the estimates from spectrophotometry to those from Gaia, and use this to assess the reliability of both catalogues and improve our distance estimates. We find that the distances to the lowest log g stars are, on average, overestimated and caution that they may not be reliable. We also find that it is likely that the Gaia random uncertainties are smaller than the reported values. As a byproduct we derive ages for the RAVE stars, many with relative uncertainties less than 20 percent. These results for 219566 RAVE sources have been made publicly available, and we encourage their use for studies that combine the radial velocities provided by RAVE with the proper motions provided by Gaia. A sample that we believe to be reliable can be found by taking only the stars with the flag notification flag_any=0

We explore the connections between stellar age, chemistry, and kinematics across a Galactocentric distance of $7.5 < R\,(\mathrm{kpc}) < 9.0$, using a sample of $\sim 12\,000$ intermediate-mass (FGK) turnoff stars observed with the RAdial Velocity Experiment (RAVE) survey. The kinematics of this sample are determined using radial velocity measurements from RAVE, and parallax and proper motion measurements from the Tycho-Gaia Astrometric Solution (TGAS). In addition, ages for RAVE stars are determined using a Bayesian method, taking TGAS parallaxes as a prior. We divide our sample into young ($0 < \tau < 3$ Gyr) and old ($8 < \tau < 13$ Gyr) populations, and then consider different metallicity bins for each of these age groups. We find significant differences in kinematic trends of young and old, metal-poor and metal-rich, stellar populations. In particular, we find a strong metallicity dependence in the mean Galactocentric radial velocity as a function of radius ($\partial {V_{\rm R}}/\partial R$) for young stars, with metal-rich stars having a much steeper gradient than metal-poor stars. For $\partial {V_{\phi}}/\partial R$, young, metal-rich stars significantly lag the LSR with a slightly positive gradient, while metal-poor stars show a negative gradient above the LSR. We interpret these findings as correlations between metallicity and the relative contributions of the non-axisymmetries in the Galactic gravitational potential (the spiral arms and the bar) to perturb stellar orbits.

We use data from the Radial Velocity Experiment (RAVE) and the Tycho-Gaia astrometric solution catalogue (TGAS) to compute the velocity fields yielded by the radial (VR), azimuthal (Vphi) and vertical (Vz) components of associated Galactocentric velocity. We search in particular for variation in all three velocity components with distance above and below the disc midplane, as well as how each component of Vz (line-of-sight and tangential velocity projections) modifies the obtained vertical structure. To study the dependence of velocity on proper motion and distance we use two main samples: a RAVE sample including proper motions from the Tycho-2, PPMXL and UCAC4 catalogues, and a RAVE-TGAS sample with inferred distances and proper motions from the TGAS and UCAC5 catalogues. In both samples, we identify asymmetries in VR and Vz. Below the plane we find the largest radial gradient to be dVR / dR = -7.01+- 0.61 km\s kpc, in agreement with recent studies. Above the plane we find a similar gradient with dVR / dR= -9.42+- 1.77 km\s kpc. By comparing our results with previous studies, we find that the structure in Vz is strongly dependent on the adopted proper motions. Using the Galaxia Milky Way model, we demonstrate that distance uncertainties can create artificial wave-like patterns. In contrast to previous suggestions of a breathing mode seen in RAVE data, our results support a combination of bending and breathing modes, likely generated by a combination of external or internal and external mechanisms.

It has been known for more than thirty years that the distribution of molecular gas in the innermost 300 parsecs of the Milky Way, the Central Molecular Zone, is strongly asymmetric. Indeed, approximately three quarters of molecular emission comes from positive longitudes, and only one quarter from negative longitudes. However, despite much theoretical effort, the origin of this asymmetry has remained a mystery. Here we show that the asymmetry can be neatly explained by unsteady flow of gas in a barred potential. We use high-resolution 3D hydrodynamical simulations coupled to a state-of-the-art chemical network. Despite the initial conditions and the bar potential being point-symmetric with respect to the Galactic Centre, asymmetries develop spontaneously due to the combination of a hydrodynamical instability known as the "wiggle instability" and the thermal instability. The observed asymmetry must be transient: observations made tens of megayears in the past or in the future would often show an asymmetry in the opposite sense. Fluctuations of amplitude comparable to the observed asymmetry occur for a large fraction of the time in our simulations, and suggest that the present is not an exceptional moment in the life of our Galaxy.

We run hydrodynamical simulations of a 2D isothermal non self-gravitating inviscid gas flowing in a rigidly rotating externally imposed potential formed by only two components: a monopole and a quadrupole. We explore systematically the effects of varying the quadrupole while keeping fixed the monopole and discuss the consequences for the interpretation of longitude-velocity diagrams in the Milky Way. We find that the gas flow can constrain the quadrupole of the potential and the characteristics of the bar that generates it. The exponential scale length of the bar must be at least $1.5\rm\, kpc$. The strength of the bar is also constrained. Our global interpretation favours a pattern speed of $\Omega=40\,\rm km s^{-1} {kpc}^{-1}$. We find that for most observational features, there exist a value of the parameters that matches each individual feature well, but is difficult to reproduce all the important features at once. Due to the intractably high number of parameters involved in the general problem, quantitative fitting methods that can run automatic searches in parameter space are necessary.

Our Galaxy's bar makes the Galaxy's potential distinctly non-axisymmetric. All orbits are affected by non-axisymmetry, and significant numbers are qualitatively changed by being trapped at a resonance with the bar. Orbital tori are used to compute these effects. Thick-disc orbits are no less likely to be trapped by corotation or a Lindblad resonance than thin-disc orbits. Perturbation theory is used to create non-axisymmetric orbital tori from standard axisymmetric tori, and both trapped and untrapped orbits are recovered to surprising accuracy. Code is added to the TorusModeller library that makes it as easy to manipulate non-axisymmetric tori as axisymmetric ones. The augmented TorusModeller is used to compute the velocity structure of the solar neighbourhood for bars of different pattern speeds and a simple action-based distribution function. The technique developed here can be applied to any non-axisymmetric potential that is stationary in a rotating from - hence also to classical spiral structure.

The present-day response of a Galactic disc stellar population to a non-axisymmetric perturbation of the potential has previously been computed through perturbation theory within the phase-space coordinates of the unperturbed axisymmetric system. Such an Eulerian linearized treatment however leads to singularities at resonances, which prevent quantitative comparisons with data. Here, we manage to capture the behaviour of the distribution function (DF) at a resonance in a Lagrangian approach, by averaging the Hamiltonian over fast angle variables and re-expressing the DF in terms of a new set of canonical actions and angles variables valid in the resonant region. We then follow the prescription of Binney (2016), assigning to the resonant DF the time average along the orbits of the axisymmetric DF expressed in the new set of actions and angles. This boils down to phase-mixing the DF in terms of the new angles, such that the DF for trapped orbits only depends on the new set of actions. This opens the way to quantitatively fitting the effects of the bar and spirals to Gaia data in terms of distribution functions in action space.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

We perform controlled N-body simulations of disc galaxies growing within live dark matter (DM) haloes to present-day galaxies that contain both thin and thick discs. We consider two types of models: a) thick-disc initial conditions to which stars on near-circular orbits are continuously added over ~10 Gyr and b) models in which the birth velocity dispersion of stars decreases continuously over the same time-scale. We show that both schemes produce double-exponential vertical profiles similar to that of the Milky Way (MW). We indicate how the spatial age structure of galaxies can be used to discriminate between scenarios. We show that the presence of a thick disc significantly alters and delays bar formation and thus makes possible models with a realistic bar and a high baryon-to-DM mass ratio in the central regions, as required by microlensing constraints. We examine how the radial mass distribution in stars and DM is affected by disc growth and non-axisymmetries. We discuss how bar buckling shapes the vertical age distribution of thin- and thick-disc stars in the bar region. The extent to which the combination of observationally motivated inside-out growth histories and cosmologically motivated dark halo properties leads to the spontaneous formation of non-axisymmetries that steer the models towards present-day MW-like galaxies is noteworthy.

We analyse disc heating and radial migration in N-body models of growing disc galaxies with thick and thin discs. Similar to thin-disc-only models, galaxies with appropriate non-axisymmetric structures reproduce observational constraints on radial disc heating in and migration to the Solar Neighbourhood (Snhd). The presence of thick discs can suppress non-axisymmetries and thus higher baryonic-to-dark matter fractions are required than in models that only have a thin disc. Models that are baryon-dominated to roughly the Solar radius R_0 are favoured, in agreement with data for the Milky Way. For inside-out growing discs, today's thick-disc stars at R_0 are dominated by outwards migrators. Whether outwards migrators are vertically hotter than non-migrators depends on the radial gradient of the thick disc vertical velocity dispersion. There is an effective upper boundary in angular momentum that thick disc stars born in the centre of a galaxy can reach by migration, which explains the fading of the high-alpha sequence outside R_0. Our models compare well to Snhd kinematics from RAVE-TGAS. For such comparisons it is important to take into account the azimuthal variation of kinematics at R ~ R_0 and biases from survey selection functions. The vertical heating of thin disc stars by giant molecular clouds is only mildly affected by the presence of thick discs. Our models predict higher vertical velocity dispersions for the oldest stars than found in the Snhd age-velocity dispersion relation, possibly because of measurement uncertainties or an underestimation of the number of old cold stars in our models.

Angle/action based distribution function (DF) models can be optimised based on how well they reproduce observations thus revealing the current matter distribution in the Milky Way. Gaia data combined with data from other surveys, e.g. the RAVE/TGAS sample, and its full selection function will greatly improve their accuracy.

Galaxy models are fundamental to exploiting surveys of our Galaxy. There is now a significant body of work on axisymmetric models. A model can be defined by giving the DF of each major class of stars and of dark matter. Then the self-consistent gravitational potential is determined. Other modelling techniques are briefly considered before an overview of some early work on non-axisymmetric models.

We characterize the selection function of RAVE using 2MASS as our underlying population, which we assume represents all stars which could have potentially been observed. We evaluate the completeness fraction as a function of position, magnitude, and color in two ways: first, on a field-by-field basis, and second, in equal-size areas on the sky. Then, we consider the effect of the RAVE stellar parameter pipeline on the final resulting catalogue, which in principle limits the parameter space over which our selection function is valid. Our final selection function is the product of the completeness fraction and the selection function of the pipeline. We then test if the application of the selection function introduces biases in the derived parameters. To do this, we compare a parent mock catalogue generated using Galaxia with a mock-RAVE catalogue where the selection function of RAVE has been applied. We conclude that for stars brighter than I = 12, between $4000 \rm K < T_{\rm eff} < 8000 \rm K$ and $0.5 < \rm{log}\,g < 5.0$, RAVE is kinematically and chemically unbiased with respect to expectations from Galaxia.

We construct a model for the Galactic globular cluster system based on a realistic gravitational potential and a distribution function (DF) analytic in the action integrals. The DF comprises disc and halo components whose functional forms resemble those recently used to describe the stellar discs and stellar halo. We determine the posterior distribution of our model parameters using a Bayesian approach. This gives us an understanding of how well the globular cluster data constrain our model. The favoured parameter values of the disc and halo DFs are similar to values previously obtained from fits to the stellar disc and halo, although the cluster halo system shows clearer rotation than does the stellar halo. Our model reproduces the generic features of the globular cluster system, namely the density profile, the mean rotation velocity. The fraction of disc clusters coincides with the observed fraction of metal-rich clusters. However, the data indicate either incompatibility between catalogued cluster distances and current estimates of distance to the Galactic Centre, or failure to identify clusters behind the bulge. As the data for our Galaxy's components increase in volume and precision over the next few years, it will be rewarding to revisit the present analysis.

Massive and diffuse haloes of hot gas (coronae) are important intermediaries between cosmology and galaxy evolution, storing mass and angular momentum acquired from the cosmic web until eventual accretion on to star-forming discs. We introduce a method to reconstruct the rotation of a galactic corona, based on its angular momentum distribution (AMD). This allows us to investigate in what conditions the angular momentum acquired from tidal torques can be transferred to star forming discs and explain observed galaxy-scale processes, such as inside-out growth and the build-up of abundance gradients. We find that a simple model of an isothermal corona with a temperature slightly smaller than virial and a cosmologically motivated AMD is in good agreement with galaxy evolution requirements, supporting hot-mode accretion as a viable driver for the evolution of spiral galaxies in a cosmological context. We predict moderately sub-centrifugal rotation close to the disc and slow rotation close to the virial radius. Motivated by the observation that the Milky Way has a relatively hot corona (T ~ 2 x 10^6 K), we also explore models with a temperature larger than virial. To be able to drive inside-out growth, these models must be significantly affected by feedback, either mechanical (ejection of low angular momentum material) or thermal (heating of the central regions). However, the agreement with galaxy evolution constraints becomes, in these cases, only marginal, suggesting that our first and simpler model may apply to a larger fraction of galaxy evolution history.

Data Release 5 (DR5) of the Radial Velocity Experiment (RAVE) is the fifth data release from a magnitude-limited (9< I < 12) survey of stars randomly selected in the southern hemisphere. The RAVE medium-resolution spectra ($R\sim7500$) covering the Ca-triplet region (8410-8795\AA) span the complete time frame from the start of RAVE observations in 2003 to their completion in 2013. Radial velocities from 520,781 spectra of 457,588 unique stars are presented, of which 255,922 stellar observations have parallaxes and proper motions from the Tycho-Gaia astrometric solution (TGAS) in Gaia DR1. For our main DR5 catalog, stellar parameters (effective temperature, surface gravity, and overall metallicity) are computed using the RAVE DR4 stellar pipeline, but calibrated using recent K2 Campaign 1 seismic gravities and Gaia benchmark stars, as well as results obtained from high-resolution studies. Also included are temperatures from the Infrared Flux Method, and we provide a catalogue of red giant stars in the dereddened color $(J-Ks)_0$ interval (0.50,0.85) for which the gravities were calibrated based only on seismology. Further data products for sub-samples of the RAVE stars include individual abundances for Mg, Al, Si, Ca, Ti, Fe, and Ni, and distances found using isochrones. Each RAVE spectrum is complemented by an error spectrum, which has been used to determine uncertainties on the parameters. The data can be accessed via the RAVE Web site or the Vizier database.

We construct a new family of models of our Galaxy in which dark matter and disc stars are both represented by distribution functions that are analytic functions of the action integrals of motion. The potential that is self-consistently generated by the dark matter, stars and gas is determined, and parameters in the distribution functions are adjusted until the model is compatible with observational constraints on the circular-speed curve, the vertical density profile of the stellar disc near the Sun, the kinematics of nearly 200 000 giant stars within 2 kpc of the Sun, and estimates of the optical depth to microlensing of bulge stars. We find that the data require a dark halo in which the phase-space density is approximately constant for actions |J| \lesssim 140 kpc km ^-1. In real space these haloes have core radii ~ 2 kpc.

The distribution of Milky Way halo blue horizontal-branch (BHB) stars is examined using action-based extended distribution functions (EDFs) that describe the locations of stars in phase space, metallicity, and age. The parameters of the EDFs are fitted using stars observed in the Sloan Extension for Galactic Understanding and Exploration-II (SEGUE-II) survey that trace the phase-space kinematics and chemistry out to ~70 kpc. A maximum a posteriori probability (MAP) estimate method and a Markov Chain Monte Carlo method are applied, taking into account the selection function in positions, distance, and metallicity for the survey. The best-fit EDF declines with actions less steeply at actions characteristic of the inner halo than at the larger actions characteristic of the outer halo, and older ages are found at smaller actions than at larger actions. In real space, the radial density profile steepens smoothly from -2 at ~2 kpc to -4 in the outer halo, with an axis ratio ~0.7 throughout. There is no indication for rotation in the BHBs, although this is highly uncertain. A moderate level of radial anisotropy is detected, with $\beta_s$ varying from isotropic to between ~0.1 and ~0.3 in the outer halo depending on latitude. The BHB data are consistent with an age gradient of -0.03 Gyr kpc$^{-1}$, with some uncertainty in the distribution of the larger ages. These results are consistent with a scenario in which older, larger systems contribute to the inner halo, whilst the outer halo is primarily comprised of younger, smaller systems.

Galaxy modelling is greatly simplified by assuming the existence of a global system of angle-action coordinates. Unfortunately, global angle-action coordinates do not exist because some orbits become trapped by resonances, especially where the radial and vertical frequencies coincide. We show that in a realistic Galactic potential such trapping occurs only on thick-disc and halo orbits (speed relative to the guiding centre >~80 km/s). We explain how the Torus Mapper code (TM) behaves in regions of phase space in which orbits are resonantly trapped, and we extend TM so trapped orbits can be manipulated as easily as untrapped ones. The impact that the resonance has on the structure of velocity space depends on the weights assigned to trapped orbits. The impact is everywhere small if each trapped orbit is assigned the phase space density equal to the time average along the orbit of the DF for untrapped orbits. The impact could be significant with a different assignment of weights to trapped orbits.

We analyse the heating of stellar discs by non axisymmetric structures and giant molecular clouds (GMCs) in N-body simulations of growing disc galaxies. The analysis resolves long-standing discrepancies between models and data by demonstrating the importance of distinguishing between measured age-velocity dispersion relations (AVRs) and the heating histories of the stars that make up the AVR. We fit both AVRs and heating histories with formulae proportional to t^beta and determine the exponents beta_R and beta_z derived from in-plane and vertical AVRs and ~beta_R and ~beta_z from heating histories. Values of beta_z are in almost all simulations larger than values of ~beta_z, whereas values of beta_R are similar to or mildly larger than values of ~beta_R. Moreover, values of beta_z (~beta_z) are generally larger than values of beta_R (~beta_R). The dominant cause of these relations is the decline over the life of the disc in importance of GMCs as heating agents relative to spiral structure and the bar. We examine how age errors and biases in solar neighbourhood surveys influence the measured AVR: they tend to decrease beta values by smearing out ages and thus measured dispersions. We compare AVRs and velocity ellipsoid shapes sigma_z/sigma_R from simulations to Solar neighbourhood data. We conclude that for the expected disc mass and dark halo structure, combined GMC and spiral/bar heating can explain the AVR of the Galactic thin disc. Strong departures of the disc mass or the dark halo structure from expectation spoil fits to the data.

We present evidence from the RAdial Velocity Experiment (RAVE) survey of chemically separated, kinematically distinct disc components in the solar neighbourhood. We apply probabilistic chemical selection criteria to separate our sample into $\alpha$-low (`thin disc') and $\alpha$-high (`thick disc') sequences. Using newly derived distances, which will be utilized in the upcoming RAVE DR5, we explore the kinematic trends as a function of metallicity for each of the disc components. For our $\alpha$-low disc, we find a negative trend in the mean rotational velocity ($V_{\mathrm{\phi}}$) as a function of iron abundance ([Fe/H]). We measure a positive gradient $\partial V_{\mathrm{\phi}}$/$\partial$[Fe/H] for the $\alpha$-high disc, consistent with results from high-resolution surveys. We also find differences between the $\alpha$-low and $\alpha$-high discs in all three components of velocity dispersion. We discuss the implications of an $\alpha$-low, metal-rich population originating from the inner Galaxy, where the orbits of these stars have been significantly altered by radial mixing mechanisms in order to bring them into the solar neighbourhood. The probabilistic separation we propose can be extended to other data sets for which the accuracy in [$\alpha$/Fe] is not sufficient to disentangle the chemical disc components a priori. For such datasets which will also have significant overlap with Gaia DR1, we can therefore make full use of the improved parallax and proper motion data as it becomes available to investigate kinematic trends in these chemical disc components.

We perform a series of controlled N-body simulations of growing disc galaxies within non-growing, live dark matter haloes of varying mass and concentration. Our initial conditions include either a low-mass disc or a compact bulge. New stellar particles are continuously added on near-circular orbits to the existing disc, so spiral structure is continuously excited. To study the effect of combined spiral and giant molecular cloud (GMC) heating on the discs we introduce massive, short-lived particles that sample a GMC mass function. An isothermal gas component is introduced for a subset of the models. We perform a resolution study and vary parameters governing the GMC population, the histories of star formation and radial scale growth. Models with GMCs and standard values for the disc mass and halo density provide the right level of self-gravity to explain the age velocity dispersion relation of the Solar neighbourhood (Snhd). GMC heating generates remarkably exponential vertical profiles with scaleheights that are radially constant and agree with observations of galactic thin discs. GMCs are also capable of significantly delaying bar formation. The amount of spiral induced radial migration agrees with what is required for the metallicity distribution of the Snhd. However, in our standard models the outward migrating populations are not hot enough vertically to create thick discs. Thick discs can form in models with high baryon fractions, but the corresponding bars are too long, the young stellar populations too hot and the discs flare considerably.

We fit an Extended Distribution Function (EDF) to K giants in the Sloan Extension for Galactic Understanding and Exploration (SEGUE) survey. These stars are detected to radii ~80 kpc and span a wide range in [Fe/H]. Our EDF, which depends on [Fe/H] in addition to actions, encodes the entanglement of metallicity with dynamics within the Galaxy's stellar halo. Our maximum-likelihood fit of the EDF to the data allows us to model the survey's selection function. The density profile of the K giants steepens with radius from a slope ~-2 to ~-4 at large radii. The halo's axis ratio increases with radius from 0.7 to almost unity. The metal-rich stars are more tightly confined in action space than the metal-poor stars and form a more flattened structure. A weak metallicity gradient ~-0.001 dex/kpc, a small gradient in the dispersion in [Fe/H] of ~0.001 dex/kpc, and a higher degree of radial anistropy in metal-richer stars result. Lognormal components with peaks at ~-1.5 and ~-2.3 are required to capture the overall metallicity distribution, suggestive of the existence of two populations of K giants. The spherical anisotropy parameter varies between 0.3 in the inner halo to isotropic in the outer halo. If the Sagittarius stream is included, a very similar model is found but with a stronger degree of radial anisotropy throughout.

We review the available methods for estimating actions, angles and frequencies of orbits in both axisymmetric and triaxial potentials. The methods are separated into two classes. Unless an orbit has been trapped by a resonance, convergent, or iterative, methods are able to recover the actions to arbitrarily high accuracy given sufficient computing time. Faster non-convergent methods rely on the potential being sufficiently close to a separable potential and the accuracy of the action estimate cannot be improved through further computation. We critically compare the accuracy of the methods and the required computation time for a range of orbits in an axisymmetric multi-component Galactic potential. We introduce a new method for estimating actions that builds on the adiabatic approximation of Sch\"onich & Binney (2012) and discuss the accuracy required for the actions, angles and frequencies using suitable distribution functions for the thin and thick discs, the stellar halo and a star stream. We conclude that for studies of the disc and smooth halo component of the Milky Way the most suitable compromise between speed and accuracy is the St\"ackel Fudge, whilst when studying streams the non-convergent methods do not offer sufficient accuracy and the most suitable method is computing the actions from an orbit integration via a generating function. All the software used in this study can be downloaded from https://github.com/jls713/tact.

Multi-object spectrographs have opened a new window on the analyses of the chemo-dynamical properties of old Milky Way stars. These analyses allow us to trace back the internal mechanisms and the external factors that have influenced the evolution of our Galaxy, and therefore understand fundamental aspects of galaxy evolution in general. Here, we present recent results from the RAdial Velocity Experiment (RAVE) and the Gaia-ESO survey. These surveys explore the Milky Way properties in different ways, in terms of sample size and selection, magnitude range, and spectral resolution. We focus here on (i) the first direct detection of evidence for radial migration within the thin disc, providing insight into the history of spiral structure of the Milky Way, and (ii) the chemo-dynamical characterisation of the metal-weak thick and thin discs, for which chemo-dynamical models still have difficulties in reproducing.

Chemodynamical models of our Galaxy that have analytic Extended Distribution Functions (EDFs) are likely to play a key role in extracting science from surveys in the era of Gaia.

We present a freely downloadable software package for modelling the dynamics of galaxies, which we call the Torus Mapper (TM). The package is based around `torus mapping', which is a non-perturbative technique for creating orbital tori for specified values of the action integrals. Given an orbital torus and a star's position at a reference time, one can compute its position at any other time, no matter how remote. One can also compute the velocities with which the star will pass through any given point and the contribution it will make to the time-averaged density there. A system of angle-action coordinates for the given potential can be created by foliating phase space with orbital tori. Such a foliation is facilitated by the ability of TM to create tori by interpolating on a grid of tori. We summarise the advantages of using TM rather than a standard time-stepper to create orbits, and give segments of code that illustrate applications of TM in several contexts, including setting up initial conditions for an N-body simulation. We examine the precision of the orbital tori created by TM and the behaviour of the code when orbits become trapped by a resonance.

Starting from the hypothesis that the Galaxy's dark halo responded adiabatically to the infall of baryons, we have constructed a self-consistent dynamical model of the Galaxy that satisfies a large number of observations, including measurements of gas terminal velocities and masers, the kinematics of a 180,000 giant stars from the RAVE survey, and star count data from the SDSS. The stellar disc and the dark halo are both specified by distribution functions (DFs) of the action integrals. The model is obtained by extending the work of Piffl Penoyre & Binney (2015} from the construction of a single model to a systematic search of model space. Whereas the model of Piffl et al violated constraints on the terminal-velocity curve, our model respects these constraints by adopting a long scale length R_d=3.66 kpc for the thin and thick discs. The model is, however, inconsistent with the measured optical depth for microlensing of bulge stars because it attributes too large a fraction of the density at R <~ 3 kpc to dark matter rather than stars. Moreover, it now seems likely that the thick disc's scale-length is significantly shorter than the model implies. Shortening this scale-length would cause the constraints from the rotation curve to be violated anew. We conclude that we can now rule out adiabatic compression of our Galaxy's dark halo.

We extend models of our Galaxy based on distribution functions (DFs) that are analytic functions of the action integrals to extended distribution functions (EDFs), which have an analytic dependence on metallicity as well. We use a simple, but physically-motivated, functional forms for the metallicity of the interstellar medium as a function of radius and time and for the star-formation rate, and a model for the diffusion of stars through phase space to suggest the required functional form of an EDF. We introduce a simple prescription for radial migration that preserves the overall profile of the disc while allowing individual stars to migrate throughout the disc. Our models explicitly consider the thin and thick discs as two distinct components separated in age. We show how an EDF can be used to incorporate realistic selection functions in models, and to construct mock catalogues of observed samples. We show that the selection function of the Geneva-Copenhagen Survey (GCS) biases in favour of young stars, which have atypically small random velocities. With the selection function taken into account our models produce good fits of the GCS data in chemo-dynamical space and the Gilmore and Reid (1983) density data. From our EDF, we predict the structure of the SEGUE G-dwarf sample. The kinematics are successfully predicted. The predicted metallicity distribution has too few stars with [Fe/H]$\simeq-0.5$ dex and too many metal-rich stars. A significant problem may be the lack of any chemical-kinematic correlations in our thick disc. We argue that EDFs will prove essential tools for the analysis of both observational data and sophisticated models of Galaxy formation and evolution.

The Gaia-ESO Survey (GES) is a large public spectroscopic survey at the European Southern Observatory Very Large Telescope. A key aim is to provide precise radial velocities (RVs) and projected equatorial velocities (v sin i) for representative samples of Galactic stars, that will complement information obtained by the Gaia astrometry satellite. We present an analysis to empirically quantify the size and distribution of uncertainties in RV and v sin i using spectra from repeated exposures of the same stars. We show that the uncertainties vary as simple scaling functions of signal-to-noise ratio (S/N) and v sin i, that the uncertainties become larger with increasing photospheric temperature, but that the dependence on stellar gravity, metallicity and age is weak. The underlying uncertainty distributions have extended tails that are better represented by Student's t-distributions than by normal distributions. Parametrised results are provided, that enable estimates of the RV precision for almost all GES measurements, and estimates of the v sin i precision for stars in young clusters, as a function of S/N, v sin i and stellar temperature. The precision of individual high S/N GES RV measurements is 0.22-0.26 km/s, dependent on instrumental configuration.

Spiral arms that emerge from the ends of a galactic bar are important in interpreting observations of our and external galaxies. It is therefore important to understand the physical mechanism that causes them. We find that these spiral arms can be understood as kinematic density waves generated by librations around underlying ballistic closed orbits. This is even true in the case of a strong bar, provided the librations are around the appropriate closed orbits and not around the circular orbits that form the basis of the epicycle approximation. An important consequence is that it is a potential's orbital structure that determines whether a bar should be classified as weak or strong, and not crude estimates of the potential's deviation from axisymmetry.

We present the identification of potential members of nearby Galactic globular clusters using radial velocities from the RAdial Velocity Experiment Data Release 4 (RAVE-DR4) survey database. Our identifications are based on three globular clusters -- NGC 3201, NGC 5139 ($\omega$ Cen) and NGC 362 -- all of which are shown to have |RV|>100 km/s. The identification of globular cluster stars in RAVE DR4 data offers a unique opportunity to test the precision and accuracy of the stellar parameters determined with the currently available Stellar Parameter Pipelines (SPPs) used in the survey, as globular clusters are ideal testbeds for the validation of stellar atmospheric parameters, abundances, distances and ages. For both NGC 3201 and $\omega$ Cen, there is compelling evidence for numerous members (> 10) in the RAVE database; in the case of NGC 362 the evidence is more ambiguous, and there may be significant foreground and/or background contamination in our kinematically-selected sample. A comparison of the RAVE-derived stellar parameters and abundances with published values for each cluster and with BASTI isochrones for ages and metallicities from the literature reveals overall good agreement, with the exception of the apparent underestimation of surface gravities for giants, in particular for the most metal-poor stars. Moreover, if the selected members are part of the main body of each cluster our results would also suggest that the distances from Binney et al. 2013, where only isochrones more metal-rich than -0.9 dex were used, are typically underestimated by ~ 40% with respect to the published distances for the clusters, while the distances from Zwitter et al. 2010 show stars ranging from 1 to ~ 6.5 kpc -- with indications of a trend toward higher distances at lower metallicities -- for the three clusters analysed in this study.

We present a new method to construct fully self-consistent equilibrium models of multi-component disc galaxies similar to the Milky Way. We define distribution functions for the stellar disc and dark halo that depend on phase space position only through action coordinates. We then use an iterative approach to find the corresponding gravitational potential. We study the adiabatic response of the initially spherical dark halo to the introduction of the baryonic component and find that the halo flattens in its inner regions with final minor-major axis ratios $q$ = 0.75 - 0.95. The extent of the flattening depends on the velocity structure of the halo particles with radially biased models exhibiting a stronger response. In this latter case, which is according to cosmological simulations the most likely one, the new density structure resembles a "dark disc" superimposed on a spherical halo. We discuss the implications of these results for our recent estimate of the local dark matter density. The velocity distribution of the dark-matter particles near the Sun is very non-Gaussian. All three principal velocity dispersions are boosted as the halo contracts, and at low velocities a plateau develops in the distribution of $v_z$. For models similar to a state-of-the-art Galaxy model we find velocity dispersions around 155 km s$^{-1}$ for $v_z$ and the tangential velocity, $v_\varphi$, and 140 - 175 km s$^{-1}$ for the in-plane radial velocity, $v_R$, depending on the anisotropy of the model.

Fluctuations in a stellar system's gravitational field cause the orbits of stars to evolve. The resulting evolution of the system can be computed with the orbit-averaged Fokker-Planck equation once the diffusion tensor is known. We present the formalism that enables one to compute the diffusion tensor from a given source of noise in the gravitational field when the system's dynamical response to that noise is included. In the case of a cool stellar disc we are able to reduce the computation of the diffusion tensor to a one-dimensional integral. We implement this formula for a tapered Mestel disc that is exposed to shot noise and find that we are able to explain analytically the principal features of a numerical simulation of such a disc. In particular the formation of narrow ridges of enhanced density in action space is recovered. As the disc's value of Toomre's $Q$ is reduced and the disc becomes more responsive, there is a transition from a regime of heating in the inner regions of the disc through the inner Lindblad resonance to one of radial migration of near-circular orbits via the corotation resonance in the intermediate regions of the disc. The formalism developed here provides the ideal framework in which to study the long-term evolution of all kinds of stellar discs.

We use a Cartesian grid to simulate the flow of gas in a barred Galactic potential and investigate the effects of varying the sound speed in the gas and the resolution of the grid. For all sound speeds and resolutions, streamlines closely follow closed orbits at large and small radii. At intermediate radii shocks arise and the streamlines shift between two families of closed orbits. The point at which the shocks appear and the streamlines shift between orbit families depends strongly on sound speed and resolution. For sufficiently large values of these two parameters, the transfer happens at the cusped orbit as hypothesised by Binney et al. over two decades ago. For sufficiently high resolutions the flow downstream of the shocks becomes unsteady. If this unsteadiness is physical, as appears to be the case, it provides a promising explanation for the asymmetry in the observed distribution of CO.

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

Using the RAdial Velocity Experiment fourth data release (RAVE DR4), and a new metallicity calibration that will be also taken into account in the future RAVE DR5, we investigate the existence and the properties of super-solar metallicity stars ([M/H] > +0.1 dex) in the sample, and in particular in the Solar neighbourhood. We find that RAVE is rich in super-solar metallicity stars, and that the local metallicity distribution function declines remarkably slowly up to +0.4 dex. Our results show that the kinematics and height distributions of the super-solar metallicity stars are identical to those of the [M/H] < 0 thin-disc giants that we presume were locally manufactured. The eccentricities of the super-solar metallicity stars indicate that half of them are on a roughly circular orbit (e < 0.15), so under the assumption that the metallicity of the interstellar medium at a given radius never decreases with time, they must have increased their angular momenta by scattering at corotation resonances of spiral arms from regions far inside the Solar annulus. The likelihood that a star will migrate radially does not seem to decrease significantly with increasing amplitude of vertical oscillations within range of oscillation amplitudes encountered in the disc.

We present an approach to approximating rapidly the actions in a general triaxial potential. The method is an extension of the axisymmetric approach presented by Binney (2012), and operates by assuming that the true potential is locally sufficiently close to some St\"ackel potential. The choice of St\"ackel potential and associated ellipsoidal coordinates is tailored to each individual input phase-space point. We investigate the accuracy of the method when computing actions in a triaxial Navarro-Frenk-White potential. The speed of the algorithm comes at the expense of large errors in the actions, particularly for the box orbits. However, we show that the method can be used to recover the observables of triaxial systems from given distribution functions to sufficient accuracy for the Jeans equations to be satisfied. Consequently, such models could be used to build models of external galaxies as well as triaxial components of our own Galaxy. When more accurate actions are required, this procedure can be combined with torus mapping to produce a fast convergent scheme for action estimation.

We present an approach to the design of distribution functions that depend on the phase-space coordinates through the action integrals. The approach makes it easy to construct a dynamical model of a given stellar component. We illustrate the approach by deriving distribution functions that self-consistently generate several popular stellar systems, including the Hernquist, Jaffe, and Navarro, Frenk and White models. We focus on non-rotating spherical systems, but extension to flattened and rotating systems is trivial. Our distribution functions are easily added to each other and to previously published distribution functions for discs to create self-consistent multi-component galaxies. The models this approach makes possible should prove valuable both for the interpretation of observational data and for exploring the non-equilibrium dynamics of galaxies via N-body simulation.

Henon sought the most general spherical potential in which the radial periods of orbits depended only on energy. He named this potential the isochrone, and discovered that it provided a good representation of data for globular clusters. He sought an explanation in terms of resonant relaxation. The role that resonant relaxation might play in globular clusters is still an open question, but the isochrone potential is guaranteed a role in dynamical astronomy because it is the most general potential in which closed-form expressions for angle-action coordinates are available. I explain how this property makes the isochrone invaluable for the powerful technique of torus mapping. I also describe flattened isochrone models, which enable us to explore a powerful general method of generating self-consistent stellar systems.

Simulations of the clustering of cold dark matter yield dark-matter halos that have central density cusps, but observations of totally dark-matter dominated dwarf spheroidal galaxies imply that they do not have cuspy central density profiles. We use analytic calculations and numerical modelling to argue that whenever stars form, central density cusps are likely to be erased. Gas that accumulates in the potential well of an initially cuspy dark-matter halo settles into a disc. Eventually the surface density of the gas exceeds the threshold for fragmentation into self-gravitating clouds. The clouds are massive enough to transfer energy to the dark-matter particles via dynamical friction on a short time-scale. The halo's central cusp is heated to form a core with central logarithmic density slope gamma=0 before stellar feedback makes its impact. Since star formation is an inefficient process, the clouds are disrupted by feedback when only a small fraction of their mass has been converted to stars, and the dark matter dominates the final mass distribution.

We use the kinematics of $\sim200\,000$ giant stars that lie within $\sim 1.5$ kpc of the plane to measure the vertical profile of mass density near the Sun. We find that the dark mass contained within the isodensity surface of the dark halo that passes through the Sun ($(6\pm0.9)\times10^{10}\,\mathrm{M_\odot}$), and the surface density within $0.9$ kpc of the plane ($(69\pm10)\,\mathrm{M_\odot\,pc^{-2}}$) are almost independent of the (oblate) halo's axis ratio $q$. If the halo is spherical, 46 per cent of the radial force on the Sun is provided by baryons, and only 4.3 per cent of the Galaxy's mass is baryonic. If the halo is flattened, the baryons contribute even less strongly to the local radial force and to the Galaxy's mass. The dark-matter density at the location of the Sun is $0.0126\,q^{-0.89}\,\mathrm{M_\odot\,pc^{-3}}=0.48\,q^{-0.89}\,\mathrm{GeV\,cm^{-3}}$. When combined with other literature results we find hints for a mildly oblate dark halo with $q \simeq 0.8$. Our value for the dark mass within the solar radius is larger than that predicted by cosmological dark-matter-only simulations but in good agreement with simulations once the effects of baryonic infall are taken into account. Our mass models consist of three double-exponential discs, an oblate bulge and a Navarro-Frenk-White dark-matter halo, and we model the dynamics of the RAVE stars in the corresponding gravitational fields by finding distribution functions $f(\mathbf{J})$ that depend on three action integrals. Statistical errors are completely swamped by systematic uncertainties, the most important of which are the distance to the stars in the photometric and spectroscopic samples and the solar distance to the Galactic centre. Systematics other than the flattening of the dark halo yield overall uncertainties $\sim 15$ per cent.

The Gaia-ESO Survey is obtaining high-quality spectroscopic data for about 10^5 stars using FLAMES at the VLT. UVES high-resolution spectra are being collected for about 5000 FGK-type stars. These UVES spectra are analyzed in parallel by several state-of-the-art methodologies. Our aim is to present how these analyses were implemented, to discuss their results, and to describe how a final recommended parameter scale is defined. We also discuss the precision (method-to-method dispersion) and accuracy (biases with respect to the reference values) of the final parameters. These results are part of the Gaia-ESO 2nd internal release and will be part of its 1st public release of advanced data products. The final parameter scale is tied to the one defined by the Gaia benchmark stars, a set of stars with fundamental atmospheric parameters. A set of open and globular clusters is used to evaluate the physical soundness of the results. Each methodology is judged against the benchmark stars to define weights in three different regions of the parameter space. The final recommended results are the weighted-medians of those from the individual methods. The recommended results successfully reproduce the benchmark stars atmospheric parameters and the expected Teff-log g relation of the calibrating clusters. Atmospheric parameters and abundances have been determined for 1301 FGK-type stars observed with UVES. The median of the method-to-method dispersion of the atmospheric parameters is 55 K for Teff, 0.13 dex for log g, and 0.07 dex for [Fe/H]. Systematic biases are estimated to be between 50-100 K for Teff, 0.10-0.25 dex for log g, and 0.05-0.10 dex for [Fe/H]. Abundances for 24 elements were derived: C, N, O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Y, Zr, Mo, Ba, Nd, and Eu. The typical method-to-method dispersion of the abundances varies between 0.10 and 0.20 dex.

The diffuse interstellar bands (DIBs) are absorption lines observed in visual and near infrared spectra of stars. Understanding their origin in the interstellar medium (ISM) is one of the oldest problems in astronomical spectroscopy, as DIBs ahave been known since 1922. In a completely new approach to understanding DIBs, we combined information from nearly 500,000 stellar spectra obtained by the massive spectroscopic survey RAVE (Radial Velocity Experiment) to produce the first pseudo three-dimensional map of the strength of the DIB at 8620~{\AA}ngstroms covering the nearest 3 kiloparsecs from the Sun, and show that it follows our independently constructed spatial distribution of extinction by interstellar dust along the Galactic plane. Despite having a similar distribution in the Galactic plane, the DIB 8620 carrier has a significantly larger vertical scale height than the dust. Even if one DIB may not represent the general DIB population, our observations outline the future direction of DIB research.

We present a method for extracting actions, angles and frequencies from an orbit's time series. The method recovers the generating function that maps an analytic phase-space torus to the torus to which the orbit is confined by simultaneously solving the constraints provided by each time step. We test the method by recovering the actions and frequencies of a triaxial St\"ackel potential, and use it to investigate the structure of orbits in a triaxial potential that has been fitted to our Galaxy's Sagittarius stream. The method promises to be useful for analysing N-body simulations. It also takes a step towards constructing distribution functions for the triaxial components of our Galaxy, such as the bar and dark halo.

We investigate the kinematic parameters of the Milky Way disc using the RAVE and GCS stellar surveys. We do this by fitting a kinematic model to the data taking the selection function of the data into account. For stars in the GCS we use all phase-space coordinates, but for RAVE stars we use only $(l,b,v_{\rm los})$. Using MCMC technique, we investigate the full posterior distributions of the parameters given the data. We investigate the `age-velocity dispersion' relation for the three kinematic components ($\sigma_R,\sigma_{\phi},\sigma_z$), the radial dependence of the velocity dispersions, the Solar peculiar motion ($U_{\odot},V_{\odot}, W_{\odot} $), the circular speed $\Theta_0$ at the Sun and the fall of mean azimuthal motion with height above the mid-plane. We confirm that the Besan\c{c}on-style Gaussian model accurately fits the GCS data, but fails to match the details of the more spatially extended RAVE survey. In particular, the Shu distribution function (DF) handles non-circular orbits more accurately and provides a better fit to the kinematic data. The Gaussian distribution function not only fits the data poorly but systematically underestimates the fall of velocity dispersion with radius. We find that correlations exist between a number of parameters, which highlights the importance of doing joint fits. The large size of the RAVE survey, allows us to get precise values for most parameters. However, large systematic uncertainties remain, especially in $V_{\odot}$ and $\Theta_0$. We find that, for an extended sample of stars, $\Theta_0$ is underestimated by as much as $10\%$ if the vertical dependence of the mean azimuthal motion is neglected. Using a simple model for vertical dependence of kinematics, we find that it is possible to match the Sgr A* proper motion without any need for $V_{\odot}$ being larger than that estimated locally by surveys like GCS.

The Gaia-ESO Survey is a large public spectroscopic survey that aims to derive radial velocities and fundamental parameters of about 10^5 Milky Way stars in the field and in clusters. Observations are carried out with the multi-object optical spectrograph FLAMES, using simultaneously the medium resolution (R~20,000) GIRAFFE spectrograph and the high resolution (R~47,000) UVES spectrograph. In this paper, we describe the methods and the software used for the data reduction, the derivation of the radial velocities, and the quality control of the FLAMES-UVES spectra. Data reduction has been performed using a workflow specifically developed for this project. This workflow runs the ESO public pipeline optimizing the data reduction for the Gaia-ESO Survey, performs automatically sky subtraction, barycentric correction and normalisation, and calculates radial velocities and a first guess of the rotational velocities. The quality control is performed using the output parameters from the ESO pipeline, by a visual inspection of the spectra and by the analysis of the signal-to-noise ratio of the spectra. Using the observations of the first 18 months, specifically targets observed multiple times at different epochs, stars observed with both GIRAFFE and UVES, and observations of radial velocity standards, we estimated the precision and the accuracy of the radial velocities. The statistical error on the radial velocities is sigma~0.4 km s^-1 and is mainly due to uncertainties in the zero point of the wavelength calibration. However, we found a systematic bias with respect to the GIRAFFE spectra (~0.9 km s^-1) and to the radial velocities of the standard stars (~0.5 kms^-1) retrieved from the literature. This bias will be corrected in the future data releases, when a common zero point for all the setups and instruments used for the survey will be established.

We present a family of self-consistent axisymmetric stellar systems that have analytic distribution functions (DFs) of the form f(J), so they depend on three integrals of motion and have triaxial velocity ellipsoids. The models, which are generalisations of Henon's isochrone sphere, have four dimensionless parameters, two determining the part of the DF that is even in L_z, and two determining the odd part of the DF (which determines the azimuthal velocity distribution). Outside their cores, the velocity ellipsoids of all models tend to point to the model's centre, and we argue that this behaviour is generic, so near the symmetry axis of a flattened model, the long axis of the velocity ellipsoid is naturally aligned with the symmetry axis and not perpendicular to it as in many published dynamical models of well-studied galaxies. By varying one of the DF's parameters, the intensity of rotation can be increased from zero up to a maximum value set by the requirement that the DF be non-negative. Since angle-action coordinates are easily computed for these models, they are ideally suited for perturbative treatments and stability analysis. They can also be used to choose initial conditions for an N-body model that starts in perfect equilibrium and to model observations of early-type galaxies. The modelling technique introduced here is readily extended to different radial density profiles, more complex kinematics, and multi-component systems. A number of important technical issues surrounding the determination of the models' observable properties are explained in two appendices.

We analyse the kinematics of ~400000 RAVE stars. We split the sample into hot and cold dwarfs, red-clump and non-clump giants. The kinematics of the clump giants are consistent with being identical with those of non-clump giants. We fit Gaussian velocity ellipsoids to the meridional-plane components of velocity of each star class and give formulae from which the shape and orientation of the velocity ellipsoid can be determined at any location. The data are consistent with the giants and the cool dwarfs sharing the same velocity ellipsoids; sigma_z rises from 21 kms in the plane to sim 55 kms at |z|=2 kpc, while sigma_r rises from 37 kms to 82 kms. At (R,z) the longest axis of one of these velocity ellipsoids is inclined to the Galactic plane by an angle ~0.8 arctan(z/R). We use a novel formula to obtain precise fits to the highly non-Gaussian distributions of v_phi components. We compare the observed velocity distributions with the predictions of a dynamical model fitted to the velocities of stars that lie within ~150 pc of the Sun and star counts towards the Galactic pole. The model accurately reproduces the non-Gaussian nature of the v_r and v_z distributions and provides excellent fits to the data for v_z at all locations. The model v_phi distributions for the cool dwarfs fit the data extremely well, while those for the hot dwarfs have displacements to low v_phi that grow with |z| from very small values near the plane. At |z|>0.5 kpc, the theoretical v_phi distributions for giants show a deficit of stars with large v_phi and the model v_r distributions are too narrow. Systematically over-estimating distances by 20 per cent introduces asymmetry into the model v_r and v_z distributions near the plane and but significantly improves the fits to the data at |z|>0.5 kpc. The quality of the fits lends credence to the assumed, disc-dominated, gravitational potential.

Probability density functions are determined from new stellar parameters for the distance moduli of stars for which the RAdial Velocity Experiment (RAVE) has obtained spectra with S/N>=10. Single-Gaussian fits to the pdf in distance modulus suffice for roughly half the stars, with most of the other half having satisfactory two-Gaussian representations. As expected, early-type stars rarely require more than one Gaussian. The expectation value of distance is larger than the distance implied by the expectation of distance modulus; the latter is itself larger than the distance implied by the expectation value of the parallax. Our parallaxes of Hipparcos stars agree well with the values measured by Hipparcos, so the expectation of parallax is the most reliable distance indicator. The latter are improved by taking extinction into account. The effective temperature absolute-magnitude diagram of our stars is significantly improved when these pdfs are used to make the diagram. We use the method of kinematic corrections devised by Schoenrich, Binney & Asplund to check for systematic errors for general stars and confirm that the most reliable distance indicator is the expectation of parallax. For cool dwarfs and low-gravity giants <pi> tends to be larger than the true distance by up to 30 percent. The most satisfactory distances are for dwarfs hotter than 5500 K. We compare our distances to stars in 13 open clusters with cluster distances from the literature and find excellent agreement for the dwarfs and indications that we are over-estimating distances to giants, especially in young clusters.

The velocity dispersions of stars near the Sun are known to increase with stellar age, but age can be difficult to determine so a proxy like the abundance of alpha elements (e.g., Mg) with respect to iron, [alpha/Fe], is used. Here we report an unexpected behavior found in the velocity dispersion of a sample of giant stars from the RAdial Velocity Experiment (RAVE) survey with high quality chemical and kinematical information, in that it decreases strongly for stars with [Mg/Fe] > 0.4 dex (i.e., those that formed in the first Gyr of the Galaxy's life). These findings can be explained by perturbations from massive mergers in the early Universe, which have affected more strongly the outer parts of the disc, and the subsequent radial migration of stars with cooler kinematics from the inner disc. Similar reversed trends in velocity dispersion are also found for different metallicity subpopulations. Our results suggest that the Milky Way disc merger history can be recovered by relating the observed chemo-kinematic relations to the properties of past merger events.

We construct new estimates on the Galactic escape speed at various Galactocentric radii using the latest data release of the Radial Velocity Experiment (RAVE DR4). Compared to previous studies we have a database larger by a factor of 10 as well as reliable distance estimates for almost all stars. Our analysis is based on the statistical analysis of a rigorously selected sample of 90 high-velocity halo stars from RAVE and a previously published data set. We calibrate and extensively test our method using a suite of cosmological simulations of the formation of Milky Way-sized galaxies. Our best estimate of the local Galactic escape speed, which we define as the minimum speed required to reach three virial radii $R_{340}$, is $533^{+54}_{-41}$ km/s (90% confidence) with an additional 5% systematic uncertainty, where $R_{340}$ is the Galactocentric radius encompassing a mean over-density of 340 times the critical density for closure in the Universe. From the escape speed we further derive estimates of the mass of the Galaxy using a simple mass model with two options for the mass profile of the dark matter halo: an unaltered and an adiabatically contracted Navarro, Frenk & White (NFW) sphere. If we fix the local circular velocity the latter profile yields a significantly higher mass than the un-contracted halo, but if we instead use the statistics on halo concentration parameters in large cosmological simulations as a constraint we find very similar masses for both models. Our best estimate for $M_{340}$, the mass interior to $R_{340}$ (dark matter and baryons), is $1.3^{+0.4}_{-0.3} \times 10^{12}$ M$_\odot$ (corresponding to $M_{200} = 1.6^{+0.5}_{-0.4} \times 10^{12}$ M$_\odot$). This estimate is in good agreement with recently published independent mass estimates based on the kinematics of more distant halo stars and the satellite galaxy Leo I.

Diffuse interstellar bands are usually observed in spectra of hot stars, where interstellar lines are rarely blended with stellar ones. The need for hot stars is a strong limitation in the number of sightlines we can observe and the distribution of sightlines in the Galaxy, as hot stars are rare and concentrated in the Galactic plane. We are introducing a new method, where interstellar lines can be observed in spectra of cool stars in large spectroscopic surveys. The method is completely automated and does not require prior knowledge of the stellar parameters. If known, the stellar parameters only reduce the computational time and are not involved in the extraction of the interstellar spectrum. The main step in extracting interstellar lines is a construction of the stellar spectrum, which is in our method done by finding other observed spectra that lack interstellar features and are otherwise very similar to the spectrum in question. Such spectra are then combined into a single stellar spectrum template, which matches the stellar component in an observed spectrum. We demonstrate the performance of this new method on a sample of 482,430 spectra observed in RAVE survey. However, many spectra have to be combined (48 on average) in order to achieve a S/N ratio high enough to measure the DIB's profile, hence limiting the spatial information about the ISM. Only one strong interstellar line is included in the RAVE spectral range, a diffuse interstellar band at 8620 \AA. We compare its equivalent width with extinction maps and with Bayesian reddening, calculated for individual stars, and provide a linear relation between the equivalent width and reddening. Separately from the introduced method, we calculate equivalent widths of the diffuse interstellar band in spectra of hot stars with known extinction and compare all three linear relations with each other and with relations from the literature.

We present the stellar atmospheric parameters (effective temperature, surface gravity, overall metallicity), radial velocities, individual abundances and distances determined for 425 561 stars, which constitute the fourth public data release of the RAdial Velocity Experiment (RAVE). The stellar atmospheric parameters are computed using a new pipeline, based on the algorithms of MATISSE and DEGAS. The spectral degeneracies and the 2MASS photometric information are now better taken into consideration, improving the parameter determination compared to the previous RAVE data releases. The individual abundances for six elements (magnesium, aluminum, silicon, titanium, iron and nickel) are also given, based on a special-purpose pipeline which is also improved compared to that available for the RAVE DR3 and Chemical DR1 data releases. Together with photometric information and proper motions, these data can be retrieved from the RAVE collaboration website and the Vizier database.

The RAVE survey, combined with proper motions and distance estimates, can be used to study in detail stellar kinematics in the extended solar neighbourhood (solar suburb). Using the red clump, we examine the mean velocity components in 3D between an R of 6 and 10 kpc and a Z of -2 to 2 kpc, concentrating on North-South differences. Simple parametric fits to the R, Z trends for VPHI and the velocity dispersions are presented. We confirm the recently discovered gradient in mean Galactocentric radial velocity, VR, finding that the gradient is more marked below the plane, with a Z gradient also present. The vertical velocity, VZ, also shows clear structure, with indications of a rarefaction-compression pattern, suggestive of wave-like behaviour. We perform a rigorous error analysis, tracing sources of both systematic and random errors. We confirm the North-South differences in VR and VZ along the line-of-sight, with the VR estimated independent of the proper motions. The complex three-dimensional structure of velocity space presents challenges for future modelling of the Galactic disk, with the Galactic bar, spiral arms and excitation of wave-like structures all probably playing a role.

Our Galaxy is a complex machine in which several processes operate simultaneously: metal-poor gas is accreted, is chemically enriched by dying stars, and then drifts inwards, surrendering its angular momentum to stars; new stars are formed on nearly circular orbits in the equatorial plane and then diffuse through orbit space to eccentric and inclined orbits; the central stellar bar surrenders angular momentum to the surrounding disc and dark halo while acquiring angular momentum from inspiralling gas; the outer parts of the disc are constantly disturbed by satellite objects, both luminous and dark, as they sweep through pericentre. We review the conceptual tools required to bring these complex happenings into focus. Our first concern must be the construction of equilibrium models of the Galaxy, for upon these hang our hopes of determining the Galaxy's mean gravitational field, which is required for every subsequent step. Ideally our equilibrium model should be formulated so that the secular evolution of the system can be modelled with perturbation theory. Such theory can be used to understand how stars diffuse through orbit space from either the thin gas disc in which we presume disc stars formed, or the debris of an accreted object, the presumed origin of many halo stars. Coupling this understanding to the still very uncertain predictions of the theory of stellar evolution and nucleosynthesis, we can finally extract a complete model of the chemodynamic evolution of our reasonably generic Galaxy. We discuss the relation of such a model to cosmological simulations of galaxy formation, which provide general guidance but cannot be relied on for quantitative detail.

Equilibrium dynamical models are essential tools for extracting science from surveys of our Galaxy. We show how models can be tested with data from a survey before the survey's selection function has been determined. We illustrate the application of this method by presenting some results for the RAVE survey. We extend our published analytic distribution functions to include chemistry and fit the chosen functional form to a combination of the Geneva--Copenhagen survey (GCS) and a sample of G-dwarfs observed at z~1.75 kpc by the SEGUE survey. By including solid dynamics we are able to predict the contribution that the thick disc/halo stars surveyed by SEGUE should make to the GCS survey. We show that the measured [Fe/H] distribution from the GCS includes many fewer stars at [Fe/H]<-0.6 than are predicted. The problem is more likely to lie in discordant abundance scales than with incorrect dynamics.

We consider the problem of determining the Galaxy's gravitational potential from a star catalogue. We show that orbit-based approaches to this problem suffer from unacceptable numerical noise deriving from the use of only a finite number of orbits. An alternative approach, which requires an ability to determine the model's phase-space density at predetermined positions and velocities, has a level of numerical noise that lies well below the intrinsic uncertainty associated with the finite size of the catalogue analysed. A catalogue of 10000 stars brighter than V=17 and distributed over the sky at b>30 degrees enables us to determine the scaleheight of the disc that contributes to the potential with an uncertainty below 20pc if the catalogue gives proper motions, line-of-sight velocities and parallaxes with errors typical of the Gaia Catalogue, rising to 36pc if only proper motions are available. The uncertainty in the disc's scalelength is significantly smaller than 0.25kpc.

Tidal streams don't, in general, delineate orbits. A stream-orbit misalignment is expected to lead to biases when using orbit-fitting to constrain models for the Galactic potential. In this first of two papers we discuss the expected magnitude of the misalignment and the resulting dangers of using orbit-fitting algorithms to constrain the potential. We summarize data for known streams which should prove useful for constraining the Galactic potential, and compute their actions in a realistic Galactic potential. We go on to discuss the formation of tidal streams in angle-action space, and explain why, in general, streams do not delineate orbits. The magnitude of the stream-orbit misalignment is quantified for a logarithmic potential and a multi-component Galactic potential. Specifically, we focus on the expected misalignment for the known streams. By introducing a two-parameter family of realistic Galactic potentials we demonstrate that assuming these streams delineate orbits can lead to order one errors in the halo flattening and halo-to-disc force ratio at the Sun. We present a discussion of the dependence of these results on the progenitor mass, and demonstrate that the misalignment is mass-independent for the range of masses of observed streams. Hence, orbit-fitting does not yield better constraints on the potential if one uses narrower, lower-mass streams.

In the first of these two papers we demonstrated that assuming streams delineate orbits can lead to order one errors in potential parameters for realistic Galactic potentials. Motivated by the need for an improvement on orbit-fitting, we now present an algorithm for constraining the Galactic potential using tidal streams without assuming that streams delineate orbits. This approach is independent of the progenitor mass so is valid for all observed tidal streams. The method makes heavy use of angle-action variables and seeks the potential which recovers the expected correlations in angle space. We demonstrate that the method can correctly recover the parameters of a simple two-parameter logarithmic potential by analysing an N-body simulation of a stream. We investigate the magnitude of the errors in observational data for which the method can still recover the correct potential and compare this to current and future errors in data. The errors in the observables of individual stars for current and near future data are shown to be too large for the direct use of this method, but when the data are averaged in bins on the sky, the resulting averaged data are accurate enough to constrain correctly the potential parameters for achievable observational errors. From pseudo-data with errors comparable to those that will be furnished in the era of Gaia (20 per cent distance errors, 1.2 mas/yr proper motion errors, and 10 km/s line-of-sight velocity errors) we recover the circular velocity, V_c=220 km/s, and the flattening of the potential, q=0.9, to be V_c=223+/-10km/s and q=0.91+/-0.09.

We show that the ultraviolet absorption features, newly discovered in HST spectra, are consistent with being formed in a layer that extends a few kpc above the disk of the Milky Way. In this interface between the disk and the Galactic corona, high-metallicity gas ejected from the disk by supernova feedback can mix efficiently with the virial-temperature coronal material. The mixing process triggers the cooling of the lower corona down to temperatures encompassing the characteristic range of the observed absorption features, producing a net supernova-driven gas accretion onto the disk at a rate of a few Msun/yr. We speculate that this mechanism explains how the hot-mode of cosmological accretion feeds star formation in galactic disks.

Using a new algorithm for estimating the actions of orbits a parametrised distribution function is automatically fitted to observational data for the solar neighbourhood. We adopt a gravitational potential that is generated by three discs (gas and both thin and thick stellar discs), a bulge and a dark halo, and fit the thin-disc component of the distribution function to the solar-neighbourhood velocity distribution from the Geneva-Copenhagen Survey. We find that the disc's vertical density profile is in good agreement with data at z<~500 pc. The thick-disc component of the distribution function is then used to extend the fit to data from Gilmore & Reid (1983) for z<~2.5 kpc. The resulting model predicts excellent fits to the profile of the vertical velocity dispersion \sigma_z(z) from the RAVE survey and to the distribution of v_\phi velocity components at |z|~1kpc from the SDSS survey. The ability of this model to predict successfully data that was not used in the fitting process suggests that the adopted gravitational potential (which is close to a maximum-disc potential) is close to the true one. We show that if another plausible potential is used, the predicted values of \sigma_z are too large. The models imply that in contrast to the thin disc, the thick disc has to be hotter vertically than radially, a prediction that it will be possible to test in the near future. When the model parameters are adjusted in an unconstrained manner, there is a tendency to produce models that predict unexpected radial variations in quantities such as scale height. This finding suggests that to constrain these models adequately one needs data that extends significantly beyond the solar cylinder. The models presented in this paper might prove useful to the interpretation of data for external galaxies that has been taken with an integral field unit.

We give an algorithm for the economical calculation of angles and actions for stars in axisymmetric potentials. We test the algorithm by integrating orbits in a realistic model of the Galactic potential, and find that, even for orbits characteristic of thick-disc stars, the errors in the actions are typically smaller than 2 percent. We describe a scheme for obtaining actions by interpolation on tabulated values that significantly accelerates the process of calculating observables quantities, such as density and velocity moments, from a distribution function.

Using the RAVE survey, we recently brought to light a gradient in the mean galactocentric radial velocity of stars in the extended solar neighbourhood. This gradient likely originates from non-axisymmetric perturbations of the potential, among which a perturbation by spiral arms is a possible explanation. Here, we apply the traditional density wave theory and analytically model the radial component of the two-dimensional velocity field. Provided that the radial velocity gradient is caused by relatively long-lived spiral arms that can affect stars substantially above the plane, this analytic model provides new independent estimates for the parameters of the Milky Way spiral structure. Our analysis favours a two-armed perturbation with the Sun close to the inner ultra-harmonic 4:1 resonance, with a pattern speed \Omega_p=18.6^{+0.3}_{-0.2} km/s/kpc and a small amplitude A=0.55 \pm 0.02% of the background potential (14% of the background density). This model can serve as a basis for numerical simulations in three dimensions, additionally including a possible influence of the galactic bar and/or other non-axisymmetric modes.

The text of lectures to the 2011 Tenerife Winter School. The School's theme was "Secular Evolution of Galaxies" and my task was to present the underlying stellar-dynamical theory. Other lecturers were speaking on the role of bars and chemical evolution, so these topics are avoided here. We start with an account of the connections between isolating integrals, quasiperiodicity and angle-action variables - these variables played a unifying role throughout the lectures. This leads on to the phenomenon of resonant trapping and how this can lead to chaos in cuspy potentials and phase-space mixing in slowly evolving potentials. Surfaces of section and frequency analysis are introduced as diagnostics of phase-space structure. Real galactic potentials include a fluctuating part that drives the system towards unattainable thermal equilibrium. Two-body encounters are only one source of fluctuations, and all fluctuations will drive similar evolution. We derive the orbit-averaged Fokker-Planck equation and relations that hold between the second-order diffusion coefficients and both the power spectrum of the fluctuations and the first-order diffusion coefficients. From the observed heating of the solar neighbourhood we show that the second-order diffusion coefficients must scale as J^{1/2}. We show that periodic spiral structure shifts angular momentum outwards, heating at the Lindblad resonances and mixing at corotation. The equation that would yield the normal modes of a stellar disc is first derived and then used to discuss the propagation of tightly-wound spiral waves. The winding up of such waves is explains why cool stellar discs are responsive systems that amplify ambient noise. An explanation is offered of why the Lin-Shu-Kalnajs dispersion relation and even global normal-mode calculations provide a very incomplete understanding of the dynamics of stellar discs.

We present a new method for detecting and correcting systematic errors in the distances to stars when both proper motions and line-of-sight velocities are available. The method, which is applicable for samples of 200 or more stars that have a significant extension on the sky, exploits correlations between the measured U, V and W velocity components that are introduced by distance errors. We deliver a formalism to describe and interpret the specific imprints of distance errors including spurious velocity correlations and shifts of mean motion in a sample. We take into account correlations introduced by measurement errors, Galactic rotation and changes in the orientation of the velocity ellipsoid with position in the Galaxy. Tests on pseudodata show that the method is more robust and sensitive than traditional approaches to this problem. We investigate approaches to characterising the probability distribution of distance errors, in addition to the mean distance error, which is the main theme of the paper. Stars with the most overestimated distances bias our estimate of the overall distance scale, leading to the corrected distances being slightly too small. We give a formula that can be used to correct for this effect. We apply the method to samples of stars from the SEGUE survey, exploring optimal gravity cuts, sample contamination, and correcting the used distance relations.

We consider what is the best way to extract science from large surveys of the Milky Way galaxy. The diversity of data gathered in these surveys, together with our position within the Galaxy, imply that science must be extracted by fitting dynamical models to the data in the space of the observables. Models based on orbital tori promise to be superior for this task than traditional types of models, such as N-body models and Schwarzschild models. A formalism that allows such models to be fitted to data is developed and tested on pseudodata of varying richness.

Accretion of fresh gas at a rate of ~ 1 M_{sun} yr^{-1} is necessary in star-forming disc galaxies, such as the Milky Way, in order to sustain their star-formation rates. In this work we present the results of a new hydrodynamic simulation supporting the scenario in which the gas required for star formation is drawn from the hot corona that surrounds the star-forming disc. In particular, the cooling of this hot gas and its accretion on to the disc are caused by the passage of cold galactic fountain clouds through the corona.

Cosmology requires at least half of the baryons in the Universe to be in the intergalactic medium, much of which is believed to form hot coronae around galaxies. Star-forming galaxies must be accreting from their coronae. HI observations of external galaxies show that they have HI halos associated with star formation. These halos are naturally modelled as ensembles of clouds driven up by supernova bubbles. These models can fit the data successfully only if clouds exchange mass and momentum with the corona. As a cloud orbits, it is ablated and forms a turbulent wake where cold high-metallicity gas mixes with hot coronal gas causing the prompt cooling of the latter. As a consequence the total mass of HI increases. This model has recently been used to model the Leiden-Argentina-Bonn survey of Galactic HI. The values of the model's parameters that are required to model NGC 891, NGC 2403 and our Galaxy show a remarkable degree of consistency, despite the very different natures of the two external galaxies and the dramatic difference in the nature of the data for our Galaxy and the external galaxies. The parameter values are also consistent with hydrodynamical simulations of the ablation of individual clouds. The model predicts that a galaxy that loses its cool-gas disc for instance through a major merger cannot reform it from its corona; it can return to steady star formation only if it can capture a large body of cool gas, for example by accreting a gas-rich dwarf. Thus the model explains how major mergers can make galaxies "red and dead."

In a disc galaxy the distribution of azimuthal components of velocity is very skew. In the past this skewness has been modelled by superposed Gaussians. We use dynamical arguments to derive an analytic formula that can be fitted to observed velocity distributions, and validate it by fits to the velocities derived from a dynamically rigorous model, and to a sample of local stars with accurate space velocities. Our formula is much easier to use than a full distribution function. It has fewer parameters than a multi-Gaussian fit, and the best-fitting model parameters give insight into the underlying disc dynamics. In particular, once the azimuthal velocities of a sample have been successfully fitted, the apparatus provides a prediction for the corresponding distribution of radial velocities vR . An effective formula like ours is invaluable when fitting to data for stars at some distance from the Sun because it enables one to make proper allowance for the errors in distance and proper motion when determining the underlying disc kinematics. The derivation of our formula elucidates the way the horizontal and vertical motions are closely intertwined, and makes it evident that no stellar population can have a scale height and vertical velocity dispersions that are simultaneously independent of radius. We show that the oscillation of a star perpendicular to the Galactic plane modifies the effective potential in which the star moves radially in such a way that the more vertical energy a star has, the larger is the mean radius of its orbit.

We use a model of the Galactic fountain to simulate the neutral-hydrogen emission of the Milky Way Galaxy. The model was developed to account for data on external galaxies with sensitive HI data. For appropriate parameter values, the model reproduces well the HI emission observed at Intermediate Velocities. The optimal parameters imply that cool gas is ionised as it is blasted out of the disc, but becomes neutral when its vertical velocity has been reduced by ~30 per cent. The parameters also imply that cooling of coronal gas in the wakes of fountain clouds transfers gas from the virial-temperature corona to the disc at ~2 Mo/yr. This rate agrees, to within the uncertainties with the accretion rate required to sustain the Galaxy's star formation without depleting the supply of interstellar gas. We predict the radial profile of accretion, which is an important input for models of Galactic chemical evolution. The parameter values required for the model to fit the Galaxy's HI data are in excellent agreement with values estimated from external galaxies and hydrodynamical studies of cloud-corona interaction. Our model does not reproduce the observed HI emission at High Velocities, consistent with High Velocity Clouds being extragalactic in origin. If our model is correct, the structure of the Galaxy's outer HI disc differs materially from that used previously to infer the distribution of dark matter on the Galaxy's outskirts.

We apply the method of Burnett & Binney (2010) for the determination of stellar distances and parameters to the internal catalogue of the Radial Velocity Experiment (Steinmetz et al. 2006). Subsamples of stars that either have Hipparcos parallaxes or belong to well-studied clusters, inspire confidence in the formal errors. Distances to dwarfs cooler than ~6000 K appear to be unbiased, but those to hotter dwarfs tend to be too small by ~10% of the formal errors. Distances to giants tend to be too large by about the same amount. The median distance error in the whole sample of 216,000 stars is 28% and the error distribution is similar for both giants and dwarfs. Roughly half the stars in the RAVE survey are giants. The giant fraction is largest at low latitudes and in directions towards the Galactic Centre. Near the plane the metallicity distribution is remarkably narrow and centred on [M/H]-0.04 dex; with increasing |z| it broadens out and its median moves to [M/H] ~ -0.5. Mean age as a function of distance from the Galactic centre and distance |z| from the Galactic plane shows the anticipated increase in mean age with |z|.

In galaxies like the Milky Way, cold (~ 10^4 K) gas ejected from the disc by stellar activity (the so-called galactic-fountain gas) is expected to interact with the virial-temperature (~ 10^6 K) gas of the corona. The associated transfer of momentum between cold and hot gas has important consequences for the dynamics of both gas phases. We quantify the effects of such an interaction using hydrodynamical simulations of cold clouds travelling through a hot medium at different relative velocities. Our main finding is that there is a velocity threshold between clouds and corona, of about 75 km/s, below which the hot gas ceases to absorb momentum from the cold clouds. It follows that in a disc galaxy like the Milky Way a static corona would be rapidly accelerated: the corona is expected to rotate and to lag, in the inner regions, by ~ 80-120 km/s with respect to the cold disc. We also show how the existence of this velocity threshold can explain the observed kinematics of the cold extra-planar gas.

We present the third data release of the RAdial Velocity Experiment (RAVE) which is the first milestone of the RAVE project, releasing the full pilot survey. The catalog contains 83,072 radial velocity measurements for 77,461 stars in the southern celestial hemisphere, as well as stellar parameters for 39,833 stars. This paper describes the content of the new release, the new processing pipeline, as well as an updated calibration for the metallicity based upon the observation of additional standard stars. Spectra will be made available in a future release. The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

Our knowledge of the Galaxy is being revolutionised by a series of photometric, spectroscopic and astrometric surveys. Already an enormous body of data is available from completed surveys, and data of ever increasing quality and richness will accrue at least until the end of this decade. To extract science from these surveys we need a class of models that can give probability density functions in the space of the observables of a survey -- we should not attempt to "invert" the data from the space of observables into the physical space of the Galaxy. Currently just one class of model has the required capability, so-called "torus models". A pilot application of torus models to understanding the structure of the Galaxy's thin and thick discs has already produced two significant results: a major revision of our best estimate of the Sun's velocity with respect to the Local Standard of Rest, and a successful prediction of the way in which the vertical velocity dispersion in the disc varies with distance from the Galactic plane.

We present here a fast code for creating a synthetic survey of the Milky Way. Given one or more color-magnitude bounds, a survey size and geometry, the code returns a catalog of stars in accordance with a given model of the Milky Way. The model can be specified by a set of density distributions or as an N-body realization. We provide fast and efficient algorithms for sampling both types of models. As compared to earlier sampling schemes which generate stars at specified locations along a line of sight, our scheme can generate a continuous and smooth distribution of stars over any given volume. The code is quite general and flexible and can accept input in the form of a star formation rate, age metallicity relation, age velocity dispersion relation and analytic density distribution functions. Theoretical isochrones are then used to generate a catalog of stars and support is available for a wide range of photometric bands. As a concrete example we implement the Besancon Milky Way model for the disc. For the stellar halo we employ the simulated stellar halo N-body models of Bullock & Johnston (2005). In order to sample N-body models, we present a scheme that disperses the stars spawned by an N-body particle, in such a way that the phase space density of the spawned stars is consistent with that of the N-body particles. The code is ideally suited to generating synthetic data sets that mimic near future wide area surveys such as GAIA, LSST and HERMES. As an application we study the prospect of identifying structures in the stellar halo with a simulated GAIA survey. We plan to make the code publicly available at http://galaxia.sourceforge.net.

Stars near the Sun oscillate both horizontally and vertically. In Paper I the coupling between these motions was modelled by determining the horizontal motion without reference to the vertical motion, and recovering the coupling by assuming that the vertical action is adiabatically conserved as the star oscillates horizontally. Here we show that, although the assumption of adiabatic invariance works well, more accurate results can be obtained by taking the vertical action into account when calculating the horizontal motion. We use orbital tori to present a simple but fairly realistic model of the Galaxy's discs in which the motion of stars is handled rigorously, without decomposing it into horizontal and vertical components. We examine the ability of the adiabatic approximation to calculate the model's observables, and find that it performs perfectly in the plane, but errs slightly away from the plane. When the new correction to the adiabatic approximation is used, the density, streaming velocity and velocity dispersions are in error by less than 10 per cent for distances up to $2.5\kpc$ from the Sun. The torus-based model reveals that at locations above the plane the long axis of the velocity ellipsoid points almost to the Galactic centre, even though the model potential is significantly flattened. This result contradicts the widespread belief that the shape of the Galaxy's potential can be strongly constrained by the orientation of velocity ellipsoid near the Sun. An analysis of orbits reveals that in a general potential the orientation of the velocity ellipsoid depends on the structure of the model's distribution function as much as on its gravitational potential, contrary to what is the case for Staeckel potentials. We argue that the adiabatic approximation will provide a valuable complement to torus-based models in the interpretation of current surveys of the Galaxy.

We present an analysis of the mechanics of thin streams, which are formed following the tidal disruption of cold, low-mass clusters in the potential of a massive host galaxy. The analysis makes extensive use of action-angle variables, in which the physics of stream formation and evolution is expressed in a particularly simple form. We demonstrate the formation of streams by considering examples in both spherical and flattened potentials, and we find that the action-space structures formed in each take on a consistent and characteristic shape. We demonstrate that tidal streams formed in realistic galaxy potentials are poorly represented by single orbits, contrary to what is often assumed. We further demonstrate that attempting to constrain the parameters of the Galactic potential by fitting orbits to such streams can lead to significant systematic error. However, we show that it is possible to predict accurately the track of streams from simple models of the action-space distribution of the disrupted cluster.

We identify a new, nearby (0.5 < d < 10 kpc) stream in data from the RAdial Velocity Experiment (RAVE). As the majority of stars in the stream lie in the constellation of Aquarius we name it the Aquarius Stream. We identify 15 members of the stream lying between 30 < l < 75 and -70< b <-50, with heliocentric line-of-sight velocities V_los~-200 km/s. The members are outliers in the radial velocity distribution, and the overdensity is statistically significant when compared to mock samples created with both the Besan\c{c}on Galaxy model and newly-developed code Galaxia. The metallicity distribution function and isochrone fit in the log g - T_eff plane suggest the stream consists of a 10 Gyr old population with [m/H]~-1.0. We explore relations to other streams and substructures, finding the stream cannot be identified with known structures: it is a new, nearby substructure in the Galaxy's halo. Using a simple dynamical model of a dissolving satellite galaxy we account for the localization of the stream. We find that the stream is dynamically young and therefore likely the debris of a recently disrupted dwarf galaxy or globular cluster. The Aquarius stream is thus a specimen of ongoing hierarchical Galaxy formation, rare for being right in the solar suburb.

Using a sample of 213,713 stars from the Radial Velocity Experiment (RAVE) survey, limited to a distance of 2 kpc from the Sun and to |z|<1 kpc, we report the detection of a velocity gradient of disc stars in the fourth quadrant, directed radially from the Galactic centre. In the direction of the Galactic centre, we apply a simple method independent of stellar proper motions and of Galactic parameters to assess the existence of this gradient in the RAVE data. This velocity gradient corresponds to |K+C| < 3 km/s/kpc, where K and C are the Oort constants measuring the local divergence and radial shear of the velocity field, respectively. In order to illustrate the effect, assuming a zero radial velocity of the Local Standard of Rest we then reconstruct the two-dimensional Galactocentric velocity maps using two different sets of proper motions and photometric distances based either on isochrone fitting or on K-band magnitudes, and considering two sets of values for the Galactocentric radius of the Sun and local circular speed. Further observational confirmation of our finding with line-of-sight velocities of stars at low latitudes, together with further modelling, should help constrain the non-axisymmetric components of the Galactic potential, including the bar, the spiral arms and possibly the ellipticity of the dark halo.

Very metal-poor stars are of obvious importance for many problems in chemical evolution, star formation, and galaxy evolution. Finding complete samples of such stars which are also bright enough to allow high-precision individual analyses is of considerable interest. We demonstrate here that stars with iron abundances [Fe/H] < -2 dex, and down to below -4 dex, can be efficiently identified within the Radial Velocity Experiment (RAVE) survey of bright stars, without requiring additional confirmatory observations. We determine a calibration of the equivalent width of the Calcium triplet lines measured from the RAVE spectra onto true [Fe/H], using high spectral resolution data for a subset of the stars. These RAVE iron abundances are accurate enough to obviate the need for confirmatory higher-resolution spectroscopy. Our initial study has identified 631 stars with [Fe/H] <= -2, from a RAVE database containing approximately 200,000 stars. This RAVE-based sample is complete for stars with [Fe/H] < -2.5, allowing statistical sample analysis. We identify three stars with [Fe/H] <= -4. Of these, one was already known to be `ultra metal-poor', one is a known carbon-enhanced metal-poor star, but we obtain [Fe/H]= -4.0, rather than the published [Fe/H]=-3.3, and derive [C/Fe] = +0.9, and [N/Fe] = +3.2, and the third is at the limit of our S/N. RAVE observations are on-going and should prove to be a rich source of bright, easily studied, very metal-poor stars.

The RAdial Velocity Experiment (RAVE) is a spectroscopic survey of the Milky Way. We use the subsample of spectra with spectroscopically determined values of stellar parameters to determine the distances to these stars. The list currently contains 235,064 high quality spectra which show no peculiarities and belong to 210,872 different stars. The numbers will grow as the RAVE survey progresses. The public version of the catalog will be made available through the CDS services along with the ongoing RAVE public data releases. The distances are determined with a method based on the work by Breddels et al.~(2010). Here we assume that the star undergoes a standard stellar evolution and that its spectrum shows no peculiarities. The refinements include: the use of either of the three isochrone sets, a better account of the stellar ages and masses, use of more realistic errors of stellar parameter values, and application to a larger dataset. The derived distances of both dwarfs and giants match within ~21% to the astrometric distances of Hipparcos stars and to the distances of observed members of open and globular clusters. Multiple observations of a fraction of RAVE stars show that repeatability of the derived distances is even better, with half of the objects showing a distance scatter of \simlt 11%. RAVE dwarfs are ~300 pc from the Sun, and giants are at distances of 1 to 2 kpc, and up to 10 kpc. This places the RAVE dataset between the more local Geneva-Copenhagen survey and the more distant and fainter SDSS sample. As such it is ideal to address some of the fundamental questions of Galactic structure and evolution in the pre-Gaia era. Individual applications are left to separate papers, here we show that the full 6-dimensional information on position and velocity is accurate enough to discuss the vertical structure and kinematic properties of the thin and thick disks.

We study the eccentricity distribution of a thick disc sample of stars observed in the Radial Velocity Experiment (RAVE) and compare it to that expected in four simulations of thick disc formation in the literature (accretion of satellites, heating of a primordial thin disc during a merger, radial migration, and gas-rich mergers), as compiled by Sales et al. (2009). We find that the distribution of our sample is peaked at low eccentricities and falls off smoothly and rather steeply to high eccentricities. This distribution is fairly robust to changes in distances, thin disc contamination, and the particular thick disc sample used. Our results are inconsistent with what is expected for the pure accretion simulation, since we find that the dynamics of local thick disc stars implies that the majority must have formed "in situ". Of the remaining models explored, the eccentricity distribution of our stars appears to be most consistent with the gas-rich merger case.

We report on the discovery of new members of nearby young moving groups, exploiting the full power of combining the RAVE survey with several stellar age diagnostic methods and follow-up high-resolution optical spectroscopy. The results include the identification of one new and five likely members of the beta Pictoris moving group, ranging from spectral types F9 to M4 with the majority being M dwarfs, one K7 likely member of the epsilon Cha group and two stars in the Tuc-Hor association. Based on the positive identifications we foreshadow a great potential of the RAVE database in progressing toward a full census of young moving groups in the solar neighbourhood.

Theories of thick disk formation can be differentiated by measurements of stellar elemental abundances. We have undertaken a study of metal-poor stars selected from the RAVE spectroscopic survey of bright stars to establish whether or not there is a significant population of metal-poor thick-disk stars ([Fe/H] <~ -1.0) and to measure their elemental abundances. In this paper, we present abundances of four alpha-elements (Mg, Si, Ca, Ti) and iron for a subsample of 212 RGB and 31 RC/HB stars from this study. We find that the [alpha/Fe] ratios are enhanced implying that enrichment proceeded by purely core-collapse supernovae. This requires that star formation in each star forming region had a short duration. The relative lack of scatter in the [alpha/Fe] ratios implies good mixing in the ISM prior to star formation. In addition, the ratios resemble that of the halo, indicating that the halo and thick disk share a similar massive star IMF. We conclude that the alpha enhancement of the metal-poor thick disk implies that direct accretion of stars from dwarf galaxies similar to surviving dwarf galaxies today did not play a major role in the formation of the thick disk.

We consider the problem of self-regulated heating and cooling in galaxy clusters and the implications for cluster magnetic fields and turbulence. Viscous heating of a weakly collisional magnetised plasma is regulated by the pressure anisotropy with respect to the local direction of the magnetic field. The intracluster medium is a high-beta plasma, where pressure anisotropies caused by the turbulent stresses and the consequent local changes in the magnetic field will trigger very fast microscale instabilities. We argue that the net effect of these instabilities will be to pin the pressure anisotropies at a marginal level, controlled by the plasma beta parameter. This gives rise to local heating rates that turn out to be comparable to the radiative cooling rates. Furthermore, we show that a balance between this heating and Bremsstrahlung cooling is thermally stable, unlike the often conjectured balance between cooling and thermal conduction. Given a sufficient (and probably self-regulating) supply of turbulent power, this provides a physical mechanism for mitigating cooling flows and preventing cluster core collapse. For observed density and temperature profiles, the assumed balance of viscous heating and radiative cooling allows us to predict magnetic-field strengths, turbulent velocities and turbulence scales as functions of distance from the centre. Specific predictions and comparisons with observations are given for several different clusters. Our predictions can be further tested by future observations of cluster magnetic fields and turbulent velocities.

We devise a new method for the detection of double-lined binary stars in a sample of the Radial Velocity Experiment (RAVE) survey spectra. The method is both tested against extensive simulations based on synthetic spectra, and compared to direct visual inspection of all RAVE spectra. It is based on the properties and shape of the cross-correlation function, and is able to recover ~80% of all binaries with an orbital period of order 1 day. Systems with periods up to 1 year are still within the detection reach. We have applied the method to 25,850 spectra of the RAVE second data release and found 123 double-lined binary candidates, only eight of which are already marked as binaries in the SIMBAD database. Among the candidates, there are seven that show spectral features consistent with the RS CVn type (solar type with active chromosphere) and seven that might be of W UMa type (over-contact binaries). One star, HD 101167, seems to be a triple system composed of three nearly identical G-type dwarfs. The tested classification method could also be applicable to the data of the upcoming Gaia mission.

A Bayesian approach to the determination of stellar distances from photometric and spectroscopic data is presented and tested both on pseudodata, designed to mimic data for stars observed by the RAVE survey, and on the real stars from the Geneva-Copenhagen survey. It is argued that this method is optimal in the sense that it brings to bear all available information and that its results are limited only by observational errors and the underlying physics of stars. The method simultaneously returns the metallicities, ages and masses of programme stars. Remarkably, the uncertainty in the output metallicity is typically 44 per cent smaller than the uncertainty in the input metallicity.

(Abridged) Aims:We develop a method for deriving distances from spectroscopic data and obtaining full 6D phase-space coordinates for the RAVE survey's second data release. Methods: We used stellar models combined with atmospheric properties from RAVE (Teff, logg and [Fe/H]) and (J-Ks) photometry from archival sources to derive absolute magnitudes. We are able to derive the full 6D phase-space coordinates for a large sample of RAVE stars. This method is tested with artificial data, Hipparcos trigonometric parallaxes and observations of the open cluster M67. Results: When we applied our method to a set of 16 146 stars, we found that 25% (4 037) of the stars have relative (statistical) distance errors of < 35%, while 50% (8 073) and 75% (12 110) have relative (statistical) errors smaller than 45% and 50%, respectively. Our various tests show that we can reliably estimate distances for main-sequence stars, but there is an indication of potential systematic problems with giant stars. For the main-sequence star sample (defined as those with log(g) > 4), 25% (1 744) have relative distance errors < 31%, while 50% (3 488) and 75% (5 231) have relative errors smaller than 36% and 42%, respectively. Our full dataset shows the expected decrease in the metallicity of stars as a function of distance from the Galactic plane. The known kinematic substructures in the U and V velocity components of nearby dwarf stars are apparent in our dataset, confirming the accuracy of our data and the reliability of our technique. We provide independent measurements of the orientation of the UV velocity ellipsoid and of the solar motion, and they are in very good agreement with previous work. Conclusions: The distance catalogue for the RAVE second data release is available at http://www.astro.rug.nl/~rave

It is argued that galaxies like ours sustain their star formation by transferring gas from an extensive corona to the star-forming disc. The transfer is effected by the galactic fountain -- cool clouds that are shot up from the plane to kiloparsec heights above the plane. The Kelvin-Helmholtz instability strips gas from these clouds. If the pressure and the the metallicity of the corona are high enough, the stripped gas causes a similar mass of coronal gas to condense in the cloud's wake. Hydrodynamical simulations of cloud-corona interaction are presented. These confirm the existence of a critical ablation rate above which the corona is condensed, and imply that for the likely parameters of the Galactic corona this rate lies near the actual ablation rate of clouds. In external galaxies trails of HI behind individual clouds will not be detectable, although the integrated emission from all such trails should be significant. Parts of the trails of the clouds that make up the Galaxy's fountain should be observable and may account for features in targeted 21-cm observations of individual high-velocity clouds and surveys of Galactic HI emission. Taken in conjunction with the known decline in the availability of cold infall with increasing cosmic time and halo mass, the proposed mechanism offers a promising explanation of the division of galaxies between the blue cloud to the red sequence in the colour-luminosity plane.

Star-forming disc galaxies such as the Milky Way need to accrete $\gsim$ 1 $M_{\odot}$ of gas each year to sustain their star formation. This gas accretion is likely to come from the cooling of the hot corona, however it is still not clear how this process can take place. We present simulations supporting the idea that this cooling and the subsequent accretion are caused by the passage of cold galactic-fountain clouds through the hot corona. The Kelvin-Helmholtz instability strips gas from these clouds and the stripped gas causes coronal gas to condense in the cloud's wake. For likely parameters of the Galactic corona and of typical fountain clouds we obtain a global accretion rate of the order of that required to feed the star formation.

We re-examine the stellar kinematics of the Solar neighbourhood in terms of the velocity of the Sun with respect to the local standard of rest. We show that the classical determination of its component V_sun in the direction of Galactic rotation via Stroemberg's relation is undermined by the metallicity gradient in the disc, which introduces a correlation between the colour of a group of stars and the radial gradients of its properties. Comparing the local stellar kinematics to a chemodynamical model which accounts for these effects, we obtain (U,V,W)_sun = (11.1 +/- 0.74, 12.24 +/- 0.47, 7.25 +/-0.37) km/s, with additional systematic uncertainties of ~ (1,2,0.5) km/s. In particular, V_sun is 7 km/s larger than previously estimated. The new values of solar motion are extremely insensitive to the metallicity gradient within the disc.

We reanalyse the measurements of parallax, proper motion, and line-of-sight velocity for 18 masers in high mass star-forming regions presented by Reid et al. (2009). We use a likelihood analysis to investigate the distance of the Sun from the Galactic centre, R_0, the rotational speed of the local standard of rest, v_0, and the peculiar velocity of the Sun, vsol, for various models of the rotation curve, and models which allow for a typical peculiar motion of the high mass star-forming regions. We find that these data are best fit by models with non-standard values for vsol or a net peculiar motion of the high mass star-forming regions. We argue that a correction to vsol is much more likely, and that these data support the conclusion of Binney (2009) that V_sol should be revised upwards from 5.2 km/s to 11 km/s. We find that the values of R_0 and v_0 that we determine are heavily dependent on the model we use for the rotation curve, with model-dependent estimates of R_0 ranging from 6.7 \pm 0.5kpc to 8.9 \pm 0.9kpc, and those of v_0 ranging from 200 \pm 20 km/s to 279 \pm 33 km/s. We argue that these data cannot be thought of as implying any particular values of R_0 or v_0. However, we find that v_0/R_0 is better constrained, lying in the range 29.9-31.6 km/s/kpc for all models but one.

The body of photometric and astrometric data on stars in the Galaxy has been growing very fast in recent years (Hipparcos/Tycho, OGLE-3, 2-Mass, DENIS, UCAC2, SDSS, RAVE, Pan Starrs, Hermes, ...) and in two years ESA will launch the Gaia satellite, which will measure astrometric data of unprecedented precision for a billion stars. On account of our position within the Galaxy and the complex observational biases that are built into most catalogues, dynamical models of the Galaxy are a prerequisite full exploitation of these catalogues. On account of the enormous detail in which we can observe the Galaxy, models of great sophistication are required. Moreover, in addition to models we require algorithms for observing them with the same errors and biases as occur in real observational programs, and statistical algorithms for determining the extent to which a model is compatible with a given body of data. JD5 reviewed the status of our knowledge of the Galaxy, the different ways in which we could model the Galaxy, and what will be required to extract our science goals from the data that will be on hand when the Gaia Catalogue becomes available.

Analytic distribution functions (DFs) for the Galactic disc are discussed. The DFs depend on action variables and their predictions for observable quantities are explored under the assumption that the motion perpendicular to the Galactic plane is adiabatically invariant during motion within the plane. A promising family of DFs is defined that has several adjustable parameters. A standard DF is identified by adjusting these parameters to optimise fits to the stellar density in the column above the Sun, and to the velocity distribution of nearby stars and stars ~1 kpc above the Sun. The optimum parameters imply a radial structure for the disc which is consistent with photometric studies of the Milky Way and similar galaxies, and that 20 per cent of the disc's luminosity comes from thick disc. The fits suggest that the value of the V component of the Sun's peculiar velocity should be revised upwards from 5.2 km/s to ~11 km/s. It is argued that the standard DF provides a significantly more reliable way to divide solar-neighbourhood stars into members of the thin and thick discs than is currently used. The standard DF provides predictions for surveys of stars observed at any distance from the Sun. It is anticipated that DFs of the type discussed here will provide useful starting points for much more sophisticated chemo-dynamical models of the Milky Way.

We re-examine the method of deprojection of proper motions, which has been used for finding the velocity ellipsoid of stars in the nearby Galaxy. This method is only legitimate if the lines of sight to the individual stars are uncorrelated with the stars' velocities. Very simple models are used to show that spurious results similar to ones recently reported are obtained when velocity dispersion decreases with galactocentric radius in the expected way. A scheme that compensates for this bias is proposed.

The Galaxy's stellar halo seems to be a tangle of disrupted systems that have been tidally stretched out into streams. Each stream approximately delineates an orbit in the Galactic force-field. In the first paper in this series we showed that all six phase-space coordinates of each point on an orbit can be reconstructed from the orbit's path across the sky and measurements of the line-of-sight velocity along the orbit. In this paper we complement this finding by showing that the orbit can also be reconstructed if we know proper motions along the orbit rather than the radial velocities. We also show that accurate proper motions of stream stars would enable distances to be determined to points on the stream that are independent of any assumption about the Galaxy's gravitational potential. Such "Galactic parallaxes" would be as fundamental as conventional trigonometric parallaxes, but measureable to distances ~70 times further.

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.

Virtually all massive galaxies, including our own, host central black holes ranging in mass from millions to billions of solar masses. The growth of these black holes releases vast amounts of energy that powers quasars and other weaker active galactic nuclei. A tiny fraction of this energy, if absorbed by the host galaxy, could halt star formation by heating and ejecting ambient gas. A central question in galaxy evolution is the degree to which this process has caused the decline of star formation in large elliptical galaxies, which typically have little cold gas and few young stars, unlike spiral galaxies.

We describe a technique that finds orbits through the Galaxy that are consistent with measurements of a tidal stream, taking into account the extent that tidal streams do not precisely delineate orbits. We show that if accurate line-of-sight velocities are measured along a well defined stream, the technique recovers the underlying orbit through the Galaxy and predicts the distances and proper motions along the stream to high precision. As the error bars on the location and velocities of the stream grow, the technique is able to find more and more orbits that are consistent with the data and the uncertainties in the predicted distances and proper motions increase. With radial-velocity data along a stream ~40deg long and <0.3deg wide on the sky accurate to ~1 km/s the precisions of the distances and tangential velocities along the stream are 4 percent and 5 km/s, respectively. The technique can be used to diagnose the Galactic potential: if circular-speed curve is actually flat, both a Keplerian potential and Phi(r) proportional to r are readily excluded. Given the correct radial density profile for the dark halo, the halo's mass can be determined to a precision of 5 percent.

CONTEXT: The RAVE spectroscopic survey for galactic structure and evolution obtains 8400-8800 Ang spectra at 7500 resolving power at the UK Schmidt Telescope using the 6dF multi-fiber positioner. More than 300,000 9<I<12 and |b|>25 deg southern stars have been observed to date. AIMS: This paper presents the first intrinsic examination of stellar spectra from the RAVE survey, aimed at evaluating their diagnostic potential for peculiar stars and at contributing to the general understanding of Luminous Blue Variables (LBVs). METHODS: We used the multi-epoch spectra for all seven LBVs observed, between 2005 and 2008, in the Large Magellanic Cloud (LMC) by the RAVE survey. RESULTS: We demonstrate that RAVE spectra possess significant diagnostic potential when applied to peculiar stars and, in particular, LBVs. The behaviour of the radial velocities for both emission and absorption lines, and the spectral changes between outburst and quiescence states are described and found to agree with evidence gathered at more conventional wavelengths. The wind outflow signatures and their variability are investigated, with multi-components detected in S Doradus. Photoionisation modelling of the rich emission line spectrum of R 127 shows evidence of a massive detached ionised shell that was ejected during the 1982-2000 outburst. Surface inhomogeneities in the nuclear-processed material, brought to the surface by heavy mass loss, could have been observed in S Doradus, even if alternative explanations are possible. We also detect the transition from quiescence to outburst state in R 71. Finally, our spectrum of R 84 offers one of the clearest views of its cool companion.

We examine the proposal that the HI "high-velocity" clouds (HVCs) surrounding the Milky Way and other disc galaxies form by condensation of the hot galactic corona via thermal instability. Under the assumption that the galactic corona is well represented by a non-rotating, stratified atmosphere, we find that for this formation mechanism to work the corona must have an almost perfectly flat entropy profile. In all other cases the growth of thermal perturbations is suppressed by a combination of buoyancy and thermal conduction. Even if the entropy profile were nearly flat, cold clouds with sizes smaller than 10 kpc could form in the corona of the Milky Way only at radii larger than 100 kpc, in contradiction with the determined distances of the largest HVC complexes. Clouds with sizes of a few kpc can form in the inner halo only in low-mass systems. We conclude that unless even slow rotation qualitatively changes the dynamics of a corona, thermal instability is unlikely to be a viable mechanism for formation of cold clouds around disc galaxies.

We use proper motions and parallaxes from the new reduction of Hipparcos data and Geneva-Copenhagen radial velocities for a complete sample of ~15000 main-sequence and subgiant stars, and new Padova isochrones to constrain the kinematics and star-formation history of the solar neighbourhood. We rederive the solar motion and the structure of the local velocity ellipsoids. When the principal velocity dispersions are assumed to increase with time as t^\beta, the index \beta is larger for \sigma_W (\beta_W~0.45) than for \sigma_U (\beta_U~0.31). For the three-dimensional velocity dispersion we obtain \beta=0.35. We exclude saturation of disc heating after ~3 Gyr as proposed by Quillen & Garnett(2000). Saturation after >~4 Gyr combined with an abrupt increase in velocity dispersion for the oldest stars cannot be excluded. For all our models the star-formation rate is declining, being a factor 2-7 lower now than it was at the beginning. Models in which the SFR declines exponentially favour very high disc ages between 11.5 and 13 Gyr and exclude ages below ~10.5 Gyr as they yield worse fits to the number density and velocity dispersion of red stars. Models in which the SFR is the sum of two declining exponentials representing the thin and thick discs favour ages between 10.5 and 12 Gyr with a lower limit of ~10.0 Gyr. Although in our models the star-formation rate peaked surprisingly early, the mean formation time of solar-neighbourhood stars is later than in ab-initio models of galaxy formation, probably on account of weaknesses in such models.

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 alpha 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 monotonic declining, its disc naturally splits into an alpha-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 sigma_z from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of sigma_z with age than predicted by theory. The metallicity gradient in the ISM 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 cutoff 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.

We present a measure of the inclination of the velocity ellipsoid at 1 kpc below the Galactic plane using a sample of red clump giants from the RAVE DR2 release. We find that the velocity ellipsoid is tilted towards the Galactic plane with an inclination of 7.3 +/-1.8 degree. We compare this value to computed inclinations for two mass models of the Milky Way. We find that our measurement is consistent with a short scale length of the stellar disc (Rd ~2 kpc) if the dark halo is oblate or with a long scale length (Rd~3 kpc) if the dark halo is prolate. Once combined with independent constraints on the flattening of the halo, our measurement suggests that the scale length is approximately halfway between these two extreme values, with a preferred range [2.5-2.7] kpc for a nearly spherical halo. Nevertheless, no model can be clearly ruled out. With the continuation of the RAVE survey, it will be possible to provide a strong constraint on the mass distribution of the Milky Way using refined measurements of the orientation of the velocity ellipsoid.

Angle-action coordinates are used to study the relic of an N-body simulation of a self-gravitating satellite galaxy that was released on a short-period orbit within the disc of the Galaxy. Satellite stars that lie within 1.5 kpc of the Sun are confined to a grid of patches in action space. As the relic phase-mixes for longer, the patches become smaller and more numerous. These patches can be seen even when the angle-action coordinates of an erroneous Galactic potential are used, but using the wrong potential displaces them. Diagnostic quantities constructed from the angle coordinates both allow the true potential to be identified, and the relic to be dated. Hence when the full phase space coordinates of large numbers of solar-neighbourhood stars are known, it should be possible to identify members of particular relics from the distribution of stars in an approximate action space. This would then open up the possibility of determining the time since the relic was disrupted and gaining better knowledge of the Galactic potential. The availability of angle-action coordinates for arbitrary potentials is the key to these developments. The paper includes a brief introduction to the torus technique used to generate them.

We have used spectra of hot stars from the RAVE Survey in order to investigate the visibility and properties of five diffuse interstellar bands previously reported in the literature. The RAVE spectroscopic survey for Galactic structure and kinematics records CCD spectra covering the 8400-8800 Ang wavelength region at 7500 resolving power. The spectra are obtained with the UK Schmidt at the AAO, equipped with the 6dF multi-fiber positioner. The DIB at 8620.4 Ang is by far the strongest and cleanest of all DIBs occurring within the RAVE wavelength range, with no interference by underlying absorption stellar lines in hot stars. It correlates so tightly with reddening that it turns out to be a reliable tool to measure it, following the relation E(B-V) = 2.72 (+/- 0.03) x E.W.(Ang), valid throughout the general interstellar medium of our Galaxy. The presence of a DIB at 8648 Ang is confirmed. Its intensity appears unrelated to reddening, in agreement with scanty and preliminary reports available in the literature, and its measurability is strongly compromised by severe blending with underlying stellar HeI doublet at 8649 Ang. The two weak DIBS at 8531 and 8572 Ang do not appear real and should actually be blends of underlying stellar lines. The very weak DIB at 8439 Ang cannot be resolved within the profile of the much stronger underlying hydrogen Paschen 18 stellar line.

We present the second data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities (RVs) and stellar atmosphere parameters of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). It is obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region (8,410--8,795 \AA) for southern hemisphere stars in the magnitude range 9<I<12. Following the first data release (Steinmetz et al. 2006) the current release doubles the sample of published RVs, now containing 51,829 RVs for 49,327 individual stars observed on 141 nights between April 11 2003 and March 31 2005. Comparison with external data sets shows that the new data collected since April 3 2004 show a standard deviation of 1.3 km/s, about twice better than for the first data release. For the first time this data release contains values of stellar parameters from 22,407 spectra of 21,121 individual stars. They were derived by a penalized \chi^2 method using an extensive grid of synthetic spectra calculated from the latest version of Kurucz models. From comparison with external data sets, our conservative estimates of errors of the stellar parameters (for a spectrum with S/N=40) are 400 K in temperature, 0.5 dex in gravity, and 0.2 dex in metallicity. We note however that the internal errors estimated from repeat RAVE observations of 822 stars are at least a factor 2 smaller. We demonstrate that the results show no systematic offsets if compared to values derived from photometry or complementary spectroscopic analyses. The data release includes proper motion and photometric measurements. It can be accessed via the RAVE webpage: http://www.rave-survey.org and through CDS.

We have tested a previous analytical estimate of the dynamical friction timescale in Modified Newtonian Dynamics (MOND) with fully non-linear N-body simulations. The simulations confirm that the dynamical friction timescale is significantly shorter in MOND than in equivalent Newtonian systems, i.e. systems with the same phase-space distribution of baryons and additional dark matter. An apparent conflict between this result and the long timescales determined for bars to slow and mergers to be completed in previous N-body simulations of MOND systems is explained. The confirmation of the short dynamical-friction timescale in MOND underlines the challenge that the Fornax dwarf spheroidal poses to the viability of MOND.

Recent years have seen the discovery of many tidal streams through the Galaxy. Relatively straightforward observations of a stream allow one to deduce three phase-space coordinates of an orbit. An algorithm is presented that reconstructs the missing phase-space coordinates from these data. The reconstruction starts from assumed values of the Galactic potential and a distance to one point on the orbit, but with noise-free data the condition that energy be conserved on the orbit enables one to reject incorrect assumptions. The performance of the algorithm is investigated when errors are added to the input data that are comparable to those in published data for the streams of Pal 5. It is found that the algorithm returns distances and proper motions that are accurate to of order one percent, and enables one to reject quite reasonable but incorrect trial potentials. In practical applications it will be important to minimize errors in the input data, and there is considerable scope for doing this.

We present evidence for cosmological gas accretion onto spiral galaxies in the local universe. The accretion is seen through its effects on the dynamics of the extra-planar neutral gas. The accretion rates that we estimate for two nearby spiral galaxies are of the order of their star formation rates. Our model shows that most of the extra-planar gas is produced by supernova feedback (galactic fountain) and only 10-20 % comes from accretion. The accreting material must have low specific angular momentum about the disc's spin axis, although the magnitude of the specific angular-momentum vector can be higher. We also explore the effects of a hot corona on the dynamics of the extra-planar gas and find that it is unlikely to be responsible for the observed kinematical pattern and the source of accreted gas. However, the interaction with the fountain flow should profoundly affect the hydrodynamics of the corona.

We analyze the distribution of G and K type stars towards the Galactic poles using RAVE and ELODIE radial velocities, 2MASS photometric star counts, and UCAC2 proper motions. The combination of photometric and 3D kinematic data allows us to disentangle and describe the vertical distribution of dwarfs, sub-giants and giants and their kinematics. We identify discontinuities within the kinematics and magnitude counts that separate the thin disk, thick disk and a hotter component. The respective scale heights of the thin disk and thick disk are 225$\pm$10 pc and 1048$\pm$36 pc. We also constrain the luminosity function and the kinematic distribution function. The existence of a kinematic gap between the thin and thick disks is incompatible with the thick disk having formed from the thin disk by a continuous process, such as scattering of stars by spiral arms or molecular clouds. Other mechanisms of formation of the thick disk such as `created on the spot' or smoothly `accreted' remain compatible with our findings.

We have searched for in-falling stellar streams on to the local Milky Way disc in the CORAVEL and RAVE surveys. The CORAVEL survey consists of local dwarf stars (Nordstrom et al. Geneva-Copenhagen survey) and local Famaey et al. giant stars. We select RAVE stars with radial velocities that are sensitive to the Galactic vertical space velocity (Galactic latitude b < -45 deg). Kuiper statistics have been employed to test the symmetry of the Galactic vertical velocity distribution functions in these samples for evidence of a net vertical flow that could be associated with a (tidal?) stream of stars with vertically coherent kinematics. In contrast to the `Field of Streams' found in the outer halo, we find that the local volumes of the solar neighbourhood sampled by the CORAVEL dwarfs (complete within ~3 x 10^-4 kpc^3), CORAVEL giants (complete within ~5 x 10^-2 kpc^3) and RAVE (5-15% complete within ~8 kpc^3) are devoid of any vertically coherent streams containing hundreds of stars. This is sufficiently sensitive to allow our RAVE sample to rule out the passing of the tidal stream of the disrupting Sagittarius (Sgr) dwarf galaxy through the solar neighbourhood. This agrees with the most recent determination of its orbit and dissociates it from the Helmi et al. halo stream. Our constraints on the absence of the Sgr stream near the Sun could prove a useful tool for discriminating between Galactic potential models. The lack of a net vertical flow through the solar neighbourhood in the CORAVEL giants and RAVE samples argues against the Virgo overdensity crossing the disc near the Sun. There are no vertical streams in the CORAVEL giants and RAVE samples with stellar densities >1.6 x 10^4 and 1.5 x 10^3 stars kpc^-3 respectively and therefore no evidence for locally enhanced dark matter.

In colour-magnitude diagrams most galaxies fall in either the ``blue cloud'' or the ``red sequence'', with the red sequence extending to significantly brighter magnitudes than the blue cloud. The bright-end of the red sequence comprises elliptical galaxies (Es) with boxy isophotes and luminosity profiles with shallow central cores, while fainter Es have disky isophotes and power-law inner surface-brightness (SB) profiles. An analysis of published data reveals that the centres of galaxies with power-law central SB profiles have younger stellar populations than the centres of cored galaxies. We argue that thermal evaporation of cold gas by virial-temperature gas plays an important role in determining these phenomena. In less massive galaxies, thermal evaporation is not very efficient, so significant amounts of cold gas can reach the galaxy centre and fill a central core with newly formed stars, consistent with the young stellar ages of the cusps of Es with power-law SB profiles. In more massive galaxies, cold gas is evaporated within a dynamical time, so star formation is inhibited, and a core in the stellar density profile produced by dissipationless dynamics cannot be refilled. The different observed properties of AGN in higher-mass and lower-mass ellipticals are also explained because in the former the central black holes invariably accrete hot gas, while in the latter they typically accrete cold gas. An important consequence of our results is that at the present time there cannot be blue, star-forming galaxies in the most massive galactic halos, consistent with the observed truncation of the blue cloud at L*. [abridged]

We analyzed the feedback of AGN jets on cooling flow clusters using three-dimensional AMR hydrodynamic simulations. We studied the interaction of the jet with the intracluster medium and creation of low X-ray emission cavities (Bubbles) in cluster plasma. The distribution of energy input by the jet into the system was quantified in its different forms, i.e. internal, kinetic and potential. We find that the energy associated with the bubbles, (pV + gamma pV/(gamma-1)), accounts for less than 10 percent of the jet energy.

We report new constraints on the local escape speed of our Galaxy. Our analysis is based on a sample of high velocity stars from the RAVE survey and two previously published datasets. We use cosmological simulations of disk galaxy formation to motivate our assumptions on the shape of the velocity distribution, allowing for a significantly more precise measurement of the escape velocity compared to previous studies. We find that the escape velocity lies within the range $498\kms < \ve < 608 \kms$ (90 per cent confidence), with a median likelihood of $544\kms$. The fact that $\ve^2$ is significantly greater than $2\vc^2$ (where $\vc=220\kms$ is the local circular velocity) implies that there must be a significant amount of mass exterior to the Solar circle, i.e. this convincingly demonstrates the presence of a dark halo in the Galaxy. For a simple isothermal halo, one can calculate that the minimum radial extent is $\sim58$ kpc. We use our constraints on $\ve$ to determine the mass of the Milky Way halo for three halo profiles. For example, an adiabatically contracted NFW halo model results in a virial mass of $1.42^{+1.14}_{-0.54}\times10^{12}M_\odot$ and virial radius of $305^{+66}_{-45}$ kpc (90 per cent confidence). For this model the circular velocity at the virial radius is $142^{+31}_{-21}\kms$. Although our halo masses are model dependent, we find that they are in good agreement with each other.

We examine the distribution of injected energy in three-dimensional, adaptive-grid simulations of the heating of cooling flows. We show that less than 10 percent of the injected energy goes into bubbles. Consequently, the energy input from the nucleus is underestimated by a factor of order 6 when it is taken to be given by PVgamma/(gamma-1), where P and V are the pressure and volume of the bubble, and gamma the ratio of principal specific heats.

In recent years evidence has accumulated that nearby spiral galaxies are surrounded by massive haloes of neutral and ionised gas. These gaseous haloes rotate more slowly than the disks and show inflow motions. They are clearly analogous to the High Velocity Clouds of the Milky Way. We show that these haloes cannot be produced by a galactic fountain process (supernova outflows from the disk) where the fountain gas conserves its angular momentum. Making this gas interact with a pre-existing hot corona does not solve the problem. These results point at the need for a substantial accretion of low angular momentum material from the IGM.

We analysed the long-term variability of four microquasars (GRS 1915+105, Cyg X-1, Cyg X-3, and Sco X-1) in radio and X rays. The results of our analysis indicate the existence of two distinct modes of energy output, which we refer to as the `coupled' mode and the `flaring' mode. The coupled mode is responsible for mildly fluctuating, flat-spectrum radio emission, coupled with the X-ray emission; the flaring mode produces powerful, steep-spectrum radio flares, with no significant counterpart in X rays. We find that the fraction of time spent by a typical microquasar in the flaring mode is similar to the fraction of quasars that are radio-loud. This is consistent with the hypothesis that radio-loudness of quasars is a function of the epoch at which the source is observed.

We present the first data release of the Radial Velocity Experiment (RAVE), an ambitious spectroscopic survey to measure radial velocities and stellar atmosphere parameters (temperature, metallicity, surface gravity) of up to one million stars using the 6dF multi-object spectrograph on the 1.2-m UK Schmidt Telescope of the Anglo-Australian Observatory (AAO). The RAVE program started in 2003, obtaining medium resolution spectra (median R=7,500) in the Ca-triplet region ($\lambda\lambda$ 8,410--8,795 \AA) for southern hemisphere stars drawn from the Tycho-2 and SuperCOSMOS catalogs, in the magnitude range 9<I<12. The first data release is described in this paper and contains radial velocities for 24,748 individual stars (25,274 measurements when including re-observations). Those data were obtained on 67 nights between 11 April 2003 to 03 April 2004. The total sky coverage within this data release is $\sim$4,760 square degrees. The average signal to noise ratio of the observed spectra is 29.5, and 80% of the radial velocities have uncertainties better than 3.4 km/s. Combining internal errors and zero-point errors, the mode is found to be 2 km/s. Repeat observations are used to assess the stability of our radial velocity solution, resulting in a variance of 2.8 km/s. We demonstrate that the radial velocities derived for the first data set do not show any systematic trend with color or signal to noise. The RAVE radial velocities are complemented in the data release with proper motions from Starnet 2.0, Tycho-2 and SuperCOSMOS, in addition to photometric data from the major optical and infrared catalogs (Tycho-2, USNO-B, DENIS and 2MASS). The data release can be accessed via the RAVE webpage: http://www.rave-survey.org.

Recent HI observations reveal that the discs of spiral galaxies are surrounded by extended gaseous haloes. This extra-planar gas reaches large distances (several kpc) from the disc and shows peculiar kinematics (low rotation and inflow). We have modelled the extra-planar gas as a continuous flow of material from the disc of a spiral galaxy into its halo region. The output of our models are pseudo-data cubes that can be directly compared to the HI data. We have applied these models to two spiral galaxies (NGC891 and NGC2403) known to have a substantial amount of extra-planar gas. Our models are able to reproduce accurately the vertical distribution of extra-planar gas for an energy input corresponding to a small fraction (<4%) of the energy released by supernovae. However they fail in two important aspects: 1) they do not reproduce the right gradient in rotation velocity; 2) they predict a general outflow of the extra-planar gas, contrary to what is observed. We show that neither of these difficulties can be removed if clouds are ionized and invisible at 21cm as they leave the disc but become visible at some point on their orbits. We speculate that these failures indicate the need for accreted material from the IGM that could provide the low angular momentum and inflow required.

The study of gas outside the plane of disc galaxies is crucial to understanding the circulation of material within a galaxy and between galaxies and the intergalactic environment. We present new HI observations of the edge-on galaxy NGC891, which show an extended halo component lagging behind the disc in rotation. We compare these results for NGC891 with other detections of gaseous haloes. Finally, we present a dynamical model for the formation of extra-planar gas.

Both microlensing surveys and radio-frequency observations of gas flow imply that the inner Milky Way is completely dominated by baryons, contrary to the predictions of standard cold dark matter (CDM) cosmology. We investigate the predictions of the Modified Newtonian Dynamics (MOND) formula for the Galaxy given the measured baryon distribution. Satisfactory fits to the observationally determined terminal-velocity curve are obtained for different choices of the MOND's interpolating function mu(x). However, with simple analytical forms of mu(x), the local circular speed v_c(R_0) can be as large as 220 km/s only for values of the parameter a_0 that are excluded by observations of NGC 3198. Only a numerically specified interpolating function can produce v_c(R_0)=220 km/s, which is therefore an upper limit in MOND, while the asymptotic velocity is predicted to be v_c(infty)=170+-5 km/s. The data are probably not consistent with the functional form of mu(x) that has been explored as a toy model in the framework of Bekenstein's covariant theory of gravity.

The low-frequency power spectra of the X-ray and radio emission from four microquasars suggest that two distinct modes of energy output are at work: (i) the `coupled' mode in which the X-ray and radio luminosities are closely coupled and vary only weakly, and (ii) the `flaring' mode, which dramatically boosts the radio luminosity but makes no impact on the X-ray luminosity. The systems are in the flaring mode only a few percent of the time. However, flares completely dominate the power spectrum of radio emission, with the consequence that sources in which the flaring mode occurs, such as GRS 1915+105 and Cyg X-3, have radio power spectra that lie more than an order of magnitude above the corresponding X-ray power spectra. Of the four microquasars for which we have examined data, in only one, Cyg X-1, is the flaring mode seemingly inactive. While Cyg X-1 is a black-hole candidate, one of the three flaring sources, Sco X-1, is a neutron star. Consequently, it is likely that both modes are driven by the accretion disk rather than black-hole spin. Radio imaging strongly suggests that the flaring mode involves relativistic jets. A typical microquasar is in the flaring mode a few percent of the time, which is similar to the fraction of quasars that are radio loud. Thus there may be no essential difference between radio-loud and radio-quiet quasars; radio loudness may simply be a function of the epoch at which the source is observed.

There is increasing evidence that AGN mechanical feedback is important in the energetics of cooling flows in galaxies and galaxy clusters. We investigate the implications of the variability of AGN mechanical luminosity L_m on observations of cooling flows and radio galaxies in general. It is natural to assume that l=ln(L_m/L_x) is a Gaussian process. Then L_m will be log-normally distributed at fixed cooling luminosity L_x, and the variance in a measure of L_m will increase with the time-resolution of the measure. We test the consistency of these predictions with existing data. These tests hinge on the power spectrum of l(t). Monitoring of Seyfert galaxies combined with estimates of the duty cycle of quasars imply flicker noise spectra, similar to those of microquasars. We combine a sample of sources in cooling flows that have cavities with the assumption that the average mechanical luminosity of the AGN equals L_x. Given that the mechanical luminosities are characterized by flicker noise, we find that their spectral amplitudes lie between the estimated amplitudes of quasars and the measured values for the radio luminosities of microquasars. The model together with the observation that powerful radio galaxies lie within a narrow range in optical luminosity, predicts the luminosity function of radio galaxies, in agreement with observations. Forthcoming radio surveys will test the prediction that the luminosity function turns over at about the smallest luminosities so far probed. [Abridged]

We have obtained Integral Field Spectroscopy for 23 spiral bulges using INTEGRAL on the William Herschel Telescope and SPIRAL on the Anglo-Australian Telescope. This is the first 2D survey directed solely at the bulges of spiral galaxies. Eleven galaxies of the sample do not have previous measurements of the stellar velocity dispersion (sigma*). These data are designed to complement our Space Telescope Imaging Spectrograph program for estimating black hole masses in the range 10^6-10^8M_sun using gas kinematics from nucleated disks. These observations will serve to derive the stellar dynamical bulge properties using the traditional Mgb and CaII triplets. We use both Cross Correlation and Maximum Penalized Likelihood to determine projected sigma* in these systems and present radial velocity fields, major axis rotation curves, curves of growth and sigma* fields. Using the Cross Correlation to extract the low order 2D stellar dynamics we generally see coherent radial rotation and irregular velocity dispersion fields suggesting that sigma* is a non-trivial parameter to estimate.

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.

We investigate the central regions of 23 spiral galaxies using archival NICMOS imaging and STIS spectroscopy. The sample is taken from our program to determine the masses of central massive black holes (MBH) in 54 nearby spiral galaxies. Stars are likely to contribute significantly to any dynamical central mass concentration that we find in our MBH program and this paper is part of a series to investigate the nuclear properties of these galaxies. We use the Nuker law to fit surface brightness profiles, derived from the NICMOS images, to look for nuclear star clusters and find possible extended sources in 3 of the 23 galaxies studied (13 per cent). The fact that this fraction is lower than that inferred from optical Hubble Space Telescope studies is probably due to the greater spatial resolution of those studies. Using R-H and J-H colors and equivalent widths of H-alpha emission (from the STIS spectra) we investigate the nature of the stellar population with evolutionary models. Under the assumption of hot stars ionizing the gas, as opposed to a weak AGN, we find that there are young stellar populations (~10-20 Myr) however these data do not allow us to determine what percentage of the total nuclear stellar population they form. Also, in an attempt to find any unknown AGN we use [N II] and [S II] line flux ratios (relative to H-alpha) and find tentative evidence for weak AGN in NGC 1300 and NGC 4536.

We present Space Telescope Imaging Spectrograph emission-line spectra of the central regions of the spiral galaxies NGC 1300 and NGC 2748. From the derived kinematics of the nuclear gas we have found evidence for central supermassive black holes in both galaxies. The estimated mass of the black hole in NGC 1300 is 6.6 (+6.3, -3.2) x 10^7 solar masses and in NGC 2748 is 4.4 (+3.5, -3.6) x 10^7 solar masses (both at the 95% confidence level). These two black hole mass estimates contribute to the poorly sampled low-mass end of the nuclear black hole mass spectrum.

Techniques for the construction of dynamical Galaxy models should be considered essential infrastructure that should be put in place before GAIA flies. Three possible modelling techniques are discussed. Although one of these seems to have significantly more potential than the other two, at this stage work should be done on all three. A major effort is needed to decide how to make a model consistent with a catalogue such as that which GAIA will produce. Given the complexity of the problem, it is argued that a hierarchy of models should be constructed, of ever increasing complexity and quality of fit to the data. The potential that resonances and tidal streams have to indicate how a model should be refined is briefly discussed.

I review evidence that galaxies form from gas that falls into potential wells cold, rather than from virialized gas, and that formation stops once an atmosphere of trapped virialized gas has accumulated. Disk galaxies do not have such atmospheres, so their formation is ongoing. During galaxy formation feedback is an efficient process, and the nuclear regions of disk galaxies blow winds. The cold infalling gas that drives continued star formation has a significant component of angular momentum perpendicular to that of the disk. Extraplanar gas has to be understood in the context set by nuclear outflows and cold skew-rotating cosmic infall.

[Abridged] We present a novel technique, dubbed FiEstAS, to estimate the underlying density field from a discrete set of sample points in an arbitrary multidimensional space. FiEstAS assigns a volume to each point by means of a binary tree. Density is then computed by integrating over an adaptive kernel. As a first test, we construct several Monte Carlo realizations of a Hernquist profile and recover the particle density in both real and phase space. At a given point, Poisson noise causes the unsmoothed estimates to fluctuate by a factor ~2 regardless of the number of particles. This spread can be reduced to about 1 dex (~26 per cent) by our smoothing procedure. [...] We conclude that our algorithm accurately measure the phase-space density up to the limit where discreteness effects render the simulation itself unreliable. Computationally, FiEstAS is orders of magnitude faster than the method based on Delaunay tessellation that Arad et al. employed, making it practicable to recover smoothed density estimates for sets of 10^9 points in 6 dimensions.

We present surface photometry for the central regions of a sample of 48 spiral galaxies (mostly unbarred and barred of types Sbc or Sc) observed with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. Surface brightness profiles were derived and modeled with a Nuker law. We also analyzed archival Wide Field Planetary Camera~2 images with a larger field of view, available for 18 galaxies in our sample. We modeled the extracted bulge surface brightness profiles with an exponential, a de Vaucouleurs or a Sersic profile. In agreement with previous studies, we find that bulges of Sbc galaxies fall into two categories: bulges well described by an exponential profile and those well described by an de Vaucouleurs profile. Only one galaxy requires the use of a more general Sersic profile to properly describe the bulge. Nuclear photometrically distinct components are found in ~55% of the galaxies. For those that we classify as star clusters based on their resolved extent we find absolute magnitudes that are brighter on average than those previously identified in spiral galaxies. This might be due to a bias in our sample toward star forming galaxies, combined with a trend for star forming galaxies to host brighter central clusters.

In cuspy atmospheres, jets driven by supermassive black holes (BHs) offset radiative cooling. The jets fire episodically, but often enough that the cuspy atmosphere does not move very far towards a cooling catastrophe in the intervals of jet inactivity. The ability of energy released on the sub-parsec scale of the BH to balance cooling on scales of several tens of kiloparsecs arises through a combination of the temperature sensitivity of the accretion rate and the way in which the radius of jet disruption varies with ambient density. Accretion of hot gas does not significantly increase BH masses, which are determined by periods of rapid BH growth and star formation when cold gas is briefly abundant at the galactic centre. Hot gas does not accumulate in shallow potential wells. As the Universe ages, deeper wells form, and eventually hot gas accumulates. This gas soon prevents the formation of further stars, since jets powered by the BH prevent it from cooling, and it mops up most cold infalling gas before many stars can form. Thus BHs set the upper limit to the masses of galaxies. The formation of low-mass galaxies is inhibited by a combination of photo-heating and supernova-driven galactic winds. Working in tandem these mechanisms can probably explain the profound difference between the galaxy luminosity function and the mass function of dark halos expected in the cold dark matter cosmology.

A naive extension to MOND of the standard computation of the two-body relaxation time Tb implies that Tb is comparable to the crossing time regardless of the number N of stars in the system. This computation is questionable in view of the non-linearity of MOND's field equation. A non-standard approach to the calculation of Tb is developed that can be extended to MOND whenever discreteness noise generates force fluctuations that are small compared to the mean-field force. It is shown that this approach yields standard Newtonian results for systems in which the mean density profile is either plane-parallel or spherical. In the plane-parallel case we find that in the deep-MOND regime Tbb scales with N as in the Newtonian case, but is shorter by the square of the factor by which MOND enhances the gravitational force over its Newtonian value for the same system. Application of these results to dwarf galaxies and groups and clusters of galaxies reveals that in MOND luminosity segregation should be far advanced in groups and clusters of galaxies, two body relaxation should have substantially modified the density profiles of galaxy groups, while objects with masses in excess of ~10M_sun should have spiralled to the centres of dwarf galaxies.(Abridged)

We determine the critical size l_crit of a filament of cold (T~10^4 K) gas that is in radiative equilibrium with X-ray emitting gas at temperatures T_out~10^6 - 10^8 K. Filaments smaller than l_crit will be rapidly evaporated, while longer ones will induce the condensation of the ambient medium. At fixed pressure P, l_crit increases as T_out^(11/4), while at fixed T_out it scales as 1/P. It scales as f^(1/2), where f is the factor by which the magnetic field depresses the thermal conductivity below Spitzer's benchmark value. For plausible values of f, l_crit is similar to the lengths of observed filaments. In a cluster such as Perseus, the value of l_crit increases by over an order of magnitude between the centre and a radius of 100 kpc. If the spectrum of seed filament lengths l is strongly falling with l, as is natural, then these results explain why filaments are only seen within a few kiloparsecs of the centres of clusters, and are not seen in clusters that have no cooling flow. We calculate the differential emission measure as a function of temperature for the interface between filaments and ambient gas of various temperatures. We discuss the implications of our results for the origin of the galaxy luminosity function.

Three-dimensional hydrodynamical simulations are used to investigate the structural stability of cooling flows that are episodically heated by jets from a central AGN. The radial profile of energy deposition is controlled by (a) the power of the jets, and (b) the pre-outburst density profile. A delay in the ignition of the jets causes more powerful jets to impact on a more centrally concentrated medium. The net effect is a sufficient increase in the central concentration of energy deposition to cause the post-outburst density profile to be less centrally concentrated than that of an identical cluster in which the outburst happened earlier and was weaker. These results suggest that the density profiles of cooling flows oscillate around an attracting profile, thus explaining why cooling flows are observed to have similar density profiles. The possibility is raised that powerful FR II systems are ones in which this feedback mechanism has broken down and a runaway growth of the source parameters has occurred.

The competition between CDM and MOND to account for the `missing mass' phenomena is asymmetric. MOND has clearly demonstrated that a characteristic acceleration $a_0$ underlies the data and understanding what gives rise to $a_0$ is an important task. The reason for MOND's success may lie in either the details of galaxy formation, or an advance in fundamental physics that reduces to MOND in a suitable limit. CDM has enjoyed great success on large scales. The theory cannot be definitively tested on small scales until galaxy formation has been understood because baryons either are, or possibly have been, dominant in all small-scale objects. MOND's predictive power is seriously undermined by its isolation from the rest of physics. In view of this isolation, the way forward is probably to treat CDM as an established theory to be used alongside relativity and electromagnetism in efforts to understand the formation and evolution of galaxies.

Active galactic nuclei are clearly heating gas in `cooling flows'. The effectiveness and spatial distribution of the heating are controversial. We use three-dimensional simulations on adaptive grids to study the impact on a cooling flow of weak, subrelativistic jets. The simulations show cavities and vortex rings as in the observations. The cavities are fast-expanding dynamical objects rather than buoyant bubbles as previously modelled, but shocks still remain extremely hard to detect with X-rays. At late times the cavities turn into overdensities that strongly excite the cluster's g-modes. These modes damp on a long timescale. Radial mixing is shown to be an important phenomenon, but the jets weaken the metallicity gradient only very near the centre. The central entropy density is modestly increased by the jets. We use a novel algorithm to impose the jets on the simulations.

Evidence is summarized that suggests that when a protogalaxy collapses, a fraction $f$ of its gas fails to heat to the virial temperature, where $f$ is large for haloes less massive than the value $M^*$ associated with $L^*$ galaxies. Stars and galaxies form only from the cool gas fraction. Hot gas is ejected from low-mass systems as in conventional semi-analytic models of galaxy formation. In high-mass systems it is retained but does not cool and form stars. Instead it builds up as a largely inert atmosphere, in which cooling is inhibited by an episodically active galactic nucleus. Cold gas frequently falls into galactic haloes. In the absence of a dense atmosphere of virial-temperature gas it builds up on nearly circular orbits and forms stars. When there is a sufficiently dense hot atmosphere, cold infalling gas tends to be ablated and absorbed by the hot atmosphere before it can form stars. The picture nicely explains away the surfeit of high-luminosity galaxies that has recently plagued semi-analytic models of galaxy formation, replacing them by systems of moderate luminosity from old stars and large X-ray luminosities from hot gas.

An estimate of the convergence radius of a simulated CDM halo is obtained under the assumption that the peak phase-space density in the system is set by discreteness effects that operate prior to relaxation. The predicted convergence radii are approximately a factor 2 larger than those estimated for numerical convergence studies. A toy model is used to study the formation of sheets of the cosmic web, from which DM haloes form later. This model demonstrates the interplay between phase mixing and violent relaxation that must also be characteristic of spherical collapse. In the limit that sheets contain arbitrarily many particles, it seems that power-law profiles are established in both distance and energy. When only a finite number of particles is employed, relaxation is prematurely terminated and the power laws are broken. In a given simulation, the sheets with the highest peak phase-space densities are those that form from the longest waves. Hence simulations with little small-scale power are expected to form the cuspiest haloes.

It is now clear that AGN heat cooling flows, largely by driving winds. The winds may contain a relativistic component that generates powerful synchrotron radiation, but it is not clear that all winds do so. The spatial and temporal stability of the AGN/cooling flow interaction are discussed. Collimation of the winds probably provides spatial stability. Temporal stability may be possible only for black holes with masses above a critical value. Both the failure of cooling flows to have adiabatic cores and the existence of X-ray cavities confirm the importance of collimated outflows. I quantify the scale of the convective flow that the AGN Hydra would need to drive if it balanced radiative inward flow by outward flow parallel to the jets. At least in Virgo any such flow must be confined to r<~20 kpc. Hydrodynamical simulations suggest that AGN outbursts cannot last longer than ~25 Myr. Data for four clusters with well studied X-ray cavities suggests that heating associated with cavity formation approximately balances radiative cooling. The role of cosmic infall and the mechanism of filament formation are briefly touched on.

We emphasise the importance of the gas entropy in studying the evolution of cluster gas evolving under the influence of radiative cooling. On this basis, we develop an analytical model for this evolution. We then show that the assumptions needed for such a model are consistent with a numerical solution of the same equations. We postulate that the passive cooling phase ends when the central gas temperature falls to very low values. It follows a phase during which an unspecified mechanism heats the cluster gas. We show that in such a scenario the small number of clusters containing gas with temperatures below about 1 keV is simply a consequence of the radiative cooling.

Using only the assumption that all interactions between particles in an accretion flow are electromagnetically mediated, it is shown that the time to establish equipartition between ions and electrons is shorter than the characteristic accretion time. Consequently, two-temperature fits to the spectra of accreting objects are unphysical, and models in which significant thermal energy is carried across the event horizon are effectively ruled out.

We present HST/STIS spectra of the Sbc spiral galaxy NGC 4041 which were used to map the velocity field of the gas in its nuclear region. We detect the presence of a compact (r~0.4" ~40 pc), high surface brightness, rotating nuclear disk co-spatial with a nuclear star cluster. The disk is characterized by a rotation curve with a peak to peak amplitude of ~40 km/s and is systematically blueshifted by 10 - 20 km/s with respect to the galaxy systemic velocity. With the standard assumption of constant mass-to-light ratio and with the nuclear disk inclination taken from the outer disk, we find that a dark point mass of 1(-0.7;+0.6) 10^7 Msun is needed to reproduce the observed rotation curve. However the observed blueshift suggests the possibility that the nuclear disk could be dynamically decoupled. Following this line of reasoning we relax the standard assumptions and find that the kinematical data can be accounted for by the stellar mass provided that either the central mass-to-light ratio is increased by a factor of ~2 or that the inclination is allowed to vary. This model results in a 3 sigma upper limit of 6 10^6 Msun on the mass of any nuclear black hole. Overall, our analysis only allows us to set an upper limit of 2 10^7 Msun on the mass of the nuclear black hole. If this upper limit is taken in conjunction with an estimated bulge B magnitude of -17.7 and with a central stellar velocity dispersion of 95 km/s, then these results are not inconsistent with both the MBH-Lsph and the MBH-sigma correlations. Constraints on black hole masses in spiral galaxies of types as late as Sbc are still very scarce and therefore the present result adds an important new datapoint to our understanding of black hole demography.

We present a simple analytical model and the results of numerical simulations supporting the idea of periodical heating of cluster gas by AGN outflows. We show why, under this assumption, we are extremely unlikely to observe clusters containing gas with a temperature below about 1 keV. We review the results from numerical hydro-simulations studying the mechanisms by which AGN outflows heat the cluster gas.

A semi-analytic model of cluster cooling flows is presented. The model assumes that episodic nuclear activity followed by radiative cooling without mass-dropout cycles the cluster gas between a relatively homogeneous, nearly isothermal post-outburst state and a cuspy configuration in which a cooling catastrophe initiates the next nuclear outburst. Fitting the model to Chandra data for the Hydra cluster, a lower limit of $284\Myr$ until the next outburst of Hydra A is derived. Density, temperature and emission-measure profiles at several times prior to the cooling catastrophe are presented. It proves possible to fit the mass $M(\sigma)$ with entropy index $P\rho^{-\gamma}$ less than $\sigma$ to a simple power-law form, which is almost invariant as the cluster cools. We show that radiative cooling automatically establishes this power-law form if the entropy index was constant throughout the cluster gas at some early epoch or after an AGN activity cycle. To high precision, the central value of $\sigma$ decreases linearly in time. The fraction of clusters in a magnitude-limited sample that have gas cooler than $T$ is calculated, and is shown to be small for $T=2\keV$. Similarly, only 1 percent of clusters in such a sample contain gas with $P\rho^{-\gamma} < 2\keV\cm^2$. Entropy production in shocks is shown to be small. The entropy that is radiated from the cluster can be replaced if a few percent of the cluster gas passes through bubbles heated during an outburst of the AGN.

We show that spiral waves in galaxy discs churn the stars and gas in a manner that largely preserves the overall angular momentum distribution and leads to little increase in random motion. Changes in the angular momenta of individual stars are typically as large as ~50% over the lifetime of the disk. The changes are concentrated around the corotation radius for an individual spiral wave, but since transient waves with a wide range of pattern speeds develop in rapid succession, the entire disk is affected. This behaviour has profound consequences for the metallicity gradients with radius in both stars and gas. The ISM is also stirred by the same mechanism. We find observational support for stirring, propose a simple model for the distribution of stars over metallicity and age, and discuss other possible consequences.

It is logically possible that early two-body relaxation in simulations of cosmological clustering influences the final structure of massive clusters. Convergence studies in which mass and spatial resolution are simultaneously increased, cannot eliminate this possibility. We test the importance of two-body relaxation in cosmological simulations with simulations in which there are two species of particles. The cases of two mass ratios, sqrt(2):1 and 4:1, are investigated. Simulations are run with both a spatially fixed softening length and adaptive softening using the publicly available codes GADGET and MLAPM, respectively. The effects of two-body relaxation are detected in both the density profiles of halos and the mass function of halos. The effects are more pronounced with a fixed softening length, but even in this case they are not so large as to suggest that results obtained with one mass species are significantly affected by two-body relaxation. The simulations that use adaptive softening are less affected by two-body relaxation and produce slightly higher central densities in the largest halos. They run about three times faster than the simulations that use a fixed softening length.

Cold, dense clouds of gas have been proposed as baryonic candidates for the dark matter in Galactic haloes, and have also been invoked in the Galactic disc as an explanation for the excess faint sub-mm 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. This presents the possibility of detecting them by looking for occultations of background stars. 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 cold clouds. We compute the rate and timescale 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.

The GAIA mission will produce an extraordinary database from which we should be able to deduce not only the Galaxy's current structure, but also much of its history, and thus cast a powerful light on the way in which galaxies in general are made up of components, and of how these formed. The database can be fully exploited only by fitting to it a sophisticated model of the entire Galaxy. Steady-state models are of fundamental importance even though the Galaxy cannot be in a steady state. A very elaborate model of the Galaxy will be required to reproduce the great wealth of detail that GAIA will reveal. A systematic approach to model-building will be required if such a model is to be successfully constructed, however. The natural strategy is to proceed through a series of models of ever increasing elaborateness, and to be guided in the specification of the next model by mismatches between the data and the current model. An approach to the dynamics of systems with steady gravitational potentials that we call the `torus programme' promises to provide an appropriate framework within which to carry out the proposed modelling programme. The basic principles of this approach have been worked out in some detail and are summarized here. Some extensions will be required before the GAIA database can be successfully confronted. Other modelling techniques that might be employed are briefly examined.

The microlensing optical depth to Baade's Window constrains the minimum total mass in baryonic matter within the Solar circle to be greater than 3.9 x 10^{10} solar masses, assuming the inner Galaxy is barred with viewing angle of roughly 20 degrees. From the kinematics of solar neighbourhood stars, the local surface density of dark matter is about 30 +/- 15 solar masses per square parsec. We construct cuspy haloes normalised to the local dark matter density and calculate the circular-speed curve of the halo in the inner Galaxy. This is added in quadrature to the rotation curve provided by the stellar and ISM discs, together with a bar sufficiently massive so that the baryonic matter in the inner Galaxy reproduces the microlensing optical depth. Such models violate the observational constraint provided by the tangent-velocity data in the inner Galaxy (typically at radii 2-4 kpc). The high baryonic contribution required by the microlensing is consistent with implications from hydrodynamical modelling and the pattern speed of the Galactic bar. We conclude that the cuspy haloes favoured by the Cold Dark Matter cosmology (and its variants) are inconsistent with the observational data on the Galaxy.

We present a computer code written in C that is designed to simulate structure formation from collisionless matter. The code is purely grid-based and uses a recursively refined Cartesian grid to solve Poisson's equation for the potential, rather than obtaining the potential from a Green's function. Refinements can have arbitrary shapes and in practice closely follow the complex morphology of the density field that evolves. The timestep shortens by a factor two with each successive refinement. It is argued that an appropriate choice of softening length is of great importance and that the softening should be at all points an appropriate multiple of the local inter-particle separation. Unlike tree and P3M codes, multigrid codes automatically satisfy this requirement. We show that at early times and low densities in cosmological simulations, the softening needs to be significantly smaller relative to the inter-particle separation than in virialized regions. Tests of the ability of the code's Poisson solver to recover the gravitational fields of both virialized halos and Zel'dovich waves are presented, as are tests of the code's ability to reproduce analytic solutions for plane-wave evolution. The times required to conduct a LCDM cosmological simulation for various configurations are compared with the times required to complete the same simulation with the ART, AP3M and GADGET codes. The power spectra, halo mass functions and halo-halo correlation functions of simulations conducted with different codes are compared.

For two decades the steady-state cooling-flow model has dominated the literature of cluster and elliptical-galaxy X-ray sources. For ten years this model has been in severe difficulty from a theoretical point of view, and it is now coming under increasing pressure observationally. For two decades the steady-state cooling-flow model has dominated the literature of cluster and elliptical-galaxy X-ray sources. For ten years this model has been in severe difficulty from a theoretical point of view, and it is now coming under increasing pressure observationally. A small number of enthusiasts have argued for a radically different interpretation of the data, but had little impact on prevailing opinion because the unsteady heating picture that they advocate is extremely hard to work out in detail. Here I explain why it is difficult to extract robust observational predictions from the heating picture. Major problems include the variability of the sources, the different ways in which a bi-polar flow can impact on X-ray emission, the weakness of synchrotron emission from sub-relativistic flows, and the sensitivity of synchrotron emission to a magnetic field that is probably highly localized.

We describe an on-going HST program aimed at determining the relationship between the nuclear black hole mass and bulge mass in spiral galaxies. We have selected a volume limited sample of 54 nearby spiral galaxies for which we already have ground based emission line rotation curves, CCD surface photometry and radio maps. We are now obtaining HST/STIS longslit observations of each of the galaxies in the sample in order to determine the nuclear Halpha rotation curve at high (~0.1") spatial resolution. We will use these data to measure the unresolved dark mass concentration at the nucleus of each object. Here we show the first results from observations of objects in the sample.

Motivated by STIS observations of more than 50 nearby galactic nuclei, we consider long-slit emission-line spectra when the slit is wider than the instrumental PSF, and the target has arbitrarily large velocity gradients. The finite width of the slit generates complex patterns in the spectra that can be misinterpreted as coming from various physically distinct nuclear components, but when interpreted correctly, they can have considerable diagnostic power. For a thin disk in circular motion around a central galactic black hole (BH), a characteristic artifact occurs in the spectrum at the outer edge of the BH's sphere of influence. It betrays the presence of a BH, and allows us to develop a new method for estimating its mass, which gives higher sensitivity to BH detection than traditional methods.

We consider long-slit emission-line spectra of galactic nuclei when the slit is wider than the instrumental PSF, and the target has large velocity gradients. The finite width of the slit generates complex distributions of brightness at a given spatial point in the measured spectrum, which can be misinterpreted as coming from additional physically distinct nuclear components. We illustrate this phenomenon for the case of a thin disc in circular motion around a nuclear black hole (BH). We develop a new method for estimating the mass of the BH that exploits a feature in the spectrum at the outer edge of the BH's sphere of influence, and therefore gives higher sensitivity to BH detection than traditional methods. Moreover, with this method we can determine the black hole mass and the inclination of the surrounding disc separately, whereas the traditional approach to black-hole estimation requires two long-slit spectra to be taken. We show that with a given spectrograph, the detectability of a BH depends on the sense of rotation of the nuclear disc. We apply our method to estimate the BH mass in M84 from a publicly available spectrum, and recover a value 4 times lower than that published previously from the same data.

In the solar-neighbourhood, older stars have larger random velocities than younger ones. It is argued that the increase in velocity dispersion with time is predominantly a gradual process rather than one induced by discrete events such as minor mergers. Ephemeral spiral arms seem to be the fundamental drivers of disk heating, although scattering by giant molecular clouds plays an important moderating role. In addition to heating the disk, spiral arms cause stars' guiding centres to diffuse radially. The speed of this diffusion is currently controversial. Data from the Hipparcos satellite has made it clear that the Galaxy is by no means in a steady state. This development enormously increases the complexity of the models required to account for the data. There are preliminary indications that we see in the local phase-space distribution the dynamical footprints of long-dissolved spiral waves.

High-quality Hipparcos data for a complete sample of nearly 12000 main-sequence and subgiant stars, together with Padua isochrones, are used to constrain the star-formation history of the solar neigbourhood and the processes that stochastically accelerate disk stars. The velocity dispersion of a coeval group of stars is found to increase with time from ~8 kms at birth as t^{0.33}. In the fits, the slope of the IMF near 1 Msun proves to be degenerate with the rate at which the star-formation rate declines. If the slope of the IMF is to lie near Salpeter's value, -2.35, the star-formation rate has to be very nearly constant. The age of the solar neighbourhood is found to be 11.2+/-0.75 Gyr with remarkably little sensitivity to variations in the assumed metallicity distribution of old disk stars. This age is only a Gyr younger than the age of the oldest globular clusters when the same isochrones and distance scale are employed. It is compatible with current indications of the redshift of luminous galaxy formation only if there is a large cosmological constant. A younger age is formally excluded because it provides a poor fit to the number density of red stars. Since this density is subject to a significantly uncertain selection function, ages as low as 9 Gyr are plausible even though they lie outside our formal error bars.

An energetic argument implies that a galaxy like the Milky Way is blowing a powerful wind that carries away most of the heavy elements currently synthesized and has impacted the IGM out to at least 180 kpc. Rich clusters of galaxies appear to be closed systems in which most heavy elements are ejected from galaxies. More supernovae are required than the yield of core-collapse SNe from a Salpeter IMF. X-ray observations imply that the IGM in groups and clusters as been strongly preheated. SNe probably cannot supply the required energy, which must come from AGN.

The paper claimed that significant radial migration of stars in a stellar disk like that of the Milky Way could not occur. We now think that while the treatment of the effects of molecular clouds was correct, the paper seriously underestimated the ability of spiral arms to shift the radii of stars that corotate with them. Consequently, it is likely that significant radial migration_is_ possible.

The Milky Way's satellites provide unique information about the density of the Galactic halo at large radii. The inclusion of even a few rather inaccurate proper motions resolves an ambiguity in older mass estimates in favour of higher values. Many of the satellites are concentrated into streams. The dynamics of the Magellanic Stream provided an early indication that the halo reaches out to beyond 100 kpc. Tidal forces between the Clouds are currently disturbing the Clouds' internal dynamics. One would expect this damage to worsen rapidly as the tidal field of the MW excites the eccentricity of the Clouds' mutual orbit. This process, which has yet to be completely modelled, is important for understanding the degree of self-lensing in searches for gravitational lensing events. The Sagittarius Dwarf galaxy very likely contributes significantly to the Galactic warp. The direction of the warp's line of nodes is incorrectly predicted by the simplest models of the Dwarf's orbit. More sophisticated models, in which a complex distribution of stripped dark matter is predicted, may be more successful.

Because we know little about the Galactic force-field away from the plane, the Galactic mass distribution is very ill-determined. I show that a microlensing survey of galaxies closer than 50 Mpc would enable us to map in three dimensions the Galactic density of stellar mass, which should be strictly less than the total mass density. A lower limit can be placed on the stellar mass needed at R<R_0 to generate the measured optical depth towards sources in the bulge. If the Galaxy is barred, this limit is lower by a factor of up to two than in the axisymmetric case. Even our limited knowledge of the Galactic force field suffices to rule out the presence of the amount of mass an axisymmetric Galaxy needs to generate the measured optical depth. Several lines of argument imply that the Galaxy is strongly barred only at R < 4 kpc, and if this is the case, even barred Galaxy models cannot generate the measured optical depth without violating some constraint on the Galactic force-field. Galactic mass models that are based on the assumption that light traces mass, for which there is significant support in the inner Galaxy, yield microlensing optical depths that are smaller than the measured value by a factor of more than 2.5.

We generalize to elliptical models the argument of Kuijken (1997), which connects the microlensing optical depth towards the Galactic bulge to the Galactic rotation curve. When applied to the latest value from the MACHO collaboration for the optical depth for microlensing of bulge sources, the argument implies that the Galactic bar cannot plausibly reconcile the measured values of the optical depth, the rotation curve and the local mass density. Either there is a problem with the interpretation of the microlensing data, or our line of sight to the Galactic centre is highly atypical in that it passes through a massive structure that wraps only a small distance around the Galactic centre.

Both theory and observation indicate that galaxies like the Milky Way accrete matter at the rate of a few Msun per year.

In the generic CDM cosmogony, dark-matter halos emerge too lumpy and centrally concentrated to host observed galactic disks. Moreover, disks are predicted to be smaller than those observed. We argue that the resolution of these problems may lie with a combination of the effects of protogalactic disks, which would have had a mass comparable to that of the inner dark halo and be plausibly non-axisymmetric, and of massive galactic winds, which at early times may have carried off as many baryons as a galaxy now contains. A host of observational phenomena, from quasar absorption lines and intracluster gas through the G-dwarf problem point to the existence of such winds. Dynamical interactions will homogenize and smooth the inner halo, and the observed disk will be the relic of a massive outflow. The inner halo expanded after absorbing energy and angular momentum from the ejected material. Observed disks formed at the very end of the galaxy formation process, after the halo had been reduced to a minor contributor to the central mass budget and strong radial streaming of the gas had died down.

Possible orbital histories of the Sgr dwarf galaxy are explored. A special-purpose N-body code is used to construct the first models of the Milky Way - Sgr Dwarf system in which both the Milky Way and the Sgr Dwarf are represented by full N-body systems and followed for a Hubble time. These models are used to calibrate a semi-analytic model of the Dwarf's orbit that enable us to explore a wider parameter space than is accessible to the N-body models. We conclude that the extant data on the Dwarf are compatible with a wide range of orbital histories. At one extreme the Dwarf initially possesses 10^{11} Solar Mass and starts from a Galactocentric distance 200 kpc. At the other extreme the Dwarf starts with 10^9 Solar Mass and Galactocentric distance 60 kpc, similar to its present apocentric distance. In all cases the Dwarf is initially dark-matter dominated and the current velocity dispersion of the Dwarf's dark matter is tightly constrained to be 21 km/s. This number is probably compatible with the smaller measured dispersion of the Dwarf's stars because of (a) the dynamical difference between dark and luminous matter, and (b) velocity anisotropy.

The amount of mass contained in low-mass objects is investigated anew. Instead of using a mass-luminosity relation to convert a luminosity function to a mass function, I predict the mass-luminosity relation from assumed mass functions and the luminosity functions of Jahreiss & Wielen (1997) and Gould et al (1997). Comparison of the resulting mass-luminosity relations with data from binary stars constrains the permissible mass functions. If the mass function is assumed to be a power law, the best fitting slope lies either side of the critical slope, -2, below which the mass in low-mass objects is divergent, depending on the luminosity function adopted. If these power-law mass functions are truncated at 0.001Msun, the contribution to the local density of stars lies between 0.016 and 0.039 Msun pc^-3, in conformity with the density measured dynamically from Hipparcos stars. If the mass function is generalized from a power law to a low-order polynomial in log(M), the mass in stars with M<0.1Msun is either negligible or strongly divergent, depending on the order of the polynomial adopted.

An extension of Schwarzschild's galaxy-building technique is presented that, for the first time, enables one to build Schwarzschild models with known distribution functions (DFs). The new extension makes it possible to combine a DF that depends only on classical integrals with orbits that respect non-classical integrals. With such a combination, Schwarzschild's orbits are used only to represent the difference between the true galaxy DF and an approximating classical DF. The new method is used to construct a dynamical model of the inner Galaxy. The model is based on an orbit library that contains 22168 regular orbits. The model aims to reproduce the three-dimensional mass density of Binney, Gerhard & Spergel (1997), which was obtained through deprojection of the COBE surface photometry, and to reproduce the observed kinematics in three windows - namely Baade's Window and two off-axis fields. The model fits essentially all the available data within the innermost 3 kpc. The axis ratio and the morphology of the projected density contours of the COBE bar are recovered to good accuracy within corotation. The kinematic quantities - the line-of-sight streaming velocity and velocity dispersion, as well as the proper motions when available - are recovered, not merely for the fitted fields, but also for three new fields. The dynamical model deviates most from the input density close to the Galactic plane just outside corotation, where the deprojection of the surface photometry is suspect. The dynamical model does not reproduce the kinematics at the most distant window, where disk contamination may be severe.

N-body simulations show that when infall reorientates the outer parts of a galactic halo by several degrees per Gyr, a self-gravitating disk that is embedded in the halo develops an integral-sign warp that is comparable in amplitude to observed warps. Studies of angular-momentum acquisition suggest that the required rate of halo reorientation is realistic for galaxies like the Milky Way.

The simplest interpretation of the microlensing events towards the Large Magellanic Cloud detected by the MACHO and EROS collaborations is that about one third of the halo of our own Milky Way galaxy exists in the form of objects of around 0.5 solar mass. There are grave problems with this interpretation. A normal stellar population of 0.5 solar mass stars should be visible. The other obvious candidate for the lenses is a population of white dwarfs. But, the precursor population must have polluted the interstellar medium with metals, in conflict with current population II abundance. Here, we propose a more conventional, but at the moment more speculative, explanation. Some of the lenses are stars in the disk of the Milky Way. They lie along the line of sight to the LMC because of warping and flaring of the Galactic disk. Depending on its scalelength and ellipticity, the disk's optical depth may lie anywhere between $0.2 \times 10^{-7}$ and $0.9 \times 10^{-7}$. Together with contributions from the LMC disk and bar and perhaps even intervening stellar contaminants, the total optical depth may match the data within the uncertainties. Microlensing towards the LMC may be telling us more about the distorted structure and stellar populations of the outer Milky Way disk than the composition of the dark halo.

We summarize recent work on the structure and dynamics of the Galactic bar and inner disk. Current work focusses on constructing a quantitative model which integrates NIR photometry, source count observations, gas kinematics, stellar dynamical observations, and microlensing. Some avenues for future research are discussed.

(shortened) From a kinematically unbiased subsample of the Hipparcos catalogue we have redetermined as a function of colour the kinematics of main-sequence stars. The stars' mean heliocentric velocity nicely follows the asymmetric drift relation, except for stars blueward of B-V=0.1. Extrapolating to zero dispersion yields for the velocity of the Sun w.r.t. the LSR in km/s: U_0=10.00+/-0.36 (radially inwards), V_0=5.23+/-0.62 (in direction of galactic rotation), and W_0=7.17+/-0.38 (vertically upwards). A plot of velocity dispersion vs. colour beautifully shows Parenago's discontinuity: the dispersion is constant for B-V>0.62 and decreases towards bluer colour. We determine the velocity-dispersion tensor sigma^2_ij as function of B-V. The mixed moments involving vertical motion are zero within the errors, while sigma^2_xy is non-zero at about (10km/s)^2 independent of colour. The resulting vertex deviations are about 20 deg for early-type stars and 10+/-4 deg for old-disc stars. The persistence of the vertex deviation to late-type stars implies that the Galactic potential is significantly non-axisymmetric at the solar radius. If spiral arms are responsible for this, they cannot be tightly wound. Except for stars bluer than B-V=0.1 the ratios of the principal velocity dispersions are 2.2 : 1.4 :1, while the absolute values increase with colour from sigma_1=20km/s at B-V=0.2 to sigma_1=38km/s at Parenago's discontinuity and beyond. These ratios imply significant heating of the disc by spiral structure and that R_0/R_d=3 to 3.5, where R_d is the scale length of the disc.

A parameterized 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 scale length R_d* of the stellar disk to R0. The disk and bulge dominate v_c(R) at R<R0 only for R_d*/R0< 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 << R0 to affect profoundly the structure of the best-fitting model at R >> R0. For example, changing the disk slightly from an exponential surface-density profile significantly changes the form of v_c(R) at R >> R0, where the disk 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 halos 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 > R0 are especially important.

We show that the usual relation between redshift and angular-diameter distance can be obtained by considering light from a source to be gravitationally lensed by material that lies in the telescope beam as it passes from source to observer through an otherwise empty universe. This derivation yields an equation for the dependence of angular diameter on redshift in an inhomogeneous universe. We use this equation to model the distribution of angular-diameter distance for redshift z=3 in a realistically clustered cosmology. The distribution is such that attempts to determine q_0 from angular-diameter distances will systematically underestimate q_0 by ~0.15, and large samples would be required to beat down the intrinsic dispersion in measured values of q_0.

We study the persistence of warps in galactic discs in the presence of massive halos. A disc is approximated by a set of massive rings, while a halo is represented by a conventional n-body simulation. We confirm the conclusion of Nelson & Tremaine (1995) that a halo responds strongly to an embedded precessing disc. This response invalidates the approximations made by in the derivation of classical `modified tilt' modes. We show that the response of the halo causes the line of nodes of a disc that starts from a modified tilt mode to wind up within a few dynamical times. We explain this finding in terms of the probable spectrum of true normal modes of a combined disc-halo system.

A straightforward determination of the circular-speed curve vc(R) of the Milky Way suggests that near the Sun, vc starts to rise approximately linearly with R. If this result were correct, the Galactic mass density would have to be independent of radius at R ~> R0. We show that the apparent linear rise in v_c arises naturally if the true circular-speed curve is about constant or gently falling at R0 < R ~< 2 R0, but most tracers that appear to be at R ~> 1.25 R0 are actually concentrated into a ring of radius ~1.6 R0.

We examine the left-right asymmetry in the cleaned COBE/DIRBE near-infrared data of the inner Galaxy and show (i) that the Galactic bar is probably not seen very nearly end-on, and (ii) that even if it is, it is not highly elongated. The assumption of constant mass-to-light ratio is used to derive simulated terminal-velocity plots for the ISM from our model luminosity distributions. By comparing these plots with observed terminal velocities we determine the mass-to-light ratio of the near-IR bulge and disk. Assuming that all this mass contributes to gravitational microlensing we compute optical depths $\tau$ for microlensing in Galactic-centre fields. For three models with bar major axis between $10\deg-25\deg$ from the Sun-Galactic Center line, the resulting optical depths in Baade's window lie in the range $0.83\times10^{-6} \lta \tau \lta 0.89\times10^{-6}$ for main-sequence stars and $1.2\times10^{-6} \lta \tau \lta 1.3\times10^{-6}$ for red-clump giants. We discuss a number of uncertainties including possible variations of the near-infrared mass-to-light ratio. We conclude that, although the values predicted from analyzing the COBE and gas velocity data are inconsistent at the $2-2.5\sigma$ level with recent observational determinations of $\tau$, we believe they should be taken seriously.

The light distribution in the inner few kiloparsecs of the Milky Way is recovered non-parametrically from a dust-corrected near-infrared COBE/DIRBE surface brightness map of the inner Galaxy. The best fits to the photometry are obtained when the Sun is assumed to lie $\sim14\pm4\pc$ below the plane. The recovered density distributions clearly show an elongated three-dimensional bulge set in a highly non-axisymmetric disk. In the favoured models, the bulge has axis ratios $1{:}0.6{:}0.4$ and semi-major axis length $\sim2\kpc$. Its nearer long axis lies in the first quadrant. The bulge is surrounded by an elliptical disk that extends to $\sim2\kpc$ on the minor axis and $\sim3.5\kpc$ on the major axis. In all models there is a local density minimum $\sim2.2\kpc$ down the minor axis. The subsequent maximum $\sim3\kpc$ down the minor axis (corresponding to $l\simeq-22\deg$ and $l\simeq 17\deg$) may be associated with the Lagrange point L$_4$. From this identification and the length of the bulge-bar, we infer a pattern speed $\Omega_b\simeq 60-70\kms\kpc^{-1}$ for the bar. Experiments in which pseudo-data derived from models with spiral structure were deprojected under the assumption that the Galaxy is either eight-fold or four-fold symmetric, indicate that the highly non-axisymmetric disks recovered from the COBE data could reflect spiral structure within the Milky Way if that structure involves density contrasts greater than $\gta 3$ at NIR wavelengths. These experiments indicate that the angle $\phi_0$ between the Sun--centre line and a major axis of the bulge lies near $20\deg$.

The only way to map the Galaxy's gravitational potential $\Phi({\bf x})$ and the distribution of matter that produces it is by modelling the dynamics of stars and gas. Observations of the kinematics of gas provide key information about gradients of $\Phi$ within the plane, but little information about the structure of $\Phi$ out of the plane. Traditional Galaxy models {\em assume}, for each of the Galaxy's components, arbitrary flattenings, which together with the components' relative masses yield the model's equipotentials. However, the Galaxy's isopotential surfaces should be {\em determined\/} directly from the motions of stars that move far from the plane. Moreover, from the kinematics of samples of such stars that have well defined selection criteria, one should be able not only to map $\Phi$ at all positions, but to determine the distribution function $f_i({\bf x},{\bf v})$ of each stellar population $i$ studied. These distribution functions will contain a wealth of information relevant to the formation and evolution of the Galaxy. An approach to fitting a wide class of dynamical models to the very heterogeneous body of available data is described and illustrated.

An algorithm is developed and tested for the problem posed by photometric observations of the bulge of the Milky Way. The latter subtends a non-trivial solid angle on the sky, and we show that this permits inversion of the projected brightness distribution under the assumption that the bulge has three orthogonal mirror planes of specified orientation. A serious error in the assumed orientation of the mirror planes should be detectable.

Axisymmetric density distributions are constructed which are invisible when viewed from a range of inclination angles $i$. By adding such distributions to a model galaxy, it can be made either disky or boxy without in any way affecting its projected image. As the inclination of a galaxy decreases from edge-on to face on, the range of `invisible' densities, the uncertainty in the deprojection, and the sensitivity of the deprojection to noise all increase. The relation between these phenomena is clarified by an analysis of Palmer's deprojection algorithm. These results imply that disk-to-bulge ratios are in principle ill-determined from photometry unless the disk is strong or the system is seen precisely edge-on. The uncertain role of third integrals in galaxies makes it unclear to what degree this indeterminacy can be resolved by kinematic studies.

A scheme is developed which enables one to trace backwards in time the cosmic density and velocity fields, and to determine accurately the current-epoch velocity field from the current-epoch density field, or vice versa. The scheme implements the idea of Giavalisco \etal\ (1993) that the principle of least action should be used to formulate gravitational instability as a two-point boundary-value problem. We argue that the Eulerian formulation of the problem is to be preferred to the Lagrangian one, on grounds of computational simplicity, of ease of interfacing with observational data, and of internal consistency at early times. The scheme is successfully tested on an exact solution in one dimension, and on currently Gaussian fields in one and two dimensions. The application of the scheme to real observational data appears to be eminently feasible, though computationally costly.

This paper contains all known errata in the book "Galactic Dynamics".