We study the late-time evolution of the central regions of two Milky Way-like simulations of galaxies formed in a cosmological context, one hosting a fast bar and the other a slow one. We find that bar length, R_b, measurements fluctuate on a dynamical timescale by up to 100%, depending on the spiral structure strength and measurement threshold. The bar amplitude oscillates by about 15%, correlating with R_b. The Tremaine-Weinberg-method estimates of the bars' instantaneous pattern speeds show variations around the mean of up to ~20%, typically anti-correlating with the bar length and strength. Through power spectrum analyses, we establish that these bar pulsations, with a period in the range ~60-200 Myr, result from its interaction with multiple spiral modes, which are coupled with the bar. Because of the presence of odd spiral modes, the two bar halves typically do not connect at exactly the same time to a spiral arm, and their individual lengths can be significantly offset. We estimated that in about 50% of bar measurements in Milky Way-mass external galaxies, the bar lengths of SBab type galaxies are overestimated by ~15% and those of SBbc types by ~55%. Consequently, bars longer than their corotation radius reported in the literature, dubbed "ultra-fast bars", may simply correspond to the largest biases. Given that the Scutum-Centaurus arm is likely connected to the near half of the Milky Way bar, recent direct measurements may be overestimating its length by 1-1.5 kpc, while its present pattern speed may be 5-10 km/s/kpc smaller than its time-averaged value.

We use the second Gaia data release in combination with the catalog of Sanders & Das (2018) to dissect the Milky Way disc in phase-space and relative ages. We confirm and report the existence of multiple velocity moving groups at low azimuthal velocities and angular momenta, below Arcturus, regularly separated by $\sim~18-20\,\rm{km s^{-1}}$ in azimuthal velocity. Such features were predicted to exist more than ten years ago from the reaction of the Milky Way to a perturbation in the disc undergoing phase-mixing. These structures appear slightly younger than their phase-space surroundings, arguing against an extra-galactic origin. We also identify in relative age many of the classical ridges in the plane of azimuthal velocity vs. Galactocentric radius, traditionally associated with resonance features. These ridges are also younger than their phase-space surroundings in accordance with predictions from recent state-of-the-art cosmological hydrodynamical simulations of Milky Way-like galaxies. We study the response of dynamically young and old stellar disc populations to resonances from an analytic model of a large bar, which, remarkably, qualitatively reproduces the trends seen in the data. Our results re-inforce the idea that the Galactic disc is currently being shaped by both internal and external perturbations, and that, while absolute isochrone ages have to be taken with great care, exploring the dynamical structure of the disc in stellar ages, especially with future asteroseismic data, will provide much stronger constraints than metallicity/abundance trends alone.

We investigate the disc-halo connection in massive (Mstar/Msun>5e10) disc galaxies from the cosmological hydrodynamical simulations EAGLE and IllustrisTNG, and compare it with that inferred from the study of HI rotation curves in nearby massive spirals from the Spitzer Photometry and Accurate Rotation Curves (SPARC) dataset. We find that discrepancies between the the simulated and observed discs arise both on global and on local scales. Globally, the simulated discs inhabit halos that are a factor ~4 (in EAGLE) and ~2 (in IllustrisTNG) more massive than those derived from the rotation curve analysis of the observed dataset. We also use synthetic rotation curves of the simulated discs to demonstrate that the recovery of the halo masses from rotation curves are not systematically biased. We find that the simulations predict dark-matter dominated systems with stellar-to-total enclosed mass ratios that are a factor of 1.5-2 smaller than real galaxies at all radii. This is an alternative manifestation of the `failed feedback problem', since it indicates that simulated halos hosting massive discs have been too inefficient at converting their baryons into stars, possibly due to an overly efficient stellar and/or AGN feedback implementation.

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.

Double-lined spectroscopic binaries (SB2s) are one of the main sources of stellar masses, as additional observations are only needed to give the inclinations of the orbital planes in order to obtain the individual masses of the components. For this reason, we are observing a selection of SB2s using the SOPHIE spectrograph at the Haute-Provence observatory in order to precisely determine their orbital elements. Our objective is to finally obtain masses with an accuracy of the order of one percent by combining our radial velocity (RV) measurements and the astrometric measurements that will come from the Gaia satellite. We present here the RVs and the re-determined orbits of 10 SB2s. In order to verify the masses we will derive from Gaia, we obtained interferometric measurements of the ESO VLTI for one of these SB2s. Adding the interferometric or speckle measurements already published by us or by others for 4 other stars, we finally obtain the masses of the components of 5 binary stars, with masses ranging from 0.51 to 2.2 solar masses, including main-sequence dwarfs and some more evolved stars whose location in the HR diagram has been estimated.

We reconsider the case for the association of Galactic globular clusters (GCs) to the tidal stream of the Sagittarius dwarf spheroidal galaxy (Sgr dSph), using Gaia DR2 data. We use RR Lyrae to trace the stream in 6D and we select clusters matching the observed stream in position and velocity. In addition to the clusters residing in the main body of the galaxy (M 54, Ter 8, Ter 7, Arp 2) we confirm the membership of Pal 12 and Whiting 1 to the portion of the trailing arm populated by stars lost during recent perigalactic passages. NGC 2419, NGC 5634 and NGC 4147 are very interesting candidates, possibly associated to more ancient wraps of the stream. We note that all these clusters, with the exception of M 54, that lies within the stellar nucleus of the galaxy, are found in the trailing arm of the stream. The selected clusters are fully consistent with the [Fe/H] vs. [Mg/Fe], [Ca/Fe] patterns and the age-metallicity relation displayed by field stars in the main body of Sgr dSph.

We present the first full six-dimensional panoramic portrait of the Sagittarius stream, obtained by searching for wide stellar streams in the Gaia DR2 dataset with the STREAMFINDER algorithm. We use the kinematic behavior of the sample to devise a selection of Gaia RR Lyrae, providing excellent distance measurements along the stream. The proper motion data are complemented with radial velocities from public surveys. We find that the global morphological and kinematic properties of the Sagittarius stream are still reasonably well reproduced by the simple Law & Majewski (2010) model (LM10), although the model overestimates the leading arm and trailing arm distances by up to $\sim 15$%. The sample newly reveals the leading arm of the Sagittarius stream as it passes into very crowded regions of the Galactic disk towards the Galactic Anticenter direction. Fortuitously, this part of the stream is almost exactly at the diametrically opposite location from the Galactic Center to the progenitor, which should allow an assessment of the influence of dynamical friction and self-gravity in a way that is nearly independent of the underlying Galactic potential model.

The density variations in thin stellar streams may encode important information on the nature of the dark matter. For instance, if dark matter aggregates into massive sub-halos, these perturbers are expected to scatter stars out of dynamically cold stellar streams, possibly leading to detectable gaps in those structures. Here we re-examine the density variations in the GD-1 stream, using Gaia DR2 astrometry, Pan-STARRS photometry, together with high precision radial velocities measured with the CFHT/ESPaDOnS and VLT/UVES instruments and complemented with public radial velocity catalogs. We show that after correcting for projection effects, the density profile exhibits high contrast periodic peaks, separated by $2.64\pm0.18$ kpc. An N-body simulation is presented that reproduces this striking morphology with simple epicyclic motion in a smooth Galactic potential. We also discuss the reliability of measuring density variations using ground-based photometric surveys, and for the particular case of GD-1 we highlight some of the artifacts present in the Gaia DR2 catalog along its track. Massive dark subhalos do not appear to be required to explain the density clumping along GD-1.

The formation and evolution of galaxies is highly dependent on the dynamics of stars and gas, which is governed by the underlying law of gravity. To investigate how the formation and evolution of galaxies takes place in Milgromian gravity (MOND), we present full hydrodynamical simulations with the Phantom of Ramses (POR) code. These are the first-ever galaxy formation simulations done in MOND with detailed hydrodynamics, including star formation, stellar feedback, radiative transfer and supernovae. These models start from simplified initial conditions, in the form of isolated, rotating gas spheres in the early Universe. These collapse and form late-type galaxies obeying several scaling relations, which was not a priori expected. The formed galaxies have a compact bulge and a disk with exponentially decreasing surface mass density profiles and scale lengths consistent with observed galaxies, and vertical stellar mass distributions with distinct exponential profiles (thin and thick disk). This work thus shows for the first time that disk galaxies with exponential profiles in both gas and stars are a generic outcome of collapsing gas clouds in MOND. These models have a slight lack of stellar angular momentum because of their somewhat compact stellar bulge, which is connected to the simple initial conditions and the negligible later gas accretion. We also analyse how the addition of more complex baryonic physics changes the main resulting properties of the models and find this to be negligibly so in the Milgromian framework.

The empirical scaling relations observed in disk galaxies remain challenging for models of galaxy formation. The most striking among these is the Mass Discrepancy-Acceleration Relation (MDAR), which encodes both a tight baryonic Tully-Fisher relation (BTFR) and the observed diversity of galaxy rotation curves through the central surface density relation (CSDR). Building on our earlier work, we propose here that the MDAR is the result of interactions between baryons and 'Baryon-Interacting Dark Matter' (BIDM), which heat up the dark matter. Following a bottom-up, hydrodynamical approach, we find that the MDAR follows if: $i)$ the BIDM equation of state approximates that of an ideal gas; $ii)$ the BIDM relaxation time is order the Jeans time; $iii)$ the heating rate is inversely proportional to the BIDM density. Remarkably, under these assumptions the set of hydrodynamical equations together with Poisson's equation enjoy an anisotropic scaling symmetry. In the BIDM-dominated regime, this gives rise to an enhanced symmetry which fully captures the low-acceleration limit of the MDAR. We then show that, assuming a cored pseudo-isothermal profile at equilibrium, this set of equations gives rise to parameters reproducing the MDAR. Specifically, in the flat part of the rotation curve the asymptotic rotational velocity matches the parametric dependence of the BTFR. Moreover, in the central region of high-surface brightness galaxies, the profile reproduces the CSDR. Finally, by studying the time-dependent approach to equilibrium, we derive a global combination of the BTFR and CSDR, which matches the expectations in low surface-brightness galaxies. The form of the heating rate also makes model-independent predictions for various cosmological observables. We argue that our scenario satisfies existing observational constraints, and, intriguingly, offers a possible explanation to the EDGES anomaly.

The outskirts of globular clusters (GCs) simultaneously retain crucial information about their formation mechanism and the properties of their host galaxy. Thanks to the advent of precision astrometry both their morphological and kinematic properties are now accessible. Here we present the first dynamical study of the outskirts of the retrograde GC NGC 3201 until twice its Jacobi radius (< 100 pc), using specifically-selected high-quality astrometric data from Gaia DR2. We report the discovery of a stellar overdensity along the South-East/North-West direction that we identify as tidal tails. The GC is characterized globally by radial anisotropy and a hint of isotropy in the outer parts, with an excess of tangential orbits around the lobes corresponding to the tidal tails, in qualitative agreement with an N-body simulation. Moreover, we measure flat velocity dispersion profiles, reaching values of $3.5\pm0.9$ km/s until beyond the Jacobi radius. While tidal tails could contribute to such a flattening, this high velocity dispersion value is in disagreement with the expectation from the sole presence of potential escapers. To explain this puzzling observation, we discuss the possibility of an accreted origin of the GC, the presence of a dark matter halo --leftover of its formation at high redshift -- and the possible effects of non-Newtonian dynamics. Our study uncovers a new path for the study of GC formation and of the properties of the Milky Way potential in the era of precision astrometry.

We report the discovery of a 70-degree long stellar stream in the Milky Way halo, which criss-crosses the well known "GD-1" stream. We show that this new stellar structure ("Kshir") and GD-1 lie at similar distance, and are remarkably correlated in kinematics. We propose several explanations for the nature of this new structure and its possible association with GD-1. However, a scenario in which these two streams were accreted onto the Milky Way within the same dark matter sub-halo seems to provide a natural explanation for their phase-space entanglement, and other complexities of this coupled-system.

We present radial velocities for five member stars of the recently discovered young (age= 100-150 Myr) stellar system Price-Whelan 1 (PW 1), that is located far away in the Galactic Halo (D~ 29 kpc, Z~ 15 kpc), and that is probably associated to the Leading Arm (LA) of the Magellanic Stream. We measure the systemic radial velocity of PW 1, V_r=275 +/- 10 km/s, significantly larger than the velocity of the LA gas in the same direction. We re-discuss the main properties and the origin of this system in the light of these new observations, computing the orbit of the system and comparing its velocity with that of the HI in its surroundings. We show that the bulk of the gas at the velocity of the stars is more than 10 deg (5~kpc) away from PW 1 and the velocity difference between the gas and the stars become larger as gas closer to the stars is considered. We discuss the possibilities that (a) the parent gas cloud was dissolved by the interaction with the Galactic gas, and (b) that the parent cloud is the high velocity cloud HVC287.5+22.5+240, lagging behind the stellar system by ~25 km/s and ~10 deg.

We present a new fully data-driven algorithm that uses photometric data from the Canada-France-Imaging-Survey (CFIS; $u$), Pan-STARRS 1 (PS1; $griz$), and Gaia ($G$) to discriminate between dwarf and giant stars and to estimate their distances and metallicities. The algorithm is trained and tested using the SDSS/SEGUE spectroscopic dataset and Gaia photometric/astrometric dataset. At [Fe/H]$<-1.2$, the algorithm succeeds in identifying more than 70% of the giants in the training/test set, with a dwarf contamination fraction below 30% (with respect to the SDSS/SEGUE dataset). The photometric metallicity estimates have uncertainties better than 0.2 dex when compared with the spectroscopic measurements. The distances estimated by the algorithm are valid out to a distance of at least $\sim 80$ kpc without requiring any prior on the stellar distribution, and have fully independent uncertainities that take into account both random and systematic errors. These advances allow us to estimate these stellar parameters for approximately 12 million stars in the photometric dataset. This will enable studies involving the chemical mapping of the distant outer disc and the stellar halo, including their kinematics using the Gaia proper motions. This type of algorithm can be applied in the Southern hemisphere to the first release of LSST data, thus providing an almost complete view of the external components of our Galaxy out to at least $\sim 80$ kpc. Critical to the success of these efforts will be ensuring well-defined spectroscopic training sets that sample a broad range of stellar parameters with minimal biases. A catalogue containing the training/test set and all relevant parameters within the public footprint of CFIS is available online.

The second data release of the Gaia mission has revealed, in stellar velocity and action space, multiple ridges, the exact origin of which is still debated. Recently, we demonstrated that a large Galactic bar with pattern speed 39 km/s/kpc does create most of the observed ridges. Among those ridges, the Hercules moving group would then be associated to orbits trapped at the co-rotation resonance of the bar. Here we show that a distinctive prediction of such a model is that the angular momentum of Hercules at the Sun's radius must significantly decrease with increasing Galactocentric azimuth, i.e. when getting closer to the major axis of the bar. We show that such a dependence of the angular momentum of trapped orbits on the azimuth would on the other hand not happen close to the outer Lindblad resonance of a faster bar, unless the orbital distribution is still far from phase-mixed, namely for a bar perturbation younger than ~ 2 Gyr. Using Gaia DR2 and Bayesian distances from the StarHorse code, and tracing the average Galactocentric radial velocity as a function of angular momentum and azimuth, we show that the Hercules angular momentum changes significantly with azimuth as expected for the co-rotation resonance of a dynamically old large bar.

In $\Lambda$CDM cosmology, to first order, galaxies form out of the cooling of baryons within the virial radius of their dark matter halo. The fractions of mass and angular momentum retained in the baryonic and stellar components of disc galaxies put strong constraints on our understanding of galaxy formation. In this work, we derive the fraction of angular momentum retained in the stellar component of spirals, $f_j$, the global star formation efficiency $f_M$, and the ratio of the asymptotic circular velocity ($V_{\rm flat}$) to the virial velocity $f_V$, and their scatter, by fitting simultaneously the observed stellar mass-velocity (Tully-Fisher), size-mass, and mass-angular momentum (Fall) relations. We compare the goodness of fit of three models: (i) where the logarithm of $f_j$, $f_M$, and $f_V$ vary linearly with the logarithm of the observable $V_{\rm flat}$; (ii) where these values vary as a double power law; and (iii) where these values also vary as a double power law but with a prior imposed on $f_M$ such that it follows the expectations from widely used abundance matching models. We conclude that the scatter in these fractions is particularly small ($\sim 0.07$ dex) and that the linear model is by far statistically preferred to that with abundance matching priors. This indicates that the fundamental galaxy formation parameters are small-scatter single-slope monotonic functions of mass, instead of being complicated non-monotonic functions. This incidentally confirms that the most massive spiral galaxies should have turned nearly all the baryons associated with their haloes into stars. We call this the failed feedback problem.

The period 2035-50 considered in the ESA Voyage long-term plan will coincide with a series of foreseeable advances in the characterization of the stellar content of the Milky Way. The Gaia mission, combined with large-scale spectroscopic surveys, is helping to build an unprecedented census in terms of the astrometric, kinematic and chemical properties of Galactic stellar populations. Within a decade, precise measurements of such properties will be available for hundreds of millions of stars. Meanwhile, time-domain surveys initiated with CoRoT and Kepler/K2 and carried on by space missions such as TESS and PLATO or ground-based projects like the LSST, will have brought asteroseismology to a high level of maturity. The combination of precise ages from asteroseismology with astrometric and spectroscopic data, on large stellar samples, is allowing Galactic archaeologists to gain new insight into the assembly history of the Milky Way. Recent breakthroughs --- based on the detection of solar-like oscillations in tens of thousands of red-giant stars --- demonstrate the potential of such approach. Therefore, we are convinced that an all-sky, high-cadence, long-duration stellar variability survey will become a scientific priority in the 2035-50 period. The Chronos concept presented here consists in a time-domain extension to Gaia. It will allow for mass and age estimates for half a million red giants within 1.7 kpc from the Sun and hence shed a new light on our understanding of the Galactic dynamics and archaeology. In terms of the targeted pulsators, Chronos will bridge the gap between PLATO and the LSST by surveying stars all the way from the subgiant branch to the early AGB. Finally, it will surpass all previous surveys capable of conducting asteroseismology in terms of the combined sky coverage and duration of the observations (2 x 3.75 months over the whole sky and >5 years in the CVZ).

The missing mass problem has not been solved decisively yet. Observations show that if gravity is to be modified, then the MOND theory is its excellent approximation on galactic scales. MOND suggests an adjustments of the laws of physics in the limit of low accelerations. Comparative simulations of interacting galaxies in MOND and Newtonian gravity with dark matter revealed two principal differences: 1) galaxies can have close flybys without ending in mergers in MOND because of weaker dynamical friction, and 2) tidal dwarf galaxies form very easily in MOND. When this is combined with the fact that many interacting galaxies are observed at high redshift, we obtain a new perspective on tidal features: they are often formed by non-merging encounters and tidal disruptions of tidal dwarf galaxies. Here we present the results from our self-consistent MOND $N$-body simulation of a close flyby of two galaxies similar to the Milky Way. It turns out that most types of the structures that are traditionally assigned to galaxy mergers can be formed by non-merging encounters, including tidal arms, bridges, streams, shells, disk warps, thick disks, and most probably also disks of satellites. The success of MOND in explaining the dynamics of galaxies hints us that this way of formation of tidal structures should be considered seriously.

Recent studies have shown that the passage of a massive satellite through the disk of a spiral galaxy can induce vertical wobbles in the disk and produce features such as in-plane rings and phase-space spirals. Here we analyze a high-resolution N-body simulation of a live stellar disk perturbed by the recent passage of a massive dwarf galaxy that induces such perturbations. We study the implications of the phase-space structures for the estimate of the matter density through traditional Jeans modelling. The dwarf satellite excites rapid time-variations in the potential, leading to a significant bias of the local matter surface density determined through such a method. In particular, while the Jeans modelling gives reasonable estimates in the most overdense regions of the disk, we show that it tends to overestimate the dynamical surface density in underdense regions. In these regions, the phase-space spiral is indeed more marked in the surface of section of height vs. vertical velocity. This prediction can be verified with future Gaia data releases. Our finding is highly relevant for future attempts at determining the dynamical surface density of the outer Milky Way disk as a function of radius. The outer disk of the Milky Way is indeed heavily perturbed, and \textit{Gaia} DR2 data have clearly shown that such phase-space perturbations are even present locally. While our results show that traditional Jeans modelling should give reliable results in overdense regions of the disk, the important biases in underdense regions call for the development of non-equilibrium methods to estimate the dynamical matter density locally and in the outer disk.

Using Jeans modelling, we calculate the velocity dispersion profile of the ultra-diffuse galaxy (UDG) Dragonfly 44 in MOND. For the nominal mass-to-light ratio from the literature and an isotropic profile, the agreement with the data is excellent near the center of the galaxy. However, in modified gravity, close to the cluster core, the gravitational environment should bring the galaxy back towards the Newtonian behavior. This success of the isolated MOND prediction for the central velocity dispersion could then mean that the galaxy is at large distances ($\gg5$ Mpc) from the cluster core, as hinted by the fact that close-by UDGs share similar velocities with a dispersion well below that of the cluster itself. There is however a two-sigma tension in the outer part of the UDG due to an increase of the observed dispersion profile with respect to the flat MOND prediction. This deviation could simply be a measurement error. Other possibilities could be, for a UDG far from the cluster, a higher-than-nominal baryonic mass with a tangentially anisotropic dispersion profile, or even a sign of a dark baryonic halo. If the UDG is closer to the cluster core, the deviation could be a sign that it is in the process of disruption.

The second data release of the Gaia mission has revealed a very rich structure in local velocity space. In terms of in-plane motions, this rich structure is also seen as multiple ridges in the actions of the axisymmetric background potential of the Galaxy. These ridges are probably related to a combination of effects from ongoing phase-mixing and resonances from the spiral arms and the bar. We have recently developed a method to capture the behaviour of the stellar phase-space distribution function at a resonance, by re-expressing it in terms of a new set of canonical actions and angles variables valid in the resonant region. Here, by properly treating the distribution function at resonances, and by using a realistic model for a slowly rotating large Galactic bar with pattern speed 39 km/s/kpc, we show that no less than six ridges in local action space can be related to resonances with the bar. Two of these at low angular momentum correspond to the corotation resonance, and can be associated to the Hercules moving group in local velocity space. Another one at high angular momentum corresponds to the outer Lindblad resonance, and can tentatively be associated to the velocity structure seen as an arch at high azimuthal velocities in Gaia data. The other ridges are associated to the 3:1, 4:1 and 6:1 resonances. The latter can be associated to the so-called 'horn' of the local velocity distribution. While it is clear that effects from spiral arms and incomplete phase-mixing related to external perturbations also play a role in shaping the complex kinematics revealed by Gaia data, the present work demonstrates that, contrary to common misconceptions, the bar alone can create multiple prominent ridges in velocity and action space.

The Milky Way, the Andromeda galaxy, and Centaurus A host flattened distributions of satellite galaxies which exhibits coherent velocity trends indicative of rotation. Comparably extreme satellite structures are very rare in cosmological LCDM simulations, giving rise to the `satellite plane problem'. As a possible explanation it has been suggested that earlier-forming, higher concentration host halos contain more flattened and kinematically coherent satellite planes. We have tested for such a proposed correlation between the satellite plane and host halo properties in the ELVIS suite of simulations. We find evidence neither for a correlation of plane flattening with halo concentration or formation time, nor for a correlation of kinematic coherence with concentration. The height of the thinnest sub-halo planes does correlate with the host virial radius and with the radial extent of the sub-halo system. This can be understood as an effect of not accounting for differences in the radial distribution of sub-halos, and selecting them from different volumes than covered by the actual observations. Being part of a halo pair like the Local Group does not result in more narrow or more correlated satellite planes either. Additionally, using the PhatELVIS simulations we show that the presence of a central galaxy potential does not favor more narrow or more correlated satellite planes, it rather leads to slightly wider planes. Such a central potential is a good approximation of the dominant effect baryonic physics in cosmological simulations has on a sub-halo population. This suggests that, in contrast to other small-scale problems, the planes of satellite galaxies issue is made worse by accounting for baryonic effects.

