Showing 1–50 of 100 results for author: Sohn, T

HSTPROMO Internal Proper Motion Kinematics of Dwarf Spheroidal Galaxies: I. Velocity Anisotropy and Dark Matter Cusp Slope of Draco

Authors: Eduardo Vitral, Roeland P. van der Marel, Sangmo Tony Sohn, Mattia Libralato, Andrés del Pino, Laura L. Watkins, Andrea Bellini, Matthew G. Walker, Gurtina Besla, Marcel S. Pawlowski, Gary A. Mamon

Abstract: We analyze four epochs of HST imaging over 18 years for the Draco dwarf spheroidal galaxy. We measure precise proper motions (PMs) for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially-resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass-anisotropy degen… ▽ More We analyze four epochs of HST imaging over 18 years for the Draco dwarf spheroidal galaxy. We measure precise proper motions (PMs) for hundreds of stars and combine these with existing line-of-sight (LOS) velocities. This provides the first radially-resolved 3D velocity dispersion profiles for any dwarf galaxy. These constrain the intrinsic velocity anisotropy and resolve the mass-anisotropy degeneracy. We solve the Jeans equations in oblate axisymmetric geometry to infer the mass profile. We find the velocity dispersion to be radially anisotropic along the symmetry axis and tangentially anisotropic in the equatorial plane, with a globally-averaged value $\overline{β_{\mathrm B}}=-0.20^{+ 0.28}_{- 0.53}$, (where $1 - β_{\mathrm B} \equiv \langle v_{\mathrm{ tan}}^2 \rangle / \langle v_{\mathrm{ rad}}^2 \rangle$ in 3D). The logarithmic dark matter (DM) density slope over the observed radial range, $Γ_{\mathrm{ dark}}$, is $-0.83^{+ 0.32}_{- 0.37}$, consistent with the inner cusp predicted in $Λ$CDM cosmology. As expected given Draco's low mass and ancient star formation history, it does not appear to have been dissolved by baryonic processes. We rule out cores larger than 487, 717, 942 pc at respective 1-, 2-, 3-$σ$ confidence, thus imposing important constraints on the self-interacting DM cross-section. Spherical models yield biased estimates for both the velocity anisotropy and the inferred slope. The circular velocity at our outermost data point (900 pc) is $24.19^{+ 6.31}_{- 2.97} \ \mathrm{km~s^{-1}}s$. We infer a dynamical distance of $75.37^{+ 4.73}_{- 4.00}$ kpc, and show that Draco has a modest LOS rotation, with $\left<v / σ\right> = 0.22 \pm 0.09$. Our results provide a new stringent test of the so-called `cusp-core' problem that can be readily extended to other dwarfs. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 35 pages, 18 figures, 3 Tables. Accepted for publication in ApJ. Journal version has better readability. Data is available at Zenodo: https://zenodo.org/records/11111113

Journal ref: The Astrophysical Journal, 970:1 (26pp), 2024 July 20

Satellite group infall into the Milky Way: exploring the Crater-Leo case with new HST proper motions

Authors: Mariana P. Júlio, Marcel S. Pawlowski, Sangmo Tony Sohn, Salvatore Taibi, Roeland P. van der Marel, Stacy S. McGaugh

Abstract: Context. Within $Λ$ Cold Dark Matter ($Λ$CDM) simulations, Milky Way-like galaxies accrete some of their satellite galaxies in groups of 3-5 members rather than individually. It was also suggested that this might be the reason behind the origin of satellite planes. Objects accreted in groups are expected to share similar specific total energy and angular momentum, and also identical orbital planes… ▽ More Context. Within $Λ$ Cold Dark Matter ($Λ$CDM) simulations, Milky Way-like galaxies accrete some of their satellite galaxies in groups of 3-5 members rather than individually. It was also suggested that this might be the reason behind the origin of satellite planes. Objects accreted in groups are expected to share similar specific total energy and angular momentum, and also identical orbital planes and directions. Aims. Looking at observations of Milky Way satellites, the dwarf galaxies Leo II, IV, V, Crater II, and the star cluster Crater 1 were proposed to be a vestige of group infall. The suggested "Crater-Leo group" shows a coherent distance gradient and all these objects align along a great circle on the sky. We use proper motion data to investigate whether the phase-space distribution of the members of the proposed group are indeed consistent with group infall. Methods. To further investigate this possibility, we use Gaia Data Release 3 (DR3) and new Hubble Space Telescope (HST) proper motions, $(μ_{α*}, μ_δ) = (-0.1921 \pm 0.0514, -0.0686 \pm 0.0523)$ mas yr$^{-1}$ for Leo IV and $(μ_{α*}, μ_δ) = (0.1186 \pm 0.1943, -0.1183 \pm 0.1704)$ mas yr$^{-1}$ for Leo V, to derive accurate orbital properties for the proposed group objects. In addition, we explore other possible members of this putative association. Results. Leo II, Leo IV, and Crater 1 show orbital properties consistent with those we predict from assuming group infall. However, our results suggest that Crater II was not accreted with the rest of the objects. If confirmed with increasingly accurate proper motions in the future, the Crater-Leo objects appear to constitute the first identified case of a cosmologically expected, typical group infall event, as opposed to the highly hierarchical Magellanic Cloud system. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A on 23-04-2024. 20 pages, 18 figures, 5 tables

Deep Hubble Space Telescope Photometry of LMC and Milky Way Ultra-Faint Dwarfs: A careful look into the magnitude-size relation

Authors: Hannah Richstein, Nitya Kallivayalil, Joshua D. Simon, Christopher T. Garling, Andrew Wetzel, Jack T. Warfield, Roeland P. van der Marel, Myoungwon Jeon, Jonah C. Rose, Paul Torrey, Anna Claire Engelhardt, Gurtina Besla, Yumi Choi, Marla Geha, Puragra Guhathakurta, Evan N. Kirby, Ekta Patel, Elena Sacchi, Sangmo Tony Sohn

Abstract: We present deep Hubble Space Telescope (HST) photometry of ten targets from Treasury Program GO-14734, including six confirmed ultra-faint dwarf galaxies (UFDs), three UFD candidates, and one likely globular cluster. Six of these targets are satellites of, or have interacted with, the Large Magellanic Cloud (LMC). We determine their structural parameters using a maximum-likelihood technique. Using… ▽ More We present deep Hubble Space Telescope (HST) photometry of ten targets from Treasury Program GO-14734, including six confirmed ultra-faint dwarf galaxies (UFDs), three UFD candidates, and one likely globular cluster. Six of these targets are satellites of, or have interacted with, the Large Magellanic Cloud (LMC). We determine their structural parameters using a maximum-likelihood technique. Using our newly derived half-light radius ($r_h$) and $V$-band magnitude ($M_V$) values in addition to literature values for other UFDs, we find that UFDs associated with the LMC do not show any systematic differences from Milky Way UFDs in the magnitude-size plane. Additionally, we convert simulated UFD properties from the literature into the $M_V-r_h$ observational space to examine the abilities of current dark matter (DM) and baryonic simulations to reproduce observed UFDs. Some of these simulations adopt alternative DM models, thus allowing us to also explore whether the $M_V-r_h$ plane could be used to constrain the nature of DM. We find no differences in the magnitude-size plane between UFDs simulated with cold, warm, and self-interacting dark matter, but note that the sample of UFDs simulated with alternative DM models is quite limited at present. As more deep, wide-field survey data become available, we will have further opportunities to discover and characterize these ultra-faint stellar systems and the greater low surface-brightness universe. △ Less

Submitted 13 February, 2024; originally announced February 2024.

Comments: 42 pages, 13 figures, 13 tables, 2 appendices; Submitted to AAS Journals

Proper Motions and Orbits of Distant Local Group Dwarf Galaxies from a combination of Gaia and Hubble Data

Authors: Paul Bennet, Ekta Patel, Sangmo Tony Sohn, Andres del Pino, Roeland P. van der Marel, Mattia Libralato, Laura L. Watkins, Antonio Aparicio, Gurtina Besla, Carme Gallart, Mark A. Fardal, Matteo Monelli, Elena Sacchi, Erik Tollerud, Daniel R. Weisz

Abstract: We have determined the proper motions (PMs) of 12 dwarf galaxies in the Local Group (LG), ranging from the outer Milky Way (MW) halo to the edge of the LG. We used HST as the first and Gaia as the second epoch using the GaiaHub software. For Leo A and Sag DIG we also used multi-epoch HST measurements relative to background galaxies. Orbital histories derived using these PMs show that two-thirds of… ▽ More We have determined the proper motions (PMs) of 12 dwarf galaxies in the Local Group (LG), ranging from the outer Milky Way (MW) halo to the edge of the LG. We used HST as the first and Gaia as the second epoch using the GaiaHub software. For Leo A and Sag DIG we also used multi-epoch HST measurements relative to background galaxies. Orbital histories derived using these PMs show that two-thirds of the galaxies in our sample are on first infall with $>$90\% certainty. The observed star formation histories (SFHs) of these first-infall dwarfs are generally consistent with infalling dwarfs in simulations. The remaining four galaxies have crossed the virial radius of either the MW or M31. When we compare their star formation (SF) and orbital histories we find tentative agreement between the inferred pattern of SF with the timing of dynamical events in the orbital histories. For Leo~I, SF activity rises as the dwarf crosses the MW's virial radius, culminating in a burst of SF shortly before pericenter ($\approx1.7$~Gyr ago). The SF then declines after pericenter, but with some smaller bursts before its recent quenching ($\approx0.3$~Gyr ago). This shows that even small dwarfs like Leo~I can hold on to gas reservoirs and avoid quenching for several Gyrs after falling into their host, which is longer than generally found in simulations. Leo~II, NGC~6822, and IC~10 are also qualitatively consistent with this SF pattern in relation to their orbit, but more tentatively due to larger uncertainties. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: 27 pages, 9 Figures, 8 Tables, Submitted to ApJ

BP3M: Bayesian Positions, Parallaxes, and Proper Motions derived from the Hubble Space Telescope and Gaia data

Authors: Kevin A. McKinnon, Andrés del Pino, Constance M. Rockosi, Miranda Apfel, Puragra Guhathakurta, Roeland P. van der Marel, Paul Bennet, Mark A. Fardal, Mattia Libralato, Sangmo Tony Sohn, Eduardo Vitral, Laura L. Watkins

Abstract: We present a hierarchical Bayesian pipeline, BP3M, that measures positions, parallaxes, and proper motions (PMs) for cross-matched sources between Hubble~Space~Telescope (HST) images and Gaia -- even for sparse fields ($N_*<10$ per image) -- expanding from the recent GaiaHub tool. This technique uses Gaia-measured astrometry as priors to predict the locations of sources in HST images, and is there… ▽ More We present a hierarchical Bayesian pipeline, BP3M, that measures positions, parallaxes, and proper motions (PMs) for cross-matched sources between Hubble~Space~Telescope (HST) images and Gaia -- even for sparse fields ($N_*<10$ per image) -- expanding from the recent GaiaHub tool. This technique uses Gaia-measured astrometry as priors to predict the locations of sources in HST images, and is therefore able to put the HST images onto a global reference frame without the use of background galaxies/QSOs. Testing our publicly-available code in the Fornax and Draco dSphs, we measure accurate PMs that are a median of 8-13 times more precise than Gaia DR3 alone for $20.5<G<21~\mathrm{mag}$. We are able to explore the effect of observation strategies on BP3M astrometry using synthetic data, finding an optimal strategy to improve parallax and position precision at no cost to the PM uncertainty. Using 1619 HST images in the sparse COSMOS field (median 9 Gaia sources per HST image), we measure BP3M PMs for 2640 unique sources in the $16<G<21.5~\mathrm{mag}$ range, 25% of which have no Gaia PMs; the median BP3M PM uncertainty for $20.25<G<20.75~\mathrm{mag}$ sources is $0.44~$mas/yr compared to $1.03~$mas/yr from Gaia, while the median BP3M PM uncertainty for sources without Gaia-measured PMs ($20.75<G<21.5~\mathrm{mag}$) is $1.16~$mas/yr. The statistics that underpin the BP3M pipeline are a generalized way of combining position measurements from different images, epochs, and telescopes, which allows information to be shared between surveys and archives to achieve higher astrometric precision than that from each catalog alone. △ Less

