We present MINESweeper, a tool to measure stellar parameters by jointly fitting observed spectra and broadband photometry to model isochrones and spectral libraries. This approach enables the measurement of spectrophotometric distances, in addition to stellar parameters such as Teff, log(g), [Fe/H], [a/Fe], and radial velocity. MINESweeper employs a Bayesian framework and can easily incorporate a variety of priors, including Gaia parallaxes. Mock data are fit in order to demonstrate how the precision of derived parameters depends on evolutionary phase and SNR. We then fit a selection of data in order to validate the model outputs. Fits to a variety of benchmark stars including Procyon, Arcturus, and the Sun result in derived stellar parameters that are in good agreement with the literature. We then fit combined spectra and photometry of stars in the open and globular clusters M92, M13, M3, M107, M71, and M67. Derived distances, [Fe/H], [a/Fe], and log(g)-Teff, relations are in overall good agreement with literature values, although there are trends between metallicity and log(g), within clusters that point to systematic uncertainties at the ~0.1 dex level. Finally, we fit a large sample of stars from the H3 Spectroscopic Survey in which high quality Gaia parallaxes are also available. These stars are fit without the Gaia parallaxes so that the geometric parallaxes can serve as an independent test of the spectrophotometric distances. Comparison between the two reveals good agreement within their formal uncertainties after accounting for the Gaia zero point uncertainties.

In the $\Lambda$CDM paradigm the Galactic stellar halo is predicted to harbor the accreted debris of smaller systems. To identify these systems, the H3 Spectroscopic Survey, combined with $Gaia$, is gathering 6D phase-space and chemical information in the distant Galaxy. Here we present a comprehensive inventory of structure within 50 kpc from the Galactic center using a sample of 5684 giants at $|b|>40^{\circ}$ and $|Z|>2$ kpc. We identify known structures including the high-$\alpha$ disk, the in-situ halo (disk stars heated to eccentric orbits), Sagittarius (Sgr), $Gaia$-Sausage-Enceladus (GSE), the Helmi Streams, Sequoia, and Thamnos. Additionally, we identify the following new structures: (i) Aleph ([Fe/H]$=-0.5$), a low eccentricity structure that rises a surprising 10 kpc off the plane, (ii, iii) Arjuna ([Fe/H]$=-1.2$) and I'itoi ([Fe/H]$<-2$), which comprise the high-energy retrograde halo along with Sequoia, and (iv) Wukong ([Fe/H]$=-1.6$), a prograde phase-space overdensity chemically distinct from GSE. For each structure we provide [Fe/H], [$\alpha$/Fe], and orbital parameters. Stars born within the Galaxy are a major component at $|Z|\sim$2 kpc ($\approx$60$\%$), but their relative fraction declines sharply to $\lesssim$5$\%$ past 15 kpc. Beyond 15 kpc, $>$80$\%$ of the halo is built by two massive ($M_{\star}\sim10^{8}-10^{9}M_{\odot}$) accreted dwarfs: GSE ([Fe/H]$=-1.2$) within 25 kpc, and Sgr ([Fe/H]$=-1.0$) beyond 25 kpc. This explains the relatively high overall metallicity of the halo ([Fe/H]$\approx-1.2$). We attribute $\gtrsim$95$\%$ of the sample to one of the listed structures, pointing to a halo built entirely from accreted dwarfs and heating of the disk.

We present new kinematic data for the Ophiuchus stellar stream. Spectra have been taken of member candidates at the MMT telescope using Hectospec, Hectochelle and Binospec, which provide more than 1800 new velocities. Combined with proper motion measurements of stars in the field by the Gaia - DR2 catalog, we have derived stream membership probabilities, resulting in the detection of more than 200 likely members. These data show the stream extends to more than three times the length shown in the discovery data. A spur to the main stream is also detected. The high resolution spectra allow us to resolve the stellar velocity dispersion, found to be $1.6 \pm 0.3 $ km/sec.

