Showing 1–26 of 26 results for author: Garavito-Camargo, N

All-Sky Kinematics of the Distant Halo: The Reflex Response to the LMC

Authors: Vedant Chandra, Rohan P. Naidu, Charlie Conroy, Nicolas Garavito-Camargo, Chervin Laporte, Ana Bonaca, Phillip A. Cargile, Emily Cunningham, Jiwon Jesse Han, Benjamin D. Johnson, Hans-Walter Rix, Yuan-Sen Ting, Turner Woody, Dennis Zaritsky

Abstract: The infall of the Large Magellanic Cloud (LMC) is predicted to displace the inner Milky Way (MW), imprinting an apparent 'reflex motion' on the observed velocities of distant halo stars. We construct the largest all-sky spectroscopic dataset of luminous red giant stars from $50-160$ kpc, including a new survey of the southern celestial hemisphere. We fit the full 6D kinematics of our data to measu… ▽ More The infall of the Large Magellanic Cloud (LMC) is predicted to displace the inner Milky Way (MW), imprinting an apparent 'reflex motion' on the observed velocities of distant halo stars. We construct the largest all-sky spectroscopic dataset of luminous red giant stars from $50-160$ kpc, including a new survey of the southern celestial hemisphere. We fit the full 6D kinematics of our data to measure the amplitude and direction of the inner MW's motion towards the outer halo. The observed velocity grows with distance such that, relative to halo stars at $100$ kpc, the inner MW is lurching at $\approx 40$ km s$^{-1}$ towards a recent location along the LMC's past orbit. Our measurements align with N-body simulations of the halo's response to a $1.8 \times 10^{11} M_\odot$ LMC on first infall, suggesting that the LMC is at least 15% as massive as the MW. Our findings highlight the dramatic disequilibrium of the MW outskirts, and will enable more accurate measurements of the total mass of our Galaxy. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 25 pages, 15 figures. Submitted to ApJ

Generating synthetic star catalogs from simulated data for next-gen observatories with py-ananke

Authors: Adrien C. R. Thob, Robyn E. Sanderson, Andrew P. Eden, Farnik Nikakhtar, Nondh Panithanpaisal, Nicolás Garavito-Camargo, Sanjib Sharma

Abstract: We find ourselves on the brink of an exciting era in observational astrophysics, driven by groundbreaking facilities like JWST, Euclid, Rubin, Roman, SKA, or ELT. Simultaneously, computational astrophysics has shown significant strides, yielding highly realistic galaxy formation simulations, thanks to both hardware and software enhancements. Bridging the gap between simulations and observations ha… ▽ More We find ourselves on the brink of an exciting era in observational astrophysics, driven by groundbreaking facilities like JWST, Euclid, Rubin, Roman, SKA, or ELT. Simultaneously, computational astrophysics has shown significant strides, yielding highly realistic galaxy formation simulations, thanks to both hardware and software enhancements. Bridging the gap between simulations and observations has become paramount for meaningful comparisons. We introduce py-ananke, a Python pipeline designed to generate synthetic resolved stellar surveys from cosmological simulations, adaptable to various instruments. Building upon its predecessor, ananke by Sanderson et al. 2020 (arXiv:1806.10564), which produced Gaia DR2 mock star surveys, the py-ananke package offers a user-friendly "plug & play" experience. The pipeline employs cutting-edge phase-space density estimation and initial mass function sampling to convert particle data into synthetic stars, while interpolating pre-computed stellar isochrone tracks for photometry. Additionally, it includes modules for estimating interstellar reddening, dust-induced extinctions, and for quantifying errors through dedicated modeling approaches. py-ananke promises to serve as a vital bridge between computational astrophysics and observational astronomy, facilitating preparations and making scientific predictions for the next generation of telescopes. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: Submitted to the Journal of Open Source Software. 12 pages (6 pages of bibliography) and 1 figure. Software repository at https://github.com/athob/py-ananke

On the co-rotation of Milky Way satellites: LMC-mass satellites induce apparent motions in outer halo tracers

Authors: Nicolas Garavito-Camargo, Adrian M. Price-Whelan, Jenna Samuel, Emily C. Cunningham, Ekta Patel, Andrew Wetzel, Kathryn V. Johnston, Arpit Arora, Robyn E. Sanderson, Lehman Garrison, Danny Horta

