Showing 1–50 of 82 results for author: Laporte, P

On the Galactic rotation curve inferred from the Jeans equations Assessing its robustness using Gaia DR3 and cosmological simulations

Authors: Orlin Koop, Teresa Antoja, Amina Helmi, Thomas M. Callingham, Chervin F. P. Laporte

Abstract: Several works have recently applied Jeans modelling to Gaia-based datasets to infer the circular velocity curve for the Milky Way. Such works have consistently found evidence for a continuous decline in the rotation curve beyond $\sim$15kpc possibly indicative of a light dark matter halo. We used Gaia DR3 RVS data, supplemented with Bayesian distances to determine the radial variation of the secon… ▽ More Several works have recently applied Jeans modelling to Gaia-based datasets to infer the circular velocity curve for the Milky Way. Such works have consistently found evidence for a continuous decline in the rotation curve beyond $\sim$15kpc possibly indicative of a light dark matter halo. We used Gaia DR3 RVS data, supplemented with Bayesian distances to determine the radial variation of the second moments of the velocity distribution for stars close to the Galactic plane. We have used these profiles to determine the rotation curve using the Jeans equations under the assumption of axisymmetry and explored how they vary with azimuth and above and below the Galactic disk plane. We have applied the same methodology to an N-body simulation of a Milky Way-like galaxy impacted by a satellite akin the Sagittarius dwarf and to the Auriga suite of cosmological simulations. We reveal evidence of disequilibrium and deviations from axisymmetry closer in. We find that the second moment of $V_R$ flattens out at $R \gtrsim 12.5$kpc, and that the second moment of $V_φ$ is different above and below the plane for $R \gtrsim 11$kpc. The simulations indicate that these features are typical of galaxies that have been perturbed by external satellites. They also suggest that the difference between the true circular velocity curve and that inferred from Jeans equations can be as high as 15$\%$, but is likely of order 10$\%$ for the Milky Way. This is of larger amplitude than the systematics inherent to Jeans equations. However, if the density of the tracer population were truncated at large radii, the erroneous conclusion of a steeply declining rotation curve can be reached. We find that steady-state axisymmetric Jeans modelling becomes less robust at large radii, indicating that particular caution is needed when interpreting the rotation curve inferred in those regions. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 20 pages, 22 figures

[X/Fe] Marks the Spot: Mapping Chemical Azimuthal Variations in the Galactic Disk with APOGEE

Authors: Zoe Hackshaw, Keith Hawkins, Carrie Filion, Danny Horta, Chervin F. P. Laporte, Chris Carr, Adrian M. Price-Whelan

Abstract: Chemical cartography of the Galactic disk provides insights to its structure and assembly history over cosmic time. In this work, we use chemical cartography to explore chemical gradients and azimuthal substructure in the Milky Way disk with giant stars from APOGEE DR17. We confirm the existence of a radial metallicity gradient in the disk of $Δ$[Fe/H]/$Δ$R $\sim -0.066 \pm 0.0004$ dex/kpc and a v… ▽ More Chemical cartography of the Galactic disk provides insights to its structure and assembly history over cosmic time. In this work, we use chemical cartography to explore chemical gradients and azimuthal substructure in the Milky Way disk with giant stars from APOGEE DR17. We confirm the existence of a radial metallicity gradient in the disk of $Δ$[Fe/H]/$Δ$R $\sim -0.066 \pm 0.0004$ dex/kpc and a vertical metallicity gradient of $Δ$[Fe/H]/$Δ$Z $\sim -0.164 \pm 0.001$ dex/kpc. We find azimuthal variations ($\pm0.1$ dex) on top of the radial metallicity gradient that have been previously established with other surveys. The APOGEE giants show strong correlations with stellar age and the intensity of azimuthal variations in iron; older stellar populations show the largest deviations from the radial metallicity gradient. Beyond iron, we show that other elements (e.g., Mg, O) display azimuthal variations at the $\pm0.05$ dex-level across the Galactic disk. We illustrate that moving into the orbit-space could help constrain the mechanisms producing these azimuthal metallicity variations. These results suggest that the spiral arms of the Galaxy are not solely responsible for azimuthal metallicity variations and other Galactic processes are at play. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 19 pages, 11 figures

Mapping the anisotropic Galactic stellar halo with Blue Horizontal Branch stars

Authors: João A. S. Amarante, Sergey E. Koposov, Chervin F. P. Laporte

Abstract: We use Legacy Survey photometric data to probe the stellar halo in multiple directions of the sky using a probabilistic methodology to identify Blue Horizontal Branch (BHB) stars. The measured average radial density profile follows a double power law in the range $ 5 < r_{gc}/{\rm kpc} < 120$, with a density break at $r_{gc}\approx20$ kpc. This description, however, falls short, depending on the c… ▽ More We use Legacy Survey photometric data to probe the stellar halo in multiple directions of the sky using a probabilistic methodology to identify Blue Horizontal Branch (BHB) stars. The measured average radial density profile follows a double power law in the range $ 5 < r_{gc}/{\rm kpc} < 120$, with a density break at $r_{gc}\approx20$ kpc. This description, however, falls short, depending on the chosen line-of-sight, with some regions showing no signature of a break in the profile and a wide range of density slopes, e.g. outer slope $-5.5 \lesssim α_{out} \lesssim -4$, pointing towards a highly anisotropic stellar halo. This explains in part the wide range of density profiles reported in the literature owing to different tracers and sky coverage. Using our detailed 3-D stellar halo density map, we quantify the shape of the Pisces overdensity associated with the transient wake response of the Galaxy's (dark) halo to the Large Magellanic Cloud (LMC). Measured in the LMC's coordinate system, Pisces stands above the background, is 60 degrees long and 25 degrees wide and aligned with the LMC's orbit. This would correspond to a wake width of $\sim 32$ kpc at $\sim 70$ kpc. We do not find a statistically significant signature of the collective response in density as previously reported in the literature measured with K giant stars, despite our larger numbers. We release the catalogue constructed in this study with 95,446 possible BHB stars and their BHB probability. △ Less

Submitted 25 June, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A. The catalogue can be found here: https://zenodo.org/records/12155488

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

Azimuthal metallicity variations, spiral structure, and the failure of radial actions based on assuming axisymmetry

Authors: Victor P. Debattista, Tigran Khachaturyants, Joao A. S. Amarante, Christopher Carr, Leandro Beraldo e Silva, Chervin F. P. Laporte

Abstract: We study azimuthal variations in the mean stellar metallicity, <[Fe/H]>, in a self-consistent, isolated simulation in which all stars form out of gas. We find <[Fe/H]> variations comparable to those observed in the Milky Way and which are coincident with the spiral density waves. The azimuthal variations are present in young and old stars and therefore are not a result of recently formed stars. Si… ▽ More We study azimuthal variations in the mean stellar metallicity, <[Fe/H]>, in a self-consistent, isolated simulation in which all stars form out of gas. We find <[Fe/H]> variations comparable to those observed in the Milky Way and which are coincident with the spiral density waves. The azimuthal variations are present in young and old stars and therefore are not a result of recently formed stars. Similar variations are present in the mean age and alpha-abundance. We measure the pattern speeds of the <[Fe/H]>-variations and find that they match those of the spirals, indicating that they are at the origin of the metallicity patterns. Because younger stellar populations are not just more [Fe/H]-rich and alpha-poor but also dynamically cooler, we expect them to more strongly support spirals, which is indeed the case in the simulation. However, if we measure the radial action, J_R, using the Stackel axisymmetric approximation, we find that the spiral ridges are traced by regions of high J_R, contrary to expectations. Assuming that the passage of stars through the spirals leads to unphysical variations in the measured J_R, we obtain an improved estimate of J_R by averaging over a 1 Gyr time interval. This time-averaged J_R is a much better tracer of the spiral structure, with minima at the spiral ridges. We conclude that the errors incurred by the axisymmetric approximation introduce correlated deviations large enough to render the instantaneous radial action inadequate for tracing spirals. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 21 pages. Submitted to MNRAS. Comments welcomed

Spectacular nucleosynthesis from early massive stars

Authors: Alexander P. Ji, Sanjana Curtis, Nicholas Storm, Vedant Chandra, Kevin C. Schlaufman, Keivan G. Stassun, Alexander Heger, Marco Pignatari, Adrian M. Price-Whelan, Maria Bergemann, Guy S. Stringfellow, Carla Frohlich, Henrique Reggiani, Erika M. Holmbeck, Jamie Tayar, Shivani P. Shah, Emily J. Griffith, Chervin F. P. Laporte, Andrew R. Casey, Keith Hawkins, Danny Horta, William Cerny, Pierre Thibodeaux, Sam A. Usman, Joao A. S. Amarante , et al. (17 additional authors not shown)

Abstract: Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931… ▽ More Stars formed with initial mass over 50 Msun are very rare today, but they are thought to be more common in the early universe. The fates of those early, metal-poor, massive stars are highly uncertain. Most are expected to directly collapse to black holes, while some may explode as a result of rotationally powered engines or the pair-creation instability. We present the chemical abundances of J0931+0038, a nearby low-mass star identified in early followup of SDSS-V Milky Way Mapper, which preserves the signature of unusual nucleosynthesis from a massive star in the early universe. J0931+0038 has relatively high metallicity ([Fe/H] = -1.76 +/- 0.13) but an extreme odd-even abundance pattern, with some of the lowest known abundance ratios of [N/Fe], [Na/Fe], [K/Fe], [Sc/Fe], and [Ba/Fe]. The implication is that a majority of its metals originated in a single extremely metal-poor nucleosynthetic source. An extensive search through nucleosynthesis predictions finds a clear preference for progenitors with initial mass > 50 Msun, making J0931+0038 one of the first observational constraints on nucleosynthesis in this mass range. However the full abundance pattern is not matched by any models in the literature. J0931+0038 thus presents a challenge for the next generation of nucleosynthesis models and motivates study of high-mass progenitor stars impacted by convection, rotation, jets, and/or binary companions. Though rare, more examples of unusual early nucleosynthesis in metal-poor stars should be found in upcoming large spectroscopic surveys. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 11 pages + 22 page appendix, accepted to ApJL

All-sky Kinematics and Chemistry of Monoceros Stellar Overdensity