The archaeological record of stars in the Milky Way opens a uniquely detailed window into the early formation and assembly of galaxies. Here we use 11,000 main-sequence turn-off stars with well-measured ages, [Fe/H], [$\alpha$/Fe], and orbits from the H3 Survey and Gaia to time the major events in the early Galaxy. Located beyond the Galactic plane, $1\lesssim |Z|/\rm kpc \lesssim4$, this sample contains three chemically distinct groups: a low metallicity population, and low-$\alpha$ and high-$\alpha$ groups at higher metallicity. The age and orbit distributions of these populations show that: 1) the high-$\alpha$ group, which includes both disk stars and the in-situ halo, has a star-formation history independent of eccentricity that abruptly truncated $8.3\pm0.1$ Gyr ago ($z\simeq1$); 2) the low metallicity population, which we identify as the accreted stellar halo, is on eccentric orbits and its star formation truncated $10.2.^{+0.2}_{-0.1}$ Gyr ago ($z\simeq2$); 3) the low-$\alpha$ population is primarily on low eccentricity orbits and the bulk of its stars formed less than 8 Gyr ago. These results suggest a scenario in which the Milky Way accreted a satellite galaxy at $z\approx2$ that merged with the early disk by $z\approx1$. This merger truncated star formation in the early high-$\alpha$ disk and perturbed a fraction of that disk onto halo-like orbits. The merger enabled the formation of a chemically distinct, low-$\alpha$ disk at $z\lesssim1$. The lack of any stars on halo-like orbits at younger ages indicates that this event was the last significant disturbance to the Milky Way disk.

Globular clusters (GCs) formed when the Milky Way experienced a phase of rapid assembly. We use the wealth of information contained in the Galactic GC population to quantify the properties of the satellite galaxies from which the Milky Way assembled. To achieve this, we train an artificial neural network on the E-MOSAICS cosmological simulations of the co-formation and co-evolution of GCs and their host galaxies. The network uses the ages, metallicities, and orbital properties of GCs that formed in the same progenitor galaxies to predict the stellar masses and accretion redshifts of these progenitors. We apply the network to Galactic GCs associated with five progenitors: {\it Gaia}-Enceladus, the Helmi streams, Sequoia, Sagittarius, and the recently discovered, `low-energy' GCs, which provide an excellent match to the predicted properties of the enigmatic galaxy `Kraken'. The five galaxies cover a narrow stellar mass range [$M_\star=(0.6{-}4.6)\times10^8~{\rm M}_\odot$], but have widely different accretion redshifts ($z_{\rm acc}=0.57{-}2.65$). All accretion events represent minor mergers, but Kraken likely represents the most major merger ever experienced by the Milky Way, with stellar and virial mass ratios of $r_{M_\star}=1$:$31^{+34}_{-16}$ and $r_{M_{\rm h}}=1$:$7^{+4}_{-2}$, respectively. The progenitors match the $z=0$ relation between GC number and halo virial mass, but have elevated specific frequencies, suggesting an evolution with redshift. Even though these progenitors likely were the Milky Way's most massive accretion events, they contributed a total mass of only $\log{(M_{\rm \star,tot}/{\rm M}_\odot)}=9.0\pm0.1$, similar to the stellar halo. This implies that the Milky Way grew its stellar mass mostly by in-situ star formation. We conclude by organising these accretion events into the most detailed reconstruction to date of the Milky Way's merger tree.

Thin stellar streams, formed from the tidal disruption of globular clusters, are important gravitational tools, sensitive to both global and small-scale properties of dark matter. The Palomar 5 stellar stream (Pal 5) is an exemplar stream within the Milky Way: Its $\sim 20^\circ$ tidal tails connect back to the progenitor cluster, and the stream has been used to study the shape, total mass, and substructure fraction of the dark matter distribution of the Galaxy. However, most details of the phase-space distribution of the stream are not fully explained, and dynamical models that use the stream for other inferences are therefore incomplete. Here we aim to measure distance and kinematic properties along the Pal 5 stream in order to motivate improved models of the system. We use a large catalog of RR Lyrae-type stars (RRLs) with astrometric data from the Gaia mission to probabilistically identify RRLs in the Pal 5 stream. RRLs are useful because they are intrinsically-luminous standard candles and their distances can be inferred with small relative precision ($\sim3\%$). By building a probabilistic model of the Pal 5 cluster and stream in proper motion and distance, we find 27 RRLs consistent with being members of the cluster (10) and stream (17). Using these RRLs, we detect gradients in distance and proper motion along the stream, and provide an updated measurement of the distance to the Pal 5 cluster using the RRLs, $d = 20.6 \pm 0.2~\textrm{kpc}$. We provide a catalog of Pal 5 RRLs with inferred membership probabilities for future modeling work.