Dark matter-baryon scaling relations in galaxies are important in order to constrain galaxy formation models. Here, we provide a modern quantitative assessment of those relations, by modelling the rotation curves of galaxies from the Spitzer Photometry and Accurate Rotation Curves (SPARC) database with the Einasto dark halo model. We focus in particular on the comparison between the original SPARC parameters, with constant mass-to-light ratios for bulges and disks, and the parameters for which galaxies follow the tightest radial acceleration relation. We show that fits are improved in the second case, and that the pure halo scaling relations also become tighter. We report that the density at the radius where the slope is -2 is strongly anticorrelated to this radius, and to the Einasto index. The latter is close to unity for a large number of galaxies, indicative of large cores. In terms of dark matter-baryon scalings, we focus on relations between the core properties and the extent of the baryonic component, which are relevant to the cusp-core transformation process. We report a positive correlation between the core size of halos with small Einasto index and the stellar disk scale-length, as well as between the averaged dark matter density within 2 kpc and the baryon-induced rotational velocity at that radius. This finding is related to the consequence of the radial acceleration relation on the diversity of rotation curve shapes, quantified by the rotational velocity at 2 kpc. While a tight radial acceleration relation slightly decreases the observed diversity compared to the original SPARC parameters, the diversity of baryon-induced accelerations at 2 kpc is sufficient to induce a large diversity, incompatible with current hydrodynamical simulations of galaxy formation, while maintaining a tight radial acceleration relation.

The possibility that ultra-diffuse galaxies lacking dark matter has recently stimulated interest to check the validity of Modified Newton Dynamics (MOND) predictions on the scale of such galaxies. It has been shown that the External Field Effect (EFE) induced by the close-by galaxy can suppress the velocity dispersion of these systems, so that they appear almost dark matter free in the Newtonian context. Here, following up on this, we are making a priori predictions for the velocity dispersion of 22 ultra-diffuse galaxies in the nearby Universe. This sample can be used to test MOND and the EFE with future follow-up measurements. We construct a catalog of nearby ultra-diffuse galaxies in galaxy group environments, and set upper and lower limits for the possible velocity dispersion allowed in MOND, taking into account possible variations in the mass-to-light ratio of the dwarf and in the distance to the galaxy group. The prediction for the velocity dispersion is made as a function of the three dimensional separation of the dwarf to its host. In 17 out of 22 cases, the EFE plays a crucial role in the prediction.

We present the kinematics of Blue Straggler (BS) stars identified in the Canada-France-Imaging-Survey (CFIS), covering 4000 deg$^2$ on the sky in the $u$-band. The BS sample, characterised through CFIS and Pan-STARRS photometry, has been kinematically decomposed into putative halo and disc populations after cross-matching with Gaia astrometry and SDSS/SEGUE/LAMOST spectroscopy. This decomposition clearly reveals the strong flaring of the outer Milky Way disc. In particular, we show that we can detect this flaring up to a vertical height of $|Z| \simeq 8$ kpc at a Galactocentric distance of $R\sim $27 kpc. While some small level of flaring is expected for extended discs built up by radial migration, we demonstrate that this observed very strong flaring of the Milky Way disc is more likely consistent with it being dynamically heated by the repeated passage of the Sagittarius dwarf spheroidal galaxy through the midplane.

We present the stellar density profile of the outer halo of the Galaxy traced over a range of Galactocentric radii from $15< R_{GC} < 220$ kpc by blue horizontal branch (BHB) stars. These stars are identified photometrically using deep $u-$band imaging from the new Canada-France-Imaging-Survey (CFIS) that reaches 24.5 mag. This is combined with $griz$ bands from Pan-STARRS 1 and covers a total of $\sim4000$ deg$^2$ of the northern sky. We present a new method to select BHB stars that has low contamination from blue stragglers and high completeness. We use this sample to measure and parameterize the three dimensional density profile of the outer stellar halo. We fit the profile using (i) a simple power-law with a constant flattening (ii) a flattening that varies as a function of Galactocentric radius (iii) a broken power law profile. We find that outer stellar halo traced by the BHB is well modelled by a broken power law with a constant flattening of $q=0.86 \pm 0.02$, with an inner slope of $\gamma=4.24 \pm 0.08$. This is much steeper than the preferred outer profile that has a slope of $\beta=3.21\pm 0.07$ after a break radius of $r_b=41.4^{+2.5}_{-2.4}$ kpc. The outer profile of the stellar halo trace by BHB stars is shallower than that recently measured using RR Lyrae, a surprising result given the broad similarity of the ages of these stellar populations.

We measure the escape speed curve of the Milky Way based on the analysis of the velocity distribution of $\sim 2850$ counter-rotating halo stars from the Gaia DR2. The distances were estimated through the StarHorse code, and only stars with distance errors smaller than 10 per cent were used in the study. The escape speed curve is measured at Galactocentric radii ranging from $\sim 5$ kpc to $\sim 10.5$ kpc. The local Galactic escape at the Sun's position is estimated to be $v_\mathrm{e}(r_\odot)=580 \pm 63~\mathrm{km~s^{-1}}$, and it rises towards the Galactic center. Defined as the minimum speed required to reach three virial radii, our estimate of the escape speed as a function of radius implies, for a Navarro-Frenk-White profile and local circular velocity of $240~\mathrm{km~s^{-1}}$, a dark matter mass $M_{200}=1.28^{+0.68}_{-0.50} \times 10^{12}~M_\odot$ and a high concentration $c_{200}=11.09^{+2.94}_{-1.79}$. Assuming the mass-concentration relation of $\Lambda$CDM, we get $M_{200}=1.55_{-0.51}^{+0.64}\times 10^{12}~M_\odot$, $c_{200}=7.93_{-0.27}^{+0.33}$, for a local circular velocity of $228~\mathrm{km~s^{-1}}$.

Before the publication of the Gaia DR2 we confirmed with RAVE and TGAS an observation recently made with the GALAH survey by Quillen ey al. concerning the Coma Berenices moving group in the Solar neighbourhood, namely that it is only present at negative Galactic latitudes. This allowed us to show that it is coherent in vertical velocity, providing a first evidence for incomplete vertical phase-mixing. We estimated for the first time from dynamical arguments that the moving group must have formed at most ~ 1.5 Gyr ago, and related this to a pericentric passage of the Sagittarius dwarf satellite galaxy. The present note is a rewritten version of the original arXiv post on this result now also including a confirmation of our finding with Gaia DR2.

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.

The observed tightness of the mass discrepancy-acceleration relation (MDAR) poses a fine-tuning challenge to current models of galaxy formation. We propose that this relation could arise from collisional interactions between baryons and dark matter (DM) particles, without the need for modification of gravity or ad hoc feedback processes. We assume that these interactions satisfy the following three conditions: (i) the relaxation time of DM particles is comparable to the dynamical time in disk galaxies; (ii) DM exchanges energy with baryons due to elastic collisions; (iii) the product between the baryon-DM cross section and the typical energy exchanged in a collision is inversely proportional to the DM number density. We present an example of a particle physics model that gives a DM-baryon cross section with the desired density and velocity dependence. Direct detection constraints require our DM particles to be either very light ($m << m_b$) or very heavy ($m >> m_b$), corresponding respectively to heating and cooling of DM by baryons. In both cases, our mechanism applies and an equilibrium configuration can in principle be reached. Here, we focus on the heavy DM/cooling case as it is technically simpler. Under these assumptions, we find that rotationally-supported disk galaxies could naturally settle to equilibrium configurations satisfying a MDAR at all radii without invoking finely tuned feedback processes. We also discuss issues related to the small scale clumpiness of baryons, as well as predictions for pressure-supported systems. We argue in particular that galaxy clusters do not follow the MDAR despite being DM-dominated because they have not reached their equilibrium configuration. Finally, we revisit existing phenomenological, astrophysical and cosmological constraints on baryon-DM interactions in light of the unusual density dependence of the cross section.

The dwarf galaxy NGC1052-DF2 has recently been identified as potentially lacking dark matter. If correct, this could be a challenge for MOND, which predicts that low surface brightness galaxies should evince large mass discrepancies. However, the correct prediction of MOND depends on both the internal field of the dwarf and the external field caused by its proximity to the giant elliptical NGC1052. Taking both into consideration under plausible assumptions, we find $\sigma_{\rm MOND} = 13.4^{+4.8}_{-3.7}\;\mathrm{km}\,\mathrm{s}^{-1}$. This is only marginally higher than the claimed 90\% upper limit on the velocity dispersion ($\sigma < 10.5\;\mathrm{km}\,\mathrm{s}^{-1}$), and compares well with the observed root mean square velocity dispersion ($\sigma = 14.3\;\mathrm{km}\,\mathrm{s}^{-1}$). We also discuss a few caveats on both the observational and theoretical side. On the theory side, the internal virialization time in this dwarf may be longer that the time scale of variation of the external field. On the observational side, the paucity of data and their large uncertainties call for further analysis of the velocity dispersion of NGC1052-DF2, to check whether it poses a challenge to MOND or is a success thereof.

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.

(abridged) The Milky Way (MW) and Andromeda (M31) galaxies possess rotating planes of satellites. Their formation has not been explained satisfactorily yet. It was suggested that the MW and M31 satellites are ancient tidal dwarf galaxies, which could explain their configuration. This suggestion gained support by an analytic backward calculation of the relative MW-M31 orbit in the MOND modified dynamics paradigm by Zhao et al. (2013) implying their close flyby 7-11 Gyr ago. Here we explore the Local Group history in MOND in more detail using a simplified first-ever self-consistent simulation. We note the features induced by the encounter in the simulation and identify their possible real counterparts. The simulation was set to approximately reproduce the observed MW and M31 masses, effective radii, separation, relative velocity and disk inclinations. We used the publicly available adaptive-mesh-refinement code Phantom of RAMSES. In the simulation, matter was transferred from the MW to M31 along a tidal tail. The encounter induced formation of several structures resembling the peculiarities of the Local Group. Most notably: 1) A rotating planar structure formed around M31 from the transferred material. It had a size similar to the observed satellite plane and was oriented edge-on to the simulated MW, just as the real one. 2) The same structure also resembled the tidal features observed around M31 by its size and morphology. 3) A warp in the MW developed with an amplitude and orientation similar to that observed. 4) A cloud of particles formed around the simulated MW, with the extent of the actual MW satellite system. The encounter did not end by merging in a Hubble time. The simulation thus demonstrated that MOND can possibly explain many peculiarities of the Local Group and, moreover, that non-merging galaxy encounters in MOND can produce tidal features in galaxies.

