Source author record

Roeland P. van der Marel

Roeland P. van der Marel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2026arXiv

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

We analyze three epochs of HST imaging over 20 years for the Sculptor dwarf spheroidal galaxy, measuring precise proper motions for 119 stars and combining them with 1760 existing line-of-sight velocities. This catalog yields the first radially-resolved 3D velocity dispersion profiles for Sculptor. We confirm mild oblate rotation, with major-axis velocities reaching $\sim 2$ km s$^{-1}$ beyond 20.0 arcmin. Using a methodology similar to that in the first paper in this series, we solve the Jeans equations in oblate axisymmetric geometry to infer the galaxy's mass profile. Our modeling reveals a significant degeneracy due to the unknown galaxy inclination, which is overlooked under spherical symmetry assumptions. This degeneracy allows acceptable fits across a range of dark matter profiles, from cuspy to cored. While we do not directly constrain the inclination with our Jeans models, higher-order line-of-sight velocity moments provide useful additional constraints: comparisons with scalefree models from de Bruijne et al. (1996) favor highly flattened (more face-on) configurations. Adopting an inclination well consistent with these comparisons ($i = 57.1$ degrees), we find, alongside radial velocity anisotropy, a dark matter density slope of $Γ_{\rm dark} = 0.29^{+0.31}_{-0.41}$ within the radial extent of the 3D velocity data, ruling out a cusp with $Γ_{\rm dark} \leq -1$ at 99.8% confidence. This confidence increases for lower inclinations and decreases drastically for nearly edge-on configurations. The results qualitatively agree with $Λ$CDM, SIDM, and Fuzzy DM scenarios that predict core formation, while our specific measurements provide quantitative constraints on the prescriptions of feedback, cross sections, or particle masses required by these models, respectively.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

The recent LMC-SMC collision: Timing and impact parameter constraints from comparison of Gaia LMC disk kinematics and N-body simulations

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 $<$250 Myr. When adopting the impact timing constraint of $\sim$140--160 Myr ago from previous studies, our results suggest that the most recent SMC encounter must have occurred with an impact parameter of $\sim$5 kpc. We also find consistent radial trends in the kinematically- and geometrically-derived disk inclination and line-of-node position angles, indicating a common origin.

preprint2021arXiv

Hunting for intermediate-mass black holes in globular clusters: an astrometric study of NGC 6441

We present an astrometric study of the proper motions (PMs) in the core of the globular cluster NGC 6441. The core of this cluster has a high density and observations with current instrumentation are very challenging. We combine ground-based, high-angular-resolution NACO@VLT images with Hubble Space Telescope ACS/HRC data and measure PMs with a temporal baseline of 15 yr for about 1400 stars in the centermost 15 arcseconds of the cluster. We reach a PM precision of $\sim$30 $μ$as yr$^{-1}$ for bright, well-measured stars. Our results for the velocity dispersion are in good agreement with other studies and extend already-existing analyses of the stellar kinematics of NGC 6441 to its centermost region never probed before. In the innermost arcsecond of the cluster, we measure a velocity dispersion of (19.1 $\pm$ 2.0) km s$^{-1}$ for evolved stars. Because of its high mass, NGC 6441 is a promising candidate for harbouring an intermediate-mass black hole (IMBH). We combine our measurements with additional data from the literature and compute dynamical models of the cluster. We find an upper limit of $M_{\rm IMBH} < 1.32 \times 10^4\,\textrm{M}_\odot$ but we can neither confirm nor rule out its presence. We also refine the dynamical distance of the cluster to $12.74^{+0.16}_{-0.15}$ kpc. Although the hunt for an IMBH in NGC 6441 is not yet concluded, our results show how future observations with extremely-large telescopes will benefit from the long temporal baseline offered by existing high-angular-resolution data.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

SMASHing the low surface brightness SMC

The periphery of the Small Magellanic Cloud (SMC) can unlock important information regarding galaxy formation and evolution in interacting systems. Here, we present a detailed study of the extended stellar structure of the SMC using deep colour-magnitude diagrams (CMDs), obtained as part of the Survey of the MAgellanic Stellar History (SMASH). Special care was taken in the decontamination of our data from MW foreground stars, including from foreground globular clusters NGC 362 and 47 Tuc. We derived the SMC surface brightness using a ``conservative'' approach from which we calculated the general parameters of the SMC, finding a staggered surface brightness profile. We also traced the fainter outskirts by constructing a stellar density profile. This approach, based on stellar counts of the oldest main sequence turn-off (MSTO) stars, uncovered a tidally disrupted stellar feature that reaches as far out as 12 degrees from the SMC centre. We also serendipitously found a faint feature of unknown origin located at $\sim 14$ degrees from the centre of the SMC and that we tentatively associated to a more distant structure. We compared our results to in-house simulations of a $1\times10^{9} M_\odot$ SMC, finding that its elliptical shape can be explained by its tidal disruption under the combined presence of the MW and the LMC. Finally, we found that the older stellar populations show a smooth profile while the younger component presents a jump in the density followed by a flat profile, confirming the heavily disturbed nature of the SMC.

preprint2020arXiv

The Orbital Histories of Magellanic Satellites Using Gaia DR2 Proper Motions

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

preprint2020arXiv

The Second Data Release of the Survey of the MAgellanic Stellar History (SMASH)

The Large and Small Magellanic Clouds (LMC and SMC) are the largest satellite galaxies of the Milky Way and close enough to allow for a detailed exploration of their structure and formation history. The Survey of the MAgellanic Stellar History (SMASH) is a community Dark Energy Camera (DECam) survey of the Magellanic Clouds using $\sim$50 nights to sample over $\sim$2400 deg$^2$ centered on the Clouds at $\sim$20% filling factor (but with contiguous coverage in the central regions) and to depths of $\sim$24th mag in $ugriz$. The primary goals of SMASH are to map out the extended stellar peripheries of the Clouds and uncover their complicated interaction and accretion history as well as to derive spatially-resolved star formation histories of the central regions and create a "movie" of their past star formation. Here we announce the second SMASH public data release (DR2), which contains all 197 fully-calibrated DECam fields including the main body fields in the central regions. The DR2 data are available through the Astro Data Lab hosted by the NSF's National Optical-Infrared Astronomy Research Laboratory. We highlight three science cases that make use of the SMASH DR2 data and will be published in the future: (1) preliminary star formation histories of the LMC; (2) the search for Magellanic star clusters using citizen scientists; and, (3) photometric metallicities of Magellanic Cloud stars using the DECam $u$-band.

preprint2019arXiv

Imprints of Evolution on the Internal Kinematics of Globular Clusters

Globular clusters are collisional systems, meaning that the stars inside them interact on timescales much shorter than the age of the Universe. These frequent interactions transfer energy between stars and set up observable trends that tell the story of a cluster's evolution. This contribution focuses on what we can learn by studying velocity anisotropy and energy equipartition in globular clusters with Hubble Space Telescope proper motions.

