The Large Magellanic Cloud (LMC) is the closest and most studied example of an irregular galaxy. Among its principal defining morphological features, its off-centred bar and single spiral arm stand out, defining a whole family of galaxies known as the Magellanic spirals (Sm). These structures are thought to be triggered by tidal interactions and possibly maintained via gas accretion. However, it is still unknown whether they are long-lived stable structures. In this work, by combining photometry that reaches down to the oldest main sequence turn-off in the colour-magnitude diagrams (CMD, up to a distance of $\sim$4.4 kpc from the LMC centre) from the SMASH survey and CMD fitting techniques, we find compelling evidence supporting the long-term stability of the LMC spiral arm, dating the origin of this structure to more than 2~Gyr ago. The evidence suggests that the close encounter between the LMC and the Small Magellanic Cloud (SMC) that produced the gaseous Magellanic Stream and its Leading Arm (LA) also triggered the formation of the LMC's spiral arm. Given the mass difference between the Clouds and the notable consequences of this interaction, we can speculate that this should have been one of their closest encounters. These results set important constraints on the timing of LMC-SMC collisions, as well as on the physics behind star formation induced by tidal encounters.

We report the discovery of two ultra-faint stellar systems found in early data from the DECam Local Volume Exploration survey (DELVE). The first system, Centaurus I (DELVE J1238-4054), is identified as a resolved overdensity of old and metal-poor stars with a heliocentric distance of ${\rm D}_{\odot} = 116.3_{-0.6}^{+0.6}$ kpc, a half-light radius of $r_h = 2.3_{-0.3}^{+0.4}$ arcmin, an age of $\tau > 12.85$ Gyr, a metallicity of $Z = 0.0002_{-0.0002}^{+0.0001}$, and an absolute magnitude of $M_V = -5.55_{-0.11}^{+0.11}$ mag. This characterization is consistent with the population of ultra-faint satellites, and confirmation of this system would make Centaurus I one of the brightest recently discovered ultra-faint dwarf galaxies. Centaurus I is detected in Gaia DR2 with a clear and distinct proper motion signal, confirming that it is a real association of stars distinct from the Milky Way foreground; this is further supported by the clustering of blue horizontal branch stars near the centroid of the system. The second system, DELVE 1 (DELVE J1630-0058), is identified as a resolved overdensity of stars with a heliocentric distance of ${\rm D}_{\odot} = 19.0_{-0.6}^{+0.5} kpc$, a half-light radius of $r_h = 0.97_{-0.17}^{+0.24}$ arcmin, an age of $\tau = 12.5_{-0.7}^{+1.0}$ Gyr, a metallicity of $Z = 0.0005_{-0.0001}^{+0.0002}$, and an absolute magnitude of $M_V = -0.2_{-0.6}^{+0.8}$ mag, consistent with the known population of faint halo star clusters. Given the low number of probable member stars at magnitudes accessible with Gaia DR2, a proper motion signal for DELVE 1 is only marginally detected. We compare the spatial position and proper motion of both Centaurus I and DELVE 1 with simulations of the accreted satellite population of the Large Magellanic Cloud (LMC) and find that neither is likely to be associated with the LMC.

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.

A number of scenarios for the formation of multiple populations in globular clusters (GCs) predict that second generation (2G) stars form in a compact and dense subsystem embedded in a more extended first-generation (1G) system. If these scenarios are accurate, a consequence of the denser 2G formation environment is that 2G binaries should be more significantly affected by stellar interactions and disrupted at a larger rate than 1G binaries. The fractions and properties of binary stars can thus provide a dynamical fingerprint of the formation epoch of multiple-population GCs and their subsequent dynamical evolution. We investigate the connection between binaries and multiple populations in five GCs, NGC 288, NGC 6121 (M 4), NGC 6352, NGC 6362, and NGC 6838 (M 71). To do this, we introduce a new method based on the comparison of Hubble Space Telescope observations of binaries in the F275W, F336W, F438W, F606W and F814W filters with a large number of simulated binaries. In the inner regions probed by our data we do not find large differences between the local 1G and the 2G binary incidences in four of the studied clusters, the only exception being M 4 where the 1G binary incidence is about three times larger than the 2G incidence. The results found are in general agreement with the results of simulations predicting significant differences in the global 1G and 2G incidences and in the local values in the clusters' outer regions but similar incidences in the inner regions. The significant difference found in M 4 is consistent with simulations with a larger fraction of wider binaries. Our analysis also provides the first evidence of mixed (1G-2G) binaries, a population predicted by numerical simulations to form in a cluster's inner regions as a result of stellar encounters during which one component of a binary is replaced by a star of a different population.

In the present work we analyzed seven globular clusters selected from their location in the Galactic bulge and with metallicity values in the range $-1.30\lesssim\rm{[Fe/H]}\lesssim-0.50$. The aim of this work is first to derive cluster ages assuming single stellar populations, and secondly, to identify the stars from first (1G) and second generations (2G) from the main sequence, subgiant and red giant branches, and to derive their age differences. Based on a combination of UV and optical filters used in this project, we apply the Gaussian mixture models to distinguish the multiple stellar populations. Applying statistical isochrone fitting, we derive self-consistent ages, distances, metallicities, and reddening values for the sample clusters. An average of $12.3\pm0.4$ Gyr was obtained both using Dartmouth and BaSTI (accounting atomic diffusion effects) isochrones, without a clear distinction between the moderately metal-poor and the more metal-rich bulge clusters, except for NGC 6717 and the inner halo NGC 6362 with $\sim 13.5$ Gyr. We derived a weighted mean age difference between the multiple populations hosted by each globular cluster of $41\pm170$ Myr adopting canonical He abundances; whereas for higher He in 2G stars, this difference reduces to $17\pm170$ Myr, but with individual uncertainties of $500$ Myr.

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) will be the first high-performance stellar coronagraph using active wavefront control for deep starlight suppression in space, providing unprecedented levels of contrast, spatial resolution, and sensitivity for astronomical observations in the optical. One science case enabled by the CGI will be taking images and(R~50)spectra of faint interplanetary dust structures present in the habitable zone of nearby sunlike stars (~10 pc) and within the snow-line of more distant ones(~20pc), down to dust density levels commensurate with that of the solar system zodiacal cloud. Reaching contrast levels below~10-7 for the first time, CGI will cross an important threshold in debris disks physics, accessing disks with low enough optical depths that their structure is dominated by transport phenomena than collisions. Hence, CGI results will be crucial for determining how exozodiacal dust grains are produced and transported in low-density disks around mature stars. Additionally, CGI will be able to measure the brightness level and constrain the degree of asymmetry of exozodiacal clouds around individual nearby sunlike stars in the optical, at the ~10x solar zodiacal emission level. This information will be extremely valuable for optimizing the observational strategy of possible future exo-Earth direct imaging missions, especially those planning to operate at optical wavelengths, such as Habitable Exoplanet Observatory (HabEx) and the Large Ultraviolet/Optical/Infrared Surveyor (LUVOIR).

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.

Understanding the evolutionary history of the Magellanic Clouds requires an in-depth exploration and characterization of the stellar content in their outer regions, which ultimately are key to tracing the epochs and nature of past interactions. We present new deep images of a shell-like over-density of stars in the outskirts of the Small Magellanic Cloud (SMC). The shell, also detected in photographic plates dating back to the fifties, is located at ~1.9 degr from the center of the SMC in the north-east direction.The structure and stellar content of this feature were studied with multi-band, optical data from the Survey of the MAgellanic Stellar History (SMASH) carried out with the Dark Energy Camera on the Blanco Telescope at Cerro Tololo Inter-American Observatory. We also investigate the kinematic of the stars in the shell using the Gaia Data Release 2. The shell is composed of a young population with an age ~ 150 Myr, with no contribution from an old population. Thus, it is hard to explain its origin as the remnant of a tidally disrupted stellar system. The spatial distribution of the young main-sequence stars shows a rich sub-structure, with a spiral arm-like feature emanating from the main shell and a separated small arc of young stars close to the globular cluster NGC 362. We find that the absolute g-band magnitude of the shell is M_{g,shell} = -10.78+/- 0.02, with a surface brightness of mu_{g,shell} = 25.81+/- 0.01 mag/arcsec^{2}. We have not found any evidence that this feature is of tidal origin or a bright part of a spiral arm-like structure. Instead, we suggest that the shell formed in a recent star formation event, likely triggered by an interaction with the Large Magellanic Cloud and/or the Milky Way, ~150 Myr ago.

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

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

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.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

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

We present the analysis of the radial distributions and kinematic properties of the multiple stellar populations (mPOPs) hosted in the globular cluster (GC) NGC 6352 as part of the Hubble Space Telescope (HST) "UV Legacy Survey of Galactic Globular Clusters" program. NGC 6352 is one of the few GCs for which the mPOP tagging in appropriate color-magnitude diagrams is clear in all evolutionary sequences. We computed high-precision stellar proper motions for the stars from the cluster's core out to 75 arcsec ($\sim$1.5 core radii, or $\sim$0.6 half-light radii). We find that, in the region explored, first- and second-generation stars share the same radial distribution and kinematic properties. Velocity dispersions, anisotropy radial profiles, differential rotation, and level of energy equipartition, all suggest that NGC 6352 is probably in an advanced evolutionary stage, and any possible difference in the structural and kinematic properties of its mPOPs have been erased by dynamical processes in the core of the cluster. We also provide an estimate of the mass of blue stragglers and of main-sequence binaries through kinematics alone. In general, in order to build a complete dynamical picture of this and other GCs, it will be essential to extend the analyses presented in this paper to the GCs' outer regions where some memories of the initial differences in the mPOP properties, and those imprinted by dynamical processes, might still be present.

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

The Wide Field Infrared Survey Telescope (WFIRST) Coronagraph Instrument (CGI) is a high-contrast imager and integral field spectrograph that will enable the study of exoplanets and circumstellar disks at visible wavelengths. Ground-based high-contrast instrumentation has fundamentally limited performance at small working angles, even under optimistic assumptions for 30m-class telescopes. There is a strong scientific driver for better performance, particularly at visible wavelengths. Future flagship mission concepts aim to image Earth analogues with visible light flux ratios of more than 10^10. CGI is a critical intermediate step toward that goal, with a predicted 10^8-9 flux ratio capability in the visible. CGI achieves this through improvements over current ground and space systems in several areas: (i) Hardware: space-qualified (TRL9) deformable mirrors, detectors, and coronagraphs, (ii) Algorithms: wavefront sensing and control; post-processing of integral field spectrograph, polarimetric, and extended object data, and (iii) Validation of telescope and instrument models at high accuracy and precision. This white paper, submitted to the 2018 NAS Exoplanet Science Strategy call, describes the status of key CGI technologies and presents ways in which performance is likely to evolve as the CGI design matures.

We explore the stellar structure of the Large Magellanic Cloud (LMC) disk using data from the Survey of the MAgellanic Stellar History (SMASH) and the Dark Energy Survey. We detect a ring-like stellar overdensity in the red clump star count map at a radius of ~6 degrees (~5.2 kpc at the LMC distance) that is continuous over ~270 degrees in position angle and is only limited by the current data coverage. The overdensity shows an amplitude up to 2.5 times higher than that of the underlying smooth disk. This structure might be related to the multiple arms found by de Vaucouleurs. We find that the overdensity shows spatial correlation with intermediate-age star clusters, but not with young (< 1 Gyr) main-sequence stars, indicating the stellar populations associated with the overdensity are intermediate in age or older. Our findings on the LMC overdensity can be explained by either of two distinct formation mechanisms of a ring-like overdensity: (1) the overdensity formed out of an asymmetric one-armed spiral wrapping around the LMC main body, which is induced by repeated encounters with the Small Magellanic Cloud (SMC) over the last Gyr, or (2) the overdensity formed very recently as a tidal response to a direct collision with the SMC. Although the measured properties of the overdensity alone cannot distinguish between the two candidate scenarios, the consistency with both scenarios suggests that the ring-like overdensity is likely a product of tidal interaction with the SMC, but not with the Milky Way halo.

We present a study of the three-dimensional (3D) structure of the Large Magellanic Cloud (LMC) using ~2.2 million red clump (RC) stars selected from the Survey of the MAgellanic Stellar History. To correct for line-of-sight dust extinction, the intrinsic RC color and magnitude and their radial dependence are carefully measured by using internal nearly dust-free regions. These are then used to construct an accurate 2D reddening map (165 square degrees with ~10 arcmin resolution) of the LMC disk and the 3D spatial distribution of RC stars. An inclined disk model is fit to the 2D distance map yielding a best-fit inclination angle i = 25.86(+0.73,-1.39) degrees with random errors of +\-0.19 degrees and line-of-nodes position angle theta = 149.23(+6.43,-8.35) degrees with random errors of +/-0.49 degrees. These angles vary with galactic radius, indicating that the LMC disk is warped and twisted likely due to the repeated tidal interactions with the Small Magellanic Cloud (SMC). For the first time, our data reveal a significant warp in the southwestern part of the outer disk starting at rho ~ 7 degrees that departs from the defined LMC plane up to ~4 kpc toward the SMC, suggesting that it originated from a strong interaction with the SMC. In addition, the inner disk encompassing the off-centered bar appears to be tilted up to 5-15 degrees relative to the rest of the LMC disk. These findings on the outer warp and the tilted bar are consistent with the predictions from the Besla et al. simulation of a recent direct collision with the SMC.

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

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.

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

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

How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of $\sim$150 $M_\odot$ . It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We investigate the kinematics of VFTS682 as obtained by Gaia and Hubble Space Telescope astrometry. We derive a projected velocity relative to the cluster of $38 \pm 17 \mathrm{km \ s^{-1}}$ (1$\sigma$ confidence interval). Although the error bars are substantial, two independent measures suggest that VFTS682 is a runaway ejected from the central cluster. This hypothesis is further supported by a variety of circumstantial clues. The central cluster is known to harbor other stars more massive than 150 $M_\odot$ of similar spectral type and recent astrometric studies on VFTS16 and VFTS72 provide direct evidence that the cluster can eject some of its most massive members, in agreement with theoretical predictions. If future data confirm the runaway nature, this would make VFTS682 the most massive runaway star known to date.

Recent work, based on data from the Hubble Space Telescope (HST) UV Legacy Survey of Galactic Globular Clusters (GCs), has revealed that all the analyzed clusters host two groups of first- (1G) and second-generation (2G) stars. In most GCs, both 1G and 2G stars host sub-stellar populations with different chemical composition. We compare multi-wavelength HST photometry with synthetic spectra to determine for the first time the average helium difference between the 2G and 1G stars in a large sample of 57 GCs and the maximum helium variation within each of them. We find that in all clusters 2G stars are consistent with being enhanced in helium with respect to 1G. The maximum helium variation ranges from less than 0.01 to more than 0.10 in helium mass fraction and correlates with both the cluster mass and the color extension of the horizontal branch (HB). These findings demonstrate that the internal helium variation is one of the main (second) parameters governing the HB morphology.

In this paper we present the astro-photometric catalogues of 56 globular clusters and one open cluster. Astrometry and photometry are mainly based on images collected within the "HST Legacy Survey of Galactic Globular Clusters: Shedding UV Light on Their Populations and Formation" (GO-13297, PI:~Piotto), and the "ACS Survey of Galactic Globular Clusters" (GO-10775, PI:~Sarajedini). For each source in the catalogues for which we have reliable proper motion we also publish a membership probability for separation of field and cluster stars. These new catalogues, which we make public in Mikulski Archive for Space Telescopes, replace previous catalogues by Paper VIII of this series.

Line-of-sight kinematic studies indicate that many Galactic globular clusters have a significant degree of internal rotation. However, three-dimensional kinematics from a combination of proper motions and line-of-sight velocities are needed to unveil the role of angular momentum in the formation and evolution of these old stellar systems. Here we present the first quantitative study of internal rotation on the plane-of-the-sky for a large sample of globular clusters using proper motions from Gaia DR2. We detect signatures of rotation in the tangential component of proper motions for 11 out of 51 clusters at a $>$3-sigma confidence level, confirming the detection reported in Gaia collaboration et al. (2018) for 8 clusters, and additionally identify 11 GCs with a 2-sigma rotation detection. For the clusters with a detected global rotation, we construct the two-dimensional rotation maps and proper motion rotation curves, and we assess the relevance of rotation with respect to random motions ($V/\sigma\sim0.08-0.51$). We find evidence of a correlation between the degree of internal rotation and relaxation time, highlighting the importance of long-term dynamical evolution in shaping the clusters current properties. This is a strong indication that angular momentum must have played a fundamental role in the earliest phases of cluster formation. Finally, exploiting the spatial information of the rotation maps and a comparison with line-of-sight data, we provide an estimate of the inclination of the rotation axis for a subset of 8 clusters. Our work demonstrates the potential of Gaia data for internal kinematic studies of globular clusters and provides the first step to reconstruct their intrinsic three-dimensional structure.

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

We present a new distance estimation for the Milky Way dwarf spheroidal satellite Sculptor obtained from multi-epoch mid-infrared observations of RR Lyrae stars. The 3.6 {\mu}m observations have been acquired with the Infrared Array Camera on board the Spitzer Space Telescope as part of the SMHASH Program. Mid-infrared light curves for 42 RRL were obtained, from which we measured Sculptor's distance modulus to be {\mu} = 19.60 $\pm$ 0.02 (statistical) $\pm$ 0.04 (photometric) mag (with $\sigma_{sys}=$ = 0.09 mag), using the 3.6 {\mu}m empirical period-luminosity relations derived from the Galactic globular cluster M4, or {\mu} = 19.57 $\pm$ 0.02 (statistical) $\pm$ 0.04 (photometric) mag (with $\sigma_{sys}=$ = 0.11 mag) using empirical relations in the same passband recently derived from the Large Magellanic Cloud globular cluster Reticulum. Both these measurements are in good agreement with values presented in previous works with Sculptor RR Lyrae stars in optical bands, and are also consistent with recent near-infrared RR Lyrae results. Best agreement with the literature is found for the latter modulus which is equivalent to a distance of d = 82 $\pm$ 1 (statistical) $\pm$ 2 (photometric) kpc (with $\sigma_{sys}=$ = 4 kpc). Finally, using a subsample of RR Lyrae stars with spectroscopic metallicities, we demonstrate that these distance estimates are not affected by metallicity effects.

A previous spectroscopic study identified the very massive O2 III star VFTS 16 in the Tarantula Nebula as a runaway star based on its peculiar line-of-sight velocity. We use the Gaia DR2 catalog to measure the relative proper motion of VFTS 16 and nearby bright stars to test if this star might have been ejected from the central cluster, R136, via dynamical ejection. We find that the position angle and magnitude of the relative proper motion (0.338 +/- 0.046 mas/yr, or approximately 80 +\- 11 km/s) of VFTS 16 are consistent with ejection from R136 approximately 1.5 +/- 0.2 Myr ago, very soon after the cluster was formed. There is some tension with the presumed age of VFTS 16 that, from published stellar parameters, cannot be greater than 0.9 +0.3/-0.2 Myr. Older ages for this star would appear to be prohibited due to the absence of He I lines in its optical spectrum, since this sets a firm lower limit on its effective temperature. The dynamical constraints may imply an unusual evolutionary history for this object, perhaps indicating it is a merger product. Gaia DR2 also confirms that another very massive star in the Tarantula Nebula, VFTS 72 (alias BI253; O2 III-V(n)((f*)), is also a runaway on the basis of its proper motion as measured by Gaia. While its tangential proper motion (0.392 +/-0.062 mas/yr or 93 +/-15 km/s) would be consistent with dynamical ejection from R136 approximately 1 Myr ago, its position angle is discrepant with this direction at the 2$\sigma$ level. From their Gaia DR2 proper motions we conclude that the two ~100 solar mass O2 stars, VFTS 16 and VFTS72, are fast runaway stars, with space velocities of around 100 km/s relative to R136 and the local massive star population. The dynamics of VFTS16 are consistent with it having been ejected from R136, and this star therefore sets a robust lower limit on the age of the central cluster of ~1.3 Myr.

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

The Hubble Tarantula Treasury Project (HTTP) has provided an unprecedented photometric coverage of the entire star-burst region of 30 Doradus down to the half Solar mass limit. We use the deep stellar catalogue of HTTP to identify all the pre--main-sequence (PMS) stars of the region, i.e., stars that have not started their lives on the main-sequence yet. The photometric distinction of these stars from the more evolved populations is not a trivial task due to several factors that alter their colour-magnitude diagram positions. The identification of PMS stars requires, thus, sophisticated statistical methods. We employ Machine Learning Classification techniques on the HTTP survey of more than 800,000 sources to identify the PMS stellar content of the observed field. Our methodology consists of 1) carefully selecting the most probable low-mass PMS stellar population of the star-forming cluster NGC 2070, 2) using this sample to train classification algorithms to build a predictive model for PMS stars, and 3) applying this model in order to identify the most probable PMS content across the entire Tarantula Nebula. We employ Decision Tree, Random Forest and Support Vector Machine classifiers to categorise the stars as PMS and Non-PMS. The Random Forest and Support Vector Machine provided the most accurate models, predicting about 20,000 sources with a candidateship probability higher than 50 percent, and almost 10,000 PMS candidates with a probability higher than 95 percent. This is the richest and most accurate photometric catalogue of extragalactic PMS candidates across the extent of a whole star-forming complex.

We present an improved data-reduction technique to obtain high-precision proper motions (PMs) of globular clusters using Hubble Space Telescope data. The new reduction is superior to the one presented in the first paper of this series for the faintest sources in very crowded fields. We choose the globular cluster NGC 362 as a benchmark to test our new procedures. We measure PMs of 117 450 sources in the field, showing that we are able to obtain a PM precision better than 10 $\mu$as yr$^{-1}$ for bright stars. We make use of this new PM catalog of NGC 362 to study the cluster's internal kinematics. We investigate the velocity-dispersion profiles of the multiple stellar populations hosted by NGC 362 and find new pieces of information on the kinematics of first- and second-generation stars. We analyze the level of energy equipartition of the cluster and find direct evidence for its post-core-collapsed state from kinematic arguments alone. We refine the dynamical mass of the blue stragglers and study possible kinematic differences between blue stragglers formed by collisions and mass transfer. We also measure no significant cluster rotation in the plane of the sky. Finally, we measure the absolute PM of NGC 362 and of the background stars belonging to the Small Magellanic Cloud, finding a good agreement with previous estimates in the literature. We make the PM catalog publicly available.

We present the detection of very extended stellar populations around the Large Magellanic Cloud (LMC) out to R~21 degrees, or ~18.5 kpc at the LMC distance of 50 kpc, as detected in the Survey of the MAgellanic Stellar History (SMASH) performed with the Dark Energy Camera on the NOAO Blanco 4m Telescope. The deep (g~24) SMASH color magnitude diagrams (CMDs) clearly reveal old (~9 Gyr), metal-poor ([Fe/H]=-0.8 dex) main-sequence stars at a distance of 50 kpc. The surface brightness of these detections is extremely low with our most distant detection having 34 mag per arcsec squared in g-band. The SMASH radial density profile breaks from the inner LMC exponential decline at ~13-15 degrees and a second component at larger radii has a shallower slope with power-law index of -2.2 that contributes ~0.4% of the LMC's total stellar mass. In addition, the SMASH densities exhibit large scatter around our best-fit model of ~70% indicating that the envelope of stellar material in the LMC periphery is highly disturbed. We also use data from the NOAO Source catalog to map the LMC main-sequence populations at intermediate radii and detect a steep dropoff in density on the eastern side of the LMC (at R~8 deg) as well as an extended structure to the far northeast. These combined results confirm the existence of a very extended, low-density envelope of stellar material with disturbed shape around the LMC. The exact origin of this structure remains unclear but the leading options include a classical accreted halo or tidally stripped outer disk material.

The parameter A+, defined as the area enclosed between the cumulative radial distribution of blue straggler stars (BSSs) and that of a reference population, is a powerful indicator of the level of BSS central segregation. As part of the Hubble Space Telescope UV Legacy Survey of Galactic globular clusters (GCs), here we present the BSS population and the determination of A+ in 27 GCs observed out to about one half-mass radius. In combination with 21 additional clusters discussed in a previous paper this provides us with a global sample of 48 systems (corresponding to \sim 32\% of the Milky Way GC population), for which we find a strong correlation between A+ and the ratio of cluster age to the current central relaxation time. Tight relations have been found also with the core radius and the central luminosity density, which are expected to change with the long-term cluster dynamical evolution. An interesting relation is emerging between A+ and the ratio of the BSS velocity dispersion relative to that of main sequence turn-off stars, which measures the degree of energy equipartition experienced by BSSs in the cluster. These results provide further confirmation that BSSs are invaluable probes of GC internal dynamics and A+ is a powerful dynamical clock.

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

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

