The vast majority of the mass in the Milky Way (MW) is in dark matter (DM); we therefore cannot directly observe the MW mass distribution, and have to use tracer populations in order to infer properties of the MW DM halo. However, MW halo tracers do not only feel the gravitational influence of the MW itself. Tracers can also be affected by MW satellites; Garavito-Camargo et al. (2019) (hereafter GC19) demonstrate that the Large Magellanic Cloud (LMC) induces a density wake in the MW DM, resulting in large scale kinematic patterns in the MW stellar halo. In this work, we use spherical harmonic expansion (SHE) of the velocity fields of simulated stellar halos in an effort to disentangle perturbations on large scales (e.g., due to the LMC itself as well as the LMC-induced DM wake) and small scales (due to substructure). Using the GC19 simulations, we demonstrate how the different terms in the SHE of the stellar velocity field reflect the different wake components, and show that these signatures are a strong function of the LMC mass. An exploration of model halos built from accreted dwarfs Bullock & Johnston (2005) suggests that stellar debris from massive, recent accretion events can produce much more power in the velocity angular power spectra than the perturbation from the LMC-induced wake. We therefore consider two models for the Sagittarius (Sgr) stream -- the most recent, massive accretion event in the MW apart from the LMC -- and find that the angular power on large scales is generally dominated by the LMC-induced wake, even when Sgr is included. We conclude that SHE of the MW stellar halo velocity field may therefore be a useful tool in quantifying the response of the MW DM halo to the LMC's infall.

$Gaia$ DR2 has provided an unprecedented wealth of information about the positions and motions of stars in our Galaxy, and has highlighted the degree of disequilibria in the disc. As we collect data over a wider area of the disc it becomes increasingly appealing to start analysing stellar actions and angles, which specifically label orbit space, instead of their current phase space location. Conceptually, while $\bar{x}$ and $\bar{v}$ tell us about the potential and local interactions, grouping in action puts together stars that have similar frequencies and hence similar responses to dynamical effects occurring over several orbits. Grouping in actions and angles refines this further to isolate stars which are travelling together through space and hence have shared histories. Mixing these coordinate systems can confuse the interpretation. For example, it has been suggested that by moving stars to their guiding radius, the Milky Way spiral structure is visible as ridge-like overdensities in the $Gaia$ data \citep{Khoperskov+19b}. However, in this work, we show that these features are in fact the known kinematic moving groups, both in the $L_z-\phi$ and the $v_{\mathrm{R}}-v_{\phi}$ planes. Using simulations we show how this distinction will become even more important as we move to a global view of the Milky Way. As an example, we show that the radial velocity wave seen in the Galactic disc in $Gaia$ and APOGEE should become stronger in the action-angle frame, and that it can be reproduced by transient spiral structure.

We present a detailed analysis of the disk and circumstellar environment of the forming O-type star AFGL 4176 mm1, placing results from the Atacama Large Millimeter/submillimeter Array (ALMA) into context with multiwavelength data. With ALMA, we detect seventeen 1.2 mm continuum sources within 5$''$ (21,000 au) of AFGL 4176 mm1. We find that mm1 has a spectral index of 3.4$\pm$0.2 across the ALMA band, with $>$87% of its 1.2 mm continuum emission from dust. The source mm2, projected 4200 au from mm1, may be a companion or a blueshifted knot in a jet. We also explore the morphological differences between the molecular lines detected with ALMA, finding 203 lines from 25 molecules, which we categorize into several morphological types. Our results show that AFGL 4176 mm1 provides an example of a forming O-star with a large and chemically complex disk, which is mainly traced by nitrogen-bearing molecules. Lines that show strong emission on the blueshifted side of the disk are predominantly oxygen-bearing, which we suggest are tracing a disk accretion shock. The molecules C$^{34}$S, H$_2$CS and CH$_{3}$CCN trace a slow wide-angle wind or dense structures in the outflow cavity walls. With the Australia Telescope Compact Array (ATCA), we detect a compact continuum source ($<$2000 $\times$ 760 au) at 1.2 cm, associated with mm1, of which $>$96% is from ionized gas. The ATCA NH$_3$(1,1) and (2,2) emission traces a large-scale (r$\sim$0.5 pc) rotating toroid with the disk source mm1 in the blueshifted part of this structure offset to the northwest.

The Galactic bulge is a tumultuous dense region of space, packed with stars separated by far smaller distances than those in the Solar neighborhood. A quantification of the frequency and proximity of close stellar encounters in this environment dictates the exchange of material, disruption of planetary orbits, and threat of sterilizing energetic events. We present estimated encounter rates for stars in the Milky Way bulge found using a combination of numerical and analytical methods. By integrating the orbits of bulge stars with varying orbital energy and angular momentum to find their positions over time, we were able to estimate how many close stellar encounters the stars should experience as a function of orbit shape. We determined that ~80% of bulge stars have encounters within 1000 AU and that half of bulge stars will have >35 such encounters, both over a gigayear. Our work has interesting implications for the long-term survivability of planets in the Galactic bulge.

We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

A major goal in the field of galaxy formation is to understand the formation of the Milky Way's disk. The first step toward doing this is to empirically describe its present state. We use the new high-dimensional dataset of 19 abundances from 27,135 red clump APOGEE stars to examine the distribution of clusters defined using abundances. We explore different dimensionality reduction techniques and implement a non-parametric agglomerate hierarchical clustering method. We see that groups defined using abundances are spatially separated, as a function of age. Furthermore, the abundance groups represent different distributions in the [Fe/H]-age plane. Ordering our clusters by age reveals patterns suggestive of the sequence of chemical enrichment in the disk over time. Our results indicate that a promising avenue to trace the details of the disk's assembly is via a full interpretation of the empirical connections we report.

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

As the diversity of people in higher education grows, Universities are struggling to provide inclusive environments that nurture the spirit of free inquiry in the presence of these differences. At the extreme, the value of diversity is under attack as a few, vocal academics use public forums to question the innate intellectual abilities of certain demographic groups. Throughout my career as an astronomer, from graduate student, through professor to department chair, I have witnessed these struggles firsthand. Exclusive cultures result in lost opportunities in the form of unfulfilled potential of all members of the institution - students, administrators and faculty alike. How to move steadily towards inclusion is an unsolved problem that hampers the advancement of knowledge itself. As every scientist knows, problem definition is an essential feature of problem solution. This article draws on insights from dynamical systems descriptions of conflict developed in the social and behavioral sciences to present a model that captures the convoluted, interacting challenges that stifle progress on this problem. This description of complexity explains the persistence of exclusive cultures and the inadequacy of quick or simple fixes. It also motivates the necessity of prolonged and multifaceted approaches to solutions. It is incumbent on our faculties to recognize the complexities in both problem and solutions, and persevere in responding to these intractable dynamics. It is incumbent on our administrations to provide the consistent structure that supports these tasks. It incumbent on all of our constituents - students, administration and faculty - to be cognizant of and responsive to these efforts.

Our multi-view metric learning framework enables robust characterization of star categories by directly learning to discriminate in a multi-faceted feature space, thus, eliminating the need to combine feature representations prior to fitting the machine learning model. We also demonstrate how to extend standard multi-view learning, which employs multiple vectorized views, to the matrix-variate case which allows very novel variable star signature representations. The performance of our proposed methods is evaluated on the UCR Starlight and LINEAR datasets. Both the vector and matrix-variate versions of our multi-view learning framework perform favorably --- demonstrating the ability to discriminate variable star categories.

We present the construction of a novel time-domain signature extraction methodology and the development of a supporting supervised pattern detection algorithm. We focus on the targeted identification of eclipsing binaries that demonstrate a feature known as the O'Connell effect. Our proposed methodology maps stellar variable observations to a new representation known as distribution fields (DFs). Given this novel representation, we develop a metric learning technique directly on the DF space that is capable of specifically identifying our stars of interest. The metric is tuned on a set of labeled eclipsing binary data from the Kepler survey, targeting particular systems exhibiting the O'Connell effect. The result is a conservative selection of 124 potential targets of interest out of the Villanova Eclipsing Binary Catalog. Our framework demonstrates favorable performance on Kepler eclipsing binary data, taking a crucial step in preparing the way for large-scale data volumes from next-generation telescopes such as LSST and SKA.

We study the effects of the local environmental density and the cosmic web environment (filaments, walls, and voids) on key properties of dark matter halos using the Bolshoi-Planck LCDM cosmological simulation. The z = 0 simulation is analysed into filaments, walls, and voids using the SpineWeb method and also the VIDE package of tools, both of which use the watershed transform. The key halo properties that we study are the specific mass accretion rate, spin parameter, concentration, prolateness, scale factor of the last major merger, and scale factor when the halo had half of its z = 0 mass. For all these properties, we find that there is no discernible difference between the halo properties in filaments, walls, or voids when compared at the same environmental density. As a result, we conclude that environmental density is the core attribute that affects these properties. This conclusion is in line with recent findings that properties of galaxies in redshift surveys are independent of their cosmic web environment at the same environmental density at z ~ 0. We also find that the local web environment of the Milky Way and the Andromeda galaxies near the centre of a cosmic wall does not appear to have any effect on the properties of these galaxies' dark matter halos except for their orientation, although we find that it is rather rare to have such massive halos near the centre of a relatively small cosmic wall.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

It is well established that the chemical structure of the Milky Way exhibits a bimodality with respect to the $\alpha$-enhancement of stars at a given [Fe/H]. This has been studied largely based on a bulk $\alpha$ abundance, computed as a summary of several individual $\alpha$-elements. Inspired by the expected subtle differences in their nucleosynthetic origins, here we probe the higher level of granularity encoded in the inter-family [Mg/Si] abundance ratio. Using a large sample of stars with APOGEE abundance measurements, we first demonstrate that there is additional information in this ratio beyond what is already apparent in [$\alpha$/Fe] and [Fe/H] alone. We then consider Gaia astrometry and stellar age estimates to empirically characterize the relationships between [Mg/Si] and various stellar properties. We find small but significant trends between this ratio and $\alpha$-enhancement, age, [Fe/H], location in the Galaxy, and orbital actions. To connect these observed [Mg/Si] variations to a physical origin, we attempt to predict the Mg and Si abundances of stars with the galactic chemical evolution model Chempy. We find that we are unable to reproduce abundances for the stars that we fit, which highlights tensions between the yield tables, the chemical evolution model, and the data. We conclude that a more data-driven approach to nucleosynthetic yield tables and chemical evolution modeling is necessary to maximize insights from large spectroscopic surveys.

Satellite galaxies are predicted to generate gravitational density wakes as they orbit within the dark matter (DM) halos of their hosts, causing their orbits to decay over time. The recent infall of the Milky Way's (MW) most massive satellite galaxy, the Large Magellanic Cloud (LMC), affords us the unique opportunity to study this process in action. In this work, we present high-resolution ($m_{dm} = 4 \times 10^4 M_{\odot}$ ) N-body simulations of the MW-LMC interaction over the past 2 Gyr. We quantify the impact of the LMC's passage on the density and kinematics of the MW's DM halo and the observability of these structures in the MW's stellar halo. The LMC is found to generate pronounced Local and Global wakes in both the DM and stellar halos, leads to both local overdensities and distinct kinematic patterns that should be observable with ongoing and future surveys. Specifically, the Global Wake will result in redshifted radial velocities of stars in the North and blueshifts in the South, at distances larger than 45 kpc. The Local Wake traces the orbital path of the LMC through the halo (50-200 kpc), resulting in a stellar overdensity with a distinct, tangential kinematic pattern that persists to the present day. The detection of the MW's halo response will constrain: the infall mass of the LMC and its orbital trajectory, the mass of the MW, and it may inform us about the nature of the dark matter particle itself.

Orbital properties of stars, computed from their six-dimensional phase space measurements and an assumed Galactic potential, are used to understand the structure and evolution of the Galaxy. Stellar actions, computed from orbits, have the attractive quality of being invariant under certain assumptions and are therefore used as quantitative labels of a star's orbit. We report a subtle but important systematic error that is induced in the actions as a consequence of local midplane variations expected for the Milky Way. This error is difficult to model because it is non-Gaussian and bimodal, with neither mode peaking on the null value. An offset in the vertical position of the Galactic midplane of $\sim15\,\text{pc}$ for a thin disk-like orbit or $\sim 120\,\text{pc}$ for a thick disk-like orbit induces a $25\%$ systematic error in the vertical action $J_z$. In FIRE simulations of Milky Way-mass galaxies, these variations are on the order of $\sim100\,\text{pc}$ at the solar circle. From observations of the mean vertical velocity variation of $\sim5\text{--}10\,\text{km}\,\text{s}^{-1}$ with radius, we estimate that the Milky Way midplane variations are $\sim60\text{--}170\,\text{pc}$, consistent with three-dimensional dust maps. Action calculations and orbit integrations, which assume the global and local midplanes are identical, are likely to include this induced error, depending on the volume considered. Variation in the local standard of rest or distance to the Galactic center causes similar issues. The variation of the midplane must be taken into account when performing dynamical analysis across the large regions of the disk accessible to Gaia and future missions.

We explore the structure of the element abundance--age--orbit distribution of the stars in the Milky Way's low-$\alpha$ disk, by (re-)deriving precise [Fe/H], [X/Fe] and ages, along with orbits, for red clump stars from the APOGEE survey. There has been a long-standing theoretical expectation and observational evidence that metallicity ([Fe/H]) and age are informative about a star's orbit, e.g. about its angular momentum and the corresponding mean Galactocentric distance or its vertical motion. Indeed, our analysis of the APOGEE data confirms that [Fe/H] or age alone can predict the stars' orbits far less well than the combination of the two. Remarkably, we find and show explicitly, that for known [Fe/H] and age, the other abundances [X/Fe] of Galactic disk stars can be predicted well (on average to 0.02 dex) across a wide range of Galactocentric radii, and therefore provide little additional information, e.g. for predicting their orbit. While the age-abundance space for metal poor stars and potentially for stars near the Galactic center is rich or complex, for the bulk of the Galaxy's low-$\alpha$ disk it is simple: [Fe/H] and age contain most information, unless [X/Fe] can be measured to 0.02, or better. Consequently, we do not have the precision with current (and likely near-future) data to assign stars to their individual (coeval) birth clusters, from which the disk is presumably formed. We can, however, place strong constraints on future models of galactic evolution, chemical enrichment and mixing.

The next decade affords tremendous opportunity to achieve the goals of Galactic archaeology. That is, to reconstruct the evolutionary narrative of the Milky Way, based on the empirical data that describes its current morphological, dynamical, temporal and chemical structures. Here, we describe a path to achieving this goal. The critical observational objective is a Galaxy-scale, contiguous, comprehensive mapping of the disk's phase space, tracing where the majority of the stellar mass resides. An ensemble of recent, ongoing, and imminent surveys are working to deliver such a transformative stellar map. Once this empirical description of the dust-obscured disk is assembled, we will no longer be operationally limited by the observational data. The primary and significant challenge within stellar astronomy and Galactic archaeology will then be in fully utilizing these data. We outline the next-decade framework for obtaining and then realizing the potential of the data to chart the Galactic disk via its stars. One way to support the investment in the massive data assemblage will be to establish a Galactic Archaeology Consortium across the ensemble of stellar missions. This would reflect a long-term commitment to build and support a network of personnel in a dedicated effort to aggregate, engineer, and transform stellar measurements into a comprehensive perspective of our Galaxy.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

Discoveries of low mass galaxy pairs and groups are increasing. Studies indicate that dwarf galaxy pairs are gas rich in the field and exhibit elevated star formation rates, suggestive of interactions. Lacking are dynamical models of observed dwarf galaxy pairs to disentangle the physical processes regulating their baryon cycles. We present new optical data and the first detailed theoretical model of an observed tidal encounter between two isolated low mass galaxies, NGC 4490 & NGC 4485. This system is an isolated analog of the Magellanic Clouds and is surrounded by a ~50 kpc extended HI envelope. We use hybrid $N$-body and test-particle simulations along with a visualization interface $Identikit$ to simultaneously reproduce the observed present-day morphology and kinematics. Our results demonstrate how repeated encounters between two dwarf galaxies can "park" baryons at very large distances, without the aid of environmental effects. Our best match to the data is an 8:1 mass ratio encounter where a one-armed spiral is induced in the NGC 4490-analog, which we postulate explains the nature of diffuse starlight presented in the new optical data. We predict that the pair will fully merge in ~370 Myr, but that the extended tidal features will continue to evolve and return to the merged remnant over ~5 Gyr. This pre-processing of baryons will affect the efficiency of gas stripping if such dwarf pairs are accreted by a massive host. In contrast, in isolated environments this study demonstrates how dwarf-dwarf interactions can create a long-lived supply of gas to the merger remnant.

The morphology of thin stellar streams can be used to test the nature of dark matter. It is therefore crucial to extend searches for globular cluster streams to other galaxies than the Milky Way. In this paper, we investigate the current and future prospects of detecting globular cluster streams in external galaxies in resolved stars (e.g. with WFIRST) and using integrated light (e.g. with HSC, LSST and Euclid). In particular, we inject mock-streams to data from the PAndAS M31 survey, and produce simulated M31 backgrounds mimicking what WFIRST will observe in M31. Additionally, we estimate the distance limit to which globular cluster streams will be observable. Our results demonstrate that for a 1 hour (1000 sec.) exposure, using conservative estimates, WFIRST should detect globular cluster streams in resolved stars in galaxies out to distances of ~3.5 Mpc (~2 Mpc). This volume contains 199 (122) galaxies of which >90% are dwarfs. With integrated light, thin streams can be resolved out to ~100 Mpc with HSC and LSST and to ~600 Mpc with WFIRST and Euclid. The low surface brightness of the streams (typically >30 mag/arcsec$^2$), however, will make them difficult to detect, unless the streams originate from very young clusters. We emphasize that if the external galaxies do not host spiral arms or galactic bars, gaps in their stellar streams provide an ideal test case for evidence of interactions with dark matter subhalos. Furthermore, obtaining a large samples of thin stellar streams can help constrain the orbital structure and hence the potentials of external halos.

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

Massive young stellar objects (MYSOs) have recently been shown to drive jets whose particles can interact with either the magnetic fields of the jet or ambient medium to emit non-thermal radiation. We report a search for non-thermal radio emission from a sample of 15 MYSOs to establish the prevalence of the emission in the objects. We used their spectra across the L-, C- and Q-bands along with spectral index maps to characterise their emission. We find that about 50% of the sources show evidence for non-thermal emission with 40% showing clear non-thermal lobes, especially sources of higher bolometric luminosity. The common or IRAS names of the sources that manifest non-thermal lobes are; V645Cyg, IRAS 22134+5834, NGC 7538 IRS 9, IRAS 23262+640, AFGL 402d and AFGL 490. All the central cores of the sources are thermal with corresponding mass-loss rates that lie in the range 3X10^{-7} to 7X10^{-6} solar masses per year. Given the presence of non-thermal lobes in some of the sources and the evidence of non-thermal emission from some spectral index maps, it seems that magnetic fields play a significant role in the jets of massive protostars. Also noted is that some of the sources show evidence of binarity and variability.

Studying our Galaxy, the Milky Way (MW), gives us a close-up view of the interplay between cosmology, dark matter, and galaxy formation. In the next decade our understanding of the MW's dynamics, stellar populations, and structure will undergo a revolution thanks to planned and proposed astrometric, spectroscopic and photometric surveys, building on recent advances by the Gaia astrometric survey. Together, these new efforts will measure three-dimensional positions and velocities and numerous chemical abundances for stars to the MW's edge and well into the Local Group, leading to a complete multidimensional view of our Galaxy. Studies of the multidimensional Milky Way beyond the Gaia frontier---from the edge of the Galactic disk to the edge of our Galaxy's dark matter halo---will unlock new scientific advances across astrophysics, from constraints on dark matter to insights into galaxy formation.

This paper investigates the gravitational trapping of HII regions predicted by steady-state analysis using radiation hydrodynamical simulations. We present idealised spherically symmetric radiation hydrodynamical simulations of the early evolution of HII regions including the gravity of the central source. As with analytic steady state solutions of spherically symmetric ionised Bondi accretion flows, we find gravitationally trapped HII regions with accretion through the ionisation front onto the source. We found that, for a constant ionising luminosity, fluctuations in the ionisation front are unstable. This instability only occurs in this spherically symmetric accretion geometry. In the context of massive star formation, the ionising luminosity increases with time as the source accretes mass. The maximum radius of the recurring HII region increases on the accretion timescale until it reaches the sonic radius, where the infall velocity equals the sound speed of the ionised gas, after which it enters a pressure-driven expansion phase. This expansion prevents accretion of gas through the ionisation front, the accretion rate onto the star decreases to zero, and it stops growing from accretion. Because of the time required for any significant change in stellar mass and luminosity through accretion our simulations keep both mass and luminosity constant and follow the evolution from trapped to expanding in a piecewise manner. Implications of this evolution of HII regions include a continuation of accretion of material onto forming stars for a period after the star starts to emit ionising radiation, and an extension of the lifetime of ultracompact HII regions.

We analyse an N-body simulation of the interaction of the Milky Way (MW) with a Sagittarius-like dSph (Sgr), looking for signatures which may be attributed to its orbital history in the phase space volume around the Sun in light of $Gaia$ DR2 discoveries. The repeated impacts of Sgr excite coupled vertical and radial oscillations in the disc which qualitatively, and to a large degree quantitatively are able to reproduce many features in the 6D $Gaia$ DR2 samples, from the median $V_{R},V_{\phi}, V_{z} $ velocity maps to the local $\delta\rho(v_{z},z)$ phase-space spiral which is a manifestation of the global disc response to coupled oscillations within a given volume. The patterns in the large-scale velocity field are well described by tightly wound spirals and vertical corrugations excited from Sgr's impacts. We show that the last pericentric passage of Sgr resets the formation of the local present-day $\delta\rho(v_z,z)$ spiral and situate its formation around 500-800 Myr. As expected $\delta\rho(vz,z)$ grows in size and decreases in woundedness as a function of radius in both the $Gaia$ DR2 data and simulations. This is the first N-body model able to explain so many of the features in the data on different scales. We demonstrate how to use the full extent of the Galactic disc to date perturbations dating from Myr to Gyr, probe the underlying potential and constrain the mass-loss history of Sgr. $\delta\rho(vz,z)$ looks the same in all stellar populations age bins down to the youngest ages which rules out a bar buckling origin.

We provide homogeneous optical (UBVRI) and near-infrared (JHK) time series photometry for 254 cluster (omega Cen, M4) and field RR Lyrae (RRL) variables. We ended up with more than 551,000 measurements. For 94 fundamental (RRab) and 51 first overtones (RRc) we provide a complete optical/NIR characterization (mean magnitudes, luminosity amplitudes, epoch of the anchor point). The NIR light curves of these variables were adopted to provide new and accurate light-curve templates for both RRc (single period bin) and RRab (three period bins) variables. The templates for the J and the H band are newly introduced, together with the use of the pulsation period to discriminate among the different RRab templates. To overcome subtle uncertainties in the fit of secondary features of the light curves we provide two independent sets of analytical functions (Fourier series, Periodic Gaussian functions). The new templates were validated by using 26 omega Cen and Bulge RRLs covering the four period bins. We found that the difference between the measured mean magnitude along the light curve and the mean magnitude estimated by using the template on a single randomly extracted phase point is better than 0.01 mag (sigma=0.04 mag). We also validated the template on variables for which at least three phase points were available, but without information on the phase of the anchor point. The accuracy of the mean magnitudes is ~0.01 mag (sigma=0.04 mag). The new templates were applied to the Large Magellanic Cloud (LMC) globular Reticulum and by using literature data and predicted PLZ relations we found true distance moduli of 18.47+-0.10+-0.03 mag (J) and 18.49+-0.09+-0.05 mag (K). We also used literature optical and mid-infrared data and we found a mean true distance modulus of 18.47+-0.02+-0.06 mag, suggesting that Reticulum is ~1 kpc closer than the LMC.

Deep observations of faint surface brightness stellar tidal streams in external galaxies with LSST are addressed in this White Paper contribution. We propose using the Wide--Fast--Deep survey that contains several nearby galaxies (at distances where the stars themselves are not resolved, i.e., beyond 20 Mpc). In the context of hierarchical galaxy formation, it is necessary to understand the prevalence and properties of tidal substructure around external galaxies based on integrated (i.e., unresolved) diffuse light. This requires collecting observations on much larger samples of galaxies than the Milky Way and M31. We will compare the observed structures to the predictions of cosmological models of galactic halo formation that inform us about the number and properties of streams around Milky Way-like galaxies. The insight gained from these comparisons will allow us to infer the properties of stream progenitors (masses, dynamics, metallicities, stellar populations). The changes in the host galaxies caused by the interactions with the dissolving companion galaxies will be another focus of our studies. We conclude by discussing synergies with WFIRST and Euclid, and also provide concrete suggestions for how the effects of scattered light could be minimized in LSST images to optimize the search for low surface brightness features, such as faint unresolved stellar tidal streams.