The $100^\circ$-long thin stellar stream in the Milky Way halo, GD-1, has an ensemble of features that may be due to dynamical interactions. Using high-resolution MMT/Hectochelle spectroscopy we show that a spur of GD-1-like stars outside of the main stream are kinematically and chemically consistent with the main stream. In the spur, as in the main stream, GD-1 has a low intrinsic radial velocity dispersion, $\sigma_{V_r}\lesssim1\,\rm km\,s^{-1}$, is metal-poor, $\rm [Fe/H]\approx-2.3$, with little $\rm [Fe/H]$ spread and some variation in $\rm [\alpha/Fe]$ abundances, which point to a common globular cluster progenitor. At a fixed location along the stream, the median radial velocity offset between the spur and the main stream is smaller than $0.5\,\rm km\,s^{-1}$, comparable to the measurement uncertainty. A flyby of a massive, compact object can change orbits of stars in a stellar stream and produce features like the spur observed in GD-1. In this scenario, the radial velocity of the GD-1 spur relative to the stream constrains the orbit of the perturber and its current on-sky position to $\approx5,000\,\rm deg^2$. The family of acceptable perturber orbits overlaps the stellar and dark-matter debris of the Sagittarius dwarf galaxy in present-day position and velocity. This suggests that GD-1 may have been perturbed by a globular cluster or an extremely compact dark-matter subhalo formerly associated with Sagittarius.

Stars that escape globular clusters form tidal tails that are predominantly shaped by the global distribution of mass in the Galaxy, but also preserve a historical record of small-scale perturbations. Using deep $grz$ photometry from DECaLS, we present highly probable members of the tidal tails associated with the disrupting globular cluster Palomar 5. These data yield the cleanest view of a stellar stream beyond $\sim20\,\rm kpc$ and reveal: (1) a wide, low surface-brightness extension of the leading tail; (2) significant density variations along the stream; and (3) sharp changes in the direction of both the leading and the trailing tail. In the fiducial Milky Way model, a rotating bar perturbs the Palomar 5 tails and can produce streams with similar width and density profiles to those observed. However, the deviations of the stream track in this simple model do not match those observed in the Palomar 5 trailing tail, indicating the need for an additional source of perturbation. These discoveries open up the possibility of measuring the population of perturbers in the Milky Way, including dark-matter subhalos, with an ensemble of stellar streams and deep photometry alone.

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.

The timing argument provides a lower limit on the mass of the Milky Way. We find, using a sample of 32 stars at $R > 60$ kpc drawn from the H3 Spectroscopic Survey and mock catalogs created from published numerical simulations, that M$_{200} > 0.91\times 10^{12}$ M$_\odot$ with 90% confidence. We recommend using this limit to refine the allowed prior mass range in more complex and sophisticated statistical treatments of Milky Way dynamics. The use of such a prior would have significantly reduced many previously published uncertainty ranges. Our analysis suggests that the most likely value of M$_{200}$ is $\sim 1.4 \times 10^{12}$ M$_\odot$, but establishing this as the Milky Way mass requires a larger sample of outer halo stars and a more complete analysis of the inner halo stars in H3. The imminent growth in the sample of outer halo stars due to ongoing and planned surveys will make this possible.

The Galactic stellar halo is predicted to have formed at least partially from the tidal disruption of accreted dwarf galaxies. This assembly history should be detectable in the orbital and chemical properties of stars. The H3 Survey is obtaining spectra for 200,000 stars, and, when combined with Gaia data, is providing detailed orbital and chemical properties of Galactic halo stars. Unlike previous surveys of the halo, the H3 target selection is based solely on magnitude and Gaia parallax; the survey therefore provides a nearly unbiased view of the entire stellar halo at high latitudes. In this paper we present the distribution of stellar metallicities as a function of Galactocentric distance and orbital properties for a sample of 4232 kinematically-selected halo giants to 100 kpc. The stellar halo is relatively metal-rich, [Fe/H]=-1.2, and there is no discernable metallicity gradient over the range $6<R_{\rm gal}<100$ kpc. However, the halo metallicity distribution is highly structured including distinct metal-rich and metal-poor components at $R_{\rm gal}<10$ kpc and $R_{\rm gal}>30$ kpc, respectively. Metal-poor stars with [Fe/H]$<-2$ are a small population of the halo at all distances and orbital categories. We associate the "in-situ" stellar halo with stars displaying thick-disk chemistry on halo-like orbits; such stars are confined to $|z|<10$ kpc. The majority of the stellar halo is resolved into discrete features in orbital-chemical space, suggesting that the bulk of the stellar halo formed from the accretion and tidal disruption of dwarf galaxies. (ABRIDGED)