Submitted 4 November, 2023; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 33 pages, 25 figures, 3 tables

JWST-TST High Contrast: Achieving direct spectroscopy of faint substellar companions next to bright stars with the NIRSpec IFU

Authors: Jean-Baptiste Ruffio, Marshall D. Perrin, Kielan K. W. Hoch, Jens Kammerer, Quinn M. Konopacky, Laurent Pueyo, Alex Madurowicz, Emily Rickman, Christopher A. Theissen, Shubh Agrawal, Alexandra Z. Greenbaum, Brittany E. Miles, Travis S. Barman, William O. Balmer, Jorge Llop-Sayson, Julien H. Girard, Isabel Rebollido, Rémi Soummer, Natalie H. Allen, Jay Anderson, Charles A. Beichman, Andrea Bellini, Geoffrey Bryden, Néstor Espinoza, Ana Glidden , et al. (11 additional authors not shown)

Abstract: The JWST NIRSpec integral field unit (IFU) presents a unique opportunity to observe directly imaged exoplanets from 3-5 um at moderate spectral resolution (R~2,700) and thereby better constrain the composition, disequilibrium chemistry, and cloud properties of their atmospheres. In this work, we present the first NIRSpec IFU high-contrast observations of a substellar companion that requires starli… ▽ More The JWST NIRSpec integral field unit (IFU) presents a unique opportunity to observe directly imaged exoplanets from 3-5 um at moderate spectral resolution (R~2,700) and thereby better constrain the composition, disequilibrium chemistry, and cloud properties of their atmospheres. In this work, we present the first NIRSpec IFU high-contrast observations of a substellar companion that requires starlight suppression techniques. We develop specific data reduction strategies to study faint companions around bright stars, and assess the performance of NIRSpec at high contrast. First, we demonstrate an approach to forward model the companion signal and the starlight directly in the detector images, which mitigates the effects of NIRSpec's spatial undersampling. We demonstrate a sensitivity to planets that are 3e-6 fainter than their stars at 1'', or 3e-5 at 0.3''. Then, we implement a reference star point spread function (PSF) subtraction and a spectral extraction that does not require spatially and spectrally regularly sampled spectral cubes. This allows us to extract a moderate resolution (R~2,700) spectrum of the faint T-dwarf companion HD 19467 B from 2.9-5.2 um with signal-to-noise ratio (S/N)~10 per resolution element. Across this wavelength range, HD~19467~B has a flux ratio varying between 1e-5-1e-4 and a separation relative to its star of 1.6''. A companion paper by Hoch et al. more deeply analyzes the atmospheric properties of this companion based on the extracted spectrum. Using the methods developed here, NIRSpec's sensitivity may enable direct detection and spectral characterization of relatively old (~1 Gyr), cool (~250 K), and closely separated (~3-5 au) exoplanets that are less massive than Jupiter. △ Less

Submitted 31 May, 2024; v1 submitted 15 October, 2023; originally announced October 2023.

Comments: Accepted to AJ. The data analysis scripts for this work are published https://github.com/jruffio/HD_19467_B (https://zenodo.org/doi/10.5281/zenodo.11391740). The main revisions of the manuscript are listed in the change history section of the readme

JWST-TST DREAMS: Quartz Clouds in the Atmosphere of WASP-17b

Authors: David Grant, Nikole K. Lewis, Hannah R. Wakeford, Natasha E. Batalha, Ana Glidden, Jayesh Goyal, Elijah Mullens, Ryan J. MacDonald, Erin M. May, Sara Seager, Kevin B. Stevenson, Jeff A. Valenti, Channon Visscher, Lili Alderson, Natalie H. Allen, Caleb I. Cañas, Knicole Colón, Mark Clampin, Néstor Espinoza, Amélie Gressier, Jingcheng Huang, Zifan Lin, Douglas Long, Dana R. Louie, Maria Peña-Guerrero , et al. (17 additional authors not shown)

Abstract: Clouds are prevalent in many of the exoplanet atmospheres that have been observed to date. For transiting exoplanets, we know if clouds are present because they mute spectral features and cause wavelength-dependent scattering. While the exact composition of these clouds is largely unknown, this information is vital to understanding the chemistry and energy budget of planetary atmospheres. In this… ▽ More Clouds are prevalent in many of the exoplanet atmospheres that have been observed to date. For transiting exoplanets, we know if clouds are present because they mute spectral features and cause wavelength-dependent scattering. While the exact composition of these clouds is largely unknown, this information is vital to understanding the chemistry and energy budget of planetary atmospheres. In this work, we observe one transit of the hot Jupiter WASP-17b with JWST's MIRI LRS and generate a transmission spectrum from 5-12 $\rmμ$m. These wavelengths allow us to probe absorption due to the vibrational modes of various predicted cloud species. Our transmission spectrum shows additional opacity centered at 8.6 $\rmμ$m, and detailed atmospheric modeling and retrievals identify this feature as SiO$_2$(s) (quartz) clouds. The SiO$_2$(s) clouds model is preferred at 3.5-4.2$σ$ versus a cloud-free model and at 2.6$σ$ versus a generic aerosol prescription. We find the SiO$_2$(s) clouds are comprised of small ${\sim}0.01$ $\rmμ$m particles, which extend to high altitudes in the atmosphere. The atmosphere also shows a depletion of H$_2$O, a finding consistent with the formation of high-temperature aerosols from oxygen-rich species. This work is part of a series of studies by our JWST Telescope Scientist Team (JWST-TST), in which we will use Guaranteed Time Observations to perform Deep Reconnaissance of Exoplanet Atmospheres through Multi-instrument Spectroscopy (DREAMS). △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 19 pages, 7 figures, accepted for publication in ApJL

RomAndromeda: The Roman Survey of the Andromeda Halo

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

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

Submitted 21 June, 2023; originally announced June 2023.

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

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

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

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

Submitted 20 June, 2023; originally announced June 2023.

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

The Near Infrared Imager and Slitless Spectrograph for the James Webb Space Telescope -- I. Instrument Overview and in-Flight Performance

Authors: Rene Doyon, C. J Willott, John B. Hutchings, Anand Sivaramakrishnan, Loic Albert, David Lafreniere, Neil Rowlands, M. Begona Vila, Andre R. Martel, Stephanie LaMassa, David Aldridge, Etienne Artigau, Peter Cameron, Pierre Chayer, Neil J. Cook, Rachel A. Cooper, Antoine Darveau-Bernier, Jean Dupuis, Colin Earnshaw, Nestor Espinoza, Joseph C. Filippazzo, Alexander W. Fullerton, Daniel Gaudreau, Roman Gawlik, Paul Goudfrooij , et al. (38 additional authors not shown)

Abstract: The Near-Infrared Imager and Slitless Spectrograph (NIRISS) is the science module of the Canadian-built Fine Guidance Sensor (FGS) onboard the James Webb Space Telescope (JWST). NIRISS has four observing modes: 1) broadband imaging featuring seven of the eight NIRCam broadband filters, 2) wide-field slitless spectroscopy (WFSS) at a resolving power of $\sim$150 between 0.8 and 2.2 $μ$m, 3) single-… ▽ More The Near-Infrared Imager and Slitless Spectrograph (NIRISS) is the science module of the Canadian-built Fine Guidance Sensor (FGS) onboard the James Webb Space Telescope (JWST). NIRISS has four observing modes: 1) broadband imaging featuring seven of the eight NIRCam broadband filters, 2) wide-field slitless spectroscopy (WFSS) at a resolving power of $\sim$150 between 0.8 and 2.2 $μ$m, 3) single-object cross-dispersed slitless spectroscopy (SOSS) enabling simultaneous wavelength coverage between 0.6 and 2.8 $μ$m at R$\sim$700, a mode optimized for exoplanet spectroscopy of relatively bright ($J<6.3$) stars and 4) aperture masking interferometry (AMI) between 2.8 and 4.8 $μ$m enabling high-contrast ($\sim10^{-3}-10^{-4}$) imaging at angular separations between 70 and 400 milliarcsec for relatively bright ($M<8$) sources. This paper presents an overview of the NIRISS instrument, its design, its scientific capabilities, and a summary of in-flight performance. NIRISS shows significantly better response shortward of $\sim2.5\,μ$m resulting in 10-40% sensitivity improvement for broadband and low-resolution spectroscopy compared to pre-flight predictions. Two time-series observations performed during instrument commissioning in the SOSS mode yield very stable spectro-photometry performance within $\sim$10% of the expected noise. The first space-based companion detection of the tight binary star AB Dor AC through AMI was demonstrated. △ Less

Submitted 5 June, 2023; originally announced June 2023.

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

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

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

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

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

The James Webb Space Telescope Mission

Authors: Jonathan P. Gardner, John C. Mather, Randy Abbott, James S. Abell, Mark Abernathy, Faith E. Abney, John G. Abraham, Roberto Abraham, Yasin M. Abul-Huda, Scott Acton, Cynthia K. Adams, Evan Adams, David S. Adler, Maarten Adriaensen, Jonathan Albert Aguilar, Mansoor Ahmed, Nasif S. Ahmed, Tanjira Ahmed, Rüdeger Albat, Loïc Albert, Stacey Alberts, David Aldridge, Mary Marsha Allen, Shaune S. Allen, Martin Altenburg , et al. (983 additional authors not shown)

Abstract: Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astrono… ▽ More Twenty-six years ago a small committee report, building on earlier studies, expounded a compelling and poetic vision for the future of astronomy, calling for an infrared-optimized space telescope with an aperture of at least $4m$. With the support of their governments in the US, Europe, and Canada, 20,000 people realized that vision as the $6.5m$ James Webb Space Telescope. A generation of astronomers will celebrate their accomplishments for the life of the mission, potentially as long as 20 years, and beyond. This report and the scientific discoveries that follow are extended thank-you notes to the 20,000 team members. The telescope is working perfectly, with much better image quality than expected. In this and accompanying papers, we give a brief history, describe the observatory, outline its objectives and current observing program, and discuss the inventions and people who made it possible. We cite detailed reports on the design and the measured performance on orbit. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted by PASP for the special issue on The James Webb Space Telescope Overview, 29 pages, 4 figures

JWST-TST Proper Motions: I. High-Precision NIRISS Calibration and Large Magellanic Cloud Kinematics

Authors: M. Libralato, A. Bellini, R. P. van der Marel, J. Anderson, S. T. Sohn, L. L. Watkins, L. Alderson, N. Allen, M. Clampin, A. Glidden, J. Goyal, K. Hoch, J. Huang, J. Kammerer, N. K. Lewis, Z. Lin, D. Long, D. Louie, R. J. MacDonald, M. Mountain, M. Peña-Guerrero, M. D. Perrin, L. Pueyo, I. Rebollido, E. Rickman , et al. (5 additional authors not shown)

Abstract: We develop and disseminate effective point-spread functions and geometric-distortion solutions for high-precision astrometry and photometry with the JWST NIRISS instrument. We correct field dependencies and detector effects, and assess the quality and the temporal stability of the calibrations. As a scientific application and validation, we study the proper motion (PM) kinematics of stars in the J… ▽ More We develop and disseminate effective point-spread functions and geometric-distortion solutions for high-precision astrometry and photometry with the JWST NIRISS instrument. We correct field dependencies and detector effects, and assess the quality and the temporal stability of the calibrations. As a scientific application and validation, we study the proper motion (PM) kinematics of stars in the JWST calibration field near the Large Magellanic Cloud (LMC) center, comparing to a first-epoch Hubble Space Telescope (HST) archival catalog with a 16-yr baseline. For stars with G~20, the median PM uncertainty is ~13 $μ$as yr$^{-1}$ (3.1 km s$^{-1}$), better than Gaia DR3 typically achieves for its very best-measured stars. We kinematically detect the known star cluster OGLE-CL LMC 407, measure its absolute PM for the first time, and show how this differs from other LMC populations. The inferred cluster dispersion sets an upper limit of 24 $μ$as yr$^{-1}$ (5.6 km s$^{-1}$) on systematic uncertainties. Red-giant-branch stars have a velocity dispersion of 33.8 $\pm$ 0.6 km s$^{-1}$, while younger blue populations have a narrower velocity distribution, but with a significant kinematical substructure. We discuss how this relates to the larger velocity dispersions inferred from Gaia DR3. These results establish JWST as capable of state-of-the-art astrometry, building on the extensive legacy of HST. This is the first paper in a series by our JWST Telescope Scientist Team (TST), in which we will use Guaranteed Time Observations to study the PM kinematics of various stellar systems in the Local Group. △ Less