preprint2016arXiv

Active compensation of aperture discontinuities for WFIRST-AFTA: analytical and numerical comparison of propagation methods and preliminary results with a WFIRST-AFTA-like pupil

The new frontier in the quest for the highest contrast levels in the focal plane of a coronagraph is now the correction of the large diffractive artifacts effects introduced at the science camera by apertures of increasing complexity. The coronagraph for the WFIRST/AFTA mission will be the first of such instruments in space with a two Deformable Mirrors wavefront control system. Regardless of the control algorithm for these multi Deformable Mirrors, they will have to rely on quick and accurate simulation of the propagation effects introduced by the out-of-pupil surface. In the first part of this paper, we present the analytical description of the different approximations to simulate these propagation effects. In Annex A, we prove analytically that, in the special case of surfaces inducing a converging beam, the Fresnel method yields high fidelity for simulations of these effects. We provide numerical simulations showing this effect. In the second part, we use these tools in the framework of the Active Compensation of Aperture Discontinuities technique (ACAD) applied to pupil geometries similar to WFIRST-AFTA. We present these simulations in the context of the optical layout of the High-contrast imager for Complex Aperture Telescopes, which will test ACAD on a optical bench. The results of this analysis show that using the ACAD method, an apodized pupil lyot coronagraph and the performance of our current deformable mirrors, we are able to obtain, in numerically simulations, a dark hole with an AFTA-like pupil. Our numerical simulation shows that we can obtain contrast better than $2.10^{-9}$ in monochromatic light and better than 3.e-8 with 10% bandwidth between 5 and 14 lambda/D.

preprint2016arXiv

Hubble Space Telescope Proper Motions of Individual Stars in Stellar Streams: Orphan, Sagittarius, Lethe, and the New "Parallel" Stream

We present a multi-epoch Hubble Space Telescope (HST) study of stellar proper motions (PMs) for four fields along the Orphan Stream. We determine absolute PMs of several individual stars per target field using established techniques that utilize distant background galaxies to define a stationary reference frame. Five Orphan Stream stars are identified in one of the four fields based on combined color-magnitude and PM information. The average PM is consistent with the existing model of the Orphan stream by Newberg et al. In addition to the Orphan stream stars, we detect stars that likely belong to other stellar streams. To identify which stellar streams these stars belong to, we examine the 2-d bulk motion of each group of stars on the sky by subtracting the PM contribution of the solar motion (which is a function of position on the sky and distance) from the observed PMs, and comparing the vector of net motion with the spatial extent of known stellar streams. By doing this, we identify candidate stars in the Sagittarius and Lethe streams, and a newly-found stellar stream at a distance of ~17 kpc, which we tentatively name the "Parallel stream". Together with our Sagittarius stream study (Sohn et al., 2015, ApJ, 803, 56), this work demonstrates that even in the Gaia era, HST will continue to be advantageous in measuring PMs of old stellar populations on a star-by-star basis, especially for distances beyond ~10 kpc.

preprint2016arXiv

Isotropic at the Break? 3D Kinematics of Milky Way Halo Stars in the Foreground of M31

We present the line-of-sight (LOS) velocities for 13 distant main sequence Milky Way halo stars with published proper motions. The proper motions were measured using long baseline (5-7 years) multi-epoch HST/ACS photometry, and the LOS velocities were extracted from deep (5-6 hour integrations) Keck II/DEIMOS spectra. We estimate the parameters of the velocity ellipsoid of the stellar halo using a Markov chain Monte Carlo ensembler sampler method. The velocity second moments in the directions of the Galactic $(l,b,$ LOS) coordinate system are $\langle v^2_l \rangle^{1/2} = 138^{+43}_{-26}$ km/s, $\langle v^2_b \rangle^{1/2} = 88^{+28}_{-17}$ km/s, and $\langle v^2_{\rm{LOS}} \rangle^{1/2} = 91^{+27}_{-14}$ km/s. We use these ellipsoid parameters to constrain the velocity anisotropy of the stellar halo. Ours is the first measurement of the anisotropy parameter $β$ using 3D kinematics outside of the solar neighborhood. We find $β=-0.3^{+0.4}_{-0.9}$, consistent with isotropy and lower than solar neighborhood $β$ measurements by 2$σ$ ($β_{SN} \sim 0.5-0.7$). We identify two stars in our sample that are likely members of the known TriAnd substructure, and excluding these objects from our sample increases our estimate of the anisotropy to $β=0.1^{+0.4}_{-1.0}$, which is still lower than solar neighborhood measurements by $1σ$. The potential decrease in $β$ with Galactocentric radius is inconsistent with theoretical predictions, though consistent with recent observational studies, and may indicate the presence of large, shell-type structure (or structures) at $r \sim 25$ kpc. The methods described in this paper will be applied to a much larger sample of stars with 3D kinematics observed through the ongoing HALO7D program.

preprint2016arXiv

Low Surface Brightness Imaging of the Magellanic System: Imprints of Tidal Interactions between the Clouds in the Stellar Periphery

We present deep optical images of the Large and Small Magellanic Clouds (LMC and SMC) using a low cost telephoto lens with a wide field of view to explore stellar substructure in the outskirts of the stellar disk of the LMC (r < 10 degrees from the center). These data have higher resolution than existing star count maps, and highlight the existence of stellar arcs and multiple spiral arms in the northern periphery, with no comparable counterparts in the South. We compare these data to detailed simulations of the LMC disk outskirts, following interactions with its low mass companion, the SMC. We consider interaction in isolation and with the inclusion of the Milky Way tidal field. The simulations are used to assess the origin of the northern structures, including also the low density stellar arc recently identified in the DES data by Mackey et al. 2015 at ~ 15 degrees. We conclude that repeated close interactions with the SMC are primarily responsible for the asymmetric stellar structures seen in the periphery of the LMC. The orientation and density of these arcs can be used to constrain the LMC's interaction history with and impact parameter of the SMC. More generally, we find that such asymmetric structures should be ubiquitous about pairs of dwarfs and can persist for 1-2 Gyr even after the secondary merges entirely with the primary. As such, the lack of a companion around a Magellanic Irregular does not disprove the hypothesis that their asymmetric structures are driven by dwarf-dwarf interactions.

preprint2015arXiv

Hubble Space Telescope proper motion (HSTPROMO) catalogs of Galactic globular clusters. II. Kinematic profiles and maps