Authors: Lais Borbolato, Hélio D. Perottoni, Silvia Rossi, Guilherme Limberg, Angeles Pérez-Villegas, Friedrich Anders, Teresa Antoja, Chervin F. P. Laporte, Helio J. Rocha-Pinto, Rafael M. Santucci

Abstract: We explore the kinematic and chemical properties of Monoceros stellar overdensity by combining data from 2MASS, WISE, APOGEE, and $\text{Gaia}$. Monoceros is a structure located towards the Galactic anticenter and close to the disk. We identified that its stars have azimuthal velocity in the range of $200 < v_φ\,{\rm(km\,s^{-1})}< 250$. Combining their kinematics and spatial distribution, we desig… ▽ More We explore the kinematic and chemical properties of Monoceros stellar overdensity by combining data from 2MASS, WISE, APOGEE, and $\text{Gaia}$. Monoceros is a structure located towards the Galactic anticenter and close to the disk. We identified that its stars have azimuthal velocity in the range of $200 < v_φ\,{\rm(km\,s^{-1})}< 250$. Combining their kinematics and spatial distribution, we designed a new method to select stars from this overdensity. This method allows us to easily identify the structure in both hemispheres and estimate their distances. Our analysis was supported by comparison with simulated data from the entire sky generated by $\texttt{Galaxia}$ code. Furthermore, we characterized, for the first time, the Monoceros overdensity in several chemical-abundance spaces. Our results confirm its similarity to stars found in the thin disk of the Galaxy and suggest an $\textit{in situ}$ formation. Furthermore, we demonstrate that the southern (Mon-S) and northern (Mon-N) regions of Monoceros exhibit indistinguishable chemical compositions. △ Less

Submitted 18 October, 2023; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: Paper accepted for publication in ApJ

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

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

Structure, Kinematics, and Observability of the Large Magellanic Cloud's Dynamical Friction Wake in Cold vs. Fuzzy Dark Matter

Authors: Hayden R. Foote, Gurtina Besla, Philip Mocz, Nicolás Garavito-Camargo, Lachlan Lancaster, Martin Sparre, Emily C. Cunningham, Mark Vogelsberger, Facundo A. Gómez, Chervin F. P. Laporte

Abstract: The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW's stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of… ▽ More The Large Magellanic Cloud (LMC) will induce a dynamical friction (DF) wake on infall to the Milky Way (MW). The MW's stellar halo will respond to the gravity of the LMC and the dark matter (DM) wake, forming a stellar counterpart to the DM wake. This provides a novel opportunity to constrain the properties of the DM particle. We present a suite of high-resolution, windtunnel-style simulations of the LMC's DF wake that compare the structure, kinematics, and stellar tracer response of the DM wake in cold DM (CDM), with and without self-gravity, vs. fuzzy DM (FDM) with $m_a = 10^{-23}$ eV. We conclude that the self-gravity of the DM wake cannot be ignored. Its inclusion raises the wake's density by $\sim 10\%$, and holds the wake together over larger distances ($\sim$ 50 kpc) than if self-gravity is ignored. The DM wake's mass is comparable to the LMC's infall mass, meaning the DM wake is a significant perturber to the dynamics of MW halo tracers. An FDM wake is more granular in structure and is $\sim 20\%$ dynamically colder than a CDM wake, but with comparable density. The granularity of an FDM wake increases the stars' kinematic response at the percent level compared to CDM, providing a possible avenue of distinguishing a CDM vs. FDM wake. This underscores the need for kinematic measurements of stars in the stellar halo at distances of 70-100 kpc. △ Less

Submitted 8 September, 2023; v1 submitted 30 June, 2023; originally announced July 2023.

Comments: 38 pages, 31 figures. Updated to published version

Journal ref: ApJ 954 163 (2023)

Exploring the diversity and similarity of radially anisotropic Milky Way-like stellar haloes: implications for disrupted dwarf galaxy searches

Authors: Matthew D. A. Orkney, Chervin F. P. Laporte, Robert J. J. Grand, Facundo A. Gómez, Freeke van de Voort, Azadeh Fattahi, Federico Marinacci, Rüdiger Pakmor, Francesca Fragkoudi, Volker Springel

Abstract: We investigate the properties of mergers comparable to the Gaia-Sausage-Enceladus (GSE) using cosmological hydrodynamical simulations of Milky Way-like galaxies. The merger progenitors span an order of magnitude in their peak stellar mass ($3\times10^8<M_{\star}/\rm{M}_{\odot}<4\times10^9$) and include both rotation and pressure-supported galaxies ($0.10<D/T<0.77$). In a minority of cases, the GSE… ▽ More We investigate the properties of mergers comparable to the Gaia-Sausage-Enceladus (GSE) using cosmological hydrodynamical simulations of Milky Way-like galaxies. The merger progenitors span an order of magnitude in their peak stellar mass ($3\times10^8<M_{\star}/\rm{M}_{\odot}<4\times10^9$) and include both rotation and pressure-supported galaxies ($0.10<D/T<0.77$). In a minority of cases, the GSE-like debris is comprised of stars from more than one merger progenitor. However, there is a close similarity in their chemodynamical properties and the triaxial shapes of their debris, and so it is not always possible to distinguish them. The merger progenitors host a variety of luminous satellites ($0-8$ with $M_{\star}>10^6\,\rm{M}_{\odot}$), but most of these do not follow the merger to low orbital energies. Between $0-1$ of these satellites may survive to $z=0$, but with no clear signatures of their past association. We show that the fraction of stars originating from GSE-like mergers is reduced for lower metallicities (reaching a minimum around $\text{[Fe/H]} = -2$), and also within $5\,$kpc of the galactic centre. Whilst these central regions are dominated by in-situ stars, the ex-situ fraction trends towards a 100 per cent asymptote when considering the most metal-poor stars ($\text{[Fe/H]}\ll-2.5$). Considering this, its near proximity, and its small volume on the sky, the Galactic centre lends itself as a prime environment in the search for the stars from the earliest galaxies, whilst avoiding contamination from GSE stars. △ Less

Submitted 7 September, 2023; v1 submitted 3 March, 2023; originally announced March 2023.

Comments: 23 pages, 18 figures, accepted by MNRAS

The phase spiral in Gaia DR3

Authors: T. Antoja, P. Ramos, B. García-Conde, M. Bernet, C. F. P. Laporte, D. Katz

Abstract: We aim to study the phase spiral in the Milky Way (MW) with Gaia DR3. We used an edge detection algorithm to find the border of the phase spiral, allowing us to robustly quantify its shape at different positions and for different selections. We calculated the time of onset of the phase-mixing by determining the different turns of the phase spiral and using the vertical frequencies from commonly us… ▽ More We aim to study the phase spiral in the Milky Way (MW) with Gaia DR3. We used an edge detection algorithm to find the border of the phase spiral, allowing us to robustly quantify its shape at different positions and for different selections. We calculated the time of onset of the phase-mixing by determining the different turns of the phase spiral and using the vertical frequencies from commonly used MW potential models. We find that the phase spiral extends down to $-1.2$ kpc in height below the plane (about 3 to 5 scale heights of the thin disc) and beyond $\pm 50$ km/s in $V_Z$. We see a secondary branch mostly at positive vertical velocities when coloured by azimuthal velocity and in the counts projection. We also find complex variations of the phase spirals with angular momentum and azimuth. All these possibly provide evidence of multiple perturbations (from different times or from different perturbers) and/or of the complexity of the phase mixing process. We detect the phase spiral from 6 to 11 kpc from the Galactic centre and find signatures of vertical asymmetries 1-2 kpc beyond this range. We measure small but clear variations with azimuth. When we determine the phase mixing times from the phase spiral at different angular momenta and using the different spiral turns (at different $Z$) we obtain inconsistent times with systematic differences (times increasing with $|L_Z|$ and with $|Z|$). Our determinations are mostly in the range of [0.3-0.9] Gyr, with an average of 0.5 Gyr. The inconsistencies do not change when using different usual potential models, different stellar distances or frequencies for different kinetic temperatures. They could stem from the inconsistency of potential models with the true MW, and from too simple modelling, in particular neglecting self-gravity, not considering the multiple perturbations and the interference with other processes. △ Less

Submitted 25 May, 2023; v1 submitted 22 December, 2022; originally announced December 2022.

Comments: version after proofs corrections

Journal ref: A&A 673, A115 (2023)

Lopsided Galaxies in a cosmological context: a new galaxy-halo connection

Authors: Silvio Varela-Lavin, Facundo A. Gómez, Patricia B. Tissera, Gurtina Besla, Nicolás Garavito-Camargo, Federico Marinacci, Chervin F. P. Laporte

Abstract: Disc galaxies commonly show asymmetric features in their morphology, such as warps and lopsidedness. These features can provide key information regarding the recent evolution of a given disc galaxy. In the nearby Universe, up to $\sim30$ percent of late-type galaxies display a global non-axisymmetric lopsided mass distribution. However, the origin of this perturbation is not well understood. In th… ▽ More Disc galaxies commonly show asymmetric features in their morphology, such as warps and lopsidedness. These features can provide key information regarding the recent evolution of a given disc galaxy. In the nearby Universe, up to $\sim30$ percent of late-type galaxies display a global non-axisymmetric lopsided mass distribution. However, the origin of this perturbation is not well understood. In this work, we study the origin of lopsided perturbations in simulated disc galaxies extracted from the TNG50 simulation of the IllustrisTNG project. We statistically explore different excitation mechanisms for this perturbation, such as direct satellite tidal interactions and distortions of the underlying dark matter distributions. We also characterize the main physical conditions that lead to lopsided perturbations. 50 percent of our sample galaxy have lopsided modes $m=1$ greater than $\sim 0.12$. We find a strong correlation between internal galaxy properties, such as central stellar surface density and disc radial extension with the strength of lopsided modes. The majority of lopsided galaxies have lower central surface densities and more extended discs than symmetric galaxies. As a result, such lopsided galaxies are less self-gravitationally cohesive, and their outer disc region is more susceptible to different types of external perturbations. However, we do not find strong evidence that tidal interactions with satellite galaxies are the main driving agent of lopsided modes. Lopsided galaxies tend to live in asymmetric dark matter halos with high spin, indicating strong galaxy-halo connections in late-type lopsided galaxies. △ Less

Submitted 14 June, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: 16 pages, 16 figures, accepted for publication in MNRAS

The Prince and the Pauper: Evidence for the early high-redshift formation of the Galactic $α$-poor disc population