Submitted 25 April, 2023; v1 submitted 28 February, 2023; originally announced March 2023.

Comments: 28 pages, 25 figures, 2 tables. Accepted for publication in ApJ. The effective point-spread-function models, the geometric-distortion solutions and a preliminary version of the code are available at the links provided in the manuscript

Classifying the satellite plane membership of Centaurus A's dwarf galaxies using orbital alignment constraints

Authors: Kosuke Jamie Kanehisa, Marcel S. Pawlowski, Oliver Müller, Sangmo Tony Sohn

Abstract: The flattened, possibly co-rotating plane of satellite galaxies around Centaurus A, if more than a fortuitous alignment, adds to the pre-existing tension between the well-studied Milky Way and M31 planes and the $Λ$CDM model of structure formation. It was recently reported that the Centaurus A satellite plane (CASP) may be rotationally supported, but a further understanding of the system's kinemat… ▽ More The flattened, possibly co-rotating plane of satellite galaxies around Centaurus A, if more than a fortuitous alignment, adds to the pre-existing tension between the well-studied Milky Way and M31 planes and the $Λ$CDM model of structure formation. It was recently reported that the Centaurus A satellite plane (CASP) may be rotationally supported, but a further understanding of the system's kinematics is elusive in the absence of full three-dimensional velocities. We constrain the transverse velocities of 27 satellites that would rotationally stabilise the Centaurus A plane, and classify the satellites by whether their possible orbits are consistent with the CASP. Five satellites are identified to be unlikely to participate in the plane, two of which are clearly non-members. Despite their previously reported line-of-sight velocity trend suggestive of a common co-rotating motion, 17 out of 22 potential CASP members are consistent with either orbital direction within both the full range of possible kinematics as well as when limiting orbits to those within the plane. On the other hand, disregarding the 5 off-plane satellites found to be inconsistent with CASP membership enhances the significance of the CASP's line-of-sight velocity trend fivefold. Our results are robust with different mass estimates of the Centaurus A halo, and the adoption of either spherical or triaxial NFW potentials. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 15 pages, 8 figures, 3 tables, accepted for publication in MNRAS

The James Webb Space Telescope Mission: Optical Telescope Element Design, Development, and Performance

Authors: Michael W. McElwain, Lee D. Feinberg, Marshall D. Perrin, Mark Clampin, C. Matt Mountain, Matthew D. Lallo, Charles-Philippe Lajoie, Randy A. Kimble, Charles W. Bowers, Christopher C. Stark, D. Scott Acton, Ken Aiello, Charles Atkinson, Beth Barinek, Allison Barto, Scott Basinger, Tracy Beck, Matthew D. Bergkoetter, Marcel Bluth, Rene A. Boucarut, Gregory R. Brady, Keira J. Brooks, Bob Brown, John Byard, Larkin Carey , et al. (104 additional authors not shown)

Abstract: The James Webb Space Telescope (JWST) is a large, infrared space telescope that has recently started its science program which will enable breakthroughs in astrophysics and planetary science. Notably, JWST will provide the very first observations of the earliest luminous objects in the Universe and start a new era of exoplanet atmospheric characterization. This transformative science is enabled by… ▽ More The James Webb Space Telescope (JWST) is a large, infrared space telescope that has recently started its science program which will enable breakthroughs in astrophysics and planetary science. Notably, JWST will provide the very first observations of the earliest luminous objects in the Universe and start a new era of exoplanet atmospheric characterization. This transformative science is enabled by a 6.6 m telescope that is passively cooled with a 5-layer sunshield. The primary mirror is comprised of 18 controllable, low areal density hexagonal segments, that were aligned and phased relative to each other in orbit using innovative image-based wavefront sensing and control algorithms. This revolutionary telescope took more than two decades to develop with a widely distributed team across engineering disciplines. We present an overview of the telescope requirements, architecture, development, superb on-orbit performance, and lessons learned. JWST successfully demonstrates a segmented aperture space telescope and establishes a path to building even larger space telescopes. △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: accepted by PASP for JWST Overview Special Issue; 34 pages, 25 figures

The Near Infrared Imager and Slitless Spectrograph for the James Webb Space Telescope -- IV. Aperture Masking Interferometry

Authors: Anand Sivaramakrishnan, Peter Tuthill, James P. Lloyd, Alexandra Z. Greenbaum, Deepashri Thatte, Rachel A. Cooper, Thomas Vandal, Jens Kammerer, Joel Sanchez-Bermudez, Benjamin J. S. Pope, Dori Blakely, Loïc Albert, Neil J. Cook, Doug Johnstone, André R. Martel, Kevin Volk, Anthony Soulain, Étienne Artigau, David Lafrenière, Chris J. Willott, Sébastien Parmentier, K. E. Saavik Ford, Barry McKernan, M. Begoña Vila, Neil Rowlands , et al. (14 additional authors not shown)

Abstract: The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, pres… ▽ More The James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) flies a 7-hole non-redundant mask (NRM), the first such interferometer in space, operating at 3-5 \micron~wavelengths, and a bright limit of $\simeq 4$ magnitudes in W2. We describe the NIRISS Aperture Masking Interferometry (AMI) mode to help potential observers understand its underlying principles, present some sample science cases, explain its operational observing strategies, indicate how AMI proposals can be developed with data simulations, and how AMI data can be analyzed. We also present key results from commissioning AMI. Since the allied Kernel Phase Imaging (KPI) technique benefits from AMI operational strategies, we also cover NIRISS KPI methods and analysis techniques, including a new user-friendly KPI pipeline. The NIRISS KPI bright limit is $\simeq 8$ W2 magnitudes. AMI (and KPI) achieve an inner working angle of $\sim 70$ mas that is well inside the $\sim 400$ mas NIRCam inner working angle for its circular occulter coronagraphs at comparable wavelengths. △ Less

Submitted 7 November, 2022; v1 submitted 31 October, 2022; originally announced October 2022.

Comments: 30 pages, 10 figures

HubPUG: Proper Motions for Local Group Dwarfs observed with HST utilizing Gaia as a Reference Frame

Authors: Jack T. Warfield, Nitya Kallivayalil, Paul Zivick, Tobias Fritz, Hannah Richstein, Sangmo Tony Sohn, Andrés del Pino, Alessandro Savino, Daniel R. Weisz

Abstract: We present the method behind HubPUG, a software tool built for recovering systemic proper motions (PMs) of Hubble Space Telescope (HST) fields with two epochs of observations by utilizing stars observed by Gaia as a foreground frame of reference. HST PM experiments have typically relied on the use of distant background galaxies or quasi-stellar objects (QSOs) as stationary sources against which to… ▽ More We present the method behind HubPUG, a software tool built for recovering systemic proper motions (PMs) of Hubble Space Telescope (HST) fields with two epochs of observations by utilizing stars observed by Gaia as a foreground frame of reference. HST PM experiments have typically relied on the use of distant background galaxies or quasi-stellar objects (QSOs) as stationary sources against which to measure PMs. Without consistent profiles, background galaxies are more difficult to centroid, but benefit on-aggregate from their large numbers. QSOs, though they can be fit with stellar point-spread functions, are sparse, with most fields containing none. Historically, the use of stars as references against which to measure PMs would have been difficult because they have individual PMs of their own. However, Gaia has now provided positions and PMs for over 1.4 billion stars, which are much more likely to be well-imaged in the fields around targets versus background sources, have predictable stellar profiles, and require less observing time per-image for good signal-to-noise. This technique allows us to utilize the power of Gaia to measure the PM of targets too faint for Gaia to observe itself. We have recovered PMs for the Milky Way satellites Sculptor and Draco with comparable uncertainties over HST-only and Gaia-only measurements, limited primarily by the current capabilities of the Gaia data. We also show the promise of this method for satellites of M31 with a new PM measurement for Andromeda VII. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Astrometric and Wavelength Calibration of the NIRSpec Instrument during Commissioning using a model-based approach

Authors: Nora Lützgendorf, Giovanna Giardino, Catarina Alves de Oliveira, Peter Zeidler, Pierre Ferruit, Peter Jakobsen, Nimisha Kumari, Timothy Rawle, Stephan M. Birkmann, Torsten Böker, Charles Proffitt, Marco Sirianni, Maurice Te Plate, Sangmo Tony Sohn

Abstract: The NIRSpec instrument for the James Webb Space Telescope (JWST) is a highly versatile near-infrared spectrograph that can be operated in various observing modes, slit apertures, and spectral resolutions. Obtaining dedicated calibration data for all possible combinations of aperture and disperser is an intractable task. We have therefore developed a procedure to derive a highly realistic model of… ▽ More The NIRSpec instrument for the James Webb Space Telescope (JWST) is a highly versatile near-infrared spectrograph that can be operated in various observing modes, slit apertures, and spectral resolutions. Obtaining dedicated calibration data for all possible combinations of aperture and disperser is an intractable task. We have therefore developed a procedure to derive a highly realistic model of the instrument's optical geometry across the entire field of view, using calibration data acquired through only a subset of NIRSpec apertures, which nevertheless allows the light paths within the spectrograph to be accurately computed for all apertures and all observing modes. This parametric instrument model thus provides the basis for the extraction of wavelength-calibrated spectra from any NIRSpec exposure, regardless of observing mode or aperture used. Optimizing the NIRSpec instrument model and deriving its final wavelength and astrometric calibration was one of the most crucial elements of the NIRSpec commissioning phase. Here, we describe the process of re-fitting the NIRSpec instrument model with in-orbit commissioning data, and present its final performance in terms of wavelength accuracy and astrometric calibration. △ Less

Submitted 10 August, 2022; originally announced August 2022.

Comments: 15 pages, to appear in Proceedings of SPIE conference "Space Telescopes and Instrumentation 2022: Optical, Infrared, and Millimeter Wave"

JWST/NIRCam Coronagraphy: Commissioning and First On-Sky Results

Authors: Julien H. Girard, Jarron Leisenring, Jens Kammerer, Mario Gennaro, Marcia Rieke, John Stansberry, Armin Rest, Eiichi Egami, Ben Sunnquist, Martha Boyer, Alicia Canipe, Matteo Correnti, Bryan Hilbert, Marshall D. Perrin, Laurent Pueyo, Remi Soummer, Marsha Allen, Howard Bushouse, Jonathan Aguilar, Brian Brooks, Dan Coe, Audrey DiFelice, David Golimowski, George Hartig, Dean C. Hines , et al. (31 additional authors not shown)

