Source author record

Benoit Famaey

Benoit Famaey appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

62works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

62 published item(s)

preprint2026arXiv

Matière noire et (ou) gravitation modifiée : une approche historique et épistémologique

The current standard model of cosmology assumes that the majority of matter in the Universe is made of dark matter, and that the latter is fundamentally different from ordinary matter. Dark matter can in principle explain the rotation of galaxies, the gravitational lensing from galaxy clusters or the appearance of the cosmic microwave background, the oldest light in the Universe. But does dark matter really exist? Here, we review the history of this concept and its implications for the formation and evolution of galaxies. We also consider the questions that remain, the limitations of the model, and present alternative theories, in particular modifications to the gravitional law that would -- perhaps -- make it possible to do without it.

preprint2022arXiv

3D hydrodynamic simulations for the formation of the Local Group satellite planes

The existence of mutually correlated thin and rotating planes of satellite galaxies around both the Milky Way (MW) and Andromeda (M31) calls for an explanation. Previous work in Milgromian dynamics (MOND) indicated that a past MW-M31 encounter might have led to the formation of these satellite planes. We perform the first-ever hydrodynamical MOND simulation of the Local Group using Phantom of RAMSES. We show that an MW-M31 encounter at $z \approx 1$, with a perigalactic distance of about 80 kpc, can yield two disc galaxies at $z=0$ oriented similarly to the observed galactic discs and separated similarly to the observed M31 distance. Importantly, the tidal debris are distributed in phase space similarly to the observed MW and M31 satellite planes, with the correct preferred orbital pole for both. The MW-M31 orbital geometry is consistent with the presently observed M31 proper motion despite this not being considered as a constraint when exploring the parameter space. The mass of the tidal debris around the MW and M31 at $z=0$ compare well with the mass observed in their satellite systems. The remnant discs of the two galaxies have realistic radial scale lengths and velocity dispersions, and the simulation naturally produces a much hotter stellar disc in M31 than in the MW. However, reconciling this scenario with the ages of stellar populations in satellite galaxies would require that a higher fraction of stars previously formed in the outskirts of the progenitors ended up within the tidal debris, or that the MW-M31 interaction occurred at $z>1$.

preprint2022arXiv

Constraining the Milky Way halo kinematics via its Linear Response to the Large Magellanic Cloud

We model the response of spherical, non-rotating Milky Way (MW) dark matter and stellar halos to the Large Magellanic Cloud (LMC) using the matrix method of linear response theory. Our computations reproduce the main features of the dark halo response from simulations. We show that these features can be well separated by a harmonic decomposition: the large scale over/underdensity in the halo (associated with its reflex motion) corresponds to the $\ell=1$ terms, and the local overdensity to the $\ell\geq2$ multipoles. Moreover, the dark halo response is largely dominated by the first order 'forcing' term, with little influence from self-gravity. This makes it difficult to constrain the underlying velocity distribution of the dark halo using the observed response of the stellar halo, but it allows us to investigate the response of stellar halo models with various velocity anisotropies: a tangential (respectively radial) halo produces a shallower (respectively stronger) response. We also show that only the local wake is responsible for these variations, the reflex motion being solely dependent on the MW potential. Therefore, we identify the structure (orientation and winding) of the in-plane quadrupolar ($m=2$) response as a potentially good probe of the stellar halo anisotropy. Finally, our method allows us to tentatively relate the wake strength and shape to resonant effects: the strong radial response could be associated with the inner Lindblad resonance, and the weak tangential one with corotation.

preprint2022arXiv

On the Effect of the Large Magellanic Cloud on the Orbital Poles of Milky Way Satellite Galaxies

The reflex motion and distortion of the Milky Way (MW) halo caused by the infall of a massive Large Magellanic Cloud (LMC) has been demonstrated to result in an excess of orbital poles of dark matter halo particles towards the LMC orbital pole. This was suggested to help explain the observed preference of MW satellite galaxies to co-orbit along the Vast Polar Structure (VPOS). We test this idea by correcting the positions and velocities of the MW satellites for the Galactocentric-distance-dependent shifts inferred from a LMC-infall simulation. While this should substantially reduce the observed clustering of orbital poles if it were mainly caused by the LMC, we instead find that the strong clustering remains preserved. We confirm the initial study's main result with our simulation of an MW-LMC-like interaction, and use it to identify two reasons why this scenario is unable to explain the VPOS: (1) the orbital pole density enhancement in our simulation is very mild (~10% within 50-250 kpc) compared to the observed enhancement (~220-300%), and (2) it is very sensitive to the specific angular momenta (AM) of the simulation particles, with higher AM particles being affected the least. Particles in simulated dark matter halos tend to follow more radial orbits (lower AM), so their orbital poles are more easily affected by small offsets in position and velocity caused by an LMC infall than objects with more tangential velocity (higher AM), such as the observed dwarf galaxies surrounding the MW. The origin of the VPOS thus remains unexplained.

