Search

We present the method behind HubPUG, a software tool built for recovering systemic proper motions (PMs) of targets (e.g. clusters or resolved galaxies) in 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.

We present analysis of the proper-motion (PM) field of the red clump stars in the Large Magellanic Cloud (LMC) disk using the Gaia Early Data Release 3 catalog. Using a kinematic model based on old stars with 3D velocity measurements, we construct the residual PM field by subtracting the center-of-mass motion and internal rotation motion components. The residual PM field reveals asymmetric patterns, including larger residual PMs in the southern disk. Comparisons between the observed residual PM field with those of five numerical simulations of an LMC analog that is subject to the tidal fields of the Milky Way and the Small Magellanic Cloud (SMC) show that the present-day LMC is not in dynamical equilibrium. We find that both the observed level of disk heating (PM residual root-mean-square of 0.057$\pm$0.002 mas yr$^{-1}$) and kinematic asymmetry are not reproduced by Milky Way tides or if the SMC impact parameter is larger than the size of the LMC disk. This measured level of disk heating provides a novel and important method to validate numerical simulations of the LMC-SMC interaction history. Our results alone put constraints on an impact parameter $\lesssim$10 kpc and impact timing $, Comment: 17 pages, 10 figures, resubmitted to AAS Journal

We use Gaia EDR3 data to identify stars associated with six classical dwarf spheroidals (Draco, Ursa Minor, Sextans, Sculptor, Fornax, Carina) at their outermost radii, beyond their nominal King stellar limiting radius. For all of the dSphs examined, we find radial velocity matches with stars residing beyond the King limiting radius and with $> 50\%$ astrometric probability (four in Draco, two in Ursa Minor, eight in Sextans, two in Sculptor, twelve in Fornax, and five in Carina), indicating that these stars are associated with their respective dwarf spheroidals (dSphs) at high probability. We compare the positions of our candidate "extra-tidal" stars with the orbital tracks of the galaxies, and identify stars, both with and without radial velocity matches, that are consistent with lying along the orbital track of the satellites. However, given the small number of candidate stars, we cannot make any conclusive statements about the significance of these spatially correlated stars. Cross matching with publicly available catalogs of RR Lyrae, we find one RR Lyrae candidate with $> 50\%$ astrometric probability outside the limiting radius in each of Sculptor and Fornax, two such candidates in Draco, nine in Ursa Minor, seven in Sextans, and zero in Carina. Follow-up spectra on all of our candidates, including possible metallicity information, will help confirm association with their respective dSphs, and could represent evidence for extended stellar halos or tidal debris around these classical dSphs., 19 pages, 9 figures, 6 tables, accepted for publication in MNRAS

We combine proper motion data from $Gaia$ EDR3 and HST with line-of-sight velocity data to study the stellar kinematics of the Omega Centauri globular cluster. Using a steady-state, axisymmetric dynamical model, we measure the distribution of both the dark and luminous mass components. Assuming both Gaussian and NFW mass profiles, depending on the dataset, we measure an integrated mass of $\lesssim 10^6$ M$_\odot$ within the Omega Centauri half-light radius for a dark component that is distinct from the luminous stellar component. For the HST and radial velocity data, models with a non-luminous mass component are strongly statistically preferred relative to a stellar mass-only model with a constant mass-to-light ratio. While a compact core of stellar remnants may account for a dynamical mass up to $\sim 5 \times 10^5$ M$_\odot$, they likely cannot explain the higher end of the range. This leaves open the possibility that this non-luminous dynamical mass component is comprised of non-baryonic dark matter. In comparison to the dark matter distributions around dwarf spheroidal galaxies, the Omega Centauri dark mass component is much more centrally concentrated. Interpreting the non-luminous mass distribution as particle dark matter, we use these results to obtain the J-factor, which sets the sensitivity to the annihilation cross section. For the datasets considered, the range of median J-factors is $\sim 10^{22} - 10^{24}$ GeV$^2$ cm$^{-5}$, which is larger than that obtained for any dwarf spheroidal galaxy., 14 pages with 8 figures. Accepted for publication in MNRAS

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., 7 pages, 3 figures, 2 tables. Accepted for publication in ApJL

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., Comment: 18 pages, 11 figures; submitted to ApJ

We present our effort to measure the proper motions of satellites in the halo of the Milky Way with mainly ground based telescopes as a precursor on what is possible with Gaia. For our first study, we used wide field optical data from the LBT combined with a first epoch of SDSS observations, on the globular cluster Palomar 5 (Pal 5). Since Pal 5 is associated with a tidal stream it is very useful to constrain the shape of the potential of the Milky Way. The motion and other properties of the Pal 5 system constrain the inner halo of the Milky Way to be rather spherical. Further, we combined adaptive optics and HST to get an absolute proper motion of the globular cluster Pyxis. Using the proper motion and the line-of-sight velocity we find that the orbit of Pyxis is rather eccentric with its apocenter at more than 100 kpc and its pericenter at about 30 kpc. The dynamics excludes an association with the ATLAS stream, the Magellanic clouds, and all satellites of the Milky Way at least down to the mass of Leo II. However, the properties of Pyxis, like metallicity and age, point to an origin from a dwarf of at least the mass of Leo II. We therefore propose that Pyxis originated from an unknown relatively massive dwarf galaxy, which is likely today fully disrupted. Assuming that Pyxis is bound to the Milky Way we derive a 68% lower limit on the mass of the Milky Way of 9.5 × 1011 M⊙.

