Showing 1–50 of 367 results for author: Ibata, R

Simulations of cluster ultra-diffuse galaxies in MOND

Authors: Srikanth T. Nagesh, Jonathan Freundlich, Benoit Famaey, Michal Bílek, Graeme Candlish, Rodrigo Ibata, Oliver Müller

Abstract: Ultra-diffuse galaxies (UDGs) in the Coma cluster have velocity dispersion profiles that are in full agreement with the predictions of Modified Newtonian Dynamics (MOND) in isolation. However, the external field effect (EFE) from the cluster seriously deteriorates this agreement. It has been suggested that this could be related to the fact that UDGs are out-of-equilibrium objects whose stars have… ▽ More Ultra-diffuse galaxies (UDGs) in the Coma cluster have velocity dispersion profiles that are in full agreement with the predictions of Modified Newtonian Dynamics (MOND) in isolation. However, the external field effect (EFE) from the cluster seriously deteriorates this agreement. It has been suggested that this could be related to the fact that UDGs are out-of-equilibrium objects whose stars have been heated by the cluster tides or that they recently fell onto the cluster on radial orbits, such that their velocity dispersion may not reflect the EFE at their instantaneous distance from the cluster center. Here, we simulate UDGs within the Coma cluster in MOND, using the Phantom of Ramses (\textsc{por}) code, and show that if UDGs are initially at equilibrium within the cluster, tides are not sufficient to increase their velocity dispersions to values as high as the observed ones. On the other hand, if they are on a first radial infall onto the cluster, they can keep high velocity dispersions without being destroyed until their first pericentric passage. We conclude that, without alterations such as a screening of the EFE in galaxy clusters or much higher baryonic masses than currently estimated, in the MOND context UDGs must be out-of-equilibrium objects on their first infall onto the cluster. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 10 Pages, 9 Figures, accepted by A&A

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

The Pristine survey -- XXVI. The very metal-poor Galaxy: Chemodynamics through the follow-up of the Pristine-Gaia synthetic catalogue

Authors: Akshara Viswanathan, Zhen Yuan, Anke Ardern-Arentsen, Else Starkenburg, Nicolas F. Martin, Kris Youakim, Rodrigo A. Ibata, Federico Sestito, Tadafumi Matsuno, Carlos Allende Prieto, Freya Barwell, Manuel Bayer, Amandine Doliva-Dolinsky, Emma Fernandez-Alvar, Pablo M. Galan-de Anta, Kiran Jhass, Nicolas Longeard, Jose Maria Arroyo-Polonio, Pol Massana, Martin Montelius, Samuel Rusterucci, Judith Santos, Guillaume F. Thomas, Sara Vitali, Wenbo Wu , et al. (5 additional authors not shown)

Abstract: The Pristine-\textit{Gaia} synthetic catalogue provides reliable photometric metallicities for $\sim$30 million FGK stars using the Pristine survey model and Gaia XP spectra. We perform the first low-to-medium-resolution spectroscopic follow-up of bright (G<15) and distant (up to 35 kpc) very and extremely metal-poor (V/EMP, [Fe/H]<-2.5) red giant branch stars from this. We use Isaac Newton Telesc… ▽ More The Pristine-\textit{Gaia} synthetic catalogue provides reliable photometric metallicities for $\sim$30 million FGK stars using the Pristine survey model and Gaia XP spectra. We perform the first low-to-medium-resolution spectroscopic follow-up of bright (G<15) and distant (up to 35 kpc) very and extremely metal-poor (V/EMP, [Fe/H]<-2.5) red giant branch stars from this. We use Isaac Newton Telescope/Intermediate Dispersion Spectrograph (INT/IDS) observations centred around the calcium triplet region ideal for V/EMP stars. We find that 76\% of our stars indeed have [Fe/H]<-2.5 with these inferred spectroscopic metallicities and only 3\% are outliers with [Fe/H] > -2.0. We report a success rate of 77\% and 38\% in finding stars with [Fe/H]<-2.5 and -3.0 respectively. This will allow for 10,000-20,000 homogeneously analysed EMP stars using the WEAVE survey follow-up of Pristine EMP candidates. We associate 20\%, 46\%, and 34\% of the stars to be confined to the disc plane, or to have inner and outer halo orbits, respectively. We also associate these V/EMP stars to known accretion events such as Gaia-Enceladus-Sausage (GES), LMS-1/Wukong, Thamnos, Helmi streams, Sagittarius, Sequoia, etc. For the stars that orbit close to the disc plane, we find that the prograde region with low vertical action is overdense with a significance of 4$σしぐま$ as compared to its retrograde counterpart. We also find three new (brightest) members of the most metal-poor stellar stream, C-19, one of which is 50$^\circ$ away from the main body of the stream. Our measured mean metallicity, velocity dispersion, and stream width are consistent with the literature, but our results favour a higher distance ($\sim$21.5 kpc) for the stream. We publish a catalogue (and 1D spectra) of 215 V/EMP stars from this spectroscopic follow-up and showcase the power of chemokinematic analysis of V/EMP end. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: Submitted to A&A. 17 pages (9 figures) + 3 pages (3 figures) in Appendix. Comments are very welcome! The catalogue and 1D spectra will be made available public after acceptance and before upon reasonable request to the first author

Scale-dependent local primordial non-Gaussianity as a solution to the $S_8$ tension

Authors: Clément Stahl, Benoit Famaey, Rodrigo Ibata, Oliver Hahn, Nicolas Martinet, Thomas Montandon

Abstract: For the last decade, several probes have pointed to a cosmological tension between the amplitude of density fluctuations extrapolated from the cosmic microwave background within the standard cosmological model and the one encapsulated by the $S_8$ parameter from large scale structure. The origin of this $S_8$ tension has not yet been elucidated and may hint at systematics in the data, unaccounted… ▽ More For the last decade, several probes have pointed to a cosmological tension between the amplitude of density fluctuations extrapolated from the cosmic microwave background within the standard cosmological model and the one encapsulated by the $S_8$ parameter from large scale structure. The origin of this $S_8$ tension has not yet been elucidated and may hint at systematics in the data, unaccounted effects from baryonic physics, or new physics beyond the standard model of cosmology. Baryonic physics may in principle provide a nonlinear solution to the tension by suppressing the matter power spectrum more strongly on nonlinear scales than is traditionally assumed. Such a solution would not worsen the Hubble tension, contrary to many other proposed solutions to the $S_8$ tension. However, no realistic baryonic feedback in hydrodynamical simulations provides the needed suppression as a function of redshift. Here, we point out that a scale-dependence of local-type primordial non-Gaussianities (PNG), with significant PNG at scales of a few Mpc, can provide the needed suppression, since such PNG can suppress the power spectrum at slightly larger scales than baryons do. We demonstrate this by devising collisionless numerical simulations of structure formation in boxes of 0.5 Gpc/h with scale-dependent local-type PNG. Our simple models show that, as a proof of principle, scale-dependent PNG, with a Gaussian random field for primodial density fluctuations on large scales and $f_{\rm NL} \simeq -300$ at $\lesssim 10$ Mpc scales, together with state-of-the-art baryonification of the matter power spectrum, can in principle solve the $S_8$ tension. The $S_8$ tension would then be a smoking-gun of non-trivial inflationary physics. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 12 pages, 7 figures, comments welcome :)

The Gaia RVS benchmark stars II. A sample of stars selected for their Gaia high radial velocity

Authors: E. Caffau, D. Katz, A. Gómez, P. Bonifacio, R. Lallement, P. Sartoretti, L. Sbordone, M. Spite, A. Mucciarelli, R. Ibata, L. Chemin, F. Thévenin, P. Panuzzo, N. Leclerc, P. François, H. -G. Ludwig, L. Monaco, M. Haywood, C. Soubiran

Abstract: The Gaia satellite has already provided the astronomical community with three data releases, and the Radial Velocity Spectrometer (RVS) on board Gaia has provided the radial velocity for 33 million stars. When deriving the radial velocity from the RVS spectra, several stars are measured to have large values. To verify the credibility of these measurements, we selected some bright stars with the mo… ▽ More The Gaia satellite has already provided the astronomical community with three data releases, and the Radial Velocity Spectrometer (RVS) on board Gaia has provided the radial velocity for 33 million stars. When deriving the radial velocity from the RVS spectra, several stars are measured to have large values. To verify the credibility of these measurements, we selected some bright stars with the modulus of radial velocity in excess of 500\ to be observed with SOPHIE at OHP and UVES at VLT. This paper is devoted to investigating the chemical composition of the stars observed with UVES. We derived atmospheric parameters using Gaia photometry and parallaxes, and we performed a chemical analysis using the code. We find that the sample consists of metal-poor stars, although none have extremely low metallicities. The abundance patterns match what has been found in other samples of metal-poor stars selected irrespective of their radial velocities. We highlight the presence of three stars with low Cu and Zn abundances that are likely descendants of pair-instability supernovae. Two stars are apparently younger than 1\,Ga, and their masses exceed twice the turn-off mass of metal-poor populations. This makes it unlikely that they are blue stragglers because it would imply they formed from triple or multiple systems. We suggest instead that they are young metal-poor stars accreted from a dwarf galaxy. Finally, we find that the star RVS721 is associated with the Gjoll stream, which itself is associated with the Globular Cluster NGC\,3201. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: Astronomy and Astrophysics - A\&A, In press

From inflation to dark matter halo profiles: the impact of primordial non-Gaussianities on the central density cusp

Authors: Clément Stahl, Nicolas Mai, Benoit Famaey, Yohan Dubois, Rodrigo Ibata

Abstract: It has recently been shown that local primordial non-Gaussianities (PNG) with significant amplitude ($|f_{\rm NL}| \sim 1000$), at small (Mpc) scales, can help in forming simulated galaxies with more disky baryonic kinematics than in the Gaussian case, while generating matter power spectra that can differ by up to 20% from the Gaussian case at non-linear scales. Here, we explore in detail the cons… ▽ More It has recently been shown that local primordial non-Gaussianities (PNG) with significant amplitude ($|f_{\rm NL}| \sim 1000$), at small (Mpc) scales, can help in forming simulated galaxies with more disky baryonic kinematics than in the Gaussian case, while generating matter power spectra that can differ by up to 20% from the Gaussian case at non-linear scales. Here, we explore in detail the consequences of such small-scale PNG on the dark matter halo profiles. We show in particular that, for negative $f_{\rm NL}$, dark matter halos formed in collisionless simulations are not always well described by the traditional Navarro-Frenk-White (NFW) profiles, as supported by their sparsity distribution. We conclude that NFW profiles are not as clear attractors for the density profiles of dark matter halos in the presence of PNG than in the case of a Gaussian contrast density field. We show how alternatives to the NFW profile can describe halos both in the Gaussian and non-Gaussian cases. From the combination of our sparsity analysis and the quality of the adjustments of the density profiles with a minimal extension to NFW, we conclude that $z=1$ halos carry the most interesting information about PNG △ Less