In 2008 it was reported that the stellar stream of the edge-on spiral NGC5907 loops twice around the galaxy, enveloping it in a giant corkscrew-like structure. Here we present imaging of this iconic object with the Dragonfly Telephoto Array, reaching a $1\sigma$ surface brightness level of $\mu_g\approx 30.5$ mag/arcsec$^2$ on spatial scales of 1' (the approximate width of the stream). We find a qualitatively different morphology from that reported in the 2008 study. The Dragonfly data do not show two loops but a single curved stream with a total length of 45' (220 kpc). The surface brightness of the stream ranges from $\mu_g \approx 27.6$ mag/arcsec$^2$ to $\mu_g\approx 28.8$ mag/arcsec$^2$, and it extends significantly beyond the region where tidal features had previously been detected. We find a density enhancement near the luminosity-weighted midpoint of the stream which we identify as the likely remnant of a nearly-disrupted progenitor galaxy. A restricted N-body simulation provides a qualitative match to the detected features. In terms of its spatial extent and stellar mass the stream is similar to Sagittarius, and our results demonstrate the efficacy of low surface brightness-optimized telescopes for obtaining maps of such large streams outside the Local Group. The census of these rare, relatively high mass events complements the census of common, low mass ones that is provided by studies of streams in the Milky Way halo.

Modern theories of galaxy formation predict that the Galactic stellar halo was hierarchically assembled from the accretion and disruption of smaller systems. This hierarchical assembly is expected to produce a high degree of structure in the combined phase and chemistry space; this structure should provide a relatively direct probe of the accretion history of our Galaxy. Revealing this structure requires precise 3D positions (including distances), 3D velocities, and chemistry for large samples of stars. The Gaia satellite is delivering proper motions and parallaxes for >1 billion stars to G~20. However, radial velocities and metallicities will only be available to G~15, which is insufficient to probe the outer stellar halo (>10 kpc). Moreover, parallaxes will not be precise enough to deliver high-quality distances for stars beyond ~10 kpc. Identifying accreted systems throughout the stellar halo therefore requires a large ground-based spectroscopic survey to complement Gaia. Here we provide an overview of the H3 Stellar Spectroscopic Survey, which will deliver precise stellar parameters and spectrophotometric distances for 200,000 stars to r=18. Spectra are obtained with the Hectochelle instrument at the MMT, which is configured for the H3 Survey to deliver resolution R~23,000 spectra covering the wavelength range 5150A-5300A. The survey is optimized for stellar halo science and therefore focuses on high Galactic latitude fields (|b|>30 deg.), sparsely sampling 15,000 sq. degrees. Targets are selected on the basis of Gaia parallaxes, enabling very efficient selection of bone fide halo stars. The survey began in the Fall of 2017 and has collected 88,000 spectra to-date. All of the data, including the derived stellar parameters, will eventually be made publicly available via the survey website: h3survey.rc.fas.harvard.edu.

We present a model for the interaction of the GD-1 stellar stream with a massive perturber that naturally explains many of the observed stream features, including a gap and an off-stream spur of stars. The model involves an impulse by a fast encounter, after which the stream grows a loop of stars at different orbital energies. At specific viewing angles, this loop appears offset from the stream track. A quantitative comparison of the spur and gap features prefers models where the perturber is in the mass range of $10^6\,\rm M_\odot$ to $10^8\,\rm M_\odot$. Orbit integrations back in time show that the stream encounter could not have been caused by any known globular cluster or dwarf galaxy with a determined orbit, and mass, size and impact-parameter arguments show that it could not have been caused by a molecular cloud in the Milky Way disk. The most plausible explanation for the gap-and-spur structure is an encounter with a dark-matter substructure, like those predicted to populate galactic halos in LCDM cosmology. However, the expected densities of LCDM subhalos in this mass range and in this part of the Milky Way are $2-3\,\sigma$ lower than the inferred density of the GD-1 perturber. This observation opens up the possibility that detailed observations of streams could measure the mass spectrum of dark-matter substructures and even identify individual substructures and their orbits in the Galactic halo.