We present kinematical analyses of 22 Galactic globular clusters using the Hubble Space Telescope proper motion (HSTPROMO) catalogues recently presented in Bellini et al. (2014). For most clusters, this is the first proper-motion study ever performed, and, for many, this is the most detailed kinematic study of any kind. We use cleaned samples of bright stars to determine binned velocity-dispersion and velocity-anisotropy radial profiles and two-dimensional velocity-dispersion spatial maps. Using these profiles, we search for correlations between cluster kinematics and structural properties. We find that: (1) more centrally-concentrated clusters have steeper radial velocity-dispersion profiles; (2) on average, at 1σconfidence in two dimensions, the photometric and kinematic centres of globular clusters agree to within ~1", with a cluster-to-cluster rms of 4" (including observational uncertainties); (3) on average, the cores of globular clusters have isotropic velocity distributions to within 1% (σ_t/σ_r = 0.992 +/- 0.005), with a cluster-to-cluster rms of 2% (including observational uncertainties); (4) clusters generally have mildly radially anisotropic velocity distributions (σ_t/σ_r ~ 0.8-1.0) near the half-mass radius, with bigger deviations from isotropy for clusters with longer relaxation times; (5) there is a relation between σ_minor/σ_major and ellipticity, such that the more flattened clusters in the sample tend to be more anisotropic, with σ_minor/σ_major ~ 0.9-1.0. Aside from these general results and correlations, the profiles and maps presented here can provide a basis for detailed dynamical modelling of individual globular clusters. Given the quality of the data, this is likely to provide new insights into a range of topics concerning globular cluster mass profiles, structure, and dynamics.

preprint2015arXiv

Hubble Space Telescope proper motion (HSTPROMO) catalogs of Galactic globular clusters. III. Dynamical distances and mass-to-light ratios

We present dynamical distance estimates for 15 Galactic globular clusters and use these to check the consistency of dynamical and photometric distance estimates. For most of the clusters, this is the first dynamical distance estimate ever determined. We extract proper-motion dispersion profiles using cleaned samples of bright stars from the Hubble Space Telescope proper-motion catalogs recently presented in Bellini et al. (2014) and compile a set of line-of-sight velocity-dispersion profiles from a variety of literature sources. Distances are then estimated by fitting spherical, non-rotating, isotropic, constant mass-to-light (M/L) dynamical models to the proper-motion and line-of-sight dispersion profiles together. We compare our dynamical distance estimates with literature photometric estimates from the Harris (1996, 2010 edition) globular cluster catalog and find that the mean fractional difference between the two types is consistent with zero at just $-1.9 \pm 1.7 \%$. This indicates that there are no significant biases in either estimation method and provides an important validation of the stellar-evolution theory that underlies photometric distance estimates. The analysis also estimates dynamical M/L ratios for our clusters; on average, the dynamically-inferred M/L ratios agree with existing stellar-population-based M/L ratios that assume a Chabrier initial mass function (IMF) to within $-8.8 \pm 6.4 \%$, implying that such an IMF is consistent with our data. Our results are also consistent with a Kroupa IMF, but strongly rule out a Salpeter IMF. We detect no correlation between our M/L offsets from literature values and our distance offsets from literature values, strongly indicating that our methods are reliable and our results are robust.

preprint2015arXiv

Hubble Space Telescope Proper Motions along the Sagittarius Stream: I. Observations and Results for Stars in Four Fields

We present a multi-epoch Hubble Space Telescope (HST) study of stellar proper motions (PMs) for four fields spanning 200 degrees along the Sagittarius (Sgr) stream: one trailing arm field, one field near the Sgr dwarf spheroidal tidal radius, and two leading arm fields. We determine absolute PMs of dozens of individual stars per field, using established techniques that use distant background galaxies as stationary reference frame. Stream stars are identified based on combined color-magnitude diagram and PM information. The results are broadly consistent with the few existing PM measurements for the Sgr galaxy and the trailing arm. However, our new results provide the highest PM accuracy for the stream to date, the first PM measurements for the leading arm, and the first PM measurements for individual stream stars; we also serendipitously determine the PM of the globular cluster NGC~6652. In the trailing-arm field, the individual PMs allow us to kinematically separate trailing-arm stars from leading-arm stars that are 360 degrees further ahead in their orbit. Also, in three of our fields we find indications that two distinct kinematical components may exist within the same arm and wrap of the stream. Qualitative comparison of the HST data to the predictions of the Law & Majewski and Penarrubia et al. N-body models show that the PM measurements closely follow the predicted trend with Sgr longitude. This provides a successful consistency check on the PM measurements, as well as on these N-body approaches (which were not tailored to fit any PM data).

preprint2015arXiv

Hydra II: a faint and compact Milky Way dwarf galaxy found in the Survey of the Magellanic Stellar History

We present the discovery of a new dwarf galaxy, Hydra II, found serendipitously within the data from the ongoing Survey of the MAgellanic Stellar History (SMASH) conducted with the Dark Energy Camera on the Blanco 4m Telescope. The new satellite is compact (r_h = 68 +/- 11 pc) and faint (M_V = -4.8 +/- 0.3), but well within the realm of dwarf galaxies. The stellar distribution of HydraII in the color-magnitude diagram is well-described by a metal-poor ([Fe/H] = -2.2) and old (13 Gyr) isochrone and shows a distinct blue horizontal branch, some possible red clump stars, and faint stars that are suggestive of blue stragglers. At a heliocentric distance of 134 +/- 10 kpc, Hydra II is located in a region of the Galactic halo that models have suggested may host material from the leading arm of the Magellanic Stream. A comparison with N-body simulations hints that the new dwarf galaxy could be or could have been a satellite of the Magellanic Clouds.

preprint2015arXiv

Ram Pressure Stripping of the Large Magellanic Cloud's Disk as a Probe of the Milky Way's Circumgalactic Medium

Recent observations have constrained the orbit and structure of the Large Magellanic Cloud (LMC), implying a well-constrained pericentric passage about the Milky Way (MW) ~ 50 Myr ago. In this scenario, the LMC's gaseous disk has recently experienced maximal ram pressure stripping, suggesting the current extent of its HI disk directly probes the medium in which it is moving. From the observed stellar and HI distributions of the system we find evidence of a truncated gas profile along the windward ``leading edge' of the LMC disk, despite a far more extended stellar component. We explore the implications of this ram pressure stripping signature, using both analytic prescriptions and full 3-dimensional hydrodynamic simulations of the LMC. Our simulations subject the system to a headwind whose velocity components correspond directly to the recent orbital history of the LMC. We vary the density of this headwind, using a variety of sampled parameters for a Beta-profile for a theoretical MW circumgalactic medium (CGM), comparing the resulting HI morphology directly to observations of the LMC HI and stellar components. This model can match the radial extent of the LMC's leading (windward) edge only in scenarios where the MW CGM density at pericentric passage is n(R = 48.2 +/- 5 kpc) = 1.1 (+.44/-.45) x 1e-4 cm^-3. The implied pericentric density proves insensitive to both the broader CGM structure and temperature profile, thus providing a model-independent constraint on the local gas density. This result imposes an important constraint on the density profile of the MW's CGM, and thus the total baryon content of the MW. From our work, assuming a Beta-profile valid to ~ Rvir, we infer a total diffuse CGM mass M(300 kpc) = 2.6 +/- 1.4 x 1e10 Msun or approximately 15% of a 1e12 Msun MW's baryonic mass budget.