The Wide Field Infrared Survey Telescope (WFIRST) was the top ranked large space mission in the 2010 New Worlds, New Horizons decadal survey, and it was formed by merging the science programs of 3 different mission concepts, including the Microlensing Planet Finder (MPF) concept (Bennett \etal\ 2010). The WFIRST science program (Spergel \etal\ 2015) consists of a general observer program, a wavefront controlled technology program, and two targeted science programs: a program to study dark energy, and a statistical census of exoplanets with a microlensing survey, which uses nearly one quarter of WFIRST's observing time in the current design reference mission. The New Worlds, New Horizons (decadal survey) midterm assessment summarizes the science case for the WFIRST exoplanet microlensing survey with this statement: "WFIRST's microlensing census of planets beyond 1 AU will perfectly complement Kepler's census of compact systems, and WFIRST will also be able to detect free-floating planets unbound from their parent stars\rlap."

We present Magellan/IMACS, Anglo-Australian Telescope/AAOmega+2dF, and Very Large Telescope/GIRAFFE+FLAMES spectroscopy of the CarinaII (Car II) & Carina III (Car III) dwarf galaxy candidates, recently discovered in the Magellanic Satellites Survey (MagLiteS). We identify 18 member stars in Car II, including 2 binaries with variable radial velocities and 2 RR Lyrae stars. The other 14 members have a mean heliocentric velocity $v_{\rm hel} = 477.2 \pm 1.2$ km/s and a velocity dispersion of $\sigma_v = 3.4^{+1.2}_{-0.8}$ km/s. Assuming Car II is in dynamical equilibrium, we derive a total mass within the half-light radius of $1.0^{+0.8}_{-0.4} \times 10^{6} M_\odot$, indicating a mass-to-light ratio of $369^{+309}_{-161} M_\odot/L_\odot$. From equivalent width measurements of the calcium triplet lines of 9 RGB stars, we derive a mean metallicity of [Fe/H] = $-2.44 \pm 0.09$ with dispersion $\sigma_{\rm [Fe/H]} = 0.22 ^{+0.10}_{-0.07}$. Considering both the kinematic and chemical properties, we conclude that Car II is a dark-matter-dominated dwarf galaxy. For Car III, we identify 4 member stars, from which we calculate a systemic velocity of $v_{\rm hel} = 284.6^{+3.4}_{-3.1}$ km/s. The brightest RGB member of Car III has a metallicity of [Fe/H] $= -1.97 \pm 0.12$. Due to the small size of the Car III spectroscopic sample, we cannot conclusively determine its nature. Although these two systems have the smallest known physical separation ($\Delta d\sim10~kpc$) among Local Group satellites, the large difference in their systemic velocities, $\sim200$ km/s, indicates that they are unlikely to be a bound pair. One or both systems are likely associated with the Large Magellanic Cloud (LMC), and may remain LMC satellites today. No statistically significant excess of $\gamma$-rays emission is found at the locations of Car II and Car III in eight years of Fermi-LAT data.

We present intermediate-resolution (R$\sim$1000) spectra in the $\sim$3500-10,000 A range of 14 globular clusters in the magellanic irregular galaxy NGC 4449 acquired with the Multi Object Double Spectrograph on the Large Binocular Telescope. We derived Lick indices in the optical and the CaII-triplet index in the near-infrared in order to infer the clusters' stellar population properties. The inferred cluster ages are typically older than $\sim$9 Gyr, although ages are derived with large uncertainties. The clusters exhibit intermediate metallicities, in the range $-1.2\lesssim$[Fe/H]$\lesssim-0.7$, and typically sub-solar [$\alpha/Fe$] ratios, with a peak at $\sim-0.4$. These properties suggest that i) during the first few Gyrs NGC 4449 formed stars slowly and inefficiently, with galactic winds having possibly contributed to the expulsion of the $\alpha$-elements, and ii) globular clusters in NGC 4449 formed relatively "late", from a medium already enriched in the products of type Ia supernovae. The majority of clusters appear also under-abundant in CN compared to Milky Way halo globular clusters, perhaps because of the lack of a conspicuous N-enriched, second-generation of stars like that observed in Galactic globular clusters. Using the cluster velocities, we infer the dynamical mass of NGC 4449 inside 2.88 kpc to be M($<$2.88 kpc)=$3.15^{+3.16}_{-0.75} \times 10^9~M_\odot$. We also report the serendipitous discovery of a planetary nebula within one of the targeted clusters, a rather rare event.

As part of the Hubble Space Telescope UV Legacy Survey of Galactic Globular Clusters, 110 parallel fields were observed with the Wide Field Channel of the Advanced Camera for Surveys, in the outskirts of 48 globular clusters, plus the open cluster NGC 6791. Totalling about $0.3$ square degrees of observed sky, this is the largest homogeneous Hubble Space Telescope photometric survey of Galalctic globular clusters outskirts to date. In particular, two distinct pointings have been obtained for each target on average, all centred at about $6.5$ arcmin from the cluster centre, thus covering a mean area of about $23\,{\rm arcmin^{2}}$ for each globular cluster. For each field, at least one exposure in both F475W and F814W filters was collected. In this work, we publicly release the astrometric and photometric catalogues and the astrometrised atlases for each of these fields.

The Hubble Space Telescope UV Legacy Survey of Galactic Globular Clusters is providing a major breakthrough in our knowledge of Globular Clusters (GCs) and their stellar populations. Among the main results, we discovered that all the studied GCs host two main discrete groups consisting of first generation (1G) and second generation (2G) stars. We exploit the multiwavelength photometry from this project to investigate, for the first time, the Red Giant Branch Bump (RGBB) of the two generations in a large sample of GCs. We identified, with high statistical significance, the RGBB of 1G and 2G stars in 26 GCs and found that their magnitude separation as a function of the filter wavelength follows comparable trends. The comparison of observations to synthetic spectra reveals that the RGBB luminosity depends on the stellar chemical composition and that the 2G RGBB is consistent with stars enhanced in He and N and depleted in C and O with respect to 1G stars. For metal-poor GCs the 1G and 2G RGBB relative luminosity in optical bands mostly depends on helium content, Y. We used the RGBB observations in F606W and F814W bands to infer the relative helium abundance of 1G and 2G stars in 18 GCs, finding an average helium enhancement $\Delta$Y=0.011$\pm$0.002 of 2G stars with respect to 1G stars. This is the first determination of the average difference in helium abundance of multiple populations in a large number of clusters and provides a lower limit to the maximum internal variation of helium in GCs.

In this second installment of the series, we look at the internal kinematics of the multiple stellar populations of the globular cluster $\omega$ Centauri in one of the parallel Hubble Space Telescope (HST) fields, located at about 3.5 half-light radii from the center of the cluster. Thanks to the over 15-year-long baseline and the exquisite astrometric precision of the HST cameras, well-measured stars in our proper-motion catalog have errors as low as $\sim 10\ \mu$as yr$^{-1}$, and the catalog itself extends to near the hydrogen-burning limit of the cluster. We show that second-generation (2G) stars are significantly more radially anisotropic than first-generation (1G) stars. The latter are instead consistent with an isotropic velocity distribution. In addition, 1G have excess systemic rotation in the plane of the sky with respect to 2G stars. We show that the six populations below the main-sequence (MS) knee identified in our first paper are associated to the five main population groups recently isolated on the upper MS in the core of cluster. Furthermore, we find both 1G and 2G stars in the field to be far from being in energy equipartition, with $\eta_{\rm 1G}=-0.007\pm0.026$ for the former, and $\eta_{\rm 2G}=0.074\pm0.029$ for the latter, where $\eta$ is defined so that the velocity dispersion $\sigma_\mu$ scales with stellar mass as $\sigma_\mu \propto m^{-\eta}$. The kinematical differences reported here can help constrain the formation mechanisms for the multiple stellar populations in $\omega$ Centauri and other globular clusters. We make our astro-photometric catalog publicly available.

In this paper we report a new estimate of the absolute proper motion (PM) of the globular cluster NGC 5139 ($\omega$ Cen) as part of the HST large program GO-14118+14662. We analyzed a field 17 arcmin South-West of the center of $\omega$ Cen and computed PMs with an epoch span of $\sim$15.1 years. We employed 45 background galaxies to link our relative PMs to an absolute reference-frame system. The absolute PM of the cluster in our field is: $(\mu_\alpha \cos\delta , \mu_\delta) = (-3.341 \pm 0.028 , -6.557 \pm 0.043)$ mas yr$^{-1}$. Upon correction for the effects of viewing perspective and the known cluster rotation, this implies that for the cluster center of mass $(\mu_\alpha \cos\delta , \mu_\delta) = (-3.238 \pm 0.028, -6.716 \pm 0.043)$ mas yr$^{-1}$. This measurement is direct and independent, has the highest random and systematic accuracy to date, and will provide an external verification for the upcoming Gaia Data Release 2. It also differs from most reported PMs for $\omega$ Cen in the literature by more than 5$\sigma$, but consistency checks compared to other recent catalogs yield excellent agreement. We computed the corresponding Galactocentric velocity, calculated the implied orbit of $\omega$ Cen in two different Galactic potentials, and compared these orbits to the orbits implied by one of the PM measurements available in the literature. We find a larger (by about 500 pc) perigalactic distance for $\omega$ Cen with our new PM measurement, suggesting a larger survival expectancy for the cluster in the Galaxy.

Stellar tidal streams provide an opportunity to study the motion and structure of the disrupting galaxy as well as the gravitational potential of its host. Streams around the Milky Way are especially promising as phase space positions of individual stars will be measured by ongoing or upcoming surveys. Nevertheless, it remains a challenge to accurately assess distances to stars farther than 10 kpc from the Sun, where we have the poorest knowledge of the Galaxy's mass distribution. To address this we present observations of 32 candidate RR Lyrae stars in the Orphan tidal stream taken as part of the Spitzer Merger History and Shape of the Galactic Halo (SMHASH) program. The extremely tight correlation between the periods, luminosities, and metallicities of RR Lyrae variable stars in the Spitzer IRAC $\mathrm{3.6 \mu m}$ band allows the determination of precise distances to individual stars; the median statistical distance uncertainty to each RR Lyrae star is $2.5\%$. By fitting orbits in an example potential we obtain an upper limit on the mass of the Milky Way interior to 60 kpc of $\mathrm{5.6_{-1.1}^{+1.2}\times 10^{11}\ M_\odot}$, bringing estimates based on the Orphan Stream in line with those using other tracers. The SMHASH data also resolve the stream in line--of--sight depth, allowing a new perspective on the internal structure of the disrupted dwarf galaxy. Comparing with N--body models we find that the progenitor had an initial dark halo mass of approximately $\mathrm{3.2 \times 10^{9}\ M_\odot}$, placing the Orphan Stream's progenitor amongst the classical dwarf spheroidals.

We present the discovery of a very faint stellar system, SMASH 1, that is potentially a satellite of the Large Magellanic Cloud. Found within the Survey of the MAgellanic Stellar History (SMASH), SMASH 1 is a compact ($r_h = 9.1^{+5.9}_{-3.4}$ pc) and very low luminosity (M_V = -1.0 +/- 0.9, $L_V=10^{2.3 +/- 0.4}$ Lsun) stellar system that is revealed by its sparsely populated main sequence and a handful of red-giant-branch candidate member stars. The photometric properties of these stars are compatible with a metal-poor ([Fe/H]=-2.2) and old (13 Gyr) isochrone located at a distance modulus of ~18.8, i.e. a distance of ~57 kpc. Situated at 11.3$^\circ$ from the LMC in projection, its 3-dimensional distance from the Cloud is ~13 kpc, consistent with a connection to the LMC, whose tidal radius is at least 16 kpc. Although the nature of SMASH 1 remains uncertain, its compactness favors it being a stellar cluster and hence dark-matter free. If this is the case, its dynamical tidal radius is only <19 pc at this distance from the LMC, and smaller than the system's extent on the sky. Its low luminosity and apparent high ellipticity ($\epsilon=0.62^{+0.17}_{-0.21}$) with its major axis pointing toward the LMC may well be the tell-tale sign of its imminent tidal demise.

We present new proper motion (PM) measurements of the dwarf spheroidal galaxies (dSphs) Draco and Sculptor using multi-epoch images obtained with the Hubble Space Telescope ACS/WFC. Our PM results have uncertainties far lower than previous measurements, even made with the same instrument. The PM results for Draco and Sculptor are (mu_W,mu_N)_Dra = (-0.0562+/-0.0099,-0.1765+/-0.0100) mas/yr and (mu_W,mu_N)_Scl = (-0.0296+/-0.0209,-0.1358 +/-0.0214) mas/yr. The implied Galactocentric velocity vectors for Draco and Sculptor have radial and tangential components: (V_rad,V_tan)_Dra = (-88.6,161.4) +/- (4.4,5.6) km/s; and (V_rad,V_tan)_Scl = (72.6,200.2) +/- (1.3,10.8) km/s. We study the detailed orbital history of both Draco and Sculptor via numerical orbit integrations. Orbital periods of Draco and Sculptor are found to be 1-2 and 2-5 Gyrs, respectively, accounting for uncertainties in the MW mass. We also study the influence of the Large Magellanic Cloud (LMC) on the orbits of Draco and Sculptor. Overall, the inclusion of the LMC increases the scatter in the orbital results. Based on our calculations, Draco shows a rather wide range of orbital parameters depending on the MW mass and inclusion/exclusion of the LMC, but Sculptor's orbit is very well constrained with its most recent pericentric approach to the MW being 0.3-0.4 Gyr ago. Our new PMs imply that the orbital trajectories of both Draco and Sculptor are confined within the Disk of Satellites (DoS), better so than implied by earlier PM measurements, and likely rule out the possibility that these two galaxies were accreted together as part of a tightly bound group.

The Large and Small Magellanic Clouds (LMC and SMC) are unique local laboratories for studying the formation and evolution of small galaxies in exquisite detail. The Survey of the MAgellanic Stellar History (SMASH) is an NOAO community DECam survey of the Clouds mapping 480 square degrees (distributed over ~2400 square degrees at ~20% filling factor) to ~24th mag in ugriz with the goal of identifying broadly distributed, low surface brightness stellar populations associated with the stellar halos and tidal debris of the Clouds. SMASH will also derive spatially-resolved star formation histories covering all ages out to large radii from the MCs that will further complement our understanding of their formation. Here, we present a summary of the survey, its data reduction, and a description of the first public Data Release (DR1). The SMASH DECam data have been reduced with a combination of the NOAO Community Pipeline, PHOTRED, an automated PSF photometry pipeline based mainly on the DAOPHOT suite, and custom calibration software. The attained astrometric precision is ~15 mas and the accuracy is ~2 mas with respect to the Gaia DR1 astrometric reference frame. The photometric precision is ~0.5-0.7% in griz and ~1% in u with a calibration accuracy of ~1.3% in all bands. The median 5 sigma point source depths in ugriz bands are 23.9, 24.8, 24.5, 24.2, 23.5 mag. The SMASH data already have been used to discover the Hydra II Milky Way satellite, the SMASH 1 old globular cluster likely associated with the LMC, and very extended stellar populations around the LMC out to R~18.4 kpc. SMASH DR1 contains measurements of ~100 million objects distributed in 61 fields. A prototype version of the NOAO Data Lab provides data access, including a data discovery tool, SMASH database access, an image cutout service, and a Jupyter notebook server with example notebooks for exploratory analysis.

We constrain the possible presence of a central black hole (BH) in the center of the Large Magellanic Cloud (LMC). This requires spectroscopic measurements over an area of order a square degree, due to the poorly known position of the kinematic center. Such measurements are now possible with the impressive field of view of the Multi Unit Spectroscopic Explorer (MUSE) on the ESO Very Large Telescope. We used the Calcium Triplet (~850nm) spectral lines in many short-exposure MUSE pointings to create a two-dimensional integrated-light line-of-sight velocity map from the ~$10^8$ individual spectra, taking care to identify and remove Galactic foreground populations. The data reveal a clear velocity gradient at an unprecedented spatial resolution of 1 arcmin$^{2}$. We fit kinematic models to arrive at a $3\sigma$ upper-mass-limit of $10^{7.1}$ M$_{Sun}$ for any central BH - consistent with the known scaling relations for supermassive black holes and their host systems. This adds to the growing body of knowledge on the presence of BHs in low-mass and dwarf galaxies, and their scaling relations with host-galaxy properties, which can shed light on theories of BH growth and host system interaction.

We report the discovery of a stellar over-density 8$^{\circ}$ north of the center of the Small Magellanic Cloud (Small Magellanic Cloud Northern Over-Density; SMCNOD) using data from the first two years of the Dark Energy Survey (DES) and the first year of the MAGellanic SatelLITEs Survey (MagLiteS). The SMCNOD is indistinguishable in age, metallicity and distance from the nearby SMC stars, being primarly composed of intermediate-age stars (6 Gyr, Z=0.001), with a small fraction of young stars (1 Gyr, Z=0.01). The SMCNOD has an elongated shape with an ellipticity of 0.6 and a size of $\sim$ 6x2 deg. It has an absolute magnitude of $M_V \cong$ -7.7, $r_h = 2.1$ kpc, and $\mu_V(r<r_h)$ = 31.2 mag arcsec$^{-2}$. We estimate a stellar mass of $\sim 10^5$ $M_{\odot}$, following a Kroupa mass function. The SMCNOD was probably removed from the SMC disk by tidal stripping, since it is located near the head of the Magellanic Stream, and the literature indicates likely recent LMC-SMC encounters. This scenario is supported by the lack of significant HI gas. Other potential scenarios for the SMCNOD origin are a transient over-density within the SMC tidal radius or a primordial SMC satellite in advanced stage of disruption.

High-precision proper motions of the globular cluster 47 Tuc have allowed us to measure for the first time the cluster rotation in the plane of the sky and the velocity anisotropy profile from the cluster core out to about 13'. These profiles are coupled with prior measurements along the line of sight and the surface-brightness profile, and fit all together with self-consistent models specifically constructed to describe quasi-relaxed stellar systems with realistic differential rotation, axisymmetry and pressure anisotropy. The best-fit model provides an inclination angle i between the rotation axis and the line-of-sight direction of 30 deg, and is able to simultaneously reproduce the full three-dimensional kinematics and structure of the cluster, while preserving a good agreement with the projected morphology. Literature models based solely on line-of-sight measurements imply a significantly different inclination angle (i=45 deg), demonstrating that proper motions play a key role in constraining the intrinsic structure of 47 Tuc. Our best-fit global dynamical model implies an internal rotation higher than previous studies have shown, and suggests a peak of the intrinsic V/sigma ratio of ~0.9 at around two half-light radii, with a non-monotonic intrinsic ellipticity profile reaching values up to 0.45. Our study unveils a new degree of dynamical complexity in 47 Tuc, which may be leveraged to provide new insights into the formation and evolution of globular clusters.

We take advantage of the exquisite quality of the Hubble Space Telescope 26-filter astro-photometric catalog of the core of Omega Cen presented in the first paper of this series and the empirical differential-reddening correction presented in the second paper in order to distill the main sequence into its constituent populations. To this end, we restrict ourselves to the five most useful filters: the magic "trio" of F275W, F336W, and F438W, along with F606W and F814W. We develop a strategy for identifying color systems where different populations stand out most distinctly, then we isolate those populations and examine them in other filters where their sub-populations also come to light. In this way, we have identified at least 15 sub-populations, each of which has a distinctive fiducial curve through our 5-dimensional photometric space. We confirm the MSa to be split into two subcomponents, and find that both the bMS and the rMS are split into three subcomponents. Moreover, we have discovered two additional MS groups: the MSd (which has three subcomponents) shares similar properties with the bMS, and the MSe (which has four subcomponents), has properties more similar to those of the rMS. We examine the fiducial curves together and use synthetic spectra to infer relative heavy-element, light-element, and Helium abundances for the populations. Our findings show that the stellar populations and star formation history of Omega Cen are even more complex than inferred previously. Finally, we provide as a supplement to the original catalog a list that identifies for each star which population it most likely is associated with.

We present a pilot study of Galactic globular cluster (GC) proper motion (PM) determinations using Gaia data. We search for GC stars in the Tycho-Gaia Astrometric Solution (TGAS) catalogue from Gaia Data Release 1 (DR1), and identify five members of NGC104 (47 Tucanae), one member of NGC5272 (M3), five members of NGC6121 (M4), seven members of NGC6397, and two members of NGC6656 (M22). By taking a weighted average of member stars, fully accounting for the correlations between parameters, we estimate the parallax (and, hence, distance) and PM of the GCs. This provides a homogeneous PM study of multiple GCs based on an astrometric catalog with small and well-controlled systematic errors, and yields random PM errors similar to existing measurements. Detailed comparison to the available Hubble Space Telescope (HST) measurements generally shows excellent agreement, validating the astrometric quality of both TGAS and HST. By contrast, comparison to ground-based measurements shows that some of those must have systematic errors exceeding the random errors. Our parallax estimates have uncertainties an order of magnitude larger than previous studies, but nevertheless imply distances consistent with previous estimates. By combining our PM measurements with literature positions, distances, and radial velocities, we measure Galactocentric space motions for the clusters and find that these also agree well with previous analyses. Our analysis provides a framework for determining more accurate distances and PMs of Galactic GCs using future Gaia data releases. This will provide crucial constraints on the near end of the cosmic distance ladder and provide accurate GC orbital histories.

We have constructed the most-comprehensive catalog of photometry and proper motions ever assembled for a globular cluster (GC). The core of $\omega$Cen has been imaged over 650 times through WFC3's UVIS and IR channels for the purpose of detector calibration. There exist from 4 to over 60 exposures through each of 26 filters, stretching continuously from F225W in the UV to F160W in the infrared. Furthermore, the 11-year baseline between these data and a 2002 ACS survey has allowed us to more than double the proper-motion accuracy and triple the number of well-measured stars compared to our previous groundbreaking effort. This totally unprecedented complete spectral coverage for over 470,000 stars within the cluster's core, from the tip of the red-giant branch down to the white dwarfs, provides the best astro-photometric observational data base yet to understand the multiple-population phenomenon in any GC. In this first paper of the series we describe in detail the data-reduction processes and deliver the astro-photometric catalog to the astronomical community.

We take advantage of the exquisite quality of the Hubble Space Telescope astro-photometric catalog of the core of wCen presented in the first paper of this series to derive a high-resolution, high-precision, high-accuracy differential-reddening map of the field. The map has a spatial resolution of 2x2 square arcsecs over a total field of view of about 4.3'x4.3'. The differential reddening itself is estimated via an iterative procedure using five distinct color-magnitude diagrams, which provided consistent results to within the 0.1% level. Assuming an average reddening value E(B-V)=0.12, the differential-reddening within the cluster's core can vary by up to +/- 10%, with a typical a standard deviation of about 4%. Our differential-reddening map is made available to the astronomical community in the form of a multi-extension FITS file. This differential-reddening map is essential for a detailed understanding of the multiple stellar populations of wCen, as presented in the next paper in this series. Moreover, it provides unique insight into the level of small spatial-scale extinction variations in the Galactic foreground.

We present a proper motion measurement for the halo globular cluster Pyxis, using HST/ACS data as the first epoch, and GeMS/GSAOI Adaptive Optics data as the second, separated by a baseline of about 5 years. This is both the first measurement of the proper motion of Pyxis and the first calibration and use of Multi-Conjugate Adaptive Optics data to measure an absolute proper motion for a faint, distant halo object. Consequently, we present our analysis of the Adaptive Optics data in detail. We obtain a proper motion of mu_alpha cos(delta)=1.09+/-0.31 mas/yr and mu_delta=0.68+/-0.29 mas/yr. From the proper motion and the line-of-sight velocity we find the orbit of Pyxis is rather eccentric with its apocenter at more than 100 kpc and its pericenter at about 30 kpc. We also investigate two literature-proposed associations for Pyxis with the recently discovered ATLAS stream and the Magellanic system. Combining our measurements with dynamical modeling and cosmological numerical simulations we find it unlikely Pyxis is associated with either system. We examine other Milky Way satellites for possible association using the orbit, eccentricity, metallicity, and age as constraints and find no likely matches in satellites down to the mass of Leo II. We propose that Pyxis probably originated in an unknown galaxy, which today is fully disrupted. Assuming that Pyxis is bound and not on a first approach, we derive a 68% lower limit on the mass of the Milky Way of 0.95*10^12 M_sun.

We search for runaway stars in the Large Magellanic Cloud (LMC) by computing the space velocities of the visually brightest stars in the LMC that are included in the Gaia TGAS proper motion catalog. We compare with predictions from stellar dynamical models to obtain (peculiar) velocities relative to their local stellar environment. Two of the 31 stars have unusually high proper motions. Of the remaining 29 stars we find that most objects in this sample have velocities in very good agreement with model predictions of a circularly rotating disk model. Indeed the excellent fit to the model implies that the TGAS uncertainty estimates are likely overestimated. The fastest outliers in this subsample contain the LBV R71 and a few other well known emission line objects though in no case do we derive velocities consistent with fast (~100 km/s) runaways. Our results imply that R 71 in particular has a moderate deviation from the local stellar velocity field (40 km/s) lending support to the proposition that this object cannot have evolved as a normal single star since it lies too far from massive star forming complexes to have arrived at its current position during its lifetime. Our findings therefore strengthen the case for this LBV being the result of binary evolution. Of the two stars with unusually high proper motions we find that one, the isolated B1.5 Ia+ supergiant Sk-67 2 (HIP 22237), is a candidate hypervelocity star, the TGAS proper motion implying a very large peculiar transverse velocity (~360 km/s) directed radially away from the LMC centre. If confirmed, for example by Gaia Data Release 2, it would imply that this massive supergiant, on the periphery of the LMC, is leaving the galaxy where it will explode as a supernova.