Recently, a new form of dark matter has been suggested to naturally reproduce the empirically successful aspects of Milgrom's law in galaxies. The dark matter particle candidates are axion-like, with masses of order eV and strong self-interactions. They Bose-Einstein condense into a superfluid phase in the central regions of galaxy halos. The superfluid phonon excitations in turn couple to baryons and mediate an additional long-range force. For a suitable choice of the superfluid equation of state, this force can mimic Milgrom's law. In this paper we develop in detail some of the main phenomenological consequences of such a formalism, by revisiting the expected dark matter halo profile in the presence of an extended baryon distribution. In particular, we show how rotation curves of both high and low surface brightness galaxies can be reproduced, with a slightly rising rotation curve at large radii in massive high surface brightness galaxies, thus subtly different from Milgrom's law. We finally point out other expected differences with Milgrom's law, in particular in dwarf spheroidal satellite galaxies, tidal dwarf galaxies, and globular clusters, whose Milgromian or Newtonian behavior depends on the position with respect to the superfluid core of the host galaxy. We also expect ultra-diffuse galaxies within galaxy clusters to have velocities slightly above the baryonic Tully-Fisher relation. Finally, we note that, in this framework, photons and gravitons follow the same geodesics, and that galaxy-galaxy lensing, probing larger distances within galaxy halos than rotation curves, should follow predictions closer to the standard cosmological model than those of Milgrom's law.

The orbital motion of non-contact double-lined spectroscopic binaries (SB2), with periods of a few tens of days to several years, holds unique accurate informations on individual stellar masses, that only long-term monitoring can unlock. The combination of radial velocity measurements from high-resolution spectrographs and astrometric measurements from high-precision interferometers allows the derivation of SB2 components masses down to the percent precision. Since 2010, we observed a large sample of SB2 with the SOPHIE spectrograph at the Observatoire de Haute-Provence, aiming at the derivation of orbital elements with sufficient accuracy to obtain masses of components with relative errors as low as 1$\%$ when the astrometric measurements of the Gaia satellite will be taken into account. In this paper we present the results from six years of observations of 14 SB2 systems with periods ranging from 33 to 4185 days. Using the TODMOR algorithm we computed radial velocities from the spectra, and then derived the orbital elements of these binary systems. The minimum masses of the 28 stellar components are then obtained with a sample average accuracy of 1.0$\pm$0.2$\,\%$. Combining the radial velocities with existing interferometric measurements, we derived the masses of the primary and secondary components of HIP 61100, HIP 95995 and HIP 101382 with relative errors for components (A,B) of respectively (2.0, 1.7)$\,\%$, (3.7, 3.7)$\,\%$, and (0.2, 0.1)$\,\%$. Using the Cesam2k stellar evolution code, we could constrain the initial He-abundance, age and metallicity for HIP 61100 and HIP 95995.

The radial velocity (RV) of a single star is easily obtained from cross-correlation of the spectrum with a template, but the treatment of double-lined spectroscopic binaries (SB2s) is more difficult. Two different approaches were applied to a set of SB2s: the fit of the cross-correlation function with two normal distributions, and the cross-correlation with two templates, derived with the TODCOR code. It appears that the minimum masses obtained through the two methods are sometimes rather different, although their estimated uncertainties are roughly equal. Moreover, both methods induce a shift in the zero point of the secondary RVs, but it is less pronounced for TODCOR. All-in-all the comparison between the two methods is in favour of TODCOR.

Accurate stellar masses are requested in order to improve our understanding of stellar interiors, but they are still rather rare. Fortunately, the forthcoming Gaia Mission will provide astrometric measurements permitting the derivation of the orbital inclinations of nearby binaries which are also observed as double-lined spectroscopic binaries (SB2s) with ground-based telescopes. A programme of radial velocity (RV) measurements was initiated in 2010 with the Sophie spectrograph of the Haute-Provence observatory in order to derive accurate SB2 orbits for a large set of stars. Therefore, combined SB2+astrometric orbits will be derived thanks to Gaia, and masses with errors around 1~\% are expected for both components. The programme includes 70 SB2s, and the accurate SB2 orbits of 24 of them were already derived. In addition, two complementary programmes devoted to southern stars or to late-type dwarf stars were also initiated with the HERMES and the CARMENES spectrographs, respectively. Interferometric measurements were obained with the VLTI/PIONIER for 7 SB2s, and were taken from other sources for 4 others. Currently, combined "visual binary" (VB) +SB2 solutions were derived for 7 binaries, leading to the masses of the components and to the parallaxes. The parallaxes from the Hipparcos 2 catalogue were corrected for orbital motion and compared to our solution, confirming the high quality of Hipparcos 2.

It has recently been proposed, by assuming that dark matter is a superfluid, that MOND-like effects can be achieved on small scales whilst preserving the success of $\Lambda$CDM on large scales. Here we aim to provide the first set of spherical models of galaxy clusters in the context of superfluid dark matter. We first outline the theoretical structure of the superfluid core and the surrounding "normal phase" dark halo of quasi-particles in thermal equlibrium. The latter should encompass the largest part of galaxy clusters. Here, we set the SfDM transition at the radius where the density and pressure of the superfluid and normal phase coincides, neglecting the effect of phonons in the suprefluid core. We then apply the theory to a sample of galaxy clusters, and directly compare the SfDM predicted mass profiles to data. We find that the superfluid formulation can reproduce the X-ray dynamical mass profile of clusters. The SfDM fits however display slight under-predictions of the gravity in the central regions which might be partly related to our neglecting of the effect of phonons in these regions. We conclude that this superfluid formulation is successful in describing galaxy clusters, but further work will be needed to determine whether the parameter choice is consistent with galaxies. Our model could be made more realistic by exploring non-sphericity and the SfDM transition condition we impose.

Kinematically cold tidal streams of globular clusters (GC) are excellent tracers of the Galactic gravitational potential at moderate Galactocentric distances, and can also be used as probes of the law of gravity on Galactic scales. Here, we compare for the first time the generation of such streams in Newtonian and Milgromian gravity (MOND). We first compute analytical results to investigate the expected shape of the GC gravitational potential in both frameworks, and we then run N-body simulations with the Phantom of Ramses code. We find that the GCs tend to become lopsided in MOND. This is a consequence of the external field effect which breaks the strong equivalence principle. When the GC is filling its tidal radius the lopsidedness generates a strongly asymmetric tidal stream. In Newtonian dynamics, such markedly asymmetric streams can in general only be the consequence of interactions with dark matter subhalos, giant molecular clouds, or interaction with the Galactic bar. In these Newtonian cases, the asymmetry is the consequence of a very large gap in the stream, whilst in MOND it is a true asymmetry. This should thus allow us in the future to distinguish these different scenarios by making deep observations of the environment of the asymmetric stellar stream of Palomar 5. Moreover, our simulations indicate that the high internal velocity dispersion of Palomar 5 for its small stellar mass would be natural in MOND.

The existence of a spatially thin, kinematically coherent Disk of Satellites (DoS) around the Milky Way (MW), is a problem that often garners vivacious debate in the literature or at scientific meetings. One of the most recent incarnations of this wrangle occurred with two papers by Maji et al, who argued that these structures "maybe a misinterpretation of the data". These claims are in stark contrast to previous works. Motivated by this and other recent publications on this problem, we discuss necessary considerations to make, observational effects to consider, and pitfalls to avoid when investigating satellite galaxy planes such as the MW's DoS. In particular, we emphasize that conclusions need to have a statistical basis including a determination of the significance of satellite alignments, observational biases must not be ignored, and measurement errors (e.g. for proper motions) need to be considered. We discuss general problems faced by attempts to determine the dynamical stability of the DoS via orbit integrations of MW satellite galaxies, and demonstrate that to interpret simulations, it is helpful to compare them with a null case of isotropically distributed satellite positions and velocities. Based on these criteria, we find that the conclusions of Maji et al. do not hold up to scrutiny, and that their hydrodynamic cosmological simulation of a single host shows no evidence for a significant kinematic coherence among the simulated satellite galaxies, in contrast to the observed MW system.

We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.

The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.

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.

Tidal streams of disrupting dwarf galaxies orbiting around their host galaxy offer a unique way to constrain the shape of galactic gravitational potentials. Such streams can be used as leaning tower gravitational experiments on galactic scales. The most well motivated modification of gravity proposed as an alternative to dark matter on galactic scales is Milgromian dynamics (MOND), and we present here the first ever N-body simulations of the dynamical evolution of the disrupting Sagittarius dwarf galaxy in this framework. Using a realistic baryonic mass model for the Milky Way, we attempt to reproduce the present-day spatial and kinematic structure of the Sagittarius dwarf and its immense tidal stream that wraps around the Milky Way. With very little freedom on the original structure of the progenitor, constrained by the total luminosity of the Sagittarius structure and by the observed stellar mass-size relation for isolated dwarf galaxies, we find reasonable agreement between our simulations and observations of this system. The observed stellar velocities in the leading arm can be reproduced if we include a massive hot gas corona around the Milky Way that is flattened in the direction of the principal plane of its satellites. This is the first time that tidal dissolution in MOND has been tested rigorously at these mass and acceleration scales.

We explore the correlations between velocity and metallicity and the possible distinct chemical signatures of the velocity over-densities of the local Galactic neighbourhood. We use the large spectroscopic survey RAVE and the Geneva Copenhagen Survey. We compare the metallicity distribution of regions in the velocity plane ($v_R,v_\phi$) with that of their symmetric counterparts ($-v_R,v_\phi$). We expect similar metallicity distributions if there are no tracers of a sub-population (e.g., a dispersed cluster, accreted stars), if the disk of the Galaxy is axisymmetric, and if the orbital effects of the spiral arms and the bar are weak. We find that the metallicity-velocity space of the solar neighbourhood is highly patterned. A large fraction of the velocity plane shows differences in the metallicity distribution when comparing symmetric $v_R$ regions. The typical differences in the median metallicity are of $0.05$ dex with a statistical significance of at least $95\%$, and with values up to $0.6$ dex. For low azimuthal velocity $v_\phi$, stars moving outwards in the Galaxy have on average higher metallicity than those moving inwards. These include stars in the Hercules and Hyades moving groups and other velocity branch-like structures. For higher $v_\phi$, the stars moving inwards have higher metallicity than those moving outwards. The most likely interpretation of the metallicity asymmetry is that it is due to the orbital effects of the bar and the radial metallicity gradient of the disk. We present a simulation that supports this idea. We have also discovered a positive gradient in $v_\phi$ with respect to metallicity at high metallicities, apart from the two known positive and negative gradients for the thick and thin disks, respectively.

It was recently proposed that the effects usually attributed to particle dark matter on galaxy scales are due to the displacement of dark energy by baryonic matter, a paradigm known as emergent gravity. This formalism leads to predictions similar to Modified Newtonian Dynamics (MOND) in spherical symmetry, but not quite identical. In particular, it leads to a well defined transition between the Newtonian and the modified gravitational regimes, a transition depending on both the Newtonian acceleration and its first derivative with respect to radius. Under the hypothesis of the applicability of this transition to aspherical systems, we investigate whether it can reproduce observed galaxy rotation curves. We conclude that the formula leads to marginally acceptable fits with strikingly low best-fit distances, low stellar mass-to-light ratios, and a low Hubble constant. In particular, some unobserved wiggles are produced in rotation curves because of the dependence of the transition on the derivative of the Newtonian acceleration, leading, even in the most favorable case, to systematically less good fits than MOND. Then, applying the predicted transition from emergent gravity in a regime where it should a priori be applicable, i.e. in spherical symmetry and outside of the bulk of matter, we show that the predictions for the secular advances of Solar System planets' perihelia are discrepant with the data by seven orders of magnitude, ruling out the present emergent gravity weak-field formula with high confidence.

We report on the detection of a small overdensity of stars in velocity space with systematically higher Galactocentric rotation velocity than the Sun by about 20 km s$^{-1}$ in the $Gaia$ Data Release 1 Tycho-Gaia astrometric solution (TGAS) data. We find these fast Galactic rotators more clearly outside of the Solar radius, compared to inside of the Solar radius. In addition, the velocity of the fast Galactic rotators is independent of the Galactocentric distance up to $R-R_{\odot}\sim0.6$ kpc. Comparing with numerical models, we qualitatively discuss that a possible cause of this feature is the co-rotation resonance of the Perseus spiral arm, where the stars in peri-centre phase in the trailing side of the Perseus spiral arm experience an extended period of acceleration owing to the torque from the Perseus arm.

The length and pattern speed of the Milky Way bar are still controversial. Photometric and spectroscopic surveys of the inner Galaxy, as well as gas kinematics, favour a long and slowly rotating bar, with corotation around a Galactocentric radius of 6 kpc. On the other hand, the existence of the Hercules stream in local velocity space favours a short and fast bar with corotation around 4 kpc. This follows from the fact that the Hercules stream looks like a typical signature of the outer Lindblad resonance of the bar. As we showed recently, reconciling this local stream with a slow bar would need to find a yet unknown alternative explanation, based for instance on the effect of spiral arms. Here, by combining the TGAS catalogue of the Gaia DR1 with LAMOST radial velocities, we show that the position of Hercules in velocity space as a function of radius in the outer Galaxy indeed varies exactly as predicted by fast bar models with a pattern speed no less than 1.8 times the circular frequency at the Sun's position.

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.

Both the three-dimensional density of red clump giants and the gas kinematics in the inner Galaxy indicate that the pattern speed of the Galactic bar could be much lower than previously estimated. Here, we show that such slow bar models are unable to reproduce the bimodality observed in local stellar velocity space. We do so by computing the response of stars in the Solar neighbourhood to the gravitational potential of slow and fast bars, in terms of their perturbed distribution function in action-angle space up to second order, as well as by identifying resonantly trapped orbits. We also check that the bimodality is unlikely to be produced through perturbations from spiral arms, and conclude that, contrary to gas kinematics, local stellar kinematics still favour a fast bar in the Milky Way, with a pattern speed of the order of almost twice (and no less than 1.8 times) the circular frequency at the Sun's position. This leaves open the question of the nature of the long flat extension of the bar in the Milky Way.

Together with interstellar turbulence, gravitation is one key player in star formation. It acts both at galactic scales in the assembly of gas into dense clouds, and inside those structures for their collapse and the formation of pre-stellar cores. To understand to what extent the large scale dynamics govern the star formation activity of galaxies, we present hydrodynamical simulations in which we generalise the behaviour of gravity to make it differ from Newtonian dynamics in the low acceleration regime. We focus on the extreme cases of interacting galaxies, and compare the evolution of galaxy pairs in the dark matter paradigm to that in the Milgromian Dynamics (MOND) framework. Following up on the seminal work by Tiret & Combes, this paper documents the first simulations of galaxy encounters in MOND with a detailed Eulerian hydrodynamical treatment of baryonic physics, including star formation and stellar feedback. We show that similar morphologies of the interacting systems can be produced by both the dark matter and MOND formalisms, but require a much slower orbital velocity in the MOND case. Furthermore, we find that the star formation activity and history are significantly more extended in space and time in MOND interactions, in particular in the tidal debris. Such differences could be used as observational diagnostics and make interacting galaxies prime objects in the study of the nature of gravitation at galactic scales.

We investigate models of the Milky Way disc taking into account simultaneously the bar and a two-armed quasi-static spiral pattern. Away from major resonance overlaps, the mean stellar radial motions in the plane are essentially a linear superposition of the isolated effects of the bar and spirals. Thus, provided the bar is strong enough, even in the presence of spiral arms, these mean radial motions are predominantly affected by the Galactic bar for large scale velocity fluctuations. This is evident when comparing the peculiar line-of-sight velocity power spectrum of our coupled models with bar-only models. However, we show how forthcoming spectroscopic surveys could disentangle bar-only non-axisymmetric models of the Galaxy from models in which spiral arms have a significant amplitude. We also point out that overlaps of low-order resonances are sufficient to enhance stellar churning within the disc, even when the spirals amplitude is kept constant. Nevertheless, for churning to be truly non-local, stronger or (more likely) transient amplitudes would be needed: otherwise the disc is actually mostly unaffected by churning in the present models. Finally, regarding vertical breathing modes, the combined effect of the bar and spirals on vertical motions is a clear non-linear superposition of the isolated effects of both components, significantly superseding the linear superposition of modes produced by each perturber separately, thereby providing an additional effect to consider when analysing the observed breathing mode of the Galactic disc in the extended Solar neighbourhood.

Assuming virial equilibrium and Newtonian dynamics, low-mass early-type galaxies have larger velocity dispersions than expected from the amount of baryons they contain. The conventional interpretation of this finding is that their dynamics is dominated by non-baryonic matter. However, there is also strong evidence that many low-mass early-type galaxies formed as tidal dwarf galaxies, which would contain almost no dark matter. Using an extensive catalogue of early-type galaxies, we therefore discuss how the internal dynamics of early-type galaxies in general can be understood by replacing the assumption of non-baryonic dark matter with two alternative assumptions. The first assumption is that Milgromian dynamics (i.e., MOND) is valid, which changes the effective gravitational force in the weak-field limit. The second assumption is that binary stars affect the observed line-of-sight velocity dispersions. Some moderate discrepancies between observed and predicted velocity dispersions remain also when these effects are implemented. Nevertheless, the observed velocity dispersions in early-type galaxies can then easily be explained without invoking the presence of non-baryonic dark matter in them, but with already documented variations of the galaxy-wide stellar initial mass function and non-equilibrium dynamics in some of the low-mass early-type galaxies.

Starting from an axisymmetric equilibrium distribution function (DF) in action space, representing a Milky Way thin disc stellar population, we use the linearized Boltzmann equation to explicitly compute the response to a three-dimensional spiral potential in terms of the perturbed DF. This DF, valid away from the main resonances, allows us to investigate a snapshot of the velocity distribution at any given point in three-dimensional configuration space. Moreover, the first order moments of the DF give rise to non-zero radial and vertical bulk flows -- namely breathing modes -- qualitatively similar to those recently observed in the extended Solar neighbourhood. We show that these analytically predicted mean stellar motions are in agreement with the outcome of test-particle simulations. Moreover, we estimate for the first time the reduction factor for the vertical bulk motions of a stellar population compared to the case of a cold fluid. Such an explicit expression for the full perturbed DF of a thin disc stellar population in the presence of spiral arms will be helpful in order to dynamically interpret the detailed information on the Milky Way disc stellar kinematics that will be provided by upcoming large astrometric and spectroscopic surveys of the Galaxy.