preprint2022arXiv

Origin of the spectacular tidal shells of galaxy NGC474

The lenticular galaxy NGC474 hosts a rich system of tidal shells and streams, some of which are exceptionally bright. Two teams recently presented spectroscopic observations of the brightest shells. These were the first shell spectra ever observed in integrated starlight. The authors studied the stellar populations of the shell, of the center of the galaxy and of its globular clusters. The precise formation scenario for the tidal features of this prominent galaxy however still remained unclear. Here, we add further clues on their formation from the radii of the shells, and we present a scenario for the formation of the tidal features that seems to be unique and explaining all available data. Shell radii are analyzed with the shell identification method, and we run self-consistent simulations of the formation of the tidal features. We consider Newtonian as well as MOND gravity. Observations suggest that the tidal features originate from the accretion of a spiral galaxy. The shell identification method yields that the merging galaxies collided first 1.3Gyr ago and then again 0.9Gyr ago, thereby forming the shells in two generations. This would also explain the young ages of stellar populations in the center of the galaxy and the young age of the globular clusters. The analytic models of shell propagation, that underlie the shell identification method, are verified by a simulation. The simulations reproduce well the observed morphology of the tidal features. The accreted spiral likely reached NGC474 nearly radially, in the plane of the sky, from the south, its rotation axis pointing toward us. It should have had a stellar mass of around 1/6 of NGC474, i.e. $10^{9.8}\,M_\odot$. It seems that all tidal features in the galaxy originate from one merger.

preprint2022arXiv

Revisiting a disky origin for the faint branch of the Sagittarius stellar stream

We investigate ways to produce the bifurcation observed in the stellar stream of the Sagittarius dwarf galaxy (Sgr). Our method consists in running $N$-body simulations of Sgr falling into the Milky Way for the last 3~Gyr, with added test particles on disk orbits that span a wide range of initial positions, energies, and angular momenta. We find that particles that end up in the faint branch are predominantly high angular momentum particles that can all originate from a single plane within the progenitor, nearly perpendicular both to the orbital plane of the progenitor and to the Milky Way stellar disk. Their original configuration at the start of the simulation corresponds to spiral features already present 3~Gyr ago, which could be, e.g., the result of a disk-like component being tidally perturbed, or the tidal tails of a satellite being disrupted within Sgr. We then run a simulation including the self-gravity of this disky component. Despite the remaining ambiguity of its origin, this disk component of the Sgr dwarf with spiral over-densities provides a first step towards a working model to reproduce the observed faint branch of the bifurcated Sgr stream.

preprint2022arXiv

The Complexity of the Cetus Stream Unveiled from the Fusion of STREAMFINDER and StarGO

We combine the power of two stream-searching tools, STREAMFINDER and StarGO applied to the Gaia EDR3 data, to detect stellar debris belonging to the Cetus stream system that forms a complex, nearly polar structure around the Milky Way. In this work, we find the southern extensions of the northern Cetus stream as the Palca stream and a new southern stream, which overlap on the sky but have different distances. These two stream wraps extend over more than $\sim100°$ on the sky ($-60°<δ<+40°$). The current N-body model of the system reproduces both wraps in the trailing arm. We also show that the Cetus system is confidently associated with the Triangulum/Pisces, Willka Yaku, and the recently discovered C-20 streams. The association with the ATLAS-Aliqa Uma stream is much weaker. All of these stellar debris are very metal-poor, comparable to the average metallicity of the southern Cetus stream with [Fe/H] = $-2.17\pm0.2$. The estimated stellar mass of the Cetus progenitor is at least $10^{5.6}\rm\,M_\odot$, compatible with Ursa Minor or Draco dwarf galaxies. The associated globular cluster with similar stellar mass, NGC 5824 very possibly was accreted in the same group infall. The multi-wrap Cetus stream is a perfect example of a dwarf galaxy that has undergone several periods of stripping, leaving behind debris at multiple locations in the halo. The full characterization of such systems is crucial to unravel the assembly history of the Milky Way and, as importantly, to provide nearby fossils to study ancient low-mass dwarf galaxies.