Authors: Matthew Raymond Gent, Philipp Eitner, Aldo Serenelli, Jennifer K. S. Friske, Sergey E. Koposov, Chervin F. P. Laporte, Tobias Buck, Maria Bergemann

Abstract: Context. The presence of [$α$/Fe]-[Fe/H] bi-modality in the Milky Way disc has animated the Galactic archaeology community since more than two decades. Aims. Our goal is to investigate the chemical, temporal, and kinematical structure of the Galactic discs using abundances, kinematics, and ages derived self-consistently with the new Bayesian framework SAPP. Methods. We employ the public Gaia-ESO s… ▽ More Context. The presence of [$α$/Fe]-[Fe/H] bi-modality in the Milky Way disc has animated the Galactic archaeology community since more than two decades. Aims. Our goal is to investigate the chemical, temporal, and kinematical structure of the Galactic discs using abundances, kinematics, and ages derived self-consistently with the new Bayesian framework SAPP. Methods. We employ the public Gaia-ESO spectra, as well as Gaia EDR3 astrometry and photometry. Stellar parameters and chemical abundances are determined for 13 426 stars using NLTE models of synthetic spectra. Ages are derived for a sub-sample of 2 898 stars, including subgiants and main-sequence stars. The sample probes a large range of Galactocentric radii, $\sim$ 3 to 12 kpc, and extends out of the disc plane to $\pm$ 2 kpc. Results. Our new data confirm the known bi-modality in the [Fe/H] - [$α$/Fe] space, which is often viewed as the manifestation of the chemical thin and thick discs. The over-densities significantly overlap in metallicity, age, and kinematics, and none of these is a sufficient criterion for distinguishing between the two disc populations. Different from previous studies, we find that the $α$-poor disc population has a very extended [Fe/H] distribution and contains $\sim$ 20$\%$ old stars with ages of up to $\sim$ 11 Gyr. Conclusions. Our results suggest that the Galactic thin disc was in place early, at look-back times corresponding to redshifts z $\sim$ 2 or more. At ages $\sim$ 9 to 11 Gyr, the two disc structures shared a period of co-evolution. Our data can be understood within the clumpy disc formation scenario that does not require a pre-existing thick disc to initiate a formation of the thin disc. We anticipate that a similar evolution can be realised in cosmological simulations of galaxy formation. △ Less

Submitted 14 November, 2023; v1 submitted 22 June, 2022; originally announced June 2022.

Comments: 14 pages, 10 figures, re-submitted to A&A

The impact of two massive early accretion events in a Milky Way-like galaxy: repercussions for the buildup of the stellar disc and halo

Authors: Matthew D. A. Orkney, Chervin F. P. Laporte, Robert J. J. Grand, Facundo A. Gómez, Freeke van de Voort, Federico Marinacci, Ruediger Pakmor, Volker Springel

Abstract: We identify and characterise a Milky Way-like realisation from the Auriga simulations with two consecutive massive mergers $\sim2\,$Gyr apart at high redshift, comparable to the reported Kraken and Gaia-Sausage-Enceladus. The Kraken-like merger ($z=1.6$, $M_{\rm Tot} = 8\times10^{10}\,$M$_{\odot}$) is gas-rich, deposits most of its mass in the inner $10\,$kpc, and is largely isotropic. The Sausage… ▽ More We identify and characterise a Milky Way-like realisation from the Auriga simulations with two consecutive massive mergers $\sim2\,$Gyr apart at high redshift, comparable to the reported Kraken and Gaia-Sausage-Enceladus. The Kraken-like merger ($z=1.6$, $M_{\rm Tot} = 8\times10^{10}\,$M$_{\odot}$) is gas-rich, deposits most of its mass in the inner $10\,$kpc, and is largely isotropic. The Sausage-like merger ($z=1.14$, $M_{\rm Tot} = 1\times10^{11}\,$M$_{\odot}$) leaves a more extended mass distribution at higher energies, and has a radially anisotropic distribution. For the higher redshift merger, the stellar mass ratio of the satellite to host galaxy is 1:3. As a result, the chemistry of the remnant is indistinguishable from contemporaneous in-situ populations, making it challenging to identify this component through chemical abundances. This naturally explains why all abundance patterns attributed so far to Kraken are in fact fully consistent with the metal-poor in-situ so-called Aurora population and thick disc. However, our model makes a falsifiable prediction: if the Milky Way underwent a gas-rich double merger at high redshift, then this should be imprinted on its star formation history with bursts about $\sim2\,$Gyrs apart. This may offer constraining power on the highest-redshift major mergers. △ Less

Submitted 7 November, 2022; v1 submitted 18 June, 2022; originally announced June 2022.

Comments: 5 pages, 5 figures, published by MNRAS letters. Comments welcome

Journal ref: Volume 517, Issue 1, November 2022, Pages L138-L142

On the Hunt for the Origins of the Orphan--Chenab Stream: Detailed Element Abundances with APOGEE and Gaia

Authors: Keith Hawkins, Adrian M. Price-Whelan, Allyson A. Sheffield, Aidan Z. Subrahimovic, Rachael L. Beaton, Vasily Belokurov, Denis Erkal, Sergey E. Koposov, Richard R. Lane, Chervin F. P. Laporte, Christian Nitschelm

Abstract: Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan--Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the… ▽ More Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan--Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the parent of the OC stream, we use astrometry from the early third data release of ESA's Gaia mission (Gaia EDR3) and radial velocity information from the SDSS-IV APOGEE survey to find up to 13 stars that are likely members of the OC stream. We use the APOGEE survey to study the chemical nature (for up to 13 stars) of the OC stream in the $α$ (O, Mg, Ca, Si, Ti, S), odd-Z (Al, K, V), Fe-peak (Fe, Ni, Mn, Co, Cr) and neutron capture (Ce) elemental groups. We find that the stars that make up the OC stream are not consistent with a mono-metallic population and have a median metallicity of --1.92~dex with a dispersion of 0.28 dex. Our results also indicate that the $α$-elements are depleted compared to the known Milky Way populations and that its [Mg/Al] abundance ratio is not consistent with second generation stars from globular clusters. The detailed chemical pattern of these stars indicates that the OC stream progenitor is very likely to be a dwarf spheroidal galaxy with a mass of ~10$^6$ M$_\odot$. △ Less

Submitted 27 May, 2022; originally announced May 2022.

Comments: 13 Pages; 4 Figures. Submitted to AAS Journals, comments welcome

GASTRO library I: the simulated chemodynamical properties of several GSE-like stellar halos

Authors: João A. S. Amarante, Victor P. Debattista, Leandro Beraldo e Silva, Chervin F. P. Laporte, Nathan Deg

Abstract: The Milky Way stellar halo contains relics of ancient mergers that tell the story of our Galaxy's formation. Some of them are identified due to their similarity in energy, actions and chemistry, referred to as the "chemodynamical space", and are often attributed to distinct merger events. It is also known that our Galaxy went through a significant merger event that shaped the local stellar halo du… ▽ More The Milky Way stellar halo contains relics of ancient mergers that tell the story of our Galaxy's formation. Some of them are identified due to their similarity in energy, actions and chemistry, referred to as the "chemodynamical space", and are often attributed to distinct merger events. It is also known that our Galaxy went through a significant merger event that shaped the local stellar halo during its first Gyr. Previous studies using $N$-body only and cosmological hydrodynamical simulations have shown that such single massive merger can produce several "signatures" in the chemodynamical space, which can potentially be misinterpreted as distinct merger events. Motivated by these, in this work we use a subset of the GASTRO, library which consists of several SPH+$N$-body models of single accretion event in a Milky Way-like galaxy. Here, we study models with orbital properties similar to the main merger event of our Galaxy and explore the implications to known stellar halo substructures. We find that: $i.$ supernova feedback efficiency influences the satellite's structure and orbital evolution, resulting in distinct chemodynamical features for models with the same initial conditions, $ii.$ very retrograde high energy stars are the most metal-poor of the accreted dwarf galaxy and could be misinterpreted as a distinct merger $iii.$ the most bound stars are more metal-rich in our models, the opposite of what is observed in the Milky Way, suggesting a secondary massive merger, and finally $iv.$ our models can reconcile other known substructures to an unique progenitor. △ Less

Submitted 17 September, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: Published in ApJ

Journal ref: ApJ 937 12 (2022)

Weighing the Galactic disk using phase-space spirals IV. Tests on a three-dimensional galaxy simulation

Authors: Axel Widmark, Jason A. S. Hunt, Chervin F. P. Laporte, Giacomo Monari

Abstract: In this fourth article on weighing the Galactic disk using the shape of the phase-space spiral, we have tested our method on a billion particle three-dimensional N-body simulation, comprised of a Milky Way like host galaxy and a merging dwarf satellite. The main purpose of this work was to test the validity of our model's fundamental assumptions that the spiral inhabits a locally static and vertic… ▽ More In this fourth article on weighing the Galactic disk using the shape of the phase-space spiral, we have tested our method on a billion particle three-dimensional N-body simulation, comprised of a Milky Way like host galaxy and a merging dwarf satellite. The main purpose of this work was to test the validity of our model's fundamental assumptions that the spiral inhabits a locally static and vertically separable gravitational potential. These assumptions might be compromised in the complex kinematic system of a disturbed three-dimensional disk galaxy; in fact, the statistical uncertainty and any potential biases related to these assumptions are expected to be amplified for this simulation, which differs from the Milky Way in that it is more strongly perturbed and has a phase-space spiral that inhabits higher vertical energies. We constructed 44 separate data samples from different spatial locations in the simulated host galaxy. Our method produced accurate results for the vertical gravitational potential of these 44 data samples, with an unbiased distribution of errors with a standard deviation of 7 %. We also tested our method under severe and unknown spatially dependent selection effects, also with robust results; this sets it apart from traditional dynamical mass measurements that are based on the assumption of a steady state and which are highly sensitive to unknown or poorly modelled incompleteness. Hence, we will be able to make localised mass measurements of distant regions in the Milky Way disk, which would otherwise be compromised by complex and poorly understood selection effects. △ Less

Submitted 11 April, 2022; v1 submitted 12 January, 2022; originally announced January 2022.

Comments: 18 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 663, A16 (2022)

Migration and Mixing in the Galactic Disc from Encounters between Sagittarius and the Milky Way

Authors: Christopher Carr, Kathryn V. Johnston, Chervin F. P. Laporte, Melissa K. Ness

