We present a Bayesian method to identify multiple (chemodynamic) stellar populations in dwarf spheroidal galaxies (dSphs) using velocity, metallicity, and positional stellar data without the assumption of spherical symmetry. We apply this method to a new Keck/DEIMOS spectroscopic survey of the Ursa Minor (UMi) dSph. We identify 892 likely members, making this the largest UMi sample with line-of-sight velocity and metallicity measurements. Our Bayesian method detects two distinct chemodynamic populations with high significance ($\ln{B}\sim33$). The metal-rich ($[{\rm Fe/H}]=-2.05\pm0.03$) population is kinematically colder (radial velocity dispersion of $\sigma_v=4.9\pm0.8 \, {\rm km \, s^{-1}}$) and more centrally concentrated than the metal-poor ($[{\rm Fe/H}]=-2.29\pm0.05$) and kinematically hotter population ($\sigma_v =11.5\pm0.9\, {\rm km \, s^{-1}}$). Furthermore, we apply the same analysis to an independent MMT/Hectochelle data set and confirm the existence of two chemodynamic populations in UMi. In both data sets, the metal-rich population is significantly flattened ($\epsilon=0.75\pm0.03$) and the metal-poor population is closer to spherical ($\epsilon=0.33_{-0.09}^{+0.12}$). Despite the presence of two populations, we are unable to robustly estimate the slope of the dynamical mass profile. We found hints for prolate rotation of order $\sim 2 \, {\rm km \, s^{-1}}$ in the MMT data set, but further observations are required to verify this. The flattened metal-rich population invalidates assumptions built into simple dynamical mass estimators, so we computed new astrophysical dark matter annihilation (J) and decay profiles based on the rounder, hotter metal-poor population and inferred $\log_{10}{(J(0.5^{\circ})/{\rm GeV^{2} \, cm^{-5}})}\approx19.1$ for the Keck data set. Our results paint a more complex picture of the evolution of Ursa Minor than previously discussed.

We study ram pressure stripping in simulated Milky Way-like halos at z>=2 from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. These simulations reach exquisite resolution in their circumgalactic medium (CGM) gas owing to FOGGIE's novel refinement scheme. The CGM of each halo spans a wide dynamic range in density and velocity over its volume---roughly 6 dex and 1000 km/s, respectively---translating into a 5 dex range in ram pressure imparted to interacting satellites. The ram pressure profiles of the simulated CGM are highly stochastic, owing to kpc-scale variations of the density and velocity fields of the CGM gas. As a result, the efficacy of ram pressure stripping depends strongly on the specific path a satellite takes through the CGM. The ram-pressure history of a single satellite is generally unpredictable and not well correlated with its approach vector with respect to the host galaxy. The cumulative impact of ram pressure on the simulated satellites is dominated by only a few short strong impulses---on average, 90% of the total surface momentum gained through ram pressure is imparted in 20% or less of the total orbital time. These results reveal an erratic mode of ram pressure stripping in Milky-Way like halos at high redshift---one that is not captured by a smooth spherically-averaged model of the circumgalactic medium.

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.

We present chemical abundances of red giant branch (RGB) stars in the dwarf spheroidal (dSph) satellite system of Andromeda (M31), using spectral synthesis of medium resolution (R $\sim 6000$) spectra obtained with the Keck II telescope and DEIMOS spectrograph via the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) survey. We coadd stars according to their similarity in photometric metallicity or effective temperature to obtain a signal-to-noise ratio (S/N) high enough to measure average [Fe/H] and [$\alpha$/Fe] abundances. We validate our method using high S/N spectra of RGB stars in Milky Way globular clusters as well as deep observations for a subset of the M31 dSphs in our sample. For this set of validation coadds, we compare the weighted average abundance of the individual stars with the abundance determined from the coadd. We present individual and coadded measurements of [Fe/H] and [$\alpha$/Fe] for stars in ten M31 dSphs, including the first [$\alpha$/Fe] measurements for And IX, XIV, XV, and XVIII. These fainter, less massive dSphs show declining [$\alpha$/Fe] relative to [Fe/H], implying an extended star formation history. In addition, these dSphs also follow the same mass-metallicity relation found in other Local Group satellites. The conclusions we infer from coadded spectra agree with those from previous measurements in brighter M31 dSphs with individual abundance measurements, as well as conclusions from photometric studies. These abundances greatly increase the number of spectroscopic measurements of the chemical composition of M31's less massive dwarf satellites, which are crucial to understanding their star formation history and interaction with the M31 system.

While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

We present new chemo--kinematics of the Hercules dwarf galaxy based on Keck II-- DEIMOS spectroscopy. Our 21 confirmed members have a systemic velocity of $v_{\mathrm{Herc}}=46.4\pm1.3$ kms$^{-1}$ and a velocity dispersion $\sigma_{v,\mathrm{Herc}}=4.4^{+1.4}_{-1.2}$ kms$^{-1}$. From the strength of the Ca II triplet, we obtain a metallicity of [Fe/H]= $-2.48\pm0.19$ dex and dispersion of $\sigma_{\rm{[Fe/H]}}= 0.63^{+0.18}_{-0.13}$ dex. This makes Hercules a particularly metal--poor galaxy, placing it slightly below the standard mass--metallicity relation. Previous photometric and spectroscopic evidence suggests that Hercules is tidally disrupting and may be on a highly radial orbit. From our identified members, we measure no significant velocity gradient. By cross--matching with the second \textit{Gaia} data release, we determine an uncertainty--weighted mean proper motion of $\mu_{\alpha}^*=\mu_{\alpha}\cos(\delta)=-0.153\pm{0.074}$ mas yr$^{-1}$, $\mu_{\delta}=-0.397\pm0.063$ mas yr$^{-1}$. This proper motion is slightly misaligned with the elongation of Hercules, in contrast to models which suggest that any tidal debris should be well aligned with the orbital path. Future observations may resolve this tension.

With a central surface brightness of $\mu_0=29.3$ mag. per sq. arcsec, and half-light radius of $r_{\rm half}=3.1^{+0.9}_{-1.1}$~kpc, Andromeda XIX (And XIX) is an extremely diffuse satellite of Andromeda. We present spectra for $\sim100$ red giant branch stars in this galaxy, plus 16 stars in a nearby stellar stream. With this exquisite dataset, we re-derive the properties of And XIX, measuring a systemic velocity of $v_r=-109.0\pm1.6$ km/s and a velocity dispersion of $\sigma_v = 7.8^{+1.7}_{-1.5}$ km/s (higher than derived in our previous work). We marginally detect a velocity gradient along the major axis of ${\rm d}v/{\rm d}r = -2.1\pm1.8$ km/s kpc$^{-1}$. We find its mass-to-light ratio is higher than galaxies of comparable stellar mass ($[M/L]_{\rm half} = 278^{+146}_{-198}M_\odot/L_\odot$), but its dynamics place it in a halo with a similar total mass to these galaxies. This could suggest that And XIX is a "puffed up" dwarf galaxy, whose properties have been altered by tidal processes, similar to its Milky Way counterpart, Antlia II. For the nearby stream, we measure $v_r=-279.2\pm3.7$ km/s, and $\sigma_v=13.8^{+3.5}_{-2.6}$ km/s. We measure its metallicity, and find it to be more metal rich than And XIX, implying that the two features are unrelated. Finally, And XIX's dynamical and structural properties imply it is a local analogue to ultra diffuse galaxies (UDGs). Its complex dynamics suggest that the masses of distant UDGs measured from velocity dispersions alone should be carefully interpreted

A Kavli foundation sponsored workshop on the theme \emph{Petabytes to Science} was held 12$^{th}$ to 14$^{th}$ of February 2019 in Las Vegas. The aim of the this workshop was to discuss important trends and technologies which may support astronomy. We also tackled how to better shape the workforce for the new trends and how we should approach education and public outreach. This document was coauthored during the workshop and edited in the weeks after. It comprises the discussions and highlights many recommendations which came out of the workshop. We shall distill parts of this document and formulate potential white papers for the decadal survey.

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

Fundamental coding and software development skills are increasingly necessary for success in nearly every aspect of astronomical and astrophysical research as large surveys and high resolution simulations become the norm. However, professional training in these skills is inaccessible or impractical for many members of our community. Students and professionals alike have been expected to acquire these skills on their own, apart from formal classroom curriculum or on-the-job training. Despite the recognized importance of these skills, there is little opportunity to develop them - even for interested researchers. To ensure a workforce capable of taking advantage of the computational resources and the large volumes of data coming in the next decade, we must identify and support ways to make software development training widely accessible to community members, regardless of affiliation or career level. To develop and sustain a technology capable astronomical and astrophysical workforce, we recommend that agencies make funding and other resources available in order to encourage, support and, in some cases, require progress on necessary training, infrastructure and policies. In this white paper, we focus on recommendations for how funding agencies can lead in the promotion of activities to support the astronomy and astrophysical workforce in the 2020s.