Abstract: In a cold and stable space environment, the James Webb Space Telescope (JWST or "Webb") reaches unprecedented sensitivities at wavelengths beyond 2 microns, serving most fields of astrophysics. It also extends the parameter space of high-contrast imaging in the near and mid-infrared. Launched in late 2021, JWST underwent a six month commissioning period. In this contribution we focus on the NIRCam… ▽ More In a cold and stable space environment, the James Webb Space Telescope (JWST or "Webb") reaches unprecedented sensitivities at wavelengths beyond 2 microns, serving most fields of astrophysics. It also extends the parameter space of high-contrast imaging in the near and mid-infrared. Launched in late 2021, JWST underwent a six month commissioning period. In this contribution we focus on the NIRCam Coronagraphy mode which was declared "science ready" on July 10 2022, the last of the 17 JWST observing modes. Essentially, this mode will allow to detect fainter/redder/colder (less massive for a given age) self-luminous exoplanets as well as other faint astrophysical signal in the vicinity of any bright object (stars or galaxies). Here we describe some of the steps and hurdles the commissioning team went through to achieve excellent performances. Specifically, we focus on the Coronagraphic Suppression Verification activity. We were able to produce firm detections at 3.35$μ$m of the white dwarf companion HD 114174 B which is at a separation of $\simeq$ 0.5" and a contrast of $\simeq$ 10 magnitudes ($10^{4}$ fainter than the K$\sim$5.3 mag host star). We compare these first on-sky images with our latest, most informed and realistic end-to-end simulations through the same pipeline. Additionally we provide information on how we succeeded with the target acquisition with all five NIRCam focal plane masks and their four corresponding wedged Lyot stops. △ Less

Submitted 31 August, 2022; v1 submitted 1 August, 2022; originally announced August 2022.

Comments: 23 pages, 18 figures, to appear in the Proceedings of the SPIE Astronomical Telescopes + Instrumentation (2022)

Kinematic Structure of the Large Magellanic Cloud Globular Cluster System from Gaia eDR3 and Hubble Space Telescope Proper Motions

Authors: Paul Bennet, Mayte Alfaro-Cuello, Andrés del Pino, Laura L. Watkins, Roeland P. van der Marel, Sangmo Tony Sohn