Abstract: Understanding the physical mechanism behind the formation of a co-rotating thin plane of satellite galaxies, like the one observed around the Milky Way (MW), has been challenging. The perturbations induced by a massive satellite galaxy, like the Large Magellanic Cloud (LMC) provide valuable insight into this problem. The LMC induces an apparent co-rotating motion in the outer halo by displacing th… ▽ More Understanding the physical mechanism behind the formation of a co-rotating thin plane of satellite galaxies, like the one observed around the Milky Way (MW), has been challenging. The perturbations induced by a massive satellite galaxy, like the Large Magellanic Cloud (LMC) provide valuable insight into this problem. The LMC induces an apparent co-rotating motion in the outer halo by displacing the inner regions of the halo with respect to the outer halo. Using the Latte suite of FIRE-2 cosmological simulations of MW-mass galaxies, we confirm that the apparent motion of the outer halo induced by the infall of a massive satellite changes the observed distribution of orbital poles of outer-halo tracers, including satellites. We quantify the changes in the distribution of orbital poles using the two-point angular correlation function and find that all satellites induce changes. However, the most massive satellites with pericentric passages between 30-100kpc induce the largest changes. The best LMC-like satellite analog shows the largest change in orbital pole distribution. The dispersion of orbital poles decreases by 20° during the first two pericentric passages. Even when excluding the satellites brought in with the LMC-like satellite, there is clustering of orbital poles. These results suggest that in the MW, the recent pericentric passage of the LMC should have changed the observed distribution of orbital poles of all other satellites. Therefore, studies of kinematically-coherent planes of satellites that seek to place the MW in a cosmological context should account for the existence of a massive satellite like the LMC. △ Less

Submitted 19 November, 2023; originally announced November 2023.

Comments: 20 pages, 10 figures. ApJ submitted, Comments are welcome

Galactoseismology in cosmological simulations: Vertical perturbations by dark matter, satellite galaxies and gas

Authors: B. García-Conde, T. Antoja, S. Roca-Fàbrega, F. Gómez, P. Ramos, N. Garavito-Camargo, MA. Gómez-Flechoso

Abstract: Only recently, complex models that include the global dynamics from dwarf satellite galaxies, dark matter halo structure, gas infalls, and stellar disk in a cosmological context became available to study the dynamics of disk galaxies such as the Milky Way (MW). We use a MW model from a high-resolution hydrodynamical cosmological simulation named GARROTXA to establish the relationship between the v… ▽ More Only recently, complex models that include the global dynamics from dwarf satellite galaxies, dark matter halo structure, gas infalls, and stellar disk in a cosmological context became available to study the dynamics of disk galaxies such as the Milky Way (MW). We use a MW model from a high-resolution hydrodynamical cosmological simulation named GARROTXA to establish the relationship between the vertical disturbances seen in its galactic disk and multiple perturbations, from the dark matter halo, satellites and gas. We calculate the bending modes in the galactic disk in the last 6 Gyr of evolution. To quantify the impact of dark matter and gas we compute the vertical acceleration exerted by these components onto the disk and compare them with the bending behavior with Fourier analysis. We find complex bending patterns at different radii and times, such as an inner retrograde mode with high frequency, as well as an outer slower retrograde mode excited at different times. The amplitudes of these bending modes are highest during the early stages of the thin disk formation and reach up to 8.5 km s-1 in the late disk evolution. We find that the infall of satellite galaxies leads to a tilt of the disk, and produces anisotropic gas accretion with subsequent star formation events, and supernovae, creating significant vertical accelerations onto the disk plane. The misalignment between the disk and the inner stellar/dark matter triaxial structure, formed during the ancient assembly of the galaxy, creates a strong vertical acceleration on the stars. We conclude that several agents trigger the bending of the stellar disk and its phase spirals in this simulation, including satellite galaxies, dark sub-halos, misaligned gaseous structures, and the inner dark matter profile, which coexist and influence each other, making it challenging to establish direct causality. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Dark matter distribution in Milky Way-analog galaxies

Authors: Natanael G. de Isídio, K. Menéndez-Delmestre, T. S. Gonçalves, M. Grossi, D. C. Rodrigues, N. Garavito-Camargo, A. Araujo-Carvalho, P. P. B. Beaklini, Y. Cavalcante-Coelho, A. Cortesi, L. H. Queiroga-Nuñez, T. Randriamampandry

