Comet 49P/Arend-Rigaux is a well known low-activity Jupiter Family comet. Despite the low activity, we have witnessed outgassing activity in 1992, 2004, and 2012. In 2012 a broad tail-like feature (PA$\sim270^\circ, \sim2.3\times10^5$ km) and a narrow jet-like feature (PA$\sim180^\circ, \sim9.3\times10^4$ km) were seen simultaneously. Using Finson-Probstein (FP) dust dynamical models we determine: grain sizes released in each event; duration of activity; when activity peaked; and velocity of the dust particles, allowing us to make comparisons between the events. We find that the tail feature in 2012 is similar to the tail in 1992 with large grains (40-4000 $\mu$m) peaking in activity near perihelion with a long outgassing duration greater than 150 days. The jet feature from 2012, however, is more similar to the 2004 event which we model with small grains (1-8 $\mu$m) with a short duration of activity ($\sim$1 month). The main difference between these two features is that the 2004 event occurs prior to perihelion, while the 2012 event is post-perihelion. We use the grain sizes from the FP models to constrain ice sublimation models. Between 1985 and 2018 we cover 6 apparitions with 26 nights of our own observations plus data from the literature and the Minor Planet Center, which together, allow us to model the heliocentric light curve. We find that the models are consistent with H$_2$O ice sublimation as the volatile responsible for driving activity over most of the active phases and a combination of H$_2$O and CO$_2$ ices are responsible for driving activity near perihelion. We measure the fractional active area over time for H$_2$O and discover that the activity decreases from an average active area of $\sim3\%$ to $\sim0.2\%$. This secular decrease in activity implies that the comet is becoming depleted of volatiles and is in the process of transitioning to a dormant or dead state.

The discovery of the first active interstellar object 2I/Borisov provides an unprecedented opportunity to study planetary formation processes in another planetary system. In particular, spectroscopic observations of 2I allow us to constrain the composition of its nuclear ices. We obtained optical spectra of 2I with the 4.2 m William Herschel and 2.5 m Isaac Newton telescopes between 2019 September 30 and October 13, when the comet was between 2.5 au and 2.4 au from the Sun. We also imaged the comet with broadband filters on 15 nights from September 11 to October 17, as well as with a CN narrow-band filter on October 18 and 20, with the TRAPPIST-North telescope. Broadband imaging confirms that the dust coma colours (B-V=0.82$\pm$0.02, V-R=0.46$\pm$0.03, R-I=0.44$\pm$0.03, B-R=1.28$\pm$0.03) are the same as for Solar System comets. We detect CN emission in all spectra and in the TRAPPIST narrow-band images with production rates between 1.6$\times10^{24}$ and 2.1$\times10^{24}$ molec/s. No other species are detected. We determine three-sigma upper limits for C$_2$, C$_3$, and OH production rates of 6$\times10^{23}$ molec/s, 3$\times10^{23}$ molec/s and 2$\times10^{27}$ molec/s, respectively, on October 02. There is no significant increase of the CN production rate or A(0)f$\rho$ during our observing period. Finally, we place a three-sigma upper limit on the Q(C$_2$)/Q(CN) ratio of 0.3 (on October 13). From this, we conclude that 2I is highly depleted in C$_2$, and may have a composition similar to Solar System carbon-chain depleted comets.

The detection of Interstellar Objects passing through the Solar System offers the promise of constraining the physical and chemical processes involved in planetary formation in other extrasolar systems. While the effect of outgassing by 1I/2017 U1 ('Oumuamua) was dynamically observed, no direct detection of the ejected material was made. The discovery of the active interstellar comet 2I/Borisov means spectroscopic investigations of the sublimated ices is possible for this object. We report the first detection of gas emitted by an interstellar comet via the near-UV emission of CN from 2I/Borisov at a heliocentric distance of $r$ = 2.7 au on 2019 September 20. The production rate was found to be Q(CN) = $(3.7\pm0.4)\times10^{24}$ s$^{-1}$, using a simple Haser model with an outflow velocity of 0.5 km s$^{-1}$. No other emission was detected, with an upper limit to the production rate of C$_2$ of $4\times10^{24}$ s$^{-1}$. The spectral reflectance slope of the dust coma over $3900$ \AA\ $< \lambda< 6000$ \AA \ is steeper than at longer wavelengths, as found for other comets. Broad band $R_c$ photometry on 2019 September 19 gave a dust production rate of $Af\rho=143\pm10$ cm. Modelling of the observed gas and dust production rates constrains the nuclear radius to $0.7-3.3$ km assuming reasonable nuclear properties. Overall, we find the gas, dust and nuclear properties for the first active Interstellar Object are similar to normal Solar System comets.

We report on the catastrophic disintegration of P/2016 G1 (PANSTARRS), an active asteroid, in April 2016. Deep images over three months show that object is constituted by a central concentration of fragments surrounded by an elongated coma, and presents previously unreported sharp arc-like and a narrow linear features. The morphology and evolution of these characteristics independently point toward a brief event on 2016 March 6. The arc and the linear feature can be reproduced by large particles on a ring, moving at ~2.5 m/s. The expansion of the ring defines a cone with a ~40deg half-opening. We propose that G1 was hit by a small object which caused its (partial or total) disruption, and that the ring corresponds to large fragments ejected during the final stages of the crater formation.}

On 2019 January 5 a streamer associated with the 4--10 km main-belt asteroid (6478)~Gault was detected by the ATLAS sky survey, a rare discovery of activity around a main-belt asteroid. Archival data from ATLAS and Pan-STARRS1 show the trail in early December 2018, but not between 2010 and January 2018. The feature has significantly changed over one month, perfectly matching predictions of pure dust dynamical evolution and changes in observing geometry for a short release of dust around 2018 October 28. Follow-up observations with HST show a second narrow trail corresponding to a brief release of dust on 2018 December 30. Both releases occurred with negligible velocity. We find the dust grains to be fairly large, with power-law size distributions in the $10^{-5} - 10^{-3}$~m range and power-law indices of $\sim -1.5$. Three runs of ground-based data find a signature of $\sim 2\,\rm h$ rotation, close to the rotational limit, suggesting that the activity is the result of landslides or reconfigurations after YORP spin-up.

We show that `Oumuamua's excited spin could be in a high energy LAM state, which implies that its shape could be far from the highly elongated shape found in previous studies. CLEAN and ANOVA algorithms are used to analyze `Oumuamua's lightcurve using 818 observations over 29.3~days. Two fundamental periodicities are found at frequencies (2.77$\pm$0.11) and (6.42$\pm$0.18)~cycles/day, corresponding to (8.67$\pm$0.34)~h and (3.74$\pm$0.11)~h, respectively. The phased data show that the lightcurve does not repeat in a simple manner, but approximately shows a double minimum at 2.77~cycles/day and a single minimum at 6.42~cycles/day. This is characteristic of an excited spin state. `Oumuamua could be spinning in either the long (LAM) or short (SAM) axis mode. For both, the long axis precesses around the total angular momentum vector with an average period of (8.67$\pm$0.34)~h. For the three LAMs we have found, the possible rotation periods around the long axis are 6.58, 13.15, or 54.48~h, with 54.48~h being the most likely. `Oumuamua may also be nutating with respective periods of half of these values. We have also found two possible SAM states where `Oumuamua oscillates around the long axis with possible periods at 13.15 and 54.48~h, the latter as the most likely. In this case any nutation will occur with the same periods. Determination of the spin state, the amplitude of the nutation, the direction of the TAMV, and the average total spin period may be possible with a direct model fit to the lightcurve. We find that `Oumuamua is "cigar-shaped"', if close to its lowest rotational energy, and an extremely oblate spheroid if close to its highest energy state for its total angular momentum.