Abstract: Stars born on near-circular orbits in spiral galaxies can subsequently migrate to different orbits due to interactions with non-axisymmetric disturbances within the disc such as bars or spiral arms. This paper extends the study of migration to examine the role of external influences using the example of the interaction of the Sagittarius dwarf galaxy (Sgr) with the Milky Way (MW). We first make im… ▽ More Stars born on near-circular orbits in spiral galaxies can subsequently migrate to different orbits due to interactions with non-axisymmetric disturbances within the disc such as bars or spiral arms. This paper extends the study of migration to examine the role of external influences using the example of the interaction of the Sagittarius dwarf galaxy (Sgr) with the Milky Way (MW). We first make impulse approximation estimates to characterize the influence of Sgr disc passages. The tidal forcing from Sgr can produce changes in both guiding radius ($ΔR_g$) and orbital eccentricity, as quantified by the maximum radial excursion, $ΔR_ {\rm max} $. These changes follow a quadrupole-like pattern across the face of the disc, with amplitude increasing with Galactocentric radius. We next examine a collisionless N-body simulation of a Sgr-like satellite interacting with a MW-like galaxy and find that Sgr's influence in the outer disc dominates over the secular evolution of orbits between disc passages. Finally, we use the same simulation to explore possible observable signatures of Sgr-induced migration by painting the simulation with different age stellar populations. We find that following Sgr disc passages, the migration it induces manifests within an annulus as an approximate quadrupole in azimuthal metallicity variations ($δ_ {\rm [Fe/H]} $), along with systematic variations in orbital eccentricity, $ΔR_ {\rm max} $. These systematic variations can persist for several rotational periods. We conclude that this combination of signatures may be used to distinguish between the different migration mechanisms shaping the chemical abundance patterns of the Milky Way's thin disc. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 17 pages, 13 figures, submitted to MNRAS

Estimating the local dark matter density in a non-axisymmetric wobbling disc

Authors: S. Sivertsson, J. I. Read, H. Silverwood, P. F. de Salas, K. Malhan, A. Widmark, C. F. P. Laporte, S. Garbari, K. Freese

Abstract: The density of dark matter near the Sun is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way's dark matter halo. Estimates to date have typically assumed that the Milky Way's stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertica… ▽ More The density of dark matter near the Sun is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way's dark matter halo. Estimates to date have typically assumed that the Milky Way's stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. We assess the impact of these assumptions on determinations of the local dark matter density by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia-Sausage-Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans-Poisson equations on our ability to correctly extract the local dark matter density from the simulated data. We find that common approximations employed in the literature - axisymmetry and a locally flat rotation curve - can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2kpc above the disc plane, implying a fractional error on the local dark matter density of order unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate, consistent with the true value within our 95% confidence intervals for realistic 20% uncertainties on the baryonic surface density of the disc. Other terms - the axial tilt, 2:nd Poisson and time dependent terms - contribute less than 10% to the local dark matter density (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis. △ Less

Submitted 5 January, 2022; originally announced January 2022.

Comments: 15 pages, 16 figues. Submitted to MNRAS. Comments welcome!

Weighing the Galactic disk using phase-space spirals III. Probing distant regions of the disk using the Gaia EDR3 proper motion sample

Authors: Axel Widmark, Chervin F. P. Laporte, Giacomo Monari

Abstract: We have applied our method for weighing the Galactic disk using phase-space spirals to the Gaia EDR3 proper motion sample. For stars in distant regions of the Galactic disk, the latitudinal proper motion has a close projection with vertical velocity, such that the phase-space spiral in the plane of vertical position and vertical velocity can be observed without requiring that all stars have availa… ▽ More We have applied our method for weighing the Galactic disk using phase-space spirals to the Gaia EDR3 proper motion sample. For stars in distant regions of the Galactic disk, the latitudinal proper motion has a close projection with vertical velocity, such that the phase-space spiral in the plane of vertical position and vertical velocity can be observed without requiring that all stars have available radial velocity information. We divided the Galactic plane into 360 separate data samples, each corresponding to an area cell in the Galactic plane in the distance range of 1.4-3.4 kpc, with an approximate cell length of 200-400 pc. Roughly half of our data samples were disqualified altogether due to severe selection effects, especially in the direction of the Galactic centre. In the remainder, we were able to infer the vertical gravitational potential by fitting an analytic model of the phase-space spiral to the data. This work is the first of its kind, in the sense that we are weighing distant regions of the Galactic disk with a high spatial resolution, without relying on the strong assumptions of axisymmetry. Post-inference, we fit a thin disk scale length of $2.2\pm 0.1$ kpc, although this value is sensitive to the considered spatial region. We see surface density variations as a function of azimuth on the order of 10-20 %, which is roughly the size of our estimated sum of potential systematic biases. With this work, we have demonstrated that our method can be used to weigh distant regions of the Galactic disk despite strong selection effects. We expect to reach even greater distances and improve our accuracy with future Gaia data releases and further improvements to our method. △ Less

Submitted 11 April, 2022; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 18 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 663, A15 (2022)

Not the Birth Cluster: the Stellar Clustering that Shapes Planetary Systems is Generated by Galactic-Dynamical Perturbations

Authors: J. M. Diederik Kruijssen, Steven N. Longmore, Mélanie Chevance, Chervin F. P. Laporte, Michal Motylinski, Benjamin W. Keller, Jonathan D. Henshaw

Abstract: Recent work has demonstrated that exoplanetary system properties correlate strongly with ambient stellar clustering in six-dimensional stellar position-velocity phase space, quantified by dividing planetary systems into sub-samples with high or low phase space densities (`overdensity' and `field' systems, respectively). We investigate the physical origins of the phase space overdensities and, ther… ▽ More Recent work has demonstrated that exoplanetary system properties correlate strongly with ambient stellar clustering in six-dimensional stellar position-velocity phase space, quantified by dividing planetary systems into sub-samples with high or low phase space densities (`overdensity' and `field' systems, respectively). We investigate the physical origins of the phase space overdensities and, thereby, which environmental mechanisms may have impacted the planetary systems. We consider the galactic-scale kinematic structure of the Milky Way observed with Gaia and show that the overdensities correspond to the well-known, kpc-scale kinematic ripples and streams in the Galactic disk, which are thought to be generated by bar and spiral arm-driven resonances and satellite galaxy passages. We also find indications that the planet demographics may vary between individual phase space overdensities, which potentially have differing physical origins and histories. Planetary systems associated with the `phase space spiral' (a recent perturbation of the Galactic disk) have a hot Jupiter-to-cold Jupiter ratio that is 10 times higher than in field systems. Finally, the hot Jupiter-to-cold Jupiter ratio within overdensities may increase with host stellar age over Gyr timescales. Because the overdensities persist for several Gyr, we argue that late-time perturbations of planetary systems most likely explain these trends, although additional perturbations at birth may contribute too. This suggests that planetary system properties are not just affected by stellar clustering in their immediate surroundings, but by galaxy-scale processes throughout their evolution. We conclude by discussing the main open questions towards understanding the diversity of physical processes that together set planetary system architectures. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 20 pages, 9 figures; ApJ submitted (September 6, 2021)

The Clustering of Orbital Poles Induced by the LMC: Hints for the Origin of Planes of Satellites

Authors: Nicolas Garavito-Camargo, Ekta Patel, Gurtina Besla, Adrian M. Price-Whelan, Facundo A. Gomez, Chervin F. P Laporte, Kathryn V. Johnston

Abstract: A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the… ▽ More A significant fraction of Milky Way (MW) satellites exhibit phase-space properties consistent with a coherent orbital plane. Using tailored N--body simulations of a spherical MW halo that recently captured a massive (1.8$\times 10^{11}$M$\odot$) LMC-like satellite, we identify the physical mechanisms that may enhance the clustering of orbital poles of objects orbiting the MW. The LMC deviates the orbital poles of MW dark matter (DM) particles from the present-day random distribution. Instead, the orbital poles of particles beyond $R\approx 50$kpc cluster near the present-day orbital pole of the LMC along a sinusoidal pattern across the sky. The density of orbital poles is enhanced near the LMC by a factor $δρ_{max}$=30\%(50\%) with respect to underdense regions, and $δρ_{iso}$=15\%(30\%) relative to the isolated MW simulation (no LMC) between 50-150 kpc (150-300 kpc). The clustering appears after the LMC's pericenter ($\approx$ 50 Myr ago, 49 kpc) and lasts for at least 1 Gyr. Clustering occurs because of three effects: 1) the LMC shifts the velocity and position of the central density of the MW's halo and disk; 2) the DM dynamical friction wake and collective response induced by the LMC changes the kinematics of particles; 3) observations of particles selected within spatial planes suffer from a bias, such that measuring orbital poles in a great circle in the sky enhances the probability of their orbital poles being clustered. This scenario should be ubiquitous in hosts that recently captured a massive satellite (at least $\approx$ 1:10 mass ratio), causing the clustering of orbital poles of halo tracers. △ Less

Submitted 16 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ. 23 pages and 16 figures. Comments are welcome

Resolving local and global kinematic signatures of satellite mergers with billion particle simulations

Authors: Jason A. S. Hunt, Ioana A. Stelea, Kathryn V. Johnston, Suroor S. Gandhi, Chervin F. P. Laporte, Jeroen Bedorf