Since its publication 1983, Milgromian dynamics (aka MOND) has been very successful in modeling the gravitational potential of galaxies from baryonic matter alone. However, the dynamical modeling has long been an unsolved issue. In particular, the setup of a stable galaxy for Milgromian N-body calculations has been a major challenge. Here, we show a way to set up disc galaxies in MOND for calculations in the PHANTOM OF RAMSES (PoR) code by L\"ughausen (2015) and Teyssier (2002). The method is done by solving the QUMOND Poisson equations based on a baryonic and a phantom dark matter component. The resulting galaxy models are stable after a brief settling period for a large mass and size range. Simulations of single galaxies as well as colliding galaxies are shown.

We present an analysis of the presence of substructures in the stellar stream of the Palomar 5 globular cluster, as derived from Sloan Digital Sky Survey data. Using a matched filter technique, we recover the positions and sizes of overdensities reported in previous studies. To explore the reality of these structures, we also create an artificial model of the stream, in which we construct a realistic background on top of which we add a perfectly smooth stream structure, taking into account the effects of photometric completeness and interstellar extinction. We find that the smooth artificial stream then shows similarly-pronounced substructures as the real structure. Interestingly, our best-fit N-body simulation does display real projected density variations linked to stellar epicyclic motions, but these become less significant when taking into account the SDSS star-count constraints. The substructures found when applying our matched filter technique to the N-body particles converted into observable stars are thus mostly unrelated to these epicyclic motions. This analysis suggests that the majority of the previously-detected substructures along the tidal tail of Palomar 5 are artefacts of observational inhomogeneities.

In anticipation of the Gaia astrometric mission, a large sample of spectroscopic binaries has been observed since 2010 with the SOPHIE spectrograph at the Haute--Provence Observatory. Our aim is to derive the orbital elements of double-lined spectroscopic binaries (SB2s) with an accuracy sufficient to finally obtain the masses of the components with relative errors as small as 1 % when the astrometric measurements of Gaia are taken into account. In this paper we present the results from five years of observations of 10 SB2 systems with periods ranging from 37 to 881 days. Using the TODMOR algorithm we computed radial velocities from the spectra, and then derived the orbital elements of these binary systems. The minimum masses of the components are then obtained with an accuracy better than 1.2 % for the ten binaries. Combining the radial velocities with existing interferometric measurements, we derived the masses of the primary and secondary components of HIP 87895 with an accuracy of 0.98 % and 1.2 % respectively.

This article is the detailed version of a paper on dark matter, dark energy, and modified gravity, published in the December 2015-January 2016 special issue of "La Recherche" (in French)

We present a dynamical measurement of the tangential motion of the Andromeda system, the ensemble consisting of the Andromeda Galaxy (M31) and its satellites. The system is modelled as a structure with cosmologically-motivated velocity dispersion and density profiles, and we show that our method works well when tested using the most massive substructures in high-resolution Lambda Cold Dark Matter simulations. Applied to the sample of 40 currently-known galaxies of this system, we find a value for the transverse velocity of 164.4 +/- 61.8 km/s (v_East = -111.5 +/- 70.2 km/s and v_North = 99.4 +/- 60.0 km/s), significantly higher than previous estimates of the proper motion of M31 itself. This result has significant implications on estimates of the mass of the Local Group, as well as on its past and future history.

We investigate the degree to which the inclusion of baryonic physics can overcome two long-standing problems of the standard cosmological model on galaxy scales: (i) the problem of satellite planes around Local Group galaxies, and (ii) the "too big to fail" problem. By comparing dissipational and dissipationless simulations, we find no indication that the addition of baryonic physics results in more flattened satellite distributions around Milky-Way-like systems. Recent claims to the contrary are shown to derive in part from a non-standard metric for the degree of flattening, which ignores the satellites' radial positions. If the full 3D positions of the satellite galaxies are considered, none of the simulations we analyse reproduce the observed flattening nor the observed degree of kinematic coherence of the Milky Way satellite system. Our results are consistent with the expectation that baryonic physics should have little or no influence on the structure of satellite systems on scales of hundreds of kiloparsecs. Claims that the "too big to fail" problem can be resolved by the addition of baryonic physics are also shown to be problematic.

In anticipation of the Gaia astrometric mission, a sample of spectroscopic binaries is being observed since 2010 with the Sophie spectrograph at the Haute--Provence Observatory. Our aim is to derive the orbital elements of double-lined spectroscopic binaries (SB2s) with an accuracy sufficient to finally obtain the masses of the components with relative errors as small as 1 % when combined with Gaia astrometric measurements. In order to validate the masses derived from Gaia, interferometric observations are obtained for three SB2s in our sample with F-K components: HIP 14157, HIP 20601 and HIP 117186. The masses of the six stellar components are derived. Due to its edge-on orientation, HIP 14157 is probably an eclipsing binary. We note that almost all the derived masses are a few percent larger than the expectations from the standard spectral-type-mass calibration and mass-luminosity relation. Our calculation also leads to accurate parallaxes for the three binaries, and the Hipparcos parallaxes are confirmed.

The DiskMass survey recently provided measurements of the vertical velocity dispersions of disk stars in a sample of nearly face-on galaxies. By setting the disk scale-heights to be equal to those of edge-on galaxies with similar scale-lengths, it was found that these disks must be sub-maximal, with surprisingly low K-band mass-to-light ratios of the order of $M_\star/L_K \simeq 0.3 M_\odot/L_\odot$. This study made use of a simple relation between the disk surface density and the measured velocity dispersion and scale height of the disk, neglecting the shape of the rotation curve and the dark matter contribution to the vertical force, which can be especially important in the case of sub-maximal disks. Here, we point out that these simplifying assumptions led to an overestimation of the stellar mass-to-light ratios. Relaxing these assumptions, we compute even lower values than previously reported for the mass-to-light ratios, with a median $M_\star/L_K \simeq 0.18 M_\odot/L_\odot$, where 14 galaxies have $M_\star/L_K < 0.11$. Invoking prolate dark matter halos made only a small difference to the derived $M_\star/L_K$, although extreme prolate halos ($q>1.5$ for the axis ratios of the potential) might help. The cross-terms in the Jeans equation are also generally negligible. These deduced K-band stellar mass-to-light ratios are even more difficult to reconcile with stellar population synthesis models than the previously reported ones.

The Modified Newtonian Dynamics (MOND) paradigm generically predicts that the external gravitational field in which a system is embedded can produce effects on its internal dynamics. In this communication, we first show that this External Field Effect can significantly improve some galactic rotation curves fits by decreasing the predicted velocities of the external part of the rotation curves. In modified gravity versions of MOND, this External Field Effect also appears in the Solar System and leads to a very good way to constrain the transition function of the theory. A combined analysis of the galactic rotation curves and Solar System constraints (provided by the Cassini spacecraft) rules out several classes of popular MOND transition functions, but leaves others viable. Moreover, we show that LISA Pathfinder will not be able to improve the current constraints on these still viable transition functions.

We present an estimate of the third integral of motion for axisymmetric three-dimensional potentials. This estimate is based on a Staeckel approximation and is explicitly written as a function of the potential. We tested this scheme for the Besancon Galactic model and two other disc-halo models and find that orbits of disc stars have an accurately conserved third quasi integral. The accuracy ranges from of 0.1% to 1% for heights varying from z = 0~kpc to z= 6 kpc and Galactocentric radii R from 5 to 15kpc. We also tested the usefulness of this quasi integral in analytic distribution functions of disc stellar populations: we show that the distribution function remains approximately stationary and that it allows to recover the potential and forces by applying Jeans equations to its moments.

We investigate the kinematic signatures induced by spiral and bar structure in a set of simulations of Milky Way-sized spiral disc galaxies. The set includes test particle simulations that follow a quasi-stationary density wave-like scenario with rigidly rotating spiral arms, and $N$-body simulations that host a bar and transient, co-rotating spiral arms. From a location similar to that of the Sun, we calculate the radial, tangential and line-of-sight peculiar velocity fields of a patch of the disc and quantify the fluctuations by computing the power spectrum from a two-dimensional Fourier transform. We find that the peculiar velocity power spectrum of the simulation with a bar and transient, co-rotating spiral arms fits very well to that of APOGEE red clump star data, while the quasi-stationary density wave spiral model without a bar does not. We determine that the power spectrum is sensitive to the number of spiral arms, spiral arm pitch angle and position with respect to the spiral arm. However, it is necessary to go beyond the line of sight velocity field in order to distinguish fully between the various spiral models with this method. We compute the power spectrum for different regions of the spiral discs, and discuss the application of this analysis technique to external galaxies.

We investigate, using the Gaia-ESO Survey internal Data-Release 2, the properties of the double sequence of the Milky Way discs (defined chemically as the high-alpha and low-alpha populations), and discuss their compatibility with discs defined by other means such as metallicity, kinematics or positions. This investigation uses two different approaches: in velocity space for stars located in the extended Solar neighbourhood, and in chemical space for stars at different ranges of Galactocentric radii and heights from the plane. The separation we find in velocity space allows us to investigate, in a novel manner, the extent in metallicity of each of the two sequences, identifying them with the two discs, without making any assumption about the shape of their metallicity distribution functions. Then, using the separation in chemical space, we characterise the spatial variation of the slopes of the [alpha/Fe] - [Fe/H] sequences for the thick and thin discs and the way in which the relative proportions of the two discs change across the Galaxy. We find that the thick disc (high-alpha sequence), extends up to [Fe/H]~ +0.2 and the thin disc (low-alpha sequence), at least down to [Fe/H]~ -0.8. Radial and vertical gradients in alpha-abundances are found for the thin disc, with mild spatial variations in its [alpha/Fe] - [Fe/H] paths, whereas for the thick disc we do not detect any such spatial variations. The small variations in the spatial [alpha/Fe] - [Fe/H] paths of the thin disc do not allow us to distinguish between formation models of this structure. On the other hand, the lack of radial gradients and [alpha/Fe] - [Fe/H] variations for the thick disc indicate that the mechanism responsible for the mixing of the metals in the young Galaxy (e.g. radial stellar migration or turbulent gaseous disc) was more efficient before the (present) thin disc started forming.

We have recently shown (Ibata et al. 2014) that pairs of satellite galaxies located diametrically opposite each other around their host possess predominantly anti-correlated velocities. This is consistent with a scenario in which $\sim 50$% of satellite galaxies belong to kinematically-coherent rotating planar structures, similar to those detected around the giant galaxies of the Local Group. Here we extend this analysis, examining the incidence of satellites of giant galaxies drawn from an SDSS photometric redshift catalog. We find that there is a $\sim 17$% overabundance ($> 3 \sigma$ significance) of candidate satellites at positions diametrically opposite a spectroscopically confirmed satellite. We show that cosmological simulations do not possess this property when the contamination is included, and that there are in fact, after subtracting contamination, 2 to 3 times more satellites diametrically opposed to a spectroscopically confirmed satellite than at $90\deg$ from it. We also examine the correlation between the satellite pair positions and the orientation of the host galaxy major axis. We find that those satellite pairs with anti-correlated velocities have a strong preference ($\sim 3:1$) to align with the major axis of the host whereas those with correlated velocities display the opposite behavior. This correlation of the satellite alignments appears to be stronger than the well-documented preference of satellites to be located close to the major axis of their host. We finally show that repeating a similar analysis to Ibata et al. (2014) with same-side satellites is generally hard to interpret, but is not inconsistent with our previous results when strong quality-cuts are applied on the sample. All these unexpected correlations strongly suggest that a substantial fraction of satellite galaxies are causally-linked in their formation and evolution.

Evidence for non-zero mean stellar velocities in the direction perpendicular to the Galactic plane has been accumulating from various recent large spectroscopic surveys. Previous analytical and numerical work has shown that a "breathing mode" of the Galactic disc, similar to what is observed in the Solar vicinity, can be the natural consequence of a non-axisymmetric internal perturbation of the disc. Here we provide a general analytical framework, in the context of perturbation theory, allowing us to compute the vertical bulk motions generated by a single internal perturber (bar or spiral pattern). In the case of the Galactic bar, we show that these analytically predicted bulk motions are well in line with the outcome of a numerical simulation. The mean vertical motions induced by the Milky Way bar are small (mean velocity of less than 1 km/sec) and cannot be responsible alone for the observed breathing mode, but they are existing. Our analytical treatment is valid close to the plane for all the non-axisymmetric perturbations of the disc that can be described by small-amplitude Fourier modes. Further work should study how the coupling of multiple internal perturbers and external perturbers is affecting the present analytical results.

This article explores the agreement between the predictions of Modified Newtonian Dynamics (MOND) and the rotation curves and stellar velocity dispersion profiles measured by the DiskMass Survey. A bulge-disk decomposition was made for each of the thirty published galaxies, and a MOND Poisson solver was used to simultaneously compute, from the baryonic mass distributions, model rotation curves and vertical velocity dispersion profiles, which were compared to the measured values. The two main free parameters, the stellar disk's mass-to-light ratio ($M/L$) and its exponential scale-height ($h_z$), were estimated by Markov Chain Monte Carlo modelling. The average best-fit K-band stellar mass-to-light ratio was $M/L \simeq 0.55 \pm 0.15$. However, to match the DiskMass Survey data, the vertical scale-heights would have to be in the range $h_z=200$ to $400$ pc which is a factor of two lower than those derived from observations of edge-on galaxies with a similar scale-length. The reason is that modified gravity versions of MOND characteristically require a larger $M/L$ to fit the rotation curve in the absence of dark matter and therefore predict a stronger vertical gravitational field than Newtonian models. It was found that changing the MOND acceleration parameter, the shape of the velocity dispersion ellipsoid, the adopted vertical distribution of stars, as well as the galaxy inclination, within any realistic range, all had little impact on these results.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

We present a study of the intrinsic deprojected ellipticity distribution of the satellite dwarf galaxies of the Andromeda galaxy, assuming that their visible components have a prolate shape, which is a natural outcome of simulations. Different possibilities for the orientation of the major axis of the prolate dwarf galaxies are tested, pointing either as close as possible to the radial direction towards the centre of Andromeda, or tangential to the radial direction, or with a random angle in the plane that contains the major axis and the observer. We find that the mean intrinsic axis ratio is ~ 1/2, with small differences depending on the assumed orientation of the population. Our deprojections also suggest that a significant fraction of the satellites, ~ 10%, are tidally disrupted remnants. We find that there is no evidence of any obvious difference in the morphology and major axis orientation between satellites that belong to the vast thin plane of co-rotating galaxies around Andromeda and those that do not belong to this structure.

We study the kinematics of a local sample of stars, located within a cylinder of 500 pc radius centered on the Sun, in the RAVE dataset. We find clear asymmetries in the $v_R$-$v_\phi$ velocity distributions of thin and thick disk stars: here are more stars moving radially outwards for low azimuthal velocities and more radially inwards for high azimuthal velocities. Such asymmetries have been previously reported for the thin disk as being due to the Galactic bar, but this is the first time that the same type of structures are seen in the thick disk. Our findings imply that the velocities of thick disk stars should no longer be described by Schwarzschild's, multivariate Gaussian or purely axisymmetric distributions. Furthermore, the nature of previously reported substructures in the thick disk needs to be revisited as these could be associated with dynamical resonances rather than to accretion events. It is clear that dynamical models of the Galaxy must fit the 3D velocity distributions of the disks, rather than the projected 1D, if we are to understand the Galaxy fully.

We review most dynamical constraints on the gravitational field of spiral galaxies in general, and of the Milky Way in particular. Such constraints are of prime importance for determining the characteristics of the putative dark matter haloes of galaxies. For the Milky Way, we review observational constraints in the inner parts (cored or cusped dark matter distribution, maximum disk or not), in the solar neighbourhood (local dark matter density) and in the outer parts (virial mass and triaxial shape of the dark matter halo). We also point out various caveats, systematic effects, and large current uncertainties. Many fundamental parameters such as the local circular velocity are poorly known, evidence for triaxiality of the dark halo is shaky, and different estimates of the virial mass as well as of the local dark matter density vary by at least a factor of two. We however argue that the current best-fit value for the local dark matter density, which should be used as a benchmark for direct dark matter detection searches, is of the order of 0.5 GeV/cm3. We also explain why alternatives to particle dark matter on galactic scales should still be very seriously considered.

We determine the Galactic potential in the solar neigbourhood from RAVE observations. We select red clump stars for which accurate distances, radial velocities, and metallicities have been measured. Combined with data from the 2MASS and UCAC catalogues, we build a sample of 4600 red clump stars within a cylinder of 500 pc radius oriented in the direction of the South Galactic Pole, in the range of 200 pc to 2000 pc distances. We deduce the vertical force and the total mass density distribution up to 2 kpc away from the Galactic plane by fitting a distribution function depending explicitly on three isolating integrals of the motion in a separable potential locally representing the Galactic one with four free parameters. Because of the deep extension of our sample, we can determine nearly independently the dark matter mass density and the baryonic disc surface mass density. We find (i) at 1kpc Kz/(2piG) = 68.5 pm 1.0 Msun/pc2, and (ii) at 2 kpc Kz/(2piG) = 96.9 pm 2.2 Msun/pc2. Assuming the solar Galactic radius at R0 = 8.5 kpc, we deduce the local dark matter density rhoDM (z=0) = 0.0143 pm 0.0011Msun pc3 = 0.542 pm 0.042 Gev/cm3 and the baryonic surface mass density Sigma = 44.4 pm 4.1 Msun/pc2 . Our results are in agreement with previously published Kz determinations up to 1 kpc, while the extension to 2 kpc shows some evidence for an unexpectedly large amount of dark matter. A flattening of the dark halo of order 0.8 can produce such a high local density in combination with a circular velocity of 240 km/s . Another explanation, allowing for a lower circular velocity, could be the presence of a secondary dark component, a very thick disc resulting either from the deposit of dark matter from the accretion of multiple small dwarf galaxies, or from the presence of an effective phantom thick disc in the context of effective galactic-scale modifications of gravity.