Abstract: Our current understanding of how dark matter (DM) is distributed within the Milky Way halo, particularly in the solar neighborhood, is based on either careful studies of the local stellar orbits, model assumptions on the global shape of the MW halo, or from direct acceleration measurements. In this work, we undertake a study of external galaxies, with the intent of providing insight to the DM dist… ▽ More Our current understanding of how dark matter (DM) is distributed within the Milky Way halo, particularly in the solar neighborhood, is based on either careful studies of the local stellar orbits, model assumptions on the global shape of the MW halo, or from direct acceleration measurements. In this work, we undertake a study of external galaxies, with the intent of providing insight to the DM distribution in MW-analog galaxies. For this, we carefully select a sample of galaxies similar to the MW, based on maximum atomic hydrogen (HI) rotational velocity ($v_{\rm max,HI}$= 200-280 km s$^{-1}$) and morphological type (Sab-Sbc) criteria. With a need for deep, highly-resolved HI, our resulting sample is composed of 5 galaxies from the VIVA and THINGS surveys. To perform our baryonic analysis, we use deep Spitzer mid-IR images at 3.6 and 4.5 $μ$m from the S$^4$G survey. Based on the dynamical three-dimensional modeling software $^{3D}$Barolo, we construct rotation curves (RCs) and derive the gas and stellar contributions from the galaxy's gaseous- and stellar-disks mass surface density profiles. Through a careful decomposition of their RCs into their baryonic (stars, gas) and DM components, we isolate the DM contribution by using an MCMC-based approach. Based on the Sun's location and the MW's R$_{25}$, we define the corresponding location of the solar neighborhood in these systems. We put forward a window for the DM density ($ρ_{dm}$ = 0.21-0.55 GeV cm$^{-3}$) at these galactocentric distances in our MW analog sample, consistent with the values found for the MW's local DM density, based on more traditional approaches found in the literature. △ Less

Submitted 4 June, 2024; v1 submitted 20 October, 2023; originally announced October 2023.

Comments: 18 pages, 5 figures, 3 tables, Accepted to ApJ

LMC-driven anisotropic boosts in stream--subhalo interactions

Authors: Arpit Arora, Nicolás Garavito-Camargo, Robyn E. Sanderson, Emily C. Cunningham, Andrew Wetzel, Nondh Panithanpaisal, Megan Barry

Abstract: Dark Matter (DM) subhalos are predicted to perturb stellar streams; stream morphologies and dynamics can constrain the mass distribution of subhalos. Using FIRE-2 simulations of Milky Way-mass galaxies, we show that presence of a Large Magellanic Cloud (LMC)--analog significantly changes stream-subhalo encounter rates. Three key factors drive these changes. First, the LMC--analog brings in many su… ▽ More Dark Matter (DM) subhalos are predicted to perturb stellar streams; stream morphologies and dynamics can constrain the mass distribution of subhalos. Using FIRE-2 simulations of Milky Way-mass galaxies, we show that presence of a Large Magellanic Cloud (LMC)--analog significantly changes stream-subhalo encounter rates. Three key factors drive these changes. First, the LMC--analog brings in many subhalos, increasing encounter rates for streams near the massive satellite by up to 20--40%. Second, the LMC--analog displaces the host from its center-of-mass (inducing reflex motion), causing a north-south asymmetry in the density and radial velocity distribution of subhalos. This asymmetry results in encounter rates varying by 50--70% across the sky at the same distance. Finally, the LMC--mass satellite induces a density wake in the host's DM halo, further boosting the encounter rates near the LMC--analog. We also explore the influence of stream orbital properties, finding a 50% increase in encounters for streams moving retrograde to the LMC--analog's orbit in the opposite hemisphere. The dependence of encounter rates on stream location and orbit has important implications for where to search for new streams with spurs and gaps in the Milky Way. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: 26 pages, 15 figures, submitted to APJ

The proto-galaxy of Milky Way-mass haloes in the FIRE simulations

Authors: Danny Horta, Emily C. Cunningham, Robyn Sanderson, Kathryn V. Johnston, Alis Deason, Andrew Wetzel, Fiona McCluskey, Nicolás Garavito-Camargo, Lina Necib, Claude-André Faucher-Giguère, Arpit Arora, Pratik J. Gandhi