Abstract: In this work we present two new $\sim10^9$ particle self-consistent simulations of the merger of a Sagittarius-like dwarf galaxy with a Milky Way-like disc galaxy. One model is a violent merger creating a thick disc, and a Gaia-Enceladus/Sausage like remnant. The other is a highly stable disc which we use to illustrate how the improved phase space resolution allows us to better examine the formati… ▽ More In this work we present two new $\sim10^9$ particle self-consistent simulations of the merger of a Sagittarius-like dwarf galaxy with a Milky Way-like disc galaxy. One model is a violent merger creating a thick disc, and a Gaia-Enceladus/Sausage like remnant. The other is a highly stable disc which we use to illustrate how the improved phase space resolution allows us to better examine the formation and evolution of structures that have been observed in small, local volumes in the Milky Way, such as the $z-v_z$ phase spiral and clustering in the $v_{\mathrm{R}}-v_φ$ plane when compared to previous works. The local $z-v_z$ phase spirals are clearly linked to the global asymmetry across the disc: we find both 2-armed and 1-armed phase spirals, which are related to breathing and bending behaviors respectively. Hercules-like moving groups are common, clustered in $v_{\mathrm{R}}-v_φ$ in local data samples in the simulation. These groups migrate outwards from the inner galaxy, matching observed metallicity trends even in the absence of a galactic bar. We currently release the best fitting `present day' merger snapshots along with the unperturbed galaxies for comparison. △ Less

Submitted 11 September, 2021; v1 submitted 13 July, 2021; originally announced July 2021.

Comments: 16 pages, 14 figures, accepted by MNRAS

Snails Across Scales: Local and Global Phase-Mixing Structures as Probes of the Past and Future Milky Way

Authors: Suroor S. Gandhi, Kathryn V. Johnston, Jason A. S. Hunt, Adrian M. Price-Whelan, Chervin F. P. Laporte, David W. Hogg

Abstract: Signatures of vertical disequilibrium have been observed across the Milky Way's disk. These signatures manifest locally as unmixed phase-spirals in $z$--$v_z$ space ("snails-in-phase") and globally as nonzero mean $z$ and $v_z$ which wraps around as a physical spiral across the $x$--$y$ plane ("snails-in-space"). We explore the connection between these local and global spirals through the example… ▽ More Signatures of vertical disequilibrium have been observed across the Milky Way's disk. These signatures manifest locally as unmixed phase-spirals in $z$--$v_z$ space ("snails-in-phase") and globally as nonzero mean $z$ and $v_z$ which wraps around as a physical spiral across the $x$--$y$ plane ("snails-in-space"). We explore the connection between these local and global spirals through the example of a satellite perturbing a test-particle Milky Way (MW)-like disk. We anticipate our results to broadly apply to any vertical perturbation. Using a $z$--$v_z$ asymmetry metric we demonstrate that in test-particle simulations: (a) multiple local phase-spiral morphologies appear when stars are binned by azimuthal action $J_φ$, excited by a single event (in our case, a satellite disk-crossing); (b) these distinct phase-spirals are traced back to distinct disk locations; and (c) they are excited at distinct times. Thus, local phase-spirals offer a global view of the MW's perturbation history from multiple perspectives. Using a toy model for a Sagittarius (Sgr)-like satellite crossing the disk, we show that the full interaction takes place on timescales comparable to orbital periods of disk stars within $R \lesssim 10$ kpc. Hence such perturbations have widespread influence which peaks in distinct regions of the disk at different times. This leads us to examine the ongoing MW-Sgr interaction. While Sgr has not yet crossed the disk (currently, $z_{Sgr} \approx -6$ kpc, $v_{z,Sgr} \approx 210$ km/s), we demonstrate that the peak of the impact has already passed. Sgr's pull over the past 150 Myr creates a global $v_z$ signature with amplitude $\propto M_{Sgr}$, which might be detectable in future spectroscopic surveys. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Journal ref: ApJ 928 80 (2022)

Weighing the Galactic disk using phase-space spirals II. Most stringent constraints to a thin dark disk using Gaia EDR3

Authors: Axel Widmark, Chervin F. P. Laporte, Pablo F. de Salas, Giacomo Monari

Abstract: Using the method that was developed in the first paper of this series, we measure the vertical gravitational potential of the Galactic disk from the time-varying structure of the phase-space spiral, using data from Gaia as well as supplementary radial velocity information from legacy spectroscopic surveys. For eleven independent data samples, we inferred gravitational potentials that were in good… ▽ More Using the method that was developed in the first paper of this series, we measure the vertical gravitational potential of the Galactic disk from the time-varying structure of the phase-space spiral, using data from Gaia as well as supplementary radial velocity information from legacy spectroscopic surveys. For eleven independent data samples, we inferred gravitational potentials that were in good agreement, despite the data samples' varied and substantial selection effects. Using a model for the baryonic matter densities, we inferred a local halo dark matter density of $0.0085 \pm 0.0039$ M$_\odot$/pc$^3 = 0.32 \pm 0.15$ GeV/cm$^3$. We were also able to place the most stringent constraint to the surface density of a thin dark disk with a scale height $\leq 50$ pc: an upper 95 % confidence limit of roughly 5 M$_\odot$/pc$^2$ (compared to previous limit of roughly 10 M$_\odot$/pc$^2$, given the same scale height). For the inferred halo dark matter density and thin dark disk surface density, the uncertainties are dominated by the baryonic model. With this level of precision, our method is highly competitive with traditional methods that rely on the assumption of a steady state. In a general sense, this illustrates that time-varying dynamical structures are not solely obstacles to dynamical mass measurements, but can also be regarded as assets containing useful information. △ Less

Submitted 30 November, 2021; v1 submitted 28 May, 2021; originally announced May 2021.

Comments: 9 figures, 20 pages

Journal ref: A&A 653, A86 (2021)

Hidden depths in the local Universe: The Stellar Stream Legacy Survey

Authors: David Martinez-Delgado, Andrew P. Cooper, Javier Roman, Annalisa Pillepich, Denis Erkal, Sarah Pearson, John Moustakas, Chervin F. P. Laporte, Seppo Laine, Mohammad Akhlaghi, Dustin Lang, Dmitry Makarov, Alejandro S. Borlaff, Giuseppe Donatiello, William J. Pearson, Juan Miro-Carretero, Jean-Charles Cuillandre, Helena Dominguez, Santi Roca-Fabrega, Carlos S. Frenk, Judy Schmidt, Maria A. Gomez-Flechoso, Rafael Guzman, Noam I. Libeskind, Arjun Dey , et al. (4 additional authors not shown)

Abstract: Mergers and tidal interactions between massive galaxies and their dwarf satellites are a fundamental prediction of the Lambda-Cold Dark Matter cosmology. These events are thought to provide important observational diagnostics of nonlinear structure formation. Stellar streams in the Milky Way and Andromeda are spectacular evidence for ongoing satellite disruption. However, constructing a statistica… ▽ More Mergers and tidal interactions between massive galaxies and their dwarf satellites are a fundamental prediction of the Lambda-Cold Dark Matter cosmology. These events are thought to provide important observational diagnostics of nonlinear structure formation. Stellar streams in the Milky Way and Andromeda are spectacular evidence for ongoing satellite disruption. However, constructing a statistically meaningful sample of tidal streams beyond the Local Group has proven a daunting observational challenge, and the full potential for deepening our understanding of galaxy assembly using stellar streams has yet to be realised. Here we introduce the Stellar Stream Legacy Survey, a systematic imaging survey of tidal features associated with dwarf galaxy accretion around a sample of ~3100 nearby galaxies within z~0.02, including about 940 Milky Way analogues. Our survey exploits public deep imaging data from the DESI Legacy Imaging Surveys, which reach surface brightness as faint as ~29 mag/arcsec^2 in the r band. As a proof of concept of our survey, we report the detection and broad-band photometry of 24 new stellar streams in the local Universe. We discuss how these observations can yield new constraints on galaxy formation theory through comparison to mock observations from cosmological galaxy simulations. These tests will probe the present-day mass assembly rate of galaxies, the stellar populations and orbits of satellites, the growth of stellar halos and the resilience of stellar disks to satellite bombardment. △ Less

Submitted 20 December, 2022; v1 submitted 13 April, 2021; originally announced April 2021.

Comments: Accepted in Astronomy & Astrophysics, 25 pages, 17 figures, 3 tables

Journal ref: A&A 671, A141 (2023)

Kinematics beats dust: unveiling nested substructure in the perturbed outer disc of the Milky Way

Authors: Chervin F. P. Laporte, Sergey E. Koposov, Vasily Belokurov

Abstract: We use $Gaia$ eDR3 data and legacy spectroscopic surveys to map the Milky Way disc substructure towards the Galactic Anticenter at heliocentric distances $d\geq10\,\rm{kpc}$. We report the discovery of multiple previously undetected new filaments embedded in the outer disc in highly extincted regions. Stars in these over-densities have distance gradients expected for disc material and move on disc… ▽ More We use $Gaia$ eDR3 data and legacy spectroscopic surveys to map the Milky Way disc substructure towards the Galactic Anticenter at heliocentric distances $d\geq10\,\rm{kpc}$. We report the discovery of multiple previously undetected new filaments embedded in the outer disc in highly extincted regions. Stars in these over-densities have distance gradients expected for disc material and move on disc-like orbits with $v_φ\sim170-230\,\rm{km\,s^{-1}}$, showing small spreads in energy. Such a morphology argues against a quiescently growing Galactic thin disc. Some of these structures are interpreted as excited outer disc material, kicked up by satellite impacts and currently undergoing phase-mixing ("feathers"). Due to the long timescale in the outer disc regions, these structures can stay coherent in configuration space over several Gyrs. We nevertheless note that some of these structures could also be folds in the perturbed disc seen in projection from the Sun's location. A full 6D phase-space characterization and age dating of these structure should help distinguish between the two possible morphologies. △ Less

Submitted 23 March, 2021; originally announced March 2021.

Comments: 5 pages, 4 figures, submitted to MNRAS

Weighing the Galactic disk using phase-space spirals I. Tests on one-dimensional simulations

Authors: Axel Widmark, Chervin F. P. Laporte, Pablo F. de Salas

Abstract: We present a new method for inferring the gravitational potential of the Galactic disk, using the time-varying structure of a phase-space spiral in the $(z,w)$-plane (where $z$ and $w$ represent vertical position and vertical velocity). Our method of inference extracts information from the shape of the spiral and disregards the bulk density distribution that is usually used to perform dynamical ma… ▽ More We present a new method for inferring the gravitational potential of the Galactic disk, using the time-varying structure of a phase-space spiral in the $(z,w)$-plane (where $z$ and $w$ represent vertical position and vertical velocity). Our method of inference extracts information from the shape of the spiral and disregards the bulk density distribution that is usually used to perform dynamical mass measurements. In this manner, it is complementary to traditional methods that are based on the assumption of a steady state. Our method consists of fitting an analytical model for the phase-space spiral to data, where the spiral is seen as a perturbation of the stellar number density in the $(z,w)$-plane. We tested our method on one-dimensional simulations, which were initiated in a steady state and then perturbed by an external force similar to that of a passing satellite. We were able to retrieve the true gravitational potentials of the simulations with high accuracy. The gravitational potential at 400-500 parsec distances from the disk mid-plane was inferred with an error of only a few percent. This is the first paper of a series in which we plan to test and refine our method on more complex simulations, as well as apply our method to Gaia data. △ Less

Submitted 29 November, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: 16 pages, 12 figures

Journal ref: A&A 650, A124 (2021)

A Holistic Review of a Galactic Interaction

Authors: Douglas Grion Filho, Kathryn V. Johnston, Eloisa Poggio, Chervin F. P. Laporte, Ronald Drimmel, Elena D'Onghia

Abstract: Our situation as occupants of the Milky Way (MW) Galaxy, bombarded by the Sagittarius dwarf galaxy, provides an intimate view of physical processes that can lead to the dynamical heating of a galactic disc. While this evolution is instigated by Sagittarius, it is also driven by the intertwined influences of the dark matter halo and the disc itself. We analyse an N-body simulation following a Sagit… ▽ More Our situation as occupants of the Milky Way (MW) Galaxy, bombarded by the Sagittarius dwarf galaxy, provides an intimate view of physical processes that can lead to the dynamical heating of a galactic disc. While this evolution is instigated by Sagittarius, it is also driven by the intertwined influences of the dark matter halo and the disc itself. We analyse an N-body simulation following a Sagittarius-like galaxy interacting with a MW-like host to disentangle these different influences during the stages of a minor merger. The accelerations in the disc plane from each component are calculated for each snapshot in the simulation, and then decomposed into Fourier series on annuli. The analysis maps quantify and compare the scales of the individual contributions over space and through time: (i) accelerations due to the satellite are only important around disc passages; (ii) the influence around these passages is enhanced and extended by the distortion of the dark matter halo; (iii) the interaction drives disc asymmetries within and perpendicular to the plane and the self-gravity of these distortions increase in importance with time eventually leading to the formation of a bar. These results have interesting implications for identifying different influences within our own Galaxy. Currently, Sagittarius is close enough to a plane crossing to search for localized signatures of its effect at intermediate radii, the distortion of the MW's dark matter halo should leave its imprint in the outer disc and the disc's own self-consistent response is sculpting the intermediate and inner disc. △ Less

Submitted 8 October, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: 19 pages, 18 figures, published in the Monthly Notices of the Royal Astronomical Society

Journal ref: Volume 507, Issue 2, October 2021

Measuring the vertical response of the Galactic disc to an infalling satellite

Authors: Eloisa Poggio, Chervin F. P. Laporte, Kathryn V. Johnston, Elena D'Onghia, Ronald Drimmel, Douglas Grion Filho

Abstract: Using N-body simulations of the Milky Way interacting with a satellite similar to the Sagittarius dwarf galaxy, we quantitatively analyse the vertical response of the Galactic disc to the satellite's repeated impacts. We approximate the vertical distortion of the Galactic disc as the sum of the first three Fourier azimuthal terms m = 0, 1 and 2, and observe their evolution in different dynamical r… ▽ More Using N-body simulations of the Milky Way interacting with a satellite similar to the Sagittarius dwarf galaxy, we quantitatively analyse the vertical response of the Galactic disc to the satellite's repeated impacts. We approximate the vertical distortion of the Galactic disc as the sum of the first three Fourier azimuthal terms m = 0, 1 and 2, and observe their evolution in different dynamical regimes of interaction. After the first interaction, the m=0 term manifests itself as outgoing ring-like vertical distortions. The m=1 term (S-shape warp) is prograde when the impacts of the satellite are more frequent, or in general close to an interaction, whereas it is slowly retrograde in the most quiescent phases. The m=2 term is typically prograde, and close to an interaction it couples with the m=1 term. Finally, we find that the vertical response of the disc can be recovered in an unbiased way using the instantaneous positions and velocities of stars in a limited volume of the Galactic disc, analogous to real data, and that the measured vertical pattern speeds have a constraining power in the context of a Milky Way-satellite interaction. △ Less

Submitted 23 November, 2020; originally announced November 2020.

Comments: 24 pages, 21 Figures. Submitted to MNRAS. We welcome comments, questions, suggestions for missing references, etc

The outer disc in shambles: blind detection of Monoceros and ACS with Gaia's astrometric sample

Authors: P. Ramos, T. Antoja, C. Mateu, F. Anders, C. F. P. Laporte, J. A. Carballo-Bello, B. Famaey, R. Ibata

Abstract: The astrometric sample of Gaia allows us to study the outermost Galactic disc, the halo and their interface. It is precisely at the very edge of the disc where the effects of external perturbations are expected to be the most noticeable. Our goal is to detect the kinematic substructure present in the halo and at the edge of the Milky Way (MW) disc, and provide observational constraints on their ph… ▽ More The astrometric sample of Gaia allows us to study the outermost Galactic disc, the halo and their interface. It is precisely at the very edge of the disc where the effects of external perturbations are expected to be the most noticeable. Our goal is to detect the kinematic substructure present in the halo and at the edge of the Milky Way (MW) disc, and provide observational constraints on their phase-space distribution. We download, one HEALpix at a time, the proper motion histogram of distant stars, to which we apply a Wavelet Transformation to reveal the significant overdensities. We then analyse the large coherent structures that appear in the sky. We reveal a sharp yet complex anticentre dominated by Monoceros (MNC) and the Anticentre Stream (ACS) in the north, which we find with an intensity comparable to the Magellanic clouds and the Sagittarius stream, and by MNC south and TriAnd at negative latitudes. Our method allows us to perform a morphological analysis of MNC and ACS, both spanning more than 100$^\circ$ in longitude, and to provide a high purity sample of giants with which we track MNC down to latitudes as low as $\sim$5$^\circ$. Their colour-magnitude diagram is consistent with extended structures at a distance of $\sim$10-11 kpc originated in the disc, with a very low ratio of RR Lyrae over M giants, and kinematics compatible with the rotation curve at those distances or only slightly slower. We present a precise characterisation of MNC and ACS, two previously known structures that our method reveals naturally, allowing us to detect them without limiting ourselves to a particular stellar type and, for the first time, using only kinematics. Our results allow future studies to model their chemo-dynamics and evolution, thus constraining some of the most influential processes that shaped the MW. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: Submitted to Astronomy and Astrophysics. Online material available upon acceptance

Journal ref: A&A 646, A99 (2021)

Exploring the pertubed Milky Way disk and the substructures of the outer disk

Authors: Y. Xu, C. Liu, H. Tian, H. J. Newberg, C. F. P. Laporte, B. Zhang, H. F. Wang, J. Li, L. C. Deng

Abstract: The recent discovery of a spiral feature in $Z-V_Z$ phase plane in the solar neighborhood implies that the Galactic disk has been remarkably affected by a dwarf galaxy passing through it some hundreds of millions of years ago. Using 429,500 LAMOST K giants stars, we show that the spiral feature exits not only in the solar vicinity; it also extends to about 15 kpc from the Galactic center, and then… ▽ More The recent discovery of a spiral feature in $Z-V_Z$ phase plane in the solar neighborhood implies that the Galactic disk has been remarkably affected by a dwarf galaxy passing through it some hundreds of millions of years ago. Using 429,500 LAMOST K giants stars, we show that the spiral feature exits not only in the solar vicinity; it also extends to about 15 kpc from the Galactic center, and then disappears beyond this radius. Moreover, we find that when the spiral features in a plot of $V_φ$ as a function of position in the $Z-V_Z$ plane, at various Galactocentric radii, are re-mapped to $R-Z$ plane, the spiral can explain well the observed asymmetric velocity substructures. This is evidence that the phase spiral features are the same as the bulk motions found in previous as well as this work. Test-particle simulations and N-body simulations show that an encounter with a dwarf galaxy a few hundred million years ago will induce a perturbation in the Galactic disk. In addition, we find that the last impact of Sgr dSph can also contribute to the flare. As a consequence of the encounter, the distribution function of disk stars at a large range of radii is imprinted by the gravitational perturbation. △ Less

Submitted 27 October, 2020; originally announced October 2020.

Comments: 31pages, 23figures, ApJ accepted

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Authors: The Cherenkov Telescope Array Consortium, :, H. Abdalla, H. Abe, F. Acero, A. Acharyya, R. Adam, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves B, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, A. Araudo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, Y. Ascasíbar, M. Ashley , et al. (474 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle… ▽ More The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology. △ Less

Submitted 26 February, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 71 pages (including affiliations and references), 13 figures, 6 tables. Accepted in JCAP; matches published version. Corresponding authors: Jonathan Biteau, Julien Lefaucheur, Humberto Martinez-Huerta, Manuel Meyer, Santiago Pita, Ievgen Vovk

Journal ref: JCAP 02 (2021) 048

Quantifying the impact of the Large Magellanic Cloud on the structure of the Milky Way's dark matter halo using Basis Function Expansions

Authors: Nicolas Garavito-Camargo, Gurtina Besla, Chervin F. P. Laporte, Adrian M. Price-Whelan, Emily C. Cunningham, Kathryn V. Johnston, Martin D. Weinberg, Facundo A. Gomez

Abstract: Indications of disequilibrium throughout the Milky Way (MW) highlight the need for compact,flexible, non-parametric descriptions of phase--space distributions of galaxies. We present a new representation of the current Dark Matter (DM) distribution and potential derived from N-body simulations of the Milky Way and Large Magellanic Cloud (LMC) system using Basis Function Expansions (BFEs). We incor… ▽ More Indications of disequilibrium throughout the Milky Way (MW) highlight the need for compact,flexible, non-parametric descriptions of phase--space distributions of galaxies. We present a new representation of the current Dark Matter (DM) distribution and potential derived from N-body simulations of the Milky Way and Large Magellanic Cloud (LMC) system using Basis Function Expansions (BFEs). We incorporate methods to maximize the physical signal in the representation. As a result, the simulations of $10^8$ DM particles representing the MW--LMC system can be described by 354 coefficients. We find that the LMC induces asymmetric perturbations (odd l, m) to the MW's halo, which are not well-described by oblate, prolate, or triaxial halos. Furthermore, the energy in high-order even modes (l,m $\geq$ 2) is similar to average triaxial halos found in cosmological simulations. As such, the response of the MW's halo to the LMC must be accounted for in order to recover the imprints of its assembly history. The LMC causes the outer halo ($\geq$ 30 kpc) to shift from the disk center of mass (COM) by $\sim$15-25 kpc at present day, manifesting as a dipole in the BFE and in the radial velocities of halo stars. The shift depends on the LMC's infall mass, the distortion of the LMC's halo and the MW halo response. Within 30 kpc, halo tracers are expected to orbit the COM of the MW's disk, regardless of LMC infall mass. The LMC's halo is also distorted by MW tides, we discuss the implications for its mass loss and the subsequent effects on current Magellanic satellites. △ Less

Submitted 4 June, 2021; v1 submitted 2 October, 2020; originally announced October 2020.

Comments: 35 pages, 25 figures, Accepted for publication in ApJ. 3d movie can be found here: https://vimeo.com/546207117

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Angüner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascasíbar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

The bar resonances and low angular momentum moving groups in the Galaxy revealed by stellar ages

Authors: Chervin F. P. Laporte, Benoit Famaey, Giacomo Monari, Vanessa Hill, Christopher Wegg, Ortwin Gerhard

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

Submitted 24 June, 2020; originally announced June 2020.

Comments: 5 pages, 3 figures, comments welcome, submitted to A&A

Journal ref: A&A 643, L3 (2020)

Quantifying the Stellar Halo's Response to the LMC's Infall with Spherical Harmonics

Authors: Emily C. Cunningham, Nicolas Garavito-Camargo, Alis J. Deason, Kathryn V. Johnston, Denis Erkal, Chervin F. P. Laporte, Gurtina Besla, Rodrigo Luger, Robyn E. Sanderson

Abstract: The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter G… ▽ More The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall. △ Less

Submitted 15 June, 2020; originally announced June 2020.

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

Fluctuations in galactic bar parameters due to bar-spiral interaction

Authors: T. Hilmi, I. Minchev, T. Buck, M. Martig, A. C. Quillen, G. Monari, B. Famaey, R. S. de Jong, C. F. P. Laporte, J. Read, J. L. Sanders, M. Steinmetz, C. Wegg

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

Submitted 30 June, 2020; v1 submitted 11 March, 2020; originally announced March 2020.

Comments: 25 p., 19 figures, Accepted to MNRAS

Assessment of the GCT prototype's optical system implementation and other key performances for the Cherenkov Telescope Array

Authors: A. Dmytriiev, L. Dangeon, G. Fasola, H. Sol, A. Zech, J. Gironnet, O. Le Blanc, J. P. Amans, G. Buchholtz, J. L. Dournaux, F. de Frondat, D. Horville, J. M. Huet, I. Jégouzo, P. Laporte

Abstract: The Cherenkov Telescope Array (CTA) project, led by an international collaboration of institutes, aims to create the world's largest next generation observatory for Very High Energy (VHE) gamma-ray astronomy. It will be devoted to observations in a wide band of energy, from a few tens of GeV to a few hundreds of TeV with Large, Medium and Small-sized telescopes. The Small-Size Telescopes (SSTs) ar… ▽ More The Cherenkov Telescope Array (CTA) project, led by an international collaboration of institutes, aims to create the world's largest next generation observatory for Very High Energy (VHE) gamma-ray astronomy. It will be devoted to observations in a wide band of energy, from a few tens of GeV to a few hundreds of TeV with Large, Medium and Small-sized telescopes. The Small-Size Telescopes (SSTs) are dedicated to the highest energy range above a few TeV and up to 300 TeV. GCT is an imaging atmospheric Cherenkov telescope (IACT) proposed for the subarray of about 70 SSTs to be installed on the Southern site of CTA in Chile. The Observatory of Paris and the National Institute for Earth Sciences and Astronomy (INSU/CNRS) have developed the mechanical structure, mirrors (aspherical lightweight aluminium segments) and control system of the GCT. The GCT is based on a Schwarzschild- Couder (S-C) dual-mirror optical design which has the advantages, compared to the current IACTs, to offer a wide field of view (~ 9°) while decreasing the cost and volume (~ 9 m x 4 m x 6 m for ~ 11 tons) of the telescope structure, as well as the camera. The prototype (pGCT) has been installed at the Meudon's site of the Observatory of Paris and was the first S-C telescope and the first CTA prototype to record VHE events on-sky in November 2015. After three years of intensive testing, pGCT has now been commissioned. This paper is a status report on the complete GCT telescope optical system and the performance it can provide for CTA. △ Less

Submitted 25 September, 2019; originally announced September 2019.

Comments: Proceeding of the 36th International Cosmic Ray Conference (ICRC2019), July 24th - August 1st, 2019, Madison, WI, U.S.A

Journal ref: PoS(ICRC2019)663, Volume 358 - 36th International Cosmic Ray Conference (ICRC2019) - GRI - Gamma Ray Indirect

Chemo-dynamical properties of the Anticenter Stream: a surviving disc fossil from a past satellite interaction

Authors: Chervin F. P. Laporte, Vasily Belokurov, Sergey E. Koposov, Martin C. Smith, Vanessa Hill

Abstract: Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly… ▽ More Using Gaia DR2, we trace the Anticenter Stream (ACS) in various stellar populations across the sky and find that it is kinematically and spatially decoupled from the Monoceros Ring. Using stars from {\sc lamost} and {\sc segue}, we show that the ACS is systematically more metal-poor than Monoceros by $0.1$ dex with indications of a narrower metallicity spread. Furthermore, the ACS is predominantly populated of old stars ($\sim 10\,\rm{Gyr}$), whereas Monoceros has a pronounced tail of younger stars ($6-10\, \rm{Gyr}$) as revealed by their cumulative age distributions. Put togehter, all of this evidence support predictions from simulations of the interaction of the Sagittarius dwarf with the Milky Way, which argue that the Anticenter Stream (ACS) is the remains of a tidal tail of the Galaxy excited during Sgr's first pericentric passage after it crossed the virial radius, whereas Monoceros consists of the composite stellar populations excited during the more extended phases of the interaction. We suggest that the ACS can be used to constrain the Galactic potential, particularly its flattening, setting strong limits on the existence of a dark disc. Importantly, the ACS can be viewed as a stand-alone fossil of the chemical enrichment history of the Galactic disc. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 6 pages, 6 Figure, submitted to MNRAS, comments welcome

The Implications of Local Fluctuations in the Galactic Midplane for Dynamical Analysis in the Gaia Era

Authors: Angus Beane, Robyn E. Sanderson, Melissa K. Ness, Kathryn V. Johnston, Douglas Grion Filho, Mordecai-Mark Mac Low, Daniel Anglés-Alcázar, David W. Hogg, Chervin F. P. Laporte

Abstract: Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important… ▽ More Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important systematic error that is induced in the actions as a consequence of local midplane variations expected for the Milky Way. This error is difficult to model because it is non-Gaussian and bimodal, with neither mode peaking on the null value. An offset in the vertical position of the Galactic midplane of $\sim15\,\text{pc}$ for a thin disk-like orbit or $\sim 120\,\text{pc}$ for a thick disk-like orbit induces a $25\%$ systematic error in the vertical action $J_z$. In FIRE simulations of Milky Way-mass galaxies, these variations are on the order of $\sim100\,\text{pc}$ at the solar circle. From observations of the mean vertical velocity variation of $\sim5\text{--}10\,\text{km}\,\text{s}^{-1}$ with radius, we estimate that the Milky Way midplane variations are $\sim60\text{--}170\,\text{pc}$, consistent with three-dimensional dust maps. Action calculations and orbit integrations, which assume the global and local midplanes are identical, are likely to include this induced error, depending on the volume considered. Variation in the local standard of rest or distance to the Galactic center causes similar issues. The variation of the midplane must be taken into account when performing dynamical analysis across the large regions of the disk accessible to Gaia and future missions. △ Less

Submitted 22 August, 2019; v1 submitted 21 May, 2019; originally announced May 2019.

Comments: 21 pages, 10 figures, 2 tables, accepted in ApJ; comments welcome; updated to match accepted version. Code used in this work available at https://github.com/gusbeane/actions_systematic

Journal ref: ApJ 883 (2019) 103

The Multidimensional Milky Way

Authors: Robyn E. Sanderson, Jeffrey L. Carlin, Emily C. Cunningham, Nicolas Garavito-Camargo, Puragra Guhathakurta, Kathryn V. Johnston, Chervin F. P. Laporte, Ting S. Li, S. Tony Sohn

Abstract: Studying our Galaxy, the Milky Way (MW), gives us a close-up view of the interplay between cosmology, dark matter, and galaxy formation. In the next decade our understanding of the MW's dynamics, stellar populations, and structure will undergo a revolution thanks to planned and proposed astrometric, spectroscopic and photometric surveys, building on recent advances by the Gaia astrometric survey.… ▽ More Studying our Galaxy, the Milky Way (MW), gives us a close-up view of the interplay between cosmology, dark matter, and galaxy formation. In the next decade our understanding of the MW's dynamics, stellar populations, and structure will undergo a revolution thanks to planned and proposed astrometric, spectroscopic and photometric surveys, building on recent advances by the Gaia astrometric survey. Together, these new efforts will measure three-dimensional positions and velocities and numerous chemical abundances for stars to the MW's edge and well into the Local Group, leading to a complete multidimensional view of our Galaxy. Studies of the multidimensional Milky Way beyond the Gaia frontier---from the edge of the Galactic disk to the edge of our Galaxy's dark matter halo---will unlock new scientific advances across astrophysics, from constraints on dark matter to insights into galaxy formation. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Submitted as a whitepaper to the 2020 decadal survey. Typos corrected in last column of Table 1

Astrophysical Tests of Dark Matter with Maunakea Spectroscopic Explorer

Authors: Ting S. Li, Manoj Kaplinghat, Keith Bechtol, Adam S. Bolton, Jo Bovy, Timothy Carleton, Chihway Chang, Alex Drlica-Wagner, Denis Erkal, Marla Geha, Johnny P. Greco, Carl J. Grillmair, Stacy Y. Kim, Chervin F. P. Laporte, Geraint F. Lewis, Martin Makler, Yao-Yuan Mao, Jennifer L. Marshall, Alan W. McConnachie, Lina Necib, A. M. Nierenberg, Brian Nord, Andrew B. Pace, Marcel S. Pawlowski, Annika H. G. Peter , et al. (5 additional authors not shown)

Abstract: We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, an… ▽ More We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics. △ Less

Submitted 9 April, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 44 pages, 19 figures. To appear as a chapter for "The Detailed Science Case for the Maunakea Spectroscopic Explorer, 2019"

Kinematics with Gaia DR2: The Force of a Dwarf

Authors: I. Carrillo, I. Minchev, M. Steinmetz, G. Monari, C. F. P. Laporte, F. Anders, A. B. A. Queiroz, C. Chiappini, A. Khalatyan, M. Martig, P. McMillan, B. X. Santiago, K. Youakim

Abstract: We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_φ$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent wor… ▽ More We use Gaia DR2 astrometric and line-of-sight velocity information combined with two sets of distances obtained with a Bayesian inference method to study the 3D velocity distribution in the Milky Way disc. We search for variations in all Galactocentric cylindrical velocity components ($V_φ$, $V_R$ and $V_z$) with Galactic radius, azimuth, and distance from the disc mid-plane. We confirm recent work showing that bulk vertical motions in the $R\text{-}z$ plane are consistent with a combination of breathing and bending modes. In the $x\text{-}y$ plane, we show that, although the amplitudes change, the structure produced by these modes is mostly invariant as a function of distance from the plane. Comparing to two different Galactic disc models, we demonstrate that the observed patterns can drastically change in short time intervals, showing the complexity of understanding the origin of vertical perturbations. A strong radial $V_R$ gradient was identified in the inner disc, transitioning smoothly from $16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $30^\circ<φ<45^\circ$ ahead of the Sun-Galactic centre line, to $-16$ km s$^{-1}$ kpc$^{-1}$ at an azimuth of $-45^\circ<φ<-30^\circ$ lagging the solar azimuth. We use a simulation with no significant recent mergers to show that exactly the opposite trend is expected from a barred potential, but overestimated distances can flip this trend to match the data. Alternatively, using an $N$-body simulation of the Sagittarius dwarf-Milky Way interaction, we demonstrate that a major recent perturbation is necessary to reproduce the observations. Such an impact may have strongly perturbed the existing bar or even triggered its formation in the last $1\text{-}2$ Gyr. △ Less

Submitted 23 August, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: 17 p., 13 fig., accepted MNRAS

Hunting for the Dark Matter Wake Induced by the Large Magellanic Cloud

Authors: Nicolas Garavito-Camargo, Gurtina Besla, Chervin F. P Laporte, Kathryn V. Johnston, Facundo A. Gómez, Laura L. Watkins

Abstract: Satellite galaxies are predicted to generate gravitational density wakes as they orbit within the dark matter (DM) halos of their hosts, causing their orbits to decay over time. The recent infall of the Milky Way's (MW) most massive satellite galaxy, the Large Magellanic Cloud (LMC), affords us the unique opportunity to study this process in action. In this work, we present high-resolution (… ▽ More Satellite galaxies are predicted to generate gravitational density wakes as they orbit within the dark matter (DM) halos of their hosts, causing their orbits to decay over time. The recent infall of the Milky Way's (MW) most massive satellite galaxy, the Large Magellanic Cloud (LMC), affords us the unique opportunity to study this process in action. In this work, we present high-resolution ($m_{dm} = 4 \times 10^4 M_{\odot}$ ) N-body simulations of the MW-LMC interaction over the past 2 Gyr. We quantify the impact of the LMC's passage on the density and kinematics of the MW's DM halo and the observability of these structures in the MW's stellar halo. The LMC is found to generate pronounced Local and Global wakes in both the DM and stellar halos, leads to both local overdensities and distinct kinematic patterns that should be observable with ongoing and future surveys. Specifically, the Global Wake will result in redshifted radial velocities of stars in the North and blueshifts in the South, at distances larger than 45 kpc. The Local Wake traces the orbital path of the LMC through the halo (50-200 kpc), resulting in a stellar overdensity with a distinct, tangential kinematic pattern that persists to the present day. The detection of the MW's halo response will constrain: the infall mass of the LMC and its orbital trajectory, the mass of the MW, and it may inform us about the nature of the dark matter particle itself. △ Less

Submitted 11 September, 2019; v1 submitted 13 February, 2019; originally announced February 2019.

Comments: 41 pages, 27 figures. Accepted to ApJ. Some terminology was changed. High-resolution images and figures can be found at https://bit.ly/2S25YzC

The total mass of the Large Magellanic Cloud from its perturbation on the Orphan stream

Authors: D. Erkal, V. Belokurov, C. F. P. Laporte, S. E. Koposov, T. S. Li, C. J. Grillmair, N. Kallivayalil, A. M. Price-Whelan, N. W. Evans, K. Hawkins, D. Hendel, C. Mateu, J. F. Navarro, A. del Pino, C. T. Slater, S. T. Sohn

Abstract: In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potentia… ▽ More In a companion paper by Koposov et al., RR Lyrae from \textit{Gaia} Data Release 2 are used to demonstrate that stars in the Orphan stream have velocity vectors significantly misaligned with the stream track, suggesting that it has received a large gravitational perturbation from a satellite of the Milky Way. We argue that such a mismatch cannot arise due to any realistic static Milky Way potential and then explore the perturbative effects of the Large Magellanic Cloud (LMC). We find that the LMC can produce precisely the observed motion-track mismatch and we therefore use the Orphan stream to measure the mass of the Cloud. We simultaneously fit the Milky Way and LMC potentials and infer that a total LMC mass of $1.38^{+0.27}_{-0.24} \times10^{11}\,\rm{M_\odot}$ is required to bend the Orphan Stream, showing for the first time that the LMC has a large and measurable effect on structures orbiting the Milky Way. This has far-reaching consequences for any technique which assumes that tracers are orbiting a static Milky Way. Furthermore, we measure the Milky Way mass within 50 kpc to be $3.80^{+0.14}_{-0.11}\times10^{11} M_\odot$. Finally, we use these results to predict that, due to the reflex motion of the Milky Way in response to the LMC, the outskirts of the Milky Way's stellar halo should exhibit a bulk, upwards motion. △ Less

Submitted 14 May, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 17 pages, 11 figures. Updated to version accepted to MNRAS after minor revision

Piercing the Milky Way: an all-sky view of the Orphan Stream

Authors: S. E. Koposov, V. Belokurov, T. S. Li, C. Mateu, D. Erkal, C. J. Grillmair, D. Hendel, A. M. Price-Whelan, C. F. P. Laporte, K. Hawkins, S. T. Sohn, A. del Pino, N. W. Evans, C. T. Slater, N. Kallivayalil, J. F. Navarro

Abstract: We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of th… ▽ More We use astrometry, broad-band photometry and variability information from the Data Release 2 of ESA's Gaia mission (GDR2) to identify members of the Orphan Stream (OS) across the whole sky. The stream is traced above and below the celestial equator and in both Galactic hemispheres, thus increasing its visible length to ~ 210 degrees equivalent to ~150 kpc in physical extent. Taking advantage of the large number of RR Lyrae stars in the OS, we extract accurate distances and proper motions across the entire stretch of the tidal debris studied. As delineated by the GDR2 RR Lyrae, the stream exhibits two prominent twists in its shape on the sky which are accompanied by changes in the tangential motion. We complement the RR Lyrae maps with those created using GDR2 Red Giants and the DECam Legacy Survey Main Sequence Turn-Off stars. The behavior of the OS track on the sky is consistent across all three tracers employed. We detect a strong non-zero motion in the across-stream direction for a substantial portion of the stream. Such a misalignment between the debris track and the streaming velocity cannot be reproduced in a static gravitational potential and signals an interaction with a massive perturber. △ Less

Submitted 11 February, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: 18 pages; accepted to MNRAS

A-type stars in the Canada-France Imaging Survey II. Tracing the height of the disc at large distances with Blue Stragglers

Authors: Guillaume F. Thomas, Chervin F. P. Laporte, Alan W. McConnachie, Benoit Famaey, Rodrigo Ibata, Nicolas F. Martin, Else Starkenburg, Raymond Carlberg, Khyati Malhan, Kim Venn

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

Submitted 4 December, 2018; v1 submitted 22 October, 2018; originally announced October 2018.

Comments: 9 pages, 6 figures, Accepted to MNRAS

Footprints of the Sagittarius dwarf galaxy in the $Gaia$ data set

Authors: Chervin F. P. Laporte, Ivan Minchev, Kathryn V. Johnston, Facundo A. Gómez

Abstract: We analyse an N-body simulation of the interaction of the Milky Way (MW) with a Sagittarius-like dSph (Sgr), looking for signatures which may be attributed to its orbital history in the phase space volume around the Sun in light of $Gaia$ DR2 discoveries. The repeated impacts of Sgr excite coupled vertical and radial oscillations in the disc which qualitatively, and to a large degree quantitativel… ▽ More We analyse an N-body simulation of the interaction of the Milky Way (MW) with a Sagittarius-like dSph (Sgr), looking for signatures which may be attributed to its orbital history in the phase space volume around the Sun in light of $Gaia$ DR2 discoveries. The repeated impacts of Sgr excite coupled vertical and radial oscillations in the disc which qualitatively, and to a large degree quantitatively are able to reproduce many features in the 6D $Gaia$ DR2 samples, from the median $V_{R},V_φ, V_{z} $ velocity maps to the local $δρ(v_{z},z)$ phase-space spiral which is a manifestation of the global disc response to coupled oscillations within a given volume. The patterns in the large-scale velocity field are well described by tightly wound spirals and vertical corrugations excited from Sgr's impacts. We show that the last pericentric passage of Sgr resets the formation of the local present-day $δρ(v_z,z)$ spiral and situate its formation around 500-800 Myr. As expected $δρ(vz,z)$ grows in size and decreases in woundedness as a function of radius in both the $Gaia$ DR2 data and simulations. This is the first N-body model able to explain so many of the features in the data on different scales. We demonstrate how to use the full extent of the Galactic disc to date perturbations dating from Myr to Gyr, probe the underlying potential and constrain the mass-loss history of Sgr. $δρ(vz,z)$ looks the same in all stellar populations age bins down to the youngest ages which rules out a bar buckling origin. △ Less

Submitted 22 February, 2019; v1 submitted 1 August, 2018; originally announced August 2018.

Comments: 21 pages, accepted in MNRAS, 20 Figures, including previously unreported new observational results (see Fig. 6 Antoja+18 spiral in stellar overdensity, Fig. 13 on $vφ-R$ ridges + new features, Fig. 14 phase-space spiral as a function of Galactocentric radius, Fig. 17 word of caution on isochrone dating of perturbations) and extended discussion of Sgr impact models incl. further predictions