We present a new kinematic model for the Small Magellanic Cloud (SMC), using data from the \gaia\ Data Release 2 catalog. We identify a sample of astrometrically well-behaved red giant (RG) stars belonging to the SMC and cross-match with publicly available radial velocity (RV) catalogs. We create a 3D spatial model for the RGs, using RR Lyrae for distance distributions, and apply kinematic models with varying rotation properties and a novel tidal expansion prescription to generate mock proper motion (PM) catalogs. When we compare this series of mock catalogs to the observed RG data, we find a combination of moderate rotation (with a magnitude of $\sim10-20$ km s$^{-1}$ at 1 kpc from the SMC center, inclination between $\sim50-80$ degrees, and a predominantly north-to-south line of nodes position angle of $\sim180$ degrees) and tidal expansion (with a scaling of $\sim10$ km s$^{-1}$ kpc$^{-1}$) is required to explain the PM signatures. The exact best-fit parameters depend somewhat on whether we assess only the PMs or include the RVs as a qualitative check, leaving some small tension remaining between the PM and RV conclusions. In either case, the parameter space preferred by our model is different both from previously inferred rotational geometries, including from the SMC H{\small I} gas and from the RG RV-only analyses, and new SMC PM analyses which conclude that a rotation signature is not detectable. Taken together this underscores the need to treat the SMC as a series of different populations with distinct kinematics., 17 pages, 11 figures, 2 tables, submitted to ApJ

We present the first proper motion measurement for an ultra-faint dwarf spheroidal galaxy, Segue 1, using SDSS and LBC data as the first and second epochs separated by a baseline of $\sim 10$ years. We obtain a motion of $\mu_{\alpha}\,\cos(\delta) = -0.37\pm0.57$ mas yr$^{-1}$ and $\mu_{\delta} =-3.39\pm0.58$ mas yr$^{-1}$. Combining this with the known line-of-sight velocity, this corresponds to a Galactocentric V$_\mathrm{rad}=84\pm9$ and V$_\mathrm{tan}=164^{+66}_{-55}$ km s$^{-1}$. Applying Milky Way halo masses between 0.8 to 1.6$\times 10^{12}$ M$_\odot$ results in an apocenter at 33.9$^{+21.7}_{-7.4}$ kpc and pericenter at 15.4$^{+10.1}_{-9.0}$ kpc from the Galactic center, indicating Segue~1 is rather tightly bound to the Milky Way. Since neither the orbital pole of Segue 1 nor its distance to the Milky Way is similar to the more massive classical dwarfs, it is very unlikely that Segue 1 was once a satellite of a massive known galaxy. Using cosmological zoom-in simulations of Milky Way-mass galaxies, we identify subhalos on similar orbits as Segue~1, which imply the following orbital properties: a median first infall 8.1$^{+3.6}_{-4.3}$ Gyrs ago, a median of 4 pericentric passages since then and a pericenter of 22.8$^{+4.7}_{-4.8}$ kpc. This is slightly larger than the pericenter derived directly from Segue 1 and Milky Way parameters, because galaxies with a small pericenter are more likely to be destroyed. Of the surviving subhalo analogs only 27\% were previously a satellite of a more massive dwarf galaxy (that is now destroyed), thus Segue 1 is more likely to have been accreted on its own., Comment: 16 pages, 7 figures, accepted by ApJ

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., 13 pages, 10 figures, 2 tables, submitted to ApJ, accepted February 8th, 2019

As new work on the proper motions (PMs) of the Large Magellanic Cloud (LMC) has come out, our view of the history of the Magellanic Clouds has evolved. We now believe they are on their first infall into the Milky Way (MW), having been tidally bound at the start of infall (though not necessarily now). Combining these observations with initial PMs of the Small Magellanic Cloud (SMC) suggests a new formation mechanism of the Magellanic Stream through the stripping of material from the SMC. However, large uncertainties remain in the exact mass of the LMC. We present a measurement of the systemic proper motions of the SMC from astrometry with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST), covering a ~3 year baseline of 30 fields with background QSOs. We find these motions to be μW = −0.82 ± 0.06 mas/yr and μN = −1.23 ± 0.07 mas/yr. Combining these measurements with previous efforts in studying the Clouds will help constrain their interactions with each other and the MW, including the mass of the LMC and the MW, as well as provide new insight into the internal dynamics of the SMC.