Methods: We observed the high-mass hot core region G351.77-0.54 with ALMA and more than 16km baselines. Results: At a spatial resolution of 18/40au (depending on the distance), we identify twelve sub-structures within the inner few thousand au of the region. The brightness temperatures are high, reaching values greater 1000K, signposting high optical depth toward the peak positions. Core separations vary between sub-100au to several 100 and 1000au. The core separations and approximate masses are largely consistent with thermal Jeans fragmentation of a dense gas core. Due to the high continuum optical depth, most spectral lines are seen in absorption. However, a few exceptional emission lines are found that most likely stem from transitions with excitation conditions above1000K. Toward the main continuum source, these emission lines exhibit a velocity gradient across scales of 100-200au aligned with the molecular outflow and perpendicular to the previously inferred disk orientation. While we cannot exclude that these observational features stem from an inner hot accretion disk, the alignment with the outflow rather suggests that it stems from the inner jet and outflow region. The highest-velocity features are found toward the peak position, and no Hubble-like velocity structure can be identified. Therefore, these data are consistent with steady-state turbulent entrainment of the hot molecular gas via Kelvin-Helmholtz instabilities at the interface between the jet and the outflow. Conclusions: Resolving this high-mass star-forming region at sub-50au scales indicates that the hierarchical fragmentation process in the framework of thermal Jeans fragmentation can continue down to the smallest accessible spatial scales. Velocity gradients on these small scales have to be treated cautiously and do not necessarily stem from disks, but may be better explained with outflow emission.

Tidal debris structures formed from disrupted satellites contain important clues about the assembly histories of galaxies. To date, studies of these structures have been hampered by reliance on by-eye identification and morphological classification which leaves their interpretation significantly uncertain. In this work we present a new machine-vision technique based on the Subspace-Constrained Mean Shift (SCMS) algorithm which can perform these tasks automatically. SCMS finds the location of the high-density `ridges' that define substructure morphology. After identification, the coefficients of an orthogonal series density estimator are used to classify points on the ridges as part of a continuum between shell-like or stream-like debris, from which a global morphological classification can be determined. We dub this procedure Subspace Constrained Unsupervised Detection of Structure (SCUDS). By applying this tool to controlled N--body simulations of minor mergers we demonstrate that the extracted classifications correspond to the well--understood underlying physics of phase mixing. The application of SCUDS to resolved stellar population data from near-future surveys will inform our understanding of the buildup of galaxies stellar halos.

Using a high resolution simulation of an isolated dwarf galaxy, accounting for multi-channel stellar feedback and chemical evolution on a star-by-star basis, we investigate how each of 15 metal species are distributed within our multi-phase interstellar medium (ISM) and ejected from our galaxy by galactic winds. For the first time, we demonstrate that the mass fraction probability distribution functions (PDFs) of individual metal species in the ISM are well described by a piecewise log-normal and power-law distribution. The PDF properties vary within each ISM phase. Hot gas is dominated by recent enrichment, with a significant power-law tail to high metal fractions, while cold gas is predominately log-normal. In addition, elements dominated by asymptotic giant branch (AGB) wind enrichment (e.g. N and Ba) mix less efficiently than elements dominated by supernova enrichment (e.g. $\alpha$ elements and Fe). This result is driven by the differences in source energetics and source locations, particularly the higher chance compared to massive stars for AGB stars to eject material into cold gas. Nearly all of the produced metals are ejected from the galaxy (only 4% are retained), but over 20% of metals dominated by AGB enrichment are retained. In dwarf galaxies, therefore, elements synthesized predominately through AGB winds should be both overabundant and have a larger spread compared to elements synthesized in either core collapse or Type Ia supernovae. We discuss the observational implications of these results, their potential use in developing improved models of galactic chemical evolution, and their generalization to more massive galaxies.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

Radio continuum surveys of the Galactic plane can find and characterize HII regions, supernova remnants (SNRs), planetary nebulae (PNe), and extragalactic sources. A number of surveys at high angular resolution (<25") at different wavelengths exist to study the interstellar medium (ISM), but no comparable high-resolution and high-sensitivity survey exists at long radio wavelengths around 21cm. We observed a large fraction of the Galactic plane in the first quadrant of the Milky Way (l=14.0-67.4deg and |b| < 1.25deg) with the Karl G. Jansky Very Large Array (VLA) in the C-configuration covering six continuum spectral windows. These data provide a detailed view on the compact as well as extended radio emission of our Galaxy and thousands of extragalactic background sources. We used the BLOBCAT software and extracted 10916 sources. After removing spurious source detections caused by the sidelobes of the synthesised beam, we classified 10387 sources as reliable detections. We smoothed the images to a common resolution of 25" and extracted the peak flux density of each source in each spectral window (SPW) to determine the spectral indices $\alpha$ (assuming $I(\nu)\propto\nu^\alpha$). By cross-matching with catalogs of HII regions, SNRs, PNe, and pulsars, we found radio counterparts for 840 HII regions, 52 SNRs, 164 PNe, and 38 pulsars. We found 79 continuum sources that are associated with X-ray sources. We identified 699 ultra-steep spectral sources ($\alpha < -1.3$) that could be high-redshift galaxies. Around 9000 of the sources we extracted are not classified specifically, but based on their spatial and spectral distribution, a large fraction of them is likely to be extragalactic background sources. More than 7750 sources do not have counterparts in the SIMBAD database, and more than 3760 sources do not have counterparts in the NED database.

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

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

We present 1,201 galaxies at $0.05<z<0.45$ that host tidal features, detected from the first $\sim\! 200$ deg$^2$ of imaging from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). All galaxies in the present sample have spectroscopic observations from the Sloan Digital Sky Survey (SDSS) spectroscopic campaigns, generating a sample of 21208 galaxies. Of these galaxies, we identify 214 shell systems and 987 stream systems. For 575 of these systems, we are additionally able to measure the $(g-i)$ colors of the tidal features. We find evidence for star formation in a subset of the streams, with the exception of streams around massive ellipticals, and find that stream host galaxies span the full range of stellar masses in our sample. Galaxies which host shells are predominantly red and massive: we find that observable shells form more frequently around ellipticals than around disc galaxies of the same stellar mass. Although the majority of the shells in our sample are consistent with being formed by minor mergers, $15\% \pm 4.4\%$ of shell host galaxies have $(g-i)$ colors as red as their host galaxy, consistent with being formed by major mergers. These "red shells" are additionally preferentially aligned with the major axis of the host galaxy, as previously predicted from simulations. We suggest that although the bulk of the observable shell population originates from fairly minor mergers, which preferentially form shells that are not aligned with the major axis of the galaxy, major mergers produce a significant number of observable shells.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object.

Stars aligned in thin stream-like features (feathers), with widths of $\delta\sim1-10^{\circ}$ and lengths as large as $\Delta l\sim180^{\circ}$, have been observed towards the Anticenter of our Galaxy and their properties mapped in abundances and phase-space. We study their origin by analysing similar features arising in an N-body simulation of a Galactic disc interacting with a Sagittarius-like dwarf spheroidal galaxy (Sgr). By following the orbits of the particles identified as contributing to feathers backwards in time, we trace their excitation to one of Sgr's previous pericentric passages. These particles initially span a large range of phase-angles but a tight range of radii, suggesting they provide a probe of populations in distinct annuli in the outer Galactic disc. The structures are long lived and persist after multiple passages on timescales of $\sim4 \,\rm{Gyrs}$. On the sky, they exhibit oscillatory motion that can be traced with a single orbit mapped over much of their full length and with amplitudes and gradients similar to those observed. We demonstrate how these properties of feathers may be exploited to measure the potential, its flattening, as well as infer the strength of recent potential perturbations.

We report on the measurement of the $^{7}$Be($n, p$)$^{7}$Li cross section from thermal to approximately 325 keV neutron energy, performed in the high-flux experimental area (EAR2) of the n_TOF facility at CERN. This reaction plays a key role in the lithium yield of the Big Bang Nucleosynthesis (BBN) for standard cosmology. The only two previous time-of-flight measurements performed on this reaction did not cover the energy window of interest for BBN, and showed a large discrepancy between each other. The measurement was performed with a Si-telescope, and a high-purity sample produced by implantation of a $^{7}$Be ion beam at the ISOLDE facility at CERN. While a significantly higher cross section is found at low-energy, relative to current evaluations, in the region of BBN interest the present results are consistent with the values inferred from the time-reversal $^{7}$Li($p, n$)$^{7}$Be reaction, thus yielding only a relatively minor improvement on the so-called Cosmological Lithium Problem (CLiP). The relevance of these results on the near-threshold neutron production in the p+$^{7}$Li reaction is also discussed.

We designed a follow-up program to find the spectroscopic properties of the Hercules-Aquila Cloud (HAC) and test scenarios for its formation. We measured the radial velocities (RVs) of 45 RR Lyrae in the southern portion of the HAC using the facilities at the MDM observatory, producing the first large sample of velocities in the HAC. We found a double-peaked distribution in RVs, skewed slightly to negative velocities. We compared both the morphology of HAC projected onto the plane of the sky and the distribution of velocities in this structure outlined by RR Lyrae and other tracer populations at different distances to N-body simulations. We found that the behaviour is characteristic of an old, well-mixed accretion event with small apo-galactic radius. We cannot yet rule out other formation mechanisms for the HAC. However, if our interpretation is correct, HAC represents just a small portion of a much larger debris structure spread throughout the inner Galaxy whose distinct kinematic structure should be apparent in RV studies along many lines of sight.

Our Galaxy is thought to have undergone an active evolutionary history dominated by star formation, the accretion of cold gas, and, in particular, mergers up to 10 gigayear ago. The stellar halo reveals rich fossil evidence of these interactions in the form of stellar streams, substructures, and chemically distinct stellar components. The impact of dwarf galaxy mergers on the content and morphology of the Galactic disk is still being explored. Recent studies have identified kinematically distinct stellar substructures and moving groups, which may have extragalactic origin. However, there is mounting evidence that stellar overdensities at the outer disk/halo interface could have been caused by the interaction of a dwarf galaxy with the disk. Here we report detailed spectroscopic analysis of 14 stars drawn from two stellar overdensities, each lying about 5 kiloparsecs above and below the Galactic plane - locations suggestive of association with the stellar halo. However, we find that the chemical compositions of these stars are almost identical, both within and between these groups, and closely match the abundance patterns of the Milky Way disk stars. This study hence provides compelling evidence that these stars originate from the disk and the overdensities they are part of were created by tidal interactions of the disk with passing or merging dwarf galaxies.

Context: The importance of magnetic fields at the onset of star formation related to the early fragmentation and collapse processes is largely unexplored today. Aims: We want to understand the magnetic field properties at the earliest evolutionary stages of high-mass star formation. Methods: The Atacama Large Millimeter Array is used at 1.3mm wavelength in full polarization mode to study the polarized emission and by that the magnetic field morphologies and strengths of the high-mass starless region IRDC18310-4. Results: The polarized emission is clearly detected in four sub-cores of the region. In general it shows a smooth distribution, also along elongated cores. Estimating the magnetic field strength via the Davis-Chandrasekhar-Fermi method and following a structure function analysis, we find comparably large magnetic field strengths between ~0.6 and 3.7mG. Comparing the data to spectral line observations, the turbulent-to-magnetic energy ratio is low, indicating that turbulence does not significantly contribute to the stability of the gas clump. A mass-to-flux ratio around the critical value 1.0 - depending on column density - indicates that the region starts to collapse which is consistent with the previous spectral line analysis of the region. Conclusions: While this high-mass region is collapsing and thus at the verge of star formation, the high magnetic field values and the smooth spatial structure indicate that the magnetic field is important for the fragmentation and collapse process. This single case study can only be the starting point for larger sample studies of magnetic fields at the onset of star formation.

Ultrafaint dwarf galaxies (UFDs) are typically assumed to have simple, stellar populations with star formation ending at reionization. Yet as the observations of these galaxies continue to improve, their star formation histories (SFHs) are revealed to be more complicated than previously thought. In this paper, we study how star formation, chemical enrichment, and mixing proceed in small, dark matter halos at early times using a high-resolution, cosmological, hydrodynamical simulation. The goals are to inform the future use of analytic models and to explore observable properties of the simulated halos in the context of UFD data. Specifically, we look at analytic approaches that might inform metal enrichment within and beyond small galaxies in the early Universe. We find that simple assumptions for modeling the extent of supernova-driven winds agree with the simulation on average whereas inhomogeneous mixing and gas flows have a large effect on the spread in simulated stellar metallicities. In the context of the UFDs, this work demonstrates that simulations can form halos with a complex SFH and a large spread in the metallicity distribution function within a few hundred Myr in the early Universe. In particular, bursty and continuous star formation are seen in the simulation and both scenarios have been argued from the data. Spreads in the simulated metallicities, however remain too narrow and too metal-rich when compared to the UFDs. Future work is needed to help reduce these discrepancies and advance our interpretation of the data.

The Monoceros Ring (also known as the Galactic Anticenter Stellar Structure) and A13 are stellar overdensities at estimated heliocentric distances of $d \sim 11$ kpc and 15 kpc observed at low Galactic latitudes towards the anticenter of our Galaxy. While these overdensities were initially thought to be remnants of a tidally-disrupted satellite galaxy, an alternate scenario is that they are composed of stars from the Milky Way (MW) disk kicked out to their current location due to interactions between a satellite galaxy and the disk. To test this scenario, we study the stellar populations of the Monoceros Ring and A13 by measuring the number of RR Lyrae and M giant stars associated with these overdensities. We obtain low-resolution spectroscopy for RR Lyrae stars in the two structures and measure radial velocities to compare with previously measured velocities for M giant stars in the regions of the Monoceros Ring and A13, to assess the fraction of RR Lyrae to M giant stars ($f_{RR:MG}$) in A13 and Mon/GASS. We perform velocity modeling on 153 RR Lyrae stars (116 in the Monoceros Ring and 37 in A13) and find that both structures have very low $f_{RR:MG}$. The results support a scenario in which stars in A13 and Mon/GASS formed in the MW disk. We discuss a possible association between Mon/GASS, A13, and the Triangulum-Andromeda overdensity based on their similar velocity distributions and $f_{RR:MG}$.

High-mass stars form in much richer environments than those associated with isolated low-mass stars, and once they reach a certain mass, produce ionised (HII) regions. The formation of these pockets of ionised gas are unique to the formation of high-mass stars (M $>8$ M$_\odot$), and present an excellent opportunity to study the final stages of accretion, which could include accretion through the HII region itself. This study of the dynamics of the gas on both sides of these ionisation boundaries in very young HII regions aims to quantify the relationship between the HII regions and their immediate environments.We present high-resolution ($\sim$ 0.5$"$) ALMA observations of nine HII regions selected from the Red MSX Source (RMS) survey with compact radio emission and bolometric luminosities greater than 10$^4$ L$_\odot$. We focus on the initial presentation of the data, including initial results from the radio recombination line H29$\alpha$, some complementary molecules, and the 256 GHz continuum emission. Of the six (out of nine) regions with H29$\alpha$ detections, two appear to have cometary morphologies with velocity gradients across them, and two appear more spherical with velocity gradients suggestive of infalling ionised gas. The remaining two were either observed at low resolution or had signals that were too weak to draw robust conclusions. We also present a description of the interactions between the ionised and molecular gas (as traced by CS (J=5-4)), often (but not always) finding theHII region had cleared its immediate vicinity of molecules. Of our sample of nine, the observations of the two clusters expected to have the youngest HII regions (from previous radio observations) are suggestive of having infalling motions in the H29$\alpha$ emission, which could be indicative of late stage accretion onto the stars despite the presence of an HII region.

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

The star formation rate in the Central Molecular Zone (CMZ) is an order of magnitude lower than predicted according to star formation relations that have been calibrated in the disc of our own and nearby galaxies. Understanding how and why star formation appears to be different in this region is crucial if we are to understand the environmental dependence of the star formation process. Here, we present the detection of a sample of high-mass cores in the CMZ's "dust ridge" that have been discovered with the Submillimeter Array as part of the CMZoom survey. These cores range in mass from ~ 50 - 2150 Msun within radii of 0.1 - 0.25 pc. All appear to be young (pre-UCHII), meaning that they are prime candidates for representing the initial conditions of high-mass stars and sub-clusters. We report that at least two of these cores ('c1' and 'e1') contain young, high-mass protostars. We compare all of the detected cores with high-mass cores in the Galactic disc and find that they are broadly similar in terms of their masses and sizes, despite being subjected to external pressures that are several orders of magnitude greater - ~ 10^8 K/cm^3, as opposed to ~ 10^5 K/cm^3. The fact that > 80% of these cores do not show any signs of star-forming activity in such a high-pressure environment leads us to conclude that this is further evidence for an increased critical density threshold for star formation in the CMZ due to turbulence.

Recent Pan-STARRS data show that the leading arm from the globular cluster Palomar 5 (Pal 5) appears shorter than the trailing arm, while simulations of Pal 5 predict similar angular extents. We demonstrate that including the spinning Galactic bar with appropriate pattern speeds in the dynamical modeling of Pal 5 can reproduce the Pan-STARRS data. As the bar sweeps by, some stream stars experience a difference in net torques near pericenter. This leads to the formation of apparent gaps along Pal 5's tidal streams and these gaps grow due to an energy offset from the rest of the stream members. We conclude that only streams orbiting far from the Galactic center or streams on retrograde orbits (with respect to the bar) can be used to unambiguously constrain dark matter subhalo interactions. Additionally, we expect that the Pal 5 leading arm debris should re-appear south of the Pan-STARRS density truncation.

We present N-body simulations of a Sagittarius-like dwarf spheroidal galaxy (Sgr) that follow its orbit about the Milky Way (MW) since its first crossing of the Galaxy's virial radius to the present day. As Sgr orbits around the MW, it excites vertical oscillations, corrugating and flaring the Galactic stellar disc. These responses can be understood by a two-phase picture in which the interaction is first dominated by torques from the wake excited by Sgr in the MW dark halo before transitioning to tides from Sgr's direct impact on the disc at late times. We show for the first time that a massive Sgr model simultaneously reproduces the locations and motions of arc-like over densities, such as the Monoceros Ring and the Triangulum Andromeda stellar clouds, that have been observed at the extremities of the disc, while also satisfying the solar neighbourhood constraints on the vertical structure and streaming motions of the disc. In additional simulations, we include the Large Magellanic Cloud (LMC) self consistently with Sgr. The LMC introduces coupling through constructive and destructive interference, but no new corrugations. In our models, the excitation of the current structure of the outer disk can be traced to interactions as far back as 6-7 Gyr ago (corresponding to $z\leq1$). Given the apparently quiescent accretion history of the MW over this timescale, this places Sgr as the main culprit behind the vertical oscillations of the disc and the last major accretion event for the Galaxy with the capacity to modulate its chemodynamical structure.

The surroundings of massive protostars constitute an accretion disc which has numerically been shown to be subject to fragmentation and responsible for luminous accretion-driven outbursts. Moreover, it is suspected to produce close binary companions which will later strongly influence the star's future evolution in the Hertzsprung-Russel diagram. We present three-dimensional gravitation-radiation-hydrodynamic numerical simulations of 100 Mo pre-stellar cores. We find that accretion discs of young massive stars violently fragment without preventing the (highly variable) accretion of gaseous clumps onto the protostars. While acquiring the characteristics of a nascent low-mass companion, some disc fragments migrate onto the central massive protostar with dynamical properties showing that its final Keplerian orbit is close enough to constitute a close massive proto-binary system, having a young high-mass and a low-mass component. We conclude on the viability of the disc fragmentation channel for the formation of such short-period binaries, and that both processes -close massive binary formation and accretion bursts- may happen at the same time. FU-Orionis-type bursts, such as observed in the young high-mass star S255IR-NIRS3, may not only indicate ongoing disc fragmentation, but also be considered as a tracer for the formation of close massive binaries - progenitors of the subsequent massive spectroscopic binaries - once the high-mass component of the system will enter the main-sequence phase of its evolution. Finally, we investigate the ALMA-observability of the disc fragments.

Deep photometric surveys of the Milky Way have revealed diffuse structures encircling our Galaxy far beyond the "classical" limits of the stellar disk. This paper reviews results from our own and other observational programs, which together suggest that, despite their extreme positions, the stars in these structures were formed in our Galactic disk. Mounting evidence from recent observations and simulations implies kinematic connections between several of these distinct structures. This suggests the existence of collective disk oscillations that can plausibly be traced all the way to asymmetries seen in the stellar velocity distribution around the Sun. There are multiple interesting implications of these findings: they promise new perspectives on the process of disk heating, they provide direct evidence for a stellar halo formation mechanism in addition to the accretion and disruption of satellite galaxies, and, they motivate searches of current and near-future surveys to trace these oscillations across the Galaxy. Such maps could be used as dynamical diagnostics in the emerging field of "Galactoseismology", which promises to model the history of interactions between the Milky Way and its entourage of satellites, as well examine the density of our dark matter halo. As sensitivity to very low surface brightness features around external galaxies increases, many more examples of such disk oscillations will likely be identified. Statistical samples of such features not only encode detailed information about interaction rates and mergers, but also about long sought-after dark matter halo densities and shapes. Models for the Milky Way's own Galactoseismic history will therefore serve as a critical foundation for studying the weak dynamical interactions of galaxies across the universe.

Thanks to modern sky surveys, over twenty stellar streams and overdensity structures have been discovered in the halo of the Milky Way. In this paper, we present an analysis of spectroscopic observations of individual stars from one such structure, "A13", first identified as an overdensity using the M giant catalog from the Two Micron All-Sky Survey. Our spectroscopic observations show that stars identified with A13 have a velocity dispersion of $\lesssim$ 40 $\mathrm{km~s^{-1}}$, implying that it is a genuine coherent structure rather than a chance super-position of random halo stars. From its position on the sky, distance ($\sim$15~kpc heliocentric), and kinematical properties, A13 is likely to be an extension of another low Galactic latitude substructure -- the Galactic Anticenter Stellar Structure (also known as the Monoceros Ring) -- towards smaller Galactic longitude and farther distance. Furthermore, the kinematics of A13 also connect it with another structure in the southern Galactic hemisphere -- the Triangulum-Andromeda overdensity. We discuss these three connected structures within the context of a previously proposed scenario that one or all of these features originate from the disk of the Milky Way.

Simulations of tidal streams show that close encounters with dark matter subhalos induce density gaps and distortions in on-sky path along the streams. Accordingly, observing disrupted streams in the Galactic halo would substantiate the hypothesis that dark matter substructure exists there, while in contrast, observing collimated streams with smoothly varying density profiles would place strong upper limits on the number density and mass spectrum of subhalos. Here, we examine several measures of stream "disruption" and their power to distinguish between halo potentials with and without substructure and with different global shapes. We create and evolve a population of 1280 streams on a range of orbits in the Via Lactea II simulation of a Milky Way-like halo, replete with a full mass range of {\Lambda}CDM subhalos, and compare it to two control stream populations evolved in smooth spherical and smooth triaxial potentials, respectively. We find that the number of gaps observed in a stellar stream is a poor indicator of the halo potential, but that (i) the thinness of the stream on-sky, (ii) the symmetry of the leading and trailing tails, and (iii) the deviation of the tails from a low-order polynomial path on-sky ("path regularity") distinguish between the three potentials more effectively. We find that globular cluster streams on low-eccentricity orbits far from the galactic center (apocentric radius ~ 30-80 kpc) are most powerful in distinguishing between the three potentials. If they exist, such streams will shortly be discoverable and mapped in high dimensions with near-future photometric and spectroscopic surveys.

OTS44 is one of only four free-floating planets known to have a disk. We have previously shown that it is the coolest and least massive known free-floating planet ($\sim$12 M$_{\rm Jup}$) with a substantial disk that is actively accreting. We have obtained Band 6 (233 GHz) ALMA continuum data of this very young disk-bearing object. The data shows a clear unresolved detection of the source. We obtained disk-mass estimates via empirical correlations derived for young, higher-mass, central (substellar) objects. The range of values obtained are between 0.07 and 0.63 M$_{\oplus}$ (dust masses). We compare the properties of this unique disk with those recently reported around higher-mass (brown dwarfs) young objects in order to infer constraints on its mechanism of formation. While extreme assumptions on dust temperature yield disk-mass values that could slightly diverge from the general trends found for more massive brown dwarfs, a range of sensible values provide disk masses compatible with a unique scaling relation between $M_{\rm dust}$ and $M_{*}$ through the substellar domain down to planetary masses.

The Ophiuchus stellar stream is peculiar: (1) its length is short given the age of its constituent stars, and (2) several probable member stars that lie close in both sky position and velocity have dispersions in these dimensions that far exceed those seen within the stream. The stream's proximity to the Galactic center suggests that the bar must have a significant influence on its dynamical history: The triaxiality and time-dependence of the bar may generate chaotic orbits in the vicinity of the stream that can greatly affect its morphology. We explore this hypothesis with models of stream formation along orbits consistent with Ophiuchus' properties in a Milky Way potential model that includes a rotating bar. We find that in all choices for the rotation parameters of the bar, orbits fit to the stream are strongly chaotic. Mock streams generated along these orbits qualitatively match the observed properties of the stream: because of chaos, stars stripped early generally form low-density, high-dispersion "fans" leaving only the most recently disrupted material detectable as a strong over-density. Our models predict that there should be more low-surface-brightness tidal debris than detected so far, likely with a complex phase-space morphology. The existence of or lack of these features around the Ophiuchus stream would provide an interesting constraint on the properties of the Milky Way bar and would help distinguish between formation scenarios for the stream. This is the first time that chaos has been used to explain the properties of a stellar stream and is the first demonstration of the dynamical importance of chaos in the Galactic halo. The existence of long, thin streams around the Milky Way---presumably formed along non- or weakly-chaotic orbits---may represent only a subset of the total population of disrupted satellites.

Aims: We resolve the small-scale structure around the high-mass hot core region G351.77-0.54 to investigate its disk and fragmentation properties. Methods: Using ALMA at 690GHz with baselines exceeding 1.5km, we study the dense gas, dust and outflow emission at an unprecedented spatial resolution of 0.06" (130AU@2.2kpc). Results: Within the inner few 1000AU, G351.77 fragments into at least four cores (brightness temperatures between 58 and 197K). The central structure around the main submm source #1 with a diameter of ~0.5" does not show additional fragmentation. While the CO(6-5) line wing emission shows an outflow lobe in the north-western direction emanating from source #1, the dense gas tracer CH3CN shows a velocity gradient perpendicular to the outflow that is indicative of rotational motions. Absorption profile measurements against the submm source #2 indicate infall rates on the order of 10^{-4} to 10^{-3}M_sun/yr which can be considered as an upper limit of the mean accretion rates. The position-velocity diagrams are consistent with a central rotating disk-like structure embedded in an infalling envelope, but they may also be influenced by the outflow. Using the CH_3CN(37_k-36_k) k-ladder with excitation temperatures up to 1300K, we derive a gas temperature map of source #1 exhibiting temperatures often in excess of 1000K. Brightness temperatures of the submm continuum never exceed 200K. This discrepancy between gas temperatures and submm dust brightness temperatures (in the optically thick limit) indicates that the dust may trace the disk mid-plane whereas the gas could be tracing a hotter gaseous disk surface layer. In addition, we conduct a pixel-by-pixel Toomre gravitational stability analysis of the central rotating structure. The derived high Q values throughout the structure confirm that this central region appears stable against gravitational instability.

We study the fragmentation of the nearest high line-mass filament, the integral shaped filament (ISF, line-mass $\sim$ 400 M$_\odot$ pc$^{-1}$) in the Orion A molecular cloud. We have observed a 1.6 pc long section of the ISF with the Atacama Large Millimetre/submillimeter Array (ALMA) at 3 mm continuum emission, at a resolution of $\sim$3" (1 200 AU). We identify from the region 43 dense cores with masses about a solar mass. 60% of the ALMA cores are protostellar and 40\% are starless. The nearest neighbour separations of the cores do not show a preferred fragmentation scale; the frequency of short separations increases down to 1 200 AU. We apply a two-point correlation analysis on the dense core separations and show that the ALMA cores are significantly grouped at separations below $\sim$17 000 AU and strongly grouped below $\sim$6 000 AU. The protostellar and starless cores are grouped differently: only the starless cores group strongly below $\sim$6 000 AU. In addition, the spatial distribution of the cores indicates periodic grouping of the cores into groups of $\sim$30 000 AU in size, separated by $\sim$50 000 AU. The groups coincide with dust column density peaks detected by Herschel. These results show hierarchical, two-mode fragmentation in which the maternal filament periodically fragments into groups of dense cores. Critically, our results indicate that the fragmentation models for lower line-mass filaments ($\sim$ 16 M$_\odot$ pc$^{-1}$) fail to capture the observed properties of the ISF. We also find that the protostars identified with Spitzer and Herschel in the ISF are grouped at separations below $\sim$17 000 AU. In contrast, young stars with disks do not show significant grouping. This suggests that the grouping of dense cores is partially retained over the protostar lifetime, but not over the lifetime of stars with disks.

We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structures. There is strong molecular line emission covering the entire MM1 region. The kinematic signatures are inconsistent with only Keplerian rotation although we propose that the shift in the emission line centroids within ~1000au of MM1-Main could hint at an underlying compact disc with Keplerian rotation. We cannot however rule out the possibility of an unresolved binary or multiple system. A putative smaller disc could be fed by the large scale spiral `feeding filament' we detect in both gas and dust emission. We also discuss the nature of the now-resolved continuum sources.