preprint2014arXiv

Local Group Proper Motion Dynamics

Our knowledge of the dynamics and masses of galaxies in the Local Group has long been limited by the fact that only line-of-sight velocities were observationally accessible. This introduces significant degeneracies in dynamical models, which can only be resolved by measuring also the velocity components perpendicular to the line of sight. However, beyond the solar neighborhood, the corresponding proper motions have generally been too small to measure. This has changed dramatically over the past decade, especially due to the angular resolution and stability available on the Hubble Space Telescope. Proper motions can now be reliably measured throughout the Local Group, as illustrated by, e.g., the work of the HSTPROMO collaboration. In this review, I summarize the importance of proper motions for Local Group science, and I describe the current and future observational approaches and facilities available to measure proper motions. I highlight recent results on various Milky Way populations (globular clusters, the bulge, the metal-poor halo, hypervelocity stars, and tidal streams), dwarf satellite galaxies, the Magellanic Clouds and the Andromeda System.

preprint2013arXiv

HST absolute Proper Motions of NGC 6681 (M70) and the Sagittarius Dwarf Spheroidal Galaxy

We have measured absolute proper motions for the three populations intercepted in the direction of the Galactic globular cluster NGC 6681: the cluster itself, the Sagittarius dwarf spheroidal galaxy and the field. For this we used Hubble Space Telescope ACS/WFC and WFC3/UVIS optical imaging data separated by a temporal baseline of 5.464 years. Five background galaxies were used to determine the zero point of the absolute-motion reference frame. The resulting absolute proper motion of NGC 6681 is ($μ_α\cosδ, μ_δ$)=($1.58\pm0.18, -4.57\pm0.16$) \masyr. This is the first estimate ever made for this cluster. For the Sgr dSph we obtain ($μ_α\cosδ, μ_δ)=(-2.54\pm0.18, -1.19\pm0.16$) \masyr, consistent with previous measurements and with the values predicted by theoretical models. The absolute proper motion of the Galaxy population in our field of view is ($μ_α\cosδ, μ_δ)=(-1.21\pm0.27, -4.39\pm0.26$) \masyr. In this study we also use background Sagittarius Dwarf Spheroidal stars to determine the rotation of the globular cluster in the plane of the sky and find that NGC 6681 is not rotating significantly:\ $v_{\rm rot}=0.82\pm1.02$ km$\,$s$^{-1}$ at a distance of 1 arcmin from the cluster center.

preprint2013arXiv

HST Proper Motions of Stars within Globular Clusters

The stable environment of space makes HST an excellent astrometric tool. Its diffraction-limited resolution allows it to distinguish and measure positions and fluxes for stars all the way to the center of most globular clusters. Apart from small changes due to breathing, its PSFs and geometric distortion have been extremely stable over its 20-year lifetime. There are now over 20 globular clusters for which there exist two or more well-separated epochs in the archive, spanning up to 10+ years. Our photometric and astrometric techniques have allowed us to measure tens of thousands of stars per cluster within one arcmin from the center, with typical proper-motion errors of ~0.02 mas/yr, which translates to ~0.8 km/s for a typical cluster. These high-quality measurements can be used to detect the possible presence of a central intermediate-mass black hole, and put constraints on its mass. In addition, they will provide a direct measurement of the cluster anisotropy and equipartition. We present preliminary results from this project, and discuss them in the context of what is already known from other techniques

preprint2013arXiv

Local Group and Star Cluster Dynamics from HSTPROMO (The Hubble Space Telescope Proper Motion Collaboration)

The Hubble Space Telescope (HST) has proven to be uniquely suited for the measurement of proper motions (PMs) of stars and galaxies in the nearby Universe. Here we summarize the main results and ongoing studies of the HSTPROMO collaboration, which over the past decade has executed some two dozen observational and theoretical HST projects on this topic. This is continuing to revolutionize our dynamical understanding of many objects, including: globular clusters; young star clusters; stars and stellar streams in the Milky Way halo; Local Group galaxies, including dwarf satellite galaxies, the Magellanic Clouds, and the Andromeda galaxy; and AGN Black Hole Jets.

preprint2013arXiv

Optical Proper Motion Measurements of the M87 Jet: New Results from the Hubble Space Telescope

We report new results from an HST archival program to study proper motions in the optical jet of the nearby radio galaxy M87. Using over 13 years of archival imaging, we reach accuracies below 0.1c in measuring the apparent velocities of individual knots in the jet. We confirm previous findings of speeds up to 4.5c in the inner 6" of the jet, and report new speeds for optical components in the outer part of the jet. We find evidence of significant motion transverse to the jet axis on the order of 0.6c in the inner jet features, and superluminal velocities parallel and transverse to the jet in the outer knot components, with an apparent ordering of velocity vectors possibly consistent with a helical jet pattern. Previous results suggested a global deceleration over the length of the jet in the form of decreasing maximum speeds of knot components from HST-1 outward, but our results suggest that superluminal speeds persist out to knot C, with large differentials in very nearby features all along the jet. We find significant apparent accelerations in directions parallel and transverse to the jet axis, along with evidence for stationary features in knots D, E, and I. These results are expected to place important constraints on detailed models of kpc-scale relativistic jets.

preprint2013arXiv

The Space Motion of Leo I: Hubble Space Telescope Proper Motion and Implied Orbit

We present the first absolute proper motion measurement of Leo I, based on two epochs of HST ACS/WFC images separated by ~5 years. The average shift of Leo I stars with respect to ~100 background galaxies implies a proper motion of (mu_W, mu_N) = (0.1140 +/- 0.0295, -0.1256 +/- 0.0293) mas/yr. The implied Galactocentric velocity vector, corrected for the reflex motion of the Sun, has radial and tangential components V_rad = 167.9 +/- 2.8 km/s and V_tan = 101.0 +/- 34.4 km/s, respectively. We study the detailed orbital history of Leo I by solving its equations of motion backward in time for a range of plausible mass models for the Milky Way and its surrounding galaxies. Leo I entered the Milky Way virial radius 2.33 +/- 0.21 Gyr ago, most likely on its first infall. It had a pericentric approach 1.05 +/- 0.09 Gyr ago at a Galactocentric distance of 91 +/- 36 kpc. We associate these time scales with characteristic time scales in Leo I's star formation history, which shows an enhanced star formation activity ~2 Gyr ago and quenching ~1 Gyr ago. There is no indication from our calculations that other galaxies have significantly influenced Leo I's orbit, although there is a small probability that it may have interacted with either Ursa Minor or Leo II within the last ~1 Gyr. For most plausible Milky Way masses, the observed velocity implies that Leo I is bound to the Milky Way. However, it may not be appropriate to include it in models of the Milky Way satellite population that assume dynamical equilibrium, given its recent infall. Solution of the complete (non-radial) timing equations for the Leo I orbit implies a Milky Way mass M_MW,vir = 3.15 (-1.36, +1.58) x 10^12 Msun, with the large uncertainty dominated by cosmic scatter. In a companion paper, we compare the new observations to the properties of Leo I subhalo analogs extracted from cosmological simulations.