Abstract: Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e., proto-Galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological… ▽ More Observational studies are finding stars believed to be relics of the earliest stages of hierarchical mass assembly of the Milky Way (i.e., proto-Galaxy). In this work, we contextualize these findings by studying the masses, ages, spatial distributions, morphology, kinematics, and chemical compositions of proto-galaxy populations from the 13 Milky Way (MW)-mass galaxies from the FIRE-2 cosmological zoom-in simulations. Our findings indicate that proto-Milky Way populations: i) can have a stellar mass range between $1\times10^{8}<\mathrm{M}_{\star}<2\times10^{10}[\mathrm{M}_{\odot}]$, a virial mass range between $3\times10^{10}<\mathrm{M}_{\star}<6\times10^{11}[\mathrm{M}_{\odot}]$, and be as young as $8 \lesssim \mathrm{Age} \lesssim 12.8$ [Gyr] ($1\lesssim z \lesssim 6$); ii) are predominantly centrally concentrated, with $\sim50\%$ of the stars contained within $5-10$ kpc; iii) on average show weak but systematic net rotation in the plane of the host's disc at $z=0$ (i.e., $0.25\lesssim\langleκ/κ_{\mathrm{disc}}\rangle\lesssim0.8$); iv) present [$α$/Fe]-[Fe/H] compositions that overlap with the metal-poor tail of the host's old disc; v) tend to assemble slightly earlier in Local Group-like environments than in systems in isolation. Interestingly, we find that ~60% of the proto-Milky Way galaxies are comprised by 1 dominant system ($1/5\lesssim$M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim4/5$) and 4-5 lower mass systems (M$_{\star}$/M$_{\star,\mathrm{proto-Milky Way}}$$\lesssim1/10$); the other ~40% are comprised by 2 dominant systems and 3-4 lower mass systems. These massive/dominant proto-Milky Way fragments can be distinguished from the lower mass ones in chemical-kinematic samples, but appear (qualitatively) indistinguishable from one another. Our results could help observational studies disentangle if the Milky Way formed from one or two dominant systems. △ Less

Submitted 11 December, 2023; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: Accepted in MNRAS. 16 pages, 12 figures and 3 tables

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)

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

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

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

Submitted 20 June, 2023; originally announced June 2023.

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

The impact of the Large Magellanic Cloud on dark matter direct detection signals

Authors: Adam Smith-Orlik, Nima Ronaghi, Nassim Bozorgnia, Marius Cautun, Azadeh Fattahi, Gurtina Besla, Carlos S. Frenk, Nicolás Garavito-Camargo, Facundo A. Gómez, Robert J. J. Grand, Federico Marinacci, Annika H. G. Peter

Abstract: We study the effect of the Large Magellanic Cloud (LMC) on the dark matter (DM) distribution in the Solar neighborhood, utilizing the Auriga magneto-hydrodynamical simulations of Milky Way (MW) analogues that have an LMC-like system. We extract the local DM velocity distribution at different times during the orbit of the LMC around the MW in the simulations. As found in previous idealized simulati… ▽ More We study the effect of the Large Magellanic Cloud (LMC) on the dark matter (DM) distribution in the Solar neighborhood, utilizing the Auriga magneto-hydrodynamical simulations of Milky Way (MW) analogues that have an LMC-like system. We extract the local DM velocity distribution at different times during the orbit of the LMC around the MW in the simulations. As found in previous idealized simulations of the MW-LMC system, we find that the DM particles in the Solar neighborhood originating from the LMC analogue dominate the high speed tail of the local DM speed distribution. Furthermore, the native DM particles of the MW in the Solar region are boosted to higher speeds as a result of a response to the LMC's motion. We simulate the signals expected in near future xenon, germanium, and silicon direct detection experiments, considering DM interactions with target nuclei or electrons. We find that the presence of the LMC causes a considerable shift in the expected direct detection exclusion limits towards smaller cross sections and DM masses, with the effect being more prominent for low mass DM. Hence, our study shows, for the first time, that the LMC's influence on the local DM distribution is significant even in fully cosmological MW analogues. △ Less

Submitted 25 October, 2023; v1 submitted 8 February, 2023; originally announced February 2023.

Comments: 35 pages, 14 figures, 3 tables, JCAP accepted version

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

On the stability of tidal streams in action space

Authors: Arpit Arora, Robyn E. Sanderson, Nondh Panithanpaisal, Emily C. Cunningham, Andrew Wetzel, Nicolás Garavito-Camargo

Abstract: In the Gaia era it is increasingly apparent that traditional static, parameterized models are insufficient to describe the mass distribution of our complex, dynamically evolving Milky Way (MW). In this work, we compare different time-evolving and time-independent representations of the gravitational potentials of simulated MW-mass galaxies from the FIRE-2 suite of cosmological baryonic simulations… ▽ More In the Gaia era it is increasingly apparent that traditional static, parameterized models are insufficient to describe the mass distribution of our complex, dynamically evolving Milky Way (MW). In this work, we compare different time-evolving and time-independent representations of the gravitational potentials of simulated MW-mass galaxies from the FIRE-2 suite of cosmological baryonic simulations. Using these potentials, we calculate actions for star particles in tidal streams around three galaxies with varying merger histories at each snapshot from 7 Gyr ago to the present day. We determine the action-space coherence preserved by each model using the Kullback-Leibler Divergence to gauge the degree of clustering in actions and the relative stability of the clusters over time. We find that all models produce a clustered action space for simulations with no significant mergers. However, a massive (mass ratio prior to infall more similar than 1:8) interacting galaxy not present in the model will result in mischaracterized orbits for stars most affected by the interaction. The locations of the action space clusters (i.e. the orbits of the stream stars) are only preserved by the time-evolving model, while the time-independent models can lose significant amounts of information as soon as 0.5--1 Gyr ago, even if the system does not undergo a significant merger. Our results imply that reverse-integration of stream orbits in the MW using a fixed potential is likely to give incorrect results if integrated longer than 0.5 Gyr into the past. △ Less