Submitted 7 April, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: 21 pages, 12 figures, JCAP accepted

Physical Symbolic Optimization

Authors: Wassim Tenachi, Rodrigo Ibata, Foivos I. Diakogiannis

Abstract: We present a framework for constraining the automatic sequential generation of equations to obey the rules of dimensional analysis by construction. Combining this approach with reinforcement learning, we built $Φふぁい$-SO, a Physical Symbolic Optimization method for recovering analytical functions from physical data leveraging units constraints. Our symbolic regression algorithm achieves state-of-the-a… ▽ More We present a framework for constraining the automatic sequential generation of equations to obey the rules of dimensional analysis by construction. Combining this approach with reinforcement learning, we built $Φふぁい$-SO, a Physical Symbolic Optimization method for recovering analytical functions from physical data leveraging units constraints. Our symbolic regression algorithm achieves state-of-the-art results in contexts in which variables and constants have known physical units, outperforming all other methods on SRBench's Feynman benchmark in the presence of noise (exceeding 0.1%) and showing resilience even in the presence of significant (10%) levels of noise. △ Less

Submitted 6 December, 2023; originally announced December 2023.

Comments: 6 pages, 2 figures, 1 table. Accepted to NeurIPS 2023, Machine Learning for Physical Sciences workshop

Class Symbolic Regression: Gotta Fit 'Em All

Authors: Wassim Tenachi, Rodrigo Ibata, Thibaut L. François, Foivos I. Diakogiannis

Abstract: We introduce 'Class Symbolic Regression' (Class SR) a first framework for automatically finding a single analytical functional form that accurately fits multiple datasets - each realization being governed by its own (possibly) unique set of fitting parameters. This hierarchical framework leverages the common constraint that all the members of a single class of physical phenomena follow a common go… ▽ More We introduce 'Class Symbolic Regression' (Class SR) a first framework for automatically finding a single analytical functional form that accurately fits multiple datasets - each realization being governed by its own (possibly) unique set of fitting parameters. This hierarchical framework leverages the common constraint that all the members of a single class of physical phenomena follow a common governing law. Our approach extends the capabilities of our earlier Physical Symbolic Optimization ($Φふぁい$-SO) framework for Symbolic Regression, which integrates dimensional analysis constraints and deep reinforcement learning for unsupervised symbolic analytical function discovery from data. Additionally, we introduce the first Class SR benchmark, comprising a series of synthetic physical challenges specifically designed to evaluate such algorithms. We demonstrate the efficacy of our novel approach by applying it to these benchmark challenges and showcase its practical utility for astrophysics by successfully extracting an analytic galaxy potential from a set of simulated orbits approximating stellar streams. △ Less

Submitted 17 June, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 10 pages, 4 figures, 1 table. Accepted for publication at ApJL

Charting the Galactic acceleration field II. A global mass model of the Milky Way from the STREAMFINDER Atlas of Stellar Streams detected in Gaia DR3

Authors: Rodrigo Ibata, Khyati Malhan, Wassim Tenachi, Anke Ardern-Arentsen, Michele Bellazzini, Paolo Bianchini, Piercarlo Bonifacio, Elisabetta Caffau, Foivos Diakogiannis, Raphael Errani, Benoit Famaey, Salvatore Ferrone, Nicolas Martin, Paola di Matteo, Giacomo Monari, Florent Renaud, Else Starkenburg, Guillaume Thomas, Akshara Viswanathan, Zhen Yuan

Abstract: We present an atlas and follow-up spectroscopic observations of 87 thin stream-like structures detected with the STREAMFINDER algorithm in Gaia DR3, of which 29 are new discoveries. Here we focus on using these streams to refine mass models of the Galaxy. Fits with a double power law halo with the outer power law slope set to $-βべーた_h=3$ yield an inner power law slope $-γがんま_h=0.97^{+0.17}_{-0.21}$, a s… ▽ More We present an atlas and follow-up spectroscopic observations of 87 thin stream-like structures detected with the STREAMFINDER algorithm in Gaia DR3, of which 29 are new discoveries. Here we focus on using these streams to refine mass models of the Galaxy. Fits with a double power law halo with the outer power law slope set to $-βべーた_h=3$ yield an inner power law slope $-γがんま_h=0.97^{+0.17}_{-0.21}$, a scale radius of $r_{0, h}=14.7^{+4.7}_{-1.0}$ kpc, a halo density flattening $q_{m, h}=0.75\pm0.03$, and a local dark matter density of $ρろー_{h, \odot}=0.0114\pm0.0007 {\rm M_\odot pc^{-3}}$. Freeing $βべーた$ yields $βべーた=2.53^{+0.42}_{-0.16}$, but this value is heavily influenced by our chosen virial mass limit. The stellar disks are found to have a combined mass of $4.20^{+0.44}_{-0.53}\times10^{10} {\rm M_\odot}$, with the thick disk contributing $12.4\pm0.7$\% to the local stellar surface density. The scale length of the thin and thick disks are $2.17^{+0.18}_{-0.08}$ kpc and $1.62^{+0.72}_{-0.13}$ kpc, respectively, while their scale heights are $0.347^{+0.007}_{-0.010}$ kpc and $0.86^{+0.03}_{-0.02}$ kpc, respectively. The virial mass of the favored model is $M_{200}=1.09^{+0.19}_{-0.14}\times 10^{12} {\rm M_\odot}$, while the mass inside of 50 kpc is $M_{R<50}=0.46\pm0.03\times 10^{12} {\rm M_\odot}$. We introduce the Large Magellanic Cloud (LMC) into the derived potential models, and fit the "Orphan" stream therein, finding a mass for the LMC that is consistent with recent estimates. Some highlights of the atlas include the nearby trailing arm of $ωおめが$-Cen, and a nearby very metal-poor stream that was once a satellite of the Sagittarius dwarf galaxy. Finally, we unambiguously detect a hot component around the GD-1 stream, consistent with it having been tidally pre-processed within its own DM subhalo. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 43 pages, 26 figures, submitted to ApJ

Exploring the impact of a rapidly decelerating bar on transforming bulge orbits into disc-like orbits

Authors: Chengdong Li, Zhen Yuan, Giacomo Monari, Nicolas F. Martin, Arnaud Siebert, Benoit Famaey, Rimpei Chiba, Georges Kordopatis, Rodrigo A. Ibata, Vanessa Hill

Abstract: The most metal-poor tail of the Milky Way ([Fe/H] $\leq$ $-$2.5) contains a population of stars with very prograde planar orbits, which is puzzling in both their origin and evolution. A possible scenario is that they are shepherded by the bar from the inner Galaxy, where many of the old and low-metallicity stars in the Galaxy are located. To investigate this scenario, we use test-particle simulati… ▽ More The most metal-poor tail of the Milky Way ([Fe/H] $\leq$ $-$2.5) contains a population of stars with very prograde planar orbits, which is puzzling in both their origin and evolution. A possible scenario is that they are shepherded by the bar from the inner Galaxy, where many of the old and low-metallicity stars in the Galaxy are located. To investigate this scenario, we use test-particle simulations with an axisymmetric background potential plus a central bar model. The test particles are generated by an extended distribution function (EDF) model based on the observational constraints of bulge stars. According to the simulation results, a bar with constant pattern speed cannot help bring stars from the bulge to the solar vicinity. In contrast, when the model includes a rapidly decelerating bar, some bulge stars can gain rotation and move outwards as they are trapped in the co-rotation regions of the bar. The resulting distribution of shepherded stars heavily depends on the present-day azimuthal angle between the bar and the Sun. The majority of the low-metallicity bulge stars driven outwards are distributed in the fourth quadrant of the Galaxy with respect to the Sun, and about 10$\%$ of them are within 6 kpc from us. Our experiments indicate that the decelerating bar perturbation can be a contributing process to explain part of the most metal-poor stars with prograde planar orbits seen in the Solar neighborhood but is unlikely to be the dominant one. △ Less

Submitted 26 February, 2024; v1 submitted 26 November, 2023; originally announced November 2023.

Comments: 12 pages, 11 figures, re-submitted to A&A with revisions

Micro galaxies in LCDM

Authors: Raphaël Errani, Rodrigo Ibata, Julio F. Navarro, Jorge Peñarrubia, Matthew G. Walker

Abstract: A fundamental prediction of the Lambda Cold Dark Matter (LCDM) cosmology is the centrally divergent cuspy density profile of dark matter haloes. Density cusps render CDM haloes resilient to tides, and protect dwarf galaxies embedded in them from full tidal disruption. The hierarchical assembly history of the Milky Way may therefore give rise to a population of micro galaxies; i.e., heavily strippe… ▽ More A fundamental prediction of the Lambda Cold Dark Matter (LCDM) cosmology is the centrally divergent cuspy density profile of dark matter haloes. Density cusps render CDM haloes resilient to tides, and protect dwarf galaxies embedded in them from full tidal disruption. The hierarchical assembly history of the Milky Way may therefore give rise to a population of micro galaxies; i.e., heavily stripped remnants of early accreted satellites, which may reach arbitrarily low luminosity. Assuming that the progenitor systems are dark matter dominated, we use an empirical formalism for tidal stripping to predict the evolution of the luminosity, size, and velocity dispersion of such remnants, tracing their tidal evolution across multiple orders of magnitude in mass and size. The evolutionary tracks depend sensitively on the progenitor distribution of stellar binding energies. We explore three cases that likely bracket most realistic models of dwarf galaxies: one where the energy distribution of the most tightly bound stars follows that of the dark matter, and two where stars are defined by either an exponential density or surface brightness profile. The tidal evolution in the size-velocity dispersion plane is quite similar for these three models, although their remnants may differ widely in luminosity. Micro galaxies are therefore best distinguished from globular clusters by the presence of dark matter; either directly, by measuring their velocity dispersion, or indirectly, by examining their tidal resilience. Our work highlights the need for further theoretical and observational constraints on the stellar energy distribution in dwarf galaxies. △ Less

Submitted 20 May, 2024; v1 submitted 24 November, 2023; originally announced November 2023.

Comments: updated to match accepted version

Journal ref: ApJ 968 89 (2024)

Could very low-metallicity stars with rotation-dominated orbits have been shepherded by the bar?

Authors: Zhen Yuan, Chengdong Li, Nicolas F. Martin, Giacomo Monari, Benoit Famaey, Arnaud Siebert, Anke Ardern-Arentsen, Federico Sestito, Guillaume F. Thomas, Vanessa Hill, Rodrigo A. Ibata, Georges Kordopatis, Else Starkenburg, Akshara Viswanathan