Proper motions (PMs) are crucial to fully understand the internal dynamics of globular clusters (GCs). To that end, the Hubble Space Telescope (HST) Proper Motion (HSTPROMO) collaboration has constructed large, high-quality PM catalogues for 22 Galactic GCs. We highlight some of our exciting recent results: the first directly-measured radial anisotropy profiles for a large sample of GCs; the first dynamical distance and mass-to-light (M/L) ratio estimates for a large sample of GCs; and the first dynamically-determined masses for hundreds of blue-straggler stars (BSSs) across a large GC sample.

The Hubble Space Telescope recently celebrated 25 years of operation. Some of the first images of extragalactic optical jets were taken by HST in the mid-1990s; with time baselines on the order of 20 years and state-of-the-art astrometry techniques, we are now able to reach accuracies in proper-motion measurements on the order of a tenth of a milliarcsecond per year. We present the results of a recent HST program to measure the kiloparsec-scale proper motions of eleven nearby optical jets with Hubble, the first sample of its kind. When paired with VLBI proper-motion measurements on the parsec scale, we are now able to map the full velocity profile of these jets from near the black hole to the final deceleration as they extend out into and beyond the host galaxy. We see convincing evidence that weak-flavor jets (i.e., FR Is) have a slowly increasing jet speed up to 100 pc from the core, where superluminal components are first seen.

The Hubble Space Telescope (HST) UV Legacy Survey of Galactic Globular Clusters (GO-13297) has been specifically designed to complement the existing F606W and F814W observations of the Advanced Camera for Surveys (ACS) Globular Cluster Survey (GO-10775) by observing the most accessible 47 of the previous survey's 65 clusters in three WFC3/UVIS filters F275W, F336W, and F438W. The new survey also adds super-solar metallicity open cluster NGC 6791 to increase the metallicity diversity. The combined survey provides a homogeneous 5-band data set that can be used to pursue a broad range of scientific investigations. In particular, the chosen UV filters allow the identification of multiple stellar populations by targeting the regions of the spectrum that are sensitive to abundance variations in C, N, and O. In order to provide the community with uniform preliminary catalogs, we have devised an automated procedure that performs high-quality photometry on the new UV observations (along with similar observations of seven other programs in the archive). This procedure finds and measures the potential sources on each individual exposure using library point-spread functions and cross-correlates these observations with the original ACS-Survey catalog. The catalog of 57 clusters we publish here will be useful to identify stars in the different stellar populations, in particular for spectroscopic follow-up. Eventually, we will construct a more sophisticated catalog and artificial-star tests based on an optimal reduction of the UV survey data, but the catalogs presented here give the community the chance to make early use of this HST Treasury survey.

We report a new ultra-faint stellar system found in Dark Energy Camera data from the first observing run of the Magellanic Satellites Survey (MagLiteS). MagLiteS J0644-5953 (Pictor II or Pic II) is a low surface brightness ({\mu} = 28.5 mag arcsec$^{-2}$ within its half-light radius) resolved overdensity of old and metal-poor stars located at a heliocentric distance of 45 kpc. The physical size (r$_{1/2}$ = 46 pc) and low luminosity (Mv = -3.2 mag) of this satellite are consistent with the locus of spectroscopically confirmed ultra-faint galaxies. MagLiteS J0644-5953 (Pic II) is located 11.3 kpc from the Large Magellanic Cloud (LMC), and comparisons with simulation results in the literature suggest that this satellite was likely accreted with the LMC. The close proximity of MagLiteS J0644-5953 (Pic II) to the LMC also makes it the most likely ultra-faint galaxy candidate to still be gravitationally bound to the LMC.

We use the Gaia data release 1 (DR1) to study the proper motion (PM) fields of the Large and Small Magellanic Clouds (LMC, SMC). This uses the Tycho-Gaia Astrometric Solution (TGAS) PMs for 29 Hipparcos stars in the LMC and 8 in the SMC. The LMC PM in the West and North directions is inferred to be $(\mu_W,\mu_N) = (-1.872 \pm 0.045, 0.224 \pm 0.054)$ mas/yr, and the SMC PM $(\mu_W,\mu_N) = (-0.874 \pm 0.066, -1.229 \pm 0.047)$ mas/yr. These results have similar accuracy and agree to within the uncertainties with existing Hubble Space Telescope (HST) PM measurements. Since TGAS uses different methods with different systematics, this provides an external validation of both data sets and their underlying approaches. Residual DR1 systematics may affect the TGAS results, but the HST agreement implies this must be below the random errors. Also in agreement with prior HST studies, the TGAS LMC PM field clearly shows the clockwise rotation of the disk, even though it takes the LMC disk in excess of $10^8$ years to complete one revolution. The implied rotation curve amplitude for young LMC stars is consistent with that inferred from line-of-sight (LOS) velocity measurements. Comparison of the PM and LOS rotation curves implies a kinematic LMC distance modulus $m-M = 18.54 \pm 0.39$, consistent but not yet competitive with photometric methods. These first results from Gaia on the topic of Local Group dynamics provide an indication of how its future data releases will revolutionize this field.

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.

Based on color-magnitude diagrams (CMDs) from the Hubble Space Telescope Hubble Tarantula Treasury Project (HTTP) survey, we present the star formation history (SFH) of Hodge~301, the oldest star cluster in the Tarantula Nebula. The HTTP photometry extends faint enough to reach, for the first time, the cluster pre-main sequence (PMS) turn-on, where the PMS joins the main sequence. Using the location of this feature, along with synthetic CMDs generated with the latest PARSEC models, we find that Hodge~301 is older than previously thought, with an age between 26.5 and 31.5 Myr. From this age, we also estimate that between 38 and 61 supernovae Type-II exploded in the region. The same age is derived from the main sequence turn-off, whereas the age derived from the post-main sequence stars is younger and between 20 and 25 Myr. Other relevant parameters are a total stellar mass of $\approx 8800\,\pm 800$M$_{\odot}$ and average reddening E(B$-$V) $\approx 0.22-0.24$ mag, with a differential reddening $\delta$E(B$-$V)$\approx 0.04$ mag.

We use high-precision photometry of red-giant-branch (RGB) stars in 57 Galactic globular clusters (GCs), mostly from the `Hubble Space Telescope (HST) UV Legacy Survey of Galactic globular clusters', to identify and characterize their multiple stellar populations. For each cluster the pseudo two-color diagram (or `chromosome map') is presented, built with a suitable combination of stellar magnitudes in the F275W, F336W, F438W and F814W filters that maximizes the separation between multiple populations. In the chromosome map of most GCs (Type I clusters), stars separate in two distinct groups that we identify with the first (1G) and the second generation (2G). This identification is further supported by noticing that 1G stars have primordial (oxygen-rich, sodium-poor) chemical composition, whereas 2G stars are enhanced in sodium and depleted in oxygen. This 1G-2G separation is not possible for a few GCs where the two sequences have apparently merged into an extended, continuous sequence. In some GCs (Type II clusters) the 1G and/or the 2G sequences appear to be split, hence displaying more complex chromosome maps. These clusters exhibit multiple SGBs also in purely optical color-magnitude diagrams, with the fainter SGB joining into a red RGB which is populated by stars with enhanced heavy-element abundance. We measure the RGB width by using appropriate colors and pseudo-colors. When the metallicity dependence is removed, the RGB width correlates with the cluster mass. The fraction of 1G stars ranges from ~8% to ~67% and anticorrelates with the cluster mass, indicating that incidence and complexity of the multiple population phenomenon both increase with cluster mass.

We present the star formation history of the extremely metal-poor dwarf galaxy DDO 68, based on our photometry with the Advanced Camera for Surveys. With a metallicity of only $12+\log(O/H)=7.15$ and a very isolated location, DDO 68 is one of the most metal-poor galaxies known. It has been argued that DDO 68 is a young system that started forming stars only $\sim 0.15$ Gyr ago. Our data provide a deep and uncontaminated optical color-magnitude diagram that allows us to disprove this hypothesis, since we find a population of at least $\sim 1$ Gyr old stars. The star formation activity has been fairly continuous over all the look-back time. The current rate is quite low, and the highest activity occurred between 10 and 100 Myr ago. The average star formation rate over the whole Hubble time is $\simeq 0.01$ M$_{\odot}$ yr$^{-1}$, corresponding to a total astrated mass of $\simeq 1.3 \times 10^8$ M$_{\odot}$. Our photometry allows us to infer the distance from the tip of the red giant branch, $D = 12.08 \pm 0.67$ Mpc; however, to let our synthetic color-magnitude diagram reproduce the observed ones we need a slightly higher distance, $D=12.65$ Mpc, or $(m-M)_0 = 30.51$, still inside the errors of the previous determination, and we adopt the latter. DDO 68 shows a very interesting and complex history, with its quite disturbed shape and a long Tail probably due to tidal interactions. The star formation history of the Tail differs from that of the main body mainly for an enhanced activity at recent epochs, likely triggered by the interaction.

We make use of the Hubble Space Telescope proper-motion catalogs derived by Bellini et al. (2014) to produce the first radial velocity-dispersion profiles sigma(R) for blue straggler stars (BSSs) in Galactic globular clusters (GCs), as well as the first dynamical estimates for the average mass of the entire BSS population. We show that BSSs typically have lower velocity dispersions than stars with mass equal to the main-sequence turnoff mass, as one would expect for a more massive population of stars. Since GCs are expected to experience some degree of energy equipartition, we use the relation sigma~M^-eta, where eta is related to the degree of energy equipartition, along with our velocity-dispersion profiles to estimate BSS masses. We estimate eta as a function of cluster relaxation from recent Monte Carlo cluster simulations by Bianchini et al. (2016b) and then derive an average mass ratio M_BSS /M_MSTO=1.50+/-0.14 and an average mass M_BSS=1.22+/-0.12 M_Sun from 598 BSSs across 19 GCs. The final error bars include any systematic errors that are random between different clusters, but not any potential biases inherent to our methodology. Our results are in good agreement with the average mass of M_BSS=1.22+/-0.06 M_Sun for the 35 BSSs in Galactic GCs in the literature with properties that have allowed individual mass determination.

High-precision kinematic studies of globular clusters require an accurate knowledge of all possible sources of contamination. Amongst other sources, binary stars can introduce systematic biases in the kinematics. Using a set of Monte Carlo cluster simulations with different concentrations and binary fractions, we investigate the effect of unresolved binaries on proper-motion dispersion profiles, treating the simulations like HST proper-motion samples. Since globular clusters evolve towards a state of partial energy equipartition, more massive stars lose energy and decrease their velocity dispersion. As a consequence, on average, binaries have a lower velocity dispersion, since they are more massive kinematic tracers. We show that, in the case of clusters with high binary fraction (initial binary fraction of 50%) and high concentration (i.e., closer to energy equipartition), unresolved binaries introduce a color-dependent bias in the velocity dispersion of main-sequence stars of the order of 0.1-0.3 km s$^{-1}$ (corresponding to $1-6$% of the velocity dispersion), with the reddest stars having a lower velocity dispersion, due to the higher fraction of contaminating binaries. This bias depends on the ability to distinguish binaries from single stars, on the details of the color-magnitude diagram and the photometric errors. We apply our analysis to the HSTPROMO data set of NGC 7078 (M15) and show that no effect ascribable to binaries is observed, consistent with the low binary fraction of the cluster. Our work indicates that binaries do not significantly bias proper-motion velocity-dispersion profiles, but should be taken into account in the error budget of kinematic analyses.

We report on the study of interstellar extinction across the Tarantula nebula (30 Doradus), in the Large Magellanic Cloud, using observations from the Hubble Tarantula Treasury Project in the 0.3 - 1.6 micron range. The considerable and patchy extinction inside the nebula causes about 3500 red clump stars to be scattered along the reddening vector in the colour-magnitude diagrams, thereby allowing an accurate determination of the reddening slope in all bands. The measured slope of the reddening vector is remarkably steeper in all bands than in the the Galactic diffuse interstellar medium. At optical wavelengths, the larger ratio of total-to-selective extinction, namely Rv = 4.5 +/- 0.2, implies the presence of a grey component in the extinction law, due to a larger fraction of large grains. The extra large grains are most likely ices from supernova ejecta and will significantly alter the extinction properties of the region until they sublimate in 50 - 100 Myr. We discuss the implications of this extinction law for the Tarantula nebula and in general for regions of massive star formation in galaxies. Our results suggest that fluxes of strongly star forming regions are likely to be underestimated by a factor of about 2 in the optical.

We report the discovery of one RR Lyrae star in the ultra--faint satellite galaxy Hydra II based on time series photometry in the g, r and i bands obtained with the Dark Energy Camera at Cerro Tololo Interamerican Observatory, Chile. The RR Lyrae star has a mean magnitude of $i = 21.30\pm 0.04$ which translates to a heliocentric distance of $151\pm 8$ kpc for Hydra II; this value is $\sim 13\%$ larger than the estimate from the discovery paper based on the average magnitude of several blue horizontal branch star candidates. The new distance implies a slightly larger half-light radius of $76^{+12}_{-10}$ pc and a brighter absolute magnitude of $M_V = -5.1 \pm 0.3$, which keeps this object within the realm of the dwarf galaxies. The pulsational properties of the RR Lyrae star ($P=0.645$ d, $\Delta g = 0.68$ mag) suggest Hydra II may be a member of the intermediate Oosterhoff or Oosterhoff II group. A comparison with other RR Lyrae stars in ultra--faint systems indicates similar pulsational properties among them, which are different to those found among halo field stars and those in the largest of the Milky Way satellites. We also report the discovery of 31 additional short period variables in the field of view (RR Lyrae, SX Phe, eclipsing binaries, and a likely anomalous cepheid). However, given their magnitudes and large angular separation from Hydra II, they must be field stars not related to Hydra II.

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.

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 $\beta$ using 3D kinematics outside of the solar neighborhood. We find $\beta=-0.3^{+0.4}_{-0.9}$, consistent with isotropy and lower than solar neighborhood $\beta$ measurements by 2$\sigma$ ($\beta_{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 $\beta=0.1^{+0.4}_{-1.0}$, which is still lower than solar neighborhood measurements by $1\sigma$. The potential decrease in $\beta$ 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.

The radio galaxy 3C 273 hosts one of the nearest and best-studied powerful quasar jets. Having been imaged repeatedly by the Hubble Space Telescope (HST) over the past twenty years, it was chosen for an HST program to measure proper motions in the kiloparsec-scale resolved jets of nearby radio-loud active galaxies. The jet in 3C 273 is highly relativistic on sub-parsec scales, with apparent proper motions up to 15$c$ observed by VLBI (Lister et al., 2013). In contrast, we find that the kpc-scale knots are compatible with being stationary, with a mean speed of $-$0.2$\pm$0.5$c$ over the whole jet. Assuming the knots are packets of moving plasma, an upper limit of 1c implies a bulk Lorentz factor $\Gamma<$2.9. This suggests that the jet has either decelerated significantly by the time it reaches the kpc scale, or that the knots in the jet are standing shock features. The second scenario is incompatible with the inverse Compton off the Cosmic Microwave Background (IC/CMB) model for the X-ray emission of these knots, which requires the knots to be in motion, but IC/CMB is also disfavored in the first scenario due to energetic considerations, in agreement with the recent finding of Meyer & Georganopoulos (2014) which ruled out the IC/CMB model for the X-ray emission of 3C 273 via gamma-ray upper limits.

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.

We present in this study a first analysis of the astrometric error budget of absolute astrometry relative to background galaxies using adaptive optics. We use for this analysis multi-conjugated adaptive optics (MCAO) images obtained with GeMS/GSAOI at Gemini South. We find that it is possible to obtain 0.3 mas reference precision in a random field with 1 hour on source using faint background galaxies. Systematic errors are correctable below that level, such that the overall error is approximately 0.4 mas. Because the reference sources are extended, we find it necessary to correct for the dependency of the PSF centroid on the used aperture size, which would otherwise cause an important bias. This effect needs also to be considered for Extremely Large Telescopes (ELTs). When this effect is corrected, ELTs have the potential to measure proper motions of dwarfs galaxies around M31 with 10 km/s accuracy over a baseline of 5 years.

The coronagraphic instrument currently proposed for the WFIRST-AFTA mission will be the first example of a space-based coronagraph optimized for extremely high contrasts that are required for the direct imaging of exoplanets reflecting the light of their host star. While the design of this instrument is still in progress, this early stage of development is a particularly beneficial time to consider the operation of such an instrument. In this paper, we review current or planned operations on the Hubble Space Telescope (HST) and the James Webb Space Telescope (JWST) with a focus on which operational aspects will have relevance to the planned WFIRST-AFTA coronagraphic instrument. We identify five key aspects of operations that will require attention: 1) detector health and evolution, 2) wavefront control, 3) observing strategies/post-processing, 4) astrometric precision/target acquisition, and 5) polarimetry. We make suggestions on a path forward for each of these items.

We present and describe the astro-photometric catalog of more than 800,000 sources found in the Hubble Tarantula Treasury Project (HTTP). HTTP is a Hubble Space Telescope (HST) Treasury program designed to image the entire 30 Doradus region down to the sub-solar (~0.5 solar masses) mass regime using the Wide Field Camera 3 (WFC3) and the Advanced Camera for Surveys (ACS). We observed 30 Doradus in the near ultraviolet (F275W, F336W), optical (F555W, F658N, F775W), and near infrared (F110W, F160W) wavelengths. The stellar photometry was measured using point-spread function (PSF) fitting across all the bands simultaneously. The relative astrometric accuracy of the catalog is 0.4 mas. The astro-photometric catalog, results from artificial star experiments and the mosaics for all the filters are available for download. Color-magnitude diagrams are presented showing the spatial distributions and ages of stars within 30 Dor as well as in the surrounding fields. HTTP provides the first rich and statistically significant sample of intermediate and low mass pre-main sequence candidates and allows us to trace how star formation has been developing through the region. The depth and high spatial resolution of our analysis highlight the dual role of stellar feedback in quenching and triggering star formation on the giant HII region scale. Our results are consistent with stellar sub-clustering in a partially filled gaseous nebula that is offset towards our side of the Large Magellanic Cloud.

We build on the evidence provided by our Legacy Survey of Galactic globular clusters (GC) to submit to a crucial test four scenarios currently entertained for the formation of multiple stellar generations in GCs. The observational constraints on multiple generations to be fulfilled are manifold, including GC specificity, ubiquity, variety, predominance, discreteness, supernova avoidance, p-capture processing, helium enrichment and mass budget. We argue that scenarios appealing to supermassive stars, fast rotating massive stars and massive interactive binaries violate in an irreparable fashion two or more among such constraints. Also the scenario appealing to AGB stars as producers of the material for next generation stars encounters severe difficulties, specifically concerning the mass budget problem and the detailed chemical composition of second generation stars. We qualitatively explore ways possibly allowing one to save the AGB scenario, specifically appealing to a possible revision of the cross section of a critical reaction rate destroying sodium, or alternatively by a more extensive exploration of the vast parameter space controlling the evolutionary behavior of AGB stellar models. Still, we cannot ensure success for these efforts and totally new scenarios may have to be invented to understand how GCs formed in the early Universe.

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.

Numerous observational studies have revealed the ubiquitous presence of multiple stellar populations in globular clusters and cast many hard challenges for the study of the formation and dynamical history of these stellar systems. In this Letter we present the results of a study of the kinematic properties of multiple populations in NGC 2808 based on high-precision Hubble Space Telescope proper-motion measurements. In a recent study, Milone et al. have identified five distinct populations (A, B, C, D, and E) in NGC 2808. Populations D and E coincide with the helium-enhanced populations in the middle and the blue main sequences (mMS and bMS) previously discovered by Piotto et al.; populations A, B, and C correspond to the redder main sequence (rMS) that in the Piotto et al. was associated with the primordial stellar population. Our analysis shows that, in the outermost regions probed (between about 1.5 and 2 times the cluster half-light radius), the velocity distribution of populations D and E is radially anisotropic (the deviation from an isotropic distribution is significant at the ~3.5-sigma level). Stars of populations D and E have a smaller tangential velocity dispersion than those of populations A, B, and C, while no significant differences are found in the radial-velocity dispersion. We present the results of a numerical simulation showing that the observed differences between the kinematics of these stellar populations are consistent with the expected kinematic fingerprint of the diffusion towards the cluster outer regions of stellar populations initially more centrally concentrated.

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.

We present measurements of positions and relative proper motions in the 30 Doradus region of the Large Magellanic Cloud (LMC). We detail the construction of a single-epoch astrometric reference frame, based on specially-designed observations obtained with the two main imaging instruments ACS/WFC and WFC3/UVIS onboard the Hubble Space Telescope (HST). Internal comparisons indicate a sub milli-arc-second (mas) precision in the positions and the presence of semi-periodic systematics with a mean amplitude of ~0.8 mas. We combined these observations with numerous archival images taken with WFPC2 and spanning 17 years. The precision of the resulting proper motions for well-measured stars around the massive cluster R 136 can be as good as ~20 microarcsec/yr, although the true accuracy of proper motions is generally lower due to the residual systematic errors. The observed proper-motion dispersion for our highest-quality measurements is ~0.1 mas/yr. Our catalog of positions and proper motions contains 86,590 stars down to V~25 and over a total area of ~70 square arcmin. We examined the proper motions of 105 relatively bright stars and identified a total of 6 candidate runaway stars. We are able to tentatively confirm the runaway status of star VFTS 285, consistent with the findings from line-of-sight velocities, and to show that this star has likely been ejected from R 136. This study demonstrates that with HST it is now possible to reliably measure proper motions of individual stars in the nearest dwarf galaxies such as the LMC.

We present a study of the recent star formation of 30 Doradus in the Large Magellanic Cloud (LMC) using the panchromatic imaging survey Hubble Tarantula Treasury Project (HTTP). In this paper we focus on the stars within 20 pc of the center of the massive ionizing cluster of 30 Doradus, NGC 2070. We recovered the star formation history by comparing deep optical and NIR color-magnitude diagrams (CMDs) with state-of-the-art synthetic CMDs generated with the latest PARSEC models, which include all stellar phases from pre-main sequence to post- main sequence. For the first time in this region we are able to measure the star formation using intermediate and low mass stars simultaneously. Our results suggest that NGC2070 experienced a prolonged activity. In particular, we find that the star formation in the region: i) exceeded the average LMC rate ~ 20 Myr ago; ii) accelerated dramatically ~ 7 Myr ago; and iii) reached a peak value 1-3 Myr ago. We did not find significant deviations from a Kroupa initial mass function down to 0.5 Msun. The average internal reddening E(B-V) is found to be between 0.3 and 0.4 mag.

In this paper we combine WFC3/UVIS F275W, F336W, and F438W data from the "UV Legacy Survey of Galactic Globular Clusters: Shedding Light on Their Populations and Formation" (GO-13297) HST Treasury program with F606W, F625W, F658N, and F814W ACS archive data for a multi-wavelength study of the globular cluster NGC 6352. In the color-magnitude and two-color diagrams obtained with appropriate combination of the photometry in the different bands we separate two distinct stellar populations and trace them from the main sequence to the subgiant, red giant, horizontal and asymptotic giant branches. We infer that the two populations differ in He by Delta Y=0.029+/-0.006. With a new method, we also estimate the age difference between the two sequences. Assuming no difference in [Fe/H] and [alpha/Fe], and the uncertainties on Delta Y, we found a difference in age between the two populations of 10+/-120 Myr. If we assume [Fe/H] and [alpha/Fe] differences of 0.02 dex (well within the uncertainties of spectroscopic measurements), the total uncertainty in the relative age rises to ~300 Myr.

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).

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.