Submitted 26 October, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Journal ref: The Astrophysical Journal, vol. 939, no. 1, 2022, p. 2

Implications of the Milky Way travel velocity for dynamical mass estimates of the Local Group

Authors: Katie Chamberlain, Adrian M. Price-Whelan, Gurtina Besla, Emily C. Cunningham, Nicolás Garavito-Camargo, Jorge Peñarrubia, Michael S. Petersen

Abstract: The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the "Timing Argument," which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of… ▽ More The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the "Timing Argument," which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history including the recent pericentric passage of the LMC, and recent work has found that the MW disk is moving with a lower bound "travel velocity" of $\sim 32~{\rm km}~{\rm s}^{-1}$ with respect to the outer stellar halo. Previous Timing Argument measurements attempt to account for this non-equilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10--12\% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume. △ Less

Submitted 30 December, 2022; v1 submitted 14 April, 2022; originally announced April 2022.

Comments: 14 pages, 3 figures; Published in ApJ; updated to match published version

Journal ref: ApJ 942 18 (2023)

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

All-Sky Dynamical Response of the Galactic Halo to the Large Magellanic Cloud

Authors: Charlie Conroy, Rohan P. Naidu, Nicolas Garavito-Camargo, Gurtina Besla, Dennis Zaritsky, Ana Bonaca, Benjamin D. Johnson

Abstract: Gravitational interactions between the Large Magellanic Cloud (LMC) and the stellar and dark matter halo of the Milky Way are expected to give rise to disequilibrium phenomena in the outer Milky Way. A local wake is predicted to trail the orbit of the LMC, while a large-scale over-density is predicted to exist across a large area of the northern Galactic hemisphere. Here we present the detection o… ▽ More Gravitational interactions between the Large Magellanic Cloud (LMC) and the stellar and dark matter halo of the Milky Way are expected to give rise to disequilibrium phenomena in the outer Milky Way. A local wake is predicted to trail the orbit of the LMC, while a large-scale over-density is predicted to exist across a large area of the northern Galactic hemisphere. Here we present the detection of both the local wake and Northern over-density (hereafter the "collective response") in an all-sky star map of the Galaxy based on 1301 stars at 60<R_gal<100 kpc. The location of the wake is in good agreement with an N-body simulation that includes the dynamical effect of the LMC on the Milky Way halo. The density contrast of the wake and collective response are both stronger in the data than in the simulation. The detection of a strong local wake is independent evidence that the Magellanic Clouds are on their first orbit around the Milky Way. The wake traces the path of the LMC, which will provide insight into the orbit of the LMC, which in turn is a sensitive probe of the mass of the LMC and the Milky Way. These data demonstrate that the outer halo is not in dynamical equilibrium, as is often assumed. The morphology and strength of the wake could be used to test the nature of dark matter and gravity. △ Less

Submitted 19 April, 2021; originally announced April 2021.

Comments: To appear in the 22 April 2021 issue of Nature

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

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

The Orbital Histories of Magellanic Satellites Using Gaia DR2 Proper Motions

Authors: Ekta Patel, Nitya Kallivayalil, Nicolas Garavito-Camargo, Gurtina Besla, Daniel R. Weisz, Roeland P. van der Marel, Michael Boylan-Kolchin, Marcel S. Pawlowski, Facundo A. Gómez