We present a kinematic and spectroscopic analysis of 38 red giant branch stars, in 7 fields, spanning the dwarf spheroidal galaxy Andromeda XXVII and the upper segment of the North West Stream. Both features are located in the outer halo of the Andromeda galaxy at a projected radius of 50-80 kpc, with the stream extending for $\sim$3$^{\circ}$ on the sky. Our data is obtained as part of the PAndAS survey and enables us to confirm that Andromeda XXVII's heliocentric distance is 827 $\pm$ 47 kpc and spectroscopic metallicity is -2.1$^{+0.4}_{-0.5}$. We also re-derive Andromeda XXVII's kinematic properties, measuring a systemic velocity = -526.1$^{+10.0}_{-11.0}${\kms} and a velocity dispersion that we find to be non-Gaussian but for which we derive a formal value of 27.0$^{+2.2}_{-3.9}${\kms}. In the upper segment of the North West Stream we measure mean values for the metallicity = -1.8$\pm$0.4, systemic velocity = -519.4 $\pm$4.0{\kms} and velocity dispersion = 10.0$\pm$4.0{\kms}. We also detect a velocity gradient of 1.7$\pm$0.3 {\kms} kpc$^{-1}$ on an infall trajectory towards M31. With a similar gradient, acting in the same direction, in the lower segment we suggest that the North West Stream is not a single structure. As the properties of the upper segment of the North West Stream and Andromeda XXVII are consistent within ~90\% confidence limits, it is likely that the two are related and plausible that Andromeda XXVII is the progenitor of this stream.

Software is a critical part of modern research, and yet there are insufficient mechanisms in the scholarly ecosystem to acknowledge, cite, and measure the impact of research software. The majority of academic fields rely on a one-dimensional credit model whereby academic articles (and their associated citations) are the dominant factor in the success of a researcher's career. In the petabyte era of astronomical science, citing software and measuring its impact enables academia to retain and reward researchers that make significant software contributions. These highly skilled researchers must be retained to maximize the scientific return from petabyte-scale datasets. Evolving beyond the one-dimensional credit model requires overcoming several key challenges, including the current scholarly ecosystem and scientific culture issues. This white paper will present these challenges and suggest practical solutions for elevating the role of software as a product of the research enterprise.

In recent years, ultrafaint dwarf (UFD) galaxies have been found through systematic searches of large optical surveys. However, the existence of Leo T, a nearby gas-rich dwarf, suggests that there could be other nearby UFDs that are optically obscured but have gas detectable at nonoptical wavelengths. With this in mind, we perform a search of the full Galactic Arecibo $L-$band Feed Array HI (GALFA-HI) survey, a radio survey which covers one-third of the sky at velocities $-650 < V_{\rm{LSR}} < +650 \; \rm{km} \; \rm{s^{-1}}$, for neutral hydrogen sources. We are able to probe regions of the sky at lower Galactic latitudes and smaller $|V_{LSR}|$ compared to previous explorations. We use the Source Finding Application (SoFiA) on GALFA-HI and select all sources with similar properties to Leo T and other local dwarf galaxies. We find 690 dwarf galaxy candidates, one of which is particularly promising and likely a new galaxy near the Galactic plane ($b=-8^{\circ}$) that is comparable in velocity width and HI-flux to other recently discovered local volume galaxies. We find we are sensitive to Leo T-like objects out to $1 \; \rm{Mpc}$ at velocities clear from background HI emission. We check each candidate's corresponding optical fields from Pan-STARRS and fit stars drawn from isochrones, but find no evidence of stellar populations. We thus find no other Leo T-like dwarfs within $500 \; \rm{kpc}$ of the Milky Way in the one-third of the sky covered by the GALFA-HI footprint and discuss our nondetection in a cosmological context.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

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

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

astroquery is a collection of tools for requesting data from databases hosted on remote servers with interfaces exposed on the internet, including those with web pages but without formal application program interfaces (APIs). These tools are built on the Python requests package, which is used to make HTTP requests, and astropy, which provides most of the data parsing functionality. astroquery modules generally attempt to replicate the web page interface provided by a given service as closely as possible, making the transition from browser-based to command-line interaction easy. astroquery has received significant contributions from throughout the astronomical community, including several significant contributions from telescope archives. astroquery enables the creation of fully reproducible workflows from data acquisition through publication. This paper describes the philosophy, basic structure, and development model of the astroquery package. The complete documentation for astroquery can be found at http://astroquery.readthedocs.io/.

We present models of the H$\beta$-emitting broad-line region (BLR) in seven Seyfert 1 galaxies from the Lick AGN (Active Galactic Nucleus) Monitoring Project 2011 sample, drawing inferences on the BLR structure and dynamics as well as the mass of the central supermassive black hole. We find that the BLR is generally a thick disk, viewed close to face-on, with preferential emission back toward the ionizing source. The dynamics in our sample range from near-circular elliptical orbits to inflowing or outflowing trajectories. We measure black hole masses of $\log_{10}(M_{\rm BH}/M_\odot) = 6.48^{+0.21}_{-0.18}$ for PG 1310$-$108, $7.50^{+0.25}_{-0.18}$ for Mrk 50, $7.46^{+0.15}_{-0.21}$ for Mrk 141, $7.58^{+0.08}_{-0.08}$ for Mrk 279, $7.11^{+0.20}_{-0.17}$ for Mrk 1511, $6.65^{+0.27}_{-0.15}$ for NGC 4593, and $6.94^{+0.14}_{-0.14}$ for Zw 229$-$015. We use these black hole mass measurements along with cross-correlation time lags and line widths to recover the scale factor $f$ used in traditional reverberation mapping measurements. Combining our results with other studies that use this modeling technique, bringing our sample size to 16, we calculate a scale factor that can be used for measuring black hole masses in other reverberation mapping campaigns. When using the root-mean-square (rms) spectrum and using the line dispersion to measure the line width, we find $\log_{10}(f_{{\rm rms},\sigma})_{\rm pred} = 0.57 \pm 0.19$. Finally, we search for correlations between $f$ and other AGN and BLR parameters and find marginal evidence that $f$ is correlated with $M_{\rm BH}$ and the BLR inclination angle, but no significant evidence of a correlation with the AGN luminosity or Eddington ratio.

Triangulum (M33) is the most massive satellite galaxy of Andromeda (M31), with a stellar mass of about $3\times10^9\; M_{\odot}$. Based on abundance matching techniques, M33's total mass at infall is estimated to be of order $10^{11}\; M_{\odot}$. $\Lambda$CDM simulations predict that M33-mass halos host several of their own low mass satellite companions, yet only one candidate M33 satellite has been detected in deep photometric surveys to date. This `satellites of satellites' hierarchy has recently been explored in the context of the dwarf galaxies discovered around the Milky Way's Magellanic Clouds in the Dark Energy Survey. Here we discuss the number of satellite galaxies predicted to reside within the virial radius ($\sim$160 kpc) of M33 based on $\Lambda$CDM simulations. We also calculate the expected number of satellite detections in N fields of data using various ground--based imagers. Finally, we discuss how these satellite population predictions change as a function of M33's past orbital history. If M33 is on its first infall into M31's halo, its proposed satellites are expected to remain bound to M33 today. However, if M33 experienced a recent tidal interaction with M31, the number of surviving satellites depends strongly on the distance achieved at pericenter due to the effects of tidal stripping. We conclude that a survey extending to $\sim$100 kpc around M33 would be sufficient to constrain its orbital history and a majority of its satellite population. In the era of WFIRST, surveys of this size will be considered small observing programs.

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

It was recently proposed that the globular cluster system of the very low surface-brightness galaxy NGC1052-DF2 is dynamically very cold, leading to the conclusion that this dwarf galaxy has little or no dark matter. Here, we show that a robust statistical measure of the velocity dispersion of the tracer globular clusters implies a mundane velocity dispersion and a poorly constrained mass-to-light ratio. Models that include the possibility that some of the tracers are field contaminants do not yield a more constraining inference. We derive only a weak constraint on the mass-to-light ratio of the system within the half-light radius (M/L_V<6.7 at the 90% confidence level) or within the radius of the furthest tracer (M/L_V<8.1 at the 90% confidence level). This limit may imply a mass-to-light ratio on the low end for a dwarf galaxy but many Local Group dwarf galaxies fall well within this contraint. With this study, we emphasize the need to reliably account for measurement uncertainties and to stay as close as possible to the data when determining dynamical masses from very small data sets of tracers.