preprint2013arXiv

The Space Motion of Leo I: The Mass of the Milky Way's Dark Matter Halo

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M_MW). Despite Leo I's large Galacto-centric space velocity (200 km/s) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M_MW > 10^{12} M_sun at 95% confidence for a variety of Bayesian priors on M_MW. In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M_MW would increase by 30%. Imposing a mass weighted LCDM prior, we find a median Milky Way virial mass of M_MW=1.6 x 10^{12} M_sun, with a 90% confidence interval of [1.0-2.4] x 10^{12} M_sun. We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.

preprint2013arXiv

Third-Epoch Magellanic Cloud Proper Motions I: HST/WFC3 data and Orbit Implications

We present proper motions for the Large & Small Magellanic Clouds (LMC & SMC) based on three epochs of \textit{Hubble Space Telescope} data, spanning a $\sim 7$ yr baseline, and centered on fields with background QSOs. The first two epochs, the subject of past analyses, were obtained with ACS/HRC, and have been reanalyzed here. The new third epoch with WFC3/UVIS increases the time baseline and provides better control of systematics. The three-epoch data yield proper motion random errors of only 1-2% per field. For the LMC this is sufficient to constrain the internal proper motion dynamics, as will be discussed in a separate paper. Here we focus on the implied center-of-mass proper motions: mu_W(LMC) = -1.910 +/- 0.020 mas/yr, mu_N(LMC) = 0.229 +/- 0.047 mas/yr, and mu_W(SMC) = -0.772 +/- 0.063 mas/yr, mu_N(SMC) = -1.117 +/- 0.061 mas/yr. We combine the results with a revised understanding of the solar motion in the Milky Way to derive Galactocentric velocities: v_{tot,LMC} = 321 +/- 24 km/s and v_{tot,SMC} = 217 +/- 26 km/s. Our proper motion uncertainties are now dominated by limitations in our understanding of the internal kinematics and geometry of the Clouds, and our velocity uncertainties are dominated by distance errors. Orbit calculations for the Clouds around the Milky Way allow a range of orbital periods, depending on the uncertain masses of the Milky Way and LMC. Periods $\lesssim 4$ Gyr are ruled out, which poses a challenge for traditional Magellanic Stream models. First-infall orbits are preferred (as supported by other arguments as well) if one imposes the requirement that the LMC and SMC must have been a bound pair for at least several Gyr.

preprint2013arXiv

Third-Epoch Magellanic Cloud Proper Motions II: The Large Magellanic Cloud Rotation Field in Three Dimensions

We present the first detailed assessment of the large-scale rotation of any galaxy based on full three-dimensional velocity measurements. We do this for the LMC by combining our HST average proper motion (PM) measurements for stars in 22 fields, with existing line-of-sight (LOS) velocity measurements for 6790 individual stars. We interpret these data with a model of circular rotation in a flat disk. The PM and LOS data paint a consistent picture of the LMC rotation and their combination yields several new insights. The PM data imply a stellar dynamical center that coincides with the HI dynamical center, and a rotation curve amplitude consistent with that inferred from LOS velocity studies. The implied disk viewing angles agree with the range of values found in the literature, but continue to indicate variations with stellar population and/or radius. Young (RSG) stars rotate faster than old (RGB/AGB) stars due to asymmetric drift. Outside the central region, the circular velocity is approximately flat at Vcirc = 91.7 +/- 18.8 km/s. This is consistent with the baryonic Tully-Fisher relation, and implies an enclosed mass M(8.7 kpc) = (1.7 +/- 0.7) x 10^10 solar masses. The virial mass is larger and depends on the full extent of the dark halo. The tidal radius is 22.3 +/- 5.2 kpc (24.0 +/- 5.6 degrees). Combination of the PM and LOS data yields kinematic distance estimates for the LMC, but these are not yet competitive with other methods.

preprint2012arXiv

Chandra Observations of the Nuclear Star Cluster and Ultraluminous X-ray Sources in NGC 2139

We report Chandra observations of the Scd galaxy NGC 2139, which is known to host a recently formed (10^7.6 yrs) nuclear star cluster. The star cluster is undetected in X-rays, with an upper bound on 0.5-7 keV luminosity of L_X < 7.1 x 10^37 erg/s. This bound implies a bolometric accretion luminosity <0.3 percent of the Eddington luminosity for a black hole with the mass (approximately 3400 M_sun) expected from extrapolation of the M-sigma relation. The lack of X-ray emission indicates that a black hole, if present, is not undergoing significant accretion at the current time. While the central cluster is undetected, the data reveal a substantial population of bright X-ray point sources elsewhere in this galaxy, with eight qualifying as ultraluminous X-ray sources with L_X > 10^39 erg/s. We use archival Hubble Space Telescope images to identify candidate optical counterparts for seven Chandra sources, which in most cases have optical luminosities and spatial profiles consistent with star clusters. Compared with other galaxies, the number of luminous X-ray sources in NGC 2139 is larger by a factor of 4 - 10 than expected based on its present star formation rate and stellar mass. This finding can be understood if NGC 2139 has concluded a burst of star formation in the recent past, and suggests that this galaxy could be important for testing the use of X-ray source populations as a chronometer of star formation history.

preprint2012arXiv

Dwarfs Gobbling Dwarfs: A Stellar Tidal Stream Around NGC 4449 and Hierarchical Galaxy Formation on Small Scales

A candidate diffuse stellar substructure was previously reported in the halo of the nearby dwarf starburst galaxy NGC 4449 by Karachentsev et al. We map and analyze this feature using a unique combination of deep integrated-light images from the Black Bird 0.5-meter telescope, and high-resolution wide-field images from the 8-meter Subaru telescope, which resolve the nebulosity into a stream of red giant branch stars, and confirm its physical association with NGC 4449. The properties of the stream imply a massive dwarf spheroidal progenitor, which after complete disruption will deposit an amount of stellar mass that is comparable to the existing stellar halo of the main galaxy. The ratio between luminosity or stellar-mass between the two galaxies is ~1:50, while the indirectly measured dynamical mass-ratio, when including dark matter, may be ~1:10-1:5. This system may thus represent a "stealth" merger, where an infalling satellite galaxy is nearly undetectable by conventional means, yet has a substantial dynamical influence on its host galaxy. This singular discovery also suggests that satellite accretion can play a significant role in building up the stellar halos of low-mass galaxies, and possibly in triggering their starbursts.

preprint2012arXiv

HST/ACS Photometry of Old Stars in NGC 1569: The Star Formation History of a Nearby Starburst