Abstract: We have determined bulk proper motions (PMs) for 31 LMC GCs from Gaia eDR3 and Hubble Space Telescope data using multiple independent analysis techniques. Combined with literature values for distances, line-of-sight velocities and existing bulk PMs, we extract full 6D phase-space information for 32 clusters, allowing us to examine the kinematics of the LMC GC system in detail. Except for two GCs (… ▽ More We have determined bulk proper motions (PMs) for 31 LMC GCs from Gaia eDR3 and Hubble Space Telescope data using multiple independent analysis techniques. Combined with literature values for distances, line-of-sight velocities and existing bulk PMs, we extract full 6D phase-space information for 32 clusters, allowing us to examine the kinematics of the LMC GC system in detail. Except for two GCs (NGC 2159 and NGC 2210) for which high velocities suggest they are not long-term members of the LMC system, the data are consistent with a flattened configuration that rotates like the stellar disk. The one-dimensional velocity dispersions are of order 30 km/s, similar to that of old stellar populations in the LMC disk. Similar to the case for Milky Way disk clusters, the velocity anisotropy is such that the dispersion is smallest in the azimuthal direction; however, alternative anisotropies cannot be ruled out due to distance uncertainties. The data are consistent with a single multi-dimensional Gaussian velocity distribution. Given the non-collisional nature of the LMC disk, this suggests that most, if not all, of the LMC GCs are formed by a single formation mechanism in the stellar disk, despite a significant spread in age and metallicity. Any accreted halo GC population is absent or far smaller in the LMC compared to the Milky Way. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted by ApJ, 19 pages, 5 figures, 4 tables

The Science Performance of JWST as Characterized in Commissioning

Authors: Jane Rigby, Marshall Perrin, Michael McElwain, Randy Kimble, Scott Friedman, Matt Lallo, René Doyon, Lee Feinberg, Pierre Ferruit, Alistair Glasse, Marcia Rieke, George Rieke, Gillian Wright, Chris Willott, Knicole Colon, Stefanie Milam, Susan Neff, Christopher Stark, Jeff Valenti, Jim Abell, Faith Abney, Yasin Abul-Huda, D. Scott Acton, Evan Adams, David Adler , et al. (601 additional authors not shown)

Abstract: This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries f… ▽ More This paper characterizes the actual science performance of the James Webb Space Telescope (JWST), as determined from the six month commissioning period. We summarize the performance of the spacecraft, telescope, science instruments, and ground system, with an emphasis on differences from pre-launch expectations. Commissioning has made clear that JWST is fully capable of achieving the discoveries for which it was built. Moreover, almost across the board, the science performance of JWST is better than expected; in most cases, JWST will go deeper faster than expected. The telescope and instrument suite have demonstrated the sensitivity, stability, image quality, and spectral range that are necessary to transform our understanding of the cosmos through observations spanning from near-earth asteroids to the most distant galaxies. △ Less

Submitted 10 April, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 5th version as accepted to PASP; 31 pages, 18 figures; https://iopscience.iop.org/article/10.1088/1538-3873/acb293

Journal ref: PASP 135 048001 (2023)

The Hubble Space Telescope UV Legacy Survey of Galactic Globular Clusters. XXIII. Proper-motion catalogs and internal kinematics

Authors: M. Libralato, A. Bellini, E. Vesperini, G. Piotto, A. P. Milone, R. P. van der Marel, J. Anderson, A. Aparicio, B. Barbuy, L. R. Bedin, L. Borsato, S. Cassisi, E. Dalessandro, F. R. Ferraro, I. R. King, B. Lanzoni, D. Nardiello, S. Ortolani, A. Sarajedini, S. T. Sohn

Abstract: A number of studies based on data collected by the $\textit{Hubble Space Telescope}$ ($\textit{HST}$) GO-13297 program "HST Legacy Survey of Galactic Globular Clusters: Shedding UV Light on Their Populations and Formation" have investigated the photometric properties of a large sample of Galactic globular clusters and revolutionized our understanding of their stellar populations. In this paper, we… ▽ More A number of studies based on data collected by the $\textit{Hubble Space Telescope}$ ($\textit{HST}$) GO-13297 program "HST Legacy Survey of Galactic Globular Clusters: Shedding UV Light on Their Populations and Formation" have investigated the photometric properties of a large sample of Galactic globular clusters and revolutionized our understanding of their stellar populations. In this paper, we expand previous studies by focusing our attention on the stellar clusters' internal kinematics. We computed proper motions for stars in 56 globular and one open clusters by combining the GO-13297 images with archival $\textit{HST}$ data. The astro-photometric catalogs released with this paper represent the most complete and homogeneous collection of proper motions of stars in the cores of stellar clusters to date, and expand the information provided by the current (and future) $\textit{Gaia}$ data releases to much fainter stars and into the crowded central regions. We also census the general kinematic properties of stellar clusters by computing the velocity-dispersion and anisotropy radial profiles of their bright members. We study the dependence on concentration and relaxation time, and derive dynamical distances. Finally, we present an in-depth kinematic analysis of the globular cluster NGC 5904. △ Less

Submitted 5 July, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: 33 pages, 19 figures, 6 tables. Accepted for publication on ApJ. Astro-photometric catalogs, velocity-dispersion values and profiles are available at https://archive.stsci.edu/hlsp/hacks

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

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

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

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

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

GaiaHub: A method for combining data from the Gaia and Hubble space telescopes to derive improved proper motions for faint stars

Authors: Andrés del Pino, Mattia Libralato, Roeland P. van der Marel, Paul Bennet, Mark A. Fardal, Jay Anderson, Andrea Bellini, Sangmo Tony Sohn, Laura L. Watkins

Abstract: We present GaiaHub, a publicly available tool that combines $Gaia$ measurements with $Hubble$ $Space$ $Telescope$ ($HST$) archival images to derive proper motions (PMs). It increases the scientific impact of both observatories beyond their individual capabilities. $Gaia$ provides PMs across the whole sky, but the limited mirror size and time baseline restrict the best PM performance to relatively… ▽ More We present GaiaHub, a publicly available tool that combines $Gaia$ measurements with $Hubble$ $Space$ $Telescope$ ($HST$) archival images to derive proper motions (PMs). It increases the scientific impact of both observatories beyond their individual capabilities. $Gaia$ provides PMs across the whole sky, but the limited mirror size and time baseline restrict the best PM performance to relatively bright stars. $HST$ can measure accurate PMs for much fainter stars over a small field, but this requires two epochs of observation which are not always available. GaiaHub yields considerably improved PM accuracy compared to $Gaia$-only measurements, especially for faint sources $(G \gtrsim 18)$, requiring only a single epoch of $HST$ data observed more than $\sim 7$ years ago (before 2012). This provides considerable scientific value especially for dynamical studies of stellar systems or structures in and beyond the Milky Way (MW) halo, for which the member stars are generally faint. To illustrate the capabilities and demonstrate the accuracy of GaiaHub, we apply it to samples of MW globular clusters (GCs) and classical dwarf spheroidal (dSph) satellite galaxies. This allows us, e.g., to measure the velocity dispersions in the plane of the sky for objects out to and beyond $\sim 100$ kpc. We find, on average, mild radial velocity anisotropy in GCs, consistent with existing results for more nearby samples. We observe a correlation between the internal kinematics of the clusters and their ellipticity, with more isotropic clusters being, on average, more round. Our results also support previous findings that Draco and Sculptor dSph galaxies appear to be radially anisotropic systems. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 24 pages, 18 figures. Accepted for publication in ApJ

Structural parameters and possible association of the Ultra-Faint Dwarfs Pegasus III and Pisces II from deep Hubble Space Telescope photometry

Authors: Hannah Richstein, Ekta Patel, Nitya Kallivayalil, Joshua D. Simon, Paul Zivick, Erik Tollerud, Tobias Fritz, Jack T. Warfield, Gurtina Besla, Roeland P. van der Marel, Andrew Wetzel, Yumi Choi, Alis Deason, Marla Geha, Puragra Guhathakurta, Myoungwon Jeon, Evan N. Kirby, Mattia Libralato, Elena Sacchi, Sangmo Tony Sohn

Abstract: We present deep Hubble Space Telescope (HST) photometry of the ultra-faint dwarf (UFD) galaxies Pegasus III (Peg III) and Pisces II (Psc II), two of the most distant satellites in the halo of the Milky Way (MW). We measure the structure of both galaxies, derive mass-to-light ratios with newly determined absolute magnitudes, and compare our findings to expectations from UFD-mass simulations. For Pe… ▽ More We present deep Hubble Space Telescope (HST) photometry of the ultra-faint dwarf (UFD) galaxies Pegasus III (Peg III) and Pisces II (Psc II), two of the most distant satellites in the halo of the Milky Way (MW). We measure the structure of both galaxies, derive mass-to-light ratios with newly determined absolute magnitudes, and compare our findings to expectations from UFD-mass simulations. For Peg III, we find an elliptical half-light radius of $a_h=1.88^{+0.42}_{-0.33}$ arcminutes ($118^{+31}_{-30}$ pc) and $M_V{=}{-4.17}^{+0.19}_{-0.22}$; for Psc II, we measure $a_h{=}1.31^{+0.10}_{-0.09}$ arcminutes ($69\pm8$ pc) and $M_V{=}{-4.28}^{+0.19}_{-0.16}$. We do not find any morphological features that indicate a significant interaction between the two has occurred, despite their close separation of only $\sim$40 kpc. Using proper motions (PMs) from Gaia early Data Release 3, we investigate the possibility of any past association by integrating orbits for the two UFDs in a MW-only and a combined MW and Large Magellanic Cloud (LMC) potential. We find that including the gravitational influence of the LMC is crucial, even for these outer-halo satellites, and that a possible orbital history exists where Peg III and Psc II experienced a close ($\sim$10-20 kpc) passage about each other just over $\sim$1 Gyr ago, followed by a collective passage around the LMC ($\sim$30-60 kpc) just under $\sim$1 Gyr ago. Considering the large uncertainties on the PMs and the restrictive priors imposed to derive them, improved PM measurements for Peg III and Psc II will be necessary to clarify their relationship. This would add to the rare findings of confirmed pairs of satellites within the Local Group. △ Less

Submitted 13 February, 2024; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 24 pages, 13 figures; 1 appendix; published in ApJ

Journal ref: The Astrophysical Journal, Volume 933, Issue 2, id.217, July 2022

Detection of Chemo-Kinematical Structures in Leo I

Authors: A. G. Alarcon Jara, M. Fellhauer, J. Simon, A. del Pino, S. W. Fu, S. T. Sohn

Abstract: Context: A variety of formation models for dwarf spheroidal (dSph) galaxies have been proposed in the literature, but generally they have not been quantitatively compared with observations. Aims: We search for chemodynamical patterns in our observational data set and compare the results with mock galaxies consisting of pure random motions, and simulated dwarfs formed via the dissolving star cluste… ▽ More Context: A variety of formation models for dwarf spheroidal (dSph) galaxies have been proposed in the literature, but generally they have not been quantitatively compared with observations. Aims: We search for chemodynamical patterns in our observational data set and compare the results with mock galaxies consisting of pure random motions, and simulated dwarfs formed via the dissolving star cluster and tidal stirring models. Methods: We made use of a new spectroscopic data set for the Milky Way dSph Leo I, combining 288 stars observed with Magellan/IMACS and existing Keck/DEIMOS data, to provide velocity and metallicity measurements for 953 Leo I member stars. We used a specially developed algorithm called Beacon to detect chemo-kinematical patterns in the observed and simulated data. Results: After analysing the Leo I data, we report the detection of 14 candidate streams of stars that may have originated in disrupted star clusters. The angular momentum vectors of these streams are randomly oriented, consistent with the lack of rotation in Leo I. These results are consistent with the predictions of the dissolving cluster model. In contrast, we find fewer candidate stream signals in mock data sets that lack coherent motions ~99% of the time. The chemodynamical analysis of the tidal stirring simulation produces streams that share a common orientation of their angular momenta, which is inconsistent with the Leo I data. Conclusions: Even though it is very difficult to distinguish which of the detected streams are real and which are only noise, we can be certain that there are more streams detected in the observational data of Leo I than expected in pure random data. △ Less

Submitted 2 March, 2023; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: 12 pages, accepted by A&A

Journal ref: A&A 672, A131 (2023)

On the Co-Orbitation of Satellite Galaxies Along the Great Plane of Andromeda: NGC 147, NGC 185, and Expectations from Cosmological Simulations

Authors: Marcel S. Pawlowski, Sangmo Tony Sohn

Abstract: Half of the satellite galaxies of Andromeda form a narrow plane termed the Great Plane of Andromeda (GPoA), and their line-of-sight velocities display correlation reminiscent of a rotating structure. Recently reported first proper motion measurements for the on-plane satellites NGC 147 and NGC 185 indicate that they indeed co-orbit along the GPoA. This provides a novel opportunity to compare the M… ▽ More Half of the satellite galaxies of Andromeda form a narrow plane termed the Great Plane of Andromeda (GPoA), and their line-of-sight velocities display correlation reminiscent of a rotating structure. Recently reported first proper motion measurements for the on-plane satellites NGC 147 and NGC 185 indicate that they indeed co-orbit along the GPoA. This provides a novel opportunity to compare the M31 satellite system to $Λ$CDM expectations. We perform the first detailed comparison of the orbital alignment of two satellite galaxies beyond the Milky Way with several hydrodynamical and dark-matter-only cosmological simulations (Illustris TNG-50, TNG-100, ELVIS, PhatELVIS), in the context of the Planes of Satellite Galaxies Problem. In line with previous works, we find that the spatial flattening and line-of-sight velocity correlation alone is already in substantial tension with $Λ$CDM, with none of the simulated analogs simultaneously reproducing both parameters. Almost none (3 to 4\%) of the simulated systems contain two satellites with orbital poles as well aligned with their satellite plane as indicated by the most-likely proper motions of NGC 147 and NGC 185. However, within current measurement uncertainties it is common (~70%) that the two best-aligned satellites of simulated systems are consistent with the orbital alignment. Yet, the chance that any two simulated on-plane satellites have as well aligned orbital poles as observed is low (~4%). We conclude that confirmation of the tight orbital alignment for these two objects via improved measurements, or the discovery of similar alignments for additional GPoA members, hold the potential to further raise the tension with $Λ$CDM expectations. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: 27 pages, 13 figures, 5 tables, accepted for publication in ApJ

Nonlinear Color-Metallicity Relations of Globular Clusters. X. Subaru/FOCAS Multi-object Spectroscopy of M87 Globular Clusters

Authors: Sooyoung Kim, Suk-Jin Yoon, Sang-Yoon Lee, Chul Chung, Sangmo Tony Sohn

Abstract: We obtained spectra of some 140 globular clusters (GCs) associated with the Virgo central cD galaxy M87 with the Subaru/FOCAS MOS mode. The fundamental properties of GCs such as age, metallicity and $α$-element abundance are investigated by using simple stellar population models. It is confirmed that the majority of M87 GCs are as old as, more metal-rich than, and more enhanced in $α$-elements tha… ▽ More We obtained spectra of some 140 globular clusters (GCs) associated with the Virgo central cD galaxy M87 with the Subaru/FOCAS MOS mode. The fundamental properties of GCs such as age, metallicity and $α$-element abundance are investigated by using simple stellar population models. It is confirmed that the majority of M87 GCs are as old as, more metal-rich than, and more enhanced in $α$-elements than the Milky Way GCs. Our high-quality, homogeneous dataset enables us to test the theoretical prediction of inflected color$-$metallicity relations (CMRs). The nonlinear-CMR hypothesis entails an alternative explanation for the widely observed GC color bimodality, in which even a unimodal metallicity spread yields a bimodal color distribution by virtue of nonlinear metallicity-to-color conversion. The newly derived CMRs of old, high-signal-to-noise-ratio GCs in M87 (the $V-I$ CMR of 83 GCs and the $M-T2$ CMR of 78 GCs) corroborate the presence of the significant inflection. Furthermore, from a combined catalog with the previous study on M87 GC spectroscopy, we find that a total of 185 old GCs exhibit a broad, unimodal metallicity distribution. The results corroborate the nonlinear-CMR interpretation of the GC color bimodality, shedding further light on theories of galaxy formation. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: 28 pages, 15 figures, 8 tables

Star Formation Histories of Ultra-Faint Dwarf Galaxies: environmental differences between Magellanic and non-Magellanic satellites?

Authors: Elena Sacchi, Hannah Richstein, Nitya Kallivayalil, Roeland van der Marel, Mattia Libralato, Paul Zivick, Gurtina Besla, Thomas M. Brown, Yumi Choi, Alis Deason, Tobias Fritz, Marla Geha, Puragra Guhathakurta, Myoungwon Jeon, Evan Kirby, Steven R. Majewski, Ekta Patel, Joshua D. Simon, Sangmo Tony Sohn, Erik Tollerud, Andrew Wetzel

Abstract: We present the color-magnitude diagrams and star formation histories (SFHs) of seven ultra-faint dwarf galaxies: Horologium 1, Hydra 2, Phoenix 2, Reticulum 2, Sagittarius 2, Triangulum 2, and Tucana 2, derived from high-precision Hubble Space Telescope photometry. We find that the SFH of each galaxy is consistent with them having created at least 80% of the stellar mass by $z\sim6$. For all galax… ▽ More We present the color-magnitude diagrams and star formation histories (SFHs) of seven ultra-faint dwarf galaxies: Horologium 1, Hydra 2, Phoenix 2, Reticulum 2, Sagittarius 2, Triangulum 2, and Tucana 2, derived from high-precision Hubble Space Telescope photometry. We find that the SFH of each galaxy is consistent with them having created at least 80% of the stellar mass by $z\sim6$. For all galaxies, we find quenching times older than 11.5 Gyr ago, compatible with the scenario in which reionization suppresses the star formation of small dark matter halos. However, our analysis also reveals some differences in the SFHs of candidate Magellanic Cloud satellites, i.e., galaxies that are likely satellites of the Large Magellanic Cloud and that entered the Milky Way potential only recently. Indeed, Magellanic satellites show quenching times about 600 Myr more recent with respect to those of other Milky Way satellites, on average, even though the respective timings are still compatible within the errors. This finding is consistent with theoretical models that suggest that satellites' SFHs may depend on their host environment at early times, although we caution that within the error bars all galaxies in our sample are consistent with being quenched at a single epoch. △ Less

Submitted 28 September, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 7 pages, 3 figures, 2 tables. Accepted for publication in ApJL

Mapping Gaia parallax systematic errors over the sky with faint Milky Way stars

Authors: Mark A. Fardal, Roeland van der Marel, Andres del Pino, Sangmo Tony Sohn

Abstract: Parallaxes measured by the Gaia mission have huge significance for astronomy, but parallaxes in Gaia DR2 are known to have systematic errors that depend on the source position and other quantities. We use the abundant information in faint Milky Way stars, along with the GOG simulation of the Gaia catalog, to probe the spatial dependence of Gaia DR2 parallax systematic errors in an empirical way. T… ▽ More Parallaxes measured by the Gaia mission have huge significance for astronomy, but parallaxes in Gaia DR2 are known to have systematic errors that depend on the source position and other quantities. We use the abundant information in faint Milky Way stars, along with the GOG simulation of the Gaia catalog, to probe the spatial dependence of Gaia DR2 parallax systematic errors in an empirical way. The parallax signal, concentrated in thick disk turnoff stars with magnitude G ~ 17, is sufficient to construct maps of the parallax systematic error over the majority of the sky. These maps show a locally regular "waffle pattern" on ~1 degree scales following Gaia scan directions, stronger linear "scar" features, and coherent variations on larger scales. The parallax bias maps also retain traces of astrophysical effects such as dust clouds. The waffle pattern, known from earlier maps of the Magellanic Clouds, extends over the entire sky; its local rms amplitude averages 15 microarcsec and varies by about a factor of two. The strength of this pattern increases by a factor ~6 from magnitude G = 13 to G = 20. Correlations with parallaxes of quasars and of stars with independent distance estimates support our bias estimates. Using similar methods, we map systematic errors in the proper motion and examine the relationship with the parallax systematics. We provide a code package to access and query our bias maps. Similar tests on the general stellar population should be useful in quantifying systematic errors in future Gaia releases. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: Accepted in AJ. 17 pages, 11 figures. Higher resolution figures available in journal. Eagerly anticipating Gaia EDR3

Revealing the Structure and Internal Rotation of the Sagittarius Dwarf Spheroidal Galaxy with Gaia and Machine Learning

Authors: Andrés del Pino, Mark A. Fardal, Roeland P. van der Marel, Ewa L. Łokas, Cecilia Mateu, Sangmo Tony Sohn

Abstract: We present a detailed study of the internal structure and kinematics of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). Using machine-learning techniques, we have combined the information provided by 3300 RR Lyrae stars, more than 2000 spectroscopically observed stars, and the Gaia second data release to derive the full phase space, i.e. 3D positions and kinematics, of more than… ▽ More We present a detailed study of the internal structure and kinematics of the core of the Sagittarius dwarf spheroidal galaxy (Sgr). Using machine-learning techniques, we have combined the information provided by 3300 RR Lyrae stars, more than 2000 spectroscopically observed stars, and the Gaia second data release to derive the full phase space, i.e. 3D positions and kinematics, of more than $1.2\times10^5$ member stars in the core of the galaxy. Our results show that Sgr has a bar structure $\sim 2.5$ kpc long, and that tidal tails emerge from its tips to form what it is known as the Sgr stream. The main body of the galaxy, strongly sheared by tidal forces, is a triaxial (almost prolate) ellipsoid with its longest principal axis of inertia inclined $43^\circ\pm6^\circ$ with respect to the plane of the sky and axis ratios of 1:0.67:0.60. Its external regions are expanding mainly along its longest principal axis, yet the galaxy conserves an inner core of about $500\times330\times300$ pc that shows no net expansion and is rotating at $v_{\rm rot} = 4.13 \pm 0.16$ ${\rm{ km \ s^{-1}}}$. The internal angular momentum of Sgr forms an angle $θ= 18^\circ\pm6^\circ$ with respect to its orbital angular momentum, meaning that the galaxy is in an inclined prograde orbit around the Milky Way. We compared our results with predictions from $N$-body models with spherical, pressure-supported progenitors and a model whose progenitor is a flattened rotating disk. Only the rotating model, based on preexisting simulations aimed at reproducing the line-of-sight velocity gradients observed in Sgr, was able to reproduce the observed properties in the core of the galaxy. △ Less