In simple models of the Milky Way, tidally disrupting satellites produce long and thin---nearly one-dimensional---stellar streams. Using astrometric data from the Gaia second data release and photometry from the Dark Energy Survey, we demonstrate that the Jhelum stream, a stellar stream in the inner halo, is a two-dimensional structure. The spatial distribution of highly probable Jhelum members reveals a dense thin component and an associated diffuse, spatially offset component. These two spatial components have indistinguishable proper motions (at $\sigma\sim1\,\rm mas\,yr^{-1}$ level) and a similar ratio of blue straggler to blue horizontal branch stars, which indicates a common origin for the two components. The best-fit orbit to the narrow component (pericenter $8\,\rm kpc$, apocenter $24\,\rm kpc$), however, does not explain the wide component of the Jhelum stream. On the other hand, an older orbital wrap of Jhelum's orbit traces the Indus stream, indicating a possible connection between these two structures and additional complexity in Jhelum's formation. Substructure in the Jhelum progenitor or precession of its tidal debris in the Milky Way potential may explain the observed structure of Jhelum. Future spectroscopic data will enable discrimination between these "nature" and "nurture" formation scenarios. Jhelum adds to the growing list of cold stellar streams that display complex morphologies and promise to reveal the dynamical history of the Milky Way.

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

It is challenging to reliably identify stars that were born together outside of actively star-forming regions and bound stellar systems. However, co-natal stars should be present throughout the Galaxy, and their demographics can shed light on the clustered nature of star formation and the dynamical state of the disk. In previous work we presented a set of simulations of the Galactic disk that followed the clustered formation and dynamical evolution of 4 billion individual stars over the last 5 Gyr. The simulations predict that a high fraction of co-moving stars with physical and 3D velocity separation of $\Delta r < 20$ pc and $\Delta v < 1.5$ km s$^{-1}$ are co-natal. In this \textit{Letter}, we use \textit{Gaia} DR2 and LAMOST DR4 data to identify and study co-moving pairs. We find that the distribution of relative velocities and separations of pairs in the data is in good agreement with the predictions from the simulation. We identify 111 co-moving pairs in the Solar neighborhood with reliable astrometric and spectroscopic measurements. These pairs show a strong preference for having similar metallicities when compared to random field pairs. We therefore conclude that these pairs were very likely born together. The simulations predict that co-natal pairs originate preferentially from high-mass and relatively young ($< 1$ Gyr) star clusters. \textit{Gaia} will eventually deliver well-determined metallicities for the brightest stars, enabling the identification of thousands of co-natal pairs due to disrupting star clusters in the solar neighborhood.

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.

The clustered nature of star formation should produce a high degree of structure in the combined phase and chemical space in the Galactic disk. To date, observed structure of this kind has been mostly limited to bound clusters and moving groups. In this paper we present a new dynamical model of the Galactic disk that takes into account the clustered nature of star formation. This model predicts that the combined phase and chemical space is rich in substructure, and that this structure is sensitive to both the precise nature of clustered star formation and the large-scale properties of the Galaxy. The model self-consistently evolves 4 billion stars over the last 5 Gyr in a realistic potential that includes an axisymmetric component, a bar, spiral arms, and giant molecular clouds (GMCs). All stars are born in clusters with an observationally-motivated range of initial conditions. As direct \textit{N}-body calculations for billions of stars is computationally infeasible, we have developed a method of initializing star cluster particles to mimic the effects of direct \textit{N}-body effects, while the actual orbit integrations are treated as test particles within the analytic potential. We demonstrate that the combination of chemical and phase space information is much more effective at identifying truly co-natal populations than either chemical or phase space alone. Furthermore, we show that co-moving pairs of stars are very likely to be co-natal if their velocity separation is $< 2$ km s$^{-1}$ and their metallicity separation is $< 0.05$ dex. The results presented here bode well for harnessing the synergies between \textit{Gaia} and spectroscopic surveys to reveal the assembly history of the Galactic disk.

Tidally-disrupted globular clusters are transformed into thin, dynamically-cold streams of stars that are extremely valuable tracers of the large- and small-scale distribution of mass in the Galaxy. Using data from the Gaia second data release combined with Pan-STARRS photometry, we present a sample of highly-probable members of the longest cold stream in the Milky Way, GD-1. The resulting map of GD-1: (1) extends the apparent length of the stream by $20^\circ$, (2) reveals plausible locations for the progenitor, (3) detects high-contrast gaps along the stream, and (4) indicates the existence of stream members perturbed off the main stream track. These discoveries are only possible because of the exquisite astrometry from Gaia, which permits a clean separation of the stream from Milky Way stars. The additional length and a proper treatment of the progenitor will aid in dynamical modeling of GD-1 for mapping the large-scale dark matter distribution. The complex morphology of the stream points to a turbulent history; detailed phase-space properties of the perturbed stream members could potentially constrain dark matter substructure in the Milky Way.