(abridged) We used HST/ACS to obtain deep V- and I-band images of NGC 1569, one of the closest and strongest starburst galaxies in the Universe. These data allowed us to study the underlying old stellar population, aimed at understanding NGC 1569's evolution over a full Hubble time. We focus on the less-crowded outer region of the galaxy, for which the color-magnitude diagram (CMD) shows predominantly a red giant branch (RGB) that reaches down to the red clump/horizontal branch feature (RC/HB). A simple stellar population analysis gives clear evidence for a more complicated star formation history (SFH) in the outer region. We derive the full SFH using a newly developed code, SFHMATRIX, which fits the CMD Hess diagram by solving a non-negative least squares problem. Our analysis shows that the relative brightnesses of the RGB tip and RC/HB, along with the curvature and color of the RGB, provide enough information to ameliorate the age-metallicity-extinction degeneracy. The distance/reddening combination that best fits the data is E(B-V) = 0.58 +/- 0.03 and D = 3.06 +/- 0.18 Mpc. Star formation began ~ 13 Gyr ago, and this accounts for the majority of the mass in the outer region. However, the initial burst was followed by a relatively low, but constant, rate of star formation until ~ 0.5-0.7 Gyr ago when there may have been a short, low intensity burst of star formation.

preprint2012arXiv

The M31 Velocity Vector. I. Hubble Space Telescope Proper Motion Measurements

We present the first proper motion measurements for the galaxy M31. We obtained new V-band imaging data with the HST ACS/WFC and WFC3/UVIS of a spheroid field near the minor axis, an outer disk field along the major axis, and a field on the Giant Southern Stream. The data provide 5-7 year time baselines with respect to pre-existing deep first-epoch observations. We measure the positions of thousands of M31 stars and hundreds of compact background galaxies in each field. High accuracy and robustness is achieved by building and fitting a unique template for each individual object. The average proper motion for each field is obtained from the average motion of the M31 stars between the epochs with respect to the background galaxies. For the three fields, the observed proper motions (mu_W,mu_N) are (-0.0458, -0.0376), (-0.0533, -0.0104), and (-0.0179,-0.0357) mas/yr, respectively. The ability to average over large numbers of objects and over the three fields yields a final accuracy of 0.012 mas/yr. The robustness of the proper-motion measurements and uncertainties are supported by the fact that data from different instruments, taken at different times and with different telescope orientations, as well as measurements of different fields, all yield statistically consistent results. Papers II and III explore the implications for our understanding of the history, future, and mass of the Local Group. (Abridged)

preprint2012arXiv

The Role of Dwarf Galaxy Interactions in Shaping the Magellanic System and Implications for Magellanic Irregulars

We present a novel pair of numerical models of the interaction history between the Large and Small Magellanic Clouds (LMC and SMC, respectively) and our Milky Way (MW) in light of recent high precision proper motions (Kallivayalil et al. 2006a,b). Given the new velocities, cosmological simulations of structure formation favor a scenario where the Magellanic Clouds (MCs) are currently on their first infall towards our Galaxy (Boylan-Kolchin et al. 2011, Busha et al. 2011). We illustrate here that the observed irregular morphology and internal kinematics of the MCs (in gas and stars) are naturally explained by interactions between the LMC and SMC, rather than gravitational interactions with the MW. This picture further supports a first infall scenario (Besla et a. 2007). In particular, we demonstrate that the Magellanic Stream, a band of HI gas trailing behind the MCs 150 degrees across the sky, can be accounted for by the action of LMC tides on the SMC before the system was accreted by the MW. We further demonstrate that the off-center, warped stellar bar of the LMC and its one-armed spiral, can be naturally explained by a recent direct collision with the SMC. Such structures are key morphological characteristics of a class of galaxies referred to as Magellanic Irregulars (de Vaucouleurs & Freeman 1972), the majority of which are not associated with massive spiral galaxies. We infer that dwarf-dwarf galaxy interactions are important drivers for the morphological evolution of Magellanic Irregulars and can dramatically affect the efficiency of baryon removal from dwarf galaxies via the formation of extended tidal bridges and tails. Such interactions are important not only for the evolution of dwarf galaxies but also have direct consequences for the buildup of baryons in our own MW, as LMC-mass systems are believed to be the dominant building blocks of MW-type halos.

preprint2011arXiv

Cluster 77 in NGC 4449: The nucleus of a satellite galaxy being transformed into a globular cluster?

We report the discovery in our HST ACS B, V, and I images of NGC 4449 of a globular cluster (GC) which appears associated with two tails of blue stars. The cluster is massive (M~1.7x10^6 M_sun) and highly flattened (ε~0.24). From the color-magnitude diagrams of the resolved stars we infer active star formation in the tails over the past ~200 Myr. In a diagram of mean projected mass density inside r_e versus total mass the cluster lies at the upper end of the GC distribution, where galaxy nuclei are. The north-west tail is associated with a concentration of HI and infrared (dust/PAHs) emission which appears as part of a much longer stream wrapping around the galaxy. These properties suggest that the cluster may be the nucleus of a former gas-rich satellite galaxy undergoing tidal disruption by NGC 4449. If so, the cluster is seen in an earlier phase compared to other suggested nuclei of disrupted galaxies such as ωCen (Milky Way) and G1 (M31).

preprint2011arXiv

Hubble Space Telescope study of resolved red giant stars in the outer halos of nearby dwarf starburst galaxies

[abridged] Aims. We observed the outer parts of NGC 1569 and NGC 4449, two of the closest and strongest dwarf starburst galaxies in the local universe, to characterize their stellar density and populations, and obtain new insights into the structure, formation, and evolution of starburst galaxies and galaxy halos. Methods. We obtained HST/WFPC2 images between 5 and 8 scale radii from the center, along the intermediate and minor axes. We performed point-source photometry to determine color magnitude diagrams of I vs. V-I. We compared the results at different radii, including also our prior HST/ACS results for more centrally located fields. Results. We detect stars in the RGB and TP-AGB (carbon star) phases in all outer fields, but not younger stars such as those present at smaller radii. The RGB star density profile is well fit by either a de Vaucouleurs profile or a power-law profile, but has more stars at large radii than a single exponential. To within the uncertainties, there are no radial gradients in the RGB color or carbon-to-RGB-star ratio at large radii. Conclusions. The galaxies have faint outer stellar envelopes that are not tidally truncated within the range of radii addressed by our study. The density profiles suggest that these are not outward extensions of the inner disks, but are instead distinct stellar halos. This agrees with other work on galaxies of similar morphology. The presence of such halos is consistent with predictions of hierarchical galaxy formation scenarios. The halos consist of intermediate-age/old stars, confirming the results of other studies that have shown the starburst phenomenon to be very centrally concentrated. There is no evidence for stellar-population age/metallicity gradients within the halos themselves.

preprint2011arXiv

The cluster population of the irregular galaxy NGC 4449 as seen by the Hubble Advanced Camera for Surveys