Abstract: The most metal-poor stars (e.g. [Fe/H] $\leq-2.5$) are the ancient fossils from the early assembly epoch of our Galaxy, very likely before the formation of the thick disc. Recent studies have shown that a non-negligible fraction of them have prograde planar orbits, which makes their origin a puzzle. It has been suggested that a later-formed rotating bar could have driven these old stars from the i… ▽ More The most metal-poor stars (e.g. [Fe/H] $\leq-2.5$) are the ancient fossils from the early assembly epoch of our Galaxy, very likely before the formation of the thick disc. Recent studies have shown that a non-negligible fraction of them have prograde planar orbits, which makes their origin a puzzle. It has been suggested that a later-formed rotating bar could have driven these old stars from the inner Galaxy outward, and transformed their orbits to be more rotation-dominated. However, it is not clear if this mechanism can explain these stars as observed in the solar neighborhood. In this paper, we explore the possibility of this scenario by tracing these stars backwards in an axisymmetric Milky Way potential with a bar perturber. We integrate their orbits backward for 6 Gyr under two bar models: one with a constant pattern speed and another one with a decelerating speed. Our experiments show that, under the constantly-rotating bar model, the stars of interest are little affected by the bar and cannot have been shepherded from a spheroidal inner Milky Way to their current orbits. In the extreme case of a rapidly decelerating bar, some of the very metal-poor stars on planar and prograde orbits can be brought from the inner Milky Way, but $\sim90\%$ of them were nevertheless already rotation-dominated ($J_φふぁい$ $\geq$ 1000 km s$^{-1}$ kpc) 6 Gyr ago. The chance of these stars having started with spheroid-like orbits with small rotation ($J_φふぁい$ $\lesssim$ 600 km s$^{-1}$ kpc) is very low ($<$ 3$\%$). We therefore conclude that, within the solar neighborhood, the bar is unlikely to have shepherded a significant fraction of inner Galaxy spheroid stars to produce the overdensity of stars on prograde, planar orbits that is observed today. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: submitted to A&A, comments are welcome

Strong constraints on the gravitational law from $Gaia$ DR3 wide binaries

Authors: Indranil Banik, Charalambos Pittordis, Will Sutherland, Benoit Famaey, Rodrigo Ibata, Steffen Mieske, Hongsheng Zhao

Abstract: We test Milgromian dynamics (MOND) using wide binary stars (WBs) with separations of $2-30$ kAU. Locally, the WB orbital velocity in MOND should exceed the Newtonian prediction by $\approx 20\%$ at asymptotically large separations given the Galactic external field effect (EFE). We investigate this with a detailed statistical analysis of $Gaia$ DR3 data on 8611 WBs within 250 pc of the Sun. Orbits… ▽ More We test Milgromian dynamics (MOND) using wide binary stars (WBs) with separations of $2-30$ kAU. Locally, the WB orbital velocity in MOND should exceed the Newtonian prediction by $\approx 20\%$ at asymptotically large separations given the Galactic external field effect (EFE). We investigate this with a detailed statistical analysis of $Gaia$ DR3 data on 8611 WBs within 250 pc of the Sun. Orbits are integrated in a rigorously calculated gravitational field that directly includes the EFE. We also allow line of sight contamination and undetected close binary companions to the stars in each WB. We interpolate between the Newtonian and Milgromian predictions using the parameter $αあるふぁ_{\rm{grav}}$, with 0 indicating Newtonian gravity and 1 indicating MOND. Directly comparing the best Newtonian and Milgromian models reveals that Newtonian dynamics is preferred at $19σしぐま$ confidence. Using a complementary Markov Chain Monte Carlo analysis, we find that $αあるふぁ_{\rm{grav}} = -0.021^{+0.065}_{-0.045}$, which is fully consistent with Newtonian gravity but excludes MOND at $16σしぐま$ confidence. This is in line with the similar result of Pittordis and Sutherland using a somewhat different sample selection and less thoroughly explored population model. We show that although our best-fitting model does not fully reproduce the observations, an overwhelmingly strong preference for Newtonian gravity remains in a considerable range of variations to our analysis. Adapting the MOND interpolating function to explain this result would cause tension with rotation curve constraints. We discuss the broader implications of our results in light of other works, concluding that MOND must be substantially modified on small scales to account for local WBs. △ Less

Submitted 24 November, 2023; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: 48 pages, 28 figures. Published in MNRAS in this form

Journal ref: Monthly Notices of the Royal Astronomical Society, volume 527, issue 3, pages 4573 - 4615 (2024)

Impact of the Galactic bar on tidal streams within the Galactic disc: The case of the tidal stream of the Hyades

Authors: Guillaume F. Thomas, Benoit Famaey, Giacomo Monari, Chervin F. P. Laporte, Rodrigo Ibata, Patrick de Laverny, Vanessa Hill, Christian Boily

Abstract: Tidal streams of disrupted clusters are routinely detected in the halo of the Milky Way. It was recently shown that tidal streams of open clusters can now also be detected within the disc. In this work, we highlight the fact that tidal streams provide a powerful new diagnostic of the non-axisymmetric disc potential and may provide a new constraint on the pattern speed of the Galactic bar. In parti… ▽ More Tidal streams of disrupted clusters are routinely detected in the halo of the Milky Way. It was recently shown that tidal streams of open clusters can now also be detected within the disc. In this work, we highlight the fact that tidal streams provide a powerful new diagnostic of the non-axisymmetric disc potential and may provide a new constraint on the pattern speed of the Galactic bar. In particular, we show how the stream-orbit misalignment for an open cluster on a quasi-circular orbit in the solar vicinity varies as a function of the position w.r.t the bar resonances. The angular shift rises beyond corotation, reaching values as high as $30^\circ$ close to the outer Lindblad resonance (OLR), then dropping again and reversing its sign beyond the OLR. We applied this mechanism to the recently detected the Hyades stream. We note that the stream would be very similar when taking a potential with no bar or with a fast pattern speed of 55 km.s$^{-1}$ kpc$^{-1}$. However, we find that the stream is different than previously detected when adopting a potential with a bar pattern speed of $39$ km.s$^{-1}$ kpc$^{-1}$. Previously detected Hyades candidate members would, on the other hand, favour a barless or a fast bar galaxy. Interestingly, the previously reported asymmetry in star counts within the leading and trailing tails of the Hyades tidal stream persists in all cases. Our study conclusively demonstrates that the effect of disc non-axisymmetries cannot be neglected when searching for tidal streams of open clusters and that current candidate members of the Hyades stream should not be trusted beyond a distance of 200 pc from the cluster. Moreover, our study allows for ideal targets to be provided for high-resolution spectroscopy follow-ups, which will enable conclusive identifications of the Hyades stream track and provide novel independent constraints on the bar pattern speed in the MW. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 16 pages, 9 figures, accepted for publication in A&A

The Pristine survey -- XXIII. Data Release 1 and an all-sky metallicity catalogue based on Gaia DR3 BP/RP spectro-photometry

Authors: Nicolas F. Martin, Else Starkenburg, Zhen Yuan, Morgan Fouesneau, Anke Arentsen, Francesca De Angeli, Felipe Gran, Martin Montelius, René Andrae, Michele Bellazzini, Paolo Montegriffo, Anna F. Esselink, Hanyuan Zhang, Kim A. Venn, Akshara Viswanathan, David S. Aguado, Giuseppina Battaglia, Manuel Bayer, Piercarlo Bonifacio, Elisabetta Caffau, Patrick Côté, Raymond Carlberg, Sébastien Fabbro, Emma Fernández Alvar, Jonay I. González Hernández , et al. (15 additional authors not shown)

Abstract: We use the spectro-photometric information of ~219 million stars from Gaia's DR3 to calculate synthetic, narrow-band, metallicity-sensitive CaHK magnitudes that mimic the observations of the Pristine survey, a survey of photometric metallicities of Milky Way stars that has been mapping more than 6,500 deg^2 of the northern sky with the CFHT since 2015. These synthetic magnitudes are used for an ab… ▽ More We use the spectro-photometric information of ~219 million stars from Gaia's DR3 to calculate synthetic, narrow-band, metallicity-sensitive CaHK magnitudes that mimic the observations of the Pristine survey, a survey of photometric metallicities of Milky Way stars that has been mapping more than 6,500 deg^2 of the northern sky with the CFHT since 2015. These synthetic magnitudes are used for an absolute re-calibration of the deeper Pristine photometry and, combined with broadband Gaia information, synthetic and Pristine CaHK magnitudes are used to estimate photometric metallicities over the whole sky. The resulting metallicity catalogue is accurate down to [Fe/H]~-3.5 and is particularly suited for the exploration of the metal-poor Milky Way ([Fe/H]<-1.0). We make available here the catalogue of synthetic CaHK_syn magnitudes for all stars with BP/RP information in Gaia DR3, as well as an associated catalogue of more than ~30 million photometric metallicities for high S/N FGK stars. This paper further provides the first public DR of the Pristine catalogue in the form of higher quality recalibrated Pristine CaHK magnitudes and photometric metallicities for all stars in common with the BP/RP information in Gaia DR3. We demonstrate that, when available, the much deeper Pristine data greatly enhances the quality of the derived metallicities, in particular at the faint end of the catalogue (G_BP>16). Combined, both catalogues include more than 2 million metal-poor star candidates as well as more than 200,000 and ~8,000 very and extremely metal-poor candidates. Finally, we show that these metallicity catalogues can be used efficiently, among other applications, for Galactic archaeology, to hunt for the most metal-poor stars, and to study how the structure of the Milky Way varies with metallicity, from the flat distribution of disk stars to the spheroid-shaped metal-poor halo. (Shortened) △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 30 pages, 24 figures, submitted to A&A. First two authors are co-first author. The CaHK photometry catalogue and the two photometric metallicity catalogues are available, before acceptance, as large compressed csv files at: https://seafile.unistra.fr/d/ee0c0f05719d4368bcbb/

Hydrodynamical simulations of galaxy formation with non-Gaussian initial conditions

Authors: Clément Stahl, Yohan Dubois, Benoit Famaey, Oliver Hahn, Rodrigo Ibata, Katarina Kraljic, Thomas Montandon

Abstract: Collisionless simulations of structure formation with significant local primordial non-Gaussianities at Mpc scales have shown that a non-Gaussian tail favouring underdensities, with a negative $f_{\rm NL}$ parameter, can significantly change the merging history of galaxy-sized dark matter halos, which then typically assemble later than in vanilla $Λらむだ$CDM. Moreover, such a small-scale negative… ▽ More Collisionless simulations of structure formation with significant local primordial non-Gaussianities at Mpc scales have shown that a non-Gaussian tail favouring underdensities, with a negative $f_{\rm NL}$ parameter, can significantly change the merging history of galaxy-sized dark matter halos, which then typically assemble later than in vanilla $Λらむだ$CDM. Moreover, such a small-scale negative $f_{\rm NL}$ could have interesting consequences for the cosmological $S_8$ tension. Here, we complement our previous work on collisionless simulations with new hydrodynamical simulations of galaxy formation in boxes of 30 Mpc/$h$, using the {\sc RAMSES} code. In particular, we show that all feedback prescriptions being otherwise identical, simulations with a negative $f_{\rm NL} \sim -1000$ on small scales, hence forming galaxies a bit later than in vanilla $Λらむだ$CDM, allow to form simulated galaxies with more disky kinematics than in the vanilla case. Therefore, such small-scale primordial non-Gaussianities could potentially help alleviate, simultaneously, tensions in cosmology and galaxy formation. These hydrodynamical simulations on small scales will need to be complemented with larger box simulations with scale-dependent non-Gaussianities, to statistically confirm these trends and explore their observational consequences in further detail. △ Less