In anticipation of the Gaia astrometric mission, a large sample of spectroscopic binaries is being observed since 2010 with the Sophie spectrograph at the Haute--Provence Observatory. Our aim is to derive the orbital elements of double-lined spectroscopic binaries (SB2s) with an accuracy sufficient to finally obtain the masses of the components with relative errors as small as 1% when the astrometric measurements of Gaia are taken into account. Simultaneously, the luminosities of the components in the Gaia photometric band G will also be obtained. Our observation program started with 200 SBs, including 152 systems that were only known as single-lined. Thanks to the high efficiency of the Sophie spectrograph, an additional component was found for 25 SBs. After rejection of 5 multiple systems, 20 new SB2s were retained, including 8 binaries with evolved primary, and their mass ratios were derived. Our final sample contains 68 SB2s, including 2 late-type giants and 10 other evolved stars.

Recent work has shown that both the Milky Way and the Andromeda galaxies possess the unexpected property that their dwarf satellite galaxies are aligned in thin and kinematically coherent planar structures. It is now important to evaluate the incidence of such planar structures in the larger galactic population, since the Local Group may not be a sufficiently representative environment. Here we report that the measurement of the velocity of pairs of diametrically opposed galaxy satellites provides a means to determine statistically the prevalence of kinematically coherent planar alignments. In the local universe (redshift $z<0.05$), we find that such satellite pairs out to a galactocentric distance of $150$ kpc are preferentially anti-correlated in their velocities (99.994% confidence level), and that the distribution of galaxies in the larger scale environment (beyond 150 kpc and up to $\approx 2$ Mpc) is strongly elongated along the axis joining the inner satellite pair ($>7\sigma$ confidence). Our finding may indicate that co-rotating planes of satellites, similar to that seen around the Andromeda galaxy, are ubiquitous in nature, while their coherent motion also suggests that they are a significant repository of angular momentum on $\sim 100$ kpc scales.

In this article we continue to test cosmological models centred on Modified Newtonian Dynamics (MOND) with light sterile neutrinos, which could in principle be a way to solve the fine-tuning problems of the standard model on galaxy scales while preserving successful predictions on larger scales. Due to previous failures of the simple MOND cosmological model, here we test a speculative model where the modified gravitational field is produced only by the baryons and the sterile neutrinos produce a purely Newtonian field (hence Modified Baryonic Dynamics). We use two component cosmological simulations to separate the baryonic N-body particles from the sterile neutrino ones. The premise is to attenuate the over-production of massive galaxy cluster halos which were prevalent in the original MOND plus light sterile neutrinos scenario. Theoretical issues with such a formulation notwithstanding, the Modified Baryonic Dynamics model fails to produce the correct amplitude for the galaxy cluster mass function for any reasonable value of the primordial power spectrum normalisation.

Both major galaxies in the Local Group host planar distributions of co-orbiting satellite galaxies, the Vast Polar Structure (VPOS) of the Milky Way and the Great Plane of Andromeda (GPoA). The $\Lambda$CDM cosmological model did not predict these features. However, according to three recent studies the properties of the GPoA and the flattening of the VPOS are common features among sub-halo based $\Lambda$CDM satellite systems, and the GPoA can be naturally explained by satellites being acquired along cold gas streams. We point out some methodological issues in these studies: either the selection of model satellites is different from that of the observed ones, or an incomplete set of observational constraints has been considered, or the observed satellite distribution is inconsistent with basic assumptions. Once these issues have been addressed, the conclusions are different: features like the VPOS and GPoA are very rare (each with probability $\lesssim 10^{-3}$, and combined probability $< 10^{-5}$) if satellites are selected from a $\Lambda$CDM simulation combined with semi-analytic modelling, and accretion along cold streams is no natural explanation of the GPoA. The origin of planar dwarf galaxy structures remains unexplained in the standard paradigm of galaxy formation.

Since its first formulation in 1983, Milgromian dynamics (MOND) has been very successful in predicting the gravitational potential of galaxies from the distribution of baryons alone, including general scaling relations and detailed rotation curves of large statistical samples of individual galaxies covering a large range of masses and sizes. Most predictions however rely on static models, and only a handful of N-body codes have been developed over the years to investigate the consequences of the Milgromian framework for the dynamics of complex evolving dynamical systems. In this work, we present a new Milgromian N-body code, which is a customized version of the RAMSES code (Teyssier 2002) and thus comes with all its features: it includes particles and gas dynamics, and importantly allows for high spatial resolution of complex systems due to the adaptive mesh refinement (AMR) technique. It further allows the direct comparison between Milgromian simulations and standard Newtonian simulations with dark matter particles. We provide basic tests of this customized code and demonstrate its performance by presenting N-body computations of dark-matter-free spherical equilibrium models as well as dark-matter-free disk galaxies in Milgromian dynamics.

Prompted by the recent successful predictions of the internal dynamics of Andromeda's satellite galaxies (McGaugh & Milgrom 2013a,b), we revisit the classical Milky Way dwarf spheroidal satellites Draco, Sculptor, Sextans, Carina, and Fornax in the framework of Milgromian dynamics (MOND). We use for the first time a Poisson solver with adaptive mesh refinement (AMR) in order to account simultaneously for the gravitational influence of the Milky Way and its satellites. This allows to rigorously model the important external field effect (EFE) of Milgromian dynamics, which can reduce the effective acceleration significantly. We make predictions on the dynamical mass-to-light ratio (Mdyn/L) expected to be measured by an observer who assumes Newtonian dynamics to be valid. We show that Milgromian dynamics predicts typical Mdyn/L ~ 10-50 Msun/Lsun. The results for the most luminous ones, Fornax and Sculptor, agree well with available velocity dispersion data. Moreover, the central power-law slopes of the dynamical masses agree exceedingly well with values inferred observationally from velocity dispersion measurements. The results for Sextans, Carina and Draco are low compared to usually quoted observational estimates, as already pointed out by Angus (2008). For Milgromian dynamics to survive further observational tests in these objects, one would thus need that either (a) previous observational findings based on velocity dispersion measurements have overestimated the dynamical mass due to, e.g., binaries and contaminant outliers, (b) the satellites are not in virial equilibrium due to the Milky Way tidal field, or (c) the specific theory used here does not describe the EFE correctly (e.g., the EFE could be practically negligible in some other theories), or a combination of (a)-(c).

The classical dwarf spheroidals (dSphs) provide a critical test for Modified Newtonian Dynamics (MOND) because they are observable satellite galactic systems with low internal accelerations and low, but periodically varying, external acceleration. This varying external gravitational field is not commonly found acting on systems with low internal acceleration. Using Jeans modelling, Carina in particular has been demonstrated to require a V-band mass-to-light ratio greater than 5, which is the nominal upper limit for an ancient stellar population. We run MOND N-body simulations of a Carina-like dSph orbiting the Milky Way to test if dSphs in MOND are stable to tidal forces over the Hubble time and if those same tidal forces artificially inflate their velocity dispersions and therefore their apparent mass-to-light ratio. We run many simulations with various initial total masses for Carina, and Galactocentric orbits (consistent with proper motions), and compare the simulation line of sight velocity dispersions (losVDs) with the observed losVDs of Walker et al. (2007). We find that the dSphs are stable, but that the tidal forces are not conducive to artificially inflating the losVDs. Furthermore, the range of mass-to-light ratios that best reproduces the observed line of sight velocity dispersions of Carina is 5.3 to 5.7 and circular orbits are preferred to plunging orbits. Therefore, some tension still exists between the required mass-to-light ratio for the Carina dSph in MOND and those expected from stellar population synthesis models. It remains to be seen whether a careful treatment of the binary population or triaxiality might reduce this tension.

In an equilibrium axisymmetric galactic disc, the mean galactocentric radial and vertical velocities are expected to be zero everywhere. In recent years, various large spectroscopic surveys have however shown that stars of the Milky Way disc exhibit non-zero mean velocities outside of the Galactic plane in both the Galactocentric radial and vertical velocity components. While radial velocity structures are commonly assumed to be associated with non-axisymmetric components of the potential such as spiral arms or bars, non-zero vertical velocity structures are usually attributed to excitations by external sources such as a passing satellite galaxy or a small dark matter substructure crossing the Galactic disc. Here, we use a three-dimensional test-particle simulation to show that the global stellar response to a spiral perturbation induces both a radial velocity flow and non-zero vertical motions. The resulting structure of the mean velocity field is qualitatively similar to what is observed across the Milky Way disc. We show that such a pattern also naturally emerges from an analytic toy model based on linearized Euler equations. We conclude that an external perturbation of the disc might not be a requirement to explain all of the observed structures in the vertical velocity of stars across the Galactic disc. Non-axisymmetric internal perturbations can also be the source of the observed mean velocity patterns.

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.

Non-axisymmetries in the Galactic potential (spiral arms and bar) induce kinematic groups such as the Hercules stream. Assuming that Hercules is caused by the effects of the Outer Lindblad Resonance of the Galactic bar, we model analytically its properties as a function of position in the Galaxy and its dependence on the bar's pattern speed and orientation. Using data from the RAVE survey we find that the azimuthal velocity of the Hercules structure decreases as a function of Galactocentric radius, in a manner consistent with our analytical model. This allows us to obtain new estimates of the parameters of the Milky Way's bar. The combined likelihood function of the bar's pattern speed and angle has its maximum for a pattern speed of $\Omega_b=(1.89 \pm 0.08) \times \Omega_0$ where $\Omega_0$ is the local circular frequency. Assuming a Solar radius of 8.05 kpc and a local circular velocity of 238 km/s, this corresponds to $\Omega_b=56\pm 2$ km/s/kpc. On the other hand, the bar's orientation $\phi_b$ cannot be constrained with the available data. In fact, the likelihood function shows that a tight correlation exists between the pattern speed and the orientation, implying that a better description of our best fit results is given by the linear relation $\Omega_b/\Omega_0=1.91 + 0.0044[\phi_b(\deg)-48]$, with standard deviation of 0.02. For example, for an angle of $\phi_b=30\deg$ the pattern speed is $54.0\pm0.5$ km/s/kpc. These results are not very sensitive to the other Galactic parameters such as the circular velocity curve or the peculiar motion of the Sun, and are robust to biases in distance.

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.

The Lambda-CDM cosmological model is succesful at reproducing various independent sets of observations concerning the large-scale Universe. This model is however currently, and actually in principle, unable to predict the gravitational field of a galaxy from it observed baryons alone. Indeed the gravitational field should depend on the relative contribution of the particle dark matter distribution to the baryonic one, itself depending on the individual assembly history and environment of the galaxy, including a lot of complex feedback mechanisms. However, for the last thirty years, Milgrom's formula, at the heart of the MOND paradigm, has been consistently succesful at predicting rotation curves from baryons alone, and has been resilient to all sorts of observational tests on galaxy scales. We show that the few individual galaxy rotation curves that have been claimed to be highly problematic for the predictions of Milgrom's formula, such as Holmberg II or NGC 3109, are actually false alarms. We argue that the fact that it is actually possible to predict the gravitational field of galaxies from baryons alone presents a challenge to the current Lambda-CDM model, and may indicate a breakdown of our understanding of gravitation and dynamics, and/or that the actual lagrangian of the dark sector is very different and richer than currently assumed. On the other hand, it is obvious that any alternative must also, in fine, reproduce the successes of the Lambda-CDM model on large scales, where this model is so well-tested that it presents by itself a challenge to any such alternative.

We use our Modified Newtonian Dynamics (MOND) cosmological particle-mesh N-body code to investigate the feasibility of structure formation in a framework involving MOND and light sterile neutrinos in the mass range 11 - 300 eV: always assuming that \Omega_{\nu_s}=0.225 for H_o=72 \kms Mpc^{-1}. We run a suite of simulations with variants on the expansion history, cosmological variation of the MOND acceleration constant, different normalisations of the power spectrum of the initial perturbations and interpolating functions. Using various box sizes, but typically with ones of length 256 Mpc/h, we compare our simulated halo mass functions with observed cluster mass functions and show that (i) the sterile neutrino mass must be larger than 30 eV to account for the low mass (M_{200}<10^{14.6} solar masses) clusters of galaxies in MOND and (ii) regardless of sterile neutrino mass or any of the variations we mentioned above, it is not possible to form the correct number of high mass (M_{200}>10^{15.1} solar masses) clusters of galaxies: there is always a considerable over production. This means that the ansatz of considering the weak-field limit of MOND together with a component of light sterile neutrinos to form structure from z ~ 200 fails. If MOND is the correct description of weak-field gravitational dynamics, it could mean that subtle effects of the additional fields in covariant theories of MOND render the ansatz inaccurate, or that the gravity generated by light sterile neutrinos (or by similar hot dark matter particles) is different from that generated by the baryons.

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.

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.

We analyse the kinematics of disc stars observed by the RAVE survey in and beyond the Solar neighbourhood.We detect significant overdensities in the velocity distributions using a technique based on the wavelet transform.We find that the main local kinematic groups are large scale features, surviving at least up to ~1 kpc from the Sun in the direction of anti-rotation, and also at ~700 pc below the Galactic plane.We also find that for regions located at different radii than the Sun, the known groups appear shifted in the velocity plane. For example, the Hercules group has a larger azimuthal velocity for regions inside the Solar circle and a lower value outside. We have also discovered a new group at (U, V) = (92,-22) km/s in the Solar neighbourhood and confirmed the significance of other previously found groups. Some of these trends detected for the first time are consistent with dynamical models of the effects of the bar and the spiral arms. More modelling is required to definitively characterise the non-axisymmetric components of our Galaxy using these groups.

The Local Group timing has been one of the first historical probes of the missing mass problem. Whilst modern cosmological probes indicate that pure baryonic dynamics is not sufficient on the largest scales, nearby galaxies and small galaxy groups persistently obey Milgrom's MOND law, which implies that dynamics at small scales is possibly entirely predicted by the baryons. Here, we investigate the Local Group timing in this context of Milgromian dynamics. Making use of the latest measured proper motions and radial velocities for Andromeda and the Magellanic clouds, we integrate backwards their orbits by making use of the Milgromian two-body equation of motion. We find that, with the currently measured proper motions and radial velocity of M31, MOND would imply that the Milky Way and M31 first moved apart via Hubble expansion after birth, but then necessarily got attracted again by the Milgromian gravitational attraction, and had a past fly-by encounter before coming to their present positions. This encounter would most probably have happened 7 to 11 Gyr ago (0.8<z<3). The absence of a dark matter halo and its associated dynamical friction is necessary for such a close encounter not to trigger a merger. Observational arguments which could exclude or favour such a past encounter would thus be very important in view of falsifying or vindicating Milgromian dynamics on the scale of the Local Group. Interestingly, the closest approach of the encounter is small enough (<55 kpc) to have severe consequences on the disk dynamics, including perhaps thick disk formation, and on the satellite systems of both galaxies. The ages of the satellite galaxies and of the young halo globular clusters, all of which form the vast polar structure around the Milky Way, are consistent with these objects having been born in this encounter.

Polar ring galaxies are ideal objects with which to study the three-dimensional shapes of galactic gravitational potentials since two rotation curves can be measured in two perpendicular planes. Observational studies have uncovered systematically larger rotation velocities in the extended polar rings than in the associated host galaxies. In the dark matter context, this can only be explained through dark halos that are systematically flattened along the polar rings. Here, we point out that these objects can also be used as very effective tests of gravity theories, such as those based on Milgromian dynamics (MOND). We run a set of polar ring models using both Milgromian and Newtonian dynamics to predict the expected shapes of the rotation curves in both planes, varying the total mass of the system, the mass of the ring with respect to the host, as well as the size of the hole at the center of the ring. We find that Milgromian dynamics not only naturally leads to rotation velocities being typically higher in the extended polar rings than in the hosts, as would be the case in Newtonian dynamics without dark matter, but that it also gets the shape and amplitude of velocities correct. Milgromian dynamics thus adequately explains this particular property of polar ring galaxies.

The Universe on large scales is well described by the Lambda-CDM cosmological model. There however remain some heavy clouds on our global understanding, especially on galaxy scales, which we review here. While some of these clouds might perhaps disappear through small compensatory adjustments of the model, such as changing the mass of the dark matter particles or accounting better for baryonic physics, others should rather be taken as strong indications that the physics of the dark sector is, at the very least, much richer and complex than currently assumed, and that our understanding of gravity and dynamics might also be at play. For instance, the empirically well-tested MOND phenomenology in galaxies, whatever its final explanation, should be understood in any model of galaxy formation and dynamics. Current alternatives to Lambda-CDM however bring with them many unsolved questions and challenges.