The inner few hundred parsecs of the Milky Way, the Central Molecular Zone (CMZ), is our closest laboratory for understanding star formation in the extreme environments (hot, dense, turbulent gas) that once dominated the universe. We present an update on the first large-area survey to expose the sites of star formation across the CMZ at high-resolution in submillimeter wavelengths: the CMZoom survey with the Submillimeter Array (SMA). We identify the locations of dense cores and search for signatures of embedded star formation. CMZoom is a three-year survey in its final year and is mapping out the highest column density regions of the CMZ in dust continuum and a variety of spectral lines around 1.3 mm. CMZoom combines SMA compact and subcompact configurations with single-dish data from BGPS and the APEX telescope, achieving an angular resolution of about 4" (0.2 pc) and good image fidelity up to large spatial scales.

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.

Currently only a small number of Milky Way (MW) stars are known to exist beyond 100 kpc from the Galactic center. Though the distribution of these stars in the outer halo is believed to be sparse, they can provide evidence of more recent accretion events than in the inner halo and help map out the MW's dark matter halo to its virial radius. We have re-examined the outermost regions of 11 existing stellar halo models with two synthetic surveys: one mimicking present-day searches for distant M giants and another mimicking RR Lyrae (RRLe) projections for LSST. Our models suggest that color and proper motion cuts currently used to select M giant candidates for follow-up successfully remove nearly all halo dwarf self-contamination and are useful for focusing observations on distant M giants, of which there are thousands to tens of thousands beyond 100 kpc in our models. We likewise expect that LSST will identify comparable numbers of RRLe at these distances. We demonstrate that several observable properties of both tracers, such as proximity of neighboring stars, proper motions, and distances (for RRLe) could help us separate different accreted structures from one another. We also discuss prospects for using ratios of M giants to RRLe as a proxy for accretion time, which in the future could provide new constraints on the recent accretion history of our Galaxy.

We present the first results from an ongoing survey for multiplicity among the bright stars using the Navy Precision Optical Interferometer (NPOI). We first present a summary of NPOI observations of known multiple systems, including the first detection of the companion of $\beta$ Scuti with precise relative astrometry, to illustrate the instrument's detection sensitivity for binaries at magnitude differences $\Delta$$m$ $\lessapprox$ 3 over the range of angular separation 3 - 860 milliarcseconds (mas). A limiting $\Delta$$m_{700}$ $\sim$ 3.5 is likely for binaries where the component spectral types differ by less than two. Model fits to these data show good agreement with published orbits, and we additionally present a new orbit solution for one of these stars, $\sigma$ Her. We then discuss early results of the survey of bright stars at $\delta$ $\geq$ -20$\deg$. This survey, which complements previous surveys of the bright stars by speckle interferometry, initially emphasizes bright stars of spectral types F0 through K2. We report observations of 41 stars of apparent visual magnitude $m_V$ $\leq$ 4.30, all having been observed on multiple nights. Analysis of these data produces fitted angular separations, position angles, and component magnitude differences for six previously known visual binaries. Three additional systems were examined as possible binaries, but no conclusive detection could be made. No evidence of close stellar companions within our detection limit of $\Delta$$m$ $\approx$ 3 was found for the remaining 32 stars observed; however, uniform-disk angular diameters are reported for 11 of the resolved stars in this last group.

Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio recombination lines as well as the continuum emission from 1 to 2GHz in full polarization over a large part of the first Galactic quadrant. Results: Covering Galactic longitudes from 14.5 to 67.4deg and latitudes between +-1.25deg, we image all of these lines and the continuum at ~20" resolution. These data allow us to study the various components of the interstellar medium (ISM): from the atomic phase, traced by the HI line, to the molecular phase, observed by the OH transitions, to the ionized medium, revealed by the cm continuum and the Halpha radio recombination lines. Furthermore, the polarized continuum emission enables magnetic field studies. In this overview paper, we discuss the survey outline and present the first data release as well as early results from the different datasets. We now release the first half of the survey; the second half will follow later after the ongoing data processing has been completed. The data in fits format (continuum images and line data cubes) can be accessed through the project web-page http://www.mpia.de/thor. Conclusions: The HI/OH/Recombination line survey of the Milky Way (THOR) opens a new window to the different parts of the ISM. It enables detailed studies of molecular cloud formation, conversion of atomic to molecular gas, and feedback from HII regions as well as the magnetic field in the Milky Way. It is highly complementary to other surveys of our Galaxy, and comparing different datasets allows us to address many open questions.

Until recently, there have been few studies of the protostellar evolution of intermediate-mass (IM) stars, which may bridge the low- and high-mass regimes. This paper aims to investigate whether the properties of an IM protostar within the IRAS 05373+2349 embedded cluster are similar to that of low- and/or high-mass protostars. We carried out Very Large Array as well as Combined Array for Research in Millimeter Astronomy continuum and 12CO(J=1-0) observations, which uncover seven radio continuum sources (VLA 1-7). The spectral index of VLA 2, associated with the IM protostar is consistent with an ionised stellar wind or jet. The source VLA 3 is coincident with previously observed H2 emission line objects aligned in the north-south direction (P.A. -20 to -12 deg), which may be either an ionised jet emanating from VLA 2 or (shock-)ionised cavity walls in the large-scale outflow from VLA 2. The position angle between VLA 2 and 3 is slightly misaligned with the large-scale outflow we map at ~5-arcsec resolution in 12CO (P.A. ~30 deg), which in the case of a jet suggests precession. The emission from the mm core associated with VLA 2 is also detected; we estimate its mass to be 12-23 Msun, depending on the contribution from ionised gas. Furthermore, the large-scale outflow has properties intermediate between outflows from low- and high-mass young stars. Therefore, we conclude that the IM protostar within IRAS 05373+2349 is phenomenologically as well as quantitatively intermediate between the low- and high-mass domains.

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 N-body and hydrodynamical simulations of the response of the Milky Way's baryonic disc to the presence of the Large Magellanic Cloud during a first infall scenario. For a fiducial galactic model reproducing the gross properties of the Galaxy, we explore a set of six initial conditions for the LMC of varying mass which all evolve to fit the measured constraints on its current position and velocity with respect to the Galactic Center. We find that the LMC can produce strong disturbances - warping of the stellar and gaseous discs - in the Galaxy, without violating constraints from the phase-space distribution of stars in the Solar Neighbourhood. All models correctly reproduce the phases of the warp and its anti-symmetrical shape about the disc's mid-plane. If the warp is due to the LMC alone, then the largest mass model is favoured ($2.5\times10^{11}\,\rm{M_{\odot}}$). Still, some quantitative discrepancies remain, including deficits in height of $\Delta Z=0.7 \, \rm{kpc}$ at $R=22 \,\rm{kpc}$ and $\Delta Z=0.7\, \rm{kpc}$ at $R=16 \,\rm{kpc}$. This suggests that even higher infall masses for the LMC's halo are allowed by the data. A comparison with the vertical perturbations induced by a heavy Sagittarius dSph model ($10^{11}\,\rm{M_{\odot}}$) suggest that positive interference with the LMC warp is expected at $R=16 \, \rm{kpc}$. We conclude that the vertical structure of the Galactic disc beyond the Solar Neighbourhood may jointly be shaped by its most massive satellites. As such, the current structure of the MW suggests we are seeing the process of disc heating by satellite interactions in action.

Filamentary structures are common morphological features of the cold, molecular interstellar medium (ISM). Recent studies have discovered massive, hundred-parsec-scale filaments that may be connected to the large-scale, Galactic spiral arm structure. Addressing the nature of these Giant Molecular Filaments (GMFs) requires a census of their occurrence and properties. We perform a systematic search of GMFs in the fourth Galactic quadrant and determine their basic physical properties. We identify GMFs based on their dust extinction signatures in near- and mid-infrared and velocity structure probed by ^{13}CO line emission. We use the ^{13}CO line emission and ATLASGAL dust emission data to estimate the total and dense gas masses of the GMFs. We combine our sample with an earlier sample from literature and study the Galactic environment of the GMFs. We identify nine GMFs in the fourth Galactic quadrant; six are located in the Centaurus spiral arm and three in inter-arm regions. Combining this sample with an earlier study using the same identification criteria in the first Galactic quadrant results in 16 GMFs, nine of which are located within spiral arms. The GMFs have sizes of 80-160 pc and ^{13}CO-derived masses between 5-90 x 10^{4} Msun. Their dense gas mass fractions are between 1.5-37%, being higher in the GMFs connected to spiral arms. We also compare the different GMF-identification methods and find that emission and extinction based techniques overlap only partially, highlighting the need to use both to achieve a complete census.

Several long, dynamically cold stellar streams have been observed around the Milky Way Galaxy, presumably formed from the tidal disruption of globular clusters. In integrable potentials---where all orbits are regular---tidal debris phase-mixes close to the orbit of the progenitor system. However, the Milky Way's dark matter halo is expected not to be fully integrable; an appreciable fraction of orbits will be chaotic. This paper examines the influence of chaos on the phase-space morphology of cold tidal streams. Streams even in weakly chaotic regions look very different from those in regular regions. We find that streams can be sensitive to chaos on a much shorter time-scale than any standard prediction (from the Lyapunov or frequency-diffusion times). For example, on a weakly chaotic orbit with a chaotic timescale predicted to be $>$1000 orbital periods ($>$1000 Gyr), the resulting stellar stream is, after just a few 10's of orbits, substantially more diffuse than any formed on a nearby but regular orbit. We find that the enhanced diffusion of the stream stars can be understood by looking at the variance in orbital frequencies of orbit ensembles centered around the parent (progenitor) orbit. Our results suggest that long, cold streams around our Galaxy must exist only on regular (or very nearly regular) orbits; they potentially provide a map of the regular regions of the Milky Way potential. This suggests a promising new direction for the use of tidal streams to constrain the distribution of dark matter around our Galaxy.

This paper presents an example where the morphology of a single stellar stream can be used to rule out a specific galactic potential form without the need for velocity information. We investigate the globular cluster Palomar5 (Pal 5), which is tidally disrupting into a cold, thin stream mapped over 22 degrees on the sky with a typical width of 0.7 degrees. We generate models of this stream by fixing Pal 5's present-day position, distance and radial velocity via observations, while allowing its proper motion to vary. In a spherical dark matter halo we easily find models that fit the observed morphology. However, no plausible Pal 5 model could be found in the triaxial potential of Law & Majewski (2010), which has been proposed to explain the properties of the Sagittarius stream. In this case, the long, thin and curved morphology of the Pal5 stream alone can be used to rule out such a potential configuration. Pal 5 like streams in this potential are either too straight, missing the curvature of the observations, or show an unusual morphology which we dub stream-fanning: a signature sensitive to the triaxiality of a potential. We conclude that the mere existence of other thin tidal streams must provide broad constraints on the orientation and shape of the dark matter halo they inhabit.

The stellar debris structures that have been discovered around the Milky Way and other galaxies are thought to be formed from the disruption of satellite stellar systems --- dwarf galaxies or globular clusters --- by galactic tidal fields. The total stellar mass in these structures is typically tiny compared to the galaxy around which they are found, and it is hence easy to dismiss them as inconsequential. However, they are remarkably useful as probes of a galaxy's history (as described in this chapter) and mass distribution (covered in a companion chapter in this volume). This power is actually a consequence of their apparent insignificance: their low contribution to the overall mass makes the physics that describes them both elegant and simple and this means that their observed properties are relatively easy to understand and interpret.

Tidal debris streams from galaxy satellites can provide insight into the dark matter distribution in halos. This is because we have more information about stars in a debris structure than about a purely random population of stars: we know that in the past they were all bound to the same dwarf galaxy; and we know that they form a dynamically cold population moving on similar orbits. They also probe a different region of the matter distribution in a galaxy than many other methods of mass determination, as their orbits take them far beyond the typical extent of those for the bulk of stars. Although conclusive results from this information have yet to be obtained, significant progress has been made in developing the methodologies for determining both the global mass distribution of the Milky Way's dark matter halo and the amount of dark matter substructure within it. Methods for measuring the halo shape are divided into "predictive methods," which predict the tidal debris properties from the progenitor satellite's mass and orbit, given an assumed parent galaxy mass distribution; and "fundamental methods," which exploit properties fundamental to the nature of tidal debris as global potential constraints. Methods for quantifying the prevalence of dark matter subhalos within halos through the analysis of the gaps left in tidal streams after these substructures pass through them are reviewed.

With the advent of digital astronomy, new benefits and new problems have been presented to the modern day astronomer. While data can be captured in a more efficient and accurate manor using digital means, the efficiency of data retrieval has led to an overload of scientific data for processing and storage. This paper will focus on the construction and application of a supervised pattern classification algorithm for the identification of variable stars. Given the reduction of a survey of stars into a standard feature space, the problem of using prior patterns to identify new observed patterns can be reduced to time tested classification methodologies and algorithms. Such supervised methods, so called because the user trains the algorithms prior to application using patterns with known classes or labels, provide a means to probabilistically determine the estimated class type of new observations. This paper will demonstrate the construction and application of a supervised classification algorithm on variable star data. The classifier is applied to a set of 192,744 LINEAR data points. Of the original samples, 34,451 unique stars were classified with high confidence (high level of probability of being the true class).

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

With the advent of digital astronomy, new benefits and new challenges have been presented to the modern day astronomer. No longer can the astronomer rely on manual processing, instead the profession as a whole has begun to adopt more advanced computational means. This paper focuses on the construction and application of a novel time-domain signature extraction methodology and the development of a supporting supervised pattern classification algorithm for the identification of variable stars. A methodology for the reduction of stellar variable observations (time-domain data) into a novel feature space representation is introduced. The methodology presented will be referred to as Slotted Symbolic Markov Modeling (SSMM) and has a number of advantages which will be demonstrated to be beneficial; specifically to the supervised classification of stellar variables. It will be shown that the methodology outperformed a baseline standard methodology on a standardized set of stellar light curve data. The performance on a set of data derived from the LINEAR dataset will also be shown

The Ophiuchus stellar stream presents a dynamical puzzle: its old stellar populations ($\sim 12$ Gyr) cannot be reconciled with (1) its orbit in a simple model for the Milky Way potential and (2) its short angular extent, both of which imply that the observed stream formed within the last $<1$ Gyr. Recent theoretical work has shown that streams on chaotic orbits may abruptly fan out near their apparent ends; stars in these fans are dispersed in both position and velocity and may be difficult to associate with the stream. Here we present the first evidence of such stream-fanning in the Ophiuchus stream, traced by four blue horizontal branch (BHB) stars beyond the apparent end of the stream. These stars stand out from the background by their high velocities ($v_{\rm los} > 230$ km s$^{-1}$) against $\sim 40$ other stars: their velocities are comparable to those of the stream, but would be exceptional if they were unrelated halo stars. Their positions and velocities are, however, inconsistent with simple extrapolation of the observed cold, high-density portion of the stream. These observations suggest that stream-fanning may be a real, observable effect and, therefore, that Ophiuchus may be on a chaotic orbit. They also show that the Ophiuchus stream is more extended and hence dynamically older than previously thought, easing the stellar population vs. dynamical age tension.

We report on observations of the hydroxyl radical (OH) within The H{\sc I}, OH Recombination line survey (THOR) pilot region. The region is bounded approximately between Galactic coordinates l=29.2 to 31.5$^\circ$ and b=-1.0 to +1.0$^\circ$ and includes the high-mass star forming region W43. We identify 103 maser sites, including 72 with 1612\,MHz masers, 42 showing masers in either of the main line transitions at 1665 and 1667\,MHz and four showing 1720\,MHz masers. Most maser sites with either main-line or 1720\,MHz emission are associated with star formation, whereas most of the 1612\,MHz masers are associated with evolved stars. We find that nearly all of the main-line maser sites are co-spatial with an infrared source, detected by GLIMPSE. We also find diffuse OH emission, as well as OH in absorption towards selected unresolved or partially resolved sites. Extended OH absorption is found towards the well known star forming complex W43 Main.

Aims: We resolve the length-scales for filament formation and fragmentation (res. <=0.1pc), in particular the Jeans length and cylinder fragmentation scale. Methods: We observed the prototypical high-mass star-forming filament IRDC18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2mm continuum and N2H+(1-0) line emission in a ten field mosaic at a spatial resolution of ~4'' (~14000AU). Results: The dust continuum emission resolves the filament into a chain of at least 12 relatively regularly spaced cores. The mean separation between cores is ~0.40(+-0.18)pc. While this is approximately consistent with the fragmentation of an infinite, isothermal, gravitationally bound gas cylinder, a high mass-to-length ratio of M/l~1000M_sun/pc requires additional turbulent and/or magnetic support against radial collapse of the filament. The N2H+(1-0) data reveal a velocity gradient perpendicular to the main filament. Although rotation of the filament cannot be excluded, the data are also consistent with the main filament being comprised of several velocity-coherent sub-filaments. Furthermore, this velocity gradient perpendicular to the filament resembles recent results toward Serpens south that are interpreted as signatures of filament formation within magnetized and turbulent sheet-like structures. Lower-density gas tracers ([CI] and C18O) reveal a similar red/blueshifted velocity structure on scales around 60'' east and west of the IRDC18223 filament. This may tentatively be interpreted as a signature of the large-scale cloud and the smaller-scale filament being kinematically coupled. We do not identify a velocity gradient along the axis of the filament. This may either be due to no significant gas flows along the filamentary axis, but it may partly also be caused by a low inclination angle of the filament with respect to the plane of the sky that could minimize such signature.

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.

We present Atacama Large Millimeter/submillimeter Array (ALMA) line and continuum observations at 1.2mm with ~0.3" resolution that uncover a Keplerian-like disk around the forming O-type star AFGL 4176. The continuum emission from the disk at 1.21 mm (source mm1) has a deconvolved size of 870+/-110 AU x 330+/-300 AU and arises from a structure ~8 M_sun in mass, calculated assuming a dust temperature of 190 K. The first-moment maps, pixel-to-pixel line modeling, assuming local thermodynamic equilibrium (LTE), and position-velocity diagrams of the CH3CN J=13-12 K-line emission all show a velocity gradient along the major axis of the source, coupled with an increase in velocity at small radii, consistent with Keplerian-like rotation. The LTE line modeling shows that where CH3CN J=13-12 is excited, the temperatures in the disk range from ~70 to at least 300 K and that the H2 column density peaks at 2.8x10^24 cm^-2. In addition, we present Atacama Pathfinder Experiment (APEX) 12CO observations which show a large-scale outflow from AFGL 4176 perpendicular to the major axis of mm1, supporting the disk interpretation. Finally, we present a radiative transfer model of a Keplerian disk surrounding an O7 star, with a disk mass and radius of 12 M_sun and 2000 AU, that reproduces the line and continuum data, further supporting our conclusion that our observations have uncovered a Keplerian disk around an O-type star.