Comet C/2017 K2 (PANSTARRS) was discovered by the Pan-STARRS1 (PS1) Survey on 2017 May 21 at a distance 16.09 au from the Sun, the second most distant discovery of an active comet. Pre-discovery images in the PS1 archive back to 2014 and additional deep CFHT images between 2013 May 10-13 showed the comet to be active at 23.75 au. We derive an upper limit to the nucleus radius of $R_N$=80 km, assuming a 4\% albedo. The spectral reflectivity of the comet surface is similar to "fresh" regions seen on comet 67P/Churyumov-Gerasimenko using the $Rosetta$ OSIRIS camera. Pre-discovery photometry combined with new data obtained with Megacam on the CFHT show that the activity is consistent with CO-ice sublimation and inconsistent with CO$_2$-ice sublimation. The ice sublimation models were run out to perihelion in 2022 at 1.8 au to predict the CO production rates, assuming that the outgassing area does not change. Assuming a canonical 4\% active surface area for water-ice sublimation, we present production rate ratios, $Q_{\rm CO}$/$Q_{\rm H2O}$, for a range of nucleus sizes. Comparing these results with other CO-rich comets we derive a lower limit to the nucleus radius of $\sim$14 km. We present predictions for $Q_{\rm CO}$ at a range of distances that will be useful for planning observations with JWST and large ground-based facilities.