Non-axisymmetric components, such as spirals and central bars, play a major role in shaping galactic discs. An important aspect of the disc secular evolution driven by these perturbers is the radial migration of stars. It has been suggested recently that migration can populate a thick-disc component from inner-disc stars with high vertical energies. Since this has never been demonstrated in simulations, we study in detail the effect of radial migration on the disc velocity dispersion and disc thickness, by separating simulated stars into migrators and non-migrators. We apply this method to three isolated barred Tree-SPH N-body galaxies with strong radial migration. Contrary to expectations, we find that as stellar samples migrate, on the average, their velocity dispersion change (by as much as 50%) in such a way as to approximately match the non-migrating population at the radius at which they arrive. We show that, in fact, migrators suppress heating in parts of the disc. To confirm the validity of our findings, we also apply our technique to three cosmological re-simulations, which use a completely different simulation scheme and, remarkably, find very similar results. We believe the inability of migration to thicken discs is a fundamental property of internal disc evolution, irrespective of the migration mechanism at work. We explain this with the approximate conservation of the (average) vertical and radial actions rather than the energy. This "action mixing" can be used to constrain the migration rate in the Milky Way: estimates of the average vertical action in observations for different populations of stars should reveal flattening with radius for older groups of stars.

We present a new method to derive orbital elements of double-lined spectroscopic binaries (SB2). The aim is to have accurate orbital parameters of a selection of SB2 in order to prepare the exploitation of astrometric Gaia observations. Combined with our results, they should allow one to measure the mass of each star with a precision of better than 1%. The new method presented here consists of using the spectra at all epochs simultaneously to derive the orbital elements without templates. It is based on a Markov chain including a new method for disentangling the spectra.

Gaia is an ambitious ESA space mission which will provide photometric and astrometric measurements with the accuracies needed to produce a kinematic census of almost one billion stars in our Galaxy. These data will revolutionize our understanding of the dynamics of the Milky Way, and our knowledge of its detailed gravitational potential and mass distribution, including the putative dark matter component and the non-axisymmetric features such as spiral arms. The Gaia mission will help to answer various currently unsettled questions by using kinematic information on both disk and halo stellar populations. Among many others: what does the rotation curve of the outer Galaxy look like? How far from axisymmetry and equilibrium is the Galaxy? What are the respective roles of hierarchical formation and secular evolution in shaping the Galaxy and its various components? Are the properties of the Galaxy in accordance with expectations from the standard model of cosmology?

We investigate the evolution of galactic disks in N-body Tree-SPH simulations. We find that disks, initially truncated at three scale-lengths, can triple their radial extent, solely driven by secular evolution. Both Type I (single exponential) and Type II (down-turning) observed disk surface-brightness profiles can be explained by our findings. We relate these results to the strong angular momentum outward transfer, resulting from torques and radial migration associated with multiple patterns, such as central bars and spiral waves of different multiplicity. We show that even for stars ending up on cold orbits, the changes in angular momentum exhibit complex structure as a function of radius, unlike the expected effect of transient spirals alone. Focussing on one of our models, we find evidence for non-linear coupling among m=1, 2, 3 and 4 density waves, where m is the pattern multiplicity. We suggest that the naturally occurring larger resonance widths at galactic radii beyond four scale-lengths may have profound consequences on the formation and location of breaks in disk density profiles, provided spirals are present at such large distances. We also consider the effect of gas inflow and show that when in-plane smooth gas accretion of ~5 M_sun/yr is included, the outer disks become more unstable, leading to a strong increase in the stellar velocity dispersion. This, in turn, causes the formation of a Type III (up-turning) profile in the old stellar population. We propose that observations of Type III surface brightness profiles, combined with an up-turn in the stellar velocity dispersions beyond the disk break, could be a signature of ongoing gas-accretion. The results of this study suggest that disk outskirts comprised of stars migrated from the inner disk would have relatively large radial velocity dispersions, and significant thickness when seen edge-on. [Abridged]

A remarkably tight relation is observed between the Newtonian gravity sourced by the baryons and the actual gravity in galaxies of all sizes. This can be interpreted as the effect of a single effective force law depending on acceleration. This is however not the case in larger systems with much deeper potential wells, such as galaxy clusters. Here we explore the possibility of an effective force law reproducing mass discrepancies in all extragalactic systems when depending on both acceleration and the depth of the potential well. We exhibit, at least at a phenomenological level, one such possible construction in the classical gravitational potential theory. Interestingly, the framework, dubbed EMOND, is able to reproduce the observed mass discrepancies in both galaxies and galaxy clusters, and to produce multi-center systems with offsets between the peaks of gravity and the peaks of the baryonic distribution.

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.

Within the framework of modified Newtonian dynamics (MOND) we investigate the kinematics of two dwarf spiral galaxies belonging to very different environments, namely KK 246 in the Local Void and Holmberg II in the M81 group. A mass model of the rotation curve of KK 246 is presented for the first time, and we show that its observed kinematics are consistent with MOND. We re-derive the outer rotation curve of Holmberg II, by modelling its HI data cube, and find that its inclination should be closer to face-on than previously derived. This implies that Holmberg II has a higher rotation velocity in its outer parts, which, although not very precisely constrained, is consistent with the MOND prediction.

A wealth of astronomical data indicate the presence of mass discrepancies in the Universe. The motions observed in a variety of classes of extragalactic systems exceed what can be explained by the mass visible in stars and gas. Either (i) there is a vast amount of unseen mass in some novel form - dark matter - or (ii) the data indicate a breakdown of our understanding of dynamics on the relevant scales, or (iii) both. Here, we first review a few outstanding challenges for the dark matter interpretation of mass discrepancies in galaxies, purely based on observations and independently of any alternative theoretical framework. We then show that many of these puzzling observations are predicted by one single relation - Milgrom's law - involving an acceleration constant (or a characteristic surface density) of the order of the square-root of the cosmological constant in natural units. This relation can at present most easily be interpreted as the effect of a single universal force law resulting from a modification of Newtonian dynamics (MOND) on galactic scales. We exhaustively review the current observational successes and problems of this alternative paradigm at all astrophysical scales, and summarize the various theoretical attempts (TeVeS, GEA, BIMOND, and others) made to effectively embed this modification of Newtonian dynamics within a relativistic theory of gravity.

Several scenarios have been suggested to explain the phase-space distribution of the Milky Way (MW) satellite galaxies in a disc of satellites (DoS). To quantitatively compare these different possibilities, a new method analysing angular momentum directions in modelled data is presented. It determines how likely it is to find sets of angular momenta as concentrated and as close to a polar orientation as is observed for the MW satellite orbital poles. The method can be easily applied to orbital pole data from different models. The observed distribution of satellite orbital poles is compared to published angular momentum directions of subhalos derived from six cosmological state-of-the-art simulations in the Aquarius project. This tests the possibility that filamentary accretion might be able to naturally explain the satellite orbits within the DoS. For the most likely alignment of main halo and MW disc spin, the probability to reproduce the MW satellite orbital pole properties turns out to be less than 0.5 per cent in Aquarius models. Even an isotropic distribution of angular momenta has a higher likelihood to produce the observed distribution. The two Via Lactea cosmological simulations give results similar to the Aquarius simulations. Comparing instead with numerical models of galaxy-interactions gives a probability of up to 90 per cent for some models to draw the observed distribution of orbital poles from the angular momenta of tidal debris. This indicates that the formation as tidal dwarf galaxies in a single encounter is a viable, if not the only, process to explain the phase-space distribution of the MW satellite galaxies.

Here we present a new particle-mesh galactic N-body code that uses the full multigrid algorithm for solving the modified Poisson equation of the Quasi Linear formulation of Modified Newtonian Dynamics (QUMOND). A novel approach for handling the boundary conditions using a refinement strategy is implemented and the accuracy of the code is compared with analytical solutions of Kuzmin disks. We then employ the code to compute the predicted rotation curves for a sample of five spiral galaxies from the THINGS sample. We generated static N-body realisations of the galaxies according to their stellar and gaseous surface densities and allowed their distances, mass-to-light ratios (M/L) and both the stellar and gas scale-heights to vary in order to estimate the best fit parameters. We found that NGC 3621, NGC 3521 and DDO 154 are well fit by MOND using expected values of the distance and M/L. NGC 2403 required a moderately larger $M/L$ than expected and NGC 2903 required a substantially larger value. The surprising result was that the scale-height of the dominant baryonic component was well constrained by the rotation curves: the gas scale-height for DDO 154 and the stellar scale-height for the others. In fact, if the suggested stellar scale-height (one-fifth the stellar scale-length) was used in the case of NGC 3621 and NGC 3521 it would not be possible to produce a good fit to the inner rotation curve. For each of the four stellar dominated galaxies, we calculated the vertical velocity dispersions which we found to be, on the whole, quite typical compared with observed stellar vertical velocity dispersions of face on spirals. We conclude that modelling the gas scale-heights of the gas rich dwarf spiral galaxies will be vital in order to make precise conclusions about MOND.

Disc non-axisymmetrc components, such as spirals and central bars, are nowadays known to play an important role in shaping galactic discs. Here we use Tree-SPH N-body simulations to examine the effect of these perturbers on two aspects: the occurrence of multiple patterns in discs and the effects of radial migration on disc thickening. We find that, in addition to a central bar, multiple spiral patterns and lopsided modes develop in all models. Interaction among these asymmetric features results in a large scale stellar migration. However, we show that, despite the strong radial mixing, discs cannot be thickened sufficiently to match observed thick discs. We relate this to the adiabatic cooling as stars migrate radially outwards. We also find that the bulge contribution to a thick-disc component for an Sa-type galaxy at ~2.5 disc scale-lengths is less than 1% and zero in the case of a Milky Way-like, Sb-type. Our findings cast doubt on the plausibility of thick disc formation via stellar radial migration.

Bars and spirals are among the main drivers of the secular evolution of galactic disks, and it is, thus, of prime importance to better understand their exact nature and dynamics. In this respect, the Milky Way provides a unique laboratory in which to study a snapshot of their detailed dynamical influence in six-dimensional phase-space. With the advent of present and future astrometric and spectroscopic surveys, such six-dimensional phase-space data can be obtained for stars in an increasingly large volume around the Sun. Here, we review the signatures of - and constraints on - disk non-axisymmetries obtained from recent stellar kinematical data sets, including (i) the detection of resonant moving groups in the solar neighbourhood, (ii) the non-zero value of the Oort constants C and K, and (iii) the detection of a radial velocity gradient of 4 km/s/kpc in the extended solar neighbourhood (d<2 kpc).

The connection between some features of the metallicity gradient in the Galactic disk, best revealed by Open Clusters and Cepheids, and the spiral structure, is explored. The step-like abrupt decrease in metallicity at 8.5 kpc (with R_0= 7.5 kpc, or at 9.5 kpc if R_0 = 8.5 kpc is adopted) is well explained by the corotation ring-shaped gap in the density of gas, which isolates the internal and external regions of the disk one from the other. This solves a long standing problem of understanding the different chemical characteristics of the inner and outer parts of the disk. The time required to build up the metallicity difference between the two sides of the step is a measure of the minimal life-time of the present grand-design spiral pattern structure, of the order of 3 Gyr. The plateaux observed on each side of the step are interpreted in terms of the large scale radial motion of the stars and of the gas flow induced by the spiral structure. The star-formation rate revealed by the density of open clusters is maximum in the Galactic radial range from 6 to 12 kpc (with an exception of a narrow gap at corotation), coinciding with the region where the 4-arms mode is allowed to exist. We argue that most of the old open clusters situated at large galactocentric radii were born in this inner region where conditions more favorable to star-formation are found. The ratio of $\alpha$-elements to Fe of the sample of Cepheids does not vary appreciably with the Galactic radius, which reveals an homogeneous history of star formation. Different arguments are given showing that usual approximations of chemical evolution models, which assume fast mixing of metallicity in the azimuthal direction and ignore the existence of the spiral arms, are a poor ones.

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.

Two incorrect arguments against MOND in elliptical galaxies could be that the equivalent circular velocity curves tend to become flat at much larger accelerations than in spiral galaxies, and that the Newtonian dark matter halos are more concentrated than in spirals. Here, we compare published scaling relations for the dark halos of elliptical galaxies to the scaling relations expected for MONDian phantom halos. We represent the baryonic content of galaxies by spherical profiles, and their corresponding MONDian phantom halos by logarithmic halos. We then derive the surface densities, central densities, and phase space densities and compare them with published scaling relations. We conclude that it is possible to get flat circular velocity curves at high acceleration in MOND, and that this happens for baryonic distributions described by Jaffe profiles in the region where the circular velocity curve is flat. Moreover, the scaling relations of dark halos of ellipticals are remarkably similar to the scaling relations of phantom halos of MOND.

The Hyades stream has long been thought to be a dispersed vestige of the Hyades cluster. However, recent analyses of the parallax distribution, of the mass function, and of the action-space distribution of stream stars have shown it to be rather composed of orbits trapped at a resonance of a density disturbance. This resonant scenario should leave a clearly different signature in the element abundances of stream stars than the dispersed cluster scenario, since the Hyades cluster is chemically homogeneous. Here, we study the metalllicty as well as the element abundances of Li, Na, Mg, Fe, Zr, Ba, La, Ce, Nd, and Eu for a random sample of stars belonging to the Hyades stream, and compare them with those of stars from the Hyades cluster. From this analysis: (i) we independently confirm that the Hyades stream cannot be solely composed of stars originating in the Hyades cluster; (ii) we show that some stars from the Hyades stream nevertheless have abundances compatible with an origin in the cluster; (iii) we emphasize that the use of Li as a chemical tag of the cluster origin of main-sequence stars is very efficient in the range 5500 - 6200 K, since the Li sequence in the Hyades cluster is very tight; (iv) we show that, while this evaporated population has a metallicity excess of ~0.2 dex w.r.t. the local thin disk population, identical to that of the Hyades cluster, the remainder of the Hyades stream population has still a metallicity excess of ~0.06 to 0.15 dex, consistent with an origin in the inner Galaxy; (v) we show that the Hyades stream can be interpreted as an inner 4:1 resonance of the spiral pattern: this then also reproduces an orbital family compatible with the Sirius stream, and places the origin of the Hyades stream up to 1 kpc inwards from the solar radius, which might explain the observed metallicity excess of the stream population.

We have recently identified a new radial migration mechanism resulting from the overlap of spiral and bar resonances in galactic disks. Here we confirm the efficiency of this mechanism in fully self-consistent, Tree-SPH simulations, as well as high-resolution pure N-body simulations. In all barred cases we clearly identify the effect of spiral-bar resonance overlap by measuring a bimodality in the changes of angular momentum in the disk, dL, whose maxima are near the bar's corotation and outer Lindblad resonance. This contrasts with the smooth distribution of dL for a simulation with no stable bar present, where strong radial migration is induced by multiple spirals. The presence of a disk gaseous component appears to increase the rate of angular momentum exchange by about 20%. The efficiency of this mechanism is such that galactic stellar disks can extend to over 10 scale-lengths within 1-3 Gyr in both Milky Way size and low-mass galaxies (circular velocity ~100 km/s). We also show that metallicity gradients can flatten in less than 1 Gyr rendering mixing in barred galaxies an order of magnitude more efficient than previously thought.

We present the analysis of 12 high-resolution galactic rotation curves from The HI Nearby Galaxy Survey (THINGS) in the context of modified Newtonian dynamics (MOND). These rotation curves were selected to be the most reliable for mass modelling, and they are the highest quality rotation curves currently available for a sample of galaxies spanning a wide range of luminosities. We fit the rotation curves with the "simple" and "standard" interpolating functions of MOND, and we find that the "simple" function yields better results. We also redetermine the value of a0, and find a median value very close to the one determined in previous studies, a0 = (1.22 +- 0.33) x 10^{-8} cm/s^2. Leaving the distance as a free parameter within the uncertainty of its best independently determined value leads to excellent quality fits for 75% of the sample. Among the three exceptions, two are also known to give relatively poor fits also in Newtonian dynamics plus dark matter. The remaining case (NGC 3198), presents some tension between the observations and the MOND fit, which might however be explained by the presence of non-circular motions, by a small distance, or by a value of a0 at the lower end of our best-fit interval, 0.9 x 10^{-8} cm/s^2. The best-fit stellar M/L ratios are generally in remarkable agreement with the predictions of stellar population synthesis models. We also show that the narrow range of gravitational accelerations found to be generated by dark matter in galaxies is consistent with the narrow range of additional gravity predicted by MOND.

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.

(abridged) Predictions of the Concordance Cosmological Model (CCM) of the structures in the environment of large spiral galaxies are compared with observed properties of Local Group galaxies. Five new most probably irreconcilable problems are uncovered. However, the Local Group properties provide hints that may lead to a solution of the above problems The DoS and bulge--satellite correlation suggest that dissipational events forming bulges are related to the processes forming phase-space correlated satellite populations. Such events are well known to occur since in galaxy encounters energy and angular momentum are expelled in the form of tidal tails, which can fragment to form populations of tidal-dwarf galaxies (TDGs) and associated star clusters. If Local Group satellite galaxies are to be interpreted as TDGs then the sub-structure predictions of CCM are internally in conflict. All findings thus suggest that the CCM does not account for the Local Group observations and that therefore existing as well as new viable alternatives have to be further explored. These are discussed and natural solutions for the above problems emerge.

We have recently identified a previously unknown radial migration mechanism resulting from the overlap of spiral and bar resonances in galactic discs (Minchev & Famaey 2010, Minchev et al. 2010). This new mechanism is much more efficient than mixing by transient spirals and its presence is unavoidable in all barred galaxies, such as our own Milky Way. The consequences of this are a strong flattening in the metallicity gradient in the disc, an extended disc profile, and the formation of a thick disc component, all taking place in only a couple of Gyr. This timescale is drastically shorter than previously expected and thus can put strong constraints on the longevity, strength and pattern speeds of the Galactic bar and Spiral Structure.