We present a study of the star cluster population in the starburst irregular galaxy NGC 4449 based on B, V, I, and Ha images taken with the Advanced Camera for Surveys on the Hubble Space Telescope. We derive the cluster properties such as size, ellipticity, and total magnitudes. Cluster ages and masses are derived fitting the observed spectral energy distributions with different population synthesis models. Our analysis is strongly affected by the age-metallicity degeneracy; however, if we assume a metallicity of ~1/4 solar, as derived from spectroscopy of HII regions, we find that the clusters have ages distributed quite continuously over a Hubble time, and they have masses from ~10^3 M_sun up to ~2 x 10^6 M_sun, assuming a Salpeters' IMF down to 0.1 M_sun. Young clusters are preferentially located in regions of young star formation, while old clusters are distributed over the whole NGC 4449 field of view, like the old stars (although we notice that some old clusters follow linear structures, possibly a reflection of past satellite accretion). The high SF activity in NGC 4449 is confirmed by its specific frequency of young massive clusters, higher than the average value found in nearby spirals and in the LMC (but lower than in other starburst dwarfs such as NGC 1705 and NGC 1569), and by the flat slope of the cluster luminosity function (dN(L_V)\propto L_V^{-1.5} dL for clusters younger than 1 Gyr). We use the upper envelope of the cluster log(mass) versus log(age) distribution to quantify cluster disruption, and do not find evidence for the high (90%) long-term infant mortality found by some studies. For the red clusters, we find correlations between size, ellipticity, luminosity and mass: brighter and more massive clusters tend to be more compact, and brighter clusters tend to be also more elliptical.

preprint2010arXiv

Simulations of the Magellanic Stream in a First Infall Scenario

Recent high precision proper motions from the Hubble Space Telescope (HST) suggest that the Large and Small Magellanic Clouds (LMC and SMC, respectively) are either on their first passage or on an eccentric long period (>6 Gyr) orbit about the Milky Way (MW). This differs markedly from the canonical picture in which the Clouds travel on a quasi-periodic orbit about the MW (period of ~2 Gyr). Without a short period orbit about the MW, the origin of the Magellanic Stream, a young (1-2 Gyr old) coherent stream of HI gas that trails the Clouds ~150 degrees across the sky, can no longer be attributed to stripping by MW tides and/or ram pressure stripping by MW halo gas. We propose an alternative formation mechanism in which material is removed by LMC tides acting on the SMC before the system is accreted by the MW. We demonstrate the feasibility and generality of this scenario using an N-body/SPH simulation with cosmologically motivated initial conditions constrained by the observations. Under these conditions we demonstrate that it is possible to explain the origin of the Magellanic Stream in a first infall scenario. This picture is generically applicable to any gas-rich dwarf galaxy pair infalling towards a massive host or interacting in isolation.

preprint2009arXiv

Mass Segregation in NGC 2298: limits on the presence of an Intermediate Mass Black Hole

[abridged] Theoretical investigations have suggested the presence of Intermediate Mass Black Holes (IMBHs, with masses in the 100-10000 Msun range) in the cores of some Globular Clusters (GCs). In this paper we present the first application of a new technique to determine the presence or absence of a central IMBH in globular clusters that have reached energy equipartition via two-body relaxation. The method is based on the measurement of the radial profile for the average mass of stars in the system, using the fact that a quenching of mass segregation is expected when an IMBH is present. Here we measure the radial profile of mass segregation using main-sequence stars for the globular cluster NGC 2298 from resolved source photometry based on HST-ACS data. The observations are compared to expectations from direct N-body simulations of the dynamics of star clusters with and without an IMBH. The mass segregation profile for NGC 2298 is quantitatively matched to that inferred from simulations without a central massive object over all the radial range probed by the observations, that is from the center to about two half-mass radii. Profiles from simulations containing an IMBH more massive than ~ 300-500 Msun (depending on the assumed total mass of NGC 2298) are instead inconsistent with the data at about 3 sigma confidence, irrespective of the IMF and binary fraction chosen for these runs. While providing a null result in the quest of detecting a central black hole in globular clusters, the data-model comparison carried out here demonstrates the feasibility of the method which can also be applied to other globular clusters with resolved photometry in their cores.

preprint2009arXiv

New Limits on an Intermediate Mass Black Hole in Omega Centauri: I. Hubble Space Telescope Photometry and Proper Motions

We analyze data from the Hubble Space Telescope's Advanced Camera for Surveys of the globular cluster Omega Cen. We construct a photometric and proper-motion catalog using the GO-9442, GO-10252, and GO-10775 data sets. The 2.5- to 4-year baseline between observations yields a catalog of some 10^5 proper motions, with 53,382 high-quality measurements in a central field. We determine the cluster center to ~1-arcsecond accuracy using two different star-count methods, and a completely independent method using 2MASS images. We also determine the kinematical center of the proper motions, which agrees with the star-count center to within its uncertainty. The proper-motion dispersion of the cluster increases gradually inwards, but there is no variation in kinematics with position within the central ~15 arcsec: there is no dispersion cusp and no stars with unusually high velocities. We measure for the first time in any globular cluster the variation in proper-motion dispersion with mass along the main sequence, and find the cluster not yet to be in equipartition. Our proper-motion results do not confirm the arguments put forward by Noyola, Gebhardt & Bergmann to suspect an intermediate-mass black hole (IMBH) in Omega Cen. In Paper II we present new dynamical models for the high-quality data presented here, with the aim of putting quantitative contraints on the mass of any possible IMBH.

preprint2009arXiv

New Limits on an Intermediate Mass Black Hole in Omega Centauri: II. Dynamical Models

We present a detailed dynamical analysis of the projected density and kinematical data available for the globular cluster Omega Cen. We solve the spherical anisotropic Jeans equation to predict the projected profiles of the RMS velocity in each of the three orthogonal coordinate directions (line of sight, proper motion radial, and proper motion tangential). We fit the models to new HST star count and proper motion data near the cluster center presented in Paper I, combined with existing ground-based measurements. We also derive and model the Gauss-Hermite moments of the observed proper motion distributions. The projected density profile is consistent with being flat near the center, with an upper limit gamma=0.07 on the central logarithmic slope. The RMS proper motion profile is also consistent with being flat near the center, and there are no unusually fast-moving stars. The models provide a good fit and yield a 1-sigma upper limit MBH < 1.2E4 solar masses on the mass of a possible intermediate-mass black hole (IMBH). The inferred upper limit corresponds to MBH/Mtot < 0.43%. We combine this with results for other clusters and discuss the implications for globular cluster IMBH demographics. Tighter limits will be needed to rule out or establish whether globular clusters follow the same black hole demographics correlations as galaxies. The arguments put forward by Noyola et al. (2008) to suspect an IMBH in Omega Cen are not confirmed by our study; the IMBH mass they suggested is firmly ruled out.

preprint2006arXiv