preprint2022arXiv

The formation of early-type galaxies through monolithic collapse of gas clouds in Milgromian gravity

Studies of stellar populations in early-type galaxies (ETGs) show that the more massive galaxies form earlier and have a shorter star formation history (SFH). In this study, we investigate the initial conditions of ETG formation. The study begins with the collapse of non-rotating post-Big-Bang gas clouds in Milgromian (MOND) gravitation. These produce ETGs with star-forming timescales (SFT) comparable to those observed in the real Universe. Comparing these collapse models with observations, we set constraints on the initial size and density of the post-Big-Bang gas clouds in order to form ETGs. The effective-radius-mass relation of the model galaxies falls short of the observed relation. Possible mechanisms for later radius expansion are discussed. Using hydrodynamic MOND simulations this work thus for the first time shows that the SFTs observed for ETGs may be a natural occurrence in the MOND paradigm. We show that different feedback algorithms change the evolution of the galaxies only to a very minor degree in MOND. The first stars have, however, formed more rapidly in the real Universe than possible just from the here studied gravitational collapse mechanism. Dark-matter-based cosmological structure formation simulations disagree with the observed SFTs at more than 5 sigma confidence.

preprint2022arXiv

The Global Dynamical Atlas of the Milky Way mergers: Constraints from Gaia EDR3 based orbits of globular clusters, stellar streams and satellite galaxies