How do galaxies move relative to one another? While we can examine the motion of dark matter subhalos around their hosts in simulations of structure formation, determining the orbits of satellites around their parent galaxies from observations is impossible except for a small number of nearby cases. In this work we outline a novel approach to probing the orbital distributions of infalling satellite galaxies using the morphology of tidal debris structures. It has long been understood that the destruction of satellites on near-radial orbits tends to lead to the formation of shells of debris, while those on less eccentric orbits produce tidal streams. We combine an understanding of the scaling relations governing the orbital properties of debris with a simple model of how these orbits phase-mix over time to produce a `morphology metric' that more rigorously quantifies the conditions required for shells to be apparent in debris structures as a function of the satellite's mass and orbit and the interaction time. Using this metric we demonstrate how differences in orbit distributions can alter the relative frequency of shells and stream structures observed around galaxies. These experiments suggest that more detailed modeling and careful comparisons with current and future surveys of low surface brightness features around nearby galaxies should be capable of actually constraining orbital distributions and provide new insights into our understanding of structure formation.

The Triangulum-Andromeda stellar clouds (TriAnd1 and TriAnd2) are a pair of concentric ring- or shell-like over-densities at large $R$ ($\approx$ 30 kpc) and $Z$ ($\approx$ -10 kpc) in the Galactic halo that are thought to have been formed from the accretion and disruption of a satellite galaxy. This paper critically re-examines this formation scenario by comparing the number ratio of RR Lyrae to M giant stars associated with the TriAnd clouds with other structures in the Galaxy. The current data suggest a stellar population for these over-densities ($f_{\rm RR:MG} < 0.38$ at 95% confidence) quite unlike any of the known satellites of the Milky Way ($f_{\rm RR:MG} \approx 0.5$ for the very largest and $f_{\rm RR:MG} >>1$ for the smaller satellites) and more like the population of stars born in the much deeper potential well inhabited by the Galactic disk ($f_{\rm RR:MG} < 0.01$). N-body simulations of a Milky-Way-like galaxy perturbed by the impact of a dwarf galaxy demonstrate that, in the right circumstances, concentric rings propagating outwards from that Galactic disk can plausibly produce similar over-densities. These results provide dramatic support for the recent proposal by Xu et al. (2015) that, rather than stars accreted from other galaxies, the TriAnd clouds could represent stars kicked-out from our own disk. If so, these would be the first populations of disk stars to be found in the Galactic halo and a clear signature of the importance of this second formation mechanism for stellar halos more generally. Moreover, their existence at the very extremities of the disk places strong constraints on the nature of the interaction that formed them.

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

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.

Observational studies of halo stars during the last two decades have placed some limits on the quantity and nature of accreted dwarf galaxy contributions to the Milky Way stellar halo by typically utilizing stellar phase-space information to identify the most recent halo accretion events. In this study we tested the prospects of using 2-D chemical abundance ratio distributions (CARDs) found in stars of the stellar halo to determine its formation history. First, we used simulated data from eleven "MW-like" halos to generate satellite template sets of 2-D CARDs of accreted dwarf satellites which are comprised of accreted dwarfs from various mass regimes and epochs of accretion. Next, we randomly drew samples of $\sim10^{3-4}$ mock observations of stellar chemical abundance ratios ([$\alpha$/Fe], [Fe/H]) from those eleven halos to generate samples of the underlying densities for our CARDs to be compared to our templates in our analysis. Finally, we used the expectation-maximization algorithm to derive accretion histories in relation to the satellite template set (STS) used and the sample size. For certain STS used we typically can identify the relative mass contributions of all accreted satellites to within a factor of 2. We also find that this method is particularly sensitive to older accretion events involving low-luminous dwarfs e.g. ultra-faint dwarfs - precisely those events that are too ancient to be seen by phase-space studies of stars and too faint to be seen by high-z studies of the early Universe. Since our results only exploit two chemical dimensions and near-future surveys promise to provide $\sim6-9$ dimensions, we conclude that these new high-resolution spectroscopic surveys of the stellar halo will allow us to recover its accretion history - and the luminosity function of infalling dwarf galaxies - across cosmic time.

Context: Jets and outflows are key ingredients in the formation of stars across the mass spectrum. In clustered regions, understanding powering sources and outflow components poses a significant problem. Aims: To understand the dynamics in the outflow(s) from a cluster in the process of forming massive stars. Methods: We use new VLA observations of the molecular gas (SiO, CS, OCS and \molec) in the massive star forming region IRAS 17233-3606 which contains a number of HII regions. We compare these observations to previously published molecular data for this source in order to get a holistic view of the outflow dynamics. Results:We find that the dynamics of the various species can be explained by a single large scale ($\sim 0.15$ pc) outflow when compared to the sizes of the HII regions, with the different morphologies of the blue and red outflow components explained with respect to the morphology of the surrounding envelope. We further find that the direction of the velocity gradients seen in OCS and \molec are suggestive of a combination of rotation and outflow motions in the warm gas surrounding the HII regions near the base of the large scale outflow. Conclusions: Our results show that the massive protostars forming within this region appear to be contributing to a single outflow on large scales. This single large scale outflow is traced by a number of different species as the outflow interacts with its surroundings. On the small scales, there appear to be multiple mechanisms contributing to the dynamics which could be a combination of either a small scale outflow or rotation with the dynamics of the large scale outflow.

In contrast to most other galaxies, star-formation rates in the Milky Way can be estimated directly from Young Stellar Objects (YSOs). In the Central Molecular Zone (CMZ) the star-formation rate calculated from the number of YSOs with 24 microns emission is up to order of magnitude higher than the value estimated from methods based on diffuse emission (such as free-free emission). Whether this effect is real or whether it indicates problems with either or both star formation rate measures is not currently known. In this paper, we investigate whether estimates based on YSOs could be heavily contaminated by more evolved objects such as main-sequence stars. We present radiative transfer models of YSOs and of main-sequence stars in a constant ambient medium which show that the main-sequence objects can indeed mimic YSOs at 24 microns. However, we show that in some cases the main-sequence models can be marginally resolved at 24 microns, whereas the YSO models are always unresolved. Based on the fraction of resolved MIPS 24 microns sources in the sample of YSOs previously used to compute the star formation rate, we estimate the fraction of misclassified "YSOs" to be at least 63%, which suggests that the star-formation rate previously determined from YSOs is likely to be at least a factor of three too high.

The massive infrared dark cloud G0.253+0.016 projected 45pc from the Galactic centre contains ~10^5Msun of dense gas whilst being mostly devoid of observed star-formation tracers. To scrutinise the physical properties, dynamics and structure of this cloud with reference to its star-forming potential, we have carried out a concerted SMA and IRAM 30m study of this cloud in dust continuum, CO isotopologues, shock tracing molecules, as well as H$_2$CO to trace the gas temperature. We detect and characterise the dust cores within G0.253+0.016 at ~1.3 mm and find that the kinetic temperature of the gas is >320K on size-scales of ~0.15 pc. Analysis of the position-velocity diagrams of our observed lines show broad linewidths and strong shock emission in the south of the cloud, indicating that G0.253+0.016 is colliding with another cloud at v(LSR)~70 km/s. We confirm via an analysis of the observed dynamics in the CMZ that it is an elongated structure, orientated with Sgr B2 closer to the Sun than Sgr A*, however our results suggest that the actual geometry may be more complex than an elliptical ring. We find that the column density PDF of G0.253+0.016 is log-normal with no discernible power-law tail, consistent with little star formation, and that its width can be explained in the framework of theory predicting the density structure of clouds created by supersonic, magnetised turbulence. We also present the delta-variance spectrum of this region, and show it is consistent with that expected for clouds with no star formation. Using G0.253+0.016 as a test-bed of the conditions required for star formation in a different physical environment to that of nearby clouds, we also conclude that there is not one column density threshold for star formation, but instead this value is dependant on the local physical conditions. [Abbrv.]

ALMA HCO+ observations of the infrared dark cloud G0.253+0.016 located in the Central Molecular Zone of the Galaxy are presented. The 89 GHz emission is area-filling, optically thick, and sub-thermally excited. Two types of filaments are seen in absorption against the HCO+ emission. Broad-line absorption filaments (BLAs) have widths of less than a few arcseconds (0.07 - 0.14 pc), lengths of 30 to 50 arcseconds (1.2 - 1.8 pc), and absorption profiles extending over a velocity range larger than 20 km/sec. The BLAs are nearly parallel to the nearby G0.18 non-thermal filaments and may trace HCO+ molecules gyrating about highly ordered magnetic fields located in front of G0.253+0.016 or edge-on sheets formed behind supersonic shocks propagating orthogonal to our line-of-sight in the foreground. Narrow-line absorption filaments (NLAs) have line-widths less than 20 km/sec. Some NLAs are also seen in absorption in other species with high optical depth such as HCN and occasionally in emission where the background is faint. The NLAs, which also trace low-density, sub-thermally excited HCO+ molecules, are mostly seen on the blueshifted side of the emission from G0.253+0.016. If associated with the surface of G0.253+0.016, the kinematics of the NLAs indicate that the cloud surface is expanding. The decompression of entrained, milli-Gauss magnetic fields may be responsible for the re-expansion of the surface layers of G0.253+0.016 as it recedes from the Galactic center following a close encounter with Sgr A.

As large-scale stellar surveys have become available over the past decade, the ability to detect and characterize substructures in the Galaxy has increased dramatically. These surveys have revealed the Triangulum-Andromeda (TriAnd) region to be rich with substructure in the distance range 20-30 kpc, and the relation of these features to each other -- if any -- remains unclear. This complex situation motivates this re-examination of the TriAnd region with a photometric and spectroscopic survey of M giants. An exploration using 2MASS photometry reveals not only the faint sequence in M giants detected by Rocha-Pinto et al. (2004) spanning the range $100^{\circ}<l<160^{\circ}$ and $-50^{\circ}<b<-15^{\circ}$ but, in addition, a second, brighter and more densely populated M giant sequence. These two sequences are likely associated with the two distinct main-sequences discovered (and labeled TriAnd1 and TriAnd2) by Martin et al. (2007) in an optical survey in the direction of M31, where TriAnd2 is the optical counterpart of the fainter RGB/AGB sequence of Rocha-Pinto et al. (2004). Here, the age, distance, and metallicity ranges for TriAnd1 and TriAnd2 are estimated by simultaneously fitting isochrones to the 2MASS RGB tracks and the optical MS/MSTO features. The two populations are clearly distinct in age and distance: the brighter sequence (TriAnd1) is younger (6-10 Gyr) and closer (distance of $\sim$ 15-21 kpc), while the fainter sequence (TriAnd2) is older (10-12 Gyr) and is at an estimated distance of $\sim$ 24-32 kpc. A comparison with simulations demonstrates that the differences and similarities between TriAnd1 and TriAnd2 can simultaneously be explained if they represent debris originating from the disruption of the same dwarf galaxy, but torn off during two distinct pericentric passages. [Abridged]

Throughout the Milky Way, molecular clouds typically appear filamentary, and mounting evidence indicates that this morphology plays an important role in star formation. What is not known is to what extent the dense filaments most closely associated with star formation are connected to the surrounding diffuse clouds up to arbitrarily large scales. How are these cradles of star formation linked to the Milky Way's spiral structure? Using archival Galactic plane survey data, we have used multiple datasets in search of large-scale, velocity-coherent filaments in the Galactic plane. In this paper, we present our methods employed to identify coherent filamentary structures first in extinction and confirmed using Galactic Ring Survey data. We present a sample of seven Giant Molecular Filaments (GMFs) that have lengths of order $\sim$100 pc, total masses of 10$^4$ - 10$^5$ M$_{\odot}$, and exhibit velocity coherence over their full length. The GMFs we study appear to be inter-arm clouds and may be the Milky Way analogues to spurs observed in nearby spiral galaxies. We find that between 2 and 12% of the total mass (above $\sim$10$^{20}$ cm$^{-2}$) is "dense" (above 10$^{22}$ cm$^{-2}$), where filaments near spiral arms in the Galactic midplane tend to have higher dense gas mass fractions than those further from the arms.

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.

The dark matter halo of the Milky Way is expected to be triaxial and filled with substructure. It is hoped that streams or shells of stars produced by tidal disruption of stellar systems will provide precise measures of the gravitational potential to test these predictions. We develop a method for inferring the Galactic potential with tidal streams based on the idea that the stream stars were once close in phase space. Our method can flexibly adapt to any form for the Galactic potential: it works in phase-space rather than action-space and hence relies neither on our ability to derive actions nor on the integrability of the potential. Our model is probabilistic, with a likelihood function and priors on the parameters. The method can properly account for finite observational uncertainties and missing data dimensions. We test our method on synthetic datasets generated from N-body simulations of satellite disruption in a static, multi-component Milky Way including a triaxial dark matter halo with observational uncertainties chosen to mimic current and near-future surveys of various stars. We find that with just four well-measured stream stars, we can infer properties of a triaxial potential with precisions of order 5-7 percent. Without proper motions we obtain 15 percent constraints on potential parameters and precisions around 25 percent for recovering missing phase-space coordinates. These results are encouraging for the eventual goal of using flexible, time-dependent potential models combined with larger data sets to unravel the detailed shape of the dark matter distribution around the Milky Way.

The Automated Planet Finder (APF) is a facility purpose-built for the discovery and characterization of extrasolar planets through high-cadence Doppler velocimetry of the reflex barycentric accelerations of their host stars. Located atop Mt. Hamilton, the APF facility consists of a 2.4-m telescope and its Levy spectrometer, an optical echelle spectrometer optimized for precision Doppler velocimetry. APF features a fixed format spectral range from 374 nm - 970 nm, and delivers a "Throughput" (resolution * slit width product) of 114,000 arc-seconds, with spectral resolutions up to 150,000. Overall system efficiency (fraction of photons incident on the primary mirror that are detected by the science CCD) on blaze at 560 nm in planet-hunting mode is 15%. First-light tests on the RV standard stars HD 185144 and HD 9407 demonstrate sub-meter per second precision (RMS per observation) held over a 3-month period. This paper reviews the basic features of the telescope, dome, and spectrometer, and gives a brief summary of first-light performance.

Near-future data from ESA's Gaia mission will provide precise, full phase-space information for hundreds of millions of stars out to heliocentric distances of ~10 kpc. This "horizon" for full phase-space measurements is imposed by the Gaia parallax errors degrading to worse than 10%, and could be significantly extended by an accurate distance indicator. Recent work has demonstrated how Spitzer observations of RR Lyrae stars can be used to make distance estimates accurate to 2%, effectively extending the Gaia, precise-data horizon by a factor of ten in distance and a factor of 1000 in volume. This Letter presents one approach to exploit data of such accuracy to measure the Galactic potential using small samples of stars associated with debris from satellite destruction. The method is tested with synthetic observations of 100 stars from the end point of a simulation of satellite destruction: the shape, orientation, and depth of the potential used in the simulation are recovered to within a few percent. The success of this simple test with such a small sample in a single debris stream suggests that constraints from multiple streams could be combined to examine the Galaxy's dark matter halo in even more detail --- a truly unique opportunity that is enabled by the combination of Spitzer and Gaia with our intimate perspective on the Galaxy.

The dynamics of the core of the Sagittarius (Sgr) dwarf spheroidal (dSph) galaxy are explored using high-resolution (R~22,500), H-band, near-infrared spectra of over 1,000 giant stars in the central 3 deg^2 of the system, of which 328 are identified as Sgr members. These data, among some of the earliest observations from the SDSS-III/Apache Point Observatory Galactic Evolution Experiment (APOGEE) and the largest published sample of high resolution Sgr dSph spectra to date, reveal a distinct gradient in the velocity dispersion of Sgr from 11-14 km/s for radii >0.8 degrees from center to a dynamical cold point of 8 km/s in the Sgr center --- a trend differing from that found in previous kinematical analyses of Sgr over larger scales that suggest a more or less flat dispersion profile at these radii. Well-fitting mass models with either cored and cusped dark matter distributions can be found to match the kinematical results, although the cored profile succeeds with significantly more isotropic stellar orbits than required for a cusped profile. It is unlikely that the cold point reflects an unusual mass distribution. The dispersion gradient may arise from variations in the mixture of populations with distinct kinematics within the dSph; this explanation is suggested (e.g., by detection of a metallicity gradient across similar radii), but not confirmed, by the present data. Despite these remaining uncertainties about their interpretation, these early test data (including some from instrument commissioning) demonstrate APOGEE's usefulness for precision dynamical studies, even for fields observed at extreme airmasses.

Tidal streams are a powerful probe of the Milky Way (MW) potential shape. In this paper, we introduce a simple test particle method to fit stream data, using a Markov Chain Monte Carlo technique to marginalise over uncertainties in the progenitor's orbit and the Milky Way halo shape parameters. Applying it to mock data of thin streams in the MW halo, we show that, even for very cold streams, stream-orbit offsets - not modelled in our simple method - introduce systematic biases in the recovered shape parameters. For the streams that we consider, and our particular choice of potential parameterisation, these errors are of order ~20% on the halo flattening parameters. However, larger systematic errors can arise for more general streams and potentials; such offsets need to be correctly modelled in order to obtain an unbiased recovery of the underlying potential. Assessing which of the known Milky Way streams are most constraining, we find NGC 5466 and Pal 5 are the most promising candidates. These form an interesting pair as their orbital planes are both approximately perpendicular to each other and to the disc, giving optimal constraints on the MW halo shape. We show that - while with current data their constraints on potential parameters are poor - good radial velocity data along the Pal 5 stream will provide constraints on qz - the flattening perpendicular to the disc. Furthermore, as discussed in a companion paper, NGC 5466 can provide rather strong constraints on the MW halo shape parameters, if the tentative evidence for a departure from the smooth orbit towards its western edge is confirmed.

We present near-infrared spectroscopy and 1 mm line and continuum observations of a recently identified site of high mass star formation likely to be located in the Central Molecular Zone near Sgr C. Located on the outskirts of the massive evolved HII region associated with Sgr C, the area is characterized by an Extended Green Object measuring ~10" in size (0.4 pc), whose observational characteristics suggest the presence of an embedded massive protostar driving an outflow. Our data confirm that early-stage star formation is taking place on the periphery of the Sgr C HII region, with detections of two protostellar cores and several knots of H2 and Brackett gamma emission alongside a previously detected compact radio source. We calculate the cores' joint mass to be ~10^3 Msun, with column densities of 1-2 x 10^24 cm-2. We show the host molecular cloud to hold ~10^5 Msun of gas and dust with temperatures and column densities favourable for massive star formation to occur, however, there is no evidence of star formation outside of the EGO, indicating that the cloud is predominantly quiescent. Given its mass, density, and temperature, the cloud is comparable to other remarkable non-star-forming clouds such as G0.253 in the Eastern CMZ.

One way to constrain the nature of the high-redshift progenitors of the Milky Way is to look at the low-metallicity stellar populations of the different Galactic components today. For example, high-resolution spectroscopy of very metal poor (VMP) stars demonstrates remarkable agreement between the distribution of [Ti/Fe] in the stellar populations of the Milky Way halo (MW) and ultra-faint dwarf (UFD) galaxies. In contrast, for the neutron capture (nc) abundance ratio distributions [(Sr,Ba)/Fe], the peak of the small UFD sample (6 stars) exhibits a signicant under-abundance relative to the VMP stars in the larger MW halo sample (~ 300 stars). We present a simple scenario that can simultaneously explain these similarities and differences by assuming: (i) that the MW VMP stars were predominately enriched by a prior generation of stars which possessed a higher total mass than the prior generation of stars that enriched the UFD VMP stars; and (ii) a much stronger mass-dependent yield (MDY) for nc-elements than for the (known) MDY for Ti. Simple statistical tests demonstrate that conditions (i) and (ii) are consistent with the observed abundance distributions, albeit without strong constraints on model parameters. A comparison of the broad constraints for these nc-MDY with those derived in the literature seems to rule out Ba production from low-mass SNs and affirms models that primarily generate yields from high-mass SN. Our scenario can be confirmed by a relatively modest (factor of ~ 3-4) increase in the number of high-resolution spectra of VMP stars in UFDs.

Recent observations suggest that abundance pattern differences exist between low metallicity stars in the Milky Way stellar halo and those in the dwarf satellite galaxies. This paper takes a first look at what role the early environment for pre-galactic star formation might have played in shaping these stellar populations. In particular, we consider whether differences in cross-pollution between the progenitors of the stellar halo and the satellites could help to explain the differences in abundance patterns. Using an N-body simulation, we find that the progenitor halos of the main halo are primarily clustered together at z=10 while the progenitors of the satellite galaxies remain on the outskirts of this cluster. Next, analytically modeled supernova-driven winds show that main halo progenitors cross-pollute each other more effectively while satellite galaxy progenitors remain more isolated. Thus, inhomogeneous cross-pollution as a result of different high-z spatial locations of each system's progenitors can help to explain observed differences in abundance patterns today. Conversely, these differences are a signature of the inhomogeneity of metal enrichment at early times.

We report the detection of a wide young hierarchical triple system where the primary has a candidate debris disc. The primary, TYC 5241-986-1 A, is a known Tycho star which we classify as a late-K star with emission in the X-ray, near and far-UV and H\alpha\ suggestive of youth. Its proper motion, photometric distance (65-105 pc) and radial velocity lead us to associate the system with the broadly defined Local Association of young stars but not specifically with any young moving group. The presence of weak lithium absorption and X-ray and calcium H and K emission support an age in the 20 to ~125 Myr range. The secondary is a pair of M4.5+-0.5 dwarfs with near and far UV and H\alpha\ emission separated by approximately 1 arcsec (~65-105 AU projected separation) which lie 145 arcsec (9200-15200 AU) from the primary. The primary has a WISE 22 micron excess and follow-up Herschel observations also detect an excess at 70 micron. The excess emissions are indicative of a 100-175 K debris disc. We also explore the possibility that this excess could be due to a coincident background galaxy and conclude that this is unlikely. Debris discs are extremely rare around stars older than 15 Myr, hence if the excess is caused by a disc this is an extremely novel system.

We simulate the evolution of a 10^9 Msun dark matter halo in a cosmological setting with an adaptive-mesh refinement code as an analogue to local low luminosity dwarf irregular and dwarf spheroidal galaxies. The primary goal of our study is to investigate the roles of reionization and supernova feedback in determining the star formation histories of low mass dwarf galaxies. We include a wide range of physical effects, including metal cooling, molecular hydrogen formation and cooling, photoionization and photodissociation from a metagalactic background, a simple prescription for self-shielding, star formation, and a simple model for supernova driven energetic feedback. We carry out simulations excluding each major effect in turn. We find that reionization is primarily responsible for expelling most of the gas in our simulations, but that supernova feedback is required to disperse the dense, cold gas in the core of the halo. Moreover, we show that the timing of reionization can produce an order of magnitude difference in the final stellar mass of the system. For our full physics run with reionization at z=9, we find a stellar mass of about 10^5 Msun at z=0, and a mass-to-light ratio within the half-light radius of approximately 130 Msun/Lsun, consistent with observed low-luminosity dwarfs. However, the resulting median stellar metallicity is 0.06 Zsun, considerably larger than observed systems. In addition, we find star formation is truncated between redshifts 4 and 7, at odds with the observed late time star formation in isolated dwarf systems but in agreement with Milky Way ultrafaint dwarf spheroidals. We investigate the efficacy of energetic feedback in our simple thermal-energy driven feedback scheme, and suggest that it may still suffer from excessive radiative losses, despite reaching stellar particle masses of about 100 Msun, and a comoving spatial resolution of 11 pc.

[Abridged] We present a medium-resolution spectroscopic survey of late-type giant stars at mid-Galactic latitudes of (30$^{\circ}<|b|<60^{\circ}$), designed to probe the properties of this population to distances of $\sim$9 kpc. Because M giants are generally metal-rich and we have limited contamination from thin disk stars by the latitude selection, most of the stars in the survey are expected to be members of the thick disk ($<$[Fe/H]$>\sim$-0.6) with some contribution from the metal-rich component of the nearby halo. Here we report first results for 1799 stars. The distribution of radial velocity (RV) as a function of l for these stars shows (1) the expected thick disk population and (2) local metal-rich halo stars moving at high speeds relative to the disk, that in some cases form distinct sequences in RV-$l$ space. High-resolution echelle spectra taken for 34 of these "RV outliers" reveal the following patterns across the [Ti/Fe]-[Fe/H] plane: seventeen of the stars have abundances reminiscent of the populations present in dwarf satellites of the Milky Way; eight have abundances coincident with those of the Galactic disk and more metal-rich halo; and nine of the stars fall on the locus defined by the majority of stars in the halo. The chemical abundance trends of the RV outliers suggest that this sample consists predominantly of stars accreted from infalling dwarf galaxies. A smaller fraction of stars in the RV outlier sample may have been formed in the inner Galaxy and subsequently kicked to higher eccentricity orbits, but the sample is not large enough to distinguish conclusively between this interpretation and the alternative that these stars represent the tail of the velocity distribution of the thick disk. Our data do not rule out the possibility that a minority of the sample could have formed from gas {\it in situ} on their current orbits.

We use the Bullock & Johnston suite of simulations to study the density profiles of L*-type galaxy stellar haloes. Observations of the Milky Way and M31 stellar haloes show contrasting results: the Milky Way has a `broken' profile, where the density falls off more rapidly beyond ~ 25 kpc, while M31 has a smooth profile out to 100 kpc with no obvious break. Simulated stellar haloes, built solely by the accretion of dwarf galaxies, also exhibit this behavior: some haloes have breaks, while others don't. The presence or absence of a break in the stellar halo profile can be related to the accretion history of the galaxy. We find that a break radius is strongly related to the build up of stars at apocentres. We relate these findings to observations, and find that the `break' in the Milky Way density profile is likely associated with a relatively early (~ 6-9 Gyr ago) and massive accretion event. In contrast, the absence of a break in the M31 stellar halo profile suggests that its accreted satellites have a wide range of apocentres. Hence, it is likely that M31 has had a much more prolonged accretion history than the Milky Way.

This paper aims to investigate the hypothesis that the embedded luminous star AFGL2591-VLA3 (2.3E5Lsun at 3.33kpc) is forming according to a scaled-up version of a low-mass star formation scenario. We present multi-configuration VLA 3.6cm and 7mm, as well as CARMA C18O and 3mm continuum observations to investigate the morphology and kinematics of the ionized gas, dust, and molecular gas around AFGL2591. We also compare our results to ancillary near-IR images, and model the SED and 2MASS image profiles of AFGL2591 using a dust continuum radiative transfer code. The observed 3.6cm images uncover for the first time that the central powering source AFGL2591-VLA3 has a compact core plus collimated jet morphology, extending 4000AU eastward from the central source with an opening angle of <10deg at this radius. However, at 7mm VLA3 does not show a jet morphology, but instead compact (<500AU) emission, some of which (<0.57 mJy of 2.9mJy) is estimated to be from dust. We determine that the momentum rate of the jet is not sufficient to ionize itself via only shocks, and thus a significant portion of the emission is instead likely created in a photoionized wind. The C18O emission uncovers dense entrained material in the outflow(s) from the young stars in the region. The main features of the SED and 2MASS images of AFGL2591-VLA3 are also reproduced by our model dust geometry of a rotationally flattened envelope with and without a disk. The above results are consistent with a picture of massive star formation similar to that seen for low-mass protostars. However, within its envelope, AFGL2591-VLA3 contains at least four other young stars, constituting a small cluster. Therefore it appears that AFGL2591-VLA3 may be able to source its accreting material from a shared gas reservoir while still exhibiting the phenomena expected during the formation of low-mass stars. (Abridged)

We have assembled a large-area spectroscopic survey of giant stars in the Sagittarius (Sgr) dwarf galaxy core. Using medium resolution (R ~15,000), multifiber spectroscopy we have measured velocities of these stars, which extend up to 12 degrees from the galaxy's center (3.7 core radii or 0.4 times the King limiting radius). From these high quality spectra we identify 1310 Sgr members out of 2296 stars surveyed distributed across 24 different fields across the Sgr core. Additional slit spectra were obtained of stars bridging from the Sgr core to its trailing tail. Our systematic, large area sample shows no evidence for significant rotation, a result at odds with the ~20 km/s rotation required as an explanation for the bifurcation seen in the Sgr tidal stream; the observed small (<= 4 km/s) velocity trend along primarily the major axis is consistent with models of the projected motion of an extended body on the sky with no need for intrinsic rotation. The Sgr core is found to have a flat velocity dispersion (except for a kinematically colder center point) across its surveyed extent and into its tidal tails, a property that matches the velocity dispersion profiles measured for other Milky Way dwarf spheroidal (dSph) galaxies. We comment on the possible significance of this observed kinematical similarity for the dynamical state of the other classical Milky Way dSphs in light of the fact that Sgr is clearly a strongly tidally disrupted system.

We distinguish between Local Group field galaxies which may have passed through the virial volume of the Milky Way, and those which have not, via a statistical compari- son against populations of dark matter haloes in the Via Lactea II (VLII) simulation with known orbital histories. Analysis of VLII provides expectations for this escaped population: they contribute 13 per cent of the galactic population between 300 and 1500 kpc from the Milky Way, and hence we anticipate that about 7 of the 54 known Local Group galaxies in that distance range are likely to be Milky Way escapees. These objects can be of any mass below that of the Milky Way, and they are expected to have positive radial velocities with respect to the Milky Way. Comparison of the radius-velocity distributions of VLII populations and measurements of Local Group galaxies presents a strong likelihood that Tucana, Cetus, NGC3109, SextansA, SextansB, Antlia, NGC6822, Phoenix, LeoT, and NGC185 have passed through the Milky Way. Most of these dwarfs have a lower HI mass fraction than the majority of dwarfs lying at similar distances to either the Milky Way or M31. Indeed, several of these galaxies - especially those with lower masses - contain signatures in their morphology, star formation history and/or gas content indicative of evolution seen in simulations of satellite/parent galactic interactions. Our results offer strong support for scenarios in which dwarfs of different types form a sequence in morphology and gas content, with evolution along the sequence being driven by interaction history.

Stellar streams provide unique probes of galactic potentials, with the longer streams normally providing the cleaner measurements. In this paper, we show an example of a short tidal stream that is particularly sensitive to the shape of the Milky Way's dark matter halo: the globular cluster tidal stream NGC 5466. This stream has an interesting deviation from a smooth orbit at its western edge. We show that such a deviation favours an underlying oblate or triaxial halo (irrespective of plausible variations in the Milky Way disc properties and the specific halo parametrisation chosen); spherical or prolate halo shapes can be excluded at a high confidence level. Therefore, more extensive data sets along the NGC 5466 tidal stream promise strong constraints on the Milky Way halo shape.

(Abridged) This paper presents the first connections made between two local features in velocity-space found in a survey of M giant stars and stellar spatial inhomogeneities on global scales. Comparison to cosmological, chemodynamical stellar halo models confirm that the M giant population is particularly sensitive to rare, recent and massive accretion events. These events can give rise to local observed velocity sequences - a signature of a small fraction of debris from a common progenitor, passing at high velocity through the survey volume, near the pericenters of their eccentric orbits. The majority of the debris is found in much larger structures, whose morphologies are more cloud-like than stream-like and which lie at the orbital apocenters. Adopting this interpretation, the full-space motions represented by the observed velocity features are derived under the assumption that the members within each sequence share a common velocity. Orbit integrations are then used to trace the past and future trajectories of these stars across the sky revealing plausible associations with large, previously-discovered, cloud-like structures. The connections made between nearby velocity structures and these distant clouds represent preliminary steps towards developing coherent maps of such giant debris systems. These maps promise to provide new insights into the origin of debris clouds, new probes of Galactic history and structure, and new constraints on the high-velocity tails of the local dark matter distribution that are essential for interpreting direct detection experiments.

The secular stellar mass-loss causes an amplification of the orbital separation in fragile, common proper motion, binary systems with separations of the order of 1000 A.U. In these systems, companions evolve as two independent coeval stars as they experience negligible mutual tidal interactions or mass transfer. We present models for how post-main sequence mass-loss statistically distorts the frequency distribution of separations in fragile binaries. These models demonstrate the expected increase in orbital seapration resulting from stellar mass-loss, as well as a perturbation of associated orbital parameters. Comparisons between our models and observations resulting from the Luyten survey of wide visual binaries, specifically those containing MS and white-dwarf pairs, demonstrate a good agreement between the calculated and the observed angular separation distribution functions.

We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Way's halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative "close pair distribution" (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $\rm r_{gc} < 20$ kpc.

In order to investigate whether massive stars form similarly to their low-mass counterparts, we have used the standard envelope plus disc geometry successfully applied to low-mass protostars to model the near-IR to sub-millimetre SED and several mid-IR images of the embedded massive star IRAS20126+4104. We have used a Monte Carlo radiative transfer dust code to model the continuum absorption, emission and scattering through two azimuthally symmetric dust geometries, the first consisting of a rotationally flattened envelope with outflow cavities, and the second which also includes a flared accretion disc. Our results show that the envelope plus disc model reproduces the observed SED and images more accurately than the model without a disc, although the latter model more closely reproduces the morphology of the mid-IR emission within a radius of 1.1" or ~1800au. We have put forward several possible causes of this discontinuity, including inner truncation of the disc due to stellar irradiation, or precession of the outflow cavity. Our best fitting envelope plus disc model has a disc radius of 9200 au. We find that it is unlikely that the outer regions of such a disc would be in hydrostatic or centrifugal equilibrium, however we calculate that the temperatures within the disc would keep it stable to fragmentation.

Dynamically cold stellar streams are ideal probes of the gravitational field of the Milky Way. This paper re-examines the question of how such streams might be used to test for the presence of "missing satellites" -the many thousands of dark-matter subhalos with masses 10^5-10^7Msolar which are seen to orbit within Galactic-scale dark-matter halos in simulations of structure formation in LCDM cosmologies. Analytical estimates of the frequency and energy scales of stream encounters indicate that these missing satellites should have a negligible effect on hot debris structures, such as the tails from the Sagittarius dwarf galaxy. However, long cold streams, such as the structure known as GD-1 or those from the globular cluster Palomar 5 (Pal 5) are expected to suffer many tens of direct impacts from missing satellites during their lifetimes. Numerical experiments confirm that these impacts create gaps in the debris' orbital energy distribution, which will evolve into degree- and sub-degree- scale fluctuations in surface density over the age of the debris. Maps of Pal 5's own stream contain surface density fluctuations on these scales. The presence and frequency of these inhomogeneities suggests the existence of a population of missing satellites in numbers predicted in the standard LCDM cosmologies.

We present here a fast code for creating a synthetic survey of the Milky Way. Given one or more color-magnitude bounds, a survey size and geometry, the code returns a catalog of stars in accordance with a given model of the Milky Way. The model can be specified by a set of density distributions or as an N-body realization. We provide fast and efficient algorithms for sampling both types of models. As compared to earlier sampling schemes which generate stars at specified locations along a line of sight, our scheme can generate a continuous and smooth distribution of stars over any given volume. The code is quite general and flexible and can accept input in the form of a star formation rate, age metallicity relation, age velocity dispersion relation and analytic density distribution functions. Theoretical isochrones are then used to generate a catalog of stars and support is available for a wide range of photometric bands. As a concrete example we implement the Besancon Milky Way model for the disc. For the stellar halo we employ the simulated stellar halo N-body models of Bullock & Johnston (2005). In order to sample N-body models, we present a scheme that disperses the stars spawned by an N-body particle, in such a way that the phase space density of the spawned stars is consistent with that of the N-body particles. The code is ideally suited to generating synthetic data sets that mimic near future wide area surveys such as GAIA, LSST and HERMES. As an application we study the prospect of identifying structures in the stellar halo with a simulated GAIA survey. We plan to make the code publicly available at http://galaxia.sourceforge.net.

The ubiquity of substructure in the stellar halo has already been demonstrated by the SDSS and 2MASS and future surveys promise to explore the halo in ever more detail. This paper examines what can be learnt from current and future photometric-databases using group-finding techniques. We compare groups recovered from a sample of M-giants from 2MASS with those found in synthetic surveys of simulated $\Lambda$CDM stellar halos and demonstrate broad consistency. We also find that these recovered groups are likely to represent the majority of high-luminosity ($\log(L/L_{Sun})>6.67$) satellites accreted within the last 10 Gyr and on orbits with apocenters within 100 kpc. However the sensitivity of the M-giant survey to accretion events that were either ancient, from low-luminosity objects or those on radial orbits is limited because of the low number of stars, bias towards high-metallicity stars and the shallow depth. We examine the extent to which these limitations are addressed by current and future surveys, in particular catalogues of main-sequence turn-off (MSTO) stars from SDSS and LSST, and of RR-Lyrae stars from LSST or PanSTARRS. The MSTO surveys are more sensitive to low-luminosity events ($\log(L/L_{Sun}~5$ or less) than the 2MASS M-giant sample, while the RR-Lyrae surveys, with superior depth, are good at detecting events on highly eccentric orbits. When combined we expect these photometric surveys to provide a comprehensive picture of the last 10 Gyr of Galactic accretion. Events older than this are too phase mixed to be discovered and would require additional information such as velocity and chemical abundance.

Four close radio sources in the International Celestial Reference Frame (ICRF) catalog were observed using phase referencing with the VLBA at 43, 23 and 8.6 GHz, and with VERA at 23 GHz over a one year period. The goal was to determine the stability of the radio cores, and to assess structure effects associated with positions in the ICRF. Although the four sources were compact at 8.6 GHz, the VLBA images at 43 GHz with 0.3-mas resolution showed that all were composed of several components. A component in each source was identified as the radio core using some or all of the following emission properties: compactness, spectral index, location at the end of the extended emission region, and stationary in the sky. Over the observing period, the relative positions between the four radio cores were constant to 0.02 mas---the phase referencing positional accuracy obtained at 23 and 43 GHz among the sources---suggesting that once a radio core is identified, it remains stationary in the sky to this accuracy. Other radio components in two of the four sources had detectable motion in the radio jet direction. Comparison of the 23 and 43 GHz VLBA images with the VLBA 8.6 GHz images and the ICRF positions suggests that some ICRF positions are dominated by a moving jet component; hence, they can be displaced up to 0.5 mas from the radio core, and may also reflect the motion of the jet component. Future astrometric efforts to determine a more accurate quasar reference frame at 23 and 43 GHz and from the VLBI2010 project are discussed, and supporting VLBA or EVN observations of ICRF sources at 43 GHz are recommended in order to determine the internal structure of the sources. A future collaboration between the radio (ICRF) and the optical frame of Gaia is discussed.

A density based hierarchical group-finding algorithm is used to identify stellar halo structures in a catalog of M-giants from the Two Micron All Sky Survey (2MASS). The intrinsic brightness of M-giant stars means that this catalog probes deep into the halo where substructures are expected to be abundant and easy to detect. Our analysis reveals 16 structures at high Galactic latitude (greater than 15 degree), of which 10 have been previously identified. Among the six new structures two could plausibly be due to masks applied to the data, one is associated with a strong extinction region and one is probably a part of the Monoceros ring. Another one originates at low latitudes, suggesting some contamination from disk stars, but also shows protrusions extending to high latitudes, implying that it could be a real feature in the stellar halo. The last remaining structure is free from the defects discussed above and hence is very likely a satellite remnant. Although the extinction in the direction of the structure is very low, the structure does match a low temperature feature in the dust maps. While this casts some doubt on its origin, the low temperature feature could plausibly be due to real dust in the structure itself. The angular position and distance of this structure encompass the Pisces overdensity traced by RR Lyraes in Stripe 82 of the Sloan Digital Sky Survey (SDSS). However, the 2MASS M-giants indicate that the structure is much more extended than what is visible with the SDSS, with the point of peak density lying just outside Stripe 82. The morphology of the structure is more like a cloud than a stream and reminiscent of that seen in simulations of satellites disrupting along highly eccentric orbits.

Measuring the full three-dimensional motions of extra-galactic objects in the Universe presents a seemingly insurmountable challenge. In this paper we investigate the application of a technique to measure tangential motion that has previously only been applied nearby within the Local Group of galaxies, to clusters of galaxies far beyond its borders. We show that mapping the mean line-of-sight motion throughout a galaxy cluster could in principle be used to detect the "perspective rotation" induced by the projection of the cluster's tangential motion into the line-of-sight. The signal will be most prominent for clusters of the largest angular extent, most symmetric intrinsic velocity distribution and surveyed with the largest number of pointings possible. We investigate the feasibility of detecting this signal using three different approaches: measuring line-of-sight motions of individual cluster members; taking spectra of intracluster gas; and mapping distortions of the Cosmic Microwave Background radiation. We conclude that future spectroscopic surveys of 1000's of members of nearby galaxy clusters hold the most promise of measuring cluster tangential motions using this technique.

We introduce the TANAMI program (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) which is monitoring an initial sample of 43 extragalactic jets located south of -30 degrees declination at 8.4 GHz and 22 GHz since 2007. All aspects of the program are discussed. First epoch results at 8.4 GHz are presented along with physical parameters derived therefrom. We present first epoch images for 43 sources, some observed for the first time at milliarcsecond resolution. Parameters of these images as well as physical parameters derived from them are also presented and discussed. These and subsequent images from the TANAMI survey are available at http://pulsar.sternwarte.uni-erlangen.de/tanami/ We obtain reliable, high dynamic range images of the southern hemisphere AGN. All the quasars and BL Lac objects in the sample have a single-sided radio morphology. Galaxies are either double-sided, single-sided or irregular. About 28% of the TANAMI sample has been detected by LAT during its first three months of operations. Initial analysis suggests that when galaxies are excluded, sources detected by LAT have larger opening angles than those not detected by LAT. Brightness temperatures of LAT detections and non-detections seem to have similar distributions. The redshift distributions of the TANAMI sample and sub-samples are similar to those seen for the bright gamma-ray AGN seen by LAT and EGRET but none of the sources with a redshift above 1.8 have been detected by LAT.

We present a summary of the observation strategy of TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry), a monitoring program to study the parsec-scale structure and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Australian Long Baseline Array (LBA) and the trans-oceanic antennas Hartebeesthoek, TIGO, and O'Higgins. TANAMI is focusing on extragalactic sources south of -30 degrees declination with observations at 8.4 GHz and 22 GHz every ~2 months at milliarcsecond resolution. The initial TANAMI sample of 43 sources has been defined before the launch of the Fermi Gamma Ray Space Telescope to include the most promising candidates for bright gamma-ray emission to be detected with its Large Area Telescope (LAT). Since November 2008, we have been adding new sources to the sample, which now includes all known radio- and gamma-ray bright AGN of the Southern Hemisphere. The combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN and the nature of relativistic jets.

The TANAMI (Tracking AGN with Austral Milliarcsecond Interferometry) and associated programs provide comprehensive radio monitoring of extragalactic gamma-ray sources south of declination -30 degrees. Joint quasi-simultaneous observations between the Fermi Gamma-ray Space Telescope and ground based observatories allow us to discriminate between competing theoretical blazar emission models. High resolution VLBI observations are the only way to spatially resolve the sub-parsec level emission regions where the high-energy radiation originates. The gap from radio to gamma-ray energies is spanned with near simultaneous data from the Swift satellite and ground based optical observatories. We present early results from the TANAMI program in the context of this panchromatic suite of observations.

TANAMI (Tracking Active Galactic Nuclei with Austral Milliarcsecond Interferometry) is a monitoring program to study the parsec-scale structures and dynamics of relativistic jets in active galactic nuclei (AGN) of the Southern Hemisphere with the Long Baseline Array and associated telescopes. Extragalactic jets south of -30 degrees declination are observed at 8.4 GHz and 22 GHz every two months at milliarcsecond resolution. The initial TANAMI sample is a hybrid radio and gamma-ray selected sample since the combination of VLBI and gamma-ray observations is crucial to understand the broadband emission characteristics of AGN.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We demonstrate that stars beyond the virial radii of galaxies may be generated by the gravitational impulse received by a satellite as it passes through the pericenter of its orbit around its parent. These stars may become energetically unbound (escaped stars), or may travel to further than a few virial radii for longer than a few Gyr, but still remain energetically bound to the system (wandering stars). Larger satellites (10-100% the mass of the parent), and satellites on more radial orbits are responsible for the majority of this ejected population. Wandering stars could be observable on Mpc scales via classical novae, and on 100 Mpc scales via SNIa. The existence of such stars would imply a corresponding population of barely-bound, old, high velocity stars orbiting the Milky Way, generated by the same physical mechanism during the Galaxy's formation epoch. Sizes and properties of these combined populations should place some constraints on the orbits and masses of the progenitor objects from which they came, providing insight into the merging histories of galaxies in general and the Milky Way in particular.

We have conducted a search for ionized gas at 3.6 cm, using the Very Large Array, towards 31 Galactic intermediate- and high-mass clumps detected in previous millimeter continuum observations. In the 10 observed fields, 35 HII regions are identified, of which 20 are newly discovered. Many of the HII regions are multiply peaked indicating the presence of a cluster of massive stars. We find that the ionized gas tends to be associated towards the millimeter clumps; of the 31 millimeter clumps observed, 9 of these appear to be physically related to ionized gas, and a further 6 have ionized gas emission within 1'. For clumps with associated ionized gas, the combined mass of the ionizing massive stars is compared to the clump masses to provide an estimate of the instantaneous star formation efficiency. These values range from a few percent to 25%, and have an average of 7 +/- 8%. We also find a correlation between the clump mass and the mass of the ionizing massive stars within it, which is consistent with a power law. This result is comparable to the prediction of star formation by competitive accretion that a power law relationship exists between the mass of the most massive star in a cluster and the total mass of the remaining stars.

Observations of the lengthy tidal streams produced by the destruction of the Sagittarius dwarf spheroidal (Sgr dSph) are capable of providing strong constraints on the shape of the Galactic gravitational potential. However, previous work, based on modeling different stream properties in axisymmetric Galactic models has yielded conflicting results: while the angular precession of the Sgr leading arm is most consistent with a spherical or slightly oblate halo, the radial velocities of stars in this arm are only reproduced by prolate halo models. We demonstrate that this apparent paradox can be resolved by instead adopting a triaxial potential. Our new Galactic halo model, which simultaneously fits all well-established phase space constraints from the Sgr stream, provides the first conclusive evidence for, and tentative measurement of, triaxiality in an individual dark matter halo. The Milky Way halo within ~ 60 kpc is best characterized by a minor/major axis ratio of the isovelocity contours c/a ~ 0.67, intermediate/major axis ratio b/a ~ 0.83, and triaxiality parameter T ~ 0.56. In this model, the minor axis of the dark halo is coincident with the Galactic X axis connecting the Sun and the Galactic Center to within ~ 15 degrees, while the major axis also lies in the Galactic plane, approximately along the Galactic Y axis.

Extensive photometric and spectroscopic surveys of the Andromeda galaxy (M31) have discovered tidal debris features throughout M31's stellar halo. We present stellar kinematics and metallicities in fields with identified substructure from our on-going SPLASH survey of M31 red giant branch stars with the DEIMOS spectrograph on the Keck II 10-m telescope. Radial velocity criteria are used to isolate members of the kinematically-cold substructures. The substructures are shown to be metal-rich relative to the rest of the dynamically hot stellar population in the fields in which they are found. We calculate the mean metallicity and average surface brightness of the various kinematical components in each field, and show that, on average, higher surface brightness features tend to be more metal-rich than lower surface brightness features. Simulations of stellar halo formation via accretion in a cosmological context are used to illustrate that the observed trend can be explained as a natural consequence of the observed dwarf galaxy mass-metallicity relation. A significant spread in metallicity at a given surface brightness is seen in the data; we show that this is due to time effects, namely the variation in the time since accretion of the tidal streams' progenitor onto the host halo. We show that in this theoretical framework a relationship between the alpha-enhancement and surface brightness of tidal streams is expected, which arises from the varying times of accretion of the progenitor satellites onto the host halo. Thus, measurements of the alpha-enrichment, metallicity, and surface brightness of tidal debris can be used to reconstruct the luminosity and time of accretion onto the host halo of the progenitors of tidal streams.

Photometric observations of 235 extragalactic objects that are potential targets for the Space Interferometry Mission (SIM) are presented. Mean B, V, R, I magnitudes at the 5% level are obtained at 1 - 4 epochs between 2005 and 2007 using the 1-m telescopes at Cerro Tololo Inter-American Observatory and Naval Observatory Flagstaff Station. Of the 134 sources which have V magnitudes in the Veron & Veron-Cetty catalog a difference of over 1.0 mag is found for the observed-catalog magnitudes for about 36% of the common sources, and 10 sources show over 3 mag difference. Our first set of observations presented here form the basis of a long-term photometric variability study of the selected reference frame sources to assist in mission target selection and to support in general QSO multi-color photometric variability studies.

Near-future surveys promise a dramatic improvement in the number and precision of astrometric, photometric and spectroscopic measurements of stars in the Milky Way's disk. We examine the impact of such surveys on our understanding of the Galaxy by "observing" particle realizations of non-axisymmetric disk distributions orbiting in an axisymmetric halo with appropriate errors and then attempting to recover the underlying potential using a Markov Chain Monte Carlo (MCMC) approach. We demonstrate that the azimuthally averaged gravitational force field in the Galactic plane--and hence, to a lesser extent, the Galactic mass distribution--can be tightly constrained over a large range of radii using a variety of types of surveys so long as the error distribution of the measurements of the parallax, proper motion and radial velocity are well-understood and the disk is surveyed globally. One advantage of our method is that the target stars can be selected non-randomly in real or apparent-magnitude space to ensure just such a global sample without biasing the results. Assuming we can always measure the line-of-sight velocity of a star with at least 1 km/s precision, we demonstrate that the force field can be determined to better than ~1% for Galactocentric radii in the range R=4-20 kpc We conclude that near-future surveys, like SIM Lite, Gaia, and VERA, will provide the first precise mapping of the gravitational force field in the region of the Galactic disk.

Abridged: The golden age of astrophysics is upon us with both grand discoveries (extra-solar planets, dark matter, dark energy) and precision cosmology. Fundamental understanding of the working of stars and galaxies is within reach, thanks to newly available precision measurements. We highlight the importance of distances and model independent distances and masses. Distances are fundamental in astrophysics and their knowledge can change our perception of phenomena dramatically: e.g., in antiquity, the Heliocentric model was rejected because the predicted stellar parallaxes were not observed. Distance measurements are directly related to the history & fate of the uni- verse as they provide 2 of 3 methods available to date the universe. The 1st method is based on the ages of stars, which can be ascertained if their lumi- nosities (distances) are accurately known. The 2nd method relies on cosmolo- gical methods. To 1st order, the age of the universe is the inverse H_0. As stressed by the previous decadal report, "the fundamental goal of ... astrophysics is to understand how the universe ... galaxies [and] stars ... formed, how they evolved, and what their destiny will be." These questions can be answered partly by micro-arcsecond astrometry: 1) Galactic archeology: a detailed reconstruction of the formation history of the Milky Way and other Local Group galaxies, 2) the very oldest stars in the Milky Way and the age of the Universe, and 3) H_0 and concordance cosmology. These goals are achievable by combining muas-arcsecond astrometry from the proposed SIM-Lite mission supplemented with ground-based spectroscopy. The results of our proposed project will force the biggest reassessment of stellar astrophysics in 50 years, which will affect most branches of astrophysics.

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

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

APOGEE is a large-scale, NIR, high-resolution (R~20,000) spectroscopic survey of Galactic stars. It is one of the four experiments in SDSS-III. Because APOGEE will observe in the H band, it will be the first survey to pierce through Galactic dust and provide a vast, uniform database of chemical abundances and radial velocities for stars across all Galactic populations (bulge, disk, and halo). The survey will be conducted with a dedicated, 300-fiber, cryogenic, spectrograph that is being built at the University of Virginia, coupled to the ARC 2.5m telescope at Apache Point Observatory. APOGEE will use a significant fraction of the SDSS-III bright time during a three-year period to observe, at high signal-to-noise ratio (S/N>100), about 100,000 giant stars selected directly from 2MASS down to a typical flux limit of H<13. The main scientific objectives of APOGEE are: (1) measuring unbiased metallicity distributions and abundance patterns for the different Galactic stellar populations, (2) studying the processes of star formation, feedback, and chemical mixing in the Milky Way, (3) surveying the dynamics of the bulge and disk, placing constraints on the nature and influence of the Galactic bar and spiral arms, and (4) using extensive chemodynamical data, particularly in the inner Galaxy, to unravel its formation and evolution.

We present a highly reliable flux-limited census of 18,949 point sources in the Galactic mid-plane that have intrinsically red mid-infrared colors. These sources were selected from the Spitzer Space Telescope GLIMPSE I and II surveys of 274 deg^2 of the Galactic mid-plane, and consist mostly of high- and intermediate- mass young stellar objects (YSOs) and asymptotic giant branch (AGB) stars. The selection criteria were carefully chosen to minimize the effects of position- dependent sensitivity, saturation, and confusion. The distribution of sources on the sky and their location in IRAC and MIPS 24 microns color-magnitude and color-color space are presented. Using this large sample, we find that YSOs and AGB stars can be mostly separated by simple color-magnitude selection criteria into approximately 50-70% of YSOs and 30-50% of AGB stars. Planetary nebulae and background galaxies together represent at most 2-3% of all the red sources. 1,004 red sources in the GLIMPSE II region, mostly AGB stars with high mass- loss rates, show significant (>0.3 mag) variability at 4.5 and/or 8.0 microns. With over 11,000 likely YSOs and over 7,000 likely AGB stars, this is to date the largest uniform census of AGB stars and high- and intermediate mass YSOs in the Milky-Way Galaxy.

This paper explores the mapping between the observable properties of a stellar halo in phase- and abundance-space and the parent galaxy's accretion history in terms of the characteristic epoch of accretion and mass and orbits of progenitor objects. The study utilizes a suite of eleven stellar halo models constructed within the context of a standard LCDM cosmology. The results demonstrate that coordinate-space studies are sensitive to the recent (0-8 Gyears ago) merger histories of galaxies (this timescale corresponds to the last few to tens of percent of mass accretion for a Milky-Way-type galaxy). Specifically, the {\it frequency, sky coverage} and {\it fraction of stars} in substructures in the stellar halo as a function of surface brightness are indicators of the importance of recent merging and of the luminosity function of infalling dwarfs. The {\it morphology} of features serves as a guide to the orbital distribution of those dwarfs. Constraints on the earlier merger history (> 8 Gyears ago) can be gleaned from the abundance patterns in halo stars: within our models, dramatic differences in the dominant epoch of accretion or luminosity function of progenitor objects leave clear signatures in the [alpha/Fe] and [Fe/H] distributions of the stellar halo - halos dominated by very early accretion have higher average [alpha/Fe], while those dominated by high luminosity satellites have higher [Fe/H]. This intuition can be applied to reconstruct much about the merger histories of nearby galaxies from current and future data sets.

The discovery of extra-solar planets is one of the greatest achievements of modern astronomy. The detection of planets with a wide range of masses demonstrates that extra-solar planets of low mass exist. In this paper we describe a mission, called Darwin, whose primary goal is the search for, and characterization of, terrestrial extrasolar planets and the search for life. Accomplishing the mission objectives will require collaborative science across disciplines including astrophysics, planetary sciences, chemistry and microbiology. Darwin is designed to detect and perform spectroscopic analysis of rocky planets similar to the Earth at mid-infrared wavelengths (6 - 20 micron), where an advantageous contrast ratio between star and planet occurs. The baseline mission lasts 5 years and consists of approximately 200 individual target stars. Among these, 25 to 50 planetary systems can be studied spectroscopically, searching for gases such as CO2, H2O, CH4 and O3. Many of the key technologies required for the construction of Darwin have already been demonstrated and the remainder are estimated to be mature in the near future. Darwin is a mission that will ignite intense interest in both the research community and the wider public.

NGC 205, a close satellite of the M31 galaxy, is our nearest example of a dwarf elliptical galaxy. Photometric and kinematic observations suggest that NGC 205 is undergoing tidal distortion from its interaction with M31. Despite earlier attempts, the orbit and progenitor properties of NGC 205 are not well known. We perform an optimized search for these unknowns by combining a genetic algorithm with restricted N-body simulations of the interaction. This approach, coupled with photometric and kinematic observations as constraints, allows for an effective exploration of the parameter space. We represent NGC 205 as a static Hernquist potential with embedded massless test particles that serve as tracers of surface brightness. We explore 3 distinct, initially stable configurations of test particles: cold rotating disk, warm rotating disk, and hot, pressure-supported spheroid. Each model reproduces some, but not all, of the observed features of NGC 205, leading us to speculate that a rotating progenitor with substantial pressure support could match all of the observables. Furthermore, plausible combinations of mass and scale length for the pressure-supported spheroid progenitor model reproduce the observed velocity dispersion profile. For all 3 models, orbits that best match the observables place the satellite 11+/-9 kpc behind M31 moving at very large velocities: 300-500 km/s on primarily radial orbits. Given that the observed radial component is only 54 km/s, this implies a large tangential motion for NGC 205, moving from the NW to the SE. These results suggest NGC 205 is not associated with the stellar arc observed to the NE of NGC 205. Furthermore, NGC 205's velocity appears to be near or greater than its escape velocity, signifying that the satellite is likely on its first M31 passage.

We present the largest and most detailed survey to date of the stellar populations in the outskirts of M31 based on the analysis of 14 deep HST/ACS pointings spanning the range 11.5-45.0 kpc. We conduct a quantitative comparison of the resolved stellar populations in these fields and identify several striking trends. The color-magnitude diagrams (CMDs), which reach ~3 magnitudes below the red clump, can be classified into two main categories based on their morphologies. `Stream-like' fields, so named for their similarity to the CMD of the giant stellar stream, are characterized by a red clump that slants bluewards at fainter magnitudes and an extended horizontal branch. On the other hand, `disk-like' fields exhibit rounder red clumps with significant luminosity width, lack an obvious horizontal branch and show evidence for recent star formation (~0.25 - 2.0 Gyr ago). We compare the spatial and line-of-sight distribution of stream-like fields with a recent simulation of the giant stream progenitor orbit and find an excellent agreement. These fields, found across the face of M31, attest to the high degree of pollution caused by this event. Disk-like material resides in the extended disk structure of M31 and is detected out to 44 kpc (projected); the uniform populations in these fields, including the ubiquitous presence of young populations, and the strong rotation reported elsewhere are most consistent with a scenario in which this structure has formed through heating and disruption of the existing thin disk, perhaps due to the impact of the giant stream progenitor. Our comparative analysis sheds new light on the likely composition of two of the ultra-deep pointings formerly presented as pure outer disk and pure halo in the literature.