We present a summary of the campaign of remote observations that supported the European Space Agency's Rosetta mission. Telescopes across the globe (and in space) followed comet 67P/Churyumov-Gerasimenko from before Rosetta's arrival until nearly the end of mission in September 2016. These provided essential data for mission planning, large-scale context information for the coma and tails beyond the spacecraft, and a way to directly compare 67P with other comets. The observations revealed 67P to be a relatively `well behaved' comet, typical of Jupiter family comets and with activity patterns that repeat from orbit-to-orbit. Comparison between this large collection of telescopic observations and the in situ results from Rosetta will allow us to better understand comet coma chemistry and structure. This work is just beginning as the mission ends -- in this paper we present a summary of the ground-based observations and early results, and point to many questions that will be addressed in future studies.

We present initial time-resolved observations of the split comet 332P/Ikeya-Murakami taken using the Hubble Space Telescope. Our images reveal a dust-bathed cluster of fragments receding from their parent nucleus at projected speeds in the range 0.06 to 3.5 m s$^{-1}$ from which we estimate ejection times from October to December 2015. The number of fragments with effective radii $\gtrsim$20 m follows a differential power law with index $\gamma$ = -3.6$\pm$0.6, while smaller fragments are less abundant than expected from an extrapolation of this power-law. We argue that, in addition to losses due to observational selection, torques from anisotropic outgassing are capable of destroying the small fragments by driving them quickly to rotational instability. Specifically, the spin-up times of fragments $\lesssim$20 m in radius are shorter than the time elapsed since ejection from the parent nucleus. The effective radius of the parent nucleus is $r_e \le$ 275 m (geometric albedo 0.04 assumed). This is about seven times smaller than previous estimates and results in a nucleus mass at least 300 times smaller than previously thought. The mass in solid pieces, $2\times10^9$ kg, is about 4% of the mass of the parent nucleus. As a result of its small size, the parent nucleus also has a short spin-up time. Brightness variations in time-resolved nucleus photometry are consistent with rotational instability playing a role in the release of fragments.

We describe 2016 January to April observations of the fragments of 332P/Ikeya-Murakami, a comet earlier observed in a 2010 October outburst (Ishiguro et al 2014). We present photometry of the fragments, and perform simulations to infer the time of breakup. We argue that the eastern-most rapidly brightening fragment ($F4$) best corresponds to the original nucleus, rather than the initial bright fragment $F1$. We compute radial and tangential non-gravitational parameters, $A_1 = (1.5 \pm 0.4) \times 10^{-8}$ AU day$^{-2}$ and $(7.2 \pm 1.9) \times 10^{-9}$ AU day$^{-2}$; both are consistent with zero at the $4\sigma$ level. Monte Carlo simulations indicate that the fragments were emitted on the outbound journey well after the 2010 outburst, with bright fragment $F1$ splitting in mid--2013 and the fainter fragments within months of the 2016 January recovery. Western fragment $F7$ is the oldest, dating from 2011. We suggest that the delayed onset of the splitting is consistent with a self-propagating crystallization of water ice.

We present detailed chemical element abundance ratios of 17 elements in three metal poor stars in the Ursa Minor dwarf spheroidal galaxy, which we combine with extant data from the literature to assess the predictions of a novel suite of galaxy chemical evolution models. The spectroscopic data were obtained with the Keck/HIRES instrument and revealed low metallicities of [Fe/H]=-2.12, -2.13 and -2.67 dex. While the most metal poor star in our sample shows an overabundance of [Mn/Fe] and other Fe-peak elements, our overall findings are in agreement with previous studies of this galaxy: elevated values of the [alpha/Fe] ratios that are similar to, or only slightly lower than, the halo values but with SN Ia enrichment at very low metallicity, as well as an enhancement of the ratio of first to second peak neutron capture elements [Y/Ba] with decreasing metallicity. The chemical evolution models which were tailored to reproduce the metallicity distribution function of the dSph, indicate that UMi had an extended star formation which lasted nearly 5 Gyr with low efficiency and are able to explain the [Y/Ba] enhancement at low metallicity for the first time. In particular, we show that the present day lack of gas is probably due to continuous loss of gas from the system, which we model as winds.

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. We find that the object is clearly active both in observations obtained in 2014 and in precovery observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggesting that its activity is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in the total brightness of the comet over the 2014 observing period, and dust modeling results showing that the dust emission persists over at least three months during both active periods, where we find start dates for emission no later than 2003 July 24+/-10 for the 2003 active period and 2014 July 28+/-10 for the 2014 active period. From serendipitous observations by the Subaru Telescope in 2004 when the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude of H_R=17.1+/-0.3, corresponding to an effective nucleus radius of r_e~1.00+/-0.15 km. The object's faintness at that time means we cannot rule out the presence of activity, and so this computed radius should be considered an upper limit. We find that 313P's orbit is intrinsically chaotic, having a Lyapunov time of T_l=12000 yr and being located near two 3-body mean-motion resonances with Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 T1 (PANSTARRS), to belong to the ~155 Myr old Lixiaohua asteroid family.

We analyze 760475 observations of 333026 main-belt objects obtained by the Pan-STARRS1 (PS1) survey telescope between 2012 May 20 and 2013 November 9, a period during which PS1 discovered two main-belt comets, P/2012 T1 (PANSTARRS) and P/2013 R3 (Catalina-PANSTARRS). PS1 comet detection procedures currently consist of the comparison of the point spread functions (PSFs) of moving objects to those of reference stars, and the flagging of objects that show anomalously large radial PSF widths. Based on the number of missed discovery opportunities among comets discovered by other observers, we estimate an upper limit comet discovery efficiency rate of ~70% for PS1. Additional analyses that could improve comet discovery yields in future surveys include linear PSF analysis, modeling of trailed stellar PSFs for comparison to trailed moving object PSFs, searches for azimuthally localized activity, comparison of point-source-optimized photometry to extended-source-optimized photometry, searches for photometric excesses in objects with known absolute magnitudes, and crowd-sourcing. Analysis of PS1 survey statistics indicates an expected fraction of 59 MBCs per 10^6 outer main-belt asteroids, and a 95% confidence upper limit of 96 MBCs per 10^6 outer main-belt asteroids. We note that more sensitive future surveys could detect many more MBCs than estimated here. We find an excess of high eccentricities (0.1 < e < 0.3) among all known MBCs relative to the background asteroid population. Theoretical calculations show that, given these eccentricities, the sublimation rate for a typical MBC is orders of magnitude larger at perihelion than at aphelion, providing a plausible physical explanation for the observed behavior of MBCs peaking in observed activity strength near perihelion.

We have calculated 90% confidence limits on the steady-state rate of catastrophic disruptions of main belt asteroids in terms of the absolute magnitude at which one catastrophic disruption occurs per year (HCL) as a function of the post-disruption increase in brightness (delta m) and subsequent brightness decay rate (tau). The confidence limits were calculated using the brightest unknown main belt asteroid (V = 18.5) detected with the Pan-STARRS1 (Pan-STARRS1) telescope. We measured the Pan-STARRS1's catastrophic disruption detection efficiency over a 453-day interval using the Pan-STARRS moving object processing system (MOPS) and a simple model for the catastrophic disruption event's photometric behavior in a small aperture centered on the catastrophic disruption event. Our simplistic catastrophic disruption model suggests that delta m = 20 mag and 0.01 mag d-1 < tau < 0.1 mag d-1 which would imply that H0 = 28 -- strongly inconsistent with H0,B2005 = 23.26 +/- 0.02 predicted by Bottke et al. (2005) using purely collisional models. We postulate that the solution to the discrepancy is that > 99% of main belt catastrophic disruptions in the size range to which this study was sensitive (100 m) are not impact-generated, but are instead due to fainter rotational breakups, of which the recent discoveries of disrupted asteroids P/2013 P5 and P/2013 R3 are probable examples. We estimate that current and upcoming asteroid surveys may discover up to 10 catastrophic disruptions/year brighter than V = 18.5.

We present the results of a search for the reactivation of active asteroid 176P/LINEAR during its 2011 perihelion passage using deep optical observations obtained before, during, and after that perihelion passage. Deep composite images of 176P constructed from data obtained between June 2011 and December 2011 show no visible signs of activity, while photometric measurements of the object during this period also show no significant brightness enhancements similar to that observed for 176P between November 2005 and December 2005 when it was previously observed to be active. An azimuthal search for dust emission likewise reveals no evidence for directed emission (i.e., a tail, as was previously observed for 176P), while a one-dimensional surface brightness profile analysis shows no indication of a spherically symmetric coma at any time in 2011. We conclude that 176P did not in fact exhibit activity in 2011, at least not on the level on which it exhibited activity in 2005, and suggest that this could be due to the devolatization or mantling of the active site responsible for its activity in 2005.

We report photometric observations for comet C/2012 S1 (ISON) obtained during the time period immediately after discovery (r=6.28 AU) until it moved into solar conjunction in mid-2013 June using the UH2.2m, and Gemini North 8-m telescopes on Mauna Kea, the Lowell 1.8m in Flagstaff, the Calar Alto 1.2m telescope in Spain, the VYSOS-5 telescopes on Mauna Loa Hawaii and data from the CARA network. Additional pre-discovery data from the Pan STARRS1 survey extends the light curve back to 2011 September 30 (r=9.4 AU). The images showed a similar tail morphology due to small micron sized particles throughout 2013. Observations at sub-mm wavelengths using the JCMT on 15 nights between 2013 March 9 (r=4.52 AU) and June 16 (r=3.35 AU) were used to search for CO and HCN rotation lines. No gas was detected, with upper limits for CO ranging between (3.5-4.5)E27 molec/s. Combined with published water production rate estimates we have generated ice sublimation models consistent with the photometric light curve. The inbound light curve is likely controlled by sublimation of CO2. At these distances water is not a strong contributor to the outgassing. We also infer that there was a long slow outburst of activity beginning in late 2011 peaking in mid-2013 January (r~5 AU) at which point the activity decreased again through 2013 June. We suggest that this outburst was driven by CO injecting large water ice grains into the coma. Observations as the comet came out of solar conjunction seem to confirm our models.

We present initial results from observations and numerical analyses aimed at characterizing main-belt comet P/2012 T1 (PANSTARRS). Optical monitoring observations were made between October 2012 and February 2013 using the University of Hawaii 2.2 m telescope, the Keck I telescope, the Baade and Clay Magellan telescopes, Faulkes Telescope South, the Perkins Telescope at Lowell Observatory, and the Southern Astrophysical Research (SOAR) telescope. The object's intrinsic brightness approximately doubles from the time of its discovery in early October until mid-November and then decreases by ~60% between late December and early February, similar to photometric behavior exhibited by several other main-belt comets and unlike that exhibited by disrupted asteroid (596) Scheila. We also used Keck to conduct spectroscopic searches for CN emission as well as absorption at 0.7 microns that could indicate the presence of hydrated minerals, finding an upper limit CN production rate of QCN<1.5x10^23 mol/s, from which we infer a water production rate of QH2O<5x10^25 mol/s, and no evidence of the presence of hydrated minerals. Numerical simulations indicate that P/2012 T1 is largely dynamically stable for >100 Myr and is unlikely to be a recently implanted interloper from the outer solar system, while a search for potential asteroid family associations reveal that it is dynamically linked to the ~155 Myr-old Lixiaohua asteroid family.

We describe the Pan-STARRS Moving Object Processing System (MOPS), a modern software package that produces automatic asteroid discoveries and identifications from catalogs of transient detections from next-generation astronomical survey telescopes. MOPS achieves > 99.5% efficiency in producing orbits from a synthetic but realistic population of asteroids whose measurements were simulated for a Pan-STARRS4-class telescope. Additionally, using a non-physical grid population, we demonstrate that MOPS can detect populations of currently unknown objects such as interstellar asteroids. MOPS has been adapted successfully to the prototype Pan-STARRS1 telescope despite differences in expected false detection rates, fill-factor loss and relatively sparse observing cadence compared to a hypothetical Pan-STARRS4 telescope and survey. MOPS remains >99.5% efficient at detecting objects on a single night but drops to 80% efficiency at producing orbits for objects detected on multiple nights. This loss is primarily due to configurable MOPS processing limits that are not yet tuned for the Pan-STARRS1 mission. The core MOPS software package is the product of more than 15 person-years of software development and incorporates countless additional years of effort in third-party software to perform lower-level functions such as spatial searching or orbit determination. We describe the high-level design of MOPS and essential subcomponents, the suitability of MOPS for other survey programs, and suggest a road map for future MOPS development.

We describe the recovery of faint Main Belt comet P/2008 R1 Garradd using several telescopes, culminating in a successful low $S/N$ recovery with the Gemini North telescope with GMOS. This recovery was a time-critical effort for a mission proposal, and had to be performed in a crowded field. We describe techniques and software tools for eliminating systematic noise artifacts and stellar residuals, bringing the final detection image statistics close to the Gaussian ideal for a median image stack, and achieving a detection sensitivity close to this theoretical optimum. The magnitude of $R_c$=26.1$\pm$0.2 with an assumed geometric albedo of 0.05 corresponds to a radius of 0.3 km. For ice to have survived in this object over the age of the solar system, it implies that the object is a more recent collisional fragment. We discuss the implications of the unexpectedly faint magnitude and nuclear size of P/2008 R1 on the survival of ice inside very small bodies.

P/2010 A2 is an object on an asteroidal orbit that was observed to have an extended tail or debris trail in January 2010. In this work, we fit the outburst of P/2010 A2 with a conical burst model, and verify previous suspicions that this was a one--time collisional event rather than an sustained cometary outburst, implying that P/2010 A2 is not a new Main Belt Comet driven by ice sublimation. We find that the best--fit cone opening angle is about 40 to 50 degrees, in agreement with numerical and laboratory simulations of cratering events. Mapping debris orbits to sky positions suggests that the distinctive arc features in the debris correspond to the same debris cone inferred from the extended dust. From the velocity of the debris, and from the presence of a velocity maximum at around 15 cm/s, we infer that the surface of A2 probably has a very low strength (<1 kPa), comparable to lunar regolith.

We present observations of comet-like main-belt object P/2010 R2 (La Sagra) obtained by Pan-STARRS 1 and the Faulkes Telescope-North on Haleakala in Hawaii, the University of Hawaii 2.2 m, Gemini-North, and Keck I telescopes on Mauna Kea, the Danish 1.54 m telescope at La Silla, and the Isaac Newton Telescope on La Palma. An antisolar dust tail is observed from August 2010 through February 2011, while a dust trail aligned with the object's orbit plane is also observed from December 2010 through August 2011. Assuming typical phase darkening behavior, P/La Sagra is seen to increase in brightness by >1 mag between August 2010 and December 2010, suggesting that dust production is ongoing over this period. These results strongly suggest that the observed activity is cometary in nature (i.e., driven by the sublimation of volatile material), and that P/La Sagra is therefore the most recent main-belt comet to be discovered. We find an approximate absolute magnitude for the nucleus of H_R=17.9+/-0.2 mag, corresponding to a nucleus radius of ~0.7 km, assuming an albedo of p=0.05. Using optical spectroscopy, we find no evidence of sublimation products (i.e., gas emission), finding an upper limit CN production rate of Q_CN<6x10^23 mol/s, from which we infer an H2O production rate of Q_H2O<10^26 mol/s. Numerical simulations indicate that P/La Sagra is dynamically stable for >100 Myr, suggesting that it is likely native to its current location and that its composition is likely representative of other objects in the same region of the main belt, though the relatively close proximity of the 13:6 mean-motion resonance with Jupiter and the (3,-2,-1) three-body mean-motion resonance with Jupiter and Saturn mean that dynamical instability on larger timescales cannot be ruled out.

Main belt asteroid (300163) 2006 VW139 (later designated P/2006 VW139) was discovered to exhibit comet-like activity by the Pan-STARRS1 survey telescope using automated point-spread-function analyses performed by PS1's Moving Object Processing System. Deep follow-up observations show both a short (\sim 10") antisolar dust tail and a longer (\sim 60") dust trail aligned with the object's orbit plane, similar to the morphology observed for another main-belt comet, P/2010 R2 (La Sagra), and other well-established comets, implying the action of a long-lived, sublimation-driven emission event. Photometry showing the brightness of the near-nucleus coma remaining constant over \sim 30 days provides further evidence for this object's cometary nature, suggesting it is in fact a main-belt comet, and not a disrupted asteroid. A spectroscopic search for CN emission was unsuccessful, though we find an upper limit CN production rate of Q_CN < 1.3x10^24 mol/s, from which we infer a water production rate of Q_H2O < 10^26 mol/s. We also find an approximately linear optical spectral slope of 7.2%/1000A, similar to other cometary dust comae. Numerical simulations indicate that P/2006 VW139 is dynamically stable for > 100 Myr, while a search for a potential asteroid family around the object reveals a cluster of 24 asteroids within a cutoff distance of 68 m/s. At 70 m/s, this cluster merges with the Themis family, suggesting that it could be similar to the Beagle family to which another main-belt comet, 133P/Elst-Pizarro, belongs.

Comet P/2010A2 LINEAR is a good candidate for membership with the Main Belt Comet family. It was observed with several telescopes (ESO NTT, La Silla; Gemini North, Mauna Kea; UH 2.2m, Mauna Kea) from 14 Jan. until 19 Feb. 2010 in order to characterize and monitor it and its very unusual dust tail, which appears almost fully detached from the nucleus; the head of the tail includes two narrow arcs forming a cross. The immediate surroundings of the nucleus were found dust-free, which allowed an estimate of the nucleus radius of 80-90m. A model of the thermal evolution indicates that such a small nucleus could not maintain any ice content for more than a few million years on its current orbit, ruling out ice sublimation dust ejection mechanism. Rotational spin-up and electrostatic dust levitations were also rejected, leaving an impact with a smaller body as the favoured hypothesis, and ruling out the cometary nature of the object. The impact is further supported by the analysis of the tail structure. Finston-Probstein dynamical dust modelling indicates the tail was produced by a single burst of dust emission. More advanced models, independently indicate that this burst populated a hollow cone with a half-opening angle alpha~40degr and with an ejection velocity v_max ~ 0.2m/s, where the small dust grains fill the observed tail, while the arcs are foreshortened sections of the burst cone. The dust grains in the tail are measured to have radii between a=1-20mm, with a differential size distribution proportional to a^(-3.44 +/- 0.08). The dust contained in the tail is estimated to at least 8x10^8kg, which would form a sphere of 40m radius. Analysing these results in the framework of crater physics, we conclude that a gravity-controlled crater would have grown up to ~100m radius, i.e. comparable to the size of the body. The non-disruption of the body suggest this was an oblique impact.

The seven known main belt comets (MBCs) have orbital characteristics of main belt asteroids yet exhibit dust ejection like comets. In order to constrain their physical and orbital properties we searched the Thousand Asteroid Light Curve Survey (TALCS; Masiero et al. 2009) for additional candidates using two diagnostics: tail and coma detection. This was the most sensitive MBC survey effort to date, extending the search from MBCs with H~18 (D~1 km) to H~21 (D~150 m). We fit each of the 924 TALCS objects to a PSF model incorporating both a coma and nuclear component to measure the fractional contribution of the coma to the total surface brightness. We determined the significance of the coma detection using the same algorithm on a sample of comparable null detections. We did not identify any MBC candidates with this technique to a sensitivity limit on the order of cometary mass loss rate of about 0.1 kg/s. Our tail detection algorithm identified statistically significant flux in a segmented annulus around the candidate object. We show that the technique can detect tail activity throughout the asteroid belt to the level of the currently known MBCs. Although we did not identify any MBC candidates with this technique, we find a statistically significant detection of faint activity in the entire ensemble of TALCS asteroids. This suggests that many main belt asteroids are active at very low levels. We set 90% upper confidence limits on the number distribution of MBCs as a function of absolute magnitude, semimajor axis, eccentricity, and inclination. There are <~ 400000 MBCs in the main belt brighter than H_V=21 (~150 m) and the MBC:MBA ratio is <~ 1:400. We further comment on the ability of observations to meaningfully constrain the snow line's location. Under some reasonable and simple assumptions we claim 85% confidence that the contemporary snow line lies beyond 2.5 AU.

We present results from a wide-field imaging campaign at the Canada-France-Hawaii Telescope to study the spectacular outburst of comet 17P/Holmes in late 2007. Using image-processing techniques we probe inside the spherical dust coma and find sixteen fragments having both spatial distribution and kinematics consistent with isotropic ejection from the nucleus. Photometry of the fragments is inconsistent with scattering from monolithic, inert bodies. Instead, each detected fragment appears to be an active cometesimal producing its own dust coma. By scaling from the coma of the primary nucleus of 17P/Holmes, assumed to be 1.7 km in radius, we infer that the sixteen fragments have maximum effective radii between ~ 10 m and ~ 100 m on UT 2007 Nov. 6. The fragments subsequently fade at a common rate of ~ 0.2 mag/day, consistent with steady depletion of ices from these bodies in the heat of the Sun. Our characterization of the fragments supports the hypothesis that a large piece of material broke away from the nucleus and crumbled, expelling smaller, icy shards into the larger dust coma around the nucleus.

Aim: We establish the mean metallicity from high-resolution spectroscopy for the recently found dwarf spheroidal galaxy Bootes I and test whether it is a common feature for ultra-faint dwarf spheroidal galaxies to show signs of inhomogeneous chemical evolution (e.g. as found in the Hercules dwarf spheroidal galaxy). Methods: We analyse high-resolution, moderate signal-to-noise spectra for seven red giant stars in the Bootes I dSph galaxy using standard abundance analysis techniques. In particular, we assume local thermodynamic equilibrium and employ spherical model atmospheres and codes that take the sphericity of the star into account when calculating the elemental abundances. Results: We confirm previous determinations of the mean metallicity of the Bootes I dwarf spheroidal galaxy to be -2.3 dex. Whilst five stars are clustered around this metallicity, one is significantly more metal-poor, at -2.9 dex, and one is more metal-rich at, -1.9 dex. Additionally, we find that one of the stars, Boo-127, shows an atypically high [Mg/Ca] ratio, indicative of stochastic enrichment processes within the dSph galaxy. Similar results have previously only been found in the Hercules and Draco dSph galaxies and appear, so far, to be unique to this type of galaxy.

We present a new suite of photometric and spectroscopic data for the faint Bootes II dwarf spheroidal galaxy candidate. Our deep photometry, obtained with the INT/WFC, suggests a distance of 46 kpc and a small half-light radius of 4.0 arcmin (56 pc), consistent with previous estimates. Follow-up spectroscopy obtained with the Gemini/GMOS instrument yielded radial velocities and metallicities. While the majority of our targets covers a broad range in velocities and metallicities, we find five stars which share very similar velocities and metallicities and which are all compatible with the colors and magnitudes of the galaxy's likely red giant branch. We interpret these as a spectroscopic detection of the Bootes II system. These stars have a mean velocity of -117 km/s, a velocity dispersion of (10.5+-7.4) km/s and a mean [Fe/H] of -1.79 dex, with a dispersion of 0.14 dex. At this metallicity, Boo II is not consistent with the stellar-mass-metallicity relation for the more luminous dwarf galaxies. Coupled with our distance estimate, its high negative systemic velocity rules out any physical connection with its projected neighbor, the Bootes I dwarf spheroidal, which has a velocity of ~+100 km/s. The velocity and distance of Bootes II coincide with those of the leading arm of Sagittarius, which passes through this region of the sky, so that it is possible that Bootes II may be a stellar system associated with the Sagittarius stream. Finally, we note that the properties of Bootes II are consistent with being the surviving remnant of a previously larger and more luminous dSph galaxy.

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

There has been a vast recent improvement in photometric and kinematic data for star clusters, Ultra Compact dwarfs, galactic nuclei, and local dSph galaxies, with Subaru contributing substantially to the photometric studies in particular. These data show that there is a bimodal distribution in half-light radii, with stable star clusters always being smaller than 35pc, while stable galaxies are always larger than 120pc. We extend the previously known observational relationships and interpret them in terms of a more fundamental pair of intrinsic properties of dark matter itself: dark matter forms cored mass distributions, with a core scale length of greater than about 100pc, and always has a maximum central mass density with a narrow range. The dark matter in dSph galaxies appears to be clustered such that there is a mean volume mass density within the stellar distribution which has the very low value of about 0.1$M_{\odot}$ pc$^{-3}$. None of the dSphs displays kinematics which require the presence of an inner cusp, while in two dSphs there is evidence that the density profile is shallow (cored) in the inner regions. The maximum central dark matter density derived is model dependent, but is likely to have a mean value (averaged over a volume of radius 10pc) of about 0.1$M_{\odot}$ pc$^{-3}$, which is 5GeV/c$^2$cm$^{-3}$). Galaxies are embedded in dark matter halos with these properties; smaller systems containing dark matter are not observed.

We have obtained high-resolution spectroscopy of ten red giants in the Carina dwarf spheroidal (dSph) with UVES at the ESO/VLT. Here we present the abundances of O,Na,Mg,Si,Ca,Ti and Fe. By comparing the iron abundances [Fe/H] with calcium triplet (CaT) metallicities we show that the empirical CaT technique yields good agreement with the high-resolution data for [Fe/H]>-2 dex, but tends to deviate at lower metallicities. We identify two metal poor stars with iron abundances of -2.72 and -2.50 dex. These stars are found to have enhanced [alpha/Fe] ratios similar to those of stars in the Milky Way halo. However, the bulk of the Carina red giants are depleted in the [alpha/Fe] abundance ratios with respect to the Galactic halo at a given metallicity. One of our targets, with a [Fe/H] of -1.5 dex, is considerably depleted in almost all of the alpha-elements by ~0.5 dex compared to the solar values. Such a low [alpha/Fe] can be produced by stochastical fluctuations in terms of an incomplete mixing of single Type Ia and II SNe events into the ISM. Our derived element ratios are consistent with the episodic and extended SF in Carina known from its color-magnitude diagram. We find a considerable star-to-star scatter in the abundance ratios. This suggests that Carina's SF history varies with position within the galaxy, with incomplete mixing. Alternatively, the SF rate is so low that the high-mass stellar IMF is sparsely populated, as statistically expected in low-mass star clusters, leading to real scatter in the resultant mass-integrated yields. Both ideas are consistent with slow stochastic SF in dissolving associations, so that one may not speak of a single SF history at a detailed level (Abridged).

We investigate a wide range of possible evolutionary histories for the recently discovered Bootes dwarf spheroidal galaxy, a Milky Way satellite. By means of N-body simulations we follow the evolution of possible progenitor galaxies of Bootes for a variety of orbits in the gravitational potential of the Milky Way. The progenitors considered cover the range from dark-matter-free star clusters to massive, dark-matter dominated outcomes of cosmological simulations. For each type of progenitor and orbit we compare the observable properties of the remnant after 10 Gyr with those of Bootes observed today. Our study suggests that the progenitor of Bootes must have been, and remains now, dark matter dominated. In general our models are unable to reproduce the observed high velocity dispersion in Bootes without dark matter. Our models do not support time-dependent tidal effects as a mechanism able to inflate significantly the internal velocity dispersion. As none of our initially spherical models is able to reproduce the elongation of Bootes, our results suggest that the progenitor of Bootes may have had some intrinsic flattening. Although the focus of the present paper is the Bootes dwarf spheroidal, these models may be of general relevance to understanding the structure, stability and dark matter content of all dwarf spheroidal galaxies.

We present the projected velocity dispersion profile for the remote (d=233kpc) Galactic dwarf spheroidal (dSph) galaxy Leo II, based on 171 discrete stellar radial velocities that were obtained from medium-resolution spectroscopy using the FLAMES/GIRAFFE spectrograph at the European Southern Observatory, Chile. The dispersion profile of those stars with good membership probabilities is essentially flat with an amplitude of 6.6+-0.7 km/s over the full radial extent of our data, which probe to the stellar boundary of this galaxy. We find no evidence of any significant apparent rotation or velocity asymmetry which suggests that tidal effects cannot be invoked to explain Leo II's properties. From basic mass modeling, employing Jeans' equation, we derive a mass out to the limiting radius of (2.7+-0.5) 10^7 Msun and a global mass to light ratio of 27-45 in solar units, depending on the adopted total luminosity. A cored halo profile and a mild amount of tangential velocity anisotropy is found to account well for Leo II's observed kinematics, although we cannot exclude the possibility of a cusped halo with radially varying velocity anisotropy. All in all, this galaxy exhibits dark matter properties which appear to be concordant with the other dSph satellites of the Milky Way, namely a halo mass profile which is consistent with a central core and a total mass which is similar to the common mass scale seen in other dSphs.

We present a synthesis of recent photometric and kinematic data for several of the most dark-matter dominated galaxies. There is a bimodal distribution in half-light radii, with stable star clusters always being smaller than $\sim30$pc, while stable galaxies are always larger than $\sim120$pc. We extend the previously known observational relationships and interpret them in terms of a more fundamental pair of intrinsic properties of dark matter itself: dark matter forms cored mass distributions, with a core scale length of greater than about 100pc, and always has a maximum central massdensity with a narrow range. The dark matter in dSph galaxies appears to be clustered such that there is a mean volume mass density within the stellar distribution which has the very low value of about 0.1$\Msun$ pc$^{-3}$ (about 5GeV/c$^2$ cm$^{-3}$). All dSphs have velocity dispersions equivalent to circular velocities at the edge of their light distributions of $\sim 15$km s$^{-1}$. In two dSphs there is evidence that the density profile is shallow (cored) in the inner regions, and so far none of the dSphs display kinematics which require the presence of an inner cusp. The maximum central dark matter density derived is model dependent, but is likely to have a mean value (averaged over a volume of radius 10pc) of $\sim0.1\Msun$ pc$^{-3}$ (about 5GeV/c$^2$ cm$^{-3}$) for our proposed cored dark mass distributions (where it is similar to the mean value), or $\sim60\Msun$ pc$^{-3}$ (about 2TeV/c$^2$ cm$^{-3}$) if the dark matter density distribution is cusped. Galaxies are embedded in dark matter halos with these properties; smaller systems containing dark matter are not observed.

We present low-resolution spectroscopy of 120 red giants in the Galactic satellite dwarf spheroidal (dSph) Leo I, obtained with the GeminiN-GMOS and Keck-DEIMOS spectrographs. We find stars with velocities consistent with membership of Leo I out to 1.3 King tidal radii. By measuring accurate radial velocities with a median measurement error of 4.6 km/s we find a mean systemic velocity of 284.2 km/s with a global velocity dispersion of 9.9 km/s. The dispersion profile is consistent with being flat out to the last data point. We show that a marginally-significant rise in the radial dispersion profile at a radius of 3' is not associated with any real localized kinematical substructure. Given its large distance from the Galaxy, tides are not likely to have affected the velocity dispersion, a statement we support from a quantitative kinematical analysis, as we observationally reject the occurrence of a significant apparent rotational signal or an asymmetric velocity distribution. Mass determinations adopting both isotropic stellar velocity dispersions and more general models yield a M/L ratio of 24, which is consistent with the presence of a significant dark halo with a mass of about 3x10^7 M_sun, in which the luminous component is embedded. This suggests that Leo I exhibits dark matter properties similar to those of other dSphs in the Local Group. Our data allowed us also to determine metallicities for 58 of the targets. We find a mildly metal poor mean of -1.31 dex and a full spread covering 1 dex. In contrast to the majority of dSphs, Leo I appears to show no radial gradient in its metallicities, which points to a negligible role of external influences in this galaxy's evolution.

We present metallicities and ages for 52 red giants in the remote Galactic dwarf spheroidal (dSph) galaxy Leo II. These stars cover the entire surface area of Leo II and are radial velocity members. We obtained medium-resolution multi-fiber spectroscopy with ESO/VLT's FLAMES spectrograph. The metallicities were determined based on the near-infrared Ca II triplet. The resulting metallicity distribution (MD) is asymmetric and peaks at [Fe/H]=-1.74 dex on the Carretta & Gratton scale. The full range in metallicities extends from -2.4 to -1.1 dex. As in other dSphs, no extremely metal-poor red giants were found. We compare Leo II's observed MD with model predictions for several other Galactic dSphs from the literature. Leo II clearly exhibits a lack of more metal poor stars, in analogy to the classical G-dwarf problem, which may indicate a comparable `K-giant problem'. Moreover, its evolution appears to have been affected by galactic winds. We use our inferred metallicities as an input parameter for isochrone fits to SDSS photometry and derive approximate ages. The resulting age-metallicity distribution covers the full age range from 2-15 Gyr on our adopted isochrone scale. During the first 7 Gyr relative to the oldest stars [Fe/H] appears to have remained almost constant. The almost constant metallicity at higher ages and a slight drop by about 0.3 dex thereafter may be indicative of rejuvenation by low metallicity gas. Overall, the age-metallicity relation appears to support the formation of Leo II from pre-enriched gas. Evidence for enrichment is seen during the recent 2-4 Gyr. Our findings support earlier photometric findings of Leo II as a galaxy with a prominent old and a dominant intermediate-age population. We do not find a significant radial metallicity gradient nor age gradient in our data.(Abridged)