We present a color-based method for identifying red giants based on Pan-STARRS grz and WISE W1 and W2 photometry. We utilize a subsample of bright stars with precise parallaxes from Gaia DR2 to verify that the color-based selection reliably separates dwarfs from giants. The selection is conservative in the sense that contamination is small (~30%) but not all giants are included (the selection primarily identifies K giants). The color-based selection can be applied to stars brighter than $W1\approx16$, more than two magnitudes fainter than techniques relying on shallower 2MASS photometry. Many streams and clouds are visible in the resulting sky maps, especially when binned by Gaia DR2 proper motions, including the Sagittarius stream, the Hercules-Aquila Cloud, the Eastern Banded Structure, Monoceros, and the Virgo Overdensity. In addition to the characterization of new and known stellar streams, we expect that this method for selecting red giants will enable detailed analysis of the diffuse stellar halo to distances exceeding 100 kpc.

Cold stellar streams---produced by tidal disruptions of clusters---are long-lived, coherent dynamical features in the halo of the Milky Way. Due to their different ages and different positions in phase space, different streams tell us different things about the Galaxy. Here we employ a Cramer--Rao (CRLB) or Fisher-matrix approach to understand the quantitative information content in eleven known streams (ATLAS, GD-1, Hermus, Kwando, Orinoco, PS1A, PS1C, PS1D, PS1E, Sangarius and Triangulum). This approach depends on a generative model, which we have developed previously, and which permits calculation of derivatives of predicted stream properties with respect to Galaxy and stream parameters. We find that in simple analytic models of the Milky Way, streams on eccentric orbits contain the most information about the halo shape. For each stream, there are near-degeneracies between dark-matter-halo properties and parameters of the bulge, the disk, and the stream progenitor, but simultaneous fitting of multiple streams will constrain all parameters at the percent level. At this precision, simulated dark matter halos deviate from simple analytic parametrizations, so we add an expansion of basis functions to give the gravitational potential more freedom. As freedom increases, the information about the halo reduces overall, and it becomes more localized to the current position of the stream. In the limit of high model freedom, a stellar stream appears to measure the local acceleration at its current position; this motivates thinking about future non-parametric approaches. The CRLB formalism also permits us to assess the value of future measurements of stellar velocities, distances, and proper motions. We show that kinematic measurements of stream stars are essential for producing competitive constraints on the distribution of dark matter, which bodes well for stream studies in the age of Gaia.

Stellar models typically use the mixing length approximation as a way to implement convection in a simplified manner. While conventionally the value of the mixing length parameter, $\alpha$, used is the solar calibrated value, many studies have shown that other values of $\alpha$ are needed to properly model stars. This uncertainty in the value of the mixing length parameter is a major source of error in stellar models and isochrones. Using asteroseismic data, we determine the value of the mixing length parameter required to properly model a set of about 450 stars ranging in $\log g$, $T_{\mathrm{eff}}$, and $\mathrm{[Fe/H]}$. The relationship between the value of $\alpha$ required and the properties of the star is then investigated. For Eddington atmosphere, non-diffusion models, we find that the value of $\alpha$ can be approximated by a linear model, in the form of $\alpha/\alpha_{\odot}=5.426 -0.101 \log (g) -1.071 \log (T_{\mathrm{eff}}) + 0.437 (\mathrm{[Fe/H]})$. This process is repeated using a variety of model physics as well as compared to previous studies and results from 3D convective simulations.