Abstract: With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest mass, ultra-faint satellite galaxies in the Milky Way's (MW) halo for the first time. Many of these faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their 6-dimensional phase space information, we calculate the orbital histories of 13 ultra-faint… ▽ More With the release of Gaia DR2, it is now possible to measure the proper motions (PMs) of the lowest mass, ultra-faint satellite galaxies in the Milky Way's (MW) halo for the first time. Many of these faint satellites are posited to have been accreted as satellites of the Magellanic Clouds (MCs). Using their 6-dimensional phase space information, we calculate the orbital histories of 13 ultra-faint satellites and five classical dwarf spheroidals in a combined MW+LMC+SMC potential to determine which galaxies are dynamically associated with the MCs. These 18 galaxies are separated into four classes: i.) long-term Magellanic satellites that have been bound to the MCs for at least the last two consecutive orbits around the MCs (Carina 2, Carina 3, Horologium 1, Hydrus 1); ii.) Magellanic satellites that were recently captured by the MCs $<$ 1 Gyr ago (Reticulum 2, Phoenix 2); iii.) MW satellites that have interacted with the MCs (Sculptor 1, Tucana 3, Segue 1); and iv.) MW satellites (Aquarius 2, Canes Venatici 2, Crater 2, Draco 1, Draco 2, Hydra 2, Carina, Fornax, Ursa Minor). Results are reported for a range of MW and LMC masses. Contrary to previous work, we find no dynamical association between Carina, Fornax, and the MCs. Finally, we determine that the addition of the SMC's gravitational potential affects the longevity of satellites as members of the Magellanic system (long-term versus recently captured), but it does not change the total number of Magellanic satellites. △ Less

Submitted 1 March, 2020; v1 submitted 6 January, 2020; originally announced January 2020.

Comments: 32 pages, 11 figures. Accepted to ApJ

The highest-speed local dark matter particles come from the Large Magellanic Cloud

Authors: Gurtina Besla, Annika Peter, Nicolas Garavito-Camargo

Abstract: Using N-body simulations of the Large Magellanic Cloud (LMC's) passage through the Milky Way (MW), tailored to reproduce observed kinematic properties of both galaxies, we show that the high-speed tail of the Solar Neighborhood dark matter distribution is overwhelmingly of LMC origin. Two populations contribute at high speeds: 1) Particles that were once bound to the LMC, and 2) MW halo particles… ▽ More Using N-body simulations of the Large Magellanic Cloud (LMC's) passage through the Milky Way (MW), tailored to reproduce observed kinematic properties of both galaxies, we show that the high-speed tail of the Solar Neighborhood dark matter distribution is overwhelmingly of LMC origin. Two populations contribute at high speeds: 1) Particles that were once bound to the LMC, and 2) MW halo particles that have been accelerated owing to the response of the halo to the recent passage of the LMC. These particles reach speeds of 700-900 km/s with respect to the Earth, above the local escape speed of the MW. The high-speed particles follow trajectories similar to the Solar reflex motion, with peak velocities reached in June. For low-mass dark matter, these high-speed particles can dominate the signal in direct-detection experiments, extending the reach of the experiments to lower mass and elastic scattering cross sections even with existing data sets. Our study shows that even non-disrupted MW satellite galaxies can leave a significant dark-matter footprint in the Solar Neighborhood. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Comments: Submitted to JCAP, comments are welcome

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

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 Extremely Luminous Quasar Survey (ELQS) in the Sloan Digital Sky Survey footprint. II. The North Galactic Cap Sample

Authors: Jan-Torge Schindler, Xiaohui Fan, Ian D. McGreer, Jinyi Yang, Feige Wang, Richard Green, Nicolas Garavito-Camargo, Yun-Hsin Huang, Christine O'Donnell, Anna Patej, Ragadeepika Pucha, Jon M. Rees, Eckhart Spalding

Abstract: We present the North Galactic Cap sample of the Extremely Luminous Quasar Survey (ELQS-N), which targets quasars with $M_{1450}<-27$ at $2.8 \leq z < 5$ in an area of $\sim7600\,\rm{deg}^2$ of the Sloan Digital Sky Survey (SDSS) footprint with $90\text{\textdegree}<\rm{RA}<270\text{\textdegree}$. Based on a near-infrared/infrared \textit{JKW2} color cut, the ELQS selection efficiently uses random… ▽ More We present the North Galactic Cap sample of the Extremely Luminous Quasar Survey (ELQS-N), which targets quasars with $M_{1450}<-27$ at $2.8 \leq z < 5$ in an area of $\sim7600\,\rm{deg}^2$ of the Sloan Digital Sky Survey (SDSS) footprint with $90\text{\textdegree}<\rm{RA}<270\text{\textdegree}$. Based on a near-infrared/infrared \textit{JKW2} color cut, the ELQS selection efficiently uses random forest methods to classify quasars and to estimate photometric redshifts; this scheme overcomes some of the difficulties of pure optical quasar selection at $z\approx3$. As a result, we retain a completeness of $>70\%$ over $z\sim3.0-5.0$ at $m_{i}\lesssim17.5$, limited toward fainter magnitudes by the depth of the Two Micron All Sky Survey (2MASS). The presented quasar catalog consists of a total of 270 objects, of which 39 are newly identified in this work with spectroscopy obtained at the Vatican Advanced Technology Telescope and the MMT $6.5\,\rm{m}$ telescope. In addition to the high completeness, which allowed us to discover new quasars in the already well-surveyed SDSS North Galactic Cap, the efficiency of our selection is relatively high at $\sim79\%$. Using 120 objects of this quasar sample we are able to extend the previously measured optical quasar luminosity function (QLF) by one magnitude toward the bright end at $2.8 \leq z \leq 4.5$. A first analysis of the QLF suggests a relatively steep bright-end slope of $β\approx-4$ for this sample. This result contrasts with previous results in the same redshift range, which find a much flatter slope around $β\sim-2.5$, but agrees with recent measurements of the bright-end slope at lower and higher redshifts. Our results constrain the bright-end slope at $z=2.8-4.5$ to $β<-2.94$ with a 99\% confidence. △ Less