The Milky Way satellite dwarf spheroidal (dSph) galaxies are the smallest dark matter dominated systems in the universe. We have underway dynamical studies of the dSph to quantify the shortest scale lengths on which Dark Matter is distributed, the range of Dark Matter central densities, and the density profile(s) of DM on small scales. Current results suggest some surprises: the central DM density profile is typically cored, not cusped, with scale sizes never less than a few hundred pc; the central densities are typically 10-20 GeV/cc; no galaxy is found with a dark mass halo less massive than ~5.10^7 M_sun. We are discovering many more dSphs, which we are analysing to test the generality of these results.

We present five new satellites of the Milky Way discovered in Sloan Digital Sky Survey (SDSS) imaging data, four of which were followed-up with either the Subaru or the Isaac Newton Telescopes. They include four probable new dwarf galaxies -- one each in the constellations of Coma Berenices, Canes Venatici, Leo and Hercules -- together with one unusually extended globular cluster, Segue 1. We provide distances, absolute magnitudes, half-light radii and color-magnitude diagrams for all five satellites. The morphological features of the color-magnitude diagrams are generally well described by the ridge line of the old, metal-poor globular cluster M92. In the last two years, a total of ten new Milky Way satellites with effective surface brightness mu_v >~ 28 mag/sq. arcsec have been discovered in SDSS data. They are less luminous, more irregular and appear to be more metal-poor than the previously-known nine Milky Way dwarf spheroidals. The relationship between these objects and other populations is discussed. We note that there is a paucity of objects with half-light radii between ~40 pc and ~ 100 pc. We conjecture that this may represent the division between star clusters and dwarf galaxies.