The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams and satellite galaxies in action ($\mathbf{J}$) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy ($\textit{E}$) to enhance the contrast between groupings. For this endeavour, we use $\textit{Gaia}$ EDR3 based measurements of $170$ globular clusters, $41$ streams and $46$ satellites to derive their $\mathbf{J}$ and $\textit{E}$. To detect groups, we use the $\texttt{ENLINK}$ software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of $N=6$ groups, including the previously known mergers $\textit{Sagittarius}$, $\textit{Cetus}$, $\textit{Gaia-Sausage/Enceladus}$, $\textit{LMS-1/Wukong}$, $\textit{Arjuna/Sequoia/I'itoi}$ and one new merger that we call $\textit{Pontus}$. All of these mergers, together, comprise $62$ objects ($\approx 25\%$ of our sample). We discuss their members, orbital properties and metallicity distributions. We find that the three most metal-poor streams of our Galaxy -- "C-19" ([Fe/H]$=-3.4$ dex), "Sylgr" ([Fe/H]$=-2.9$ dex) and "Phoenix" ([Fe/H]$=-2.7$ dex) -- are associated with $\textit{LMS-1/Wukong}$; showing it to be the most metal-poor merger. The global dynamical atlas of Milky Way mergers that we present here provides a present-day reference for galaxy formation models.

preprint2022arXiv

The global stability of M33 in MOND

The dynamical stability of disk galaxies is sensitive to whether their anomalous rotation curves are caused by dark matter halos or Milgromian Dynamics (MOND). We investigate this by setting up a MOND model of M33. We first simulate it in isolation for 6 Gyr, starting from an initial good match to the rotation curve (RC). Too large a bar and bulge form when the gas is too hot, but this is avoided by reducing the gas temperature. A strong bar still forms in 1 Gyr, but rapidly weakens and becomes consistent with the observed weak bar. Previous work showed this to be challenging in Newtonian models with a live dark matter halo, which developed strong bars. The bar pattern speed implies a realistic corotation radius of 3 kpc. However, the RC still rises too steeply, and the central line of sight velocity dispersion (LOSVD) is too high. We then add a constant external acceleration field of $8.4 \times 10^{-12}$ m/s$^2$ at $30^\circ$ to the disk as a first order estimate for the gravity exerted by M31. This suppresses buildup of material at the centre, causing the RC to rise more slowly and reducing the central LOSVD. Overall, this simulation bears good resemblance to several global properties of M33, and highlights the importance of including even a weak external field on the stability and evolution of disk galaxies. Further simulations with a time-varying external field, modeling the full orbit of M33, will be needed to confirm its resemblance to observations.

preprint2022arXiv

The Pristine survey -- XVII. The C-19 stream is dynamically hot and more extended than previously thought

The C-19 stream is the most metal poor stellar system ever discovered, with a mean metallicity $[Fe/H] = -3.38\pm0.06$. Its low metallicity dispersion ($σ_{\rm [Fe/H]}$ $<$ 0.18 at the 95\% confidence level) as well as variations in sodium abundances strongly suggest a globular cluster origin. In this work, we use VLT/UVES spectra of seven C-19 stars to derive more precise velocity measurements for member stars, and to identify two new members with radial velocities and metallicities consistent with the stream's properties. One of these new member stars is located 30 $°$ away from the previously identified body of C-19, implying that the stream is significantly more extended than previously known and that more members likely await discovery. In the main part of C-19, we measure a radial velocity dispersion $σ_v$ = 6.2$^{+2.0}_{-1.4}$ km s$^{-1}$ from nine members, and a stream width of 0.56$°\pm0.08°$, equivalent to $\sim$158 pc at a heliocentric distance of 18 kpc. These confirm that C-19 is comparatively hotter, dynamically, than other known globular cluster streams and shares the properties of faint dwarf galaxy streams. On the other hand, the variations in the Na abundances of the three newly observed bright member stars, the variations in Mg and Al for two of them, and the normal Ba abundance of the one star where it can be measured provide further evidence for a globular cluster origin. The tension between the dynamical and chemical properties of C-19 suggests that its progenitor experienced a complex birth environment or disruption history.

preprint2021arXiv

The Cen A galaxy group: dynamical mass and missing baryons

The nearby elliptical galaxy Cen A is surrounded by a flattened system of dwarf satellite galaxies with coherent motions. Using a novel Bayesian approach, we measure the mean rotation velocity $v_{\rm rot}$ and velocity dispersion $σ_{\rm int}$ of the satellite system. We find $v_{\rm rot}/σ_{\rm int} \simeq 0.7$ indicating that the satellite system has non-negligible rotational support. Using Jeans' equations, we measure a circular velocity of 258 km s$^{-1}$ and a dynamical mass of $1.2\times 10^{13}$ M$_\odot$ within 800 kpc. In a $Λ$CDM cosmological context, we find that the Cen A group has a baryon fraction $M_{\rm b}/M_{200}\simeq0.035$ and is missing $\sim$77$\%$ of the cosmologically available baryons. Consequently, Cen A should have a hot intergalactic medium with a mass of $\sim$8$\times$10$^{11}$ M$_\odot$, which is more than $\sim$20 times larger than current X-ray estimates. Intriguingly, The whole Cen A group lies on the baryonic Tully-Fisher relation defined by individual rotationally supported galaxies, as expected in Milgromian dynamics (MOND) with no need of missing baryons.

preprint2020arXiv

A panoramic landscape of the Sagittarius stream in Gaia DR2 revealed with the STREAMFINDER spyglass

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.

preprint2020arXiv

Baryon-Interacting Dark Matter: heating dark matter and the emergence of galaxy scaling relations

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.

preprint2020arXiv

Detection of Strong Epicyclic Density Spikes in the GD-1 Stellar Stream. An Absence of Evidence for the Influence of Dark Matter Subhalos?

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.

preprint2020arXiv

The formation of exponential disk galaxies in MOND

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.

preprint2020arXiv

The impact of the halo spin-concentration relation on disc scaling laws

Galaxy scaling laws, such as the Tully-Fisher, mass-size and Fall relations, can provide extremely useful clues on our understanding of galaxy formation in a cosmological context. Some of these relations are extremely tight and well described by one single parameter (mass), despite the theoretical existence of secondary parameters such as spin and concentration, which are believed to impact these relations. In fact, the residuals of these scaling laws appear to be almost uncorrelated with each other, posing significant constraints on models where secondary parameters play an important role. Here, we show that a possible solution is that such secondary parameters are correlated amongst themselves, in a way that removes correlations in observable space. In particular, we focus on how the existence of an anti-correlation between the dark matter halo spin and its concentration -- which is still debated in simulations -- can weaken the correlation of the residuals of the Tully-Fisher and mass-size relations. Interestingly, using simple analytic galaxy formation models, we find that this happens only for a relatively small portion of the parameter space that we explored, which suggests that this idea could be used to derive constraints to galaxy formation models that are still unexplored.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- I: Survey Description, Spectra and Radial Velocities

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.

preprint2020arXiv

The Sixth Data Release of the Radial Velocity Experiment (RAVE) -- II: Stellar Atmospheric Parameters, Chemical Abundances and Distances

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.

preprint2019arXiv

Phase-space correlation in stellar streams of the Milky Way halo: The clash of Kshir and GD-1

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.

preprint2016arXiv

Exploring the reality of density substructures in the Palomar 5 stellar stream

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.

preprint2016arXiv

Masses of the components of SB2 binaries observed with Gaia. III. Accurate SB2 orbits for 10 binaries and masses of HIP 87895

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.

preprint2016arXiv

Modelling the Galactic disc: perturbed distribution functions in the presence of spiral arms

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.

preprint2016arXiv

Simulating disk galaxies and interactions in Milgromian dynamics

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üghausen (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.

preprint2016arXiv

Star formation triggered by galaxy interactions in modified gravity

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.

preprint2016arXiv

Staying away from the bar: the local dynamical signature of slow and fast bars in the Milky Way

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.

preprint2016arXiv

The effects of bar-spiral coupling on stellar kinematics in the Galaxy

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.

preprint2015arXiv

Dark Matter in the Milky Way

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.

preprint2015arXiv

Dynamical measurement of the stellar surface density of face-on galaxies

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.

preprint2015arXiv

Eppur si muove: Positional and kinematic correlations of satellite pairs in the low Z universe

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 σ$ 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°$ 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.

preprint2015arXiv

Mass models of disk galaxies from the DiskMass Survey in MOND

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.

preprint2015arXiv

On the persistence of two small-scale problems in ΛCDM

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.

preprint2015arXiv

Spiral and bar driven peculiar velocities in Milky Way sized galaxy simulations

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.

preprint2015arXiv

The intrinsic ellipticity of dwarf spheroidal galaxies: constraints from the Andromeda system

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.

preprint2015arXiv

The vertical effects of disc non-axisymmetries from perturbation theory: the case of the Galactic bar

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.

preprint2014arXiv

A census of the expected properties of classical Milky Way dwarfs 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).

preprint2014arXiv

Co-orbiting satellite galaxy structures are still in conflict with the distribution of primordial dwarf galaxies

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 $Λ$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 $Λ$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 $Λ$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.

preprint2014arXiv

Galactic kinematics and dynamics from RAVE stars

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.

preprint2014arXiv

Modified Baryonic Dynamics: two-component cosmological simulations with light sterile neutrinos

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.

preprint2014arXiv

N-body simulations of the Carina dSph in MOND

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.

preprint2014arXiv

Phantom of RAMSES (POR): A new Milgromian dynamics N-body code

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.

preprint2014arXiv

Radial and vertical flows induced by galactic spiral arms: likely contributors to our "wobbly Galaxy''

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.

preprint2014arXiv

Velocity anti-correlation of diametrically opposed galaxy satellites in the low redshift universe

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σ$ 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.

preprint2013arXiv

Challenges for Lambda-CDM and MOND

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.

preprint2013arXiv

Cosmological simulations in MOND: the cluster scale halo mass function with light sterile neutrinos

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 Ω_{ν_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.

preprint2013arXiv

Local Group timing in Milgromian dynamics. A past Milky Way-Andromeda encounter at z>0.8

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.

preprint2013arXiv

The MOND phenomenology

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.

preprint2012arXiv

A QUMOND galactic N-body code I: Poisson solver and rotation curve fitting

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.

preprint2012arXiv

Gaia and the dynamics of the Galaxy

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?

preprint2012arXiv

Modified Newtonian Dynamics (MOND): Observational Phenomenology and Relativistic Extensions

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.

preprint2012arXiv

Unifying all mass discrepancies with one effective gravity law?

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.

preprint2011arXiv

Modeling disc non-axisymmetries: multiple patterns, radial migration, and thick discs

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.

preprint2011arXiv

Remarks on the properties of elliptical galaxies in modified Newtonian dynamics

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.

preprint2010arXiv

Comparing different realizations of modified Newtonian dynamics: virial theorem and elliptical shells

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.

preprint2010arXiv

Lopsidedness of cluster galaxies in modified gravity

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.

preprint2009arXiv

Current data on the globular cluster Palomar 14 are not inconsistent with MOND

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.

preprint2009arXiv

Equilibrium configurations of 11eV sterile neutrinos in MONDian galaxy clusters

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

preprint2009arXiv

The Galactic potential and the asymmetric distribution of hypervelocity stars

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.

preprint2009arXiv

Universality of galactic surface densities within one dark halo scale-length

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.