Submitted 15 February, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 37 pages, 26 figures, 2 tables. Accepted for publication in the Astrophysical Journal

2D kinematics of massive stars near the Galactic Center

Authors: Mattia Libralato, Daniel J. Lennon, Andrea Bellini, Roeland van der Marel, Simon J. Clark, Francisco Najarro, Lee R. Patrick, Jay Anderson, Luigi R. Bedin, Paul A. Crowther, Selma E. de Mink, Christopher J. Evans, Imants Platais, Elena Sabbi, Sangmo Tony Sohn

Abstract: The presence of massive stars (MSs) in the region close to the Galactic Center (GC) poses several questions about their origin. The harsh environment of the GC favors specific formation scenarios, each of which should imprint characteristic kinematic features on the MSs. We present a 2D kinematic analysis of MSs in a GC region surrounding Sgr A* based on high-precision proper motions obtained with… ▽ More The presence of massive stars (MSs) in the region close to the Galactic Center (GC) poses several questions about their origin. The harsh environment of the GC favors specific formation scenarios, each of which should imprint characteristic kinematic features on the MSs. We present a 2D kinematic analysis of MSs in a GC region surrounding Sgr A* based on high-precision proper motions obtained with the Hubble Space Telescope. Thanks to a careful data reduction, well-measured bright stars in our proper-motion catalogs have errors better than 0.5 mas yr$^{-1}$. We discuss the absolute motion of the MSs in the field and their motion relative to Sgr A*, the Arches and the Quintuplet. For the majority of the MSs, we rule out any distance further than 3-4 kpc from Sgr A* using only kinematic arguments. If their membership to the GC is confirmed, most of the isolated MSs are likely not associated with either the Arches or Quintuplet clusters or Sgr A*. Only a few MSs have proper motions suggesting they are likely members of the Arches cluster, in agreement with previous spectroscopic results. Line-of-sight radial velocities and distances are required to shed further light on the origin of most of these massive objects. We also present an analysis of other fast-moving objects in the GC region, finding no clear excess of high-velocity escaping stars. We make our astro-photometric catalogs publicly available. △ Less

Submitted 13 November, 2020; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 28 pages, 33 figures, 9 tables. Accepted for publication in MNRAS. [v3: Fixed bibliography]

HST Proper Motions of NGC 147 and NGC 185: Orbital Histories and Test of Dynamically Coherent Andromeda Satellite Plane

Authors: Sangmo Tony Sohn, Ekta Patel, Mark A. Fardal, Gurtina Besla, Roeland P. van der Marel, Marla Geha, Puragra Guhathakurta

Abstract: We present the first proper motion (PM) measurements for the dwarf elliptical galaxies NGC 147 and NGC 185, two satellite galaxies of M31, using multi-epoch HST imaging data with time baselines of $\sim 8$ years. For each galaxy, we take an error-weighted average of measurements from HST ACS/WFC and WFC3/UVIS to determine the PMs. Our final results for the PMs are… ▽ More We present the first proper motion (PM) measurements for the dwarf elliptical galaxies NGC 147 and NGC 185, two satellite galaxies of M31, using multi-epoch HST imaging data with time baselines of $\sim 8$ years. For each galaxy, we take an error-weighted average of measurements from HST ACS/WFC and WFC3/UVIS to determine the PMs. Our final results for the PMs are $(μ_\mathrm{W}, μ_\mathrm{N})_\mathrm{N147} = (-0.0232, 0.0378) \pm (0.0143, 0.0146)\,\mathrm{mas}\,\mathrm{yr}^{-1}$ for NGC 147, and $(μ_\mathrm{W}, μ_\mathrm{N})_\mathrm{N185} = (-0.0242, 0.0058) \pm (0.0141, 0.0147)\,\mathrm{mas}\,\mathrm{yr}^{-1}$ for NGC 185. The 2-dimensional direction of motion for NGC 147 about M31 is found to be aligned with its tidal tails. The 3-d positions and velocities of both galaxies are transformed into a common M31-centric coordinate system to study the detailed orbital histories of the combined M31+NGC 147+NGC 185 system via numerical orbit integration. We find that NGC 147 (NGC 185) had its closest passage to M31 0.3-0.5~Gyr ($\gtrsim 1.6$~Gyr) within the past 6 Gyr at distances of $\sim 70$ kpc (70-260 kpc). The pericentric times of NGC 147/NGC 185 correlate qualitatively well with the presence/absence of tidal tails seen around the galaxies. Our PMs show that the orbital poles of NGC 147, and also NGC 185 albeit to a lesser degree, agree within the uncertainties with the normal of the Great Plane of Andromeda (GPoA). These are the first measurements of the 3-d angular momentum vector of any satellite identified as original GPoA members. Our results strengthen the hypothesis that the GPoA may be a dynamically coherent entity. We revisit previous claims that NGC 147 and NGC 185 are binary galaxies and conclude that it is very unlikely the two galaxies were ever gravitationally bound to each other. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: 20 pages, 9 figures, 4 tables, accepted for publication in ApJ

A kinematic view of NGC 1261: structural parameters, internal dispersion, absolute proper motion and Blue Straggler Stars

Authors: S. Raso, M. Libralato, A. Bellini, F. R. Ferraro, B. Lanzoni, M. Cadelano, C. Pallanca, E. Dalessandro, G. Piotto, J. Anderson, S. T. Sohn

Abstract: We constructed a Hubble Space Telescope (HST) astro-photometric catalog of the central region of the Galactic globular cluster NGC 1261. This catalog, complemented with Gaia DR2 data sampling the external regions, has been used to estimate the structural parameters of the system (i.e., core, half-mass, tidal radii and concentration) from its resolved star density profile. We computed high-precisio… ▽ More We constructed a Hubble Space Telescope (HST) astro-photometric catalog of the central region of the Galactic globular cluster NGC 1261. This catalog, complemented with Gaia DR2 data sampling the external regions, has been used to estimate the structural parameters of the system (i.e., core, half-mass, tidal radii and concentration) from its resolved star density profile. We computed high-precision proper motions thanks to multi-epoch HST data and derived the cluster velocity dispersion profile in the plane of the sky for the innermost region, finding that the system is isotropic. The combination with line-of-sight information collected from spectroscopy in the external regions provided us with the cluster velocity dispersion profile along the entire radial extension. We also measured the absolute proper motion of NGC 1261 using a few background galaxies as a reference. The radial distribution of the Blue Straggler Star population shows that the cluster is in a low/intermediate phase of dynamical evolution. △ Less

Submitted 20 April, 2020; originally announced April 2020.

Comments: 19 pages, 12 figures, 2 tables; accepted for publication in ApJ

Astro2020 Science White Paper: Construction of an L* Galaxy: the Transformative Power of Wide Fields for Revealing the Past, Present and Future of the Great Andromeda System

Authors: Karoline M. Gilbert, Erik J. Tollerud, Jay Anderson, Rachael L. Beaton, Eric F. Bell, Alyson Brooks, Thomas M. Brown, James Bullock, Jeffrey L. Carlin, Michelle Collins, Andrew Cooper, Denija Crnojevic, Julianne Dalcanton, Andres del Pino, Richard D'Souza, Ivanna Escala, Mark Fardal, Andreea Font, Marla Geha, Puragra Guhathakurta, Evan Kirby, Geraint F. Lewis, Jennifer L. Marshall, Nicolas F. Martin, Kristen McQuinn , et al. (12 additional authors not shown)

Abstract: The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover,… ▽ More The Great Andromeda Galaxy (M31) is the nexus of the near-far galaxy evolution connection and a principal data point for near-field cosmology. Due to its proximity (780 kpc), M31 can be resolved into individual stars like the Milky Way (MW). Unlike the MW, we have the advantage of a global view of M31, enabling M31 to be observed with techniques that also apply to more distant galaxies. Moreover, recent evidence suggests that M31 may have survived a major merger within the last several Gyr, shaping the morphology of its stellar halo and triggering a starburst, while leaving the stellar disk largely intact. The MW and M31 thus provide complementary opportunities for in-depth studies of the disks, halos, and satellites of L* galaxies. Our understanding of the M31 system will be transformed in the 2020s if they include wide field facilities for both photometry (HST-like sensitivity and resolution) and spectroscopy (10-m class telescope, >1 sq. deg. field, highly multiplexed, R~ 3000 to 6000). We focus here on the power of these facilities to constrain the past, present, and future merger history of M31, via chemo-dynamical analyses and star formation histories of phase-mixed stars accreted at early times, as well as stars in surviving tidal debris features, M31's extended disk, and intact satellite galaxies that will eventually be tidally incorporated into the halo. This will yield an unprecedented view of the hierarchical formation of the M31 system and the subhalos that built it into the L* galaxy we observe today. △ Less

Submitted 1 April, 2019; originally announced April 2019.

Comments: Submitted as a science white paper to the Astro2020 Decadal Survey

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

Astro2020 Science White Paper: Science at the edges: internal kinematics of globular clusters' external fields

Authors: A. Bellini, M. Libralato, J. Anderson, D. Bennett, A. Calamida, S. Casertano, S. M. Fall, B. S. Gaudi, P. Guhathakurta, S. Ho, J. Lu, S. Malhotra, P. Melchior, E. Nelan, J. Rhodes, R. E. Sanderson, M. Shao, S. T. Sohn, E. Vesperini, R. P. van der Marel

Abstract: The outer regions of globular clusters can enable us to answer many fundamental questions concerning issues ranging from the formation and evolution of clusters and their multiple stellar populations to the study of stars near and beyond the hydrogen-burning limit and to the dynamics of the Milky Way. The outskirts of globular clusters are still uncharted territories observationally. A very effici… ▽ More The outer regions of globular clusters can enable us to answer many fundamental questions concerning issues ranging from the formation and evolution of clusters and their multiple stellar populations to the study of stars near and beyond the hydrogen-burning limit and to the dynamics of the Milky Way. The outskirts of globular clusters are still uncharted territories observationally. A very efficient way to explore them is through high-precision proper motions of low-mass stars over a large field of view. The Wide Field InfraRed Survey Telescope (WFIRST) combines all these characteristics in a single telescope, making it the best observational tool to uncover the wealth of information contained in the clusters' outermost regions. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: 8 pages, 3 figures

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

The course of the Orphan Stream in the Northern Galactic Hemisphere traced with Gaia DR2

Authors: Mark A. Fardal, Roeland P. van der Marel, Sangmo Tony Sohn, Andres del Pino Molina

Abstract: The Orphan Stream is one of the most prominent tidal streams in the Galactic halo. Using data on red giants, RR Lyrae, and horizontal branch stars from Gaia and other surveys, we determine the proper motion of the Orphan Stream over a path of more than 90 degrees on the sky. We also provide updated tracks for the sky position, distance, and radial velocity of the stream. Our tracks in these latter… ▽ More The Orphan Stream is one of the most prominent tidal streams in the Galactic halo. Using data on red giants, RR Lyrae, and horizontal branch stars from Gaia and other surveys, we determine the proper motion of the Orphan Stream over a path of more than 90 degrees on the sky. We also provide updated tracks for the sky position, distance, and radial velocity of the stream. Our tracks in these latter dimensions mostly agree with previous results. However, there are significant corrections to the earlier distance and latitude tracks as the stream approaches the Galactic disk. Stream stars selected with three-dimensional kinematics display a very tight red giant sequence. Concordantly, we find that applying a proper motion cut removes the most metal-rich stars from earlier spectroscopic samples of stream stars, though a significant dispersion remains indicating a dwarf galaxy origin. The deceleration of the stream towards its leading end suggests a circular velocity of ~200 km/s at a galactocentric radius ~30 kpc, consistent with other independent evidence. However, the track of the stream departs significantly from an orbit; the spatial track does not point along the same direction as the velocity vector, and it exhibits a lateral wiggle that is unlikely to match any reasonable orbit. The low metallicity and small dispersion of the stream in the various coordinates point to a progenitor with a relatively low dynamical mass ~$10^8$ Msun. △ Less