In this Letter, we study a localized stellar overdensity in the constellation of Ursa Major, first identified in Sloan Digital Sky Survey (SDSS) data and subsequently followed up with Subaru imaging. Its color-magnitude diagram (CMD) shows a well-defined sub-giant branch, main sequence and turn-off, from which we estimate a distance of ~30 kpc and a projected size of ~250 x 125 pc. The CMD suggests a composite population with some range in metallicity and/or age. Based on its extent and stellar population, we argue that this is a previously unknown satellite galaxy of the Milky Way, hereby named Ursa Major II (UMa II) after its constellation. Using SDSS data, we find an absolute magnitude of M_V \~ -3.8, which would make it the faintest known satellite galaxy. UMa II's isophotes are irregular and distorted with evidence for multiple concentrations; this suggests that the satellite is in the process of disruption.

We surveyed 1.75 square degrees of sky near Neptune to an R-band 50% detection efficiency of 25.8 mags (corresponding to radii of about 17 km for an assumed albedo of 0.04). We discovered one new outer satellite, Psamathe (S/2003 N1), about 20 km in radius and having a distant retrograde orbit and moderate eccentricity. Until 2003 Neptune was only known to have two satellites which exhibited orbital signatures indicative of capture. Both of these, Triton and Nereid, are unusual when compared to the irregular satellites of other giant planets. With recent discoveries of four additional satellites by Holman et al. (2004) it is now apparent that Neptune has a distant ``normal'' irregular satellite system in which the satellites have radii and orbital properties similar to those of the satellites of other giant planets. We find that the satellite size distribution at Neptune is not well determined given the few objects known to date, being especially sensitive to the inclusion of Triton and Nereid in the sample. Finally, we note that Psamathe and S/2002 N4 have similar semi-major axes, inclinations and eccentricities. They may be fragments of a once larger satellite.