Submitted 6 September, 2023; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 15 pages, 8 figures, JCAP accepted

Do old globular clusters in low mass galaxies disprove modified gravity?

Authors: Michal Bílek, Hongsheng Zhao, Benoit Famaey, Srikanth T. Nagesh, Françoise Combes, Oliver Müller, Michael Hilker, Pavel Kroupa, Rodrigo Ibata

Abstract: The controversy "dark matter vs. modified gravity" constitutes a major topic of discussion. It was proposed that dynamical friction could be used to discriminate between the two alternatives. Analytic calculations indicate that, with modified gravity, globular clusters (GCs) of low-mass galaxies experience much stronger dynamical friction than in the equivalent system with Newtonian gravity and da… ▽ More The controversy "dark matter vs. modified gravity" constitutes a major topic of discussion. It was proposed that dynamical friction could be used to discriminate between the two alternatives. Analytic calculations indicate that, with modified gravity, globular clusters (GCs) of low-mass galaxies experience much stronger dynamical friction than in the equivalent system with Newtonian gravity and dark matter. As a result, in modified gravity the old GCs of low mass galaxies should have already settled in the centers of the galaxies. This is not observed. Here we report on our efforts to verify the analytic results by self-consistent simulations with the MOND-type (modified Newtonian dynamics) gravity. The core stalling mechanism, that was not considered in the analytic calculations, prevents GCs to settle in centers of ultra-diffuse galaxies. For isolated dwarf galaxies, which are gas-rich objects, supernova explosions prevent the GCs from settling. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 4 pages, 2 figures, accepted for publication in the proceedings of IAU Symposium 379: Dynamical Masses of Local Group Galaxies

Disk galaxies are self-similar: the universality of the HI-to-Halo mass ratio for isolated disks

Authors: Marie Korsaga, Benoit Famaey, Jonathan Freundlich, Lorenzo Posti, Rodrigo Ibata, Christian Boily, Katarina Kraljic, D. Esparza-Arredondo, C. Ramos Almeida, Jean Koulidiati

Abstract: Observed scaling relations in galaxies between baryons and dark matter global properties are key to shed light on the process of galaxy formation and on the nature of dark matter. Here, we study the scaling relation between the neutral hydrogen (HI) and dark matter mass in isolated rotationally-supported disk galaxies at low redshift. We first show that state-of-the-art galaxy formation simulation… ▽ More Observed scaling relations in galaxies between baryons and dark matter global properties are key to shed light on the process of galaxy formation and on the nature of dark matter. Here, we study the scaling relation between the neutral hydrogen (HI) and dark matter mass in isolated rotationally-supported disk galaxies at low redshift. We first show that state-of-the-art galaxy formation simulations predict that the HI-to-dark halo mass ratio decreases with stellar mass for the most massive disk galaxies. We then infer dark matter halo masses from high-quality rotation curve data for isolated disk galaxies in the local Universe, and report on the actual universality of the HI-to-dark halo mass ratio for these observed galaxies. This scaling relation holds for disks spanning a range of 4 orders of magnitude in stellar mass and 3 orders of magnitude in surface brightness. Accounting for the diversity of rotation curve shapes in our observational fits decreases the scatter of the HI-to-dark halo mass ratio while keeping it constant. This finding extends the previously reported discrepancy for the stellar-to-halo mass relation of massive disk galaxies within galaxy formation simulations to the realm of neutral atomic gas. Our result reveals that isolated galaxies with regularly rotating extended HI disks are surprisingly self-similar up to high masses, which hints at mass-independent self-regulation mechanisms that have yet to be fully understood. △ Less

Submitted 5 July, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 14 pages, 4 figures. Accepted for publication in ApJL

Reionisation time fields reconstruction from 21 cm signal maps

Authors: Julien Hiegel, Emilie Thélie, Dominique Aubert, Jonathan Chardin, Nicolas Gillet, Pierre Galois, Nicolas Mai, Pierre Ocvirk, Rodrigo Ibata

Abstract: During the Epoch of reionisation, the intergalactic medium is reionised by the UV radiation from the first generation of stars and galaxies. One tracer of the process is the 21 cm line of hydrogen that will be observed by the Square Kilometre Array (SKA) at low frequencies, thus imaging the distribution of ionised and neutral regions and their evolution. To prepare for these upcoming observations,… ▽ More During the Epoch of reionisation, the intergalactic medium is reionised by the UV radiation from the first generation of stars and galaxies. One tracer of the process is the 21 cm line of hydrogen that will be observed by the Square Kilometre Array (SKA) at low frequencies, thus imaging the distribution of ionised and neutral regions and their evolution. To prepare for these upcoming observations, we investigate a deep learning method to predict from 21 cm maps the reionisation time field (treion(r)), i.e. the time at which each location has been reionised. treion(r) encodes the propagation of ionisation fronts in a single field, gives access to times of local reionisation or to the extent of the radiative reach of early sources. Moreover it gives access to the time evolution of ionisation on the plane of sky, when such evolution is usually probed along the line-of-sight direction. We trained a convolutional neural network (CNN) using simulated 21 cm maps and reionisation times fields produced by the simulation code 21cmFAST . We also investigate the performance of the CNN when adding instrumental effects. Globally, we find that without instrumental effects the 21 cm maps can be used to reconstruct the associated reionisation times field in a satisfying manner: the quality of the reconstruction is dependent on the redshift at which the 21 cm observation is being made and in general it is found that small scale (<10cMpc/h) features are smoothed in the reconstructed field, while larger scale features are well recovered. When instrumental effects are included, the scale dependance of reconstruction is even further pronounced, with significant smoothing on small and intermediate scales. △ Less

Submitted 2 October, 2023; v1 submitted 2 July, 2023; originally announced July 2023.

Journal ref: A&A 679, A125 (2023)

RomAndromeda: The Roman Survey of the Andromeda Halo

Authors: Arjun Dey, Joan Najita, Carrie Filion, Jiwon Jesse Han, Sarah Pearson, Rosemary Wyse, Adrien C. R. Thob, Borja Anguiano, Miranda Apfel, Magda Arnaboldi, Eric F. Bell, Leandro Beraldo e Silva, Gurtina Besla, Aparajito Bhattacharya, Souradeep Bhattacharya, Vedant Chandra, Yumi Choi, Michelle L. M. Collins, Emily C. Cunningham, Julianne J. Dalcanton, Ivanna Escala, Hayden R. Foote, Annette M. N. Ferguson, Benjamin J. Gibson, Oleg Y. Gnedin , et al. (28 additional authors not shown)