The Seyfert 1 galaxy Arp 151 was monitored as part of three reverberation mapping campaigns spanning $2008-2015$. We present modeling of these velocity-resolved reverberation mapping datasets using a geometric and dynamical model for the broad line region (BLR). By modeling each of the three datasets independently, we infer the evolution of the BLR structure in Arp 151 over a total of seven years and constrain the systematic uncertainties in non-varying parameters such as the black hole mass. We find that the BLR geometry of a thick disk viewed close to face-on is stable over this time, although the size of the BLR grows by a factor of $\sim 2$. The dynamics of the BLR are dominated by inflow and the inferred black hole mass is consistent for the three datasets, despite the increase in BLR size. Combining the inference for the three datasets yields a black hole mass and statistical uncertainty of $\log_{10}($M$_{\rm BH}/\rm{M}_{\odot})=6.82^{+0.09}_{-0.09}$ with a standard deviation in individual measurements of 0.13 dex.

We present an approach for comparing the detections and non-detections of Local Group (LG) dwarf galaxies in large HI surveys to the predictions of a suite of n-body simulations of the LG. This approach depends primarily on a set of empirical scaling relations to connect the simulations to the observations, rather than making strong theoretical assumptions. We then apply this methodology to the Galactic Arecibo L-band Feed Array HI Compact Cloud Catalog (GALFA-HI CCC), and compare it to the ELVIS suite of simulations. This approach reveals a strong tension between the na\"ive results of the model and the observations: while there are no LG dwarfs in the GALFA-HI CCC, the simulations predict $\sim 10$. Applying a simple model of reionization can resolve this tension by preventing low-mass halos from forming gas. However, and if this effect operates as expected, the observations provide a constraint on the mass scale of dwarf galaxy that reionization impacts. Combined with the observed properties of Leo T, the halo virial mass scale at which reionization impacts dwarf galaxy gas content is constrained to be $\sim 10^{8.5} M_\odot$, independent of any assumptions about star formation.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

We present astroplan - an open source, open development, Astropy affiliated package for ground-based observation planning and scheduling in Python. astroplan is designed to provide efficient access to common observational quantities such as celestial rise, set, and meridian transit times and simple transformations from sky coordinates to altitude-azimuth coordinates without requiring a detailed understanding of astropy's implementation of coordinate systems. astroplan provides convenience functions to generate common observational plots such as airmass and parallactic angle as a function of time, along with basic sky (finder) charts. Users can determine whether or not a target is observable given a variety of observing constraints, such as airmass limits, time ranges, Moon illumination/separation ranges, and more. A selection of observation schedulers are included which divide observing time among a list of targets, given observing constraints on those targets. Contributions to the source code from the community are welcome.

We present the velocity dispersion of red giant branch (RGB) stars in M31's halo, derived by modeling the line of sight velocity distribution of over 5000 stars in 50 fields spread throughout M31's stellar halo. The dataset was obtained as part of the SPLASH (Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo) Survey, and covers projected radii of 9 to 175 kpc from M31's center. All major structural components along the line of sight in both the Milky Way (MW) and M31 are incorporated in a Gaussian Mixture Model, including all previously identified M31 tidal debris features in the observed fields. The probability an individual star is a constituent of M31 or the MW, based on a set of empirical photometric and spectroscopic diagnostics, is included as a prior probability in the mixture model. The velocity dispersion of stars in M31's halo is found to decrease only mildly with projected radius, from 108 km/s in the innermost radial bin (8.2 to 14.1 kpc) to $\sim 80$ to 90 km/s at projected radii of $\sim 40$ to 130 kpc, and can be parameterized with a power-law of slope $-0.12\pm 0.05$. The quoted uncertainty on the power-law slope reflects only the precision of the method, although other sources of uncertainty we consider contribute negligibly to the overall error budget.

We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python.

We present the survey strategy and early results of the "Satellites Around Galactic Analogs" (SAGA) Survey. The SAGA Survey's goal is to measure the distribution of satellite galaxies around 100 systems analogous to the Milky Way down to the luminosity of the Leo I dwarf galaxy ($ M_r < -12.3 $). We define a Milky Way analog based on $K$-band luminosity and local environment. Here, we present satellite luminosity functions for 8 Milky Way analog galaxies between 20 to 40 Mpc. These systems have nearly complete spectroscopic coverage of candidate satellites within the projected host virial radius down to $ r_o < 20.75 $ using low redshift $gri$ color criteria. We have discovered a total of 25 new satellite galaxies: 14 new satellite galaxies meet our formal criteria around our complete host systems, plus 11 additional satellites in either incompletely surveyed hosts or below our formal magnitude limit. Combined with 13 previously known satellites, there are a total of 27 satellites around 8 complete Milky Way analog hosts. We find a wide distribution in the number of satellites per host, from 1 to 9, in the luminosity range for which there are five Milky Way satellites. Standard abundance matching extrapolated from higher luminosities predicts less scatter between hosts and a steeper luminosity function slope than observed. We find that the majority of satellites (26 of 27) are star-forming. These early results indicate that the Milky Way has a different satellite population than typical in our sample, potentially changing the physical interpretation of measurements based only on the Milky Way's satellite galaxies.

We present a chemodynamical analysis of the Leo~V dwarf galaxy, based on Keck II DEIMOS spectra of 8 member stars. We find a systemic velocity for the system of $\langle v_r\rangle = 170.9^{+ 2.1}_{-1.9}$kms$^{-1}$, and barely resolve a velocity dispersion for the system, with $\sigma_{vr} = 2.3^{+3.2}_{-1.6}$kms$^{-1}$, consistent with previous studies of Leo~V. The poorly resolved dispersion means we are unable to adequately constrain the dark matter content of Leo~V. We find an average metallicity for the dwarf of [Fe/H]$ = -2.48\pm0.21$, and measure a significant spread in the iron abundance of its member stars, with $-3.1\le$[Fe/H]$\le-1.9$ dex, which cleanly identifies Leo~V as a dwarf galaxy that has been able to self-enrich its stellar population through extended star formation. Owing to the tentative photometric evidence for tidal substructure around Leo~V, we also investigate whether there is any evidence for tidal stripping or shocking of the system within its dynamics. We measure a significant velocity gradient across the system, of $\frac{{\rm d}v}{{\rm d}\chi} = -4.1^{+2.8}_{-2.6}$kms$^{-1}$ per arcmin (or $\frac{{\rm d}v}{{\rm d}\chi} = -71.9^{+50.8}_{-45.6}$kms$^{-1}$~kpc$^{-1}$), which points almost directly toward the Galactic centre. We argue that Leo~V is likely a dwarf on the brink of dissolution, having just barely survived a past encounter with the centre of the Milky Way.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

We spectroscopically identify a sample of carbon stars in the satellites and halo of M31 using moderate-resolution optical spectroscopy from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo survey. We present the photometric properties of our sample of 41 stars, including their brightness with respect to the tip of the red giant branch (TRGB) and their distributions in various color-color spaces. This analysis reveals a bluer population of carbon stars fainter than the TRGB and a redder population of carbon stars brighter than the TRGB. We then apply principal component analysis to determine the sample's eigenspectra and eigencoefficients. Correlating the eigencoefficients with various observable properties reveals the spectral features that trace effective temperature and metallicity. Putting the spectroscopic and photometric information together, we find the carbon stars in the satellites and halo of M31 to be minimally impacted by dust and internal dynamics. We also find that while there is evidence to suggest that the sub-TRGB stars are extrinsic in origin, it is also possible that they are are particularly faint members of the asymptotic giant branch.

The ultra-faint satellite galaxy Hercules has a strongly elongated and irregular morphology with detections of tidal features up to 1.3 deg (3 kpc) from its center. This suggests that Hercules may be dissolving under the Milky Way's gravitational influence, and hence could be a tidal stream in formation rather than a bound, dark-matter dominated satellite. Using Bayesian inference in combination with N-body simulations, we show that Hercules has to be on a very eccentric orbit (epsilon~0.95) within the Milky Way in this scenario. On such an orbit, Hercules "explodes" as a consequence of the last tidal shock at pericenter 0.5 Gyr ago. It is currently decelerating towards apocenter of its orbit with a velocity of V=157 km/s -- of which 99% is directed radially outwards. Due to differential orbital precession caused by the non-spherical nature of the Galactic potential, its debris fans out nearly perpendicular to its orbit. This explains why Hercules has an elongated shape without showing a distance gradient along its main body: it is in fact a stream that is significantly broader than it is long. In other words, it is moving perpendicular to its apparent major axis. In this scenario, there is a spike in the radial velocity profile created by the dominant debris component that formed through the last pericenter passage. This is similar to kinematic substructure that is observed in the real Hercules. Modeling a satellite on such a highly eccentric orbit is strongly dependent on the form of the Galactic potential. We therefore propose that detailed kinematic investigation of Hercules and other exploding satellite candidates can yield strong constraints on the potential of the Milky Way.