Submitted 14 December, 2018; originally announced December 2018.

Comments: 14 pages, 10 figures, 1 table; submitted to MNRAS

MSC Class: The course of the Orphan Stream in the Northern Galactic Hemisphere traced with Gaia DR2

The Proper Motion Field Along the Magellanic Bridge: a New Probe of the LMC-SMC Interaction

Authors: Paul Zivick, Nitya Kallivayalil, Gurtina Besla, Sangmo Tony Sohn, Roeland P. van der Marel, Andrés del Pino, Sean T. Linden, Tobias K. Fritz, J. Anderson

Abstract: We present the first detailed kinematic analysis of the proper motions (PMs) of stars in the Magellanic Bridge, from both the \textit{Gaia} Data Release 2 catalog and from \textit{Hubble Space Telescope} Advanced Camera for Surveys data. For the \textit{Gaia} data, we identify and select two populations of stars in the Bridge region, young main sequence (MS) and red giant stars. The spatial locati… ▽ More We present the first detailed kinematic analysis of the proper motions (PMs) of stars in the Magellanic Bridge, from both the \textit{Gaia} Data Release 2 catalog and from \textit{Hubble Space Telescope} Advanced Camera for Surveys data. For the \textit{Gaia} data, we identify and select two populations of stars in the Bridge region, young main sequence (MS) and red giant stars. The spatial locations of the stars are compared against the known H {\small I} gas structure, finding a correlation between the MS stars and the H {\small I} gas. In the \textit{Hubble Space Telescope} fields our signal comes mainly from an older MS and turn-off population, and the proper motion baselines range between $\sim 4$ and 13 years. The PMs of these different populations are found to be consistent with each other, as well as across the two telescopes. When the absolute motion of the Small Magellanic Cloud is subtracted out, the residual Bridge motions display a general pattern of pointing away from the Small Magellanic Cloud towards the Large Magellanic Cloud. We compare in detail the kinematics of the stellar samples against numerical simulations of the interactions between the Small and Large Magellanic Clouds, and find general agreement between the kinematics of the observed populations and a simulation in which the Clouds have undergone a recent direct collision. △ Less

Submitted 18 February, 2019; v1 submitted 22 November, 2018; originally announced November 2018.

Comments: 13 pages, 10 figures, 2 tables, submitted to ApJ, accepted February 8th, 2019

HALO7D II: The Halo Velocity Ellipsoid and Velocity Anisotropy with Distant Main Sequence Stars

Authors: Emily C. Cunningham, Alis J. Deason, Robyn E. Sanderson, Sangmo Tony Sohn, Jay Anderson, Puragra Guhathakurta, Constance M. Rockosi, Roeland P. van der Marel, Sarah R. Loebman, Andrew Wetzel

Abstract: The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use… ▽ More The Halo Assembly in Lambda-CDM: Observations in 7 Dimensions (HALO7D) dataset consists of Keck II/DEIMOS spectroscopy and Hubble Space Telescope-measured proper motions of Milky Way (MW) halo main sequence turnoff stars in the CANDELS fields. In this paper, the second in the HALO7D series, we present the proper motions for the HALO7D sample. We discuss our measurement methodology, which makes use of a Bayesian mixture modeling approach for creating the stationary reference frame of distant galaxies. Using the 3D kinematic HALO7D sample, we estimate the parameters of the halo velocity ellipsoid, $\langle v_φ \rangle, σ_r, σ_φ, σ_θ$, and the velocity anisotropy $β$. Using the full HALO7D sample, we find $β=0.63 \pm 0.05$ at $\langle r \rangle =24$ kpc. We also estimate the ellipsoid parameters for our sample split into three apparent magnitude bins; the posterior medians for these estimates of $β$, while consistent with one another, increase as a function of mean sample distance. Finally, we estimate $β$ in each of the individual HALO7D fields. We find that the velocity anisotropy $β$ can vary from field to field, which suggests that the halo is not phase mixed at $\langle r \rangle =24$ kpc. We explore the $β$ variation across the skies of two stellar halos from the \textit{Latte} suite of FIRE-2 simulations, finding that both simulated galaxies show $β$ variation over a similar range to the variation observed across the four HALO7D fields. The accretion histories of the two simulated galaxies result in different $β$ variation patterns; spatially mapping $β$ is thus a way forward in characterizing the accretion history of the Galaxy. △ Less

Submitted 29 October, 2018; originally announced October 2018.

Comments: 23 pages, 13 figures. Submitted to ApJ; comments welcome

First Gaia Dynamics of the Andromeda System: DR2 Proper Motions, Orbits, and Rotation of M31 and M33

Authors: Roeland P. van der Marel, Mark A. Fardal, Sangmo Tony Sohn, Ekta Patel, Gurtina Besla, Andrés del Pino-Molina, Johannes Sahlmann, Laura L. Watkins

Abstract: The 3D velocities of M31 and M33 are important for understanding the evolution and cosmological context of the Local Group. Their most massive stars are detected by Gaia, and we use Data Release 2 (DR2) to determine the galaxy proper motions (PMs). We select galaxy members based on, e.g., parallax, PM, color-magnitude-diagram location, and local stellar density. The PM rotation of both galaxies is… ▽ More The 3D velocities of M31 and M33 are important for understanding the evolution and cosmological context of the Local Group. Their most massive stars are detected by Gaia, and we use Data Release 2 (DR2) to determine the galaxy proper motions (PMs). We select galaxy members based on, e.g., parallax, PM, color-magnitude-diagram location, and local stellar density. The PM rotation of both galaxies is confidently detected, consistent with the known line-of-sight rotation curves: $V_{\rm rot} = -206\pm86$ km s$^{-1}$ (counter-clockwise) for M31, and $V_{\rm rot} = 80\pm52$ km s$^{-1}$ (clockwise) for M33. We measure the center-of-mass PM of each galaxy relative to surrounding background quasars in DR2. This yields that $(μ_{α*},μ_δ)$ equals $(65 \pm 18 , -57 \pm 15)$ $μ$as yr$^{-1}$ for M31, and $(31 \pm 19 , -29 \pm 16)$ $μ$as yr$^{-1}$ for M33. In addition to the listed random errors, each component has an additional residual systematic error of 16 $μ$as yr$^{-1}$. These results are consistent at 0.8$σ$ and 1.0$σ$ with the (2 and 3 times higher-accuracy) measurements already available from Hubble Space Telescope (HST) optical imaging and VLBA water maser observations, respectively. This lends confidence that all these measurements are robust. The new results imply that the M31 orbit towards the Milky Way is somewhat less radial than previously inferred, $V_{\rm tan, DR2+HST} = 57^{+35}_{-31}$ km s$^{-1}$, and strengthen arguments that M33 may be on its first infall into M31. The results highlight the future potential of Gaia for PM studies beyond the Milky Way satellite system. △ Less

Submitted 28 January, 2019; v1 submitted 10 May, 2018; originally announced May 2018.

Comments: 15 pages, 7 figures, ApJ, in press

The Missing Satellites of the Magellanic Clouds? Gaia Proper Motions of the Recently Discovered Ultra-Faint Galaxies

Authors: Nitya Kallivayalil, Laura Sales, Paul Zivick, Tobias K. Fritz, Andrés Del Pino, Sangmo Tony Sohn, Gurtina Besla, Roeland P. van der Marel, Julio F. Navarro, Elena Sacchi

Abstract: According to LCDM theory, hierarchical evolution occurs on all mass scales, implying that satellites of the Milky Way should also have companions. The recent discovery of ultra-faint dwarf galaxy candidates in close proximity to the Magellanic Clouds provides an opportunity to test this theory. We present proper motion (PM) measurements for 13 of the 32 new dwarf galaxy candidates using Gaia data… ▽ More According to LCDM theory, hierarchical evolution occurs on all mass scales, implying that satellites of the Milky Way should also have companions. The recent discovery of ultra-faint dwarf galaxy candidates in close proximity to the Magellanic Clouds provides an opportunity to test this theory. We present proper motion (PM) measurements for 13 of the 32 new dwarf galaxy candidates using Gaia data release 2. All 13 also have radial velocity measurements. We compare the measured 3D velocities of these dwarfs to those expected at the corresponding distance and location for the debris of an LMC analog in a cosmological numerical simulation. We conclude that 4 of these galaxies (Hor1, Car2, Car3 and Hyi1) have come in with the Magellanic Clouds, constituting the first confirmation of the type of satellite infall predicted by LCDM. Ret2, Tuc2 and Gru1 have velocity components that are not consistent within 3 sigma of our predictions and are therefore less favorable. Hya2 and Dra2 could be associated with the LMC and merit further attention. We rule out Tuc3, Cra2, Tri2 and Aqu2 as potential members. Of the dwarfs without measured PMs, 5 of them are deemed unlikely on the basis of their positions and distances alone as being too far from the orbital plane expected for LMC debris (Eri2, Ind2, Cet2, Cet3 and Vir1). For the remaining sample, we use the simulation to predict PMs and radial velocities, finding that Phx2 has an overdensity of stars in DR2 consistent with this PM prediction. △ Less

Submitted 29 October, 2018; v1 submitted 3 May, 2018; originally announced May 2018.

Comments: 14 pages, 7 figures, 3 tables. Accepted to ApJ

Journal ref: The Astrophysical Journal, Volume 867, Issue 1, article id. 19, 10 pp. (2018)

Gaia DR2 Proper Motions of Dwarf Galaxies within 420 kpc: Orbits, Milky Way Mass, Tidal Influences, Planar Alignments, and Group Infall

Authors: T. K. Fritz, G. Battaglia, M. S. Pawlowski, N. Kallivayalil, R. van der Marel, T. S. Sohn, C. Brook, G. Besla

Abstract: A proper understanding of the Milky Way (MW) dwarf galaxies in a cosmological context requires knowledge of their 3D velocities and orbits. However, proper motion (PM) measurements have generally been of limited accuracy and available only for more massive dwarfs. We therefore present a new study of the kinematics of the MW dwarf galaxies. We use the Gaia DR2 for those dwarfs that have been spectr… ▽ More A proper understanding of the Milky Way (MW) dwarf galaxies in a cosmological context requires knowledge of their 3D velocities and orbits. However, proper motion (PM) measurements have generally been of limited accuracy and available only for more massive dwarfs. We therefore present a new study of the kinematics of the MW dwarf galaxies. We use the Gaia DR2 for those dwarfs that have been spectroscopically observed in the literature. We derive systemic PMs for 39 galaxies and galaxy candidates out to 420 kpc, and generally find good consistency for the subset with measurements available from other studies. We derive the implied Galactocentric velocities, and calculate orbits in canonical MW halo potentials of "low" ($0.8 \times 10^{12} M_\odot$) and "high" mass ($1.6 \times 10^{12} M_\odot$). Comparison of the distributions of orbital apocenters and 3D velocities to the halo virial radius and escape velocity, respectively, suggests that the satellite kinematics are best explained in the high-mass halo. Tuc III, Crater II, and additional candidates have orbital pericenters small enough to imply significant tidal influences. Relevant to the missing satellite problem, the fact that fewer galaxies are observed to be near apocenter than near pericenter implies that there must be a population of distant dwarf galaxies yet to be discovered. Of the 39 dwarfs: 12 have orbital poles that do not align with the MW plane of satellites (given reasonable assumptions about its intrinsic thickness); 10 have insufficient PM accuracy to establish whether they align; and 17 satellites align, of which 11 are co-orbiting and (somewhat surprisingly, in view of prior knowledge) 6 are counter-orbiting. Group infall might have contributed to this, but no definitive association is found for the members of the Crater-Leo group. △ Less

Submitted 29 August, 2018; v1 submitted 2 May, 2018; originally announced May 2018.