We use GAIA DR2 proper motions of the RIOTS4 field OB stars in the Small Magellanic Cloud (SMC) to study the kinematics of runaway stars. The data reveal that the SMC Wing has a systemic peculiar motion relative to the SMC Bar of (v_RA, v_Dec) = (62 +/-7, -18+/-5) km/s and relative radial velocity +4.5 +/- 5.0 km/s. This unambiguously demonstrates that these two regions are kinematically distinct: the Wing is moving away from the Bar, and towards the Large Magellanic Cloud with a 3-D velocity of 64 +/- 10 km/s. This is consistent with models for a recent, direct collision between the Clouds. We present transverse velocity distributions for our field OB stars, confirming that unbound runaways comprise on the order of half our sample, possibly more. Using eclipsing binaries and double-lined spectroscopic binaries as tracers of dynamically ejected runaways, and high-mass X-ray binaries (HMXBs) as tracers of runaways accelerated by supernova kicks, we find significant contributions from both populations. The data suggest that HMXBs have lower velocity dispersion relative to dynamically ejected binaries, consistent with the former corresponding to less energetic supernova kicks that failed to unbind the components. Evidence suggests that our fast runaways are dominated by dynamical, rather than supernova, ejections., Comment: Accepted to ApJ Letters. 10 pages, 4 figures

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 $\mu_{W}$ = $-0.82$ $\pm$ 0.02 (random) $\pm$ 0.10 (systematic) mas yr$^{-1}$ and $\mu_{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., Comment: 13 pages, 12 figures, 3 tables, accepted to ApJ

We use properties of void populations identified in $N$-body simulations to forecast the ability of upcoming galaxy surveys to differentiate models of f(R) gravity from \lcdm~cosmology. We analyze multiple simulation realizations, which were designed to mimic the expected number densities, volumes, and redshifts of the upcoming Euclid satellite and a lower-redshift ground-based counterpart survey, using the public {\tt VIDE} toolkit. We examine void abundances, ellipicities, radial density profiles, and radial velocity profiles at redshifts 1.0 and 0.43. We find that stronger f(R) coupling strengths eliminates small voids and produces voids up to $\sim 20\%$ larger in radius, leading to a significant tilt in the void number function. Additionally, under the influence of modified gravity, voids at all scales tend to be measurably emptier with correspondingly higher compensation walls. The velocity profiles reflect this, showing increased outflows inside voids and increased inflows outside voids. Using the void number function as an example, we forecast that future surveys can constrain the modified gravity coupling strength to $\sim 3 \times 10^{-5}$ using voids., 9 pages, 5 figures, accepted by MNRAS

We present VIDE, the Void IDentification and Examination toolkit, an open-source Python/C++ code for finding cosmic voids in galaxy redshift surveys and N-body simulations, characterizing their properties, and providing a platform for more detailed analysis. At its core, VIDE uses a substantially enhanced version of ZOBOV (Neyinck 2008) to calculate a Voronoi tessellation for estimating the density field and a performing a watershed transform to construct voids. Additionally, VIDE provides significant functionality for both pre- and post-processing: for example, vide can work with volume- or magnitude-limited galaxy samples with arbitrary survey geometries, or dark matter particles or halo catalogs in a variety of common formats. It can also randomly subsample inputs and includes a Halo Occupation Distribution model for constructing mock galaxy populations. VIDE uses the watershed levels to place voids in a hierarchical tree, outputs a summary of void properties in plain ASCII, and provides a Python API to perform many analysis tasks, such as loading and manipulating void catalogs and particle members, filtering, plotting, computing clustering statistics, stacking, comparing catalogs, and fitting density profiles. While centered around ZOBOV, the toolkit is designed to be as modular as possible and accommodate other void finders. VIDE has been in development for several years and has already been used to produce a wealth of results, which we summarize in this work to highlight the capabilities of the toolkit. VIDE is publicly available at http://bitbucket.org/cosmicvoids/vide public and http://www.cosmicvoids.net., Comment: 11 pages, 8 figures, minor revisions and expanded discussion from referee comments, Astronomy & Computing accepted, code available at http://bitbucket.org/cosmicvoids/vide_public and http://www.cosmicvoids.net

We use properties of void populations identified in N-body simulations to forecast the ability of upcoming galaxy surveys to differentiate models of f(R) gravity from \Lambda CDM cosmology. We analyze simulations designed to mimic the densities, volumes, and clustering statistics of upcoming surveys, using the public VIDE toolkit. We examine void abundances as a basic probe at redshifts 1.0 and 0.4. We find that stronger f(R) coupling strengths produce voids up to ~20% larger in radius, leading to a significant shift in the void number function. As an initial estimate of the constraining power of voids, we use this change in the number function to forecast a constraint on the coupling strength of $\Delta f_{R_{0}} = 10^{-5}$., Comment: To appear as proceedings of the IAU Symposium 308 "The Zeldovich Universe, Genesis and Growth of the Cosmic Web"