We have chosen the name of GYES, one of the mythological giants with one hundred arms, offspring of Gaia and Uranus, for our instrument study of a multifibre spectrograph for the prime focus of the Canada-France-Hawaii Telescope. Such an instrument could provide an excellent ground-based complement for the Gaia mission and a northern complement to the HERMES project on the AAT. The CFHT is well known for providing a stable prime focus environment, with a large field of view, which has hosted several imaging instruments, but has never hosted a multifibre spectrograph. Building upon the experience gained at GEPI with FLAMES-Giraffe and X-Shooter, we are investigating the feasibility of a high multiplex spectrograph (about 500 fibres) over a field of view 1 degree in diameter. We are investigating an instrument with resolution in the range 15000 to 30000, which should provide accurate chemical abundances for stars down to 16th magnitude and radial velocities, accurate to 1 km/s for fainter stars. The study is led by GEPI-Observatoire de Paris with a contribution from Oxford for the study of the positioner. The financing for the study comes from INSU CSAA and Observatoire de Paris. The conceptual study will be delivered to CFHT for review by October 1st 2010.

While it has long been known that a large number of short-lived transient spirals can cause stellar migration, here we report that another mechanism is also effective at mixing disks of barred galaxies. The resonance overlap of the bar and spiral structure induces a nonlinear response leading to a strong redistribution of angular momentum in the disk. We find that, depending on the amplitudes of the perturbers, the changes in angular momentum, dL, could occur up to an order of magnitude faster than in the case of recurrent spirals. The signature of this mechanism is a bimodality in dL with maxima near the bar's corotation and its outer Lindblad resonance; this is independent of the properties of the spiral structure. For parameters consistent with the Milky Way the disk mixes in about 3 Gyr and the stellar velocity dispersion increases with time in a manner roughly consistent with observations. This new mechanism could account for both the observed age-velocity relation and the absence of age-metallicity relation in the solar neighborhood. Spiral-bar interaction could also explain observations showing that strongly barred galaxies have weaker metallicity gradients than weakly barred or non-barred galaxies.

We point out an interesting theoretical prediction for elliptical galaxies residing inside galaxy clusters in the framework of modified Newtonian dynamics (MOND), that could be used to test this paradigm. Apart from the central brightest cluster galaxy, other galaxies close enough to the centre experience a strong gravitational influence from the other galaxies of the cluster. This influence manifests itself only as tides in standard Newtonian gravity, meaning that the systematic acceleration of the centre of mass of the galaxy has no consequence. However, in the context of MOND, a consequence of the breaking of the strong equivalence principle is that the systematic acceleration changes the own self-gravity of the galaxy. We show here that, in this framework, initially axisymmetric elliptical galaxies become lopsided along the external field's direction, and that the centroid of the galaxy, defined by the outer density contours, is shifted by a few hundreds parsecs with respect to the densest point.

We present HI observations performed at the GMRT of the nearby dwarf galaxy NGC 1560. This Sd galaxy is well-known for a distinct "wiggle" in its rotation curve. Our new observations have twice the resolution of the previously published HI data. We derived the rotation curve by taking projection effects into account, and we verified the derived kinematics by creating model datacubes. This new rotation curve is similar to the previously published one: we confirm the presence of a clear wiggle. The main differences are in the innermost ~100 arcsec of the rotation curve, where we find slightly (<~ 5 km/s) higher velocities. Mass modelling of the rotation curve results in good fits using the core-dominated Burkert halo (which however does not reproduce the wiggle), bad fits using the a Navarro, Frenk & White halo, and good fits using MOND (Modified Newtonian Dynamics), which also reproduces the wiggle.

There exists several modified gravity theories designed to reproduce the empirical Milgrom's formula (MOND). Here we derive analytical results in the context of the static weak-field limit of two of them (BIMOND, leading for a given set of parameters to QUMOND, and TeVeS). In this limit, these theories are constructed to give the same force field for spherical symmetry, but their predictions generally differ out of it. However, for certain realizations of these theories (characterized by specific choices for their free functions), the binding potential-energy of a system is increased, compared to its Newtonian counterpart, by a constant amount independent of the shape and size of the system. In that case, the virial theorem is exactly the same in these two theories, for the whole gravity regime and even outside of spherical symmetry, although the exact force fields are different. We explicitly show this for the force field generated by the two theories inside an elliptical shell. For more general free functions, the virial theorems are however not identical in these two theories. We finally explore the consequences of these analytical results for the two-body force.

Certain types of globular clusters have the very important property that the predictions for their kinematics in the Newtonian and modified Newtonian dynamics (MOND) contexts are divergent. Here, we caution the recent claim that the stellar kinematics data (using 17 stars) of the globular cluster Palomar 14 are inconsistent with MOND. We compare the observations to the theoretical predictions using a Kolmogorov-Smirnov test, which is appropriate for small samples. We find that, with the currently available data, the MOND prediction for the velocity distribution can only be excluded with a very low confidence level, clearly insufficient to claim that MOND is falsified.

It was recently discovered that the mean dark matter surface density within one dark halo scale length - the radius within which the volume density profile of dark matter remains approximately flat - is constant across a wide range of galaxies. This scaling relation holds for galaxies spanning a luminosity range of 14 magnitudes and the whole Hubble sequence. Here we report that the luminous matter surface density is also constant within one scale length of the dark halo. This means that the gravitational acceleration generated by the luminous component in galaxies is always the same at this radius. Although the total luminous-to-dark matter ratio is not constant, within one halo scale length it is constant. Our finding can be interpreted as a close correlation between the enclosed surface densities of luminous and dark matter in galaxies.

(Abridged) Modified Newtonian Dynamics (MOND) can fit a broad range of galaxy kinematic data, but struggles with clusters of galaxies. MONDian clusters need dark matter, and here we test the 11 eV sterile neutrino - used to fit the first three acoustic peaks of the cosmic microwave background - by investigating their equilibrium distributions in 30 groups and clusters over a wide range of temperatures. Knowing both the sterile neutrino density and velocity dispersion tells us the Tremaine-Gunn phase-space limit at all radii. We find that all 30 systems serendipitously reach the Tremaine-Gunn limit by the centre, which means a portion of the dynamical mass must always be covered by the brightest cluster galaxy. Interestingly, the typical fitted K-band mass-to-light ratio is unity and at most 1.2, which is very consistent - although leaving no margin for error - with stellar population synthesis models. Amidst the sample there are several special cases including the Coma cluster (for which dark matter was first proposed), NGC 720 (where geometrical evidence for dark matter was found) and the bullet cluster (where dark matter - of some kind - in clusters was directly proven to exist). We demonstrate that 11 eV sterile neutrinos are unlikely to influence spiral galaxy rotation curves, as they don't influence even some very massive early-types (NGC 4125 and NGC 6482). Finally, we conclude that it is intriguing that the minimum mass of sterile neutrino particle that can match the cosmic microwave background is identical to the minimum mass found here to be consistent with equilibrium configurations of MONDian clusters of galaxies.

After summarizing the respective merits of the Cold Dark Matter (CDM) and Modified Newtonian Dynamics (MOND) paradigms in various stellar systems, we investigate the possibility that a non-standard interaction between baryonic and dark matter could reproduce the successes of CDM at extragalactic scales while making baryonic matter effectively obey the MOND field equation in spiral galaxies.

We look for observational signatures that could discriminate between Newtonian and modified Newtonian (MOND) dynamics in the Milky Way, in view of the advent of large astrometric and spectroscopic surveys. Indeed, a typical signature of MOND is an apparent disk of "phantom" dark matter, which is uniquely correlated with the visible disk-density distribution. Due to this phantom dark disk, Newtonian models with a spherical halo have different signatures from MOND models close to the Galactic plane. The models can thus be differentiated by measuring dynamically (within Newtonian dynamics) the disk surface density at the solar radius, the radial mass gradient within the disk, or the velocity ellipsoid tilt angle above the Galactic plane. Using the most realistic possible baryonic mass model for the Milky Way, we predict that, if MOND applies, the local surface density measured by a Newtonist will be approximately 78 Msun/pc2 within 1.1 kpc of the Galactic plane, the dynamically measured disk scale-length will be enhanced by a factor of 1.25 with respect to the visible disk scale-length, and the local vertical tilt of the velocity ellipsoid at 1 kpc above the plane will be approximately 6 degrees. None of these tests can be conclusive for the present-day accuracy of Milky Way data, but they will be of prime interest with the advent of large surveys such as GAIA.

In recent years several hypervelocity stars (HVSs) have been observed in the halo of our Galaxy. Such HVSs have possibly been ejected from the Galactic center and then propagated in the Galactic potential up to their current position. The recent survey for candidate HVSs show an asymmetry in the kinematics of candidate HVSs (position and velocity vectors), where more outgoing stars than ingoing stars (i.e. positive Galactocentric velocities vs. negative ones) are observed. We show that such kinematic asymmetry, which is likely due to the finite lifetime of the stars and Galactic potential structure, could be used in a novel method to probe and constrain the Galactic potential, identify the stellar type of the stars in the survey and estimate the number of HVSs. Kinematics-independent identification of the stellar types of the stars in such surveys (e.g. spectroscopic identification) could further improve these results. We find that the observed asymmetry between ingoing and outgoing stars favors specific Galactic potential models. It also implies a lower limit of ~54+-8 main sequence HVSs in the survey sample (>=648+-96 in the Galaxy), assuming that all of the main sequence stars in the survey originate from the Galactic center. The other stars in the survey are likely to be hot blue horizontal branch stars born in the halo rather than stars ejected from the Galactic center.

This paper is the third one in a series devoted to studying the properties of binaries involving M giants. We use a new set of orbits to construct the first (e-logP) diagram of an extensive sample of M giant binaries, to obtain their mass-function distribution, and to derive evolutionary constraints for this class of binaries and related systems. The orbital properties of binaries involving M giants were analysed and compared with those of related families of binaries (K giants, post-AGB stars, barium stars, Tc-poor S stars). The orbital elements of post-AGB stars and M giants are not different, which may very indicate that, for the considered sample of post-AGB binaries, the post-AGB star left the AGB at quite an early stage (M4 or so). Neither are the orbital elements of post-mass-transfer binaries like barium stars very different from those of M giants, suggesting that the mass transfer did not alter the orbital elements much, contrary to current belief. Finally, we show that binary systems with e < 0.4 log P - 1 (with periods expressed in days) are predominantly post-mass-transfer systems, because (i) the vast majority of barium and S systems match this condition, and (ii) these systems have companion masses peaking around 0.6 solar mass, as expected for white dwarfs. The latter property has been shown to hold as well for open-cluster binaries involving K giants, for which a lower bound on the companion mass may easily be set.

This paper is the second one in a series devoted to the study of properties of binaries involving M giants. The binary frequency of field M giants is derived and compared with the binary fraction of K giants. Diagrams of the CORAVEL spectroscopic parameter Sb (measuring the average line-width) vs. radial-velocity standard deviation for our samples are used to define appropriate binarity criteria. These then serve to extract the binarity fraction among the M giants. Comparison is made to earlier data on K giants binarity frequency. The Sb parameter is discussed in relation to global stellar parameters and the Sb vs. stellar radius relation is used to identify fast rotators. We find that the spectroscopic binary detection rate among field M giants, in a sample with a low number of velocity measurements (~2), unbiased toward earlier known binaries, is 6.3%. This is less than half of the analogous rate for field K giants, likely resulting from a real difference. This difference originates in the greater difficulty of finding binaries among M giants because of their smaller orbital velocity amplitudes and larger intrinsic jitter and in the different distributions of K and M giants in the eccentricity-period diagram. A larger detection rate was obtained in a smaller M giant sample with more radial velocity measurements per object: 11.1% confirmed plus 2.7% possible binaries. The CORAVEL spectroscopic parameter Sb was found to correlate better with the stellar radius than with either luminosity or effective temperature separately. Two outliers of the Sb vs. stellar radius relation, HD 190658 and HD 219654, have been recognized as fast rotators. The rotation is companion-induced, as both objects turn out to be spectroscopic binaries.

This paper is a follow-up of the vast effort to collect radial velocity data for stars belonging to the Hipparcos survey. We aim at extending the orbital data available for binaries with M giant primaries. The data will be used in the companion papers of this series to (i) derive the binary frequency among M giants and compare it to that of K giants (Paper II), and (ii) analyse the eccentricity-period diagram and the mass-function distribution (Paper III). Keplerian solutions are fitted to radial-velocity data. However, for several stars, no satisfactory solution could be found, despite the fact that the radial-velocity standard deviation is larger than the instrumental error, because M giants suffer from intrinsic radial-velocity variations due to pulsations. We show that these intrinsic radial-velocity variations can be linked with both the average spectral-line width and the photometric variability. We present an extensive collection of spectroscopic orbits for M giants, with 12 new orbits, plus 17 from the literature. Moreover, to illustrate the fact that the large radial-velocity jitter present in Mira and semi-regular variables may easily be confused with orbital variations, we also present examples of pseudo-orbital variations (in S UMa, X Cnc and possibly in HD 115521, a former IAU radial-velocity standard). Because of this difficulty, M giants involving Mira variables were excluded from our monitored sample. We finally show that the majority of M giants detected as X-ray sources are actually binaries.

Observations of galaxies suggest a one-to-one analytic relation between the inferred gravity of dark matter at any radius and the enclosed baryonic mass, a relation summarized by Milgrom's law of modified Newtonian dynamics (MOND). However, present-day covariant versions of MOND usually require some additional fields contributing to the geometry, as well as an additional hot dark matter component to explain cluster dynamics and cosmology. Here, we envisage a slightly more mundane explanation, suggesting that dark matter does exist but is the source of MOND-like phenomenology in galaxies. We assume a canonical action for dark matter, but also add an interaction term between baryonic matter, gravity, and dark matter, such that standard matter effectively obeys the MOND field equation in galaxies. We show that even the simplest realization of the framework leads to a model which reproduces some phenomenological predictions of cold dark matter (CDM) and MOND at those scales where these are most successful. We also devise a more general form of the interaction term, introducing the medium density as a new order parameter. This allows for new physical effects which should be amenable to observational tests in the near future. Hence, this very general framework, which can be furthermore related to a generalized scalar-tensor theory, opens the way to a possible unification of the successes of CDM and MOND at different scales.

Intercluster filaments negligibly contribute to the weak lensing signal in general relativity (GR), $\gamma_{N}\sim 10^{-4}-10^{-3}$. In the context of relativistic modified Newtonian dynamics (MOND) introduced by Bekenstein, however, a single filament inclined by $\approx 45^\circ$ from the line of sight can cause substantial distortion of background sources pointing towards the filament's axis ($\kappa=\gamma=(1-A^{-1})/2\sim 0.01$); this is rigorous for infinitely long uniform filaments, but also qualitatively true for short filaments ($\sim 30$Mpc), and even in regions where the projected matter density of the filament is equal to zero. Since galaxies and galaxy clusters are generally embedded in filaments or are projected on such structures, this contribution complicates the interpretation of the weak lensing shear map in the context of MOND. While our analysis is of mainly theoretical interest providing order-of-magnitude estimates only, it seems safe to conclude that when modeling systems with anomalous weak lensing signals, e.g. the "bullet cluster" of Clowe et al., the "cosmic train wreck" of Abell 520 from Mahdavi et al., and the "dark clusters" of Erben et al., filamentary structures might contribute in a significant and likely complex fashion. On the other hand, our predictions of a (conceptual) difference in the weak lensing signal could, in principle, be used to falsify MOND/TeVeS and its variations.

Strong gravitational lensing by galaxies in MOdified Newtonian Dynamics (MOND) has until now been restricted to spherically symmetric models. These models were able to account for the size of the Einstein ring of observed lenses, but were unable to account for double-imaged systems with collinear images, as well as four-image lenses. Non-spherical models are generally cumbersome to compute numerically in MOND, but we present here a class of analytic non-spherical models that can be applied to fit double-imaged and quadruple-imaged systems. We use them to obtain a reasonable MOND fit to ten double-imaged systems, as well as to the quadruple-imaged system Q2237+030 which is an isolated bulge-disc lens producing an Einstein cross. However, we also find five double-imaged systems and three quadruple-imaged systems for which no reasonable MOND fit can be obtained with our models. We argue that this is mostly due to the intrinsic limitation of the analytic models, even though the presence of small amounts of additional dark mass on galaxy scales in MOND is also plausible.

Although very successful in explaining the observed conspiracy between the baryonic distribution and the gravitational field in spiral galaxies without resorting to dark matter (DM), the modified Newtonian dynamics (MOND) paradigm still requires DM in X-ray bright systems. Here, to get a handle on the distribution and importance of this DM, and thus on its possible form, we deconstruct the mass profiles of 26 X-ray emitting systems in MOND, with temperatures ranging from 0.5 to 9 keV. Initially we compute the MOND dynamical mass as a function of radius, then subtract the known gas mass along with a component of galaxies which includes the cD galaxy with $M/L_K=1$. Next we test the compatibility of the required DM with ordinary massive neutrinos at the experimental limit of detection ($m_{\nu}=2$ eV), with density given by the Tremaine-Gunn limit. Even by considering that the neutrino density stays constant and maximal within the central 100 or 150 kpc (which is the absolute upper limit of a possible neutrino contribution there), we show that these neutrinos can never account for the required DM within this region. The natural corollary of this finding is that, whereas clusters (T $\ga$ 3 keV) might have most of their mass accounted for if ordinary neutrinos have a 2 eV mass, groups (T $\lsim$ 2 keV) cannot be explained by a 2 eV neutrino contribution. This means that, for instance, cluster baryonic dark matter (CBDM, Milgrom 2007) or even sterile neutrinos would present a more satisfactory solution to the problem of missing mass in MOND X-ray emitting systems.