We use the first Gaia data release, combined with RAVE and APOGEE spectroscopic surveys, to investigate the origin of halo stars within <~3 kpc from the Sun. We identify halo stars kinematically, as moving with a relative speed of at least 220 km/s with respect to the local standard of rest. These stars are in general more metal-poor than the disk, but surprisingly, half of our halo sample is comprised of stars with [Fe/H]>-1. The orbital directions of these metal-rich halo stars are preferentially aligned with the disk rotation, in sharp contrast with the isotropic orbital distribution of the more metal-poor halo stars. We find similar properties in the Latte cosmological zoom-in simulation of a Milky Way-like galaxy from the FIRE project. In Latte, metal-rich halo stars formed primarily inside of the solar circle, while lower-metallicity halo stars preferentially formed at larger distances (extending beyond the virial radius). This suggests that metal-rich halo stars in the Solar neighborhood in fact formed in-situ within the Galactic disk rather than having been accreted from satellite systems. These stars, currently on halo-like orbits, therefore have likely undergone substantial radial migration/heating.

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.

We announce the discovery of the Aquarius~2 dwarf galaxy, a new distant satellite of the Milky Way, detected on the fringes of the VST ATLAS and the SDSS surveys. The object was originally identified as an overdensity of Red Giant Branch stars, but chosen for subsequent follow-up based on the presence of a strong Blue Horizontal Branch, which was also used to measure its distance of $\sim 110$ kpc. Using deeper imaging from the IMACS camera on the 6.5m Baade and spectroscopy with DEIMOS on Keck, we measured the satellite's half-light radius $5.1\pm 0.8$ arcmin, or $\sim 160$ pc at this distance, and its stellar velocity dispersion of $5.4^{+3.4}_{-0.9}$ km s$^{-1}$. With $\mu=30.2$ mag arcsec$^{-2}$ and $M_V=-4.36$, the new satellite lies close to two important detection limits: one in surface brightness; and one in luminosity at a given distance, thereby making Aquarius~2 one of the hardest dwarfs to find.

Tidal streams of globular clusters are ideal tracers of the Galactic gravitational potential. Compared to the few known, complex and diffuse dwarf-galaxy streams, they are kinematically cold, have thin morphologies and are abundant in the halo of the Milky Way. Their coldness and thinness in combination with potential epicyclic substructure in the vicinity of the stream progenitor turns them into high-precision scales. With the example of Palomar 5, we demonstrate how modeling of a globular cluster stream allows us to simultaneously measure the properties of the disrupting globular cluster, its orbital motion, and the gravitational potential of the Milky Way.

Using the example of the tidal stream of the Milky Way globular cluster Palomar 5 (Pal 5), we demonstrate how observational data on streams can be efficiently reduced in dimensionality and modeled in a Bayesian framework. Our approach combines detection of stream overdensities by a Difference-of-Gaussians process with fast streakline models, a continuous likelihood function built from these models, and inference with MCMC. By generating $\approx10^7$ model streams, we show that the geometry of the Pal 5 debris yields powerful constraints on the solar position and motion, the Milky Way and Pal 5 itself. All 10 model parameters were allowed to vary over large ranges without additional prior information. Using only SDSS data and a few radial velocities from the literature, we find that the distance of the Sun from the Galactic Center is $8.30\pm0.25$ kpc, and the transverse velocity is $253\pm16$ km/s. Both estimates are in excellent agreement with independent measurements of these quantities. Assuming a standard disk and bulge model, we determine the Galactic mass within Pal 5's apogalactic radius of 19 kpc to be $(2.1\pm0.4)\times10^{11}$ M$_\odot$. Moreover, we find the potential of the dark halo with a flattening of $q_z = 0.95^{+0.16}_{-0.12}$ to be essentially spherical within the radial range that is effectively probed by Pal 5. We also determine Pal 5's mass, distance and proper motion independently from other methods, which enables us to perform vital cross-checks. We conclude that with more observational data and by using additional prior information, the precision of this method can be significantly increased.