HST/STIS Spectra of Nuclear Star Clusters in Spiral Galaxies: Dependence of Age and Mass On Hubble Type

(Abridged) We study the nuclear star clusters in spiral galaxies of various Hubble types using spectra obtained with STIS on-board HST. We observed the nuclear clusters in 40 galaxies, selected from two previous HST/WFPC2 imaging surveys. The spectra provide a better separation of cluster light from underlying galaxy light than is possible with ground-based spectra. To infer the star formation history, metallicity and dust extinction, we fit weighted superpositions of single-age stellar population templates to the spectra. The luminosity-weighted age ranges from 10 Myrs to 10 Gyrs. The stellar populations of NCs are generally best fit as a mixture of populations of different ages. This indicates that NCs did not form in a single event, but instead they had additional star formation long after the oldest stars formed. On average, the sample clusters in late-type spirals have a younger luminosity-weighted mean age than those in early-type spirals (log(age/yr) = 8.37+/-0.25 vs. 9.23+/-0.21). The average cluster masses are smaller in late-type spirals than in early-type spirals (log(M/Msun) = 6.25+/-0.21 vs. 7.63+/-0.24), and exceed the masses typical of globular clusters. The cluster mass correlates strongly with both the Hubble type of the host galaxy and the luminosity of its bulge. The latter correlation has the same slope as the well-known correlation between supermassive black hole mass and bulge luminosity. The properties of both nuclear clusters and black holes are therefore intimately connected to the properties of the host galaxy.

preprint2006arXiv

Resolved Massive Star Clusters in the Milky Way and its Satellites: Brightness Profiles and a Catalogue of Fundamental Parameters

[Abridged]: We present a database of structural and dynamical properties for 153 spatially resolved star clusters (50 "young massive clusters" and 103 old globulars) in the Milky Way, the Large and Small Magellanic Clouds, and the Fornax dwarf spheroidal. This database complements and extends others in the literature, such as those of Harris, and Mackey & Gilmore. By fitting a number of models to the clusters' density profiles, we derive various characteristic surface brightnesses and radii; central potentials, velocity dispersions, and escape velocities; total luminosities, masses, and binding energies; phase-space densities and relaxation timescales; and ``kappa-space'' parameters. Population-synthesis models are used to predict intrinsic (B-V) colors, reddenings, and V-band mass-to-light ratios for the same 153 clusters plus another 63 globulars in the Milky Way, and we compare these predictions to the observed quantities where available. These results are intended to serve as the basis for future investigations of structural correlations and the fundamental plane of massive star clusters, including especially comparisons between the systemic properties of young and old clusters. We also address the question of what structural model fits each cluster best, and argue that the extended halos known to characterize many Magellanic Cloud clusters may be examples of the generic envelope structure of self-gravitating star clusters, not just transient features associated strictly with young age.

preprint2006arXiv

Stellar populations in the nuclei of late-type spiral galaxies

(Abridged) As part of an ongoing effort to study the stellar nuclei of very late-type, bulge-less spirals, we present results from a high-resolution spectroscopic survey of nine such nuclear star clusters, undertaken with VLT/UVES. We fit the spectra with population synthesis models and measure Lick-type indices to determine mean luminosity-weighted ages, which range from 4.1*10^7 to 1.1*10^10 years and are insensitive to assumed metallicity or internal extinction. The average metallicity of nuclear clusters in late-type spirals is slightly sub-solar (<Z> = 0.015) but shows significant scatter. The fits also show that the nuclear cluster spectra are best described by a mix of several generations of stars. This is supported by the fact that only models with composite stellar populations yield mass-to-light ratios that match those obtained from dynamical measurements. The last star formation episode was on average 34 Myr ago, while all clusters experienced some star formation in the last 100 Myr. We thus conclude that the nuclear clusters undergo repeated episodes of star formation. The robustness with respect to possible contamination from the underlying galaxy disk is demonstrated by comparison to spectra obtained with HST/STIS. Combining these results with those from Walcher et al. (2005), we have thus shown that the stellar nuclei of these bulge-less galaxies are massive and dense star clusters that form stars recurrently until the present day. This unique set of properties is likely due to the central location of these clusters in their host galaxies.

preprint1996arXiv

Scale-free dynamical models for galaxies: flattened densities in spherical potentials

This paper presents two families of phase-space distribution functions (DFs) that generate scale-free spheroidal mass densities in scale-free spherical potentials. The `case I' DFs are anisotropic generalizations of the flattened f(E,L_z) model, which they include as a special case. The `case II' DFs generate flattened constant-anisotropy models, in which the constant ratio of rms tangential to radial motion is characterized by Binney's parameter beta. The models can describe the outer parts of galaxies and the density cusp structure near a central black hole, but also provide general insight into the dynamical properties of flattened systems. The dependence of the intrinsic and projected properties on the model parameters and the inclination is described. The observed ratio of the rms projected line-of-sight velocities on the projected major and minor axes of elliptical galaxies is best fit by the case II models with beta > 0. These models also predict non-Gaussian velocity profile shapes consistent with existing observations. The distribution functions are used to model the galaxies NGC 2434 (E1) and NGC 3706 (E4), for which stellar kinematical measurements out to two effective radii indicate the presence of dark halos (Carollo et al.). The velocity profile shapes of both galaxies can be well fit by radially anisotropic case II models with a spherical logarithmic potential. This contrasts with the f(E,L_z) models studied previously, which require flattened dark halos to fit the data.

preprint1994arXiv

Dynamics of the nuclear gas & dust disc in the E4 radio galaxy NGC 7052

We present high spatial resolution ground-based broad-band imaging, H-alpha +[NII] narrow-band imaging and long-slit spectroscopy for the E4 radio galaxy NGC 7052, which has a nuclear dust disc. We detect ionized gas with a LINER spectrum, residing also in a nuclear disc. The gas rotates rapidly and the emission line widths increase towards the nucleus. The images are well fit by an axisymmetric model with the gas and dust in a disc of 1.5" radius (340 pc) in the equatorial plane of the stellar body, viewed at an inclination of 70 degrees. We assume the gas to be on circular orbits in the equatorial plane, with in addition a local turbulent velocity dispersion. The circular velocity is calculated from the combined gravitational potential of the stars and a possible nuclear black hole. The observed gas rotation curve is well fit, either with or without a black hole. Turbulent velocities >300 km/s must be present at radii <0.5" to fit the observed nuclear line widths. Models with a black hole and no turbulence can also fit the line widths, but these models cannot fit the observed nuclear line shapes. Models with both a black hole and gas turbulence can fit the data well, but the black hole is not required by the data, and if present, its mass must be < 5 x 10^8 solar masses. This upper limit is 5 times smaller than the black hole mass inferred for M87 from HST data. Our modelling and analysis techniques will be useful also for the study of the kinematics of nuclear gas discs in other galaxies. Such data will yield important information on the presence of massive black holes.