We present observations of the Pisces A and B galaxies with the Advanced Camera for Surveys on the Hubble Space Telescope. Photometry from these images clearly resolve a Red Giant Branch for both objects, demonstrating that they are nearby dwarf galaxies. We describe a Bayesian inferential approach to determining the distance to these galaxies using the magnitude of the tip of the RGB, and then apply this approach to these galaxies. We also provide the full probability distributions for parameters derived using this approach. This reveals the distance to these galaxies as $5.64^{+0.13}_{-0.15} \, {\rm Mpc}$ and $8.89^{+0.75}_{-0.85} \, {\rm Mpc}$ for Pisces A and B, respectively, placing both within the Local Volume but not the Local Group. We estimate the star formation histories of these galaxies, which suggests that they have recently undergone an increase in their star formation rates. Together these yield luminosities for Pisces A and B of $M_V=-11.57^{+0.06}_{-0.05}$ and $-12.9 \pm 0.2$, respectively, and estimated stellar masses of $\log(M_*/M_{\odot})= 7.0^{+0.4}_{-1.7}$ and $7.5^{+0.3}_{-1.8}$. We further show that these galaxies are likely at the boundary between nearby voids and higher-density filamentary structure. This suggests that they are entering a higher-density region from voids, where they would have experienced delayed evolution, consistent with their recent increased star formation rates. If this is indeed the case, they are useful for study as proxies of the galaxies that later evolved into typical Local Group satellite galaxies.

Context. Pisces A & Pisces B are the only two galaxies found via optical imaging and spectroscopy out of 22 HI clouds identified in the GALFAHI survey as dwarf galaxy candidates. Aims. Derive the HI content and kinematics of Pisces A & B. Methods. Aperture synthesis HI observations using the seven dish Karoo Array Telescope (KAT-7), which is a pathfinder instrument for MeerKAT, the South African precursor to the mid-frequency Square Kilometre Array (SKA-MID). Results. The small rotation velocities detected of ~5 km/sec and ~10 km/sec in Pisces A & B respectively, and their HI content show that they are really dwarf irregular galaxies (dIrr). Despite that small rotation component, it is more the random motions ~9-11 km/sec that provide most of the gravitational support, especially in the outer parts. The study of their kinematics, especially the strong gradients of random motions, suggest that those two dwarf galaxies are not yet in equilibrium. Conclusions. These HI rich galaxies may be indicative of a large population of dwarfs at the limit of detectability. However, such gas-rich dwarf galaxies will most likely never be within the virial radius of MW-type galaxies and become sub-halo candidates. Systems such as Pisces A & B are more likely to be found at a few Mpc.s from MW-type galaxies.