We present the deepest optical photometry for any dwarf elliptical (dE) galaxy based on Hubble Space Telescope ACS observations of the Local Group dE galaxies NGC 147 and NGC 185. The resulting F606W and F814W color-magnitude diagrams are the first to reach below the main sequence turnoff in a dE galaxy, allowing us to determine full star formation histories in these systems. The ACS fields are located ~1.5 effective radii from the galaxy center to avoid photometric crowding. While our ACS pointings in both dEs show unambiguous evidence for old and intermediate age stars, the mean age in NGC 147 is ~ 4 Gyr younger as compared to NGC 185. In NGC 147, only 40% of stars were in place 12.5 Gyrs ago (z~5), with the bulk of the remaining stellar population forming between 5 to 7 Gyr. In contrast, 70% of stars were formed in NGC 185 field more than12.5 Gyr ago with the majority of the remaining population forming between 8 to 10 Gyr. Star formation ceased in both ACS fields at least 3 Gyr ago. Previous observations in the central regions of NGC 185 show evidence for star formation as recent as 100 Myr ago and a strong metallicity gradient with radius. We suggest that the orbit of NGC 185 has a larger pericenter as compared to NGC 147, allowing it to preserve radial gradients and maintain a small central reservoir of recycled gas. We interpret the inferred differences in star formation histories to imply an earlier infall time into the M31 environment for NGC 185 as compared to NGC 147.

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\sigma 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% (\sigma_t/\sigma_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 (\sigma_t/\sigma_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 \sigma_minor/\sigma_major and ellipticity, such that the more flattened clusters in the sample tend to be more anisotropic, with \sigma_minor/\sigma_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.

We present high-precision multi-band photometry for the globular cluster (GC) M2. We combine the analysis of the photometric data obtained from the Hubble Space Telescope UV Legacy Survey of Galactic GCs GO-13297, with chemical abundances by Yong et al.(2014), and compare the photometry with models in order to analyze the multiple stellar sequences we identified in the color-magnitude diagram (CMD). We find three main stellar components, composed of metal-poor, metal-intermediate, and metal-rich stars (hereafter referred to as population A, B, and C, respectively). The components A and B include stars with different $s$-process element abundances. They host six sub-populations with different light-element abundances, and exhibit an internal variation in helium up to Delta Y~0.07 dex. In contrast with M22, another cluster characterized by the presence of populations with different metallicities, M2 contains a third stellar component, C, which shows neither evidence for sub-populations nor an internal spread in light-elements. Population C does not exhibit the typical photometric signatures that are associated with abundance variations of light elements produced by hydrogen burning at hot temperatures. We compare M2 with other GCs with intrinsic heavy-element variations and conclude that M2 resembles M22, but it includes an additional stellar component that makes it more similar to the central region of the Sagittarius galaxy, which hosts a GC (M54) and the nucleus of the Sagittarius galaxy itself.

We present the first study of high-precision internal proper motions (PMs) in a large sample of globular clusters, based on Hubble Space Telescope (HST) data obtained over the past decade with the ACS/WFC, ACS/HRC, and WFC3/UVIS instruments. We determine PMs for over 1.3 million stars in the central regions of 22 clusters, with a median number of ~60,000 stars per cluster. These PMs have the potential to significantly advance our understanding of the internal kinematics of globular clusters by extending past line-of-sight (LOS) velocity measurements to two- or three-dimensional velocities, lower stellar masses, and larger sample sizes. We describe the reduction pipeline that we developed to derive homogeneous PMs from the very heterogeneous archival data. We demonstrate the quality of the measurements through extensive Monte-Carlo simulations. We also discuss the PM errors introduced by various systematic effects, and the techniques that we have developed to correct or remove them to the extent possible. We provide in electronic form the catalog for NGC 7078 (M 15), which consists of 77,837 stars in the central 2.4 arcmin. We validate the catalog by comparison with existing PM measurements and LOS velocities, and use it to study the dependence of the velocity dispersion on radius, stellar magnitude (or mass) along the main sequence, and direction in the plane of the sky (radial/tangential). Subsequent papers in this series will explore a range of applications in globular-cluster science, and will also present the PM catalogs for the other sample clusters.

In this paper we describe a new UV-initiative HST project (GO-13297) that will complement the existing F606W and F814W database of the ACS Globular Cluster (GC) Treasury by imaging most of its clusters through UV/blue WFC3/UVIS filters F275W, F336W and F438W. This "magic trio" of filters has shown an uncanny ability to disentangle and characterize multiple-population (MP) patterns in GCs in a way that is exquisitely sensitive to C, N, and O abundance variations. Combination of these passbands with those in the optical also gives the best leverage for measuring helium enrichment. The dozen clusters that had previously been observed in these bands exhibit a bewildering variety of MP patterns, and the new survey will map the full variance of the phenomenon. The ubiquity of multiple stellar generations in GCs has made the formation of these cornerstone objects more intriguing than ever; GC formation and the origin of their MPs have now become one and the same problem. In the present paper we will describe the data base and our data reduction strategy, as well as the uses we intend to make of the final photometry, astrometry, and proper motions. We will also present preliminary color-magnitude diagrams from the data so far collected. These diagrams also draw on data from GO-12605 and GO-12311, which served as a pilot project for the present GO-13297.

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.

The James Webb Space Telescope (JWST) Optical Simulation Testbed (JOST) is a tabletop workbench to study aspects of wavefront sensing and control for a segmented space telescope, including both commissioning and maintenance activities. JOST is complementary to existing optomechanical testbeds for JWST (e.g. the Ball Aerospace Testbed Telescope, TBT) given its compact scale and flexibility, ease of use, and colocation at the JWST Science & Operations Center. We have developed an optical design that reproduces the physics of JWST's three-mirror anastigmat using three aspheric lenses; it provides similar image quality as JWST (80% Strehl ratio) over a field equivalent to a NIRCam module, but at HeNe wavelength. A segmented deformable mirror stands in for the segmented primary mirror and allows control of the 18 segments in piston, tip, and tilt, while the secondary can be controlled in tip, tilt and x, y, z position. This will be sufficient to model many commissioning activities, to investigate field dependence and multiple field point sensing & control, to evaluate alternate sensing algorithms, and develop contingency plans. Testbed data will also be usable for cross-checking of the WFS&C Software Subsystem, and for staff training and development during JWST's five- to ten-year mission.

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 ($\mu_{\alpha}\cos\delta, \mu_{\delta}$)=($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 ($\mu_{\alpha}\cos\delta, \mu_{\delta})=(-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 ($\mu_{\alpha}\cos\delta, \mu_{\delta})=(-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.

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.

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.

It is widely believed that globular clusters evolve over many two-body relaxation times toward a state of energy equipartition, so that velocity dispersion scales with stellar mass as \sigma ~ m^{-\eta} with \eta = 0.5. We show that this is incorrect, using direct N-body simulations with a variety of realistic IMFs and initial conditions. No simulated system ever reaches a state close to equipartition. Near the center, the luminous main-sequence stars reach a maximum \eta_{max} ~ 0.15 \pm 0.03. At large times, all radial bins convergence on an asymptotic value \eta_{\infty} ~ 0.08 \pm 0.02. The development of this "partial equipartition" is strikingly similar across our simulations, despite the range of initial conditions employed. Compact remnants tend to have higher \eta than main-sequence stars (but still \eta < 0.5), due to their steeper (evolved) mass function. The presence of an intermediate-mass black hole (IMBH) decreases \eta, consistent with our previous findings of a quenching of mass segregation under these conditions. All these results can be understood as a consequence of the Spitzer instability for two-component systems, extended by Vishniac to a continuous mass spectrum. Mass segregation (the tendency of heavier stars to sink toward the core) has often been studied observationally, but energy equipartition has not. Due to the advent of high-quality proper motion datasets from the Hubble Space Telescope, it is now possible to measure \eta. Detailed data-model comparisons open up a new observational window on globular cluster dynamics and evolution. Comparison of our simulations to Omega Cen observations yields good agreement, confirming that globular clusters are not generally in energy equipartition. Modeling techniques that assume equipartition by construction (e.g., multi-mass Michie-King models) are approximate at best.

We report a sensitive X-ray search for the proposed intermediate mass black hole (IMBH) in the massive Galactic cluster, Omega Centauri (NGC 5139). Combining Chandra X-ray Observatory data from Cycles 1 and 13, we obtain a deep (~291 ks) exposure of the central regions of the cluster. We find no evidence for an X-ray point source near any of the cluster's proposed dynamical centers, and place an upper limit on the X-ray flux from a central source of f_X(0.5-7.0 keV) <= 5.0x10^-16 erg cm^-2 s^-1, after correcting for absorption. This corresponds to an unabsorbed X-ray luminosity of L_X(0.5-7.0 keV) <= 1.6x10^30 erg s^-1, for a cluster distance of 5.2 kpc, Galactic column density N_H = 1.2x10^21 cm^-2, and powerlaw spectrum with Gamma = 2.3. If a ~10^4 M_sun IMBH resides in the cluster's core, as suggested by some stellar dynamical studies, its Eddington luminosity would be L_Edd ~10^42 erg s^-1. The new X-ray limit would then establish an Eddington ratio of L_X/L_Edd <~ 10^-12, a factor of ~10 lower than even the quiescent state of our Galaxy's notoriously inefficient supermassive black hole Sgr A*, and imply accretion efficiencies as low as eta <~ 10^-6 - 10^-8. This study leaves open three possibilities: either Omega Cen does not harbor an IMBH or, if an IMBH does exist, it must experience very little or very inefficient accretion.

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.

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.

The Hubble Tarantula Treasury Project (HTTP) is an ongoing panchromatic imaging survey of stellar populations in the Tarantula Nebula in the Large Magellanic Cloud that reaches into the sub-solar mass regime (< 0.5 Mo). HTTP utilizes the capability of HST to operate the Advanced Camera for Surveys (ACS) and the Wide Field Camera 3 (WFC3) in parallel to study this remarkable region in the near-ultraviolet, optical, and near-infrared spectral regions, including narrow band H$\alpha$ images. The combination of all these bands provides a unique multi-band view. The resulting maps of the stellar content of the Tarantula Nebula within its main body provide the basis for investigations of star formation in an environment resembling the extreme conditions found in starburst galaxies and in the early Universe. Access to detailed properties of individual stars allows us to begin to reconstruct the evolution of the stellar skeleton of the Tarantula Nebula over space and time with parcsec-scale resolution. In this first paper we describe the observing strategy, the photometric techniques, and the upcoming data products from this survey and present preliminary results obtained from the analysis of the initial set of near-infrared observations.

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.

We have derived the star formation history (SFH) of the blue compact dwarf galaxy IZw18 through comparison of deep HST/ACS data with synthetic color magnitude diagrams. A statistical analysis was implemented for the identification of the best-fit SFH and relative uncertainties. We confirm that IZw18 is not a truly young galaxy, having started forming stars earlier than ~1 Gyr ago, and possibly at epochs as old as a Hubble time. In IZw18's main body we infer a lower limit of ~2 x 10^{6} M_sun for the mass locked-up in old stars. IZw18' s main body has been forming stars very actively during the last ~10 Myr, with an average star formation rate (SFR) as high as ~1 M_sun/yr (or ~2 x 10^{-5} M_sun yr^{-1} pc^{-2}). On the other hand, the secondary body was much less active at these epochs, in agreement with the absence of significant nebular emission. The high current SFR can explain the very blue colors and the high ionized gas content in IZw18, resembling primeval galaxies in the early Universe. Detailed chemical evolution models are required to quantitatively check whether the SFH from the synthetic CMDs can explain the low measured element abundances, or if galactic winds with loss of metals are needed.

Based on long baseline (5-7 years) multi-epoch HST/ACS photometry, used previously to measure the proper motion of M31, we present the proper motions (PMs) of 13 main-sequence Milky Way halo stars. The sample lies at an average distance of r ~24 kpc from the Galactic center, with a root-mean-square spread of 6 kpc. At this distance, the median PM accuracy is 5 km/s. We devise a maximum likelihood routine to determine the tangential velocity ellipsoid of the stellar halo. The velocity second moments in the directions of the Galactic (l,b) system are < vl^2 >^{1/2} = 123 (+29, -23) km/s, and < vb^2 >^{1/2} = 83 (+24, -16) km/s. We combine these results with the known line-of-sight second moment, < vlos^2 >^{1/2} = 105 \pm 5$ km/s, at this < r > to study the velocity anisotropy of the halo. We find approximate isotropy between the radial and tangential velocity distributions, with anisotropy parameter beta = 0.0 (+0.2, -0.4). Our results suggest that the stellar halo velocity anisotropy out to r ~ 30 kpc is less radially biased than solar neighborhood measurements. This is opposite to what is expected from violent relaxation, and may indicate the presence of a shell-type structure at r ~ 24 kpc. With additional multi-epoch HST data, the method presented here has the ability to measure the transverse kinematics of the halo for more stars, and to larger distances. This can yield new improved constraints on the stellar halo formation mechanism, and the mass of the Milky Way.

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

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.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

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)

We study the future orbital evolution and merging of the MW-M31-M33 system, using a combination of collisionless N-body simulations and semi-analytic orbit integrations. Monte-Carlo simulations are used to explore the consequences of varying the initial phase-space and mass parameters within their observational uncertainties. The observed M31 transverse velocity implies that the MW and M31 will merge t = 5.86 (+1.61-0.72) Gyr from now, after a first pericenter at t = 3.87 (+0.42-0.32) Gyr. M31 may (probability p=41%) make a direct hit with the MW (defined here as a first-pericenter distance less than 25 kpc). Most likely, the MW and M31 will merge first, with M33 settling onto an orbit around them. Alternatively, M33 may make a direct hit with the MW first (p=9%), or M33 may get ejected from the Local Group (p=7%). The MW-M31 merger remnant will resemble an elliptical galaxy. The Sun will most likely (p=85%) end up at larger radius from the center of the MW-M31 merger remnant than its current distance from the MW center, possibly further than 50 kpc (p=10%). The Sun may (p=20%) at some time in the next 10 Gyr find itself moving through M33 (within 10 kpc), but while dynamically still bound to the MW-M31 merger remnant. The arrival and possible collision of M31 (and possibly M33) with the MW is the next major cosmic event affecting the environment of our Sun and solar system that can be predicted with some certainty. (Abridged)

We determine the velocity vector of M31 with respect to the Milky Way and use this to constrain the mass of the Local Group, based on HST proper-motion measurements presented in Paper I. We construct N-body models for M31 to correct the measurements for the contributions from stellar motions internal to M31. We also estimate the center-of-mass motion independently, using the kinematics of satellite galaxies of M31 and the Local Group. All estimates are mutually consistent, and imply a weighted average M31 heliocentric transverse velocity of (v_W,v_N) = (-125.2+/-30.8, -73.8+/-28.4) km/s. We correct for the reflex motion of the Sun using the most recent insights into the solar motion within the Milky Way. This implies a radial velocity of M31 with respect to the Milky Way of V_rad = -109.3+/-4.4 km/s, and a tangential velocity V_tan = 17.0 km/s (<34.3 km/s at 1-sigma confidence). Hence, the velocity vector of M31 is statistically consistent with a radial (head-on collision) orbit towards the Milky Way. We revise prior estimates for the Local Group timing mass, including corrections for cosmic bias and scatter. Bayesian combination with other mass estimates yields M_LG = M_MW(vir) + M_M31(vir) = (3.17 +/- 0.57) x 10^12 solar masses. The velocity and mass results imply at 95% confidence that M33 is bound to M31, consistent with expectation from observed tidal deformations. (Abridged)

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.

(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.

The cosmological distance ladder crucially depends on classical Cepheids (with P=3-80 days), which are primary distance indicators up to 33 Mpc. Within this volume, very few SNe Ia have been calibrated through classical Cepheids, with uncertainty related to the non-linearity and the metallicity dependence of their period-luminosity (PL) relation. Although a general consensus on these effects is still not achieved, classical Cepheids remain the most used primary distance indicators. A possible extension of these standard candles to further distances would be important. In this context, a very promising new tool is represented by the ultra-long period (ULP) Cepheids (P \geq 80 days), recently identified in star-forming galaxies. Only a small number of ULP Cepheids have been discovered so far. Here we present and analyse the properties of an updated sample of 37 ULP Cepheids observed in galaxies within a very large metallicity range of 12+log(O/H) from ~7.2 to 9.2 dex. We find that their location in the colour(V-I)-magnitude diagram as well as their Wesenheit (V-I) index-period (WP) relation suggests that they are the counterparts at high luminosity of the shorter-period (P \leq 80 days) classical Cepheids. However, a complete pulsation and evolutionary theoretical scenario is needed to properly interpret the true nature of these objects. We do not confirm the flattening in the studied WP relation suggested by Bird et al. (2009). Using the whole sample, we find that ULP Cepheids lie around a relation similar to that of the LMC, although with a large spread (~0.4 mag).

As part of the SPLASH survey of the Andromeda galaxy (M31) and its neighbors, we have obtained Keck/DEIMOS spectra of the compact elliptical (cE) satellite M32. This is the first resolved-star kinematical study of any cE galaxy. In contrast to previous studies that extended out to r<30"~1Re~100pc, we measure the rotation curve and velocity dispersion profile out to r~250" and higher order Gauss-Hermite moments out to r~70". We achieve this by combining integrated-light spectroscopy at small radii (where crowding/blending are severe) with resolved stellar spectroscopy at larger radii, using spatial and kinematical information to statistically account for M31 contamination. The rotation curve and velocity dispersion profile extend well beyond the radius (r~150") where the isophotes are distorted. Unlike NGC 205, another close dwarf companion of M31, M32's kinematic are regular and symmetric and do not show obvious sharp gradients across the region of isophotal elongation and twists. We interpret M32's kinematics using three-integral axisymmetric dynamical equilibrium models constructed using Schwarzschild's orbit superposition technique. Models with a constant M/L can fit the data remarkably well. However, since such a model requires an increasing tangential anisotropy with radius, invoking the presence of an extended dark halo may be more plausible. Such an extended dark halo is definitely required to bind a half-dozen fast-moving stars observed at the largest radii, but these stars may not be an equilibrium component of M32. The observed regularity of the stellar kinematics, as well as the possible detection of an extended dark halo, are unexpected if M31 tides are significant at large radii. While these findings by themselves do not rule out tidal models for cE formation, they suggest that tidal stripping may not be as significant for shaping cE galaxies as has often been argued.

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.

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.

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 (\epsilon~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 \omega Cen (Milky Way) and G1 (M31).

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.

The evolutionary properties and spatial distribution of I Zwicky 18 stellar populations are analyzed by means of HST/ACS deep and accurate photometry. The comparison of the resulting Colour- Magnitude diagrams with stellar evolution models indicates that stars of all ages are present in all the system components, including ob jects possibly up to 13 Gyr old, intermediate age stars and very young ones. The Colour-Magnitude diagrams show evidence of thermally pulsing asymptotic giant branch and carbon stars. Classical and ultra-long period Cepheids, as well as long period variables have been measured. About 20 ob jects could be unresolved star clusters, and are mostly concentrated in the North-West (NW) portion of the Main Body (MB). If interpreted with simple stellar population models, these ob jects indicate a particularly active star formation over the past hundred Myr in IZw 18. The stellar spatial distribution shows that the younger ones are more centrally concentrated, while old and intermediate age stars are distributed homogeneously over the two bodies, although more easily detectable at the system periphery. The oldest stars are best visible in the Secondary Body (SB) and in the South East (SE) portion of the MB, where crowding is less severe, but are present also in the rest of the MB, although measured with larger uncertainties. The youngest stars are a few Myr old, are located predominantly in the MB and mostly concentrated in its NW portion. The SE portion of the MB appears to be in a similar, but not as young evolutionary stage as the NW, while the SB stars are older than at least 10 Myr. There is then a sequence of decreasing age of the younger stars from the Secondary Body to the SE portion of the MB to the NW portion. All our results suggest that IZw18 is not atypical compared to other BCDs.

[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.

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.

Classical Cepheids are primary distance indicators playing a fundamental role in the calibration of the extragalactic distance scale. The possible dependence of their characteristic Period-Luminosity (PL) relation on chemical composition is still debated in the literature, and the behaviour of these pulsators at very low metallicity regimes is almost unexplored. In order to derive constraints on the application of the Period-Luminosity relation at low metal abundances, we investigate the properties of the few ultra-low metallicity (Z ~ 0.0004) Cepheids recently discovered in the Blue Compact Dwarf galaxy IZw18. To this purpose we have computed an updated and extended set of nonlinear convective models for Z=0.0004 and Y=0.24, spanning a wide range of stellar masses, and taking into account the evolutionary constraints for selected luminosity levels. As a result we are able to predict the topology of the instability strip, the variations of all the relevant quantities along the pulsation cycle, including the morphology of the light curves, the theoretical Period-Luminosity-Color, Period-Wesenheit and Period-Luminosity relations at such a low metallicity. For each of these relations we provide the appropriate coefficients for fundamental mode pulsators with Z=0.0004. By comparing these results with the properties of more metal rich Cepheids we find that the synthetic PL relations for Z=0.0004 are steeper than at higher Z, but similar to the Z=0.004 ones, thus suggesting a leveling off of the metallicity effect towards the lowest Zs.

Variable stars have been identified for the first time in the very metal-poor Blue Compact dwarf galaxy IZw18, using deep multi-band (F606W, F814W)time-series photometry obtained with the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope (HST). We detected 34 candidate variable stars in the galaxy. We classify three of them as Classical Cepheids, with periods of 8.71, 125.0 and 130.3 days, respectively, and other two as long period variables with periodicities longer than a hundred days. These are the lowest metallicity Classical Cepheids known so far, thus providing the opportunity to explore and fit models of stellar pulsation fo Classical Cepheids at previously inaccessible metallicities. The period distribution of the confirmed Cepheids is markedly different from what is seen in other nearby galaxies, which is likely related to the star bursting nature of IZw18. By applying to the 8.71 days Cepheid theoretical Wesenheit (V,I) relations based on new pulsation models of Classical Cepheids specifically computed for the extremely low metallicity of this galaxy (Z=0.0004, Y=0.24), we estimate the distance modulus of IZw18 to be mu_0= 31.4pm0.2 D=19.0^{+1.8}_{-1.7}Mpc) for canonical models of Classical Cepheids, and of 31.2pm0.2 mag (D=17.4^{+1.6}_{-1.6}Mpc) using over luminous models. The theoretical modeling of the star's light curves provides mu_0=31.4pm0.1 mag, D=19.0^{+0.9}_{-0.9} Mpc, in good agreement with the results from the theoretical Wesenheit relations. These pulsation distances bracket the distance of 18.2pm1.5Mpc inferred by Aloisi et al. (2007) using the galaxy's Red Giant Branch Tip.

We present kinematic and metallicity profiles for the M31 dwarf elliptical (dE) satellite galaxies NGC 147 and NGC 185. The profiles represent the most extensive spectroscopic radial coverage for any dE galaxy, extending to a projected distance of eight half-light radii (8 r_eff = 14'). We achieve this coverage via Keck/DEIMOS multislit spectroscopic observations of 520 and 442 member red giant branch stars in NGC 147 and NGC 185, respectively. In contrast to previous studies, we find that both dEs have significant internal rotation. We measure a maximum rotational velocity of 17+/-2 km/s for NGC 147 and 15+/-5 km/s for NGC 185. The velocity dispersions decrease gently with radius with an average dispersion of 16+/-1 km/s for NGC 147 and 24+/-1 km/s for NGC 185. Both dEs have internal metallicity dispersions of 0.5 dex, but show no evidence for a radial metallicity gradient. We construct two-integral axisymmetric dynamical models and find that the observed kinematical profiles cannot be explained without modest amounts of non-baryonic dark matter. We measure central mass-to-light ratios of ML_V = 4.2+/-0.6 and ML_V = 4.6+/-0.6 for NGC 147 and NGC 185, respectively. Both dE galaxies are consistent with being primarily flattened by their rotational motions, although some anisotropic velocity dispersion is needed to fully explain their observed shapes. The velocity profiles of all three Local Group dEs (NGC 147, NGC 185 and NGC 205) suggest that rotation is more prevalent in the dE galaxy class than previously assumed, but is often manifest only at several times the effective radius. Since all dEs outside the Local Group have been probed to only inside the effective radius, this opens the door for formation mechanisms in which dEs are transformed or stripped versions of gas-rich rotating progenitor galaxies.

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.

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.

Establishing or ruling out, either through solid mass measurements or upper limits, the presence of intermediate-mass black holes (IMBHs) at the centers of star clusters would profoundly impact our understanding of problems ranging from the formation and long-term dynamical evolution of stellar systems, to the nature of the seeds and the growth mechanisms of supermassive black holes. While there are sound theoretical arguments both for and against their presence in today's clusters, observational studies have so far not yielded truly conclusive IMBH detections nor upper limits. We argue that the most promising approach to solving this issue is provided by the combination of measurements of the proper motions of stars at the centers of Galactic globular clusters and dynamical models able to take full advantage of this type of data set. We present a program based on HST observations and recently developed tools for dynamical analysis designed to do just that.

Ground-based adaptive optics (AO) in the infrared has made exceptional advances in approaching space-like image quality at higher collecting area. Optical-wavelength applications are now also growing in scope. We therefore provide here a comparison of the pros and cons of observational capabilities from the ground and from space at optical wavelengths. With an eye towards the future, we focus on the comparison of a ~30m ground-based telescope with an 8-16m space-based telescope. We review the current state-of-the-art in AO, and summarize the expected future improvements in image quality, field of view, contrast, and low-wavelength cut-off. We discuss the exciting advances in extreme AO for exoplanet studies and explore what the theoretical limitations in achievable contrast might be. Our analysis shows that extreme AO techniques face both fundamental and technological hurdles to reach the contrast of 1E-10 necessary to study an Earth-twin at 10 pc. Based on our assessment of the current state-of-the-art, the future technology developments, and the inherent difficulty of observing through a turbulent atmosphere, we conclude that there will continue to be a strong complementarity between observations from the ground and from space at optical wavelengths in the coming decades. There will continue to be subjects that can only be studied from space, including imaging and (medium-resolution) spectroscopy at the deepest magnitudes, and the exceptional-contrast observations needed to characterize terrestrial exoplanets and search for biomarkers.

[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.

The fact that dark matter (DM), thus far, has revealed itself only on scales of galaxies and larger, again thrusts onto astrophysics the opportunity and the responsibility to confront the age old mystery "What is the nature of matter?" By deriving basic data on the nature of DM - e.g., mass of its particle(s), present mean temperature, distribution in galaxies and other structures in the universe, and capacity for dissipational collapse - we will be uncovering the properties of the dominant species of matter in the universe and significantly extending the standard models of particle physics. Determining the mass of the DM particle to an order of magnitude would help to sort out the particle family to which it (or they) belongs. Beyond mass, there are issues of stability. The DM particle may be unstable with a measurable half-life, or it may become unstable after absorbing a certain amount of energy from collisions. In both cases it would contribute to the present hot dark matter component. Some key parameters of DM can most accurately be measured in the very nearby universe because DM dominates the mass in the outer Milky Way (MW), in other galaxies in the Local Group, and in the Local Group in its entirety. The presence and amount of DM can be quantified by study of dynamical processes observable in fine detail within these entities. Precise measurements of 3-D velocities for stars, coherent star streams, and stars in satellite stellar systems out to the edge of the Galaxy can reveal "what is the shape, orientation, density law, and lumpiness of the dark matter halo" as well as "what is the total mass of the Galaxy?"

The concordance Lambda Cold Dark Matter (Lambda-CDM) model for the formation of structure in the Universe, while remarkably successful at describing observations of structure on large scales, continues to be challenged by observations on galactic scales. Fortunately, CDM models and their various proposed alternatives make a rich variety of testable predictions that make the Local Group and our own Milky Way Galaxy key laboratories for exploring dark matter (DM) in this regime. However, some of the most definitive tests of local DM require microarcsecond astrometry of faint sources, an astrometric regime that is a unique niche of SIM Lite. This chapter explores the important and distinct contributions that can be made by SIM Lite in the exploration of galaxy dynamics and DM on galaxy scales and that have cosmological consequences. Key areas of potential SIM Lite exploration include (1) measuring the shape, orientation, density law, and lumpiness of the dark halo of the Milky Way and other nearby galaxies, (2) determining the orbits of Galactic satellites, which may be representatives of late infall from the hierarchical formation of the Milky Way, (3) ascertaining the distribution of angular momentum and orbital anisotropy of stars and globular clusters to the outer reaches of the Galactic halo, dynamical properties that hold clues to the early hierarchical formation of the Galaxy, (4) measuring the physical nature of DM by placing strong constraints on the phase space density in the cores of nearby dSph galaxies, and (5) reconstructing the dynamical history of the Local Group through the determination of orbits and masses of its constituent galaxies.

In HST Cycles 11 and 13 we obtained two epochs of ACS/HRC data for fields in the Magellanic Clouds centered on background quasars. We used these data to determine the proper motions of the LMC and SMC to better than 5% and 15% respectively. The results had a number of unexpected implications for the Milky Way-LMC-SMC system. The implied three-dimensional velocities were larger than previously believed and close to the escape velocity in a standard 10^12 solar mass Milky Way dark halo, implying that the Clouds may be on their first passage. Also, the relative velocity between the LMC and SMC was larger than expected, leaving open the possibility that the Clouds may not be bound to each other. To further verify and refine our results we requested an additional epoch of data in Cycle 16 which is being executed with WFPC2/PC due to the failure of ACS. We present the results of an ongoing analysis of these WFPC2 data which indicate good consistency with the two-epoch results.

We review our understanding of the kinematics of the LMC and the SMC, and their orbit around the Milky Way. The line-of-sight velocity fields of both the LMC and SMC have been mapped with high accuracy using thousands of discrete traces, as well as HI gas. The LMC is a rotating disk for which the viewing angles have been well-established using various methods. The disk is elliptical in its disk plane. The disk thickness varies depending on the tracer population, with V/sigma ranging from 2-10 from the oldest to the youngest population. For the SMC, the old stellar population resides in a spheroidal distribution with considerable line-of-sight depth and low V/sigma. Young stars and HI gas reside in a more irregular rotating disk. Mass estimates based on the kinematics indicate that each Cloud is embedded in a dark halo. Proper motion measurements with HST show that both galaxies move significantly more rapidly around the Milky Way than previously believed. This indicates that for a canonical 10^12 solar mass Milky Way the Clouds are only passing by us for the first time. Although a higher Milky Way mass yields a bound orbit, this orbit is still very different from what has been previously assumed in models of the Magellanic Stream. Hence, much of our understanding of the history of the Magellanic System and the formation of the Magellanic Stream may need to be revised. The accuracy of the proper motion data is insufficient to say whether or not the LMC and SMC are bound to each other, but bound orbits do exist within the proper motion error ellipse.

The HST proper motion (PM) measurements of the Clouds have severe implications for their interaction history with the Milky Way (MW) and with each other. The Clouds are likely on their first passage about the MW and the SMC's orbit about the LMC is better described as quasi-periodic rather than circular. Binary L/SMC orbits that satisfy observational constraints on their mutual interaction history (e.g. the formation of the Magellanic Bridge during a collision between the Clouds ~300 Myr ago) can be located within 1 sigma of the mean PMs. However, these binary orbits are not co-located with the Magellanic Stream (MS) when projected on the plane of the sky and the line-of-sight velocity gradient along the LMC's orbit is significantly steeper than that along the MS. These combined results ultimately rule out a purely tidal origin for the MS: tides are ineffective without multiple pericentric passages and can neither decrease the velocity gradient nor explain the offset stream in a polar orbit configuration. Alternatively, ram pressure stripping of an extended gaseous disk may naturally explain the deviation. The offset also suggests that observations of the little-explored region between RA 21h and 23h are crucial for characterizing the full extent of the MS.

Close to its center, the bulgeless galaxy NGC 2139 hosts a star cluster that is younger and less massive than any actual nuclear star cluster (NC) studied so far. We have measured the H-alpha velocity field around the photometric center of this galaxy using the VLT ARGUS integral field unit and GIRAFFE spectrograph in order to constrain different proposed theories of NC formation. We observe that the best-fit kinematic center and the candidate NC appear to be separated by 2.8 arcsec (320 pc). Indeed, the kinematic center also is offset from the galaxy's photometric center and a possible bar or extended region of star formation in which the young cluster resides, implying that this galaxy is not in dynamic equilibrium. The H-alpha flux map also reveals other regions of strong star formation in the possible bar. These observations suggest that a nascent NC is forming away from the kinematic center of NGC 2139 which may come to rest there on a time scale of a few 100 Myr.

We present deep HST ACS/WFC photometry of the dwarf irregular galaxy NGC 1569, one of the closest and strongest nearby starburst galaxies. These data allow us, for the first time, to unequivocally detect the tip of the red giant branch and thereby determine the distance to NGC 1569. We find that this galaxy is 3.36 +/- 0.20 Mpc away, considerably farther away than the typically assumed distance of 2.2 +/- 0.6 Mpc. Previously thought to be an isolated galaxy due to its shorter distance, our new distance firmly establishes NGC 1569 as a member of the IC 342 group of galaxies. The higher density environment may help explain the starburst nature of NGC 1569, since starbursts are often triggered by galaxy interactions. On the other hand, the longer distance implies that NGC 1569 is an even more extreme starburst galaxy than previously believed. Previous estimates of the rate of star formation for stars younger than ~ 1 Gyr become stronger by more than a factor of 2. Stars older than this were not constrained by previous studies. The dynamical masses of NGC 1569's three super star clusters, which are already known as some of the most massive ever discovered, increase by ~53% to 6-7x10^5 solar masses.

In many theoretical scenarios it is expected that intermediate-mass black holes (IMBHs, with masses M ~ 100-10000 solar masses) reside at the centers of some globular clusters. However, observational evidence for their existence is limited. Several previous numerical investigations have focused on the impact of an IMBH on the cluster dynamics or brightness profile. Here we instead present results from a large set of direct N-body simulations including single and binary stars. These show that there is a potentially more detectable IMBH signature, namely on the variation of the average stellar mass between the center and the half-light radius. We find that the existence of an IMBH quenches mass segregation and causes the average mass to exhibit only modest radial variation in collisionally relaxed star clusters. This differs from when there is no IMBH. To measure this observationally requires high resolution imaging at the level of that already available from the Hubble Space Telescope (HST) for the cores of a large sample of galactic globular clusters. With a modest additional investment of HST time to acquire fields around the half-light radius, it will be possible to identify the best candidate clusters to harbor an IMBH. This test can be applied only to globulars with a half-light relaxation time less than or equal to 1 Gyr, which is required to guarantee efficient energy equipartition due to two-body relaxation.

We use the spatial information of our previously published VLT/FORS2 absorption line spectroscopy to measure mean stellar velocity and velocity dispersion profiles of 25 field early-type galaxies at a median redshift z=0.97 (full range 0.6<z<1.2). This provides the first detailed study of early-type galaxy rotation at these redshifts. From surface brightness profiles from HST imaging we calculate two-integral oblate axisymmetric Jeans equation models for the observed kinematics. Fits to the data yield for each galaxy the degree of rotational support and the mass-to-light ratio M/L_Jeans. S0 and Sa galaxies are generally rotationally supported, whereas elliptical galaxies rotate less rapidly or not at all. Down to M(B)=-19.5 (corrected for luminosity evolution), we find no evidence for evolution in the fraction of rotating early-type (E+S0) galaxies between z=1 (63+/-11%) and the present (61+/-5%). We interpret this as evidence for little or no change in the field S0 fraction with redshift. We compare M/L_Jeans with M/L_vir inferred from the virial theorem and globally averaged quantities and assuming homologous evolution. There is good agreement for non-rotating (mostly E) galaxies. However, for rotationally supported galaxies (mostly S0) M/L_Jeans is on average ~40% higher than M/L_vir. We discuss possible explanations and the implications for the evolution of M/L between z=1 and the present and its dependence on mass.

We present spectroscopic evidence for a low-luminosity, low-excitation active galactic nucleus (AGN) in NGC 1042, powered by an intermediate-mass black hole. These findings are significant in that the AGN is coincident with a compact star cluster known to reside in the nucleus, thus providing an example where the two types of central mass concentration coexist. The existence of a central black hole is additionally remarkable in that NGC 1042 lacks a stellar bulge. Objects such as NGC 1042 may have an important role in testing theories for the genesis of massive black holes in galaxy nuclei, and the extent to which they are in symbiosis with the larger stellar host.

(Abridged). We present new high spatial resolution narrowband imaging observations of extraplanar diffuse ionized gas (eDIG) in four late-type, actively star forming edge-on spirals, obtained with ACS on-board HST. Our F658N (H-alpha) observations reveal a multitude of structures on both small and large scales. Whereas all four galaxies have been studied with ground-based telescopes before, here the small scale structure of the extended emission line gas is presented for the first time at a spatial resolution of 0.05", corresponding to 5.0 pc at the mean distance to the target galaxies. The eDIG morphology is very different for all four targets, probably as a result of their different levels of star formation activity. We find that the morphology of the eDIG, in particular the break-up of diffuse emission into filaments in galaxy halos, shows a strong dependence on the level of star formation activity per unit area, and eDIG can be arranged into a morphological sequence. NGC4634 and NGC5775 have the highest SF rate per unit area in our sample and the observed morphology suggests that the break-up of the smooth eDIG layer into individual resolved filaments occurs only above a certain threshold of SF activity per unit area. Combined with ground-based data for samples that span a larger range of galaxy mass our results indicate that the gravitational potential also plays an important role in the eDIG morphology. In low-mass galaxies the gas can be expelled due to shallow gravitational potentials more easily and couple with strong star formation driven outflows on a local scale. This is in contrast to the more massive galaxies, which show smooth eDIG layers, unless they are powered by a superwind, as in the case of nucleated starburst galaxies.

(ABRIDGED) We present a new Schwarzschild orbit-superposition code designed to model discrete datasets composed of velocities of individual kinematic tracers in a dynamical system. This constitutes an extension of previous implementations that can only address continuous data in the form of (the moments of) velocity distributions, thus avoiding potentially important losses of information due to data binning. Furthermore, the code can handle any combination of available velocity components, i.e., only line-of-sight velocities, only proper motions, or a combination of both. It can also handle a combination of discrete and continuous data. The code finds the distribution function (DF, a function of the three integrals of motion E, Lz, and I3) that best reproduces the available kinematic and photometric observations in a given axisymmetric gravitational potential. The fully numerical approach ensures considerable freedom on the form of the DF f(E,Lz,I3). This allows a very general modeling of the orbital structure, thus avoiding restrictive assumptions about the degree of (an)isotropy of the orbits. We describe the implementation of the discrete code and present a series of tests of its performance based on the modeling of simulated datasets generated from a known DF. We find that the discrete Schwarzschild code recovers the original orbital structure, M/L ratios, and inclination of the input datasets to satisfactory accuracy, as quantified by various statistics. The code will be valuable, e.g., for modeling stellar motions in Galactic globular clusters, and those of individual stars, planetary nebulae, or globular clusters in nearby galaxies. This can shed new light on the total mass distributions of these systems, with central black holes and dark matter halos being of particular interest.

We present photometry with the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST) of stars in the Magellanic starburst galaxy NGC 4449. The galaxy has been imaged in the F435W (B), F555W (V) and F814W (I) broad-band filters, and in the F658N (Halpha) narrow-band filter. Our photometry includes ~300,000 objects in the (B, V) color-magnitude diagram (CMD) down to V < 28, and ~400,000 objects in the (V, I) CMD, down to I < 27 . A subsample of ~200,000 stars has been photometrized in all the three bands simultaneously. The features observed in the CMDs imply a variety of stellar ages up to at least 1 Gyr, and possibly as old as a Hubble time. The spatial variation of the CMD morphology and of the red giant branch colors point toward the presence of an age gradient: young and intermediate-age stars tend to be concentrated toward the galactic center, while old stars are present everywhere. The spatial variation in the average luminosity of carbon stars suggests that there is not a strong metallicity gradient (< 0.2 dex). Also, we detect an interesting resolved star cluster on the West side of the galaxy, surrounded by a symmetric tidal or spiral feature consisting of young stars. The positions of the stars in NGC 4449 younger than 10 Myr are strongly correlated with the Halpha emission. We derive the distance of NGC 4449 from the tip of the red giant branch to be D=3.82 \pm 0.27 Mpc. This result is in agreement with the distance that we derive from the luminosity of the carbon stars.

We present several different statistical methods to determine the transverse velocity vector of M31. The underlying assumptions are that the M31 satellites on average follow the motion of M31 through space, and that the galaxies in the outer parts of the Local Group on average follow the motion of the Local Group barycenter through space. We apply the methods to the line-of-sight velocities of 17 M31 satellites, to the proper motions of the 2 satellites M33 and IC 10, and to the line-of-sight velocities of 5 galaxies near the Local Group turn around radius, respectively. This yields 4 independent but mutually consistent determinations of the heliocentric M31 transverse velocities in the West and North directions, with weighted averages <v_W> = -78+/-41 km/s and <v_N> = -38+/-34 km/s. The Galactocentric tangential velocity of M31 is 42 km/s, with 1-sigma confidence interval V_tan <= 56 km/s. The implied M31-Milky Way orbit is bound if the total Local Group mass M exceeds 1.72^{+0.26}_{-0.25}x10^{12} solar masses. If indeed bound, then the timing argument combined with the known age of the Universe implies that M = 5.58^{+0.85}_{-0.72}x10^{12} solar masses. This is on the high end of the allowed mass range suggested by cosmologically motivated models for the individual structure and dynamics of M31 and the Milky Way, respectively. It is therefore possible that the timing mass is an overestimate of the true mass, especially if one takes into account recent results from the Millennium Simulation that show that there is also a theoretical uncertainty of 41 percent (Gaussian dispersion) in timing mass estimates. The M31 transverse velocity implies that M33 is in a tightly bound orbit around M31. This may have led to some tidal deformation of M33. It will be worthwhile to search for observational evidence of this.

We discuss the near-infrared properties of the nuclei in the 11 merging galaxies of the Toomre sequence, based on high spatial resolution J, H, and K imaging data using NICMOS onboard the Hubble Space Telescope (HST). The observations are less affected by dust extinction than our previous HST/WFPC2 observations and offer higher spatial resolution than existing ground-based near-IR data. We see a marginal trend for the nuclei to become bluer with advancing merger stage, which we attribute to a dispersal of dust at late times in the merging process. Our data also indicate a statistically significant trend for the nuclei in the sequence to become more luminous, within an aperture of fixed physical size and after correcting for dust extinction, with advancing merger stage. We derive K-band surface brightness profiles for those nuclei for which the morphology allows a meaningful isophotal analysis, and fit the profiles with a ``Nuker law'' for comparison with other samples of galaxies observed with HST. The majority of the nuclei have steep profiles that can be characterized as power-law type. In general, the Toomre sequence galaxies tend to have steeper profiles and higher central luminosity surface densities than E/S0's. We derive V-K color profiles for the nuclei to further address this possibility, but find that the large amounts of dust extinction complicate their interpretation. Overall, our results are consistent with the generic predictions of N-body simulations of spiral galaxy mergers. If left to evolve and fade for several Gyrs, it is possible that the properties of the Toomre sequence nuclei would resemble theproperties of the nuclei of normal E/S0 galaxies. Our results therefore support the view that mergers of spiral galaxies can lead to the formation of early-type galaxies (Abridged).

We present new V and I-band HST/ACS photometry of I Zw 18, the most metal-poor blue compact dwarf (BCD) galaxy in the nearby universe. It has been argued in the past that I Zw 18 is a very young system that started forming stars only < 500 Myr ago, but other work has hinted that older (> 1 Gyr) red giant branch (RGB) stars may also exist. Our new data, once combined with archival HST/ACS data, provide a deep and uncontaminated optical color-magnitude diagram (CMD) that now strongly indicates an RGB. The RGB tip (TRGB) magnitude yields a distance modulus (m-M)_0 = 31.30 +/- 0.17, i.e., D = 18.2 +/- 1.5 Mpc. The time-series nature of our observations allows us to also detect and characterize for the first time three classical Cepheids in I~Zw~18. The time-averaged Cepheid <V> and <I> magnitudes are compared to the VI reddening-free Wesenheit relation predicted from new non-linear pulsation models specifically calculated at the metallicity of I Zw 18. For the one bona-fide classical Cepheid with a period of 8.63 days this implies a distance modulus (m-M)_0 = 31.42 +/- 0.26. The other two Cepheids have unusually long periods (125.0 and 129.8 d) but are consistent with this distance. The coherent picture that emerges is that I Zw 18 is older and farther away than previously believed. This rules out the possibility that I Zw 18 is a truly primordial galaxy formed recently (z < 0.1) in the local universe.

Recent proper motion measurements of the Large and Small Magellanic Clouds (LMC and SMC, respectively) by Kallivayalil et al (2006a,b) suggest that the 3D velocities of the Clouds are substantially higher (~100 km/s) than previously estimated and now approach the escape velocity of the Milky Way (MW). Previous studies have also assumed that the Milky Way can be adequately modeled as an isothermal sphere to large distances. Here we re-examine the orbital history of the Clouds using the new velocities and a LCDM-motivated MW model with virial mass Mvir = 1e12 Msun (e.g. Klypin et al 2002). We conclude that the Clouds are either currently on their first passage about the MW or, if the MW can be accurately modeled by an isothermal sphere to distances >200 kpc (i.e., Mvir > 2e12 Msun), that their orbital period and apogalacticon distance must be a factor of two larger than previously estimated, increasing to 3 Gyr and 200 kpc, respectively. A first passage scenario is consistent with the fact that the LMC and SMC appear to be outliers when compared to other satellite galaxies of the MW: they are irregular in appearance and are moving faster. We discuss the implications of this orbital analysis for our understanding of the star formation history, the nature of the warp in the MW disk and the origin of the Magellanic Stream (MS), a band of HI gas trailing the LMC and SMC that extends roughly 100 degrees across the sky. Specifically, as a consequence of the new orbital history of the Clouds, the origin of the MS may not be explainable by current tidal and ram pressure stripping models.

We present new Space Telescope Imaging Spectrograph (STIS) observations of three spiral galaxies, NGC 4303, NGC 3310 and NGC 4258. The bright optical emission lines H$\alpha$ $\lambda$ $6564 \AA$, [NII] $\lambda$$\lambda$ $6549,6585 \AA$ and [SII] $\lambda$$\lambda$ $ 6718,6732 \AA$ were used to study the kinematics of the ionized gas in the nuclear region of each galaxy with a $\sim 0.07\arcsec$ spatial resolution. In NGC 3310, the observed gas kinematics is well matched by a circularly rotating disk model but we are only able to set an upper limit to the BH mass which, taking into account the allowed disk inclinations, varies in the range $5.0 \times 10^{6} - 4.2 \times 10^{7} M_{\odot}$ at the 95% confidence level. In NGC 4303 the kinematical data require the presence of a BH with mass $M_{BH}=(5.0)^{+0.87}_{-2.26}\times 10^{6}M_{\odot}$ (for a disk inclination $i=70$ deg).In NGC 4258, the observed kinematics require the presence of a black hole with $M_{BH}= (7.9)^{+6.2}_{-3.5} \times 10^{7}M_{\odot}$ ($i=60$ deg). This result is in good agreement with the published value $(3.9 \pm 0.1) \times 10^{7} M_{\odot}$, derived from $H_{2}O$-maser observations. Our attempt at measuring BH masses in these 3 late type Sbc spiral galaxies has shown that these measurements are very challenging and at the limit of the highest spatial resolution currently available. Nonetheless our estimates are in good agreement with the scaling relations between black holes and their host spheroids suggesting that (i) they are reliable and (ii) black holes in spiral galaxies follows the same scaling relations as those in more massive early-type galaxies. A crucial test for the gas kinematical method, the correct recovery of the known BH mass in NGC 4258, has been successful. [abridged]

HST observations have revealed that compact sources exist at the centers of many, maybe even most, galaxies across the Hubble sequence. These sources are called "nuclei" or also "nuclear star clusters" (NCs), given that their structural properties and position in the fundamental plane are similar to those of globular clusters. Interest in NCs increased recently due to the independent and contemporaneous finding of three groups (Rossa et al. for spiral galaxies; Wehner & Harris for dE galaxies; and Cote et al. for elliptical galaxies) that NC masses obey similar scaling relationships with host galaxy properties as do supermassive black holes. Here we summarize the results of our group on NCs in spiral galaxies. We discuss the implications for our understanding of the formation and evolution of NCs and their possible connection to supermassive black holes.

We present results from new deep HST/ACS photometry of I Zw 18, the most metal-poor blue compact dwarf galaxy in the nearby universe. It has been previously argued that this is a very young system that started forming stars only < 500 Myr ago, but other work has hinted that older (> 1 Gyr) red giant branch (RGB) stars may exist in this galaxy. Our deeper data indeed reveal evidence for an RGB. Underlying old (> 1 Gyr) populations are therefore present in even the most metal-poor systems, implying that star formation started at z > 0.1. The RGB tip (TRGB) magnitude and the properties of Cepheid variables identified from our program indicate that I Zw 18 is farther away (D = 19.0 +/- 1.8 Mpc) than previously believed.

(Abridged) Spiral galaxies often have extended outflows that permeate beyond the region of the disk. Such outflows have been seen both in starburst galaxies, actively star forming galaxies and galaxies with an AGN. In the latter galaxies it is unknown whether the large-scale outflows are driven by star formation activity or purely by the active nucleus. The aim of our investigation is to study the frequency of extended minor-axis outflows in edge-on Seyfert galaxies to investigate the role of the AGN, the circumnuclear environment and star formation activity within the disk regions, and their importance for IGM enrichment on large scales. Narrowband imaging in two different ionizational stages (H-alpha and [OIII]) was performed to attempt a discrimination between processes associated with the active nucleus and those connected to star forming activity within the disk. The H-alpha morphology of the Seyfert galaxies is usually complex, but only in three out of 14 galaxies did we find evidence for minor axis disk outflows. At the sensitivity of our observations [OIII] emission is generally detected only in the nuclear region. Overall, our results show that extraplanar emission of similar brightness and extent as in the previously known cases of NGC3079 and NGC4388 is not common in Seyfert galaxies of otherwise similar properties. Comparison with our previous results shows that for nearby edge-on spiral galaxies star formation may be a more powerful mechanism for producing DIG than AGN activity. While in general AGN activity undoubtedly plays some role in driving minor-axis outflows, this probably requires higher AGN luminosities than are encountered in our small distance-limited sample.

We present new Keck spectroscopy of early-type galaxies in three galaxy clusters at z~0.5. We focus on the fundamental plane (FP) relation, and combine the kinematics with structural parameters determined from HST images. The galaxies obey clear FP relations, which are offset from the FP of the nearby Coma cluster due to passive evolution of the stellar populations. The z~0.5 data are combined with published data for 11 additional clusters at 0.18<z<1.28, to determine the evolution of the mean M/L(B) ratio of cluster galaxies with masses M>10^11 M_sun, as implied by the FP. We find dlog(M/L(B))/dz = -0.555+-0.042, stronger evolution than was previously inferred from smaller samples. The observed evolution depends on the luminosity-weighted mean age of the stars in the galaxies, the initial mass function (IMF), selection effects due to progenitor bias, and other parameters. Assuming a normal IMF but allowing for various other sources of uncertainty we find z* = 2.01+-0.20 for the luminosity-weighted mean star formation epoch. The main uncertainty is the slope of the IMF in the range 1-2 Solar masses: we find z* = 4.0 for a top-heavy IMF with slope x=0. The M/L(B) ratios of the cluster galaxies are compared to those of recently published samples of field early-type galaxies at 0.32<z<1.14. Assuming that progenitor bias and the IMF do not depend on environment we find that the present-day age of stars in massive field galaxies is 4.1 +- 2.0 % (~0.4 Gyr) less than that of stars in massive cluster galaxies, consistent with most, but not all, previous studies of local and distant early-type galaxies. This relatively small age difference is surprising in the context of expectations from ``standard'' hierarchical galaxy formation models. [ABRIDGED]

We study M/L evolution of early-type galaxies using dynamical modeling of resolved internal kinematics. This makes fewer assumptions than Fundamental Plane (FP) studies and provides a powerful new approach for studying galaxy evolution. We focus on the sample of 25 galaxies in clusters at z=0.5 modeled in Paper I. For comparison we compile and homogenize M/L literature data for 60 nearby galaxies that were modeled in comparable detail. The nearby sample obeys log(M/L)_B = Z + S log(sigma_eff/[200 km/s]), with Z = 0.896 +/- 0.010, S = 0.992 +/- 0.054, and sigma_eff the effective velocity dispersion. The z=0.5 sample follows a similar relation but with lower zeropoint. The implied M/L evolution is Delta log(M/L) / Delta z = -0.457 +/- 0.046(random) +/- 0.078(systematic), consistent with passive evolution following high-redshift formation. This agrees with the FP results for this sample by van Dokkum & van der Marel. This confirms that FP evolution tracks M/L evolution, which is an important verification of the assumptions that underly FP studies. However, while we find more FP evolution for galaxies of low sigma_eff (or low mass), the dynamical M/L evolution instead shows little trend with sigma_eff. We argue that this difference can be plausibly attributed to a combination of two effects: (a) evolution in structural galaxy properties other than M/L; and (b) the neglect of rotational support in studies of FP evolution. The results leave the question open whether the low-mass galaxies in the sample have younger population ages than the high-mass galaxies. This highlights the general importance in the study of population ages for complementing dynamical measurements with broad-band colors or spectroscopic population diagnostics.

We present spatially resolved stellar rotation velocity and velocity dispersion profiles form Keck/LRIS absorption-line spectra for 25 galaxies, mostly visually classified ellipticals, in three clusters at z=0.5. We interpret the kinematical data and HST photometry using oblate axisymmetric two-integral f(E,Lz) dynamical models based on the Jeans equations. This yields good fits, provided that the seeing and observational characteristics are carefully modeled. The fits yield for each galaxy the dynamical M/L and a measure of the galaxy rotation rate. Paper II addresses the implied M/L evolution. Here we study the rotation-rate evolution by comparison to a sample of local elliptical galaxies of similar present-day luminosity. The brightest galaxies in the sample all rotate too slowly to account for their flattening, as is also observed at z=0. But the average rotation rate is higher at z=0.5 than locally. This may be due to a higher fraction of misclassified S0 galaxies (although this effect is insufficient to explain the observed strong evolution of the cluster S0 fraction with redshift). Alternatively, dry mergers between early-type galaxies may have decreased the average rotation rate over time. It is unclear whether such mergers are numerous enough in clusters to explain the observed trend quantitatively. Disk-disk mergers may affect the comparison through the so-called progenitor bias, but this cannot explain the direction of the observed rotation-rate evolution. Additional samples are needed to constrain possible environmental dependencies and cosmic variance in galaxy rotation rates. Either way, studies of the internal stellar dynamics of distant galaxies provide a valuable new approach for exploring galaxy evolution.

We present a measurement of the systemic proper motion of the Small Magellanic Cloud (SMC) made using the Advanced Camera for Surveys (ACS) on the \textit{Hubble Space Telescope} (\textit{HST}). We tracked the SMC's motion relative to 4 background QSOs over a baseline of approximately 2 years. The measured proper motion is : $\mu_W = -1.16 \pm 0.18 \masyr, \mu_N = -1.17 \pm 0.18 \masyr$. This is the best measurement yet of the SMC's proper motion. We combine this new result with our prior estimate of the proper motion of the Large Magellanic Cloud (LMC) from the same observing program to investigate the orbital evolution of both Clouds over the past 9 Gyr. The current relative velocity between the Clouds is $105 \pm 42 \kms$. Our investigations of the past orbital motions of the Clouds in a simple model for the dark halo of the Milky Way imply that the Clouds could be unbound from each other. However, our data are also consistent with orbits in which the Clouds have been bound to each other for approximately a Hubble time. Smaller proper motion errors and better understanding of the LMC and SMC masses would be required to constrain their past orbital history and their bound vs. unbound nature unambiguously. The new proper motion measurements should be sufficient to allow the construction of improved models for the origin and properties of the Magellanic Stream. In turn, this will provide new constraints on the properties of the Milky Way dark halo.

We have created a website, called "Black Holes: Gravity's Relentless Pull", which explains the physics and astronomy of black holes for a general audience. The site emphasizes user participation and is rich in animations and astronomical imagery. It won the top prize of the 2005 Pirelli INTERNETional Awards competition for the best communication of science and technology using the internet. This article provides a brief overview of the site. The site starts with an opening animation that introduces the basic concept of a black hole. The user is then invited to embark on a journey from a backyard view of the night sky to a personal encounter with a singularity. This journey proceeds through three modules, which allow the user to: find black holes in the night sky; travel to a black hole in an animated starship; and explore a black hole from up close. There are also five "experiments" that allow the user to: create a black hole; orbit around a black hole; weigh a black hole; drop a clock into a black hole; or fall into a black hole. The modules and experiments offer goal-based scenarios tailored for novices and children. The site also contains an encyclopedia of frequently asked questions and a detailed glossary that are targeted more at experts and adults. The overall result is a website where scientific knowledge, learning theory, and fun converge. Despite its focus on black holes, the site also teaches many other concepts of physics, astronomy and scientific thought. The site aims to instill an appreciation for learning and an interest in science, especially in the younger users. It can be used as an aid in teaching introductory astronomy at the undergraduate level.

We present a systematic analysis of archival HST WFPC2 ``Association'' data sets that correlate with the Chandra positions of a set of 44 ultra-luminous X-ray sources (ULXs) of nearby galaxies. We have improved the Chandra-HST relative astrometry whenever possible. Disparate numbers of potential ULX counterparts are found, and in some cases none are found. The lack of or low number of counterparts in some cases may be due to insufficient depth in the WFPC2 images. Particularly in late-type galaxies, the HST image in the ULX region was often complex or crowded. We therefore address various scenarios for the nature of the ULX since it is not known which, if any, of the sources found are true counterparts. The optical luminosities of the sources are typically in the range 10^4-6 L_sun. In several cases color information is available, with the colors roughly tending to be more red in early-type galaxies. This suggests that, in general, the (potential) counterparts found in early-type galaxies are likely to be older stellar populations, and are probably globular clusters. Several early-type galaxy counterparts have blue colors, which may be due to younger stellar populations in the host galaxies, however these could also be background sources. In spiral galaxies the sources may also be due to localized structure in the disks rather than bound stellar systems. Alternatively some of the counterparts in late-type galaxies may be isolated supergiant stars. The observed X-ray/optical flux ratio is diluted by the optical emission of the cluster in cases where the system is an X-ray binary in a cluster, particularly in the case of a low-mass X-ray binaries in old cluster. (abridged)

[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.

(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.

(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.

(abridged) We analyze and model HST /STIS observations of a sample of 27 galaxies; 16 Fanaroff & Riley Type I radio galaxies and 11 (more) normal early-type galaxies. We focus here on what can be learned from the nuclear velocity dispersion (line width) of the gas as a complement to the many studies dealing with gas rotation velocities. We find that the dispersion in a STIS aperture of ~0.1''-0.2'' generally exceeds the large-scale stellar velocity dispersion of the galaxy. This is qualitatively consistent with the presence of central BHs, but raises the question whether the excess gas dispersion is of gravitational or non-gravitational origin and whether the implied BH masses are consistent with our current understanding of BH demography(as predicted by the M-sigma relation between BH mass and stellar velocity dispersion). To address this we construct dynamical models for the gas, both thin disk models and models with more general axis ratios and velocity anisotropies. For the normal galaxies the nuclear gas dispersions are adequately reproduced assuming disks around BHs with masses that follow the M-sigma relation. In contrast, the gas dispersions observed for the radio galaxies generally exceed those predicted by any of the models. We attribute this to the presence of non-gravitational motions in the gas that are similar to or larger than the gravitational motions. The non- gravitational motions are presumably driven by the active galactic nucleus (AGN), but we do not find a relation between the radiative output of the AGN and the non-gravitational dispersion. It is not possible to uniquely determine the BH mass for each galaxy from its nuclear gas dispersion. However, for the sample as a whole the observed dispersions do not provide evidence for significant deviations from the M-sigma relation.

We present results from an HST/ACS imaging study of the metal-poor blue compact dwarf galaxy SBS 1415+437. It has been argued previously that this is a very young galaxy that started to form stars only less than 100 Myr ago. However, we find that the optical color-magnitude diagram prominently reveals asymptotic giant branch and red giant branch (RGB) stars. The brightness of the RGB tip yields a distance of 13.6 Mpc. The color of the RGB implies that its stars must be older than approximately 1.3 Gyr, with the exact age depending on the assumed metallicity and dust extinction. The number of RGB stars implies that most of the stellar mass resides in this evolved population. In view of these and other HST results for metal-poor galaxies it seems that the local Universe simply may not contain any galaxies that are currently undergoing their first burst of star formation.

We present a measurement of the systemic proper motion of the Large Magellanic Cloud (LMC) from astrometry with the High Resolution Camera (HRC) of the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope (HST). We observed LMC fields centered on 21 background QSOs that were discovered from their optical variability in the MACHO database. The QSOs are distributed homogeneously behind the central few degrees of the LMC. With 2 epochs of HRC data and a ~2 year baseline we determine the proper motion of the LMC to better than 5% accuracy: mu_W = -2.03 +/- 0.08 mas/yr; mu_N = 0.44 +/- 0.05 mas/yr. This is the most accurate proper motion measurement for any Milky Way satellite thus far. When combined with HI data from the Magellanic Stream this should provide new constraints on both the mass distribution of the Galactic Halo and models of the Stream.

We report stellar velocity dispersion measurements for a sample of 28 AGN host galaxies including our previous work. Using the mass-dispersion ($M_{\bullet}-\sigma$) and the fundamental plane relations, we estimate the black hole mass for a sample of 66 BL Lac objects and investigate the role of black hole mass in the energetics of BL Lac objects. The black hole mass range for different BL Lac spectral types is similar, $10^{7} < M_{\bullet} < 4 \times 10^{9}$. Neither X-ray nor radio luminosity correlates with black hole mass. Low-frequency-peaked BL Lac objects have higher Eddington ratios on average, because of either more beaming or higher intrinsic power. For the black hole mass range $3 \times 10^{7} < M_{\bullet} < 10^{9}$, the radio luminosity of BL Lac objects and flat-spectrum radio quasars spans over 4 orders of magnitude, with BL Lac objects being low-power AGNs. We also investigate the evolution of host galaxies for 39 AGNs out to $z \approx 0.5$ with measuredstellar velocity dispersions. Comparing the mass-to-light ratio evolution in the observed frame with population synthesis models, we find that single burst star formation models with $z_{form} = 1.4^{+0.9}_{-0.2} $ are consistent with the observations. From our $z_{form}=1.4$ model, we estimated the intrinsic mass-to-light ratio evolution in the Cousins $R$ band, $\Delta log (M/L)/ \Delta z = -0.502 \pm 0.08$, consistent with that of normal early type galaxies.

We present evidence for a counter-rotating core in the low-luminosity (M_B = -18.2) elliptical galaxy NGC 770 based on internal stellar kinematic data. This counter-rotating core is unusual as NGC 770 is not the primary galaxy in the region and it lies in an environment with evidence of on-going tidal interactions. We discovered the counter-rotating core via single-slit Keck/ESI echelle spectroscopy; subsequent integral field spectroscopy was obtained with the Gemini/GMOS IFU. The counter-rotating region has a peak rotation velocity of 21 km/s as compared to the main galaxy's rotation speed of greater than 45 km/s in the opposite direction. The counter-rotating region extends to a radius of 4'' (0.6 kpc), slightly smaller than the half-light radius of the galaxy which is 5.3'' (0.8 kpc) and is confined to a disk whose scale height is less than 0.8'' (0.1 kpc). We compute an age and metallicity of the inner counter-rotating region of 3 +/- 0.5 Gyr and [Fe/H] = 0.2 +/- 0.2 dex, based on Lick absorption-line indices. The lack of other large galaxies in this region limits possible scenarios for the formation of the counter-rotating core. We discuss several scenarios and favor one in which NGC 770 accreted a small gas-rich dwarf galaxy during a very minor merging event. If this scenario is correct, it represents one of the few known examples of merging between two dwarf-sized galaxies.

We present Space Telescope Imaging Spectrograph emission-line spectra of the central regions of the spiral galaxies NGC 1300 and NGC 2748. From the derived kinematics of the nuclear gas we have found evidence for central supermassive black holes in both galaxies. The estimated mass of the black hole in NGC 1300 is 6.6 (+6.3, -3.2) x 10^7 solar masses and in NGC 2748 is 4.4 (+3.5, -3.6) x 10^7 solar masses (both at the 95% confidence level). These two black hole mass estimates contribute to the poorly sampled low-mass end of the nuclear black hole mass spectrum.

In the last decade star clusters have been found in the centers of spiral galaxies across all Hubble types. We here present a spectroscopic study of the exceptionally bright (10^6 - 10^8 Lsun) but compact (Re ~ 5 pc) nuclear star clusters in very late type spirals with UVES at the VLT. We find the velocity dispersions of the nine clusters in our sample to range from 13 to 34 km/s. Using photometric data from the HST/WFPC2 and spherically symmetric dynamical models we determine masses between 8*10^5 and 6*10^7 Msun. The mass to light ratios range from 0.2 to 1.5 in the I band. This indicates a young mean age for most clusters, in agreement with previous studies. Given their high masses and small sizes we find that nuclear clusters are among the objects with the highest mean surface density known (up to 10^5 Msun pc^-2). From their dynamical properties we infer that, rather than small bulges, the closest structural kin of nuclear clusters appear to be massive compact star clusters. This includes such different objects as globular clusters, "super star clusters", ultra compact dwarf galaxies and the nuclei of dwarf elliptical galaxies. It is a challenge to explain why, despite the wildly different current environments, all different types of massive star clusters share very similar and structural properties. A possible explanation links UCDs and massive globular clusters to nuclear star clusters through stripping of nucleated dwarf galaxies in a merger event. The extreme properties of this type of clusters would then be a consequence of their location in the centers of their respective host galaxies.

We measured the stellar velocity dispersions of 15 active galactic nucleus (AGN) host galaxies at redshifts as high as $\sim 0.34$. Combining these with published velocity dispersion measurements from the literature, we study the Fundamental Plane of AGN host galaxies and its evolution. BL Lac hosts and radio galaxies seem to lie on the same Fundamental Plane as normal early-type galaxies. The evolution of the mass-to-light ratio of AGN host galaxies shows a similar trend to that observed in normal early-type galaxies, consistent with single-burst passive evolution models with formation redshifts $z \gtrsim 1$. The lack of a significant difference between normal and AGN host galaxies in the Fundamental plane supports the "Grand Unification" picture wherein AGNs are a transient phase in the evolution of normal galaxies. The black hole masses of BL Lac objects and radio galaxies, derived using the mass -- dispersion relation, are similar. The black hole mass is independent of BL Lac type. The local black hole mass -- host galaxy luminosity relation of our sub-sample at $z < 0.1$ is similar to that of local normal and radio galaxies, but is less well defined at higher redshift due to the luminosity evolution of the host galaxies.

We present surface photometry for the central regions of a sample of 48 spiral galaxies (mostly unbarred and barred of types Sbc or Sc) observed with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. Surface brightness profiles were derived and modeled with a Nuker law. We also analyzed archival Wide Field Planetary Camera~2 images with a larger field of view, available for 18 galaxies in our sample. We modeled the extracted bulge surface brightness profiles with an exponential, a de Vaucouleurs or a Sersic profile. In agreement with previous studies, we find that bulges of Sbc galaxies fall into two categories: bulges well described by an exponential profile and those well described by an de Vaucouleurs profile. Only one galaxy requires the use of a more general Sersic profile to properly describe the bulge. Nuclear photometrically distinct components are found in ~55% of the galaxies. For those that we classify as star clusters based on their resolved extent we find absolute magnitudes that are brighter on average than those previously identified in spiral galaxies. This might be due to a bias in our sample toward star forming galaxies, combined with a trend for star forming galaxies to host brighter central clusters.

I review our understanding of the structure and kinematics of the Large Magellanic Cloud (LMC), with a particular focus on recent results. This is an important topic, given the status of the LMC as a benchmark for studies of microlensing, tidal interactions, stellar populations, and the extragalactic distance scale. I address the observed morphology and kinematics of the LMC; the angles under which we view the LMC disk; its in-plane and vertical structure; the LMC self-lensing contribution to the total microlensing optical depth; the LMC orbit around the Milky Way; and the origin and interpretation of the Magellanic Stream. Our understanding of these topics is evolving rapidly, in particular due to the many large photometric and kinematic datasets that have become available in the last few years. It has now been established that: the LMC is considerably elongated in its disk plane; the LMC disk is thicker than previously believed; the LMC disk may have warps and twists; the LMC may have a pressure-supported halo; the inner regions of the LMC show unexpected complexities in their vertical structure; and precession and nutation of the LMC disk plane contribute measurably to the observed line-of-sight velocity field. However, many open questions remain and more work is needed before we can expect to converge on a fully coherent structural, dynamical and evolutionary picture that explains all observed features of the LMC.

We summarize some recent results from our observational campaign to study the central regions of spiral galaxies of late Hubble type (Scd - Sm). These disk-dominated, bulgeless galaxies have an apparently uneventful merger history. The evolution of their nuclei thus holds important constraints on the mechanisms that are responsible for bulge formation and nuclear activity in spiral galaxies. We discuss the structural properties, stellar populations, and dynamical masses of the compact, luminous star cluster that is found in the nuclei of most late-type spiral galaxies. Although preliminary, our results strongly indicate that many galaxies of our sample experience repeated periods of nuclear star formation. While the exact mechanism that leads to the required high gas densities in the galaxy nucleus remains unclear, results from our recent CO survey of late-type spirals demonstrate that in most cases, the central kpc contains enough molecular gas to support repetitive nuclear starbursts.

We present results from an HST study of the morphology and kinematics of NGC 6240. This merging galaxy with a double nucleus is one of the nearest and best-studied ultraluminous infrared galaxies. HST resolves both nuclei into seperate components. The distance between the northern and southern optical/near-infrared components is greater than that observed in radio and X-ray studies, arguing that even in K-band we may not be seeing all the way through the dust to the true nuclei. The ionized gas does not display rotation around either of the nuclei, or equilibrium motion in general. There is a strong velocity gradient between the nuclei, similar to what is seen in CO data. There is no such gradient in our stellar kinematics. The velocity dispersion of the gas is larger than expected for a cold disk. We also map and model the emission-line velocity field at an off-nuclear position where a steep velocity gradient was previously detected in ground-based data. Overall, the data indicate that line-of-sight projection effects, dust absorption, non-equilibrium merger dynamics, and the possible influence of a wind may be playing an important role in the observed kinematics. Chandra observations of hard X-rays have shown that both of the nuclei contain an Active Galactic Nucleus (AGN). The HST data show no clear sign of the two AGNs: neither continuum nor narrow-band imaging shows evidence for unresolved components in the nuclei, and there are no increased emission line widths or rapid rotation near the nuclei. This underscores the importance of X-ray data for identifying AGNs in highly dust-enshrouded environments.

We investigate the structural properties of nuclear star clusters in late-type spiral galaxies. More specifically, we fit analytical models to HST images of 39 nuclear clusters in order to determine their effective radii after correction for the instrumental point spread function. We use the results of this analysis to compare the luminosities and sizes of nuclear star clusters to those of other ellipsoidal stellar systems, in particular the Milky Way globular clusters. Our nuclear clusters have a median effective radius of r_e = 3.5 pc, with 50% of the sample falling between 2.4 pc < r_e < 5.0 pc. This narrow size distribution is statistically indistinguishable from that of Galactic globular clusters, even though the nuclear clusters are on average 4 magnitudes brighter than the old globulars. We discuss some possible interpretations of this result. From a comparison of nuclear cluster luminosities with various properties of their host galaxies, we confirm that more luminous galaxies harbor more luminous nuclear clusters. It remains unclear whether this correlation mainly reflects the influence of galaxy size, mass, and/or star formation rate. Since the brighter galaxies in our sample typically have stellar disks with a higher central surface brightness, nuclear cluster luminosity also correlates with this property of their hosts. On the other hand, we find no evidence for a correlation between the presence of a nuclear star cluster and the presence of a large-scale stellar bar.

We report on the properties of nuclear regions in the Toomre Sequence of merging galaxies, based on imaging data gathered with the Hubble Space Telescope WFPC2 camera. We have imaged the 11 systems in the proposed evolutionary merger sequence in the F555W and F814W broad-band filters, and in H-alpha+[NII] narrow-band filters. The broad-band morphology of the nuclear regions varies from non-nucleated starburst clumps through dust-covered nuclei to a nucleated morphology. There is no unambiguous trend in the morphology with merger stage. The emission-line morphology is extended beyond the nucleus in most cases, but centrally concentrated (within 1 kpc) emission-line gas can be seen in the four latest-stage merger systems. We have quantified the intrinsic luminosity densities and colors within the inner 100 pc and 1 kpc of each identified nucleus. We find little evidence for a clear trend in nuclear properties along the merger sequence, other than a suggestive rise in the nuclear luminosity density in the most evolved members of the sequence. The lack of clear trends in nuclear properties is likely due both to the effects of obscuration and geometry, as well as the physical variety of galaxies included in the Toomre Sequence.

Current understanding of the secular evolution of galactic disks suggests that this process is dominated by two or more heating mechanisms, which increase the random motions of stars in the disk. In particular, the gravitational influence of giant molecular clouds and irregularities in the spiral potential have been proposed to explain the observed velocity dispersions in the solar neighborhood. Each of these mechanisms acts on different components of the stellar velocities, which affects the ratio of the vertical and radial components of the stellar velocity dispersion since the relative strengths of giant molecular clouds and spiral irregularities vary with Hubble type. A study of this ratio as function of Hubble type has the potential to provide strong constraints on disk heating mechanisms. We present major and minor axis stellar kinematics for four spiral galaxies of Hubble type from Sa to Sbc, and use the data to infer the ratios sigma_z/sigma_R in the galaxy disks. The results combined with two galaxies studied previously and with Milky Way data show that the ratio is generally in the range 0.5 - 0.8. There is a marginally significant trend of decreasing ratio with advancing Hubble type, consistent with the predictions of disk heating theories. However, the errors on individual measurements are large, and the absence of any trend is consistent with the data at the 1-sigma level. As a byproduct of our study, we find that three of the four galaxies in our sample have a central drop in their stellar line-of-sight velocity dispersion, a phenomenon that is increasingly observed in spiral galaxies. [ABRIDGED]

We present spatially-resolved internal kinematics and stellar chemical abundances for a sample of dwarf elliptical (dE) galaxies in the Virgo Cluster observed with Keck/ESI. We find that 4 out of 17 dEs have major axis rotation velocities consistent with rotational flattening, while the remaining dEs have no detectable major axis rotation. Despite this difference in internal kinematics, rotating and non-rotating dEs are remarkably similar in terms of their position in the Fundamental Plane, morphological structure, stellar populations, and local environment. We present evidence for faint underlying disks and/or weak substructure in a fraction of both rotating and non-rotating dEs, but a comparable number of counter-examples exist for both types which show no evidence of such structure. Absorption-line strengths were determined based on the Lick/IDS system (Hbeta, Mgb, Fe5270, Fe5335) for the central region of each galaxy. We find no difference in the line-strength indices, and hence stellar populations, between rotating and non-rotating dE galaxies. The best-fitting mean age and metallicity for our 17 dE sample are 5 Gyr and Fe/H = -0.3 dex, respectively, with rms spreads of 3 Gyr and 0.1 dex. The majority of dEs are consistent with solar alpha/Fe abundance ratios. By contrast, the stellar populations of classical elliptical galaxies are, on average, older, more metal rich, and alpha-enhanced relative to our dE sample. The local environments of both dEs types appear to be diverse in terms of their proximity to larger galaxies in real or velocity space within the Virgo Cluster. Thus, rotating and non-rotating dEs are remarkably similar in terms of their structure, stellar content, and local environments, presenting a significant challenge to theoretical models of their formation. (abridged)

We present the results of the analysis of a set of medium resolution spectra, obtained by the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope, of the emission line gas present in the nuclei of a complete sample of 21 nearby, early-type galaxies with radio jets (the UGC FR-I Sample). For each galaxy nucleus we present spectroscopic data in the region of H-alpha and the dervived kinematics. We find that in 67% of the nuclei the gas appears to be rotating and, with one exception, the cases where rotation is not seen are either face on or have complex central morphologies. We find that in 62% of the nuclei the fit to the central spectrum is improved by the inclusion of a broad component. The broad components have a mean velocity dispersion of 1349 +/- 345 km\s and are redshifted from the narrow line components (assuming an origin in H-alpha) by 486 +/- 443 km\s.

The majority of nearby early-type galaxies contains detectable amounts of emission-line gas at their centers. The emission-line ratios and gas kinematics potentially form a valuable diagnostic of the nuclear activity and gravitational potential well. The observed central gas velocity dispersion often exceeds the stellar velocity dispersion. This could be due to either the gravitational potential of a black hole or turbulent shocks in the gas. Here we try to discriminate between these two scenarios.

This paper reviews the subject of intermediate-mass black holes (IMBHs) with masses between those of "stellar-mass" and "super-massive" black holes. The existence of IMBHs is a real possibility: they might plausibly have formed as remnants of the first generation of stars (Population III), as the result of dense star cluster evolution, or as part of the formation process of super-massive black holes. Their cosmic mass density could exceed that of super-massive black holes (Omega = 10^{-5.7}) and observations do not even rule out that they may account for all of the baryonic dark matter in the Universe (Omega = 10^{-1.7}). Unambiguous detections of individual IMBHs currently do not exist, but there are observational hints from studies of microlensing events, "ultra-luminous" X-ray sources, and centers of nearby galaxies and globular clusters. Gravitational wave experiments will soon provide another method to probe their existence. IMBHs have potential importance for several fields of astrophysics and are likely to grow as a focus of research attention.

We investigate the stellar disk properties of a sample of 19 nearby spiral galaxies with low inclination and late Hubble type (Scd or later). We combine our high-resolution HST I-band observations with existing ground-based optical images to obtain surface brightness profiles that cover a high dynamic range of galactic radius. Most of these galaxies contain a nuclear star cluster, as discussed in a separate paper. The main goal of the present work is to constrain the properties of stellar bulges at these extremely late Hubble types. We find that the surface brightness profiles of the latest-type spirals are complex, with a wide range in shapes. We have sorted our sample in a sequence, starting with ``pure'' disk galaxies (approximately 30% of the sample). These galaxies have exponential stellar disks that extend inwards to within a few tens of pc from the nucleus, where the light from the nuclear cluster starts to dominate. They appear to be truly bulge-less systems. Progressing along the sequence, the galaxies show increasingly prominent deviations from a simple exponential disk model on kpc scales. Traditionally, such deviations have prompted ``bulge-disk'' decompositions. Indeed, the surface brightness profiles of these galaxies are generally well fit by adding a second (exponential) bulge component. However, we find that most surface brightness profiles can be fit equally well (or better) with a single Sersic-type R^{1/n} profile over the entire radial range of the galaxy, without requiring a separate ``bulge'' component. We warn in a general sense against identification of bulges solely on the basis of single-band surface brightness profiles. (abridged)

We present HST/STIS spectra of the Sbc spiral galaxy NGC 4041 which were used to map the velocity field of the gas in its nuclear region. We detect the presence of a compact (r~0.4" ~40 pc), high surface brightness, rotating nuclear disk co-spatial with a nuclear star cluster. The disk is characterized by a rotation curve with a peak to peak amplitude of ~40 km/s and is systematically blueshifted by 10 - 20 km/s with respect to the galaxy systemic velocity. With the standard assumption of constant mass-to-light ratio and with the nuclear disk inclination taken from the outer disk, we find that a dark point mass of 1(-0.7;+0.6) 10^7 Msun is needed to reproduce the observed rotation curve. However the observed blueshift suggests the possibility that the nuclear disk could be dynamically decoupled. Following this line of reasoning we relax the standard assumptions and find that the kinematical data can be accounted for by the stellar mass provided that either the central mass-to-light ratio is increased by a factor of ~2 or that the inclination is allowed to vary. This model results in a 3 sigma upper limit of 6 10^6 Msun on the mass of any nuclear black hole. Overall, our analysis only allows us to set an upper limit of 2 10^7 Msun on the mass of the nuclear black hole. If this upper limit is taken in conjunction with an estimated bulge B magnitude of -17.7 and with a central stellar velocity dispersion of 95 km/s, then these results are not inconsistent with both the MBH-Lsph and the MBH-sigma correlations. Constraints on black hole masses in spiral galaxies of types as late as Sbc are still very scarce and therefore the present result adds an important new datapoint to our understanding of black hole demography.

We used the HST WFPC2 to obtain I-band images of the centers of 81 brightest cluster galaxies (BCGs), drawn from a volume-limited sample of nearby BCGs. The images show a rich variety of morphological features, including multiple or double nuclei, dust, stellar disks, point source nuclei, and central surface brightness depressions. High resolution surface brightness profiles could be inferred for 60 galaxies. Of those, 88% have well-resolved cores. Twelve percent of the BCG sample lacks a well-resolved core; all but one of these BCGs have ``power-law'' profiles. Some of these galaxies have higher luminosities than any power-law galaxy identified by Faber et al. (1997), and have physical upper limits on the break radius well below the values observed for core galaxies of the same luminosity. These results support the idea that the central structure of early-type galaxies is bimodal in its physical properties, but also suggest that there exist high luminosity galaxies with power-law profiles (or unusually small cores). The BCGs in the latter category tend to fall at the low end of the BCG luminosity function and tend to have low values of the quantity alpha (the logarithmic slope of the metric luminosity as a function of radius, at 10 kpc). Since theoretical calculations have shown that the luminosities and alpha values of BCGs grow with time as a result of accretion, this suggests a scenario in which elliptical galaxies evolve from power-law profiles to core profiles through accretion and merging. This is consistent with theoretical scenarios that invoke the formation of massive black hole binaries during merger events (Abridged).

Keck/ESI spectroscopy is presented for six Virgo Cluster dwarf elliptical (dE) galaxies. The mean line-of-sight velocity and velocity dispersion are resolved as a function of radius along the major axis of each galaxy, nearly doubling the total number of dEs with spatially-resolved stellar kinematics. None of the observed objects shows evidence of strong rotation: upper limits on v_rot/sigma are well below those expected for rotationally-flattened objects. Such limits place strong constraints on dE galaxy formation models. Although these galaxies continue the trend of low rotation velocities observed in Local Group dEs, they are in contrast to recent observations of large rotation velocities in slightly brighter cluster dEs. Using WFPC2 surface photometry and spherically-symmetric dynamical models, we determine global mass-to-light ratios 3 < M/L_V < 6. These ratios are comparable to those expected for an intermediate-age stellar population and are broadly consistent with the (V-I) colors of the galaxies. This implies that these dEs do not have a significant dark matter component inside an effective radius. Central black holes more massive than 10^7 M_sun can be ruled out. For the 5 nucleated dEs in our sample, we determine kinematic and photometric properties for the central nucleus separately from the underlying host dE galaxy. These nuclei are as bright or brighter than the most luminous Galactic globular clusters and lie near the region of Fundamental Plane space occupied by globular clusters. In this space, the Virgo dE galaxies lie in the same general region as Local Group and other nearby dEs, although non-rotating dEs appear to have a slightly higher mean mass and mass-to-light ratio than their rotating counterparts.

We analyze HST/STIS spectra (see Paper I) of the central region of the dense globular cluster M15. We infer the velocities of 64 individual stars, two-thirds of which have their velocity measured for the first time. This triples the number of stars with measured velocities in the central 1 arcsec of M15 and doubles the number in the central 2 arcsec. Combined with existing ground-based data we obtain the radial profiles of the projected kinematical quantities. The RMS velocity sigma_RMS rises to 14 km/s in the central few arcsec, somewhat higher than the values of 10-12 km/s inferred previously from ground-based data. To interpret the results we construct dynamical models based on the Jeans equation, which imply that M15 must have a central concentration of non-luminous material. If this is due to a single black hole, then its mass is M_BH = (3.9 +/- 2.2) x 10^3 solar masses. This is consistent with the relation between M_BH and sigma_RMS that has been established for galaxies. Also, the existence of intermediate-mass black holes in globular clusters is consistent with several scenarios for globular cluster evolution proposed in the literature. Therefore, these results may have important implications for our understanding of the evolution of globular clusters, the growth of black holes, the connection between globular cluster and galaxy formation, and the nature of the recently discovered `ultra-luminous' X-ray sources in nearby galaxies. Instead of a single black hole, M15 could have a central concentration of dark remnants (e.g., neutron stars) due to mass segregation. However, the best-fitting Fokker-Planck models that have previously been constructed for M15 do not predict a central mass concentration that is sufficient to explain the observed kinematics.[ABRIDGED]

It has been reported that there is an error in the figure in Dull et al. (1997, D97) that shows the radial M/L profile in Fokker-Planck models of M15. We discuss how this modifies the interpretation of our kinematical data. These imply the existence of a dark and compact mass component near the center of M15, either a single black hole (BH) or a collection of dark remnants that have sunk to the cluster center due to mass segregation. We previously showed that the latter interpretation is in conflict with the D97 M/L profile, which supported the BH interpretation. We repeat our analysis here with the corrected D97 profile. Models without a BH are now found to be statistically acceptable (within 1-sigma), although inclusion of a BH still provides a marginally better fit. It does not necessarily follow that dark remnants are now the preferred interpretation of the data. The D97 models, as well as N-body models by Baumgardt et al. (2002), assume that all neutron stars are retained during cluster evolution. This conflicts with predictions of the neutron-star retention rate (typically below 10 per cent) based on pulsar kick velocities. The presence of a single BH therefore continues to be a viable interpretation of the data. The best fit BH mass with the corrected D97 M/L profile is M_BH = 1.7^{+2.7}_{-1.7} x 10^3 solar masses, and with a constant M/L it is M_BH = 3.2^{+2.2}_{-2.2} x 10^3 solar masses. A model that includes both neutron star escape and mass segregation would probably yield a value between these numbers. This agrees with the correlation between velocity dispersion and BH mass inferred for galaxies. However, with the presently available models and data it is neither uniquely implied nor ruled out that M15 has an intermediate-mass BH.

In this series of two papers we describe a project with STIS on HST to measure the line-of-sight velocities of stars in the central few arcsec of the dense globular cluster M15. This first paper focuses on the observations and reduction of the data. We `scanned' the central region of M15 spectroscopically by consecutively placing an 0.1 arcsec slit at 18 adjacent positions. Careful and complicated data reduction and analysis were required. The spectral pixel size exceeds the velocity dispersion of M15, which puts the project at the limit of what is feasible with STIS. We applied corrections for the following effects: drifts in the STIS wavelength scale during an HST orbit; the orbital velocity component of HST along the line-of-sight to the cluster, and its variations during the HST orbit; and the apparent wavelength shift that is perceived for a star that is not centered in the slit. To make these corrections we present a detailed study of the HST/STIS point-spread function and line-spread function. We also present a catalog with positions and U, B and V magnitudes of 31,983 stars in M15, created from existing HST/WFPC2 data. We extract 19,200 one-dimensional spectra from the STIS data, each for a different position in M15, and each with a velocity scale accurate to better than 2.5 km/s. We develop an algorithm that co-adds the spectra for individual apertures and use it to extract spectra of individual stars with minimum blending and maximum S/N. In Paper II we use these spectra to extract the line-of-sight velocities of individual stars. The velocities constrain the central structure, dynamics and mass distribution of the cluster, and provide evidence for the existence of an intermediate-mass black hole in the center of the cluster.[ABRIDGED]

(abridged) We investigate the kinematics of the central gas disk of the radio-loud elliptical galaxy NGC 4335, derived from HST/STIS long-slit spectroscopic observations of Halpha+[NII] along 3 parallel slit positions. The observed mean velocities are consistent with a rotating thin disk. We model the gas disk in the customary way. This sets a 3 sigma upper limit of 10^8 Msun on black hole mass, Mbh. The velocity dispersion at r <0.5'' is in excess of that predicted by the thin rotating disk model. This does not invalidate the model, if the excess dispersion is caused by localized turbulent motion in addition to bulk circular rotation. However, if instead the dispersion is caused by the BH potential then a mass Mbh ~ 6x10^8 Msun is inferred by modeling the central gas dispersion as due to an isotropic spherical distribution of collisionless gas cloudlets. The stellar kinematics for NGC 4335 are derived from a ground-based (WHT/ISIS) long-slit observation along the galaxy major axis. A two-integral model of the stellar dynamics yields Mbh >= 3x10^9 Msun. However, there is reason to believe that this model overestimates Mbh. Reported correlations between black hole mass and inner stellar velocity dispersion sigma predict Mbh to be >= 5.4x10^8 Msun in NGC 4335. If our standard thin disk modeling of the gas kinematics is valid, then NGC 4335 has an unusually low Mbh for its velocity dispersion. If, on the other hand, this approach is flawed, and provides an underestimate of Mbh, then black hole masses for other galaxies derived from HST gas kinematics with the same assumptions should be treated with caution.

The E3 giant elliptical galaxy IC1459 is the prototypical galaxy with a fast counterrotating stellar core. We obtained one HST/STIS long-slit spectrum along the major axis of this galaxy and CTIO spectra along five position angles. We present self-consistent three-integral axisymmetric models of the stellar kinematics, obtained with Schwarzschild's numerical orbit superposition method. We study the dynamics of the kinematically decoupled core (KDC) in IC1459 and we find it consists of stars that are well-separated from the rest of the galaxy in phase space. The stars in the KDC counterrotate in a disk on orbits that are close to circular. We estimate that the KDC mass is ~0.5% of the total galaxy mass or ~3*10^9 Msun. We estimate the central black hole mass M_BH of IC1459 independently from both its stellar and its gaseous kinematics. Some complications probably explain why we find rather discrepant BH masses with the different methods. The stellar kinematics suggest that M_BH = (2.6 +/- 1.1)*10^9 Msun (3 sigma error). The gas kinematics suggests that M_BH ~ 3.5*10^8 Msun if the gas is assumed to rotate at the circular velocity in a thin disk. If the observed velocity dispersion of the gas is assumed to be gravitational, then M_BH could be as high as ~1.0*10^9 Msun. These different estimates bracket the value M_BH = (1.1 +/- 0.3)*10^9 Msun predicted by the M_BH-sigma relation. It will be an important goal for future studies to assess the reliability of black hole mass determinations with either technique. This is essential if one wants to interpret the correlation between the BH mass and other global galaxy parameters (e.g. velocity dispersion) and in particular the scatter in these correlations (believed to be only ~0.3 dex). [Abridged]

We used HST WFPC2 images to identify six early-type galaxies with surface- brightness profiles that decrease inward over a limited range of radii near their centers. The implied luminosity density profiles of these galaxies have local minima interior to their core break radii. NGC 3706 harbors a high surface brightness ring of starlight with radius ~20 pc. Its central structure may be related to that in the double-nucleus galaxies M31 and NGC 4486B. NGC 4406 and NGC 6876 have nearly flat cores that on close inspection are centrally depressed. Colors for both galaxies imply that this is not due to dust absorption. The surface brightness distributions of both galaxies are consistent with stellar tori that are more diffuse than the sharply defined system in NGC 3706. The remaining three galaxies are the brightest cluster galaxies in A260, A347, and A3574. Color information is not available for these objects, but they strongly resemble NGC 4406 and NGC 6876 in their cores. The thin ring in NGC 3706 may have formed dissipatively. The five other galaxies resemble the endpoints of some simulations of the merging of two gas-free stellar systems, each harboring a massive nuclear black hole. In one version of this scenario, diffuse stellar tori are produced when stars initially bound to one black hole are tidally stripped away by the second black hole. Alternatively, some inward-decreasing surface-brightness profiles may reflect the ejection of stars from a core during the hardening of the binary black hole created during the merger.

We present dynamical models of the nearby compact elliptical galaxy M32, using high quality kinematical measurements, obtained with the integral-field spectrograph SAURON mounted on the William Herschel Telescope on La Palma. We also include STIS data obtained by Joseph et al. We find a best-fit black hole mass of (2.5 +/- 0.5) million solar masses and a stellar I-band mass-to-light ratio of (1.85 +/- 0.15) in solar units. For the first time, we are also able to constrain the inclination along which M32 is observed to (70 +/- 5) degrees. Combined with an averaged observed flattening of 0.73, this corresponds to an intrinsic flattening of approximately 0.68 +/- 0.03. These tight constraints are mainly caused by the use of integral-field data. We show this quantitatively by comparing with models that are constrained by multiple slits only. We show the phase-space distribution and intrinsic velocity structure of the best-fit model and investigate the effect of regularisation on the orbit distribution.

We derive general expressions for the LMC velocity field which we fit to kinematical data for 1041 carbon stars. We demonstrate that all previous studies of LMC kinematics have made unnecessary over-simplifications that have led to incorrect estimates of important structural parameters. We compile and improve LMC proper motion estimates to support our analysis. We find that the kinematically determined position angle of the line of nodes is 129.9 +/- 6.0 deg. The LMC inclination changes at a rate di/dt = -103 +/- 61 deg/Gyr, a result of precession and nutation induced by Milky Way tidal torques. The LMC rotation curve V(R) has amplitude 49.8 +/- 15.9 km/s, 40% lower than what has previously (and incorrectly) been inferred from e.g. HI. The dynamical center of the carbon stars is consistent with the center of the bar and the center of the outer isophotes, but not with the HI kinematical center. The enclosed mass inside 8.9 kpc is (8.7 +/- 4.3) x 10^9 M_sun, more than half of which is due to a dark halo. The LMC has a larger vertical thickness than has traditionally been believed. Its V/sigma is less than the value for the Milky Way thick disk. We discuss the implications for the LMC self-lensing optical depth. We determine the LMC velocity and orbit in the Galactocentric rest frame and find it to be consistent with the range of velocities that has been predicted by models for the Magellanic Stream. The Milky Way dark halo must have mass >4.3 x 10^{11} M_sun and extent >39 kpc for the LMC to be bound. We predict the LMC proper motion velocity field, and discuss techniques for kinematical distance estimation. [ABRIDGED]

We present new HST I-band images of a sample of 77 nearby, late-type spiral galaxies with low inclination. The main purpose of this catalog is to study the frequency and properties of nuclear star clusters. In 59 galaxies of our sample, we have identified a distinct, compact (but resolved), and dominant source at or very close to their photocenter. In many cases, these clusters are the only prominent source within a few kpc from the galaxy nucleus. We present surface brightness profiles, derived from elliptical isophote fits, of all galaxies for which the fit was successful. We use the fitted isophotes at radii larger than 2" to check whether the location of the cluster coincides with the photocenter of the galaxy, and confirm that in nearly all cases, we are truly dealing with ``nuclear'' star clusters. From analytical fits to the surface brightness profiles, we derive the cluster luminosities after subtraction of the light contribution from the underlying galaxy disk and/or bulge.

Stars and their kinematics provide one of the tools available for studies of the shapes of galaxies and their halos. In this review I focus on two specific applications: the shape of the Milky Way dark halo and the shape of the LMC disk. The former is constrained by a variety of observations, but an accurate determination of the axial ratio q_DH remains elusive. A very flattened Milky Way dark halo with q_DH < 0.4 is ruled out, and values q_DH > 0.7 appear most consistent with the data. Near-IR surveys have revealed that the LMC disk is not approximately circular, as long believed, but instead has an axial ratio of 0.7 in the disk plane. The elongation is perpendicular to the Magellanic Stream, indicating that it is most likely due to the tidal force of the Milky Way. Equilibrium dynamical modeling of galaxies is important for many applications. At the same time, detailed studies of tidal effects and tidal streams have the potential to improve our understanding of both the Milky Way dark halo and the structure of satellite galaxies such as the LMC.

Keck spectroscopy is presented for four dwarf elliptical galaxies in the Virgo Cluster. At this distance, the mean velocity and velocity dispersion are well resolved as a function of radius between 100 to 1000 pc, allowing a clear separation between nuclear and surrounding galaxy light. We find a variety of dispersion profiles for the inner regions of these objects, and show that none of these galaxies is rotationally flattened.

We present a detailed study of the viewing angles of the LMC disk plane. We find that our viewing direction differs considerably from the commonly accepted values, which has important implications for the structure of the LMC. The discussion is based on an analysis of spatial variations in the apparent magnitude of features in the near-IR color-magnitude diagrams extracted from the DENIS and 2MASS surveys. Sinusoidal brightness variations with a peak-to-peak amplitude of approximately 0.25 mag are detected as function of position angle, for both AGB and RGB stars. This is naturally interpreted as the result of distance variations, due to one side of the LMC plane being closer to us than the opposite side. The best fitting geometric model of an inclined plane yields an inclination angle i = 34.7 +/- 6.2 degrees and line-of-nodes position angle Theta = 122.5 +/- 8.3 degrees. There is tentative evidence that the LMC disk plane may be warped. Traditional methods to estimate the position angle of the line of nodes have used either the major axis position angle Theta_maj of the spatial distribution of tracers on the sky, or the position angle Theta_max of the line of maximum gradient in the velocity field, given that for a circular disk Theta_maj = Theta_max = Theta. The present study does not rely on the assumption of circular symmetry, and is considerably more accurate than previous studies of its kind. We find that the actual position angle of the line of nodes differs considerably from both Theta_maj and Theta_max, for which measurements have fallen in the range 140-190 degrees. This indicates that the intrinsic shape of the LMC disk is not circular, but elliptical, as discussed further in Paper II. [Abridged]

I construct a near-IR star count map of the LMC and demonstrate, using the viewing angles derived in Paper I, that the LMC is intrinsically elongated. I argue that this is due to the tidal force from the Milky Way. The near-IR data from the 2MASS and DENIS surveys are used to create a star count map of RGB and AGB stars, which is interpreted through ellipse fitting. The radial number density profile is approximately exponential with a scale-length 1.3-1.5 kpc. However, there is an excess density at large radii that may be due to the tidal effect of the Milky Way. The position angle and ellipticity profile converge to PA_maj = 189.3 +/- 1.4 degrees and epsilon = 0.199 +/- 0.008 for r > 5 deg. At large radii there is a drift of the center of the star count contours towards the near side of the plane, which can be undrestood as due to viewing perspective. The fact that PA_maj differes from the line of nodes position angle Theta = 122.5 +/- 8.3 (cf. Paper I) indicates that the LMC disk is not circular, but has an intrinsic ellipticity of 0.31. The LMC is elongated in the general direction of the Galactic center, and is elongated perpendicular to the Magellanic Stream and the velocity vector of the LMC center of mass. This suggests that the elongation of the LMC has been induced by the tidal force of the Milky Way. The position angle of the line of nodes differs from the position angle Theta_max of the line of maximum line of sight velocity gradient: Theta_max - Theta = 20-60 degrees. This could be due to: (a) streaming along non-circular orbits in the elongated disk; (b) uncertainties in the transverse motion of the LMC center of mass; (c) precession and nutation of the LMC disk as it orbits the Milky Way (expected on theoretical grounds). [Abridged]

We present first results from an HST WFPC2 imaging and STIS spectroscopy program to investigate the structural and star forming properties in the nuclear regions of the Toomre Sequence of merging galaxies. Here we discuss V-band, I-band and H-alpha images of the nuclei. We comment briefly on the connection between the nuclear morphology of the ionized gas and the merger stage.

We present initial results of our analysis of line emission produced in gas disks found at the centers of a sample of nearby, radio galaxies with radio jets. We obtained data using STIS (The Space Telescope Imaging Spectrograph) at three parallel slit positions on the nucleus of each galaxy. This allows us to map the H$\alpha$ + [NII] flux, the gas radial velocity and the velocity dispersion. We find evidence of rotating disks in 11 of the sample galaxies and we can not currently rule out a rotating disk model for the remaining eight. For rotating systems, we find that the minimum central enclosed mass is greater than or similar to the predicted black hole mass based on ground-based stellar velocity dispersions. By modeling the gas dynamics we will go on to constrain the masses of the black holes. We will also investigate the properties of the gas disks themselves, giving us an insight into fueling, ionization mechanisms and the structure of the central regions.

WFPC2 images and STIS spectroscopic observations are presented of the double nucleus in the merger system NGC 6240. We find that: (a) the kinematics of the ionized gas is similar to that of the molecular gas, despite a different morphology; (b) the gaseous and stellar kinematics are quite different, suggesting an early merger stage; (c) neither the gaseous nor the stellar kinematics show an obvious sign of the supermassive black hole believed to be responsible for the X-ray emission of NGC 6240; and (d) the steep off-nuclear velocity gradient is not due to a 10E11 solar mass black hole, in contrast to earlier suggestions.

This review summarizes some aspects of the central kiloparsec scale structure of galaxies, and in particular spiral galaxies, elliptical galaxies and merger remnants. The focus is on results from optical and near-IR imaging and spectroscopy, with emphasis on recent work with the Hubble Space Telescope.

We describe an on-going HST program aimed at determining the relationship between the nuclear black hole mass and bulge mass in spiral galaxies. We have selected a volume limited sample of 54 nearby spiral galaxies for which we already have ground based emission line rotation curves, CCD surface photometry and radio maps. We are now obtaining HST/STIS longslit observations of each of the galaxies in the sample in order to determine the nuclear Halpha rotation curve at high (~0.1") spatial resolution. We will use these data to measure the unresolved dark mass concentration at the nucleus of each object. Here we show the first results from observations of objects in the sample.

We have observed a complete sample of 21 nearby (D < 70h^{-1} Mpc) Fanaroff & Riley Type I galaxies with HST/WFPC2 and detected dust disks and lanes in 19 of them. The radio jets are roughly perpendicular to the dust which is used to constrain the Doppler boosting factors of the radio jet and cores. The VLBA core flux correlates with the central Halpha+[NII] flux which might indicate that the VLBA core is dominated by an isotropic component. Twelve galaxies show nuclear optical sources. We discuss various possible origins for this emission.

We have obtained 1-2 A resolution optical Echellette spectra of the nuclear star cluster in the nearby starburst galaxy NGC 4449. The light is clearly dominated by a very young (6-10 Myr) population of stars. For our age dating, we have used recent population synthesis models to interpret the observed equivalent width of stellar absorption features such as the HI Balmer series and the CaII triplet around 8500 A. We also compare the observed spectrum of the nuclear cluster to synthesized spectra of simple stellar populations of varying ages. All these approaches yield a consistent cluster age. Metallicity estimates based on the relative intensities of various ionization lines yield no evidence for significant enrichment in the center of this low mass galaxy: the metallicity of the nuclear cluster is about one fourth of the solar value, in agreement with independent estimates for the disk material of NGC 4449.

We present a precise determination of the apparent magnitude of the tip of the red giant branch (TRGB) in the I (0.8 micron), J (1.25 micron), and K_S (2.15 micron) bands from the luminosity function of a sample of data extracted from the DENIS catalogue towards the Magellanic Clouds (Cioni et al. 2000). From the J and Ks magnitudes we derive bolometric magnitudes m_bol. We present a new algorithm for the determination of the TRGB magnitude, which we describe in detail and test extensively using Monte-Carlo simulations. We note that any method that searches for a peak in the first derivative (used by most authors) or the second derivative (used by us) of the observed luminosity function does not yield an unbiased estimate for the actual magnitude of the TRGB discontinuity. We stress the importance of correcting for this bias, which is not generally done. We combine the results of our algorithm with theoretical predictions to derive the distance modulus of the Magellanic Clouds. We obtain m-M = 18.55 (0.04 formal, 0.08 systematic) for the Large Magellanic Cloud (LMC), and m-M = 18.99 (0.03 formal, 0.08 systematic) for the Small Magellanic Cloud (SMC). These are among the most accurate determinations of these quantities currently available, which is a direct consequence of the large size of our sample and the insensitivity of near infrared observations to dust extinction.

The peculiar elliptical galaxy IC 1459 (M_V = -21.19, D = 16.5 Mpc) has a fast counterrotating stellar core, stellar shells and ripples, a blue nuclear point source and strong radio core emission. We present results of a detailed HST study of IC 1459, and in particular its central gas disk, aimed a constraining the central mass distribution. We obtained WFPC2 narrow-band imaging centered on the Halpha+[NII] emission lines to determine the flux distribution of the gas emission at small radii, and we obtained FOS spectra at six aperture positions along the major axis to sample the gas kinematics. We construct different dynamical models for the Halpha+[NII] and Hbeta kinematics that include a supermassive black hole, and in which the stellar mass distribution is constrained by the observed surface brightness distribution and ground-based stellar kinematics. All models are consistent with a black hole mass in the range Mbh=1-4 x 10^8 Msun, and models without a black hole are always ruled out at high confidence.

In the context of the CNOC1 cluster survey, redshifts were obtained for galaxies in 16 clusters. The resulting sample is ideally suited for an analysis of the internal velocity and mass distribution of clusters. Previous analyses of this dataset used the Jeans equation to model the projected velocity dispersion profile. However, the results of such an analysis always yield a strong degeneracy between the mass density profile and the velocity dispersion anisotropy profile. Here we analyze the full (R,v) dataset of galaxy positions and velocities in an attempt to break this degeneracy. We build an `ensemble cluster' from the individual clusters under the assumption that they form a homologous sequence. To interpret the data we study a one-parameter family of spherical models with different constant velocity dispersion anisotropy. The best-fit model is sought using a variety of statistics, including the overall likelihood of the dataset. Although the results of our analysis depend slightly on which statistic is used to judge the models, all statistics agree that the best-fit model is close to isotropic. This result derives primarily from the fact that the observed grand-total velocity histogram is close to Gaussian, which is not expected to be the case for a strongly anisotropic model. The best-fitting models have a mass-to-number-density ratio that is approximately independent of radius over the range constrained by the data. They also have a mass-density profile that is consistent with the dark matter halo profile advocated by Navarro, Frenk & White, in terms of both the profile shape and the characteristic scale length. This adds important new weight to the evidence that clusters do indeed follow this proposed universal mass density profile. [Abridged]

Globular clusters could harbor massive central black holes (BHs), just as galaxies do. So far, no unambiguous detection of a massive BH has been reported for any globular cluster. However, the dense core-collapsed cluster M15 seems to be a good candidate. I review the available photometric and kinematic data for this cluster. Both are consistent with a BH of approximately 2000 solar masses, although such a BH is not unambiguously required by the data. I discuss some ongoing studies with Keck and HST which should shed more light on this issue in the coming years.

We present preliminary results of an HST/NICMOS program to image merger remnants in the J, H and K bands. The nuclear brightness profiles for most sample galaxies are similar to those typical for elliptical galaxies, but some (including the well-studied NGC 3921 and 7252) have an unusually high luminosity density at small radii. This is consistent with the prediction of N-body simulations that gas flows to the center during a merger and forms new stars.

We describe an improved, practical method for constructing galaxy models that match an arbitrary set of observational constraints, without prior assumptions about the phase-space distribution function (DF). Our method is an extension of Schwarzschild's orbit superposition technique. As in Schwarzschild's original implementation, we compute a representative library of orbits in a given potential. We then project each orbit onto the space of observables, consisting of position on the sky and line-of-sight velocity, while properly taking into account seeing convolution and pixel binning. We find the combination of orbits that produces a dynamical model that best fits the observed photometry and kinematics of the galaxy. A key new element of this work is the ability to predict and match to the data the full line-of-sight velocity profile shapes. A dark component (such as a black hole and/or a dark halo) can easily be included in the models. We have tested our method, by using it to reconstruct the properties of a two-integral model built with independent software. The test model is reproduced satisfactorily, either with the regular orbits, or with the two-integral components. This paper mainly deals with the technical aspects of the method, while applications to the galaxies M32 and NGC 4342 are described elsewhere (van der Marel et al., Cretton & van den Bosch). (abridged)

We have used the NASA Infra-Red Telescope Facility (IRTF) to observe the nuclear stellar cluster in the nearby, face-on, giant Scd spiral IC 342. From high resolution (R = 21500) spectra at the 12CO (2-0) bandhead at 2.3 micron we derive a line-of-sight stellar velocity dispersion sigma = (33 +- 3) km/s. To interpret this observation we construct dynamical models based on the Jeans equation for a spherical system. The light distribution of the cluster is modeled using an isophotal analysis of an HST V-band image from the HST Data Archive, combined with new ground-based K-band imaging. Under the assumption of an isotropic velocity distribution, the observed kinematics imply a K-band mass-to-light ratio M/L_K = 0.05, and a cluster mass M ~ 6 times 10^6 Msun. We model the mass-to-light ratio with the `starburst99' stellar population synthesis models of Leitherer and collaborators, and infer a best-fitting cluster age in the range 63-630 Myears. Although this result depends somewhat on a number of uncertainties in the modeling (e.g., the assumed extinction along the line-of-sight towards the nucleus, the IMF of the stellar population model, and the velocity dispersion anisotropy of the cluster), none of these can be plausibly modified to yield a significantly larger age. We discuss the implications of this result on possible scenarios for the frequency of nuclear starbursts and their impact on secular evolution of spiral galaxy nuclei. As a byproduct of our analysis, we infer that IC 342 cannot have any central black hole more massive than 0.5 million solar masses. This is ~ 6 times less massive than the black hole inferred to exist in our Galaxy, consistent with the accumulating evidence that galaxies with less massive bulges harbor less massive black holes.

HST observations show that the surface brightness profiles of early-type galaxies have central cusps. I summarize here the results of van der Marel (1999), which show that the observed characteristics of these cusps are consistent with the hypothesis that all early-type galaxies have central black holes (BHs) that grew adiabatically in homogeneous isothermal cores. The models suggest a roughly linear correlation between BH mass and V-band galaxy luminosity, log M_BH = -1.83 + log L in solar units, similar to the relation suggested by kinematical BH detections in nearby galaxies.

The surface brightness profiles of early-type galaxies have central cusps. Two characteristic profile types are observed with HST: `core' profiles have a break at a resolved radius and logarithmic cusp slope gamma < 0.3 inside that radius; `power-law' profiles have no clear break and gamma > 0.3. With few exceptions, galaxies with M_V < -22 have core profiles, and galaxies with M_V > -20.5 have power-law profiles. Both profile types occur in galaxies with -22 < M_V < -20.5. We show that these results are consistent with the hypothesis that: (i) all early-type galaxies have black holes (BHs) that grew adiabatically in homogeneous isothermal cores; and (ii) these `progenitor' cores followed scaling relations similar to those of the fundamental plane. The models studied here are the ones first proposed by Young. Models with BH masses and progenitor cores that obey established scaling relations predict (at Virgo) that galaxies with M_V < -21.2 have core profiles and galaxies with M_V > -21.2 have power-law profiles. This reproduces both the sense and the absolute magnitude of the observed transition. Intrinsic scatter in BH and galaxy properties can explain why both types of galaxies are observed around the transition magnitude. The observed bimodality in cusp slopes may be due to a bimodality in M_bh/L, with rapidly rotating disky galaxies having larger M_bh/L than slowly rotating boxy galaxies. Application to individual galaxies with HST photometry yields a roughly linear correlation between BH mass and V-band galaxy luminosity, log M_bh = -1.83 + log L (solar units). This agrees with the average relation for nearby galaxies with kinematically determined BH masses, and also with predictions from quasar statistics (shortened abstract).

We present an HST study of the nuclear region of the E4 radio galaxy NGC 7052, which has a nuclear disk of dust and gas. The WFPC2 was used to obtain B, V and I broad-band images and an H_alpha+[NII] narrow-band image. The FOS was used to obtain H_alpha+[NII] spectra along the major axis, using a 0.26 arcsec diameter circular aperture. The observed rotation velocity of the ionized gas is V = 155 +/- 17 km/s at r = 0.2 arcsec from the nucleus. The Gaussian dispersion of the emission lines increases from sigma = 70 km/s at r=1 arcsec, to sigma = 400 km/s on the nucleus. To interpret the gas kinematics we construct axisymmetric models in which the gas and dust reside in a disk in the equatorial plane of the stellar body. It is assumed that the gas moves on circular orbits, with an intrinsic velocity dispersion due to turbulence. The circular velocity is calculated from the combined gravitational potential of the stars and a possible nuclear black hole (BH). Models without a BH predict a rotation curve that is shallower than observed (V_pred = 92 km/s at r = 0.2 arcsec), and are ruled out at > 99% confidence. Models with a BH of 3.3^{+2.3}_{-1.3} x 10^8 solar masses provide an acceptable fit. NGC 7052 can be added to the list of active galaxies for which HST spectra of a nuclear gas disk provide evidence for the presence of a central BH. The BH masses inferred for M87, M84, NGC 6251, NGC 4261 and NGC 7052 span a range of a factor 10, with NGC 7052 falling on the low end. By contrast, the luminosities of these galaxies are identical to within 25%. Any relation between BH mass and luminosity, as suggested by independent arguments, must therefore have a scatter of at least a factor 10.

The distribution of black hole masses in galaxies is constrained by photometric and kinematic studies of individual galaxies, and by the properties of the quasar population. I review our understanding of these topics, present new results of adiabatic black hole growth models for HST photometry of elliptical galaxies with brightness profiles of the `core' type, and discuss the implications of ground-based stellar kinematical data. It is not yet possible to uniquely determine the black hole mass distribution, but the available evidence is not inconsistent with a picture in which: (i) a majority of galaxies has black holes; (ii) there is a correlation (with large scatter) between black holes mass M_BH and spheroid luminosity L_sph, approximately of the form M_BH = 10^{-2} L_sph (in solar B-band units); and (iii) the black holes formed in a quasar phase through mass accretion with an efficiency of approximately 0.05.

We present multi-colour (U, V and I) photometry obtained with the second Wide Field and Planetary Camera (WFPC2) on the Hubble Space Telescope (HST), and spectra taken with the William Herschel Telescope (WHT) and the HST Faint Object Spectrograph (FOS), for the early-type galaxies NGC 4342 and NGC 4570. These galaxies are intermediate between ellipticals and lenticulars, and they both have a small nuclear stellar disc in addition to their main outer disc. Colour images reveal no colour differences between the nuclear discs and the bulges. Comparison of the U-V and V-I colours with stellar population models indicates that the central regions of both galaxies are of intermediate age (approx. 8 Gyr) and of high metallicity. The long-slit WHT spectra and the FOS spectra (of four times higher spatial resolution) reveal that both galaxies are rapidly rotating, and they both have velocity dispersions that increase strongly towards the centre. The FOS spectra of NGC 4342 indicate a central velocity dispersion of 420 km/s, higher than the 320 km/s measured from the WHT spectra. Also, the nuclear rotation gradient measured with the FOS is steeper than that measured with the WHT; it reaches V = 200 km/s at 0.25''. The rapid stellar motions seen in the centre of NGC 4342 suggest a large central mass concentration, possibly a massive black hole. The kinematics of the more massive NGC 4570 are less spectacular, with a central velocity dispersion of 250 km/s and a central rotation curve that reaches only 60 km/s at 0.25''.

Axisymmetric dynamical models are constructed for the E3 galaxy M32 to interpret high spatial resolution stellar kinematical HST data. Models are studied with two-integral phase-space distribution functions, and with fully general three- integral distribution functions. The latter are built using a new extension of Schwarzschild's orbit superposition approach. Models are constructed for inclinations of 90 and 55 degrees. No model without a nuclear dark object can fit the combined ground-based and HST data, independent of the dynamical structure of M32. Models with a nuclear dark object of 3.4 x 10^6 solar masses do provide an excellent fit. The inclined models provide the best fit, but the inferred dark mass does not depend sensitively on the assumed inclination. The models that best fit the data are not two-integral models, but like two-integral models they are azimuthally anisotropic. An extended dark object can fit the data only if its half-mass radius is r_h < 0.08 arcsec (=0.26 pc), implying a central dark matter density exceeding 1 x 10^8 solar masses / pc^3. This density is high enough to rule out most plausible alternatives to a massive black hole. The dynamically inferred dark mass is identical to that suggested by existing adiabatic black hole growth models for HST photometry of M32. The low activity of M32 implies either that only a very small fraction of the gas that is shed by evolving stars is accreted onto the black hole, or that accretion proceeds at very low efficiency, e.g. in an advection-dominated mode. (shortened version)

We study the stability of stellar dynamical equilibrium models for M32. Kinematic observations show that M32 has a central black hole of 3x10^6 solar masses, and a phase-space distribution function that is close to the `two-integral' form f=f(E,L_z). M32 is also rapidly rotating; 85-90% of the stars have the same sense of rotation around the symmetry axis. Previous work has shown that flattened, rapidly rotating two-integral models can be bar-unstable. We have performed N-body simulations to test whether this is the case for M32. Particle realizations with N=512,000 were studied for two representative inclinations, i=90 (edge-on) and i=55, corresponding to intrinsic axial ratios of q=0.73 and q=0.55, respectively. The time evolution of the models was calculated with a `self-consistent field' code on a Cray T3D parallel supercomputer. We find both models to be dynamically stable. This implies that they provide a physically meaningful description of M32, and that the inclination of M32 (and hence its intrinsic flattening) cannot be strongly constrained through stability arguments. Previous work on the stability of f(E,L_z) models has shown that the bar-mode is the only possibly unstable mode for systems rounder than q=0.3, and that the likelihood for this mode to be unstable increases with flattening and rotation rate. The f(E,L_z) models studied for M32 are stable, and M32 has a higher rotation rate than nearly all other elliptical galaxies. This suggests that f(E,L_z) models constructed to fit data for real elliptical galaxies will generally be stable for q>0.55, and possibly for flatter systems as well.

We have obtained spectra through small apertures centered on the nuclear region and major axis of M32, with the HST/FOS. A detailed analysis and reduction of the data is presented, including: (i) new calibrations and modeling of the FOS aperture sizes, PSF and line-spread-functions; (ii) determination of the aperture positioning for each observation; and (iii) accurate wavelength calibration, template matching, and kinematical analysis of the spectra. This yields measurements of the stellar rotation velocities and velocity dispersions with five times higher spatial resolution than the best available ground-based data. The inferred velocities provide the highest angular resolution stellar kinematical data obtained to date for any stellar system. The observations show a steeper rotation curve and higher central velocity dispersion than ground-based data. The rotation velocity is observed to be 30 km/s at 0.1" from the nucleus; roughly twice the value measured from the ground at this distance. The nuclear dispersion measured through the smallest FOS aperture (0.068" square) is 156 +/- 10 km/s. The average of four independent dispersion measurements at various positions inside the central 0.1" is 126 km/s, with a RMS scatter of 21 km/s. The nuclear dispersion measured from the ground is only 85-95 km/s, whereas the dispersion outside the central arcsec is only 45-55 km/s. These results significantly strengthen previous arguments for the presence of a massive nuclear black hole in M32. Detailed dynamical models are presented in a series of companion papers.

We describe a powerful technique to model and interpret the stellar line-of-sight velocity profiles of galaxies. Following Schwarzschild's modeling approach, a representative library of orbits is calculated in a given potential; then the non-negative superposition of these orbits is determined to fit best a given set of observational constraints. Our implementation incorporates several new features: (i) we calculate velocity profiles and represent them by a Gauss-Hermite series. This allows us to constrain the orbital anisotropy in the fit. (ii) we take into account the error on each observational constraint to obtain an objective chi2 measure for the quality-of-fit. Only projected, observable quantities are included in the fit, and aperture binning and seeing convolution of the data are properly taken into account. This scheme is valid for any geometry, but here we focus on spherical geometry and the issue of dark halos around elliptical galaxies. We model radially extended velocity profiles of the E0 galaxy NGC 2434, and find that constant M/L models are clearly ruled out, regardless of the orbital anisotropy. To study how much dark matter is needed, we considered a sequence of cosmologically motivated `star+halo' potentials, which are specified by the stellar mass-to-light ratio Gamma and the characteristic halo velocity, V_200 (from Navarro et al. 1996). The star+halo models provide an excellent fit to the data, with Gamma=3.35+-0.25 (in B-band solar units) and V_200=450+-100km/s. The best-fitting potential has a circular velocity Vc that is constant (at ~300km/s) to within 10% between 0.2--3 effective radii. In NGC 2434 roughly half of the mass within an effective radius appears to be dark.

Massive black holes are thought to reside at the centres of many galaxies, where they power quasars and active galactic nuclei. But most galaxies are quiescent, indicating that any central massive black hole present will be starved of fuel and therefore detectable only through its gravitational influence on the motions of the surrounding stars. M32 is a nearby, quiescent elliptical galaxy in which the presence of a black hole has been suspected; however, the limited resolution of the observational data and the restricted classes of models used to interpret this data have made it difficult to rule out alternative explanations, such as models with an anisotropic stellar velocity distribution and no dark mass or models with a central concentration of dark objects (for example, stellar remnants or brown dwarfs). Here we present high-resolution optical HST spectra of M32, which show that the stellar velocities near the centre of this galaxy exceed those inferred from previous ground-based observations. We use a range of general dynamical models to determine a central dark mass concentration of (3.4 +/- 1.6) x 10^6 solar masses, contained within a region only 0.3 pc across. This leaves a massive black hole as the most plausible explanation of the data, thereby strengthening the view that such black holes exist even in quiescent galaxies.

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.

Massive nuclear black holes of 10^6 - 10^9 solar masses are believed to be responsible for the the energy production in quasars and active galaxies, and are thought to be present in many quiescent galaxies as well. Dynamical evidence for this can be sought by studying the dynamics of gas and stars in galactic nuclei at high spatial resolution. This paper reviews the current evidence, with emphasis on some recent developments and ongoing projects. The evidence from water masers and gas kinematics in the active galaxies NGC 4258 and M87 is compelling. In quiescent galaxies only stellar kinematics are generally available. One well-studied case is M32. Stellar dynamical f(E,L_z) models with a few million solar mass black hole fit the ground-based kinematical data remarkably well. N-body simulations of an edge-on f(E,L_z) model for M32 show that this model is stable. HST spectra should soon provide new and improved constraints on the presence of black holes in quiescent galaxies.

A common target acquisition strategy with the Faint Object Spectrograph (FOS) on board the Hubble Space Telescope (HST) is to use a sequence of so-called ACQ/PEAK stages. Standard sequences and exposure times are given in the FOS Instrument Handbook (Keyes et al.~1995). The success of a given strategy is difficult to predict a priori for targets other than point sources. This paper describes software that simulates ACQ/PEAK acquisitions in a Monte-Carlo manner, for a variety of targets. The software is available from the FOS Instrument Scientists or from the author. It allows a proper quantitative assessment to be made of whether an acquisition will succeed, what exposure times should be used, and how the efficiency can be optimized.

A simple algorithm is presented for the calculation of the projected line-of-sight velocity profiles (VPs) of non-rotating anisotropic spherical stellar dynamical models with a phase-space distribution function of the Osipkov-Merritt type. The velocity distribution in these models is isotropic inside the anisotropy radius r_a and becomes increasingly radially biased at larger radii. VP shape parameters are presented for a family of models in which the luminous mass density has an r^{-gamma} power-law cusp at small radii and an r^{-4} power-law fall-off at large radii. Self-consistent models and models in which the luminous matter is embedded in a dark halo are discussed. The effects of changes in the cusp slope gamma and in the anisotropy radius r_a are documented, and the area in the (gamma,r_a)-plane that contains physical models is delineated. The shapes of the VPs of the models show a considerable (and observable) variation with projected galactocentric radius. These models will be useful for interpreting the data on the VP shapes of elliptical galaxies that are now becoming available.

We present measurements of the shape of the stellar line-of-sight velocity distribution out to two effective radii along the major axes of the four elliptical galaxies NGC 2434, 2663, 3706, and 5018. The velocity dispersion profiles are flat or decline gently with radius. We compare the data to the predictions of f=f(E,L_z) axisymmetric models with and without dark matter. Strong tangential anisotropy is ruled out at large radii. We conclude from our measurements that massive dark halos must be present in three of the four galaxies, while for the fourth galaxy (NGC 2663) the case is inconclusive.

The contour integral method of Hunter & Qian is applied to axisymmetric galaxy models in which the distribution function (DF) is of the form f=f(E,L_z), where E and L_z are the classical integrals of motion in an axisymmetric potential. A practical way to construct the unique even part of the DF for such systems is presented. It is applied to models, both oblate and prolate, in which the mass density is stratified on similar concentric spheroids. The spheroids with scale-free densities are discussed in detail. These provide useful approximations to the behaviour of more realistic models in the limit of small and large radii. The self-consistent case is treated, as well as the case in which there are additional contributions to the potential from a central black hole or dark halo. The range of axis ratios and density profile slopes is determined for which spheroidal power--law cusps with a central black hole have a physical two--integral DF. More generally, the f(E,L_z) DFs are discussed for a set of spheroidal `(alpha,beta)-models', characterized by a power-law density cusp with slope alpha at small radii, and a power-law density fall-off with slope [alpha + 2 beta] at large radii. As an application, the DF is constructed for the (alpha,beta) model with a 1.8 x 10^6 solar mass black hole used by van der Marel et al. to interpret their high spatial resolution spectroscopic data for M32. The results confirm that the model fits the observed line-of-sight velocity profiles remarkably well. The model is used to calculate predictions for future spectroscopic observations with the HST.

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.

The dynamical evidence for black holes (BHs) in galactic nuclei is reviewed, with emphasis on recent improvements in spatial resolution, methods for analyzing galaxy spectra and dynamical modeling. M31, M32 and M87 are discussed in some detail.