We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compare, in 6D phase space, an "observed" stream to models created in trial analytic potentials. We analyze a large sample of streams evolved in the Via Lactea II (VL2) simulation, which represents a realistic Galactic halo potential. The recovered potential parameters are in agreement with the best fit to the global, present-day VL2 potential. However, merely assuming an analytic potential limits the dark matter halo mass measurement to an accuracy of 5 to 20%, depending on the choice of analytic parametrization. Collectively, mass estimates using streams from our sample reach this fundamental limit, but individually they can be highly biased. Individual streams can both under- and overestimate the mass, and the bias is progressively worse for those with smaller perigalacticons, motivating the search for tidal streams at galactocentric distances larger than 70 kpc. We estimate that the assumption of a static and smooth dark matter potential in modeling of the GD-1 and Pal5-like streams introduces an error of up to 50% in the Milky Way mass estimates.

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {\it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$\propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($\propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $\sim 10^{-9.3} h^2 {\rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}\msunh < M_h < 10^{12.3}\msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged]

We present evidence for a new Milky Way stellar tidal stream in the direction of the Andromeda and Triangulum (M31 and M33) galaxies. Using a matched-filter technique, we search the Sloan Digital Sky Survey Data Release 8 by creating stellar density maps which probe the Milky Way halo at distances between 8 and 40 kpc. A visual search of these maps recovers all of the major known stellar streams, as well as a new stream in the direction of M31/M33 which we name the Triangulum stream. The stream spans 0.2 deg by 12 deg on the sky, or 75 pc by 5.5 kpc in physical units with a best fitting distance of 26+/-4 kpc. The width of the stream is consistent with being the tidal remnant of a globular cluster. A color magnitude diagram of the stream region shows an overdensity which, if identified as a main sequence turn-off, corresponds to an old (~12 Gyr) and metal-poor ([Fe/H] \sim -1.0 dex) stellar population. Future kinematic studies of this and similar cold streams will provide tight constraints on the shape of the Galactic gravitational potential.

Stellar models generally use simple parametrizations to treat convection. The most widely used parametrization is the so-called "Mixing Length Theory" where the convective eddy sizes are described using a single number, \alpha, the mixing-length parameter. This is a free parameter, and the general practice is to calibrate \alpha using the known properties of the Sun and apply that to all stars. Using data from NASA's Kepler mission we show that using the solar-calibrated \alpha is not always appropriate, and that in many cases it would lead to estimates of initial helium abundances that are lower than the primordial helium abundance. Kepler data allow us to calibrate \alpha for many other stars and we show that for the sample of stars we have studied, the mixing-length parameter is generally lower than the solar value. We studied the correlation between \alpha and stellar properties, and we find that \alpha increases with metallicity. We therefore conclude that results obtained by fitting stellar models or by using population-synthesis models constructed with solar values of \alpha are likely to have large systematic errors. Our results also confirm theoretical expectations that the mixing-length parameter should vary with stellar properties.

We use the Eighth Data Release of Sloan Digital Sky Survey (SDSS DR8) catalog with its additional sky coverage of the southern Galactic hemisphere, to measure the extent and study the nature of the Virgo Overdensity (VOD; Juric et al. 2008). The data show that the VOD extends over no less than 2000 deg^2, with its true extent likely closer to 3000 deg^2. We test whether the VOD can be attributed to a tilt in the stellar halo ellipsoid with respect to the plane of the Galactic disk and find that the observed symmetry of the north-south Galactic hemisphere star counts excludes this possibility. We argue that the Virgo Overdensity, in spite of its wide area and cloud-like appearance, is still best explained by a minor merger. Its appearance and position is qualitatively similar to a near perigalacticon merger event and, assuming that the VOD and the Virgo Stellar Stream share the same progenitor, consistent with the VSS orbit determined by Casetti-Dinescu et al. (2009).

We study hydrogenated cations of two polycyclic hydrocarbon molecules as models of hydrogenated organic species that form in the interstellar medium. Optical spectra of the hydrogenated naphthalene cation Hn-C10H8^+ for n=1,2, and 10, as well as of the astrobiologically interesting hydrogenated proflavine cation Hn-C13H11N3^+ for n=1 and 14, are calculated. The pseudopotential time dependent density functional theory is used. It is found that the fully hydrogenated proflavine cation H14-C13H11N3^+ shows a broad spectrum where positions of individual lines are almost lost. The positions of lines, their shapes, and intensities change in hydronaphthalene and hydroproflavine cations showing that hydrogen additions induce substantially different optical spectra in a comparison to base polycyclic hydrocarbon cations. One calculated line in the visible spectrum of H10-C10H8^+, and one in the visible spectrum of H-C13H11N3^+ are close to the measured diffuse interstellar bands. We also present positions of near-UV lines.

Motivated by possible astrophysical and biological applications we calculate visible and near UV spectral lines of proflavine (C13H11N3, 3,6-diaminoacridine) in vacuum, as well as its anion, cation, and dication. The pseudopotential density functional and time-dependent density functional methods are used. We find a good agreement in spectral line positions calculated by two real-time propagation methods and the Lanczos chain method. Spectra of proflavine and its ions show characteristic UV lines which are good candidates for a detection of these molecules in interstellar space and various biological processes.