We compute the Milky Way potential in different cold dark matter (CDM) based models, and compare these with the modified Newtonian dynamics (MOND) framework. We calculate the axis ratio of the potential in various models, and find that isopotentials are less spherical in MOND than in CDM potentials. As an application of these models, we predict the escape velocity as a function of the position in the Galaxy. This could be useful in comparing with future data from planned or already-underway kinematic surveys (RAVE, SDSS, SEGUE, SIM, GAIA or the hypervelocity stars survey). In addition, the predicted escape velocity is compared with the recently measured high proper motion velocity of the Large Magellanic Cloud (LMC). To bind the LMC to the Galaxy in a MOND model, while still being compatible with the RAVE-measured local escape speed at the Sun's position, we show that an external field modulus of less than $0.03 a_0$ is needed.

We present the analysis of 23 published rotation curves of disk galaxies belonging to the Ursa Major group of galaxies, with kinematics free of irregularities. The rotation curves are analysed in the context of MOND (Modified Newtonian Dynamics). We add an extra component to the rotation curve fits, in addition to the stellar and gaseous disks: a speculative halo of constant density made of, e.g., neutrinos, which would solve the bulk of the problem currently faced by MOND on rich galaxy clusters scales. We find that this additional unseen mass density is poorly constrained (as expected a priori, given that a neutrino halo never dominates the kinematics), but we also find that the best-fit value is non-zero: rho = 3.8 x 10^{-27} g/cm^3, and that a zero-density is marginally excluded with 87% confidence; also, the 95% confidence upper limit for the density is rho = 9.6 x 10^{-27} g/cm^3. These limits are slightly above the expectations from the Tremaine-Gunn phase space constraints on ordinary 2 eV neutrinos, but in accordance with the maximum density expected for one or two species of 5 eV sterile neutrinos.

We investigate the nature of the classical low-velocity structures in the local velocity field, i.e. the Pleiades, Hyades and Sirius moving groups. After using a wavelet transform to locate them in velocity space, we study their relation with the open clusters kinematically associated with them. By directly comparing the location of moving group stars in parallax space to the isochrones of the embedded clusters, we check whether, within the observational errors on the parallax, all moving group stars could originate from the on-going evaporation of the associated cluster. We conclude that, in each moving group, the fraction of stars making up the velocity-space overdensity superimposed on the background is higher than the fraction of stars compatible with the isochrone of the associated cluster. These observations thus favour a dynamical (resonant) origin for the Pleiades, Hyades and Sirius moving groups.

We summarize the status of Modified Newtonian Dynamics (MOND) in galaxy clusters. The observed acceleration is typically larger than the acceleration threshold of MOND in the central regions, implying that some dark matter is necessary to explain the mass discrepancy there. A plausible resolution of this issue is that the unseen mass in MOND is in the form of ordinary neutrinos with masses just below the experimentally detectable limit. In particular, we show that the lensing mass reconstructions of the clusters 1E0657-56 (the bullet cluster) and Cl0024+17 (the ring) do not pose a new challenge to this scenario. However, the mass discrepancy for cool X-ray emitting groups, in which neutrinos cannot cluster, pose a more serious problem, meaning that dark baryons could present a more satisfactory solution to the problem of unseen mass in MOND clusters.

We investigate how different models that have been proposed for solving the dark matter problem can fit the velocity dispersion observed around elliptical galaxies, on either a small scale (~ 20kpc) with stellar tracers, such as planetary nebulae, or large scale (~ 200kpc) with satellite galaxies as tracers. Predictions of Newtonian gravity, either containing pure baryonic matter, or embedded in massive cold dark matter (CDM) haloes, are compared with predictions of the modified gravity of MOND. The standard CDM model has problems on a small scale, and the Newtonian pure baryonic model has difficulties on a large scale, while a fit with MOND is possible on both scales.

The recent SDSS measured velocity distribution of satellite galaxies has been modelled in the context of MOND. We show that even when the extra constraint of adhering to the projected satellite number density profile is added, the two line of sight (los) velocity dispersion profiles presented in Klypin & Prada (2007) can be matched simply with a radially varying anisotropy. Interestingly, the anisotropies required to fit the los velocity dispersions are remarkably similar to the anisotropies generated by dissipationless collapse simulations in MOND. The mass-to-light ratios of the two host galaxies used are sensible and positivity of the distribution function is satisfied.

Within the cold dark matter (CDM) framework tidal dwarf galaxies (TDGs) cannot contain dark matter, so the recent results by Bournaud et al. (2007) that 3 rotating TDGs do show significant evidence for being dark matter dominated is inconsistent with the current concordance cosmological theory unless yet another dark matter component is postulated. We confirm that the TDG rotation curves are consistent with Newtonian dynamics only if either an additional dark matter component is postulated, or if all 3 TDGs happen to be viewed nearly edge-on, which is unlikely given the geometry of the tidal debris. We also find that the observed rotation curves are very naturally explained without any free parameters within the modified Newtonian dynamics (MOND) framework if inclinations are adopted as derived by Bournaud et al. We explore different inclination angles and two different assumptions about the external field effect. The results do not change significantly, and we conclude therefore that Newtonian dynamics has severe problems while MOND does exceedingly well in explaining the observed rotation curves of the 3 TDGs studied by Bournaud et al.

The self-binding energy and stability of a galaxy in MOND-based gravity are curiously decreasing functions of its center of mass acceleration towards neighbouring mass concentrations. A tentative indication of this breaking of the Strong Equivalence Principle in field galaxies is the RAVE-observed escape speed in the Milky Way. Another consequence is that satellites of field galaxies will move on nearly Keplerian orbits at large radii (100 - 500 kpc), with a declining speed below the asymptotically constant naive MOND prediction. But consequences of an environment-sensitive gravity are even more severe in clusters, where member galaxies accelerate fast: no more Dark-Halo-like potential is present to support galaxies, meaning that extended axisymmetric disks of gas and stars are likely unstable. These predicted reappearance of asymptotic Keplerian velocity curves and disappearance of "stereotypic galaxies" in clusters are falsifiable with targeted surveys.

We present a new test of modified Newtonian dynamics (MOND) on galactic scales, based on the escape speed in the solar neighbourhood. This test is independent from other empirical successes of MOND at reproducing the phenomenology of galactic rotation curves. The galactic escape speed in MOND is entirely determined by the baryonic content of the Galaxy and the external field in which it is embedded. We estimate that the external field in which the Milky Way must be embedded to produce the observed local escape speed of 550 km/s is of the order of a_0/100, where a_0 is the dividing acceleration scale below which gravity is boosted in MOND. This is compatible with the external gravitational field actually acting on the Milky Way.

Observations of spiral galaxies strongly support a one-to-one analytical relation between the inferred gravity of dark matter at any radius and the enclosed baryonic mass. It is baffling that baryons manage to settle the dark matter gravitational potential in such a precise way, leaving no "messy" fingerprints of the merging events and "gastrophysical" feedbacks expected in the history of a galaxy in a concordance Universe. This correlation of gravity with baryonic mass can be interpreted from several non-standard angles, especially as a modification of gravity called TeVeS, in which no galactic dark matter is needed. In this theory, the baryon-gravity relation is captured by the dieletric-like function mu of Modified Newtonian Dynamics (MOND), controlling the transition from 1/r^2 attraction in the strong gravity regime to 1/r attraction in the weak regime. Here, we study this mu-function in detail. We investigate the observational constraints upon it from fitting galaxy rotation curves, unveiling the degeneracy between the stellar mass-to-light ratio and the mu-function as well as the importance of the sharpness of transition from the strong to weak gravity regimes. We also numerically address the effects of non-spherical baryon geometry in the framework of non-linear TeVeS, and exhaustively examine how the mu-function connects with the free function of that theory. In that regard, we exhibit the subtle effects and wide implications of renormalizing the gravitational constant. We finally present a discontinuity-free transition between quasi-static galaxies and the evolving Universe for the free function of TeVeS, inevitably leading to a return to 1/r^2 attraction at very low accelerations in isolated galaxies.

We develop a new method to predict the density associated with weak lensing maps of (un)relaxed clusters in a range of theories interpolating between GR and MOND (General Relativity and Modified Newtonian Dynamics). We apply it to fit the lensing map of the bullet merging cluster 1E0657-56, in order to constrain more robustly the nature and amount of collisionless matter in clusters {\it beyond} the usual assumption of spherical equilibrium (Pointecouteau & Silk 2005) and the validity of GR on cluster scales (Clowe et al. 2006). Strengthening the proposal of previous authors we show that the bullet cluster is dominated by a collisionless -- most probably non-baryonic -- component in GR as well as in MOND, a result consistent with the dynamics of many X-ray clusters. Our findings add to the number of known pathologies for a purely baryonic MOND, including its inability to fit the latest data from the Wilkinson Microwave Anisotropy Probe. A plausible resolution of all these issues and standard issues of Cold Dark Matter with galaxy rotation curves is the "marriage" of MOND with ordinary hot neutrinos of 2eV. This prediction is just within the GR-independent maximum of neutrino mass from current $\beta$-decay experiments, and is falsifiable by the Karlsruhe Tritium Neutrino (KATRIN) experiment by 2009. Issues of consistency with strong lensing arcs and the large relative velocity of the two clusters comprising the bullet cluster are also addressed.

The nature of the Hyades stream, or Hyades moving group, is a long-standing question of Galactic Astronomy. Is this stream caused by non-axisymmetric perturbations of the galactic potential, such as transient or quasi-stationary spiral waves, or by the on-going evaporation of the Hyades cluster? A simple observational test has been designed to answer that question. Using the Geneva-Copenhagen survey of F and G dwarfs, we compare the mass distribution and metallicity of the stream to those of field disk stars. If the Hyades stream is composed of stars trapped at resonance, its mass distribution should obey the present-day mass function (PDMF) of the disk, and its metallicity should reflect its origin in the inner regions of the Galaxy. On the other hand, if it is an evaporated cluster, we expect a different mass distribution, depending on the inital mass function (IMF) of the cluster, and on the proportion of evaporated stars as a function of mass. We find that extreme conditions have to be adopted for the selective evaporation and IMF of the cluster to make the observed mass ditribution of the stream only roughly consistent (at a one-sigma level) with the coeval evaporated cluster scenario. The observed mass distribution is in much better agreement with the PDMF of the field. We also note that the peculiar metallicity of the stream is inconsistent with that of a field population from the solar neighbourhood trapped in the primordial cluster during its formation process and subsequently evaporated. These observations thus favour a resonant origin for the Hyades stream.

A proper test of Modified Newtonian Dynamics (MOND) in systems of non-trivial geometries depends on modelling subtle differences in several versions of its postulated theories. This is especially important for lensing and dynamics of barely virialised galaxy clusters with typical gravity of scale $\sim a_0 \sim 1\AA{\rm s}^{-2}$. The original MOND formula, the classical single field modification of the Poisson equation, and the multi-field general relativistic theory of Bekenstein (TeVeS) all lead to different predictions as we stray from spherical symmetry. In this paper, we study a class of analytical MONDian models for a system with a semi-Hernquist baryonic profile. After presenting the analytical distribution function of the baryons in spherical limits, we develop orbits and gravitational lensing of the models in non-spherical geometries. In particular, we can generate a multi-centred baryonic system with a weak lensing signal resembling that of the merging galaxy cluster 1E 0657-56 with a bullet-like light distribution. We finally present analytical scale-free highly non-spherical models to show the subtle differences between the single field classical MOND theory and the multi-field TeVeS theory.

With [Fe/H] = -3.5, CS 30322-023 is the most metal-poor star to exhibit a clear s-process signature and the most metal-poor ``lead star'' known. CS 30322-023 is also remarkable in having the lowest surface gravity (log g <= -0.3) among the metal-poor stars studied to date. The available evidence indicates that this star is presently a thermally-pulsing asymptotic giant branch (TP-AGB) star, with no strong indication of binarity thus far (although a signal of period 192 d is clearly present in the radial-velocity data, this is likely due to pulsation of the stellar envelope). We show that low-mass TP-AGB stars are not expected to be exceedingly rare in a magnitude-limited sample such as the HK survey, because their high luminosities make it possible to sample them over a very large volume. The strong N overabundance and the low 12C/13C ratio (4) in this star is typical of the operation of the CN cycle. Coupled with a Na overabundance and the absence of a strong C overabundance, this pattern seems to imply that hot-bottom burning operated in this star, which should then have a mass of at least 2 Msun. However, the luminosity associated with this mass would put the star at a distance of about 50 kpc, in the outskirts of the galactic halo. We explore alternative scenarios in which the observed abundance pattern results from some mixing mechanism yet to be identified occurring in a single low-metallicity 0.8 Msun AGB star, or from pollution by matter from an intermediate-mass AGB companion which has undergone hot-bottom burning. We stress, however, that our abundances may be subject to uncertainties due to NLTE or 3D granulation effects which were not taken into consideration.

The phenomena customly called Dark Matter or Modified Newtonian Dynamics (MOND) have been argued by Bekenstein (2004) to be the consequences of a covariant scalar field, controlled by a free function (related to the MOND interpolating function) in its Lagrangian density. In the context of this relativistic MOND theory (TeVeS), we examine critically the interpolating function in the transition zone between weak and strong gravity. Bekenstein's toy model produces too gradually varying functions and fits rotation curves less well than the standard MOND interpolating function. However, the latter varies too sharply and implies an implausible external field effect (EFE). These constraints on opposite sides have not yet excluded TeVeS, but made the zone of acceptable interpolating functions narrower. An acceptable "toy" Lagrangian density function with simple analytical properties is singled out for future studies of TeVeS in galaxies. We also suggest how to extend the model to solar system dynamics and cosmology, and compare with strong lensing data (see also astro-ph/0509590).

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 true nature of the Hyades and Sirius superclusters is still an open question. In this contribution, we confront Eggen's hypothesis that they are cluster remnants with the results of a kinematic analysis of more than 6000 K and M giants in the solar neighbourhood. This analysis includes new radial velocity data from a large survey performed with the Coravel spectrometer, complemented by Hipparcos parallaxes and Tycho-2 proper motions (Famaey et al. 2004). A maximum-likelihood method, based on a bayesian approach, has been applied to the data, in order to make full use of all the available data (including less precise parallaxes) and to derive the properties of the different kinematic subgroups. Two such subgroups can be identified with the Hyades and Sirius superclusters. Stars belonging to them span a very wide range of age, which is difficult to account for in Eggen's scenario. These groups are thus most probably "dynamical streams" related to the dynamical perturbation by spiral waves rather than to cluster remnants. In this scenario, the Hyades and Ursa Major clusters just happen to be in the Hyades and Sirius streams, which are purely dynamical features that have nothing to do with the remnants of more massive primordial clusters. This mechanism could be the key to understanding the presence of an old metal-rich population, and of many exoplanetary systems in our neighbourhood. Moreover, a strong spiral pattern seems to be needed in order to yield such prominent streams. Since spiral structure is usually baryonic, this would leave very little room for dark matter. This may be an indication that the era of the dark-matter paradigm explaining the dynamics of the Galaxy may come to an end, and is being superseded by modified gravity.

The availability of the Hipparcos Catalogue triggered many kinematic and dynamical studies of the solar neighbourhood. Nevertheless, those studies generally lacked the third component of the space velocities, i.e., the radial velocities. This work presents the kinematic analysis of 5952 K and 739 M giants in the solar neighbourhood which includes for the first time radial velocity data from a large survey performed with the CORAVEL spectrovelocimeter. It also uses proper motions from the Tycho-2 catalogue, which are expected to be more accurate than the Hipparcos ones. The UV-plane constructed from these data for the stars with precise parallaxes reveals a rich small-scale structure, with several clumps corresponding to the Hercules stream, the Sirius moving group, and the Hyades and Pleiades superclusters. A maximum-likelihood method, based on a bayesian approach, has been applied to the data, in order to make full use of all the available stars and to derive the kinematic properties of these subgroups. Isochrones in the Hertzsprung-Russell diagram reveal a very wide range of ages for stars belonging to these groups. These groups are most probably related to the dynamical perturbation by transient spiral waves rather than to cluster remnants. A possible explanation for the presence of young clusters in the same area of the UV-plane is that they have been put there by the spiral wave associated with their formation, while the kinematics of the older stars of our sample has also been disturbed by the same wave. The term "dynamical stream" for the kinematic groups is thus more appropriate than the traditional term "supercluster" since it involves stars of different ages, not born at the same place nor at the same time.

We present a set of three-component Stackel potentials defined by five parameters and designed to model the Milky Way. We review the fundamental constraints that any model of the Milky Way must satisfy, including the most recent ones derived from Hipparcos data, and we study how the parameters of the presented potentials can vary in order to match these constraints. Five different valid potentials are presented and analyzed in detail: they are designed to be confronted with kinematical surveys in the future, by the construction of three-integral analytic distribution functions.

We present new equilibrium component distribution functions that depend on three analytic integrals in a Stackel potential, and that can be used to model stellar discs of galaxies. These components are generalizations of two-integral ones and can thus provide thin discs in the two-integral approximation. Their most important properties are the partly analytical expression for their moments, the disc-like features of their configuration space densities (exponential decline in the galactic plane and finite extent in the vertical direction) and the anisotropy of their velocity dispersions. We further show that a linear combination of such components can fit a van der Kruit disc.

We present new tools to establish axisymmetric equilibrium models of the Milky Way. The models we wish to establish are pairs (V,F) where V is the gravitational potential generated by the whole mass distribution including the dark matter, and F is the distribution function in phase space for late-type tracer stars. We present a set of Stackel potentials that fit some fundamental parameters of the Milky Way (mass density in the solar neighbourhood and Oort constants). Then we define new component distribution functions that can be combined with these potentials in order to reproduce kinematical data like those that will be provided by GAIA.