Abstract: As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δでるたt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the ha… ▽ More As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($Δでるたt\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$σしぐま$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own. △ Less

Submitted 21 June, 2023; originally announced June 2023.

Comments: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

NANCY: Next-generation All-sky Near-infrared Community surveY

Authors: Jiwon Jesse Han, Arjun Dey, Adrian M. Price-Whelan, Joan Najita, Edward F. Schlafly, Andrew Saydjari, Risa H. Wechsler, Ana Bonaca, David J Schlegel, Charlie Conroy, Anand Raichoor, Alex Drlica-Wagner, Juna A. Kollmeier, Sergey E. Koposov, Gurtina Besla, Hans-Walter Rix, Alyssa Goodman, Douglas Finkbeiner, Abhijeet Anand, Matthew Ashby, Benedict Bahr-Kalus, Rachel Beaton, Jayashree Behera, Eric F. Bell, Eric C Bellm , et al. (184 additional authors not shown)

Abstract: The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GAL… ▽ More The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Submitted to the call for white papers for the Roman Core Community Survey (June 16th, 2023), and to the Bulletin of the AAS

An end-to-end strategy for recovering a free-form potential from a snapshot of stellar coordinates

Authors: Wassim Tenachi, Rodrigo Ibata, Foivos I. Diakogiannis

Abstract: New large observational surveys such as Gaia are leading us into an era of data abundance, offering unprecedented opportunities to discover new physical laws through the power of machine learning. Here we present an end-to-end strategy for recovering a free-form analytical potential from a mere snapshot of stellar positions and velocities. First we show how auto-differentiation can be used to capt… ▽ More New large observational surveys such as Gaia are leading us into an era of data abundance, offering unprecedented opportunities to discover new physical laws through the power of machine learning. Here we present an end-to-end strategy for recovering a free-form analytical potential from a mere snapshot of stellar positions and velocities. First we show how auto-differentiation can be used to capture an agnostic map of the gravitational potential and its underlying dark matter distribution in the form of a neural network. However, in the context of physics, neural networks are both a plague and a blessing as they are extremely flexible for modeling physical systems but largely consist in non-interpretable black boxes. Therefore, in addition, we show how a complementary symbolic regression approach can be used to open up this neural network into a physically meaningful expression. We demonstrate our strategy by recovering the potential of a toy isochrone system. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: 4 pages, 2 figures. Accepted for publication in the International Astronomical Union Proceedings Series

The Hubble Space Telescope Survey of M31 Satellite Galaxies II. The Star Formation Histories of Ultra-Faint Dwarf Galaxies

Authors: A. Savino, D. R. Weisz, E. D. Skillman, A. Dolphin, A. A. Cole, N. Kallivayalil, A. Wetzel, J. Anderson, G. Besla, M. Boylan-Kolchin, T. M. Brown, J. S. Bullock, M. L. M. Collins, M. C. Cooper, A. J. Deason, A. L. Dotter, M. Fardal, A. M. N. Ferguson, T. K. Fritz, M. C. Geha, K. M. Gilbert, P. Guhathakurta, R. Ibata, M. J. Irwin, M. Jeon , et al. (12 additional authors not shown)

Abstract: We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find th… ▽ More We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50\% of their stellar mass by $z=5$ (12.6~Gyr ago), similar to known UFDs around the MW, but that 10-40\% of their stellar mass formed at later times. We uncover one remarkable UFD, \A{XIII}, which formed only 10\% of its stellar mass by $z=5$, and 75\% in a rapid burst at $z\sim2-3$, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This "young" UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least massive MW UFDs ($M_*(z=5) \lesssim 5\times10^4 M_{\odot}$) are likely quenched by reionization, whereas more massive M31 UFDs ($M_*(z=5) \gtrsim 10^5 M_{\odot}$) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs. △ Less

Submitted 13 September, 2023; v1 submitted 22 May, 2023; originally announced May 2023.

Comments: 18 pages, 14 figures, 5 appendices, accepted for publication in ApJ

Deep symbolic regression for physics guided by units constraints: toward the automated discovery of physical laws

Authors: Wassim Tenachi, Rodrigo Ibata, Foivos I. Diakogiannis

Abstract: Symbolic Regression is the study of algorithms that automate the search for analytic expressions that fit data. While recent advances in deep learning have generated renewed interest in such approaches, the development of symbolic regression methods has not been focused on physics, where we have important additional constraints due to the units associated with our data. Here we present $Φふぁい$-SO, a P… ▽ More Symbolic Regression is the study of algorithms that automate the search for analytic expressions that fit data. While recent advances in deep learning have generated renewed interest in such approaches, the development of symbolic regression methods has not been focused on physics, where we have important additional constraints due to the units associated with our data. Here we present $Φふぁい$-SO, a Physical Symbolic Optimization framework for recovering analytical symbolic expressions from physics data using deep reinforcement learning techniques by learning units constraints. Our system is built, from the ground up, to propose solutions where the physical units are consistent by construction. This is useful not only in eliminating physically impossible solutions, but because the "grammatical" rules of dimensional analysis restrict enormously the freedom of the equation generator, thus vastly improving performance. The algorithm can be used to fit noiseless data, which can be useful for instance when attempting to derive an analytical property of a physical model, and it can also be used to obtain analytical approximations to noisy data. We test our machinery on a standard benchmark of equations from the Feynman Lectures on Physics and other physics textbooks, achieving state-of-the-art performance in the presence of noise (exceeding 0.1%) and show that it is robust even in the presence of substantial (10%) noise. We showcase its abilities on a panel of examples from astrophysics. △ Less

Submitted 9 October, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 29 pages, 9 figures, 11 tables. Accepted for publication at ApJ

Journal ref: ApJ 959 99 (2023)

The PAndAS View of the Andromeda Satellite System. IV Global properties

Authors: Amandine Doliva-Dolinsky, Nicolas F. Martin, Zhen Yuan, Alessandro Savino, Daniel R. Weisz, Annette M. N. Ferguson, Rodrigo A. Ibata, Stacy Y. Kim, Geraint F. Lewis, Alan W. McConnachie, Guillaume F. Thomas

Abstract: We build a statistical framework to infer the global properties of the satellite system of the Andromeda galaxy (M31) from the properties of individual dwarf galaxies located in the Pan-Andromeda Archaelogical Survey (PAndAS) and the previously determined completeness of the survey. Using forward modeling, we infer the slope of the luminosity function of the satellite system, the slope of its spat… ▽ More We build a statistical framework to infer the global properties of the satellite system of the Andromeda galaxy (M31) from the properties of individual dwarf galaxies located in the Pan-Andromeda Archaelogical Survey (PAndAS) and the previously determined completeness of the survey. Using forward modeling, we infer the slope of the luminosity function of the satellite system, the slope of its spatial density distribution, and the size-luminosity relation followed by the dwarf galaxies. We find that the slope of the luminosity function is $βべーた=-1.5\pm0.1$. Combined with the spatial density profile, it implies that, when accounting for survey incompleteness, M31 hosts $92_{-26}^{+19}$ dwarf galaxies with $M_\textrm{V}<-5.5$ and a sky-projected distance from M31 between 30 and 300kpc. We conclude that many faint or distant dwarf galaxies remain to be discovered around Andromeda, especially outside the PAndAS footprint. Finally, we use our model to test if the higher number of satellites situated in the hemisphere facing the Milky Way could be explained simply by the detection limits of dwarf galaxy searches. We rule this out at $>99.9\%$ confidence and conclude that this anisotropy is an intrinsic feature of the M31 satellite system. The statistical framework we present here is a powerful tool to robustly constrain the properties of a satellite system and compare those across hosts, especially considering the upcoming start of the Euclid or Rubin large photometric surveys that are expected to uncover a large number of dwarf galaxies in the Local Volume. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: Submitted to ApJ - 12 pages, 6 figures, 2 tables

Local Group Dwarf Galaxy Detection Limit in the CSST survey

Authors: Han Qu, Zhen Yuan, Amandine Doliva-Dolinsky, Nicolas F. Martin, Xi Kang, Chengliang Wei, Guoliang Li, Yu Luo, Jiang Chang, Chao-Wei Tsai, Zhou Fan, Rodrigo Ibata

Abstract: We predict the dwarf galaxy detection limits for the upcoming Chinese Space Station Telescope (CSST) survey that will cover 17,500 deg$^{2}$ of the sky with a wide field of view of 1.1 deg$^2$. The point-source depth reaches 26.3 mag in the $g$ band and 25.9 mag in the $i$ band. Constructing mock survey data based on the designed photometric bands, we estimate the recovery rate of artificial dwarf… ▽ More We predict the dwarf galaxy detection limits for the upcoming Chinese Space Station Telescope (CSST) survey that will cover 17,500 deg$^{2}$ of the sky with a wide field of view of 1.1 deg$^2$. The point-source depth reaches 26.3 mag in the $g$ band and 25.9 mag in the $i$ band. Constructing mock survey data based on the designed photometric bands, we estimate the recovery rate of artificial dwarf galaxies from mock point-source photometric catalogues. The detection of these artificial dwarf galaxies is strongly dependent on their distance, magnitude and size, in agreement with searches in current surveys. We expect CSST to enable the detection of dwarf galaxies with $M_V = -3.0$ and $μみゅー_{250} = 32.0$ mag/arcsec$^2$ (surface-brightness limit for a system of half-light radius $r_{\rm h}$ = 250 pc at 400 kpc, and $M_V = -4.9$ and $μみゅー_{250} = 30.5$ mag/arcsec$^2$ around the Andromeda galaxy. Beyond the Local Group, the CSST survey will achieve $M_V = -5.8$, and $μみゅー_{250}$ = 29.7 mag/arcsec$^2$ in the distance range of 1--2 Mpc, opening up an exciting discovery space for faint field dwarf galaxies. With its optical bands, wide survey footprint, and space resolution, CSST will undoubtedly expand our knowledge of low-mass dwarf galaxies to an unprecedented volume. △ Less

Submitted 27 April, 2023; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: 11 pages, 10 figures. Accepted to MNRAS

The metallicity distribution in the core of the Sagittarus dwarf spheroidal: minimising the metallicity biases

Authors: Alice Minelli, Michele Bellazzini, Alessio Mucciarelli, Piercarlo Bonifacio, Rodrigo Ibata, Donatella Romano, Lorenzo Monaco, Elisabetta Caffau, Emanuele Dalessandro, Raffaele Pascale

Abstract: We present metallicity and radial velocity for 450 bonafide members of the Sagittarius dwarf spheroidal (Sgr dSph) galaxy, measured from high resolution (R~18000) FLAMES@VLT spectra. The targets were carefully selected (a) to sample the core of the main body of Sgr dSph while avoiding contamination from the central stellar nucleus, and (b) to prevent any bias on the metallicity distribution, by se… ▽ More We present metallicity and radial velocity for 450 bonafide members of the Sagittarius dwarf spheroidal (Sgr dSph) galaxy, measured from high resolution (R~18000) FLAMES@VLT spectra. The targets were carefully selected (a) to sample the core of the main body of Sgr dSph while avoiding contamination from the central stellar nucleus, and (b) to prevent any bias on the metallicity distribution, by selecting targets based on their Gaia parallax and proper motions. All the targets selected in this way were confirmed as radial velocity members. We used this sample to derive the first metallicity distribution of the core of the Sgr dSph virtually unaffected by metallicity biases. The observed distribution ranges from [Fe/H]~ -2.3 to [Fe/H]~ 0.0, with a strong, symmetric and relatively narrow peak around [Fe/H]~ -0.5 and a weak, extended metal-poor tail, with only 13.8 +/- 1.9% of the stars having [Fe/H]< -1.0. We confirm previous evidence of correlations between chemical and kinematical properties of stars in the core of Sgr. In our sample stars with [Fe/H]>= -0.6 display a lower velocity dispersion and a higher rotation amplitude than those with [Fe/H]< -0.6, confirming previous suggestions of a disk/halo structure for the progenitor of the system. △ Less

Submitted 12 November, 2022; originally announced November 2022.

Comments: 11 pages, 10 figures, 2 tables, accepted by A&A

Journal ref: A&A 669, A54 (2023)

Dark matter halo cores and the tidal survival of Milky Way satellites

Authors: Raphaël Errani, Julio F. Navarro, Jorge Peñarrubia, Benoit Famaey, Rodrigo Ibata

Abstract: The cuspy central density profiles of cold dark matter (CDM) haloes make them highly resilient to disruption by tides. Self-interactions between dark matter particles, or the cycling of baryons, may result in the formation of a constant-density core which would make haloes more susceptible to tidal disruption. We use N-body simulations to study the evolution of NFW-like "cored" subhaloes in the ti… ▽ More The cuspy central density profiles of cold dark matter (CDM) haloes make them highly resilient to disruption by tides. Self-interactions between dark matter particles, or the cycling of baryons, may result in the formation of a constant-density core which would make haloes more susceptible to tidal disruption. We use N-body simulations to study the evolution of NFW-like "cored" subhaloes in the tidal field of a massive host, and identify the criteria and timescales for full tidal disruption. Our results imply that the survival of Milky Way satellites places constraints on the sizes of dark matter cores. Indeed, we find that no subhaloes with cores larger than 1 per cent of their initial NFW scale radius can survive for a Hubble time on orbits with pericentres <10 kpc. A satellite like Tucana 3, with pericentre ~3.5 kpc, must have a core size smaller than ~2 pc to survive just three orbital periods on its current orbit. The core sizes expected in self-interacting dark matter (SIDM) models with a velocity-independent cross section of 1 cm^2/g seem incompatible with ultra-faint satellites with small pericentric radii, such as Tuc 3, Seg 1, Seg 2, Ret 2, Tri 2, and Wil 1, as these should have fully disrupted if accreted on to the Milky Way >10 Gyr ago. These results suggest that many satellites have vanishingly small core sizes, consistent with CDM cusps. The discovery of further Milky Way satellites on orbits with small pericentric radii would strengthen these conclusions and allow for stricter upper limits on the core sizes. △ Less

Submitted 26 December, 2022; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 13 pages, 13 figures, minor edits to match accepted version

Journal ref: MNRAS 519, 384 (2023)

Exploring the effects of primordial non-Gaussianity at galactic scales

Authors: Clément Stahl, Thomas Montandon, Benoit Famaey, Oliver Hahn, Rodrigo Ibata

Abstract: While large scale primordial non-Gaussianity is strongly constrained by present-day data, there are no such constraints at Mpc scales. Here we investigate the effect of significant small-scale primordial non-Gaussianity on structure formation and the galaxy formation process with collisionless simulations: specifically, we explore four different types of non-Gaussianities. All of these prescriptio… ▽ More While large scale primordial non-Gaussianity is strongly constrained by present-day data, there are no such constraints at Mpc scales. Here we investigate the effect of significant small-scale primordial non-Gaussianity on structure formation and the galaxy formation process with collisionless simulations: specifically, we explore four different types of non-Gaussianities. All of these prescriptions lead to a distinct and potentially detectable feature in the matter power spectrum around the non-linear scale. The feature might have interesting consequences for the $S_8$ tension. We then show in particular that a negatively-skewed distribution of the potential random field, hence positively skewed in terms of overdensities, with $f_{\rm NL}$ of the order of 1000 at these scales, implies that typical galaxy-sized halos reach half of their present-day mass at an earlier stage and have a quieter merging history at $z<3$ than in the Gaussian case. Their environment between 0.5 and 4 virial radii at $z=0$ is less dense than in the Gaussian case. This quieter history and less dense environment has potentially interesting consequences in terms of the formation of bulges and bars. Moreover, we show that the two most massive subhalos around their host tend to display an interesting anti-correlation of velocities, indicative of kinematic coherence. All these hints will need to be statistically confirmed in larger-box simulations with scale-dependent non-Gaussian initial conditions, followed by hydrodynamical zoom-in simulations to explore the detailed consequences of small-scale non-Gaussianities on galaxy formation. △ Less

Submitted 21 January, 2023; v1 submitted 29 September, 2022; originally announced September 2022.

Comments: 25 pages, 9 figures, 3 tables, JCAP accepted

Andromeda XXV -- a dwarf galaxy with a low central dark matter density

Authors: Emily J. E. Charles, Michelle L. M. Collins, R. Michael Rich, Justin I. Read, Stacy Y. Kim, Rodrigo A. Ibata, Nicolas F. Martin, Scott C. Chapman, Eduardo Balbinot, Daniel R. Weisz

Abstract: Andromeda (And) XXV has previously been reported as a dwarf spheroidal galaxy (dSph) with little-to-no dark matter. However, the uncertainties on this result were significant. In this study, we double the number of member stars and re-derive the kinematics and mass of And XXV. We find that And XXV has a systemic velocity of $νにゅー_\mathrm{r}=-107.7\pm1.0 \mathrm{~km s}^{-1}$ and a velocity dispersion… ▽ More Andromeda (And) XXV has previously been reported as a dwarf spheroidal galaxy (dSph) with little-to-no dark matter. However, the uncertainties on this result were significant. In this study, we double the number of member stars and re-derive the kinematics and mass of And XXV. We find that And XXV has a systemic velocity of $νにゅー_\mathrm{r}=-107.7\pm1.0 \mathrm{~km s}^{-1}$ and a velocity dispersion of $σしぐま_νにゅー=4.5\pm1.0\mathrm{~km s}^{-1}$. With this better constrained velocity dispersion, we derive a mass contained within the half-light radius of $M(r< r_\mathrm{h})=6.9^{+3.2}_{-2.8}\times10^6\mathrm{~M}_\odot$. This mass corresponds to a mass-to-light ratio of $\mathrm{[M/L]}_\mathrm{r_\mathrm{h}}=37^{+17}_{-15}\mathrm{~M}_\odot/\mathrm{L}_\odot$, demonstrating, for the first time, that And XXV has an unambiguous dark matter component. We also measure the metallicity of And XXV to be $\mathrm{[Fe/H]}=-1.9\pm0.1$$\mathrm{~}$dex, which is in agreement with previous results. Finally, we extend the analysis of And XXV to include mass modelling using GravSphere. We find that And XXV has a low central dark matter density, $ρろー_\mathrm{DM}(150\mathrm{pc})= 2.7^{+1.8}_{-1.6}\times10^7\mathrm{~M}_\odot\mathrm{kpc}^{-3}$, making And XXV a clear outlier when compared to other Local Group (LG) dSphs of the similar stellar mass. In a companion paper, we will explore whether some combination of dark matter cusp-core transformations and/or tides can explain And XXV's low density. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 13 pages, 8 figures (7 main, 1 appendix). Submitted to MNRAS

Dipolar dark matter simulations on galaxy scales with the RAMSES code

Authors: Clément Stahl, Benoit Famaey, Guillaume Thomas, Yohan Dubois, Rodrigo Ibata

Abstract: We numerically explore on galaxy scales the Dipolar dark matter (DM) model based on the concept of gravitational polarization. This DM model has been proposed as a natural way to reproduce observed tight galactic scaling relations such as the baryonic Tully-Fisher relation and the Radial Acceleration Relation. We present a customized version of the \texttt{RAMSES} code including for the first time… ▽ More We numerically explore on galaxy scales the Dipolar dark matter (DM) model based on the concept of gravitational polarization. This DM model has been proposed as a natural way to reproduce observed tight galactic scaling relations such as the baryonic Tully-Fisher relation and the Radial Acceleration Relation. We present a customized version of the \texttt{RAMSES} code including for the first time the dynamics of this Dipolar DM in $N$-body simulations. As a first application of this code, we check that we recover an equilibrium configuration that had been found analytically, where a low density Dipolar DM halo is at rest with respect to its central galaxy, recovering the aforementioned scaling relations. A characteristic signature of this equilibrium model is that it harbours a dynamical instability with a characteristic time depending on the Dipolar DM halo density, which we recover numerically. This represents a first step towards more involved simulations needed to test this framework, ranging from galaxy interactions to structure formation. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 9 pages, 5 figures, 1 table, MNRAS accepted

Typhon: a polar stream from the outer halo raining through the Solar neighborhood

Authors: Wassim Tenachi, Pierre-Antoine Oria, Rodrigo Ibata, Benoit Famaey, Zhen Yuan, Anke Arentsen, Nicolas Martin, Akshara Viswanathan

Abstract: We report on the discovery in the Gaia DR3 astrometric and spectroscopic catalog of a new polar stream that is found as an over-density in action space. This structure is unique as it has an extremely large apocenter distance, reaching beyond 100 kpc, and yet is detected as a coherent moving structure in the Solar neighborhood with a width of $\sim 4$ kpc. A sub-sample of these stars that was fort… ▽ More We report on the discovery in the Gaia DR3 astrometric and spectroscopic catalog of a new polar stream that is found as an over-density in action space. This structure is unique as it has an extremely large apocenter distance, reaching beyond 100 kpc, and yet is detected as a coherent moving structure in the Solar neighborhood with a width of $\sim 4$ kpc. A sub-sample of these stars that was fortuitously observed by LAMOST has a mean spectroscopic metallicity of $\langle {\rm [Fe/H]}\rangle = -1.60^{+0.15}_{-0.16}$ dex and possesses a resolved metallicity dispersion of $σしぐま({\rm [Fe/H]}) = 0.32^{+0.17}_{-0.06}$ dex. The physical width of the stream, the metallicity dispersion and the vertical action spread indicate that the progenitor was a dwarf galaxy. The existence of such a coherent and highly radial structure at their pericenters in the vicinity of the Sun suggests that many other dwarf galaxy fragments may be lurking in the outer halo. △ Less

Submitted 29 February, 2024; v1 submitted 21 June, 2022; originally announced June 2022.

Comments: 7 pages, 4 figures, 1 table. Accepted for publication at ApJ. (Table typo corrected)

Journal ref: ApJL 935 L22 (2022)

Antaeus: a retrograde group of tidal debris in the Milky Way's disk plane

Authors: Pierre-Antoine Oria, Wassim Tenachi, Rodrigo Ibata, Benoit Famaey, Zhen Yuan, Anke Arentsen, Nicolas Martin, Akshara Viswanathan

Abstract: We present the discovery of a wide retrograde moving group in the disk plane of the Milky Way using action-angle coordinates derived from the \textit{Gaia} DR3 catalog. The structure is identified from a sample of its members that are currently almost at the pericenter of their orbit and are passing through the Solar neighborhood. The motions of the stars in this group are highly correlated, indic… ▽ More We present the discovery of a wide retrograde moving group in the disk plane of the Milky Way using action-angle coordinates derived from the \textit{Gaia} DR3 catalog. The structure is identified from a sample of its members that are currently almost at the pericenter of their orbit and are passing through the Solar neighborhood. The motions of the stars in this group are highly correlated, indicating that the system is probably not phase mixed. With a width of at least 1.5 kpc and with a probable intrinsic spread in metallicity, this structure is most likely the wide remnant of a tidal stream of a disrupted ancient dwarf galaxy (age $\sim 12$ Gyr, $\langle {\rm [Fe/H]} \rangle \sim -1.74$). The structure presents many similarities (e.g. in energy, angular momentum, metallicity, and eccentricity) with the Sequoia merging event. However, it possesses extremely low vertical action $J_z$ which makes it unique even amongst Sequoia dynamical groups. As the low $J_z$ may be attributable to dynamical friction, we speculate that the these stars may be the remnants of the dense core of the Sequoia progenitor. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 6 pages, 4 figures, 1 table. Submitted to ApJL

Journal ref: ApJL 936 L3 (2022)

The Hubble Space Telescope Survey of M31 Satellite Galaxies I. RR Lyrae-based Distances and Refined 3D Geometric Structure

Authors: Alessandro Savino, Daniel R. Weisz, Evan D. Skillman, Andrew Dolphin, Nitya Kallivayalil, Andrew Wetzel, Jay Anderson, Gurtina Besla, Michael Boylan-Kolchin, James S. Bullock, Andrew A. Cole, Michelle L. M. Collins, M. C. Cooper, Alis J. Deason, Aaron L. Dotter, Mark Fardal, Annette M. N. Ferguson, Tobias K. Fritz, Marla C. Geha, Karoline M. Gilbert, Puragra Guhathakurta, Rodrigo Ibata, Michael J. Irwin, Myoungwon Jeon, Evan Kirby , et al. (11 additional authors not shown)

Abstract: We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days… ▽ More We measure homogeneous distances to M31 and 38 associated stellar systems ($-$16.8$\le M_V \le$ $-$6.0), using time-series observations of RR Lyrae stars taken as part of the Hubble Space Telescope Treasury Survey of M31 Satellites. From $>700$ orbits of new/archival ACS imaging, we identify $>4700$ RR Lyrae stars and determine their periods and mean magnitudes to a typical precision of 0.01 days and 0.04 mag. Based on Period-Wesenheit-Metallicity relationships consistent with the Gaia eDR3 distance scale, we uniformly measure heliocentric and M31-centric distances to a typical precision of $\sim20$ kpc (3%) and $\sim10$ kpc (8%), respectively. We revise the 3D structure of the M31 galactic ecosystem and: (i) confirm a highly anisotropic spatial distribution such that $\sim80$% of M31's satellites reside on the near side of M31; this feature is not easily explained by observational effects; (ii) affirm the thin (rms $7-23$ kpc) planar "arc" of satellites that comprises roughly half (15) of the galaxies within 300 kpc from M31; (iii) reassess physical proximity of notable associations such as the NGC 147/185 pair and M33/AND XXII; and (iv) illustrate challenges in tip-of-the-red-giant branch distances for galaxies with $M_V > -9.5$, which can be biased by up to 35%. We emphasize the importance of RR Lyrae for accurate distances to faint galaxies that should be discovered by upcoming facilities (e.g., Rubin Observatory). We provide updated luminosities and sizes for our sample. Our distances will serve as the basis for future investigation of the star formation and orbital histories of the entire known M31 satellite system. △ Less

Submitted 12 September, 2022; v1 submitted 6 June, 2022; originally announced June 2022.

Comments: 32 pages, 15 figures, accepted for publication on ApJ

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

Authors: Pierre-Antoine Oria, Rodrigo Ibata, Pau Ramos, Benoit Famaey, Raphaël Errani

Abstract: 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 p… ▽ More 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. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 8 pages, 5 figures. Accepted for publication in ApJL

The PAndAS View of the Andromeda Satellite System. III. Dwarf galaxy detection limits

Authors: Amandine Doliva-Dolinsky, Nicolas F. Martin, Guillaume F. Thomas, Annette M. N. Ferguson, Rodrigo A. Ibata, Geraint F. Lewis, Dougal Mackey, Alan W. McConnachie, Zhen Yuan

Abstract: We determine the detection limits of the search for dwarf galaxies in the Pan-Andromeda Archaeological Survey (PAndAS) using the algorithm developed by the PAndAS team. The recovery fractions of artificial dwarf galaxies are, as expected, a strong function of physical size and luminosity and, to a lesser extent, distance. We show that these recovery fractions vary strongly with location in the sur… ▽ More We determine the detection limits of the search for dwarf galaxies in the Pan-Andromeda Archaeological Survey (PAndAS) using the algorithm developed by the PAndAS team. The recovery fractions of artificial dwarf galaxies are, as expected, a strong function of physical size and luminosity and, to a lesser extent, distance. We show that these recovery fractions vary strongly with location in the surveyed area because of varying levels of contamination from both the Milky Way foreground stars and the stellar halo of Andromeda. We therefore provide recovery fractions that are a function of size, luminosity, and location within the survey on a scale of 1 square degree. Overall, the effective surface brightness for a 50-percent detection rate range between 28 and 30 mag per square arcsecond. This is in line with expectations for a search that relies on photometric data that are as deep as the PAndAS survey. The derived detection limits are an essential ingredient on the path to constraining the global properties of Andromeda's system of satellite dwarf galaxies and, more broadly, to provide constraints on dwarf galaxy formation and evolution in a cosmological context. △ Less

Submitted 5 May, 2022; originally announced May 2022.

Comments: 17 pages, 10 figures, accepted in ApJ

The Pristine Inner Galaxy Survey (PIGS) IV: A photometric metallicity analysis of the Sagittarius dwarf spheroidal galaxy

Authors: Sara Vitali, Anke Arentsen, Else Starkenburg, Paula Jofré, Nicolas F. Martin, David S. Aguado, Raymond Carlberg, Jonay I. González Hernández, Rodrigo Ibata, Georges Kordopatis, Khyati Malhan, Pau Ramos, Federico Sestito, Zhen Yuan, Sven Buder, Geraint F. Lewis, Zhen Wan, Daniel B. Zucker

Abstract: We present a comprehensive metallicity analysis of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using $Pristine\,CaHK$ photometry. We base our member selection on $Gaia$ EDR3 astrometry applying a magnitude limit at $G_{0} = 17.3$, and our population study on the metallicity-sensitive photometry from the $Pristine$ Inner Galaxy Survey (PIGS). Working with photometric metallicities instead of… ▽ More We present a comprehensive metallicity analysis of the Sagittarius dwarf spheroidal galaxy (Sgr dSph) using $Pristine\,CaHK$ photometry. We base our member selection on $Gaia$ EDR3 astrometry applying a magnitude limit at $G_{0} = 17.3$, and our population study on the metallicity-sensitive photometry from the $Pristine$ Inner Galaxy Survey (PIGS). Working with photometric metallicities instead of spectroscopic metallicities allows us to cover an unprecedented large area ($\sim 100$ square degrees) of the dwarf galaxy, and to study the spatial distribution of its members as function of metallicity with little selection effects. Our study compares the spatial distributions of a metal-poor population of 9719 stars with [Fe/H] $< -1.3$ and a metal rich one of 30115 stars with [Fe/H] $> -1.0$. The photometric Sgr sample also allows us to assemble the largest sample of 1150 very metal-poor Sgr candidates ([Fe/H] $< -2.0$). By investigating and fitting the spatial properties of the metal-rich and metal-poor population, we find a negative metallicity gradient which extends up to 12 degrees from the Sgr center (or $\sim 5.5$ kpc at the distance of Sgr), the limit of our footprint. We conclude that the relative number of metal-poor stars increases in the outer areas of the galaxy, while the central region is dominated by metal-rich stars. These finding suggest an outside-in formation process and are an indication of the extended formation history of Sgr, which has been affected by the tidal interaction between Sgr and the Milky Way. △ Less

Submitted 4 October, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: The paper was accepted on 20 September by MNRAS

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

Authors: Indranil Banik, Ingo Thies, Roy Truelove, Graeme Candlish, Benoit Famaey, Marcel S. Pawlowski, Rodrigo Ibata, Pavel Kroupa

Abstract: 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 RAMS… ▽ More 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$. △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: 33 pages, 18 figures, 10 tables. Published in the Monthly Notices of the Royal Astronomical Society in this form

Journal ref: Monthly Notices of the Royal Astronomical Society, volume 513, issue 1, pages 129 - 158 (2022)

The Pristine survey XVIII: C-19: Tidal debris of a dark matter-dominated globular cluster?

Authors: Raphaël Errani, Julio F. Navarro, Rodrigo Ibata, Nicolas Martin, Zhen Yuan, David S. Aguado, Piercarlo Bonifacio, Elisabetta Caffau, Jonay I. González Hernández, Khyati Malhan, Rubén Sánchez-Janssen, Federico Sestito, Else Starkenburg, Guillaume F. Thomas, Kim A. Venn

Abstract: The recently discovered C-19 stellar stream is a collection of kinematically associated metal-poor stars in the halo of the Milky Way lacking an obvious progenitor. The stream spans an arc of ~15 degrees in the sky, and orbit-fitting suggests an apocentric distance of ~20 kpc and a pericentre of ~10 kpc. The narrow metallicity dispersion of stars with available spectra, together with light element… ▽ More The recently discovered C-19 stellar stream is a collection of kinematically associated metal-poor stars in the halo of the Milky Way lacking an obvious progenitor. The stream spans an arc of ~15 degrees in the sky, and orbit-fitting suggests an apocentric distance of ~20 kpc and a pericentre of ~10 kpc. The narrow metallicity dispersion of stars with available spectra, together with light element abundance variations, suggests a globular cluster (GC) origin. The observed metallicity ([Fe/H] ~ -3.4), however, is much lower than that of any known GC. In addition, the width and velocity dispersion of the stream are similar to those expected from disrupting dwarf galaxies, and substantially larger than the tidal debris of GCs able to disrupt on C-19's orbit. We propose here an unconventional model where the C-19 progenitor is a dark matter-dominated stellar system with GC-like abundance patterns. We use N-body simulations to show that the tidal disruption of a ~100 pc King-model stellar component embedded in a ~20 km/s cuspy cold dark matter halo yields debris consistent with C-19's observed width and velocity dispersion. The stellar component of the progenitor is fully disrupted, and is spread over two distinct streams; one corresponding to C-19 and another possibly hiding behind the Galactic plane. If such companion stream were found, it would suggest that dark matter-dominated dwarfs may also develop GC-like enrichment patterns, a finding that would inform our theoretical understanding of the formation of multiple populations in GCs and dwarf galaxies alike. △ Less

Submitted 21 June, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 9 pages, 9 figures. minor edits to match accepted version

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

Authors: Zhen Yuan, Nicolas F. Martin, Rodrigo A. Ibata, Elisabetta Caffau, Piercarlo Bonifacio, Lyudmila I. Mashonkina, Raphaël Errani, Amandine Doliva-Dolinsky, Else Starkenburg, Kim A. Venn, Anke Arentsen, David S. Aguado, Michele Bellazzini, Benoit Famaey, Morgan Fouesneau, Jonay I. González Hernández, Pascale Jablonka, Carmela Lardo, Khyati Malhan, Julio F. Navarro, Rubén Sánchez Janssen, Federico Sestito, Guillaume F. Thomas, Akshara Viswanathan, Sara Vitali

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 16 May, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 7 pages, 6 figures, accepted to MNRAS. The C-19 member list is published via https://github.com/zyuan-astro/C-19

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

Authors: Khyati Malhan, Rodrigo A. Ibata, Sanjib Sharma, Benoit Famaey, Michele Bellazzini, Raymond G. Carlberg, Richard D'Souza, Zhen Yuan, Nicolas F. Martin, Guillaume F. Thomas

Abstract: 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, stella… ▽ More 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. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: Published in ApJ. Key results in Table 3. The paper contains 38 pages, 17 figures, 4 Tables

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

Authors: Simon Rozier, Benoit Famaey, Arnaud Siebert, Giacomo Monari, Christophe Pichon, Rodrigo Ibata

Abstract: 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 (asso… ▽ More 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. △ Less

Submitted 14 January, 2022; originally announced January 2022.

Comments: 24 pages, 13 figures. Submitted to ApJ

Rubin-Euclid Derived Data Products: Initial Recommendations

Authors: Leanne P. Guy, Jean-Charles Cuillandre, Etienne Bachelet, Manda Banerji, Franz E. Bauer, Thomas Collett, Christopher J. Conselice, Siegfried Eggl, Annette Ferguson, Adriano Fontana, Catherine Heymans, Isobel M. Hook, Éric Aubourg, Hervé Aussel, James Bosch, Benoit Carry, Henk Hoekstra, Konrad Kuijken, Francois Lanusse, Peter Melchior, Joseph Mohr, Michele Moresco, Reiko Nakajima, Stéphane Paltani, Michael Troxel , et al. (95 additional authors not shown)

Abstract: This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum… ▽ More This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology. △ Less

Submitted 13 October, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Report of the Rubin-Euclid Derived Data Products Working Group, 78 pages, 11 figures

New constraints on the dark matter density profiles of dwarf galaxies from proper motions of globular cluster streams

Authors: Khyati Malhan, Monica Valluri, Katherine Freese, Rodrigo A. Ibata

Abstract: The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter. Recently, in Malhan et al. (2021), we employed N-body simulations to show that the cuspy cold dark matter (CDM) subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripp… ▽ More The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter. Recently, in Malhan et al. (2021), we employed N-body simulations to show that the cuspy cold dark matter (CDM) subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripping within their parent halos prior to accretion onto the Milky Way. We previously found that clusters that are accreted within cuspy subhalos produce streams with larger physical widths and higher dispersions in line-of-sight velocity and angular momentum than streams that are accreted within cored subhalos. Here, we use the same suite of simulations to demonstrate that the dispersion in the tangential velocities of streams ($σしぐま_{v_\mathrm{Tan}}$) is also sensitive to the central DM density profiles of their parent dwarfs and GCs that were accreted from: cuspy subhalos produce streams with larger $σしぐま_{v_\mathrm{Tan}}$ than those accreted inside cored subhalos. Using Gaia EDR3 observations of multiple GC streams we compare their $σしぐま_{v_\mathrm{Tan}}$ values with simulations. The measured $σしぐま_{v_\mathrm{Tan}}$ values are consistent with both an ``in situ'' origin and with accretion inside cored subhalos of $M\sim 10^{8-9}M_{\odot}$ (or very low-mass cuspy subhalos of mass $\sim 10^8M_{\odot}$). Despite the large current uncertainties in $σしぐま_{v_\mathrm{Tan}}$, we find a low probability that any of the progenitor GCs were accreted from cuspy subhalos of $M_{\rm subhalo}\buildrel > \over \sim$ $10^9 M_{\odot}$. The uncertainties on Gaia tangential velocity measurements are expected to decrease in future and will allow for stronger constraints on subhalo DM density profiles. △ Less

Submitted 29 November, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: 14 pages, 5 figures, accepted for publication in ApJ Letters

The Pristine survey -- XVI. The metallicity of 26 stellar streams around the Milky Way detected with the STREAMFINDER in Gaia EDR3

Authors: Nicolas F. Martin, Rodrigo A. Ibata, Else Starkenburg, Zhen Yuan, Khyati Malhan, Michele Bellazzini, Akshara Viswanathan, David Aguado, Anke Arentsen, Piercarlo Bonifacio, Ray Carlberg, Jonay I. González Hernández, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Carmela Lardo, Alan W. McConnachie, Julio Navarro, Rubén Sánchez-Janssen, Federico Sestito, Guillaume F. Thomas, Kim A. Venn, Sara Vitali, Karina T. Voggel

Abstract: We use the photometric metallicities provided by the panoramic Pristine survey to study the veracity and derive the metallicities of the numerous stellar streams found by the application of the STREAMFINDER algorithm to the Gaia EDR3 data. All 26 streams present in Pristine show a clear metallicity distribution function, which provides an independent check of the reality of these structures, suppo… ▽ More We use the photometric metallicities provided by the panoramic Pristine survey to study the veracity and derive the metallicities of the numerous stellar streams found by the application of the STREAMFINDER algorithm to the Gaia EDR3 data. All 26 streams present in Pristine show a clear metallicity distribution function, which provides an independent check of the reality of these structures, supporting the reliability of STREAMFINDER in finding streams and the power of Pristine to measure precise metallicities. We further present 6 candidate structures with coherent phase-space and metallicity signals that are very likely streams. The majority of studied streams are very metal-poor (14 structures with [Fe/H]<-2.0) and include 3 systems with [Fe/H]<-2.9 (C-11, C-19, and C-20). These streams could be the closest debris of low-luminosity dwarf galaxies or may have originated from globular clusters of significantly lower metallicity than any known current Milky Way globular cluster. Our study shows that the promise of the Gaia data for Galactic Archeology studies can be substantially strengthened by quality photometric metallicities, allowing us to peer back into the earliest epochs of the formation of our Galaxy and its stellar halo constituents. △ Less

Submitted 25 August, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 25 pages, 29 figures, accepted for publication in MNRAS. Updated version has only minor changes compared to initial submission. The metallicity of most streams have marginally changed in the updated version

A stellar stream remnant of a globular cluster below the metallicity floor

Authors: Nicolas F. Martin, Kim A. Venn, David S. Aguado, Else Starkenburg, Jonay I. González Hernández, Rodrigo A. Ibata, Piercarlo Bonifacio, Elisabetta Caffau, Federico Sestito, Anke Arentsen, Carlos Allende Prieto, Raymond G. Carlberg, Sébastien Fabbro, Morgan Fouesneau, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Carmela Lardo, Khyati Malhan, Lyudmila I. Mashonkina, Alan W. McConnachie, Julio F. Navarro, Rubén Sánchez Janssen, Guillaume F. Thomas, Zhen Yuan , et al. (1 additional authors not shown)

Abstract: Stellar ejecta gradually enrich the gas out of which subsequent stars form, making the least chemically enriched stellar systems direct fossils of structures formed in the early universe. Although a few hundred stars with metal content below one thousandth of the solar iron content are known in the Galaxy, none of them inhabit globular clusters, some of the oldest known stellar structures. These s… ▽ More Stellar ejecta gradually enrich the gas out of which subsequent stars form, making the least chemically enriched stellar systems direct fossils of structures formed in the early universe. Although a few hundred stars with metal content below one thousandth of the solar iron content are known in the Galaxy, none of them inhabit globular clusters, some of the oldest known stellar structures. These show metal content of at least ~0.2 percent of the solar metallicity ([Fe/H] > -2.7). This metallicity floor appears universal and it has been proposed that proto-galaxies that merge into the galaxies we observe today were simply not massive enough to form clusters that survived to the present day. Here, we report the discovery of a stellar stream, C-19, whose metallicity is less than 0.05 per cent the solar metallicity ([Fe/H]=-3.38 +/- 0.06 (stat.) +/- 0.20 (syst.)). The low metallicity dispersion and the chemical abundances of the C-19 stars show that this stream is the tidal remnant of the most metal-poor globular cluster ever discovered, and significantly below the purported metallicity floor: clusters with significantly lower metallicities than observed today existed in the past and contributed their stars to the Milky Way halo. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: 14 pages, 6 figures, 5 tables. Nature, accepted

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

Authors: Zhen Yuan, Khyati Malhan, Federico Sestito, Rodrigo A. Ibata, Nicolas F. Martin, Jiang Chang, Ting S. Li, Elisabetta Caffau, Piercarlo Bonifacio, Michele Bellazzini, Yang Huang, Karina Voggel, Nicolas Longeard, Anke Arentsen, Amandine Doliva-Dolinsky, Julio Navarro, Benoit Famaey, Else Starkenburg, David S. Aguado

Abstract: 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 differen… ▽ More 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. △ Less

Submitted 8 April, 2022; v1 submitted 10 December, 2021; originally announced December 2021.

Comments: accepted to ApJ. The N-body model of the Cetus system is published at https://zenodo.org/record/5771585

The Sagittarius stream in Gaia eDR3 and the origin of the bifurcations

Authors: P. Ramos, T. Antoja, Z. Yuan, A. Arentsen, P. -A. Oria, B. Famaey, R. Ibata, C. Mateu, J. A. Carballo-Bello

Abstract: The Sagittarius dwarf spheroidal (Sgr) is a dissolving galaxy being tidally disrupted by the Milky Way (MW). Its stellar stream still poses serious modelling challenges, which hinders our ability to use it effectively as a probe of the MW gravitational potential. Our goal is to construct the best possible sample of stars with which we can advance our understanding of the Sgr-MW interaction, focusi… ▽ More The Sagittarius dwarf spheroidal (Sgr) is a dissolving galaxy being tidally disrupted by the Milky Way (MW). Its stellar stream still poses serious modelling challenges, which hinders our ability to use it effectively as a probe of the MW gravitational potential. Our goal is to construct the best possible sample of stars with which we can advance our understanding of the Sgr-MW interaction, focusing on the characterisation of the bifurcations. We improve on previous methods based on the use of the wavelet transform to systematically search for the kinematic signature of the Sgr stream in the Gaia data. We refine our selection via the use of a clustering algorithm on the statistical properties of the colour-magnitude diagrams. Our final sample contains > 700k candidate stars, 3x larger than previous samples. We have been able to detect the bifurcation of the stream in both the northern and southern hemispheres, requiring four branches to fully describe this system. We present the detailed proper motion distribution of the trailing arm as a function of Lambda showing the presence of a sharp edge (on the side of the small proper motions) beyond which there are no Sgr stars. We also characterise the correlation between kinematics and distance. The chemical analysis of our sample shows a significant difference between the faint and bright branches. We provide analytical descriptions for the proper motion trends as well as for the sky distribution of the four branches of the stream. We interpret the bifurcations as the misaligned overlap of the material stripped at the antepenultimate pericentre (faint branches) with the stars ejected at the penultimate pericentre (bright branch). The source of this misalignment is still unknown but we argue that models with some internal rotation in the progenitor, at least during the time of stripping of the stars in the faint branches, are worth exploring. △ Less

Submitted 3 December, 2021; originally announced December 2021.

Comments: Submitted to A&A on the 03/12/2021. Online content available upon publication

Journal ref: A&A 666, A64 (2022)

Structure and kinematics of tidally limited satellite galaxies in LCDM

Authors: Raphaël Errani, Julio F. Navarro, Rodrigo Ibata, Jorge Peñarrubia

Abstract: We use N-body simulations to model the tidal evolution of dark matter-dominated dwarf spheroidal galaxies embedded in cuspy Navarro-Frenk-White subhalos. Tides gradually peel off stars and dark matter from a subhalo, trimming it down according to their initial binding energy. This process strips preferentially particles with long orbital times, and comes to an end when the remaining bound particle… ▽ More We use N-body simulations to model the tidal evolution of dark matter-dominated dwarf spheroidal galaxies embedded in cuspy Navarro-Frenk-White subhalos. Tides gradually peel off stars and dark matter from a subhalo, trimming it down according to their initial binding energy. This process strips preferentially particles with long orbital times, and comes to an end when the remaining bound particles have crossing times shorter than a fraction of the orbital time at pericentre. The properties of the final stellar remnant thus depend on the energy distribution of stars in the progenitor subhalo, which in turn depends on the initial density profile and radial segregation of the initial stellar component. The stellar component may actually be completely dispersed if its energy distribution does not extend all the way to the subhalo potential minimum, although a bound dark remnant may remain. These results imply that 'tidally-limited' galaxies, defined as systems whose stellar components have undergone substantial tidal mass loss, neither converge to a unique structure nor follow a single tidal track, as claimed in earlier work. On the other hand, tidally limited dwarfs do have characteristic sizes and velocity dispersions that trace directly the characteristic radius ($r_{max}$) and circular velocity ($V_{max}$) of the subhalo remnant. This result places strong upper limits on the size of satellites whose unusually low velocity dispersions are often ascribed to tidal effects. In particular, the large size of kinematically-cold 'feeble giant' satellites like Crater 2 or Antlia 2 cannot be explained as due to tidal effects alone in the Lambda Cold Dark Matter scenario. △ Less

Submitted 17 March, 2022; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 18 pages, 21 figures. Edited to match accepted version

Journal ref: MNRAS 511, 6001-6018 (2022)

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

Authors: Marcel S. Pawlowski, Pierre-Antoine Oria, Salvatore Taibi, Benoit Famaey, Rodrigo Ibata

Abstract: 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… ▽ More 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. △ Less

Submitted 3 May, 2022; v1 submitted 9 November, 2021; originally announced November 2021.

Comments: 13 pages, 8 figures. Accepted for publication in ApJ