Submitted 6 August, 2018; v1 submitted 8 June, 2018; originally announced June 2018.

Comments: 23 pages, 13 figures, submitted to ApJ

Modelling the gas kinematics of an atypical Lyman-alpha emitting compact dwarf galaxy

Authors: Jaime E. Forero-Romero, Max Gronke, Maria Camila Remolina-Gutiérrez, Nicolás Garavito-Camargo, Mark Dijkstra

Abstract: Star-forming Compact Dwarf Galaxies (CDGs) resemble the expected pristine conditions of the first galaxies in the Universe and are the best systems to test models on primordial galaxy formation and evolution. Here we report on one of such CDGs, Tololo 1214-277, which presents a broad, single peaked, highly symmetric Ly$α$ emission line that had evaded theoretical interpretation so far. In this pap… ▽ More Star-forming Compact Dwarf Galaxies (CDGs) resemble the expected pristine conditions of the first galaxies in the Universe and are the best systems to test models on primordial galaxy formation and evolution. Here we report on one of such CDGs, Tololo 1214-277, which presents a broad, single peaked, highly symmetric Ly$α$ emission line that had evaded theoretical interpretation so far. In this paper we reproduce for the first time these line features with two different physically motivated kinematic models: an interstellar medium composed by outflowing clumps with random motions and an homogeneous gaseous sphere undergoing solid body rotation. The multiphase model requires a clump velocity dispersion of $54.3\pm 0.6$ km s$^{-1}$ with outflows of $54.3\pm 5.1$ km s$^{-1}$, while the bulk rotation velocity is constrained to be $348^{+75}_{-48}$ km s$^{-1}$. We argue that the results from the multiphase model provide a correct interpretation of the data. In that case the clump velocity dispersion implies a dynamical mass of $2\times 10^{9}$ M$_{\odot}$, ten times its baryonic mass. If future kinematic maps of Tololo 1214-277 confirm the velocities suggested by the multiphase model, it would provide additional support to expect such kinematic state in primordial galaxies, opening the opportunity to use the models and methods presented in this paper to constrain the physics of star formation and feedback in the early generation of Ly-$α$ emitting galaxies. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: 10 pages, 5 figures, 2 tables. Accepted for publication in MNRAS

The Influence of Sagittarius and the Large Magellanic Cloud on the Milky Way Galaxy

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

Abstract: We present N-body simulations of a Sagittarius-like dwarf spheroidal galaxy (Sgr) that follow its orbit about the Milky Way (MW) since its first crossing of the Galaxy's virial radius to the present day. As Sgr orbits around the MW, it excites vertical oscillations, corrugating and flaring the Galactic stellar disc. These responses can be understood by a two-phase picture in which the interaction… ▽ More We present N-body simulations of a Sagittarius-like dwarf spheroidal galaxy (Sgr) that follow its orbit about the Milky Way (MW) since its first crossing of the Galaxy's virial radius to the present day. As Sgr orbits around the MW, it excites vertical oscillations, corrugating and flaring the Galactic stellar disc. These responses can be understood by a two-phase picture in which the interaction is first dominated by torques from the wake excited by Sgr in the MW dark halo before transitioning to tides from Sgr's direct impact on the disc at late times. We show for the first time that a massive Sgr model simultaneously reproduces the locations and motions of arc-like over densities, such as the Monoceros Ring and the Triangulum Andromeda stellar clouds, that have been observed at the extremities of the disc, while also satisfying the solar neighbourhood constraints on the vertical structure and streaming motions of the disc. In additional simulations, we include the Large Magellanic Cloud (LMC) self consistently with Sgr. The LMC introduces coupling through constructive and destructive interference, but no new corrugations. In our models, the excitation of the current structure of the outer disk can be traced to interactions as far back as 6-7 Gyr ago (corresponding to $z\leq1$). Given the apparently quiescent accretion history of the MW over this timescale, this places Sgr as the main culprit behind the vertical oscillations of the disc and the last major accretion event for the Galaxy with the capacity to modulate its chemodynamical structure. △ Less