We have used data from the Sloan Digital Sky Survey (SDSS) Data Release 5 to explore the overall structure and substructure of the stellar halo of the Milky Way using about 4 million color-selected main sequence turn-off stars. We fit oblate and triaxial broken power-law models to the data, and found a `best-fit' oblateness of the stellar halo 0.5<c/a<0.8, and halo stellar masses between Galactocentric radii of 1 and 40kpc of (3.7+/-1.2)x10^8 M_sun. The density profile of the stellar halo is approximately r^{-3}; it is possible that the power law slope is shallower inside 20kpc and steeper outside that radius. Yet, we found that all smooth and symmetric models were very poor fits to the distribution of stellar halo stars because the data exhibit a great deal of spatial substructure. We quantified deviations from a smooth oblate/triaxial model using the RMS of the data around the model profile on scales >~100pc, after accounting for the (known) contribution of Poisson uncertainties. The fractional RMS deviation of the actual stellar distribution from any smooth, parameterized halo model is >~40%: hence, the stellar halo is highly structured. We compared the observations with simulations of galactic stellar halos formed entirely from the accretion of satellites in a cosmological context by analysing the simulations in the same way as the data. While the masses, overall profiles, and degree of substructure in the simulated stellar halos show considerable scatter, the properties and degree of substructure in the Milky Way's halo match well the properties of a `typical' stellar halo built exclusively out of the debris from disrupted satellite galaxies. Our results therefore point towards a picture in which an important fraction of the Milky Way's stellar halo has been accreted from satellite galaxies.