Comments: 18 pages, 7 figures, 4 tables, accepted for publication in A&A

Journal ref: A&A 619, A103 (2018)

Evidence for an Intermediate-Mass Milky Way from Gaia DR2 Halo Globular Cluster Motions

Authors: Laura L. Watkins, Roeland P. van der Marel, Sangmo Tony Sohn, N. Wyn Evans

Abstract: We estimate the mass of the Milky Way (MW) within 21.1 kpc using the kinematics of halo globular clusters (GCs) determined by Gaia. The second Gaia data release (DR2) contained a catalogue of absolute proper motions (PMs) for a set of Galactic GCs and satellite galaxies measured using Gaia DR2 data. We select from the catalogue only halo GCs, identifying a total of 34 GCs spanning… ▽ More We estimate the mass of the Milky Way (MW) within 21.1 kpc using the kinematics of halo globular clusters (GCs) determined by Gaia. The second Gaia data release (DR2) contained a catalogue of absolute proper motions (PMs) for a set of Galactic GCs and satellite galaxies measured using Gaia DR2 data. We select from the catalogue only halo GCs, identifying a total of 34 GCs spanning $2.0 < r < 21.1$ kpc, and use their 3D kinematics to estimate the anisotropy over this range to be $β= 0.46^{+0.15}_{-0.19}$, in good agreement, though slightly lower than, a recent estimate for a sample of halo GCs using HST PM measurements further out in the halo. We then use the Gaia kinematics to estimate the mass of the MW inside the outermost GC to be $M(< 21.1 \mathrm{kpc}) = 0.21^{+0.04}_{-0.03} 10^{12} \mathrm{M_\odot}$, which corresponds to a circular velocity of $v_\mathrm{circ}(21.1 \mathrm{kpc}) = 206^{+19}_{-16}$ km/s. The implied virial mass is $M_\mathrm{virial} = 1.28^{+0.97}_{-0.48} 10^{12} \mathrm{M_\odot}$. The error bars encompass the uncertainties on the anisotropy and on the density profile of the MW dark halo, and the scatter inherent in the mass estimator we use. We get improved estimates when we combine the Gaia and HST samples to provide kinematics for 46 GCs out to 39.5 kpc: $β= 0.52^{+0.11}_{-0.14}$, $M(< 39.5 \mathrm{kpc}) = 0.42^{+0.07}_{-0.06} 10^{12} \mathrm{M_\odot}$, and $M_\mathrm{virial} = 1.54^{+0.75}_{-0.44} 10^{12} \mathrm{M_\odot}$. We show that these results are robust to potential substructure in the halo GC distribution. While a wide range of MW virial masses have been advocated in the literature, from below $10^{12} \mathrm{M_\odot}$ to above $2 \times 10^{12}\mathrm{M_\odot}$, these new data imply that an intermediate mass is most likely. △ Less

Submitted 8 February, 2019; v1 submitted 30 April, 2018; originally announced April 2018.

Comments: 16 pages, 9 figures, 3 tables, accepted for publication in ApJ (updated to match accepted version)

HST astrometry in the 30 Doradus region: II. Runaway stars from new proper motions in the Large Magellanic Cloud

Authors: I. Platais, D. J. Lennon, R. P. van der Marel, A. Bellini, E. Sabbi, L. L. Watkins, S. T. Sohn, N. R. Walborn, L. R. Bedin, C. J. Evans, S. E. de Mink, H. Sana, A. Herrero, N. Langer, P. Crowther

Abstract: We present a catalog of relative proper motions for 368,787 stars in the 30 Doradus region of the Large Magellanic Cloud (LMC), based on a dedicated two-epoch survey with the Hubble Space Telescope (HST) and supplemented with proper motions from our pilot archival study. We demonstrate that a relatively short epoch difference of 3 years is sufficient to reach a $\sim$0.1 mas yr$^{-1}$ level of pre… ▽ More We present a catalog of relative proper motions for 368,787 stars in the 30 Doradus region of the Large Magellanic Cloud (LMC), based on a dedicated two-epoch survey with the Hubble Space Telescope (HST) and supplemented with proper motions from our pilot archival study. We demonstrate that a relatively short epoch difference of 3 years is sufficient to reach a $\sim$0.1 mas yr$^{-1}$ level of precision or better. A number of stars have relative proper motions exceeding a 3-sigma error threshold, representing a mixture of Milky Way denizens and 17 potential LMC runaway stars. Based upon 183 VFTS OB-stars with the best proper motions, we conclude that none of them move faster than $\sim$0.3 mas yr$^{-1}$ in each coordinate -- equivalent to $\sim$70 km s$^{-1}$. Among the remaining 351 VFTS stars with less accurate proper motions, only one candidate OB runaway can be identified. We rule out any OB star in our sample moving at a tangential velocity exceeding $\sim$120 km s$^{-1}$. The most significant result of this study is finding 10 stars over wide range of masses, which appear to be ejected from the massive star cluster R136 in the tangential plane to angular distances from $35^{\prime\prime}$ out to $407^{\prime\prime}$, equivalent to 8-98 pc. The tangential velocities of these runaways appear to be correlated with apparent magnitude, indicating a possible dependence on the stellar mass. △ Less

Submitted 23 April, 2018; originally announced April 2018.

Comments: 45 pages (in referee format), 12 figures, 3 tables. Submitted to AJ. Comments are welcome

Connecting the Milky Way potential profile to the orbital timescales and spatial structure of the Sagittarius Stream

Authors: Mark A. Fardal, Roeland P. van der Marel, David R. Law, Sangmo Tony Sohn, Branimir Sesar, Nina Hernitschek, Hans-Walter Rix

Abstract: Recent maps of the halo using RR Lyrae from Pan-STARRS1 have clearly depicted the spatial structure of the Sagittarius stream. These maps show the leading and trailing stream apocenters differ in galactocentric radius by a factor of two, and also resolve substructure in the stream at these apocenters. Here we present dynamical models that reproduce these features of the stream in simple Galactic p… ▽ More Recent maps of the halo using RR Lyrae from Pan-STARRS1 have clearly depicted the spatial structure of the Sagittarius stream. These maps show the leading and trailing stream apocenters differ in galactocentric radius by a factor of two, and also resolve substructure in the stream at these apocenters. Here we present dynamical models that reproduce these features of the stream in simple Galactic potentials. We find that debris at the apocenters must be dynamically young, in the sense of being stripped off in the last two pericentric passages, while the Sagittarius dwarf is currently experiencing a third passage. The ratio of apocenters is sensitive to both dynamical friction and the outer slope of the Galactic rotation curve. These dependences can be understood with simple regularities connecting the apocentric radii, circular velocities, and orbital period of the progenitor. The effect of dynamical friction on the stream can be constrained using substructure within the leading apocenter. Our ensembles of models are not intended as statistically proper fits to the stream. Nevertheless, out of the range of models we consider, we consistently find the mass within 100 kpc to be $\sim 7 \times 10^{11} \, M_{\odot}$, with a nearly flat rotation curve between 50 and 100 kpc. This points to a more extended Galactic halo than assumed in some current models. As in previous work, we find prolate or triaxial halos ease agreement with the track of the leading stream. We display the behavior of our models in various observational spaces and characterize the substructure expected within the stream. In particular, the young trailing stream visible near trailing apocenter should exhibit a tight trend of velocity with distance separate from the older debris, and we suggest that this will serve as an especially useful probe of the outer Galactic potential. △ Less

Submitted 13 April, 2018; originally announced April 2018.

Comments: Submitted to MNRAS

The Proper Motion Field of the Small Magellanic Cloud: Kinematic Evidence for its Tidal Disruption

Authors: Paul Zivick, Nitya Kallivayalil, Roeland P. van der Marel, Gurtina Besla, Sean T. Linden, Szymon Kozłowski, Tobias K. Fritz, C. S. Kochanek, J. Anderson, Sangmo Tony Sohn, Marla C. Geha, Charles R. Alcock

Abstract: We present a new measurement of the systemic proper motion of the Small Magellanic Cloud (SMC), based on an expanded set of 30 fields containing background quasars and spanning a $\sim$3 year baseline, using the \textit{Hubble Space Telescope} (\textit{HST}) Wide Field Camera 3. Combining this data with our previous 5 \textit{HST} fields, and an additional 8 measurements from the \textit{Gaia}-Tyc… ▽ More We present a new measurement of the systemic proper motion of the Small Magellanic Cloud (SMC), based on an expanded set of 30 fields containing background quasars and spanning a $\sim$3 year baseline, using the \textit{Hubble Space Telescope} (\textit{HST}) Wide Field Camera 3. Combining this data with our previous 5 \textit{HST} fields, and an additional 8 measurements from the \textit{Gaia}-Tycho Astrometric Solution Catalog, brings us to a total of 43 SMC fields. We measure a systemic motion of $μ_{W}$ = $-0.82$ $\pm$ 0.02 (random) $\pm$ 0.10 (systematic) mas yr$^{-1}$ and $μ_{N}$ = $-1.21$ $\pm$ 0.01 (random) $\pm$ 0.03 (systematic) mas yr$^{-1}$. After subtraction of the systemic motion, we find little evidence for rotation, but find an ordered mean motion radially away from the SMC in the outer regions of the galaxy, indicating that the SMC is in the process of tidal disruption. We model the past interactions of the Clouds with each other based on the measured present-day relative velocity between them of $103 \pm 26$ km s$^{-1}$. We find that in 97\% of our considered cases, the Clouds experienced a direct collision $147 \pm 33$ Myr ago, with a mean impact parameter of $7.5 \pm 2.5$ kpc. △ Less

Submitted 7 November, 2018; v1 submitted 11 April, 2018; originally announced April 2018.

Comments: 13 pages, 12 figures, 3 tables, accepted to ApJ

Journal ref: The Astrophysical Journal, Volume 864, Issue 1, article id. 55, 13 pp. (2018)

Absolute HST Proper Motion (HSTPROMO) of Distant Milky Way Globular Clusters: Galactocentric Space Velocities and the Milky Way Mass

Authors: Sangmo Tony Sohn, Laura L. Watkins, Mark A. Fardal, Roeland P. van der Marel, Alis J. Deason, Gurtina Besla, Andrea Bellini

Abstract: We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing li… ▽ More We present Hubble Space Telescope (HST) absolute proper motion (PM) measurements for 20 globular clusters (GCs) in the Milky Way (MW) halo at Galactocentric distances $R_{\rm GC} \approx 10-100$ kpc, with median per-coordinate PM uncertainty 0.06 mas yr$^{-1}$. Young and old halo GCs do not show systematic differences in their 3D Galactocentric velocities, derived from combination with existing line-of-sight velocities. We confirm the association of Arp 2, Pal 12, Terzan 7, and Terzan 8 with the Sagittarius (Sgr) stream. These clusters and NGC 6101 have tangential velocity $V_{\rm tan} > 290$ km s$^{-1}$, whereas all other clusters have $V_{\rm tan} < 200$ km s$^{-1}$. NGC 2419, the most distant GC in our sample, is also likely associated with the Sgr stream, whereas NGC 4147, NGC 5024, and NGC 5053 definitely are not. We use the distribution of orbital parameters derived using the 3D velocities to separate halo GCs that either formed within the MW or were accreted. We also assess the specific formation history of e.g. Pyxis and Terzan 8. We constrain the MW mass via an estimator that considers the full 6D phase-space information for 16 of the GCs from $R_{\rm GC} = 10$ to 40 kpc. The velocity dispersion anisotropy parameter $β= 0.609^{+0.130}_{-0.229}$. The enclosed mass $M (<39.5 \rm{kpc}) = 0.61^{+0.18}_{-0.12} \times 10^{12}$ M$_{\odot}$, and the virial mass $M_\rm{vir} = 2.05^{+0.97}_{-0.79} \times 10^{12}$ M$_{\odot}$, are consistent with, but on the high side among recent mass estimates in the literature. △ Less

Submitted 14 June, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 20 pages, 10 figures, accepted for publication in ApJ