Submitted 6 October, 2017; originally announced October 2017.

Comments: 24 pages, 28 Figures (incl Appendix), submitted to MNRAS, comments welcome

Response of the Milky Way's disc to the Large Magellanic Cloud in a first infall scenario

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

Abstract: We present N-body and hydrodynamical simulations of the response of the Milky Way's baryonic disc to the presence of the Large Magellanic Cloud during a first infall scenario. For a fiducial galactic model reproducing the gross properties of the Galaxy, we explore a set of six initial conditions for the LMC of varying mass which all evolve to fit the measured constraints on its current position an… ▽ More We present N-body and hydrodynamical simulations of the response of the Milky Way's baryonic disc to the presence of the Large Magellanic Cloud during a first infall scenario. For a fiducial galactic model reproducing the gross properties of the Galaxy, we explore a set of six initial conditions for the LMC of varying mass which all evolve to fit the measured constraints on its current position and velocity with respect to the Galactic Center. We find that the LMC can produce strong disturbances - warping of the stellar and gaseous discs - in the Galaxy, without violating constraints from the phase-space distribution of stars in the Solar Neighbourhood. All models correctly reproduce the phases of the warp and its anti-symmetrical shape about the disc's mid-plane. If the warp is due to the LMC alone, then the largest mass model is favoured ($2.5\times10^{11}\,\rm{M_{\odot}}$). Still, some quantitative discrepancies remain, including deficits in height of $ΔZ=0.7 \, \rm{kpc}$ at $R=22 \,\rm{kpc}$ and $ΔZ=0.7\, \rm{kpc}$ at $R=16 \,\rm{kpc}$. This suggests that even higher infall masses for the LMC's halo are allowed by the data. A comparison with the vertical perturbations induced by a heavy Sagittarius dSph model ($10^{11}\,\rm{M_{\odot}}$) suggest that positive interference with the LMC warp is expected at $R=16 \, \rm{kpc}$. We conclude that the vertical structure of the Galactic disc beyond the Solar Neighbourhood may jointly be shaped by its most massive satellites. As such, the current structure of the MW suggests we are seeing the process of disc heating by satellite interactions in action. △ Less

Submitted 16 August, 2016; originally announced August 2016.

Comments: 13 pages, 13 figures, submitted to MNRAS, comments welcome

The impact of gas bulk rotation on the lyman-alpha line

Authors: Juan N. Garavito-Camargo, Jaime E. Forero-Romero, Mark Dijkstra

Abstract: We present results of radiative transfer calculations to measure the impact of gas bulk rotation on the morphology of the Lyman $α$ emission line in distant galaxies. We model a galaxy as a sphere with an homogeneous mixture of dust and hydrogen at a constant temperature. These spheres undergo solid-body rotation with maximum velocities in the range 0-300 \kms and neutral hydrogen optical depths i… ▽ More We present results of radiative transfer calculations to measure the impact of gas bulk rotation on the morphology of the Lyman $α$ emission line in distant galaxies. We model a galaxy as a sphere with an homogeneous mixture of dust and hydrogen at a constant temperature. These spheres undergo solid-body rotation with maximum velocities in the range 0-300 \kms and neutral hydrogen optical depths in the range $τ_{\rm H}=10^{5}-10^{7}$. We consider two types of source distributions in the sphere: central and homogeneous. Our main result is that rotation introduces a dependence of the line morphology with viewing angle and rotational velocity. Observations with a line of sight parallel to the rotation axis yield line morphologies similar to the static case. For lines of sight perpendicular to the rotation axis both the intensity at the line center and the line width increase with rotational velocity. Along the same line of sight, the line becomes single peaked at rotational velocities close to half the line width in the static case. Notably, we find that rotation does not induce any spatial anisotropy in the integrated line flux, the escape fraction or the average number of scatterings. This is because Lyman α scattering through a rotating solid-body proceeds identical as in the static case. The only difference is the doppler shift from the different regions in the sphere that move with respect to the observer. This allows us to derive an analytic approximation for the viewing-angle dependence of the emerging spectrum, as a function of rotational velocity. △ Less

Submitted 17 September, 2014; v1 submitted 11 August, 2014; originally announced August 2014.

Comments: 13 pages, 14 figures. Accepted to ApJ

Journal ref: ApJ, 795, 120, 2014