The connection between galaxies and dark matter halos is often inferred from data using probabilistic models, such as the Halo Occupation Distribution (HOD). Conventional HOD formulations assume that only halo mass governs the galaxy-halo connection. Violations of this assumption, known as galaxy assembly bias, threaten the HOD program. We introduce decorated HODs, a new, flexible class of models designed to account for assembly bias. Decorated HODs minimally expand the parameter space and maximize the independence between traditional and novel HOD parameters. We use decorated HODs to quantify the influence of assembly bias on clustering and lensing statistics. For SDSS-like samples, the impact of assembly bias on galaxy clustering can be as large as a factor of two on r ~ 200 kpc scales and ~15% in the linear regime. Assembly bias can either enhance or diminish clustering on large scales, but generally increases clustering on scales r <~ 1 Mpc. We performed our calculations with Halotools, an open-source, community-driven python package for studying the galaxy-halo connection (http://halotools.readthedocs.org). We conclude by describing the use of decorated HODs to treat assembly bias in otherwise conventional likelihood analyses.

In the Local Group, nearly all of the dwarf galaxies (M_star < 10^9 M_sun) that are satellites within 300 kpc (the virial radius) of the Milky Way (MW) and Andromeda (M31) have quiescent star formation and little-to-no cold gas. This contrasts strongly with comparatively isolated dwarf galaxies, which are almost all actively star-forming and gas-rich. This near dichotomy implies a rapid transformation of satellite dwarf galaxies after falling into the halos of the MW or M31. We combine the observed quiescent fractions for satellites of the MW and M31 with the infall times of satellites from the Exploring the Local Volume in Simulations (ELVIS) suite of cosmological zoom-in simulations to determine the typical timescales over which environmental processes within the MW/M31 halos remove gas and quench star formation in low-mass satellite galaxies. The quenching timescales for satellites with M_star < 10^8 M_sun are short, < 2 Gyr, and quenching is more rapid at lower M_star. These satellite quenching timescales can be 1 - 2 Gyr longer if one includes the time that satellites were environmentally preprocessed by low-mass groups prior to MW/M31 infall. We compare with quenching timescales for more massive satellites from previous works to synthesize the nature of satellite galaxy quenching across the observable range of M_star = 10^{3-11} M_sun. The satellite quenching timescale increases rapidly with satellite M_star, peaking at ~9.5 Gyr for M_star ~ 10^9 M_sun, and the timescale rapidly decreases at higher M_star to < 5 Gyr at M_star > 5 x 10^9 M_sun. Overall, galaxies with M_star ~ 10^9 M_sun, similar to the Magellanic Clouds, exhibit the longest quenching timescales, regardless of environmental or internal mechanisms.

We report on an informal survey about the use of software in the worldwide astronomical community. The survey was carried out between December 2014 and February 2015, collecting responses from 1142 astronomers, spanning all career levels. We find that all participants use software in their research. The vast majority of participants, 90%, write at least some of their own software. Even though writing software is so wide-spread among the survey participants, only 8% of them report that they have received substantial training in software development. Another 49% of the participants have received "little" training. The remaining 43% have received no training. We also find that astronomers' software stack is fairly narrow. The 10 most popular tools among astronomers are (from most to least popular): Python, shell scripting, IDL, C/C++, Fortran, IRAF, spreadsheets, HTML/CSS, SQL and Supermongo. Across all participants the most common programing language is Python ($67\pm 2\%$), followed by IDL ($44\pm 2\%$), C/C++ ($37\pm 2\%$) and Fortran ($28\pm 2\%$). IRAF is used frequently by $24\pm 1\%$ of participants. We show that all trends are largely independent of career stage, area of research and geographic location.

We report five Local Volume dwarf galaxies (two of which are presented here for the first time) uncovered during a comprehensive archival search for optical counterparts to ultra-compact high velocity clouds (UCHVCs). The UCHVC population of HI clouds are thought to be candidate gas-rich, low mass halos at the edge of the Local Group and beyond, but no comprehensive search for stellar counterparts to these systems has been presented. Careful visual inspection of all publicly available optical and ultraviolet imaging at the position of the UCHVCs revealed six blue, diffuse counterparts with a morphology consistent with a faint dwarf galaxy beyond the Local Group. Optical spectroscopy of all six candidate dwarf counterparts show that five have an H$\alpha$-derived velocity consistent with the coincident HI cloud, confirming their association, the sixth diffuse counterpart is likely a background object. The size and luminosity of the UCHVC dwarfs is consistent with other known Local Volume dwarf irregular galaxies. The gas fraction ($M_{HI}/M_{star}$) of the five dwarfs are generally consistent with that of dwarf irregular galaxies in the Local Volume, although ALFALFA-Dw1 (associated with ALFALFA UCHVC HVC274.68+74.70$-$123) has a very high $M_{HI}/M_{star}$$\sim$40. Despite the heterogenous nature of our search, we demonstrate that the current dwarf companions to UCHVCs are at the edge of detectability due to their low surface brightness, and that deeper searches are likely to find more stellar systems. If more sensitive searches do not reveal further stellar counterparts to UCHVCs, then the dearth of such systems around the Local Group may be in conflict with $\Lambda$CDM simulations.

We present the results of a deep study of the isolated dwarf galaxies Andromeda XXVIII and Andromeda XXIX with Gemini/GMOS and Keck/DEIMOS. Both galaxies are shown to host old, metal-poor stellar populations with no detectable recent star formation, conclusively identifying both of them as dwarf spheroidal galaxies (dSphs). And XXVIII exhibits a complex horizontal branch morphology, which is suggestive of metallicity enrichment and thus an extended period of star formation in the past. Decomposing the horizontal branch into blue (metal poor, assumed to be older) and red (relatively more metal rich, assumed to be younger) populations shows that the metal rich are also more spatially concentrated in the center of the galaxy. We use spectroscopic measurements of the Calcium triplet, combined with the improved precision of the Gemini photometry, to measure the metallicity of the galaxies, confirming the metallicity spread and showing that they both lie on the luminosity-metallicity relation for dwarf satellites. Taken together, the galaxies exhibit largely typical properties for dSphs despite their significant distances from M31. These dwarfs thus place particularly significant constraints on models of dSph formation involving environmental processes such as tidal or ram pressure stripping. Such models must be able to completely transform the two galaxies into dSphs in no more than two pericentric passages around M31, while maintaining a significant stellar populations gradient. Reproducing these features is a prime requirement for models of dSph formation to demonstrate not just the plausibility of environmental transformation but the capability of accurately recreating real dSphs.

In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the observations, reductions and measurements, and data products from the spectroscopic campaign. The reduced spectra were fitted with a multicomponent model in order to isolate the contributions of various continuum and emission-line components. We present light curves of broad emission lines and the AGN continuum, and measurements of the broad H-beta line widths in mean and root-mean square (rms) spectra. For the most highly variable AGNs we also measured broad H-beta line widths and velocity centroids from the nightly spectra. In four AGNs exhibiting the highest variability amplitudes, we detect anticorrelations between broad H-beta width and luminosity, demonstrating that the broad-line region "breathes" on short timescales of days to weeks in response to continuum variations. We also find that broad H-beta velocity centroids can undergo substantial changes in response to continuum variations; in NGC 4593 the broad H-beta velocity shifted by ~250 km/s over a one-month duration. This reverberation-induced velocity shift effect is likely to contribute a significant source of confusion noise to binary black hole searches that use multi-epoch quasar spectroscopy to detect binary orbital motion. We also present results from simulations that examine biases that can occur in measurement of broad-line widths from rms spectra due to the contributions of continuum variations and photon-counting noise.

We report the discovery of two dwarf galaxies, Pisces A and B, from a blind 21 cm HI search. These were the only two galaxies found via optical imaging and spectroscopy of 22 HI clouds identified in the GALFA-HI survey as dwarf galaxy candidates. They have properties consistent with being in the Local Volume ($<10$ Mpc), and one has resolved stellar populations such that it may be on the outer edge of the Local Group ($\sim 1 \, {\rm Mpc}$ from M31). While the distance uncertainty makes interpretation ambiguous, these may be among the faintest starforming galaxies known. Additionally, rough estimates comparing these galaxies to $\Lambda$CDM dark matter simulations suggest consistency in number density, implying that dark matter halos likely to host these galaxies are primarily HI-rich. The galaxies may thus be indicative of a large population of dwarfs at the limit of detectability that are comparable to the faint satellites of the Local Group. Because they are outside the influence of a large dark matter halo to alter their evolution, these galaxies can provide critical anchors to dwarf galaxy formation models.

We present alpha element to iron abundance ratios, [$\alpha$/Fe], for four stars in the outer stellar halo of the Andromeda Galaxy (M 31). The stars were identified as high-likelihood field halo stars by Gilbert et al. (2012) and lie at projected distances between 70 and 140 kpc from M 31's center. These are the first alpha abundances measured for a halo star in a galaxy beyond the Milky Way. The stars range in metallicity between [Fe/H]= -2.2 and [Fe/H]= -1.4. The sample's average [$\alpha$/Fe] ratio is +0.20+/-0.20. The best-fit average value is elevated above solar which is consistent with rapid chemical enrichment from Type II supernovae. The mean [$\alpha$/Fe] ratio of our M31 outer halo sample agrees (within the uncertainties) with that of Milky Way inner/outer halo stars that have a comparable range of [Fe/H].

We present the metallicity distribution of red giant branch (RGB) stars in M31's stellar halo, derived from photometric metallicity estimates for over 1500 spectroscopically confirmed RGB halo stars. The stellar sample comes from 38 halo fields observed with the Keck/DEIMOS spectrograph, ranging from 9 to 175 kpc in projected distance from M31's center, and includes 52 confirmed M31 halo stars beyond 100 kpc. While a wide range of metallicities is seen throughout the halo, the metal-rich peak of the metallicity distribution function becomes significantly less prominent with increasing radius. The metallicity profile of M31's stellar halo shows a continuous gradient from 9 to ~100 kpc, with a magnitude of -0.01 dex/kpc. The stellar velocity distributions in each field are used to identify stars that are likely associated with tidal debris features. The removal of tidal debris features does not significantly alter the metallicity gradient in M31's halo: a gradient is maintained in fields spanning 10 to 90 kpc. We analyze the halo metallicity profile, as well as the relative metallicities of stars associated with tidal debris features and the underlying halo population, in the context of current simulations of stellar halo formation. We argue that the large scale gradient in M31's halo implies M31 accreted at least one relatively massive progenitor in the past, while the field to field variation seen in the metallicity profile indicates that multiple smaller progenitors are likely to have contributed substantially to M31's outer halo.

Using the Sloan Digital Sky Survey, we examine the quenching of satellite galaxies around isolated Milky Way-like hosts in the local Universe. We find that the efficiency of satellite quenching around isolated galaxies is low and roughly constant over two orders of magnitude in satellite stellar mass ($M_{*}$ = $10^{8.5}-10^{10.5} \, M_{\odot}$), with only $\sim~20\%$ of systems quenched as a result of environmental processes. While largely independent of satellite stellar mass, satellite quenching does exhibit clear dependence on the properties of the host. We show that satellites of passive hosts are substantially more likely to be quenched than those of star-forming hosts, and we present evidence that more massive halos quench their satellites more efficiently. These results extend trends seen previously in more massive host halos and for higher satellite masses. Taken together, it appears that galaxies with stellar masses larger than about $10^{8}~M_{\odot}$ are uniformly resistant to environmental quenching, with the relative harshness of the host environment likely serving as the primary driver of satellite quenching. At lower stellar masses ($< 10^{8}~M_{\odot}$), however, observations of the Local Group suggest that the vast majority of satellite galaxies are quenched, potentially pointing towards a characteristic satellite mass scale below which quenching efficiency increases dramatically.

We present alpha to iron abundance ratios for 226 individual red giant branch stars in nine dwarf galaxies of the Andromeda (M 31) satellite system. The abundances are measured from the combined signal of Mg, Si, Ca, and Ti lines in Keck/DEIMOS medium-resolution spectra. This constitutes the first large sample of alpha abundance ratios measured in the M 31 satellite system. The dwarf galaxies in our sample exhibit a variety of alpha abundance ratios, with the average values in each galaxy ranging from approximately solar ([alpha/Fe] ~ +0.0) to alpha-enhanced ([alpha/Fe] ~ +0.5). These variations do not show a correlation with internal kinematics, environment, or stellar density. We confirm radial gradients in the iron abundance of two galaxies out of the five with sufficient data (NGC 185 and And II). There is only tentative evidence for an alpha abundance radial gradient in NGC 185. We homogeneously compare our results to the Milky Way classical dwarf spheroidals, finding evidence for wider variation in average alpha abundance. In the absence of chemical abundances for the M 31 stellar halo, we compare to the Milky Way stellar halo. A stellar halo comprised of disrupted M31 satellites is too metal-rich and inconsistent with the Milky Way halo alpha abundance distribution even if considering only satellites with predominantly old stellar populations. The M 31 satellite population provides a second system in which to study chemical abundances of dwarf galaxies and reveals a wider variety of abundance patterns than the Milky Way.

We present global metallicity properties, metallicity distribution functions (MDFs) and radial metallicity profiles for the six most luminous M31 dwarf galaxy satellites: M32, NGC 205, NGC 185, NGC 147, Andromeda VII, and Andromeda II. The results presented are the first spectroscopic MDFs for dwarf systems surrounding a host galaxy other than the Milky Way. Our sample consists of individual metallicity measurements for 1243 red giant branch (RGB) member stars spread across these six systems. We determine metallicities based on the strength of the Ca II triplet lines using the empirical calibration of Carrera et al.(2013) which is calibrated over the metallicity range -4 < [Fe/H] <+0.5. We find that these M31 satellites lie on the same luminosity-metallicity relationship as the Milky Way dwarf satellites. We do not find a trend between the internal metallicity spread and galaxy luminosity, contrary to previous studies. The MDF widths of And II and And VII are similar to the MW dwarfs of comparable luminosity, however, our four brightest M31 dwarf are more luminous than any of the MW dwarf spheroidals and have broader MDFs. The MDFs of our six M31 dwarfs are consistent with the leaky box model of chemical evolution, although our metallicity errors allow a wide range of evolution models. We find a significant radial gradient in metallicity in only two of our six systems, NGC 185 and Andromeda II, and flat radial metallicity gradients in the rest of our sample with no observed correlation between rotational support and radial metallicity gradients. While the average properties and radial trends of the M31 dwarf galaxies agree with MW counterparts at similar luminosity, the detailed MDFs are different, particularly at the metal-rich end.

We compare the kinematics of M31's satellite galaxies to the mass profiles of the subhaloes they are expected to inhabit in LCDM. We consider the most massive subhaloes of an approximately M31-sized halo, following the assumption of a monotonic galaxy luminosity-to-subhalo mass mapping. While this abundance matching relation is consistent with the kinematic data for galaxies down to the luminosity of the bright satellites of the Milky Way and M31, it is not consistent with kinematic data for fainter dwarf galaxies (those with L <~ 10^8 Lsun). Comparing the kinematics of M31's dwarf Spheroidal (dSph) satellites to the subhaloes reveals that M31's dSph satellites are too low density to be consistent with the subhaloes' mass profiles. A similar discrepancy has been reported between Milky Way dSphs and their predicted subhaloes, the "too big to fail" problem. By contrast, total mass profiles of the dwarf Elliptical (and similarly bright) satellites are consistent with the subhaloes. However, they suffer from large systematic uncertainties in their dark matter content because of substantial (and potentially dominant) contributions from baryons within their half-light radii.

We examine the star formation properties of bright (~0.1 L*) satellites around isolated ~L* hosts in the local Universe using spectroscopically confirmed systems in the Sloan Digital Sky Survey DR7. Our selection method is carefully designed with the aid of N-body simulations to avoid groups and clusters. We find that satellites are significantly more likely to be quenched than a stellar mass-matched sample of isolated galaxies. Remarkably, this quenching occurs only for satellites of hosts that are themselves quenched: while star formation is unaffected in the satellites of star-forming hosts, satellites around quiescent hosts are more than twice as likely to be quenched than stellar-mass matched field samples. One implication of this is that whatever shuts down star formation in isolated, passive L* galaxies also plays at least an indirect role in quenching star formation in their bright satellites. The previously-reported tendency for "galactic conformity" in color/morphology may be a by-product of this host-specific quenching dichotomy. The S\'ersic indices of quenched satellites are statistically identical to those of field galaxies with the same specific star formation rates, suggesting that environmental and secular quenching give rise to the same morphological structure. By studying the distribution of pairwise velocities between the hosts and satellites, we find dynamical evidence that passive host galaxies reside in dark matter halos that are ~45% more massive than those of star-forming host galaxies of the same stellar mass. We emphasize that even around passive hosts, the mere fact that galaxies become satellites does not typically result in star formation quenching: we find that only ~30% of ~0.1 L* galaxies that fall in from the field are quenched around passive hosts, compared with ~0% around star forming hosts.

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

The prominent broad Fe II emission blends in the spectra of active galactic nuclei have been shown to vary in response to continuum variations, but past attempts to measure the reverberation lag time of the optical Fe II lines have met with only limited success. Here we report the detection of Fe II reverberation in two Seyfert 1 galaxies, NGC 4593 and Mrk 1511, based on data from a program carried out at Lick Observatory in Spring 2011. Light curves for emission lines including H-beta and Fe II were measured by applying a fitting routine to decompose the spectra into several continuum and emission-line components, and we use cross-correlation techniques to determine the reverberation lags of the emission lines relative to V-band light curves. In both cases the measured lag (t_cen) of Fe II is longer than that of H-beta, although the inferred lags are somewhat sensitive to the choice of Fe II template used in the fit. For spectral decompositions done using the Fe II template of Veron-Cetty et al. (2004), we find t_cen(Fe II)/t_cen(H-beta) = 1.9+-0.6 in NGC 4593 and 1.5+-0.3 in Mrk 1511. The detection of highly correlated variations between Fe II and continuum emission demonstrates that the Fe II emission in these galaxies originates in photoionized gas, located predominantly in the outer portion of the broad-line region.

We present Keck/DEIMOS spectroscopy of resolved stars in the M31 satellites And XXVIII & And XXIX. We show that these are likely self-bound galaxies based on 18 and 24 members in And XXVIII & And XXIX, respectively. And XXVIII has a systemic velocity of -331.1 +/- 1.8 km/s and velocity dispersion of 4.9+/-1.6 km/s, implying a mass-to-light ratio (within r_1/2) of ~ 44 +/- 41. And XXIX has a systemic velocity of -194.4 +/- 1.5 km/s and velocity dispersion of 5.7+/-1.2 km/s, implying a mass-to-light ratio (within r_1/2) of ~ 124 +/- 72. The internal kinematics and implied masses of And XXVIII & And XXIX are similar to dwarf spheroidals (dSphs) of comparable luminosities, implying that these objects are dark matter-dominated dwarf galaxies. Despite the large projected distances from their host (380 and 188 kpc), the kinematics of these dSph suggest that they are bound M31 satellites.

We present the surface brightness profile of M31's stellar halo out to a projected radius of 175 kpc. The surface brightness estimates are based on confirmed samples of M31 red giant branch stars derived from Keck/DEIMOS spectroscopic observations. A set of empirical spectroscopic and photometric M31 membership diagnostics is used to identify and reject foreground and background contaminants. This enables us to trace the stellar halo of M31 to larger projected distances and fainter surface brightnesses than previous photometric studies. The surface brightness profile of M31's halo follows a power-law with index -2.2 +/- 0.2 and extends to a projected distance of at least ~175 kpc (~ 2/3 of M31's virial radius), with no evidence of a downward break at large radii. The best-fit elliptical isophotes have b/a=0.94 with the major axis of the halo aligned along the minor axis of M31's disk, consistent with a prolate halo, although the data are also consistent with M31's halo having spherical symmetry. The fact that tidal debris features are kinematically cold is used to identify substructure in the spectroscopic fields out to projected radii of 90 kpc, and investigate the effect of this substructure on the surface brightness profile. The scatter in the surface brightness profile is reduced when kinematically identified tidal debris features in M31 are statistically subtracted; the remaining profile indicates a comparatively diffuse stellar component to M31's stellar halo exists to large distances. Beyond 90 kpc, kinematically cold tidal debris features can not be identified due to small number statistics; nevertheless, the significant field-to-field variation in surface brightness beyond 90 kpc suggests that the outermost region of M31's halo is also comprised to a significant degree of stars stripped from accreted objects.

In order to understand the formation and evolution of dIm galaxies, one needs to understand their three-dimensional structure. We present measurements of the stellar velocity dispersion in NGC 1569, a nearby post-starburst dIm galaxy. The stellar vertical velocity dispersion, $\sigma_{\rm z}$, coupled with the maximum rotational velocity derived from \ion{H}{1} observations, $V_{\rm max}$, gives a measure of how kinematically hot the galaxy is, and, therefore, indicates its structure. We conclude that the stars in NGC 1569 are in a thick disk with a $V_{\rm max} / \sigma_{\rm z}$ = 2.4 $\pm$ 0.7. In addition to the structure, we analyze the ionized gas kinematics from \ion{O}{3} observations along the morphological major axis. These data show evidence for outflow from the inner starburst region and a potential expanding shell near supermassive star cluster (SSC) A. When compared to the stellar kinematics, the velocity dispersion of the stars increase in the region of SSC A supporting the hypothesis of an expanding shell. The stellar kinematics closely follow the motion of the gas. Analysis of high resolution \ion{H}{1} data clearly reveals the presence of an \ion{H}{1} cloud that appears to be impacting the eastern edge of NGC 1569. Also, an ultra-dense \ion{H}{1} cloud can be seen extending to the west of the impacting \ion{H}{1} cloud. This dense cloud is likely the remains of a dense \ion{H}{1} bridge that extended through what is now the central starburst area. The impacting \ion{H}{1} cloud was the catalyst for the starburst, thus turning the dense gas into stars over a short timescale, $\sim$ 1 Gyr. We performed a careful study of the spectral energy distribution using infrared, optical, and ultraviolet photometry producing a state-of-the-art mass model for the stellar disk. This mass modeling shows that stars dominate the gravitational potential in the inner 1 kpc.

We present kinematical profiles and metallicity for the M31 dwarf spheroidal (dSph) satellite galaxy Andromeda II (And II) based on Keck DEIMOS spectroscopy of 531 red giant branch stars. Our kinematical sample is among the largest for any M31 satellite and extends out to two effective radii (r_eff = 5.3' = 1.1 kpc). We find a mean systemic velocity of -192.4+-0.5 km/s and an average velocity dispersion of sigma_v = 7.8+-1.1 km/s. While the rotation velocity along the major axis of And II is nearly zero (<1 km/s), the rotation along the minor axis is significant with a maximum rotational velocity of v_max=8.6+-1.8 km/s. We find a kinematical major axis, with a maximum rotational velocity of v_max=10.9+-2.4 km/s, misaligned by 67 degrees to the isophotal major axis. And II is thus the first dwarf galaxy with evidence for nearly prolate rotation with a v_max/sigma_v = 1.1, although given its ellipticity of epsilon = 0.10, this object may be triaxial. We measured metallicities for a subsample of our data, finding a mean metallicity of [Fe/H] = -1.39+- 0.03 dex and an internal metallicity dispersion of 0.72+-0.03 dex. We find a radial metallicity gradient with metal-rich stars more centrally concentrated, but do not observe a significant difference in the dynamics of two metallicity populations. And II is the only known dwarf galaxy to show minor axis rotation making it a unique system whose existence offers important clues on the processes responsible for the formation of dSphs.

Vista Variables in The Via Lactea (VVV) is an ESO variability survey that is performing observations in near infrared bands (ZYJHKs) towards the Galactic bulge and part of the disk with the completeness limits at least 3 mag deeper than 2MASS. In the present work, we searched in the VVV survey data for background galaxies near the Galactic plane using ZYJHKs photometry that covers 1.636 square degrees. We identified 204 new galaxy candidates by analyzing colors, sizes, and visual inspection of multi-band (ZYJHKs) images. The galaxy candidates colors were also compared with the predicted ones by star counts models considering a more realistic extinction model at the same completeness limits observed by VVV. A comparison of the galaxy candidates with the expected one by Milennium simulations is also presented. Our results increase the number density of known galaxies behind the Milky Way by more than one order of magnitude. A catalog with galaxy properties including ellipticity, Petrosian radii and ZYJHKs magnitudes is provided, as well as comparisons of the results with other surveys of galaxies towards Galactic plane.

The Andromeda galaxy (M31) shows many tidal features in its halo, including the Giant Southern Stream (GSS) and a sharp ledge in surface density on its western side (the W Shelf). Using DEIMOS on the Keck telescope, we obtain radial velocities of M31's giant stars along its NW minor axis, in a radial range covering the W Shelf and extending beyond its edge. In the space of velocity versus radius, the sample shows the wedge pattern expected from a radial shell, which is detected clearly here for the first time. This confirms predictions from an earlier model of formation of the GSS, which proposed that the W Shelf is a shell from the third orbital wrap of the same tidal debris stream that produces the GSS, with the main body of the progenitor lying in the second wrap. We calculate the distortions in the shelf wedge pattern expected from its outward expansion and angular momentum, and show that these effects are echoed in the data. In addition, a hot, relatively smooth spheroid population is clearly present. We construct a bulge-disk-halo N-body model that agrees with surface brightness and kinematic constraints, and combine it with a simulation of the GSS. From the contrasting kinematic signatures of the hot spheroid and shelf components, we decompose the observed stellar metallicity distribution into contributions from each component using a non-parametric mixture model. The shelf component's metallicity distribution matches previous observations of the GSS superbly, further strengthening the evidence they are connected and bolstering the case for a massive progenitor of this stream.

We present a resolved-star spectroscopic survey of 15 dwarf spheroidal (dSph) satellites of the Andromeda Galaxy (M31) as part of the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) project. We filter foreground contamination from Milky Way (MW) stars, noting that MW substructure is evident in this contaminant sample. We also filter M31 halo field giant stars, and identify the remainder as probable dSph members. We then use these members to determine the kinematical properties of the dSphs. For the first time, we confirm that And XVIII, XXI, and XXII show kinematics consistent with bound, dark matter-dominated galaxies. From the velocity dispersions for the full sample of dSphs we determine masses, which we combine with the size and luminosity of the galaxies to produce mass-size-luminosity scaling relations. With these scalings we determine that the M31 dSphs are fully consistent with the MW dSphs, suggesting that the well-studied MW satellite population provides a fair sample for broader conclusions. We also estimate dark matter halo masses of the satellites, and find that there is no sign that the luminosity of these galaxies depends on their dark halo mass, a result consistent with what is seen for MW dwarfs. Two of the M31 dSphs (And XV, XVI) have estimated maximum circular velocities smaller than 12 km/s (to 1sigma), which likely places them within the lowest mass dark matter halos known to host stars (along with Bootes I of the MW). Finally, we use the systemic velocities of the M31 satellites to estimate the total mass of the M31 halo: within 139 kpc, M31's mass is 8 +4.1 - 3.7 x 10^11 corresponding to a virial mass for M31's dark matter halo of 1.2 +0.9 -0.7 x 10^12, consistent with previous results.

We present dynamical modeling of the broad line region (BLR) in the Seyfert 1 galaxy Mrk 50 using reverberation mapping data taken as part of the Lick AGN Monitoring Project (LAMP) 2011. We model the reverberation mapping data directly, constraining the geometry and kinematics of the BLR, as well as deriving a black hole mass estimate that does not depend on a normalizing factor or virial coefficient. We find that the geometry of the BLR in Mrk 50 is a nearly face-on thick disk, with a mean radius of 9.6(+1.2,-0.9) light days, a width of the BLR of 6.9(+1.2,-1.1) light days, and a disk opening angle of 25\pm10 degrees above the plane. We also constrain the inclination angle to be 9(+7,-5) degrees, close to face-on. Finally, the black hole mass of Mrk 50 is inferred to be log10(M(BH)/Msun) = 7.57(+0.44,-0.27). By comparison to the virial black hole mass estimate from traditional reverberation mapping analysis, we find the normalizing constant (virial coefficient) to be log10(f) = 0.78(+0.44,-0.27), consistent with the commonly adopted mean value of 0.74 based on aligning the M(BH)-{\sigma}* relation for AGN and quiescent galaxies. While our dynamical model includes the possibility of a net inflow or outflow in the BLR, we cannot distinguish between these two scenarios.

The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck/DEIMOS spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For the first time, we detect significant spheroid rotation (v_rot ~ 50 km/s) beyond R_proj ~ 5 kpc. The velocity dispersion decreases from about 140 km/s at R_proj = 7 kpc to 120 km/s at R_proj = 14 kpc, consistent to 2 sigma with existing measurements and models. We calculate the probability that a given star is a member of the spheroid and find that the spheroid has a significant presence throughout the spatial extent of our sample. Lastly, we show that the flattening of the spheroid is due to velocity anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our measurements are much farther out (2 - 14 r_eff) than for the comparison samples.

Supernova (SN) 2009ig was discovered 17 hours after explosion by the Lick Observatory Supernova Search, promptly classified as a normal Type Ia SN (SN Ia), peaked at V = 13.5 mag, and was equatorial, making it one of the foremost supernovae for intensive study in the last decade. Here, we present ultraviolet (UV) and optical observations of SN 2009ig, starting about 1 day after explosion until around maximum brightness. Our data include excellent UV and optical light curves, 25 premaximum optical spectra, and 8 UV spectra, including the earliest UV spectrum ever obtained of a SN Ia. SN 2009ig is a relatively normal SN Ia, but does display high-velocity ejecta - the ejecta velocity measured in our earliest spectra (v ~ -23,000 km/s for Si II 6355) is the highest yet measured in a SN Ia. The spectral evolution is very dramatic at times earlier than 12 days before maximum brightness, but slows after that time. The early-time data provide a precise measurement of 17.13 +/- 0.07 days for the SN rise time. The optical color curves and early-time spectra are significantly different from template light curves and spectra used for light-curve fitting and K-corrections, indicating that the template light curves and spectra do not properly represent all Type Ia supernovae at very early times. In the age of wide-angle sky surveys, SNe like SN 2009ig that are nearby, bright, well positioned, and promptly discovered will still be rare. As shown with SN 2009ig, detailed studies of single events can provide significantly more information for testing systematic uncertainties related to SN Ia distance estimates and constraining progenitor and explosion models than large samples of more distant SNe.

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

Like many galaxies of its size, the Milky Way is a disk with prominent spiral arms rooted in a central bar, although our knowledge of its structure and origin is incomplete. Traditional attempts to understand the Galaxy's morphology assume that it has been unperturbed by major external forces. Here we report simulations of the response of the Milky Way to the infall of the Sagittarius dwarf galaxy (Sgr), which results in the formation of spiral arms, influences the central bar and produces a flared outer disk. Two ring-like wrappings emerge towards the Galactic anti-Center in our model that are reminiscent of the low- latitude arcs observed in the same area of the Milky Way. Previous models have focused on Sgr itself to reproduce the dwarf's orbital history and place associated constraints on the shape of the Milky Way gravitational potential, treating the Sgr impact event as a trivial influence on the Galactic disk. Our results show that the Milky Way's morphology is not purely secular in origin and that low-mass minor mergers predicted to be common throughout the Universe probably have a similarly important role in shaping galactic structure.

We use a volume-limited spectroscopic sample of isolated galaxies in the Sloan Digital Sky Survey (SDSS) to investigate the frequency and radial distribution of luminous (M_r <~ -18.3) satellites like the Large Magellanic Cloud (LMC) around ~L* Milky Way analogs and compare our results object-by-object to LCDM predictions based on abundance matching in simulations. We show that 12% of Milky Way-like galaxies host an LMC-like satellite within 75 kpc (projected), and 42 % within 250 kpc (projected). This implies ~10% have a satellite within the distance of the LMC, and ~40% of L* galaxies host a bright satellite within the virialized extent of their dark matter halos. Remarkably, the simulation reproduces the observed frequency, radial dependence, velocity distribution, and luminosity function of observed secondaries exceptionally well, suggesting that LCDM provides an accurate reproduction of the observed Universe to galaxies as faint as L~10^9 Lsun on ~50 kpc scales. When stacked, the observed projected pairwise velocity dispersion of these satellites is sigma~160 km/s, in agreement with abundance-matching expectations for their host halo masses. Finally, bright satellites around L* primaries are significantly redder than typical galaxies in their luminosity range, indicating that environmental quenching is operating within galaxy-size dark matter halos that typically contain only a single bright satellite. This redness trend is in stark contrast to the Milky Way's LMC, which is unusually blue even for a field galaxy. We suggest that the LMC's discrepant color might be further evidence that it is undergoing a triggered star-formation event upon first infall.

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

We examine scaling relations of dispersion-supported galaxies over more than eight orders of magnitude in luminosity by transforming standard fundamental plane parameters into a space of mass (M1/2), radius (r1/2), and luminosity (L1/2). We find that from ultra-faint dwarf spheroidals to giant cluster spheroids, dispersion-supported galaxies scatter about a one-dimensional "fundamental curve" through this MRL space. The weakness of the M1/2-L1/2 slope on the faint end may imply that potential well depth limits galaxy formation in small galaxies, while the stronger dependence on L1/2 on the bright end suggests that baryonic physics limits galaxy formation in massive galaxies. The mass-radius projection of this curve can be compared to median dark matter halo mass profiles of LCDM halos in order to construct a virial mass-luminosity relationship (Mvir-L) for galaxies that spans seven orders of magnitude in Mvir. Independent of any global abundance or clustering information, we find that (spheroidal) galaxy formation needs to be most efficient in halos of Mvir ~ 10^12 Msun and to become inefficient above and below this scale. Moreover, this profile matching technique is most accurate at the high and low luminosity extremes (where dark matter fractions are highest) and is therefore quite complementary to statistical approaches that rely on having a well-sampled luminosity function. We also consider the significance and utility of the scatter about this relation, and find that in the dSph regime observational errors are almost at the point where we can explore the intrinsic scatter in the luminosity-virial mass relation. Finally, we note that purely stellar systems like Globular Clusters and Ultra Compact Dwarfs do not follow the fundamental curve relation. This allows them to be easily distinguished from dark-matter dominated dSph galaxies in MRL space. (abridged)

We predict that there is a population of low-luminosity dwarf galaxies orbiting within the halo of the Milky Way that have surface brightnesses low enough to have escaped detection in star-count surveys. The overall count of stealth galaxies is sensitive to the presence (or lack) of a low-mass threshold in galaxy formation. These systems have luminosities and stellar velocity dispersions that are similar to those of known ultrafaint dwarf galaxies but they have more extended stellar distributions (half light radii greater than about 100 pc) because they inhabit dark subhalos that are slightly less massive than their higher surface brightness counterparts. As a result, the typical peak surface brightness is fainter than 30 mag per square arcsec. One implication is that the inferred common mass scale for Milky Way dwarfs may be an artifact of selection bias. If there is no sharp threshold in galaxy formation at low halo mass, then ultrafaint galaxies like Segue 1 represent the high-mass, early forming tail of a much larger population of objects that could number in the hundreds and have typical peak circular velocities of about 8 km/s and masses within 300 pc of about 5 million solar masses. Alternatively, if we impose a low-mass threshold in galaxy formation in order to explain the unexpectedly high densities of the ultrafaint dwarfs, then we expect only a handful of stealth galaxies in the halo of the Milky Way. A complete census of these objects will require deeper sky surveys, 30m-class follow-up telescopes, and more refined methods to identify extended, self-bound groupings of stars in the halo.

As part of the WIYN High Image Quality Indiana Irvine (WHIQII) survey, we present 123 spectra of emission-line galaxies, selected on intermediate redshift (.4<z<.8) galaxies with blue colors that appear physically compact. The sample includes 15 true Luminous Compact Blue Galaxies (LCBGs) and an additional 27 slightly less extreme emission-line systems. These galaxies represent a highly evolving class that may play an important role in the decline of star formation since z~1, but their exact nature and evolutionary pathways remain a mystery. Here, we use emission lines to determine metallicities and ionization parameters, constraining their intrinsic properties and state of star formation. Some LCBG metallicities are consistent with a "bursting dwarf" scenario, while a substantial fraction of others are not, further confirming that LCBGs are a highly heterogeneous population but are broadly consistent with the intermediate redshift field. In agreement with previous studies, we observe overall evolution in the luminosity-metallicity relation at intermediate redshift. Our sample, and particularly the LCBGs, occupy a region in the empirical R23-O32 plane that differs from luminous local galaxies and is more consistent with dwarf Irregulars at the present epoch, suggesting that cosmic "downsizing" is observable in even the most fundamental parameters that describe star formation. These properties for our sample are also generally consistent with lying between local galaxies and those at high redshift, as expected by this scenario. Surprisingly, our sample exhibits no detectable correlation between compactness and metallicity, strongly suggesting that at these epochs of rapid star formation, the morphology of compact star-forming galaxies is largely transient.

Supernovae have been confirmed to redshift z ~ 1.7 for type Ia (thermonuclear detonation of a white dwarf) and to z ~ 0.7 for type II (collapse of the core of the star). The subclass type IIn supernovae are luminous core-collapse explosions of massive stars and, unlike other types, are very bright in the ultraviolet, which should enable them to be found optically at redshifts z ~ 2 and higher. In addition, the interaction of the ejecta with circumstellar material creates strong, long-lived emission lines that allow spectroscopic confirmation of many events of this type at z ~ 2 for 3 - 5 years after explosion. Here we report three spectroscopically confirmed type IIn supernovae, at redshifts z = 0.808, 2.013 and 2.357, detected in archival data using a method designed to exploit these properties at z ~ 2. Type IIn supernovae directly probe the formation of massive stars at high redshift. The number found to date is consistent with the expectations of a locally measured stellar initial mass function, but not with an evolving initial mass function proposed to explain independent observations at low and high redshift.

We correct the observed Milky Way satellite luminosity function for luminosity bias using published completeness limits for the Sloan Digital Sky Survey DR5. Assuming that the spatial distribution of Milky Way satellites tracks the subhalos found in the Via Lactea LCDM N-body simulation, we show that there should be between ~300 and ~600 satellites within 400 kpc of the Sun that are brighter than the faintest known dwarf galaxies, and that there may be as many as ~ 1000, depending on assumptions. By taking into account completeness limits, we show that the radial distribution of known Milky Way dwarfs is consistent with our assumption that the full satellite population tracks that of subhalos. These results alleviate the primary worries associated with the so-called missing satellites problem in CDM. We show that future, deep wide-field surveys such as SkyMapper, the Dark Energy Survey (DES), PanSTARRS, and the Large Synoptic Survey Telescope (LSST) will deliver a complete census of ultra-faint dwarf satellites out to the Milky Way virial radius, offer new limits on the free-streaming scale of dark matter, and provide unprecedented constraints on the low-luminosity threshold of galaxy formation.