The status of kinematic observations in Local Group dwarf spheroidal galaxies (dSphs) is reviewed. Various approaches to the dynamical modelling of these data are discussed and some general features of dSph dark matter haloes based on simple mass models are presented.

The history of the Milky Way Galaxy is written in the properties of its stellar populations. Here we analyse stars observed as part of surveys of local dwarf spheroidal galaxies, but which from their kinematics are highly probable to be non-members. The selection function -- designed to target metal-poor giants in the dwarf galaxies, at distances of ~100kpc -- includes F-M dwarfs in the Milky Way, at distances of up to several kpc. Thestars whose motions are analysed here lie in the cardinal directions of Galactic longitude l ~ 270 and l ~ 90, where the radial velocity is sensitive to the orbital rotational velocity. We demonstrate that the faint F/G stars contain a significant population with V_phi ~ 100km/s, similar to that found by a targeted, but limited in areal coverage, survey of thick-disk/halo stars by Gilmore, Wyse & Norris (2002). This value of mean orbital rotation does not match either the canonical thick disk or the stellar halo. We argue that this population, detected at both l ~ 270 and l ~ 90, has the expected properties of `satellite debris' in the thick-disk/halo interface, which we interpret as remnants of the merger that heated a pre-existing thin disk to form the thick disk.

There are two main tidal effects which can act on the Local Group dwarf spheroidals (dSphs): tidal stripping and tidal shocking. Using N-body simulations, we show that tidal stripping always leads to flat or rising projected velocity dispersions beyond a critical radius; it is $\sim 5$ times more likely, when averaging over all possible projection angles, that the cylindrically averaged projected dispersion will rise, rather than be flat. In contrast, the Local Group dSphs, as a class, show flat or falling projected velocity dispersions interior to $\sim 1$ kpc. This argues for tidal stripping being unimportant interior to $\sim 1$ kpc for most of the Local Group dSphs observed so far. These results have important implications for the formation of the dSphs and for cosmology. A tidal origin for the formation of these Local Group dSphs (in which they contain no dark matter) is strongly disfavoured. In the cosmological context, a naive solution to the missing satellites problem is to allow only the most massive substructure dark matter halos around the Milky Way to form stars. It is possible for dSphs to reside within these halos ($\sim 10^{10}$M$_\odot$) and have their velocity dispersions lowered through the action of tidal shocks, but only if they have a central density core in their dark matter, rather than a cusp. dSphs can reside within cuspy dark matter halos if their halos are less massive ($\sim 10^{9}$M$_\odot$) and therefore have smaller central velocity dispersions initially (abridged).

We present an improved analytic calculation for the tidal radius of satellites and test our results against N-body simulations. The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and {\it the orbit of the star within the satellite}. We demonstrate that this last point is critical and suggest using {\it three tidal radii} to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically. Over short times ($\simlt 1-2$ Gyrs $\sim 1$ satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.

The Carina dwarf spheroidal (dSph) galaxy is the only galaxy of this type that shows clearly episodic star formation separated by long pauses. Here we present metallicities for 437 radial velocity members of this Galactic satellite. We obtained medium-resolution spectroscopy with the multi-object spectrograph FLAMES at the ESO VLT. Our target red giants cover the entire projected surface area of Carina. Our spectra are centered at the near-infrared Ca triplet, which is a well-established metallicity indicator for old and intermediate-age red giants. The resulting data sample provides the largest collection of spectroscopically derived metallicities for a Local Group dSph to date. Four of our likely radial velocity members of Carina lie outside of this galaxy's nominal tidal radius, supporting earlier claims of the possible existence of such stars beyond the main body of Carina. We find a mean metallicity of [Fe/H]=-1.7 dex. The formal full width at half maximum of the metallicity distribution is 0.92 dex, while the full range of metallicities spans ~-3.0<[Fe/H]<0.0 dex. The metallicity distribution might be indicative of several subpopulations. There appears to be a mild radial gradient such that more metal-rich populations are more centrally concentrated, matching a similar trend for an increasing fraction of intermediate-age stars. This as well as the photometric colors of the more metal-rich red giants suggest that Carina exhibits an age-metallicity relation. Indeed the age-metallicity degeneracy seems to conspire to form a narrow red giant branch despite the considerable spread in metallicity and wide range of ages. The metallicity distribution is not well-matched by a simple closed-box model of chemical evolution, but requires models that take into account also infall and outflows. (Abridged)

We present photometric and kinematic information obtained by measuring 197 planetary nebulae (PNs) discovered in the flattened Fornax elliptical galaxy NGC 1344 (also known as NGC 1340) with an on-band, off-band, grism + on-band filter technique. We build the PN luminosity function (PNLF) and use it to derive a distance modulus m-M=31.4, slightly smaller than, but in good agreement with, the surface brightness fluctuation distance. The PNLF also provides an estimate of the specific PN formation rate: 6x10^-12 PNs per year per solar luminosity. Combining the positional information from the on-band image with PN positions measured on the grism + on-band image, we can measure the radial velocities of 195 PNs, some of them distant more than 3 effective radii from the center of NGC 1344. We complement this data set with stellar kinematics derived from integrated spectra along the major and minor axes, and parallel to the major axis of NGC 1344. The line-of-sight velocity dispersion profile indicates the presence of a dark matter halo around this galaxy.

The recently discovered Ursa Major dwarf spheroidal (dSph) galaxy candidate is about five to eight times less luminous than the faintest previously known dSphs And IX, Draco, and UMi. In this Letter, we present velocity measurements of seven color-magnitude selected UMajor candidate stars. Two of them are apparent non-members based on metallicity and velocity, and the remaining five stars yield a systemic heliocentric velocity of v=-52.45 +/-4.27 km/s and a central line of sight velocity dispersion of 9.3 (+11.7 -1.2) km/s, with 95% confidence that the dispersion is >6.5 km/s. Assuming that UMajor is in dynamical equilibrium, it is clearly dark matter dominated, and cannot be a purely stellar system like a globular cluster. It has an inferred central mass-to-light ratio of M/L~500 and, based on our studies of other dSphs, may possess a much larger total mass to light ratio. UMajor is unexpectedly massive for its low luminosity -- indeed, UMajor appears to be the most dark-matter dominated galaxy yet discovered. The presence of so much dark matter in UMajor immediately suggests that it may be a member of the missing population of low-mass galaxies predicted by the Cold Dark Matter (CDM) paradigm. Given the weak correlation between dSph mass and luminosity, it is entirely likely that a population of dark dwarfs surrounds our Galaxy.

We present the results from a suite of N-body simulations of the tidal stripping of two-component dwarf galaxies comprising some stars and dark matter. We show that recent kinematic data from the local group dwarf spheroidal (dSph) galaxies suggests that dSph galaxies must be sufficiently massive ($10^9 - 10^{10}$M$_\odot$) that tidal stripping is of little importance for the stars. We discuss the implications of these massive dSph galaxies for cosmology and galaxy formation.

We present metallicities for 487 red giants in the Carina dwarf spheroidal (dSph) galaxy that were obtained from FLAMES low-resolution Ca triplet (CaT) spectroscopy. We find a mean [Fe/H] of -1.91 dex with an intrinsic dispersion of 0.25 dex, whereas the full spread in metallicities is at least one dex. The analysis of the radial distribution of metallicities reveals that an excess of metal poor stars resides in a region of larger axis distances. These results can constrain evolutionary models and are discussed in the context of chemical evolution in the Carina dSph.

We present a deep optical survey of Uranus' Hill sphere for small satellites. The Subaru 8-m telescope was used to survey about 3.5 square degrees with a 50% detection efficiency at limiting red magnitude m = 26.1 mag. This magnitude corresponds to objects that are about 7 km in radius (assuming an albedo of 0.04). We detected (without prior knowledge of their positions) all previously known outer satellites and discovered two new irregular satellites (S/2001 U2 and S/2003 U3). The two inner satellites Titania and Oberon were also detected. One of the newly discovered satellites (S/2003 U3) is the first known irregular prograde of the planet. The population, size distribution and orbital parameters of Uranus' irregular satellites are remarkably similar to the irregular satellites of gas giant Jupiter. Both have shallow size distributions (power law indices q~2 for radii > 7 km) with no correlation between the sizes of the satellites and their orbital parameters. However, unlike those of Jupiter, Uranus' irregular satellites do not appear to occupy tight distinct dynamical groups in semi-major axis versus inclination phase space. Two groupings in semi-major axis versus eccentricity phase space appear to be statistically significant.

We surveyed the Hill sphere of Mars for irregular satellites. Our search covered nearly the entire Hill Sphere, but scattered light from Mars excluded the inner few arcminutes where the satellites Phobos and Deimos reside. No new satellites were found to an apparent limiting red magnitude of 23.5, which corresponds to radii of about 0.09 km using an albedo of 0.07.

We present the line of sight radial velocity dispersion profile of the Sextans dwarf spheroidal galaxy (dSph), based on a sample of 88 stars extending to about 1 degree (about 1.5 kpc). Like the Draco and Ursa Minor dSphs, Sextans shows some evidence of a fall-off in the velocity dispersion at large projected radii, with significance p=0.96. Surprisingly, the dispersion at the very centre of Sextans is close to zero (with significance p=0.96). We present evidence which suggests that this latter change in the stellar kinematics coincides with changes in the stellar populations within the dSph. We discuss possible scenarios which could lead to a kinematically and photometrically distinct population at the centre of Sextans.

We present projected velocity dispersion profiles for the Draco and Ursa Minor (UMi) dwarf spheroidal galaxies based on 207 and 162 discrete stellar velocities, respectively. Both profiles show a sharp decline in the velocity dispersion outside ~30 arcmin (Draco) and ~40 arcmin (UMi). New, deep photometry of Draco reveals a break in the light profile at ~25 arcmin. These data imply the existence of a kinematically cold population in the outer parts of both galaxies. Possible explanations of both the photometric and kinematic data in terms of both equilibrium and non-equilibrium models are discussed in detail. We conclude that these data challenge the picture of dSphs as simple, isolated stellar systems.

The nearby Ursa Minor dwarf spheroidal (UMi dSph) is one of the most dark matter dominated galaxies known, with a central mass to light ratio roughly equal to 70. Somewhat anomalously, it appears to contain morphological substructure in the form of a second peak in the stellar number density. It is often argued that this substructure must be transient because it could not survive for the > 10 Gyr age of the system, given the crossing time implied by UMi's 8.8 km/s internal velocity dispersion. In this paper, however, we present evidence that the substructure has a cold kinematical signature, and argue that UMi's clumpiness could indeed be a primordial artefact. Using numerical simulations, we demonstrate that substructure is incompatible with the cusped dark matter haloes predicted by the prevailing Cold Dark Matter (CDM) paradigm, but is consistent with an unbound stellar cluster sloshing back and forth within the nearly harmonic potential of a cored dark matter halo. Thus CDM appears to disagree with observation at the least massive, most dark matter dominated end of the galaxy mass spectrum.

We present the first clear evidence for an extended dark matter halo in the Draco dwarf spheroidal galaxy based on a sample of new radial velocities for 159 giant stars out to large projected radii. Using a two parameter family of halo models spanning a range of density profiles and velocity anisotropies, we are able to rule out (at about the 2.5 sigma confidence level) haloes in which mass follows light. The data strongly favor models in which the dark matter is significantly more extended than the visible dwarf. However, haloes with harmonic cores larger than the light distribution are also excluded. When combined with existing measurements of the proper motion of Draco, our data strongly suggest that Draco has not been tidally truncated within ~1 kpc. We also show that the rising velocity dispersion at large radii represents a serious problem for modified gravity (MOND).

This paper introduces a new two-parameter family of dwarf spheroidal (dSph) galaxy models. The density distribution has a Plummer profile and falls like the inverse fourth power of distance in projection, in agreement with the star-count data. The first free parameter controls the velocity anisotropy, the second controls the dark matter content. The dark matter distribution can be varied from one extreme of mass-follows-light through a near-isothermal halo with flat rotation curve to the other extreme of an extended dark halo with harmonic core. This family of models is explored analytically in some detail -- the distribution functions, the intrinsic moments and the projected moments are all calculated. For the nearby Galactic dSphs, samples of hundreds of discrete radial velocities are becoming available. A technique is developed to extract the anisotropy and dark matter content from such data sets by maximising the likelihood function of the sample of radial velocities. This is constructed from the distribution function and corrected for observational errors and the effects of binaries. Tests on simulated data sets show that samples of 1000 discrete radial velocities are ample to break the degeneracy between mass and anisotropy in the nearby dSphs. Interesting constraints can already be placed on the distribution of the dark matter with samples of 160 radial velocities (the size of the present-day data set for Draco).

We present stellar radial velocity data for the Draco dwarf spheroidal galaxy obtained using the AF2/WYFFOS instrument combination on the William Herschel Telescope. Our dataset consists of 186 member stars, 159 of which have good quality velocities, extending to a magnitude V=19.5 with a mean velocity precision of 2 km/s. We find statistically strong evidence of a rising velocity dispersion consistent with a dark matter halo with gently rising rotation curve. There is a <2 sigma signature of rotation about the long axis, inconsistent with tidal disruption being the source of the rising dispersion. By comparing our dataset with earlier velocities, we find that Draco probably has a binary distribution and fraction comparable to those in the solar neighbourhood. We apply a novel maximum likelihood algorithm and fit the velocity data to a two parameter spherical model with an adjustable dark matter content and velocity anisotropy. Draco is best fit by a weakly tangentially anisotropic distribution of stellar orbits in a dark matter halo with a very slowly rising rotation law. We are able to rule out both a mass-follows-light distribution and an extended halo with a harmonic core at the 2.5 to 3 sigma significance level, depending on the details of our assumptions about Draco's stellar binary population. Our modelling lends support to the idea that the dark matter in dwarf spheroidals is distributed in the form of massive, nearly isothermal haloes.