We present the results of an experiment to image the interacting binary star beta Lyrae with data from the Navy Prototype Optical Interferometer (NPOI), using a differential phase technique to correct for the effects of the instrument and atmosphere on the interferometer phases. We take advantage of the fact that the visual primary of beta Lyrae and the visibility calibrator we used are both nearly unresolved and nearly centrally symmetric, and consequently have interferometric phases near zero. We used this property to detect and correct for the effects of the instrument and atmosphere on the phases of beta Lyrae and to obtain differential phases in the channel containing the Halpha emission line. Combining the Halpha-channel phases with information about the line strength, we recovered complex visibilities and imaged the Halpha emission using standard radio interferometry methods. We find that the results from our differential phase technique are consistent with those obtained from a more-standard analysis using squared visibilities (V^2's). Our images show the position of the Halpha emitting regions relative to the continuum photocenter as a function of orbital phase and indicate that the major axis of the orbit is oriented along p.a.=248.8+/-1.7 deg. The orbit is smaller than previously predicted, a discrepancy that can be alleviated if we assume that the system is at a larger distance from us, or that the contribution of the stellar continuum to the Halpha channel is larger than estimated. Finally, we also detected a differential phase signal in the channels containing HeI emission lines at 587.6 and 706.5nm, with orbital behavior different from that of the Halpha, indicating that it originates from a different part of this interacting system.

During the last few years, considerable effort has been directed towards large-scale (>> $1 Billion US) missions to detect and characterize earth-like planets around nearby stars, such as the Terrestrial Planet Finder Interferometer (TPF-I) and Darwin missions. However, technological and budgetary issues as well as shifting science priorities will likely prevent these missions from entering Phase A until the next decade. The secondary eclipse technique using the Spitzer Space Telescope has been used to directly measure the temperature and emission spectrum of extrasolar planets. However, only a small fraction of known extrasolar planets are in transiting orbits. Thus, a simplified nulling interferometer, which produces an artificial eclipse or occultation, and operates in the near- to mid-infrared (e.g. ~ 3 to 8 or 10 microns), can characterize the atmospheres of this much larger sample of the known but non-transiting exoplanets. Many other scientific problems can be addressed with a system like this, including imaging debris disks, active galactic nuclei, and low mass companions around nearby stars. We discuss the rationale for a probe-scale mission in the $600-800 Million range, which we name here as the Small Prototype Planet Finding Interferometer (SPPFI).

We report the results of N-body simulations of disrupting satellites aimed at exploring whether the observed features of dSphs can be accounted for with simple, mass-follows-light (MFL) models including tidal disruption. As a test case, we focus on the Carina dwarf Spheroidal (dSph), which presently has the most extensive data at large radius. We find that previous N-body, MFL simulations of dSphs did not sufficiently explore the parameter space of satellite mass, density and orbital shape to find adequate matches to Galactic dSph systems, whereas with a systematic survey of parameter space we are able to find tidally disrupting, MFL satellite models that rather faithfully reproduce Carina's velocity profile, velocity dispersion profile and projected density distribution over its entire sampled radius. The successful MFL model satellites have very eccentric orbits, currently favored by Cold Dark Matter (CDM) models, and central velocity dispersions that still yield an accurate representation of the bound mass and observed central M/L~40 of Carina, despite inflation of the velocity dispersion outside the dSph core by unbound debris. Our survey of parameter space also allows us to address a number of commonly held misperceptions of tidal disruption and its observable effects on dSph structure and dynamics. The simulations suggest that even modest tidal disruption can have a profound effect on the observed dynamics of dSph stars at large radii. Satellites that are well-described by tidally disrupting MFL models could still be fully compatible with L-CDM if for example they represent a later stage in the evolution of luminous subhalos.

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

Recent observations have revealed that high surface-brightness, metal-rich debris is present over large regions of the Andromeda (M31) stellar halo. We use a set of numerical models to determine whether extended metal-rich debris is expected to exist in galaxy halos formed in a hierarchical LambdaCDM universe. We identify tidal debris in the simulations according to the current best surface brightness detection limits in M31 and demonstrate that bright features in the outer halo regions usually originate from single satellites, with very little contamination from other sources due to the low probability of tidal streams from two overlapping accretion events. In contrast, high-surface brightness features in the inner halo often originate from multiple progenitors. We also compare the age and metallicity distribution of the debris with the well-mixed stellar halos in which they reside. We find that high-surface brightness tidal debris is produced almost exclusively by relatively high mass progenitors (Mstar ~ 10^7- 10^9 Msun) and thus is expected to be of moderate- to high-metallicity. Similarly, in our models the smooth inner halo is expected to be metal-rich as this region has been built up mainly from massive satellites. Our results imply that the stellar populations of substructure observed around external galaxies with current techniques should not resemble old and metal-poor dwarf spheroidal satellites, nor the underlying component of the stellar halo.

We report on interferometric and radial-velocity observations of the double-lined 51-d period binary (A) component of the quadruple pre-main sequence (PMS) system V773 Tau. With these observations we have estimated preliminary visual and physical orbits of the V773 Tau A subsystem. Among other parameters, our orbit model includes an inclination of 66.0 $\pm$ 2.4 deg, and allows us to infer the component dynamical masses and system distance. In particular we find component masses of 1.54 $\pm$ 0.14 and 1.332 $\pm$ 0.097 M$_{\sun}$ for the Aa (primary) and Ab (secondary) components respectively. Our modeling of the subsystem component spectral energy distributions finds temperatures and luminosities consistent with previous studies, and coupled with the component mass estimates allows for comparison with PMS stellar models in the intermediate-mass range. We compare V773 Tau A component properties with several popular solar-composition models for intermediate-mass PMS stars. All models predict masses consistent to within 2-sigma of the dynamically determined values, though some models predict values that are more consistent than others.

The past two Decadal Surveys in Astronomy and Astrophysics recommended the completion of a space-based interferometry mission, known today as SIM PlanetQuest, for its unique ability to detect and characterize nearby rocky planets (Bahcall 1991, McKee & Taylor 2001), as well as contributions to a broad range of problems in astrophysics. Numerous committees of the National Research Council as well as NASA Roadmaps have similarly highlighted SIM as the one technology that offers detection and characterization of rocky planets around nearby stars and which is technically ready. To date, SIM remains the only program with the capability of detecting and confirming rocky planets in the habitable zones of nearby solar-type stars. Moreover, SIM measures masses and three-dimensional orbits of habitable planets around nearby stars (within 25 pc); these are the only stars for which follow-up by other techniques is feasible, such as space-based spectroscopy, ground-based interferometry, and of course TPF.

We present a wide-field (4.5 deg^2) photometric and spectroscopic survey of the Leo I dwarf spheroidal (dSph) galaxy to explore its extended morphology and dynamics. As in previous papers in this series, we take advantage of photometry in the M, T_2, and DDO51 filter system to select LeoI red giant branch star candidates, and, so far, this selection technique has proven 100% reliable in selecting actual Leo I members among more than 100 M < 21.5 Leo I giant candidates having previous or new Keck DEIMOS spectroscopy to a radius >1.3 times the limiting radius of the fitted, central King profile. The two-dimensional distribution of all similarly-selected Leo I giant candidates is well fitted by a central single-component King profile of limiting radius 13.3 arcmin, but many giant stars are found outside this newly derived King limiting radius. The density profile thus shows a break at a major axis radial distance of ~10 arcmin produced by an excess of stars at and beyond the King limiting radius (spectroscopically confirmed to be made of true Leo I members), and primarily along the major axis of the main body of the rather elongated satellite. This spatial configuration, a rather flat velocity dispersion profile and an asymmetric radial velocity (RV) distribution among the Leo I members at large radii together support a picture where Leo I has been tidally disrupted on at least one, but at most two, perigalactic passages of a massive Local Group member. (abridged)

We present deep HST/ACS observations of the G1 clump, a distinct stellar overdensity lying at ~30 kpc along the south-western major axis of M31 close to the G1 globular cluster (Ferguson et al. 2002). Our well-populated colour-magnitude diagram reaches ~7 magnitudes below the red giant branch tip with 90% completeness, and allows the detection of various morphological features which can be used to derive detailed constraints on the age and metallicity of the constituent stellar population. We find that the colour-magnitude diagram is best described by a population with a large age range (>= 10 Gyr) and a relatively high mean metallicity [M/H]= -0.4. The spread in metallicity is constrained to be <=0.5 dex. The star formation rate in this region has declined over time, with the bulk of the stellar mass having formed >6 Gyr ago. Nonetheless, a non-negligible mass fraction (approximately 10%) of the population has formed in the last 2 Gyr. We discuss the nature of the G1 Clump in light of these new stellar population constraints and argue that the combination of stellar content and physical size make it unlikely that the structure is the remnant of an accreted dwarf galaxy. Instead, the strong similarity between the stellar content of the G1 Clump and that of the M31 outer disk suggests the substructure is a fragment of the outer disk, perhaps torn off from the main body during a past accretion/merger event; this interpretation is consistent with extant kinematical data. If this interpretation is correct, our analysis of the stellar content provides further evidence that the outskirts of large disk galaxies have been in place for a significant time.

We have used the Very Large Array (VLA), linked with the Pie Town Very Long Baseline Array antenna, to determine astrometric positions of 46 radio stars in the International Celestial Reference Frame (ICRF). Positions were obtained in the ICRF directly through phase referencing of the stars to nearby ICRF quasars whose positions are accurate at the 0.25 mas level. Radio star positions are estimated to be accurate at the 10 mas level, with position errors approaching a few milli-arcseconds for some of the stars observed. Our measured positions were combined with previous measurements taken from as early as 1978 to obtain proper motion estimates for all 46 stars with average uncertainties of ~1.7 mas/yr. We compared our radio star positions and proper motions with the Hipparcos Catalogue data, and find consistency in the reference frames produced by each data set on the 1-sigma level, with errors of ~2.7 mas per axis for the reference frame orientation angles at our mean epoch of 2003.78. No significant spin is found between our radio data frame and the Hipparcos Celestial Reference Frame (HCRF) with largest rotation rates of +0.55 and -0.41 mas/yr around the x and z axes, respectively, with 1-sigma errors of 0.36 mas/yr. Thus, our results are consistent with a non-rotating Hipparcos frame with respect to the ICRF.

Observations of the Mercury-Manganese star Phi Herculis with the Navy Prototype Optical Interferometer (NPOI) conclusively reveal the previously unseen companion in this single-lined binary system. The NPOI data were used to predict a spectral type of A8V for the secondary star Phi Her B. This prediction was subsequently confirmed by spectroscopic observations obtained at the Dominion Astrophysical Observatory. Phi Her B is rotating at 50 +/-3 km/sec, in contrast to the 8 km/sec lines of Phi Her A. Recognizing the lines from the secondary permits one to separate them from those of the primary. The abundance analysis of Phi Her A shows an abundance pattern similar to those of other HgMn stars with Al being very underabundant and Sc, Cr, Mn, Zn, Ga, Sr, Y, Zr, Ba, Ce, and Hg being very overabundant.

A new large-area Washington M,T_2+DDO51 filter survey of more than 10 deg^2 around the Carina dSph galaxy reveals a spectroscopically confirmed power law radial density "break" population of Carina giant stars extending several degrees beyond the central King profile. Magellan telescope MIKE spectroscopy establishes the existence of Carina stars to at least 4.5 times its central King limiting radius, r_lim and primarily along Carina's major axis. To keep these stars bound to the dSph would require a global Carina mass-to-light ratio of M/L > 6,300 M/L_sun. The MIKE velocities, supplemented with ~950 additional Carina field velocities from archived VLT+GIRAFFE spectra with r<=r_lim, demonstrate a nearly constant Carina velocity dispersion to just beyond r = r_lim, and both a rising velocity dispersion and a velocity shear at still larger radii. Together, the observational evidence suggests that the discovered extended Carina population represents tidal debris from the dSph. Of 65 giant candidates at large angular radii from the Carina center for which MIKE spectra have been obtained 94% are associated either with Carina or a second, newly discovered diffuse, but strongly radial velocity-coherent (velocity dispersion of 9.8 km s^-1), foreground halo system. The fifteen stars in this second, retrograde velocity population have (1) a mean metallicity ~1 dex higher than that of Carina, and (2) colors and magnitudes consistent with the red clump of the Large Magellanic Cloud (LMC). Additional spectroscopy of giant star candidates in fields linking Carina and the LMC shows a smooth velocity gradient between the LMC and the retrograde Carina moving group. We conclude that we have found Magellanic stars almost twice as far (22 deg) from the LMC center than previously known.

Large galaxies are thought to form hierarchically, from the accretion and disruption of many smaller galaxies. Such a scenario should naturally lead to galactic phase-space distributions containing some degree of substructure. We examine the errors in mass estimates of galaxies and their dark halos made using the projected phase-space distribution of a tracer population (such as a globular cluster system or planetary nebulae) due to falsely assuming that the tracers are distributed randomly. The level of this uncertainty is assessed by applying a standard mass estimator to samples drawn from 11 random realizations of galaxy halos containing levels of substructure consistent with current models of structure formation. We find that substructure will distort our mass estimates by up to ~20% - a negligible error compared to statistical and measurement errors in current derivations of masses for our own and other galaxies. However, this represents a fundamental limit to the accuracy of any future mass estimates made under the assumption that the tracer population is distributed randomly, regardless of the size of the sample or the accuracy of the measurements.

Vega, the second brightest star in the northern hemisphere, serves as a primary spectral type standard. While its spectrum is dominated by broad hydrogen lines, the narrower lines of the heavy elements suggested slow to moderate rotation, giving confidence that the ground-based calibration of its visibile spectrum could be safely extrapolated into the ultraviolet and near-infrared (through atmosphere models), where it also serves as the primary photometric calibrator. But there have been problems: the star is too bright compared to its peers and it has unusually shaped absorption line profiles, leading some to suggest that it is a distorted, rapidly rotating star seen pole-on. Here we report optical interferometric observations of Vega which detect the asymmetric brightness distribution of the bright, slightly offset polar axis of a star rotating at 93% of breakup speed. In addition to explaining the unusual brightness and line shape pecularities, this result leads to the prediction of an excess of near-infrared emission compared to the visible, in agreement with observations. The large temperature differences predicted across its surface call into question composition determinations, adding uncertainty to Vega's age and opening the possibility that its debris disk could be substantially older than previously thought.

[Abridged] Motivated by upcoming data from astrometric and spectroscopic surveys of the Galaxy, we explore the chemical abundance properties and phase-space distributions in hierarchically-formed stellar halo simulations set in a LambdaCDM Universe. Our sample of Milky-Way type stellar halo simulations result in average metallicities that range from [Fe/H] = -1.3 to -0.9, with the most metal poor halos resulting from accretion histories that lack destructive mergers with massive (metal rich) satellites. Our stellar halo metallicities increase with stellar halo mass. The slope of the [Fe/H]-stellar mass trend mimics that of the satellite galaxies that were destroyed to build the halos, implying that the relation propagates hierarchically. All simulated halos contain a significant fraction of old stellar populations accreted more than 10 Gyr ago and in a few cases, some intermediate age populations exist. In contrast with the Milky Way, many of our simulated stellar halos contain old stellar populations which are metal rich, originating in the early accretion of massive satellites. We suggest that the (metal rich) stellar halo of M31 falls into this category, while the more metal poor halo of the Milky Way is lacking in early massive accretion events. Interestingly, our hierarchically-formed stellar halos often have non-negligible metallicity gradients in both [Fe/H] and [alpha/Fe]. These gradients extend a few tens of kpc, and can be as large as 0.5 dex in [Fe/H] and 0.2 dex in [alpha/Fe]. Finally, we find that chemical abundances can act as a rough substitute for time of accretion of satellite galaxies. We propose a criterion for identifying tidal streams spatially by selecting stars with [alpha/Fe] ratios below solar.

(Abridged) The hierarchical formation scenario for the stellar halo requires the accretion and disruption of dwarf galaxies, yet low-metallicity halo stars are enriched in alpha-elements compared to similar, low-metallicity stars in dwarf spheroidal (dSph) galaxies. We address this primary challenge for the hierarchical formation scenario for the stellar halo by combining chemical evolution modelling with cosmologically-motivated mass accretion histories for the Milky Way dark halo and its satellites. We demonstrate that stellar halo and dwarf galaxy abundance patterns can be explained naturally within the LCDM framework. Our solution relies fundamentally on the LCDM model prediction that the majority of the stars in the stellar halo were formed within a few relatively massive, ~5 x 10^10 Msun, dwarf irregular (dIrr)-size dark matter halos, which were accreted and destroyed ~10 Gyr in the past. These systems necessarily have short-lived, rapid star formation histories, are enriched primarily by Type II supernovae, and host stars with enhanced [a/Fe] abundances. In contrast, dwarf spheroidal galaxies exist within low-mass dark matter hosts of ~10^9 Msun, where supernovae winds are important in setting the intermediate [a/Fe] ratios observed. Our model includes enrichment from Type Ia and Type II supernovae as well as stellar winds, and includes a physically-motivated supernovae feedback prescription calibrated to reproduce the local dwarf galaxy stellar mass - metallicity relation. We use representative examples of the type of dark matter halos we expect to host a destroyed ``stellar halo progenitor'' dwarf, a surviving dIrr, and a surviving dSph galaxy, and show that their derived abundance patterns, stellar masses, and gas masses are consistent with those observed for each type of system.

We report successful fitting of a Roche model, with a surface temperature gradient following the von Zeipel gravity darkening law, to observations of Altair made with the Navy Prototype Optical Interferometer. We confirm the claim by Ohishi, Nordgren, & Hutter that Altair displays an asymmetric intensity distribution due to rotation, the first such detection in an isolated star. Instrumental effects due to the high visible flux of this first magnitude star appear to be the limiting factor in the accuracy of this fit, which nevertheless indicates that Altair is rotating at 0.90+/-0.02 of its breakup (angular) velocity. Our results are consistent with the apparent oblateness found by van Belle et al. and show that the true oblateness is significantly larger owing to an inclination of the rotational axis of ~64 degrees to the line of sight. Of particular interest, we conclude that instead of being substantially evolved as indicated by its classification, A7 VI-V, Altair is only barely off the ZAMS and represents a good example of the difficulties rotation can introduce in the interpretation of this part of the HR diagram.

We present a cosmologically motivated model for the hierachical formation of the stellar halo which includes a semi-analytic treatment of galactic chemical enrichment coupled to numerical simulations that track the orbital evolution and tidal disruption of satellites. A major motivating factor in this investigation is the observed systematic difference between the chemical abundances of stars in satellite galaxies and those in the Milky Way halo. Specifically, for the same [Fe/H] values, stars in neighboring satellite galaxies display significantly lower [alpha/Fe] ratios than stars in the halo. We find that the observed chemical abundance patterns are a natural outcome of the process of hierarchical assembly of the Galaxy. This result follows because the stellar halo in this context is built up from satellite galaxies accreted early on (more than 8-9 Gyr ago) and enriched in alpha-elements produced in Type II supernovae. In contrast, satellites which survive today are typically accreted late (within the last 4-5 Gyr) with nearly solar [alpha/Fe] values as a result of contributions from both Type II and Type Ia supernovae. We use our model to investigate the abundance distribution functions (using both [Fe/H] and [alpha/Fe] ratios) for stars in the halo and in surviving satellites. Our results suggest that the shapes and peaks in the abundance distribution functions provide a direct probe of the accretion histories of galaxies.

If the favored hierarchical cosmological model is correct, then the Milky Way system should have accreted ~100-200 luminous satellite galaxies in the past \~12 Gyr. We model this process using a hybrid semi-analytic plus N-body approach which distinguishes explicitly between the evolution of light and dark matter in accreted satellites. This distinction is essential to our ability to produce a realistic stellar halo, with mass and density profile much like that of our own Galaxy, and a surviving satellite population that matches the observed number counts and structural parameter distributions of the satellite galaxies of the Milky Way. Our model stellar halos have density profiles which typically drop off with radius faster than those of the dark matter. They are assembled from the inside out, with the majority of mass (~80%) coming from the \~15 most massive accretion events. The satellites that contribute to the stellar halo have median accretion times of ~9 Gyr in the past, while surviving satellite systems have median accretion times of ~5 Gyr in the past. This implies that stars associated with the inner halo should be quite different chemically from stars in surviving satellites and also from stars in the outer halo or those liberated in recent disruption events. We briefly discuss the expected spatial structure and phase space structure for halos formed in this manner. Searches for this type of structure offer a direct test of whether cosmology is indeed hierarchical on small scales.

[abridged] A "missing data" algorithm is developed to model (ie, reconstruct) the three-dimensional velocity distribution function of a sample of stars using data (velocity measurements) every one of which has one dimension unmeasured (the radial direction). It also accounts for covariant measurement uncertainties on the tangential velocity components. The algorithm is applied to tangential velocities measured in a kinematically unbiased sample of 11,865 stars taken from the Hipparcos catalog. The local stellar velocity distribution function of each of a set of 20 color-selected subsamples is modeled as a mixture of two three-dimensional Gaussian ellipsoids of arbitrary relative responsibility. In the fitting, one Gaussian (the "halo") is fixed at the known mean velocity and velocity variance tensor of the Galaxy halo, and the other (the "disk") is allowed to take arbitrary mean and arbitrary variance tensor. The mean and variance tensor (commonly the "velocity ellipsoid") of the disk velocity distribution are both found to be strong functions of stellar color, with long-lived populations showing larger velocity dispersion, slower mean rotation velocity, and smaller vertex deviation than short-lived populations. The local standard of rest (LSR) is inferred in the usual way and the Sun's motion relative to the LSR is found to be (U,V,W)_{\odot}=(10.1,4.0,6.7)+/-(0.5,0.8,0.2) km/s. Artificial data sets are made and analyzed, with the same error properties as the Hipparcos data, to demonstrate that the analysis is unbiased. The results are shown to be insensitive to the assumption that the velocity distributions are Gaussian.

We present the first results from our survey of stellar substructure in the outskirts of M31 using the Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope. We discuss the stellar populations associated with five prominent stellar overdensities discovered during the course of our panoramic ground-based imaging survey with the Isaac Newton Telescope Wide-Field Camera (INT WFC); a sixth pointing targets a region of `clean' halo. The colour-magnitude diagrams (CMDs), which contain between ~10,000--90,000 stars and reach four magnitudes below the horizontal branch, reveal clear variations in morphology between most fields, indicating that the age and/or metallicity mix of stars is not constant at large radius. This directly confirms the existence of large-scale population inhomogeneities within the halo of M31 and lends further support to the notion that M31 has formed, at least in part, through satellite accretions. We find a striking similarity between the populations of the giant stellar stream and those of another overdensity, the NE shelf, which lies north-east of the galaxy center. If these overdensities are associated with the same population, then the difference in their red clump magnitudes implies the NE shelf lies in front of the stream by several tens of kpc, in good agreement with recent orbit calculations for the stream progenitor.

We examine the nature of the progenitor of the giant stellar stream in M31 using as constraints new radial velocity measurements of stream red giant stars (presented in the companion paper by Guhathakurta et al. 2005, astro-ph/0406145) along with other M31 data sets available in the literature. We find that the observations are best fit by orbits that are highly eccentric and close to edge-on, with apo- to peri-center ratios of order 25 - 30, and with apocenters at or only slightly beyond the southern edge of the current data. Among these orbits, we are able to find a few that plausibly connect the stream with the northern spur or with the low-surface-brightness feature of similar high metallicity as the stream (originally reported by Ferguson et al. 2002) to the east of M31's center. In the latter case, if the connection is real, then the eastern debris should lie well in front of M31 near the apocenter of the orbit. Both the width of the debris and velocity dispersion measurements imply a rough lower limit on the mass of the progenitor of 10^8 M_sun. We use this limit and our orbits to discuss which of M31's satellites could be plausibly associated with the stream. In addition, we predict that the width of the stream should increase beyond the southern edge of the current data around the apocenter of the orbit and that the line-of-sight velocity dispersion should exhibit significant variations along the stream.

M giants recovered from the Two Micron All-Sky Survey (2MASS) have recently been used to map the position and velocity distributions of tidal debris from the Sagittarius (Sgr) dwarf spheroidal galaxy entirely around the Galaxy. We compare this data set to both test particle orbits and N-body simulations of satellite destruction run within a variety of rigid Milky Way potentials and find that the mass of the Milky Way within 50 kpc of its center should be 3.8-5.6 x 10^11 Msun in order for any Sgr orbit to simultaneously fit the velocity gradient in the Sgr trailing debris and the apocenter of the Sgr leading debris. Orbital pole precession of young debris and leading debris velocities in regions corresponding to older debris provide contradictory evidence in favor of oblate/prolate Galactic halo potentials respectively, leading us to conclude that the orbit of Sgr has evolved over the past few Gyr. Based upon the velocity dispersion and width along the trailing tidal stream we estimate the current bound mass of Sgr to be M_Sgr = 2 - 5 x 10^8 Msun independant of the form of the Galactic potential; this corresponds to a range of mass to light ratios (M/L)_Sgr = 14 - 36 (M/L)_Sun for the Sgr core. Models with masses in this range best fit the apocenter of leading Sgr tidal debris when they orbit with a radial period of roughly 0.85 Gyr and have periGalactica and apoGalactica of about 15 kpc and 60 kpc respectively. These distances will scale with the assumed distance to the Sgr dwarf and the assumed depth of the Galactic potential. The density distribution of debris along the orbit in these models is consistent with the M giant observations, and debris at all orbital phases where M giants are obviously present is younger (i.e. was lost more recently from the satellite) than the typical age of a Sgr M giant star.

M giants selected from the Two Micron All Sky Survey (2MASS) have been used to trace streams of tidal debris apparently associated with the Sagittarius dwarf spheroidal galaxy (Sgr) that entirely encircle the Galaxy. While the Sgr M giants are generally aligned with a single great circle on the sky, we measure a difference of 10.4 +- 2.6 degrees between the mean orbital poles of the great circles that best fit debris leading and trailing Sgr, which can be attributed to the precession of Sgr's orbit over the range of phases explored by the data set. Simulations of the destruction of Sgr in potentials containing bulge, disk and halo components best reproduce this level of precession along the same range of orbital phases if the potential contours of the halo are only slightly flattened, with the ratio between the axis length perpendicular to and in the disk in the range q = 0.90-0.95 (corresponding to isodensity contours with q_\rho ~ 0.83 - 0.92). Oblate halos are strongly preferred over prolate (q_\rho > 1) halos, and flattenings in the potential of q <= 0.85 (q_\rho <= 0.75) and q >= 1.05 (q_\rho >= 1.1) are ruled out at the 3-sigma level. More extreme values of q <= 0.80 (q_\rho <= 0.6) and q >= 1.25 (q_\rho >= 1.6) are ruled out at the 7-sigma and 5-sigma levels respectively. These constraints will improve as debris with larger separation in orbital phase can be found.

We present results from a large spectroscopic survey of M31 red giants using the Keck telescope/DEIMOS. Photometric pre-screening, based on the 100A-wide DDO51 band centered on the Mgb/MgH feature, was used to select spectroscopic targets. Red giant candidates were targeted in a small field on M31's giant southern tidal stream at a projected distance of 31kpc from the galaxy center. We isolate a clean sample of 68 giants by removing contaminants (foreground Galactic dwarfs and background galaxies) using spectroscopic, imaging, and photometric methods. About 65% of the M31 stars are found to be members of the stream, while the rest appear to be members of the general halo population. The mean (heliocentric) radial velocity of the stream in our field is -458 km/s, or -158 km/s relative to M31's systemic velocity, in good agreement with recent measurements at other stream locations. The intrinsic velocity dispersion of the stream is constrained to be 15_{-15}^{+8} km/s (90% confidence limits). The companion paper by Font et al. (2004, astro-ph/0406146) discusses possible orbits, implications of the coldness of the stream, and progenitor satellite properties. The kinematics (and perhaps [Fe/H] distribution) of our halo sample indicate that it is different from other M31 halo samples; this may be an indication of substructure in the halo. The stream seems to have a higher mean [Fe/H] than the halo, -0.51 vs -0.74 dex, and a smaller [Fe/H] spread. The stream's high metallicity implies that its progenitor must have been a luminous dwarf galaxy. The CaII triplet strengths of the M31 giants are generally consistent with photometric estimates of their metallicity (derived by fitting RGB fiducials in the color-magnitude diagram). There is indirect evidence of intermediate-age stars in the stream.

We present moderate resolution (~6 km/s) spectroscopy of 284 M giant candidates selected from the Two Micron All Sky Survey photometry. Radial velocities (RVs) are presented for stars mainly in the south, with a number having positions consistent with association to the trailing tidal tail of the Sagittarius (Sgr) dwarf galaxy. The latter show a clear RV trend with orbital longitude, as expected from models of the orbit and destruction of Sgr. A minimum 8 kpc width of the trailing stream about the Sgr orbital midplane is implied by verified RV members. The coldness of this stream (dispersion ~10 km/s) provides upper limits on the combined contributions of stream heating by a lumpy Galactic halo and the intrinsic dispersion of released stars, which is a function of the Sgr core mass. The Sgr trailing arm is consistent with a Galactic halo containing one dominant, LMC-like lump, however some lumpier halos are not ruled out. An upper limit to the total M/L of the Sgr core is 21 in solar units. A second structure that roughly mimics expectations for wrapped, leading Sgr arm debris crosses the trailing arm in the Southern Hemisphere; however, this may also be an unrelated tidal feature. Among the <13 kpc M giants toward the South Galactic Pole are some with large RVs that identify them as halo stars, perhaps part of the Sgr leading arm near the Sun. The positions and RVs of Southern Hemisphere M giants are compared with those of southern globular clusters potentially stripped from the Sgr system and support for association of Pal 2 and Pal 12 with Sgr debris is found. Our discussion includes description of a masked-filtered cross-correlation methodology that achieves better than 1/20 of a resolution element RVs in moderate resolution spectra.

If the favored hierarchical cosmological model is correct, then the Milky Way system should have accreted and subsequently tidally destroyed ~100 low-mass galaxies in the past ~12 Gyr. We model this process using a hybrid semi-analytic plus N-body approach and show that the disrupted systems lead naturally to stellar halos with masses and density profiles much like the stellar halo of our own Galaxy. We present predictions for the properties of stellar halos and show that ours is likely dominated by substructure beyond \~50kpc and more spatially smooth within that radius. The average stellar halo density profile is expected to drop off with radius more quickly than that of the dark matter because the stellar halo is formed from the most tightly bound material in infalling systems, while the majority of the accreted dark matter is stripped and deposited at larger radii. We argue that stars associated with the inner halo should be quite different chemically from stars in surviving satellites and also different from stars in the outer halo or those liberated in recent disruption events. We discuss how deep halo surveys and chemical probes may be useful tools for uncovering evidence of accretion. Searches of this kind offer a direct test of whether cosmology is indeed hierarchical on small scales.

2MASS has provided a three-dimensional map of the >360 degree, wrapped tidal tails of the Sagittarius (Sgr) dwarf spheroidal galaxy, as traced by M giant stars. With the inclusion of radial velocity data for stars along these tails, strong constraints exist for dynamical models of the Milky Way-Sgr interaction. N-body simulations of Sgr disruption with model parameters spanning a range of initial conditions (e.g., Sgr mass and orbit, Galactic rotation curve, halo flattening) are used to find parameterizations that match almost every extant observational constraint of the Sgr system. We discuss the implications of the Sgr data and models for the orbit, mass and M/L of the Sgr bound core as well as the strength, flattening, and lumpiness of the Milky Way potential.

N-body simulations are used to model the tidal disruption of the Sagittarius (Sgr) dwarf galaxy with constraints set by the positions and velocities of M giants in the Sgr tidal arms recently revealed by the Two Micron All-Sky Survey (2MASS). The simulated Sgr dwarf is placed on a variety of orbits within a Milky Way potential parameterized by variable circular velocities, halo flattenings and radial profiles. Two hundred separate test particle orbits have been used to explore a wide range of model Milky Way potentials and dwarf galaxy characteristics. The family of models is delimited by the data to a relatively narrow allowed range of parameters, and then input into N-body simulations. We present our best-fitting model, and discuss the orbital period, apoGalacticon distance, current space velocity, mass-to-light ratio, and other characteristics of the Sgr dwarf. In addition, we discuss the implications of this model for the flattening of the Galactic halo.

We present the results of multisite observations spanning two years on the pre--main-sequence (PMS) star V351 Ori. A total of around 180 hours of observations over 29 nights have been collected, allowing us to measure five different periodicities, most likely related to the delta Scuti variability of V351 Ori. Comparison with the predictions of linear nonadiabatic radial pulsation models put stringent constraints on the stellar parameters and indicate that the distance to V351 Ori is intermediate between the lower limit measured by Hipparcos (210 pc) and that of the Orion Nebula (450 pc). However, radial pulsation models are unable to reproduce all of the observed frequencies with a single choice of (M, L, and Te), suggesting the presence of additional nonradial modes.

We have used the Very Large Array, linked with the Pie Town Very Long Baseline Array antenna, to determine astrometric positions of 19 radio stars in the International Celestial Reference Frame (ICRF). The positions of these stars were directly linked to the positions of distant quasars through phase referencing observations. The positions of the ICRF quasars are known to 0.25 mas, thus providing an absolute reference at the angular resolution of our radio observations. Average values for the errors in our derived positions for all sources were 13 mas and 16 mas in R.A. and declination respectively, with accuracies approaching 1-2 mas for some of the stars observed. Differences between the ICRF positions of the 38 quasars, and those measured from our observations showed no systematic offsets, with mean values of -0.3 mas in R.A. and -1.0 mas in declination. Standard deviations of the quasar position differences of 17 mas and 11 mas in R.A. and declination respectively, are consistent with the mean position errors determined for the stars. Our measured positions were combined with previous Very Large Array measurements taken from 1978-1995 to determine the proper motions of 15 of the stars in our list. With mean errors of approximately 1.6 mas/yr, the accuracies of our proper motions approach those derived from Hipparcos, and for a few of the stars in our program, are better than the Hipparcos values. Comparing the positions of our radio stars with the Hipparcos catalog, we find that at the epoch of our observations, the two frames are aligned to within formal errors of approximately 3 mas. This result confirms that the Hipparcos frame is inertial at the expected level.

We describe first results of a spectroscopic probe of selected fields from the Grid Giant Star Survey. Multifiber spectroscopy of several hundred stars in a strip of eleven fields along delta approximately -17^{circ}, in the range 12 <~ alpha <~ 17 hours, reveals a group of 8 giants that have kinematical characteristics differing from the main field population, but that as a group maintain coherent, smoothly varying distances and radial velocities with position across the fields. Moreover, these stars have roughly the same abundance, according to their MgH+Mgb absorption line strengths. Photometric parallaxes place these stars in a semi-loop structure, arcing in a contiguous distribution between 5.7 and 7.9 kpc from the Galactic center. The spatial, kinematical, and abundance coherence of these stars suggests that they are part of a diffuse stream of tidal debris, and one roughly consistent with a wrapped, leading tidal arm of the Sagittarius dwarf spheroidal galaxy.

Recent observations of surface brightness distributions of both Milky Way and M31 satellite galaxies have revealed many instances of sudden changes or breaks in the slope of the surface brightness profiles (at some break radius, r_break). These breaks are often accompanied by increasingly elliptical isophotes and sometimes by isophote twisting. We investigate the hypothesis of a tidal origin for these features by applying the same ellipse fitting techniques that are used on observed galaxies to numerical simulations of the destruction of satellites, represented by spherical, single-component systems. We examine how observed quantities such as r_break, ellipticity and position angle of the fitted ellipses and amplitude of extra-break population vary with the satellite's orbital eccentricity and phase, and our viewpoint relative to the orbit. We also look at orbit and viewpoint dependence of the rate of change of the latter three quantities with radius. We find that there are trends with orbital phase and eccentricity in all observed quantities, many of which are preserved through a wide variety of viewing angles. In particular, a generic feature of all simulations is a depletion zone just interior to an excess zone, regions in which the surface brightness is lower and higher, respectively, than the initial profile. A clear interpretation of any individual image, however, is likely to be hampered by the dependence of the observable features on these multiple parameters. For example, breaks can be excited by several physical processes and can occur well within the bound satellite population. Nevertheless, we do find we can place loose constraints on the tidal radius, mass loss rate, orbital type and phase of the satellite, and nature of breaks using photometric data alone.

We investigate the interaction history of the M31 sub-group of galaxies by comparing surface photometry of its two nearest satellites, M32 and NGC 205, with N-body simulations of satellite destruction. The recent discovery of a giant stream in M31's outer halo, apparently pointed towards M32 and NGC 205, makes such an investigation particularly relevant. The observational component of this study is based on 1.7 x 5.0 degree B- and I- band CCD mosaic images centered on M31 and covering both satellites. The simulation component of this project is based on the analysis of single-component, spherical model satellites that are tidally disrupted through interactions with their parent galaxy. Generic features of the simulations include an excess in the surface brightness profile at large radii, a depletion zone at intermediate radii, and isophotal elongation and twists that are coincident with breaks in the brightness profile. The satellites M32 and NGC 205 display most of these features consistently across the B and I bands, strongly suggestive of tidal interaction and probable stripping by M31. We discuss what these observed features can tell us about the satellites' orbital parameters and histories. Specifically, M32 is found to be on a highly eccentric orbit and away from pericenter. Investigating M32's unusual combination of high surface brightness and low luminosity (the hallmark of compact ellipticals), we conclude on empirical grounds that M32 is not likely to be the residual core of a tidally-stripped normal elliptical galaxy as has been suggested, but rather that its precursor was intrinsically compact.

Surface photometry of the M31 satellites M32 and NGC 205 is compared to numerical simulations of satellite destruction to constrain orbital parameters and the interaction history of the M31 subgroup. Our analysis reveals the following preliminary results: (1) Generic features of tidal disruption in the simulations include an extended ``extra-tidal'' excess region and an inner depletion zone, both of which are observed in M32 and NGC 205; (2) M32 is likely to be on a highly eccentric orbit well away from pericenter; (3) Surface brightness and luminosity evolution estimates for M32, the prototypical compact elliptical galaxy, imply that it is not simply the residual core of a tidally-stripped normal elliptical galaxy, but was instead formed in a truncated state.

CDM simulations predict that there are hundreds of lumps of with masses greater than 10e7 solar masses in the Milky Way halo. However, we know of only a dozen dwarf satellites close to this mass. Are these lumps simply lacking in stars or is there a fundamental flaw in our most popular cosmology? By studying the tidal debris of known satellites we can potentially address this question. In this paper, we quantify the the effects of the dark matter lumps on tidal tails. The lumps scatter stars in the tidal tails from their original orbits producing a distinctive signature. We simulate debris evolution in smooth and lumpy halos potentials and use our simulations to motivate and test a statistical measure of the degree of scattering apparent in the angular position and radial velocity measurement of debris stars - the ``scattering index''. We find that the scattering index can in general distinguish between the levels of substructure predicted by CDM cosmologies and smooth Milky Way models, but that the sensitivity of the debris depends on the orientation of the parent satellite's orbit relative to the largest lumps orbits. We apply our results to the carbon star stream associated with the Sagittarius dwarf galaxy (Sgr) and find that these stars appear to be more scattered than we expect for debris orbiting in a smooth halo. However, the degree of scattering is entirely consistent with that expected due to the influence of the Large Magellanic Cloud, which is on an orbit that intersects Sgr's own. We conclude that the current data is unable to constrain CDM models. Nevertheless, our study suggests that future data sets of debris stars associated with other Milky Way satellites could provide strong constraints on CDM models.

In numerous studies of the outer Galactic halo some evidence for accretion has been found. If the outer halo did form in part or wholly through merger events, we might expect to find coherent streams of stars and globular clusters following similar orbits as their parent objects, which are assumed to be present or former Milky Way dwarf satellite galaxies. We present a study of this phenomenon by assessing the likelihood of potential descendant ``dynamical families'' in the outer halo. We conduct two analyses: one that involves a statistical analysis of the spatial distribution of all known Galactic dwarf satellite galaxies (DSGs) and globular clusters, and a second, more specific analysis of those globular clusters and DSGs for which full phase space dynamical data exist. In both cases our methodology is appropriate only to members of descendant dynamical families that retain nearly aligned orbital poles today. Since the Sagittarius dwarf (Sgr) is considered a paradigm for the type of merger/tidal interaction event for which we are searching, we also undertake a case study of the Sgr system and identify several globular clusters that may be members of its extended dynamical family. (ABRIDGED)

We present direct measurements of the limb-darkened intensity profiles of the late-type giant stars HR5299, HR7635, and HR8621 obtained with the Navy Prototype Optical Interferometer (NPOI) at the Lowell Observatory. A triangle of baselines with lengths of 18.9 m, 22.2 m, and 37.5 m was used. We utilized squared visibility amplitudes beyond the first minimum, as well as triple amplitudes and phases in up to 10 spectral channels covering a wavelength range of ~650 nm to ~850 nm. We find that our data can best be described by featureless symmetric limb-darkened disk models while uniform disk and fully darkened disk models can be rejected. We derive high-precision angular limb-darkened diameters for the three stars of 7.44 mas +/- 0.11 mas, 6.18 mas +/- 0.07 mas, and 6.94 mas +/- 0.12 mas, respectively. Using the HIPPARCOS parallaxes, we determine linear limb-darkened radii of 114 R$_\odot \pm $13 R$_\odot$, 56 R$_\odot \pm $4 R$_\odot$, and 98 R$_\odot \pm $9 R$_\odot$, respectively. We compare our data to a grid of Kurucz stellar model atmospheres, with them derive the effective temperatures and surface gravities without additional information, and find agreement with independent estimates derived from empirical calibrations and bolometric fluxes. This confirms the consistency of model predictions and direct observations of the limb-darkening effect.

We present light curves for 149 sources monitored with the Green Bank Interferometer. The light curves are at two radio frequencies (approximately 2.5 and 8.2 GHz) and range from 3 to 15 yrs in length, covering the interval 1979--1996, and have a typical sampling of one flux density measurement every 2 days. We have used these light curves to conduct various variability analysis (rms flux density variations and autoregressive, integrated, moving average modeling) of these sources. We find suggestive, though not unambiguous evidence, that these sources have a common, broadband mechanism for intrinsic variations, in agreement with previous studies of a subset of these source. We also find that the sources generally display a short-term variability (~ 10 d) that arises from radio-wave scattering in an extended medium. These conclusions extend those of Fiedler et al. (1987) who used a sub-sample of these data. The primary motivation for this monitoring program was the identification of extreme scattering events. In an effort to identify ESEs in a systematic manner, we have taken the wavelet transform of the light curves. We find 15 events in the light curves of 12 sources that we classify as probable ESEs. However, we also find that five ESEs previously identified from these data do not survive our wavelet selection criteria. Future identification of ESEs will probably continue to rely on both visual and systematic methods. Instructions for obtaining the data are also presented.

We examine the detectability and interpretation of debris trails caused by satellite disruption in external galaxies using semi-analytic approximations for the dependence of streamer length, width and surface brightness on satellite and primary galaxy characteristics. The semi-analytic method is tested successfully against N-body simulations and then applied to three representative astronomical applications. First, we show how streamer properties can be used to estimate mass-to-light ratios $\Upsilon$ and streamer ages of totally disrupted satellites, and apply the method to the stellar arc in NGC 5907. Second, we discuss how the lack of observed tidal debris around a satellite can provide an upper limit on its mass-loss rate, and, as an example, derive the implied limits on mass-loss rates for M32 and NGC 205 around Andromeda. Finally, we point out that a statistical analysis of streamer properties might be applied to test and refine cosmological models of hierarchical galaxy formation, and use the predicted debris from a standard $\Lambda$CDM realization to test the feasibility of such a study. Using the Local Group satellites and the few known examples of debris trails in the Galaxy and in external systems, we estimate that the best current techniques could characterize the brightest ($R < 29$ mag/ arcsec$^{2}$) portions of the youngest (3 dynamical periods) debris streamers. If systematics can be controlled, planned large-aperture telescopes such as CELT and OWL may allow fainter trails to be detected routinely and thus used for statistical studies such as those required for tests of galaxy formation.

The launching of up to four astrometric missions in the next decade will enable us to make measurements of stars in tidal streamers from Galactic satellites of sufficient accuracy to place strong constraints on the mass distribution in the Milky Way. In this paper we simulate observations of debris populations in order to assess the required properties of any data set chosen to implement this experiment. We apply our results to find the desired target properties of stars (e.g., accuracy of velocity and proper motion measurements) associated with the dwarf spheroidal satellites of the Milky Way.

As part of our survey for substructure in the Milky Way halo as traced by giant stars, and to look for tidal stellar debris in the halo viable for measurement of the Galactic mass potential with SIM, we explore the distribution of stars beyond the nominal tidal radius of, but still associated with, the Carina dwarf spheroidal galaxy. We make use of a photometric technique described in Majewski et al. (1999) to identify giant star candidates at the distance and metallicity of Carina across the entire extent of the photometric survey covering some 2.2 deg^2 on and around Carina. These Carina-associated giant candidates are identified by a combination of (1) their M-DDO51 colors, and (2) by locations in the color-magnitude diagram commensurate with the Carina RGB in the core of the galaxy. The density distribution of the extratidal giant candidates is similar to that found from previous statistical star-counting approaches. However, we are now able to pinpoint actual, remotely-situated Carina stars individually. Because we can exclude foreground dwarf stars, our approach allows greater sensitivity and the ability to map the detailed two-dimensional distribution of extended Carina populations to much larger radii, with smaller aperture telescopes, than other techniques. Moreover, we identify candidate lists of widely displaced Carina-associated stars bright enough for spectroscopic studies of large-scale dynamical and metallicity properties of the system, and for astrometric study by SIM. We obtained spectroscopy for three such ``extratidal'' stars and from their radial velocities conclude that all three are associated with Carina.

ESA's Global Astrometric Interferometer for Astrophysics (GAIA) holds the promise of mapping out the detailed phase space structure of the Galactic halo by providing unprecedented annual proper motion and parallax of $1-10\muas$ astrometric accuracy (Gilmore et al. 1998). Unlike NASA's Space Interferometry Mission (SIM), which will achieve similar accuracies but is a pointed instrument, GAIA will be able to construct a global catalogue of the halo. Here we study proper motions of giant branch stars in a tidal debris torn from a small satellite system in the halo. We follow the evolution of a cold stream on a polar orbit between 8-50 kpc in a variety of histories of the Galactic potential, and observe the bright ($V<18$mag) members of the debris tail with GAIA accuracy. We simulate effects due to the growing or flipping of the Galactic disk over the past 4 Gyrs or the perturbation from a massive accreted lump such as the progenitor of the Magellanic Clouds. Our simulations suggest that the results of Johnston, Zhao, Spergel & Hernquist (1999) and Helmi, Zhao & de Zeeuw (1999) for static Galactic potentials are likely to be largely generalizable to realistic time-dependent potentials: a tidal debris remains cold in spite of evolution and non-axial symmetry of the potential. GAIA proper motion measurements of debris stars might be used to probe both Galactic structure and Galactic history. We also study several other factors influencing our ability to identify streams, including accuracy of radial velocity and parallax data from GAIA, and contamination from random field stars. We conclude that nearby, cold streams could be detected with GAIA if these cousins of the Sagittarius stream exist.

FAME is designed to perform an all-sky, astrometric survey with unprecedented accuracy. It will create a rigid astrometric catalog of 4x10^7 stars with 5 < m_V < 15. For bright stars, 5 < m_V < 9, FAME will determine positions and parallaxes accurate to < 50 microarcseconds, with proper motion errors < 50 microarcseconds/year. For fainter stars, 9 < m_V < 15, FAME will determine positions and parallaxes accurate to < 500 microarcseconds, with proper motion errors < 500 microarcseconds/year. It will also collect photometric data on these 4 x 10^7 stars in four Sloan DSS colors.

As part of the Selected Areas Starcount Survey (SASS), a CCD survey to V>21, we have obtained VI photometry of two fields at b=+-40 aligned roughly with an extrapolation of the major axis of the Sagittarius (Sgr) dwarf spheroidal galaxy. Comparison of the color-magnitude diagram (CMD) for some of these fields to the CMDs of fields reflected about the Galactic l=0 meridian reveals an excess of stars at V_o=17.85 and 0.9 < (V-I)_o < 1.1 in the (l,b) = (11,-40) field. The excess stars have colors consistent with the Sgr red clump, and deeper CMD imaging in these locations shows evidence of a main-sequence turnoff at V=21 with the main-sequence extending to the limit of our data (V=24). The surface brightnesses we derive from either the potential excess of red clump stars or the apparent excess of MSTO stars are consistent with each other and with the results of other surveys at this latitude. No similar excess appears in our northern Galactic hemisphere fields near l=353, b=+41 field. We have obtained spectroscopy of all 30 candidate red clump stars in the range 0.9 < (V-I)_o < 1.1 and 17.75 < V_o < 17.95. The radial velocity distribution of the stars, while dissimilar from expectations of Galactic structure models, does not show a contribution by stars near the Galactocentric radial velocity seen in other studies near the Sgr core. It is difficult to reconcile a photometric result that is consistent with other explorations of the Sagittarius stream with a radial velocity distribution that is apparently inconsistent. In a companion paper (Johnston et al. 1999), we discuss how some of the discrepancies are resolved if ur potential Sgr detection corresponds to a different Sgr tidal streamer than that detected by most other surveys.

We present a comparison of semi-analytic models of the phase-space structure of tidal debris with observations of stars associated with the Sagittarius dwarf galaxy (Sgr). We find that many features in the data can be explained by these models. The properties of stars 10-15 degrees away from the center of Sgr --- in particular, the orientation of material perpendicular to Sgr's orbit (c.f. Alard 1996) and the kink in the velocity gradient (Ibata et al 1997) --- are consistent with those expected for unbound material stripped during the most recent pericentric passage ~50 Myrs ago. The break in the slope of the surface density seen by Mateo, Olszewski & Morrison (1998) at ~ b=-35 can be understood as marking the end of this material. However, the detections beyond this point are unlikely to represent debris in a trailing streamer, torn from Sgr during the immediately preceding passage ~0.7 Gyrs ago, but are more plausibly explained by a leading streamer of material that was lost more that 1 Gyr ago and has wrapped all the way around the Galaxy. The observations reported in Majewski et al (1999) also support this hypothesis. We determine debris models with these properties on orbits that are consistent with the currently known positions and velocities of Sgr in Galactic potentials with halo components that have circular velocities v_circ=140-200 km/s. The best match to the data is obtained in models where Sgr currently has a mass of ~10^9 M_sun and has orbited the Galaxy for at least the last 1 Gyr, during which time it has reduced its mass by a factor of 2-3, or luminosity by an amount equivalent to ~10% of the total luminosity of the Galactic halo. These numbers suggest that Sgr is rapidly disrupting and unlikely to survive beyond a few more pericentric passages.

Future astrometric satellites, such as SIM (NASA's Space Interferometric Mission) and GAIA (ESA's Global Astrometric Interferometer for Astrophysics), hold the promise of mapping out the detailed phase space structure of the Galactic halo by providing unprecedented annual proper motion and parallax of $1-10\mu$as astrometric accuracy. Here we show that proper motions of hundred or so giant branch stars in a tidal debris torn from a small satellite (a $10^{5-7}L_\odot$ Galactic dwarf galaxy or globular cluster) in the halo is sensitive to the current Galactic potential and its past evolution. We follow the evolution of a cold (velocity dispersion of 10 km/s) stream on a nearby (between 8-50 kpc) polar orbit in a variety of histories of the potential of the Galaxy, and observe the bright ($V<18$mag) members of the debris tail with GAIA accuracy. We simulate effects due to the growing or flipping of the Galactic disk over the past 4 Gyrs or the perturbation from a massive accreted lump such as the progenitor of the Magellanic Clouds. We study various factors influencing our ability to identify streams, including contamination from field stars, accuracy of radial velocity and distance data and evolution and non-axial symmetry of the potential. Our simulations suggest that nearby, cold streams could be detected with GAIA if these cousins of the Sagittarius stream exist. Results of Johnston, Zhao, Spergel & Hernquist (1999) and Helmi, Zhao & de Zeeuw (1999) for static Galactic potentials are likely to be largely generalizable to moderately time-evolving potentials. SIM and GAIA measurements of debris stars might be used to probe both Galactic structure and Galactic history.

We explore the use of tidal streams from Galactic satellites to recover the potential of the Milky Way. Our study is motivated both by the discovery of the first lengthy stellar stream in the halo (\cite{it98}) and by the prospect of measuring proper motions of stars brighter than 20th magnitude in such a stream with an accuracy of $\sim 4\mu as/$yr, as will be possible with the Space Interferometry Mission (SIM). We assume that the heliocentric radial velocities of these stars can be determined from supporting ground-based spectroscopic surveys, and that the mass and phase-space coordinates of the Galactic satellite with which they are associated will also be known to SIM accuracy. Using results from numerical simulations as trial data sets, we find that, if we assume the correct form for the Galactic potential, we can predict the distances to the stars as a consequence of the narrow distribution of energy expected along the streams. We develop an algorithm to evaluate the accuracy of any adopted potential by requiring that the satellite and stars recombine within a Galactic lifetime when their current phase-space coordinates are integrated backwards. When applied to a four-dimensional grid of triaxial logarithmic potentials, with varying circular velocities, axis ratios and orientation of the major-axis in the disk plane, the algorithm can recover the parameters used for the Milky Way in a simulated data set to within a few percent using only 100 stars in a tidal stream.

If the Galaxy formed hierarchically through the accretion of smaller satellite galaxies we might hope to find signatures of this in the halo's phase-space distribution. I review theoretical ideas about what form these signatures should take, concentrating in particular on the characteristics of debris from the tidal disruption of satellites less than the size of the Large Magellanic Clouds and at distances greater than 10 kpc. Under these circumstances, tidal debris can remain aligned in streams of stars close to the satellite's original orbit for the lifetime of the Galaxy and we can understand the evolution of its structure analytically. This understanding has many applications to observations and I will discuss a few examples: reconstructing the properties of a long-dead satellite from its debris; measuring the mass loss rate from Galactic satellites through the density of extra-tidal stars; and using observations of motions of tidal stream stars to measure the mass and size of the Milky Way.

The luminous high-redshift (z=2.17) quasar S5 0836+710 has been observed at 5GHz with the VSOP. We compare the properties of three images obtained from the observation: a low-resolution ground array image (dynamic range 4600:1), a full-resolution VSOP image (900:1), and an image made with only the space baselines (200:1). The space baselines alone are sufficient for a reliable recovery of the source structure, within the limits of the achieved spatial sampling of the visibility data. The curved jet ridge line observed in the images can be described by Kelvin-Helmholtz instabilities developing in a relativistic outflow with the Mach number of about 6. This description holds on the scales of up to 700h^-1 pc, and is shown to be consistent with variable apparent speeds observed in the jet.

Using CGRO/BATSE hard X-ray (HXR) data and GHz radio monitoring data from the Green Bank Interferometer (GBI), we have performed a long term study ($\sim$ 1800 days) of the unusual X-ray binary Cyg X-3 resulting in the discovery of a remarkable relationship between these two wavelength bands. We find that, during quiescent radio states, the radio flux is strongly anticorrelated with the intensity of the HXR emission. The relationship switches to a correlation with the onset of major radio flaring activity. During major radio flaring activity the HXR drops to a very low intensity during quenching in the radio and recovers during the radio flare. Injection of plasma into the radio jets of Cyg X-3 occurs during changes in the HXR emission and suggests that disk-related and jet-related components are responsible for the high energy emission.

We present the results of a study that uses numerical simulations to interpret observations of tidally disturbed satellites around the Milky Way. When analysing the simulations from the viewpoint of an observer, we find a break in the slope of the star count and velocity dispersion profiles in our models at the location where unbound stars dominate. We conclude that `extra-tidal' stars and enhanced velocity dispersions observed in the outskirts of Galactic satellites are due to contamination by stellar debris from the tidal interaction with the Milky Way. However, a significant bound population can exist beyond the break radius and we argue that it should not be identified with the tidal radius of the satellite. We also develop and test a method for determining the mass loss rate from a Galactic satellite using its extra-tidal population. We apply this method to observations of globular clusters and dwarf spheroidal satellites of the Milky Way, and conclude that a significant fraction of both satellite systems are likely be destroyed within the next Hubble time. Finally, we demonstrate that this mass loss estimate allows us to place some limits on the initial mass function (IMF) of stars in a cluster from the radial dependence of its present day mass function (PDMF).

We develop a semi-analytic method for determining the phase-space population of tidal debris along the orbit of a disrupting satellite galaxy and illustrate its use with a number of applications. We use this method to analyze Zhao's proposal that the microlensing events towards the Large Magellanic Cloud (LMC) might be explained by an appropriately placed tidal streamer, and find that his scenarios lead either to unacceptably high overdensities (10 -- 100%) in faint star counts (apparent magnitudes 17.5 -- 20.5) away from the Galactic plane or short timescales for the debris to disperse (10^8 years). We predict that the tidal streamers from the LMC and the Sagittarius dwarf galaxy currently extend over more than $2\pi$ in azimuth along their orbits. Assuming that each satellite has lost half of its primordial mass, we find that the streamers will have overdensities in faint star counts of 10 -- 100% and < 1% respectively, and conclude that this mass loss rate is unlikely for the LMC, but possible for Sagittarius. If the Galaxy has accreted one hundred $10^5-10^6 M_{\odot}$ objects (comparable to its current population of globular clusters) at distances of 20 -- 100 kpc during its lifetime then 10% of the sky will now be covered by tidal streamers.

The role that minor mergers have played in the formation and structure of the Milky Way is still an open question, about which there is much debate. We use numerical simulations to explore the evolution of debris from a tidally disrupted satellite, with the aim of developing a method that can be used to identify and quantify signatures of accretion in a survey of halo stars. For a Milky Way with a spherical halo, we find that debris from minor mergers can remain aligned along great circles throughout the lifetime of the Galaxy. We exploit this result to develop the method of Great Circle Cell Counts (GC3), which we test by applying it to artificially constructed halo distributions. Our results suggest that if as few as 1\% of the stars in a halo survey are accreted from the disruption of a single subsystem smaller than the Large Magellanic Cloud, GC3 can recover the great circle associated with this debris. The dispersion in GC3 can also be used to detect the presence of structure characteristic of accretion in distributions containing a much smaller percentage of material accreted from any single satellite.

We report multifrequency observations of the $\gamma$-ray blazar 0528+134 with the Effelsberg 100-m telescope, the IRAM 30-m telescope, and the NRL Green Bank Interferometer. 0528+134 underwent a major radio and mm outburst in 1993 a few month after a very strong outburst in high energy $\gamma$-rays. A similar behaviour is indicated for the weaker $\gamma$-ray outburst in May 1991 and the mm activity in the beginning of 1992. This result fits nicely to the general finding that blazars are bright in $\gamma$-rays preferrentially at the beginning of a radio outburst. 0528+134 was observed with a global VLBI array at 22$\,$GHz in November 1992. Within the IRIS-S and EUROPE geodetic VLBI observing campaigns the source was also observed regularly at 8.4$\,$GHz. A new VLBI component which is not present in data at 22 $\,$GHz from 1991 also exhibits superluminal motion. Backextra\-polating in time we find that the new component was expelled from the core in the first half of 1991, near the first $\gamma$-ray outburst in Mai 1991. Newest VLBI data indicate the appearance of another new component in 1994. Although this new component has yet been found only at two epochs, it appears to move superluminally and may have been expelled from the core between fall 1992 and summer 1993, i.e. near the time of the second gamma ray outburst and at the beginning of a strong radio outburst.

In this pilot investigation, precise optical positions in the FK5 system are presented for a set of 16 compact extragalactic radio sources, which will be part of the future radio--optical reference frame. The 0.9 m KPNO and CTIO telescopes equipped with 2K CCD's have been used for this project. The astrometric properties of these instruments are investigated in detail. New techniques of using wide field CCD observations for astrometry in general are developed. An internal precision of 5 to 31 mas in position per single exposure is found, depending on the brightness of the object. The tie to the primary optical reference system is established by photographic astrometry using dedicated astrographs on both hemispheres. An accuracy of $\approx 30$ mas per source is estimated for the multi--step reduction procedure when based on the future Hipparcos catalog, while the FK5--based positions suffer from system errors of 100 to 200 mas as compared to the radio positions. This work provides a contribution to the international effort to link the Hipparcos instrumental coordinate system to the quasi--inertial VLBI radio reference frame. Precise radio and optical astrometry of a large sample of compact extragalactic sources will also contribute to the astrophysics of these objects by comparing the respective centers of emission at the optical and radio wavelengths.

We report on multifrequency radio observations made with the Effelsberg 100-m telescope, the IRAM 30-m telescope and the Green Bank Interferometer between 1992 and 1994 of the $\gamma$-ray quasar 0528+134. We present a new VLBI based map of 0518+134 at 22 GHz with sub-mas angular resolution observed in November 1992, which shows a one-sided core jet structure of $\sim$ 5 mas length. One component, seen also previously at 8.4 GHz & 22 GHz, shows apparent superluminal motion with $\beta_{app}$=4.4$\pm$1.7. In summer 1993 we observed a major outburst at all frequencies higher than a few GHz, which peaked a few months after a strong outburst in high energy $\gamma$-rays and showed the canonical time evolution of a cooling and/or expanding electron distribution. Our data indicate that the outburst in the EGRET range originated very close to the central object of the AGN and that a remnant of this outburst moved further outward in the jet until it became optically thin at radio frequencies after a few months. During the flare in July 1993 we observed with the Effelsberg 100-m telescope an unusually strong decrease of the flux density by about 50% at 4.75 GHz and 10.55 GHz and slightly less at 2.695 GHz. This behaviour is also seen in the monitoring data at 2.25 GHz and 8.3 GHz taken with the Green Bank Interferometer (NRL-GBI). An extreme scattering event by a small dense plasma cloud in the line of sight is able to match the observed time lag in the lightcurves if we take into account the mas-structure of the source and different spectra of the components on the basis of their brightness in the VLBI maps.

Numerical simulations of dwarf spheroidal galaxies undergoing several close encounters with the Milky Way are described. By comparing our models to observed properties of the recently discovered dwarf galaxy in Sagittarius (Sgr), we discuss implications of our results for the formation and evolution of the Milky Way system. We find that existing observations are not sufficient to allow us to place precise limits on either the orbit or the initial state of the dwarf. Debris from the ongoing tidal stripping of the Sagittarius galaxy are expected to form moving groups in the halo of the Galaxy and the discovery of such stars would strongly constrain the history and dynamical state of the dwarf. Furthermore, if Sgr is presently being disrupted, we predict that its remains will be detectable as a moving group in the halo for more than 1 Gyr. Thus, if similar accretion events have occurred in the recent history of the Galaxy, their aftereffects may still be observable.

The radio observations are based on dual frequency Mark--III VLBI data from 1979 until the end of 1993. Currently 436 sources define the celestial inertial frame. A comparison to the JPL94R01 catalog is made. Position vs. time plots of selected sources show significant apparent non linear proper motions. The status of the project and recent results from optical observations are presented.

The Sachs-Wolfe effect is known to produce large angular scale fluctuations in the Cosmic Microwave Background Radiation (CMBR) due to gravitational potential fluctuations. We show how the angular correlation function of the CMBR can be expressed explicitly in terms of the mass autocorrelation function $\xi (r)$ in the Universe. We derive analytic expressions for the angular correlation function and its multipole moments in terms of integrals over $\xi (r)$ or its second moment, $J_3 (r)$, which does not need to satisfy the sort of integral constraint that $\xi (r)$ must. We derive similar expressions for bulk flow velocity in terms of $\xi$ and $J_3$. One interesting result that emerges directly from this analysis is that, for angles $\theta$, there is a substantial contribution to the correlation function from a wide range of distance $r$ and that the radial shape of this contribution does not vary greatly with angle.