We suggest that the Virgo Overdensity (VOD) of stars in the stellar halo is the result of a radial dwarf galaxy merger that we call the Virgo Radial Merger. Because the dwarf galaxy passed very near to the Galactic center, the debris has a large range of energies but nearly zero $L_z$ angular momentum. The debris appears to extend from 5 to 50 kpc from the Sun in the Virgo region. We connect different parts of this merger debris to the Perpendicular and Parallel Streams (the Virgo Stellar Stream is associated with either or both of these streams), the Hercules-Aquila Cloud (HAC), and possibly the Eridanus Phoenix Overdensity (EriPhe). This radial merger can explain the majority of the observed moving groups of RR Lyrae and blue horizontal branch stars that have previously been identified in Virgo. This merger also produces debris in the Solar neighborhood similar to that identified as the Gaia-Enceladus or it Gaia-sausage merger. Orbits are provided for components of the Virgo Radial Merger progenitor and for debris that appears to be related to the Cocytos Stream, which was also recovered in the Virgo region.

We use LAMOST DR4 M giants combined with Gaia DR2 proper motions and ALLWISE photometry to obtain an extremely pure sample of Sagittarius (Sgr) stream stars. Using TiO5 and CaH spectral indices as an indicator, we selected out a large sample of M giant stars from M dwarf stars in LAMOST DR4 spectra. Considering the position, distance, proper motion and the angular momentum distribution, we obtained 164 pure Sgr stream stars. We find that the trailing arm has higher energy than the leading arm in same angular momentum. The trailing arm we detected extends to a heliocentric distance of $\sim 130$ kpc at $\tilde\Lambda_{\odot}\sim 170^{\circ}$, which is consistent with the feature found in RR Lyrae in \citet{sgr2017}. Both of these detections of Sgr, in M giants and in RR Lyrae, imply that the Sgr stream may contain multiple stellar populations with a broad metallicity range.

Stellar halo substructures were identified using statistical photometric parallax of blue main sequence turnoff stars from fourteen Sloan Digital Sky Survey stripes in the north Galactic cap. Four structures are consistent with previous measurements of the Sagittarius dwarf tidal stream: the leading tail, the "bifurcated" stream, the trailing tail, and Stream C. The stellar overdensity in Virgo, about 15 kpc from the Sun, could arise from the crossing of the Parallel Stream and a new, candidate stream dubbed the Perpendicular Stream. The data suggests the presence of a wide stream near NGC 5466, with a distance of 5 to 15 kpc. Measurements of the flattening of the smooth stellar halo from the fourteen stripes average q=0.58, with a dispersion of 0.04.

We propose and test a method for applying statistical photometric parallax to main sequence turn off (MSTO) stars in the Sloan Digital Sky Survey (SDSS). Using simulated data, we show that if our density model is similar to the actual density distribution of our data, we can reliably determine the density model parameters of three major substructures in the Milky Way halo using the computational resources available on MilkyWay@home (a twenty parameter fit). We fit the stellar density in SDSS stripe 19 with a smooth stellar spheroid component and three major streams. One of these streams is consistent with the Sagittarius tidal stream at $21.1$ kpc away, one is consistent with the trailing tail of the Sagittarius tidal stream in the north Galactic cap at $48$ kpc away, and one is possibly part of the Virgo Overdensity at $6$ kpc away. We find the one sigma widths of these three streams to be $1.0$ kpc, $17.6$ kpc, and $6.1$ kpc, respectively. The width of the trailing tail is extremely wide ($41$ kpc full width at half maximum). This large width could have implications for the shape of the Milky Way dark matter halo. The width of the Virgo Overdensity-like structure is consistent with what we might expect for a "cloud"-like structure; analysis of additional stripes of data are needed to outline the full extent of this structure and confirm its association with the Virgo Overdensity.

NIRSPEC is a 1-5 micron echelle spectrograph in use on the Keck II Telescope since 1999. The spectrograph is capable of both moderate (R~2,000) and high (R~25,000) resolution observations and has been a workhorse instrument across many astronomical fields, from planetary science to extragalactic observations. In the latter half of 2018, we will upgrade NIRSPEC to improve the sensitivity and stability of the instrument and increase its lifetime. The major components of the upgrade include replacing the spectrometer and slit-viewing camera detectors with Teledyne H2RG arrays and replacing all transputer-based electronics. We present detailed design, testing, and analysis of the upgraded instrument, including the finalized optomechanical design of the new 1-5 micron slit-viewing camera, final alignment and assembly of the science array, electronics systems, and updated software design.

We report on the design and performance of the Keck Cosmic Web Imager (KCWI), a general purpose optical integral field spectrograph that has been installed at the Nasmyth port of the 10 m Keck II telescope on Mauna Kea, HI. The novel design provides blue-optimized seeing-limited imaging from 350-560 nm with configurable spectral resolution from 1000 - 20000 in a field of view up to 20"x33". Selectable volume phase holographic (VPH) gratings and high performance dielectric, multilayer silver and enhanced aluminum coatings provide end-to-end peak efficiency in excess of 45% while accommodating the future addition of a red channel that will extend wavelength coverage to 1 micron. KCWI takes full advantage of the excellent seeing and dark sky above Mauna Kea with an available nod-and-shuffle observing mode. The instrument is optimized for observations of faint, diffuse objects such as the intergalactic medium or cosmic web. In this paper, a detailed description of the instrument design is provided with measured performance results from the laboratory test program and ten nights of on-sky commissioning during the spring of 2017. The KCWI team is lead by Caltech and JPL (project management, design and implementation) in partnership with the University of California at Santa Cruz (camera optical and mechanical design) and the W. M. Keck Observatory (observatory interfaces).

IRIS is a near-infrared (0.84 to 2.4 microns) integral field spectrograph and wide-field imager being developed for first light with the Thirty Meter (TMT). It mounts to the advanced optics (AO) system NFIRAOS and has integrated on-instrument wavefront sensors (OIWFS) to achieve diffraction-limited spatial resolution at wavelengths longer than 1 micron. With moderate spectral resolution (R ~4,000 - 8,000) and large bandpass over a continuous field of view, IRIS will open new opportunities in virtually every area of astrophysical science. It will be able to resolve surface features tens of kilometers across Titan, while also mapping the distant galaxies at the scale of an individual star forming region. This paper summarizes the entire design and capabilities, and includes the results from the nearly completed preliminary design phase.

The InfraRed Imaging Spectrograph (IRIS) will be a first-light client instrument for the Narrow Field Infrared Adaptive Optics System (NFIRAOS) on the Thirty Meter Telescope. IRIS includes three configurable tip/tilt (TT) or tip/tilt/focus (TTF) On-Instrument Wavefront Sensors (OIWFS). These sensors are positioned over natural guide star (NGS) asterisms using movable polar-coordinate pick-off arms (POA) that patrol an approximately 2-arcminute circular field-of-view (FOV). The POAs are capable of colliding with one another, so an algorithm for coordinated motion that avoids contact is required. We have adopted an approach in which arm motion is evaluated using the gradient descent of a scalar potential field that includes an attractive component towards the goal configuration (locations of target stars), and repulsive components to avoid obstacles (proximity to adjacent arms). The resulting vector field is further modified by adding a component transverse to the repulsive gradient to avoid problematic local minima in the potential. We present path planning simulations using this computationally inexpensive technique, which exhibit smooth and efficient trajectories.

IRIS (InfraRed Imaging Spectrograph) is the diffraction-limited first light instrument for the Thirty Meter Telescope (TMT) that consists of a near-infrared (0.84 to 2.4 $\mu$m) imager and integral field spectrograph (IFS). The IFS makes use of a lenslet array and slicer for spatial sampling, which will be able to operate in 100's of different modes, including a combination of four plate scales from 4 milliarcseconds (mas) to 50 mas with a large range of filters and gratings. The imager will have a field of view of 34$\times$34 arcsec$^{2}$ with a plate scale of 4 mas with many selectable filters. We present the preliminary design of the data reduction system (DRS) for IRIS that need to address all of these observing modes. Reduction of IRIS data will have unique challenges since it will provide real-time reduction and analysis of the imaging and spectroscopic data during observational sequences, as well as advanced post-processing algorithms. The DRS will support three basic modes of operation of IRIS; reducing data from the imager, the lenslet IFS, and slicer IFS. The DRS will be written in Python, making use of open-source astronomical packages available. In addition to real-time data reduction, the DRS will utilize real-time visualization tools, providing astronomers with up-to-date evaluation of the target acquisition and data quality. The quicklook suite will include visualization tools for 1D, 2D, and 3D raw and reduced images. We discuss the overall requirements of the DRS and visualization tools, as well as necessary calibration data to achieve optimal data quality in order to exploit science cases across all cosmic distance scales.

The Thirty Meter Telescope (TMT) first light instrument IRIS (Infrared Imaging Spectrograph) will complete its preliminary design phase in 2016. The IRIS instrument design includes a near-infrared (0.85 - 2.4 micron) integral field spectrograph (IFS) and imager that are able to conduct simultaneous diffraction-limited observations behind the advanced adaptive optics system NFIRAOS. The IRIS science cases have continued to be developed and new science studies have been investigated to aid in technical performance and design requirements. In this development phase, the IRIS science team has paid particular attention to the selection of filters, gratings, sensitivities of the entire system, and science cases that will benefit from the parallel mode of the IFS and imaging camera. We present new science cases for IRIS using the latest end-to-end data simulator on the following topics: Solar System bodies, the Galactic center, active galactic nuclei (AGN), and distant gravitationally-lensed galaxies. We then briefly discuss the necessity of an advanced data management system and data reduction pipeline.

We test a particular theory of dark matter in which dark matter axions form ring "caustics" in the plane of the Milky Way against actual observations of Milky Way stars. According to this theory, cold, collisionless dark matter particles with angular momentum flow in and out of the Milky Way on sheets. These flows form caustic rings (at the positions of the rings, the density of the flow is formally infinite) at the locations of closest approach to the Galactic center. We show that the caustic ring dark matter theory reproduces a roughly logarithmic halo, with large perturbations near the rings. We show that the theory can reasonably match the observed rotation curve of the Milky Way. We explore the effects of the caustic rings on dwarf galaxy tidal disruption using N-body simulations. In particular, simulations of the Sagittarius dwarf galaxy tidal disruption in a caustic ring halo potential match observations of the trailing tidal tail as far as 90 kpc from the Galactic center; they do not, however, match the leading tidal tail. None of the caustic ring, NFW, or triaxial logarithmic halos fit all of the data. The source code for calculating the acceleration due to a caustic ring halo has been made publicly available in the NEMO Stellar Dynamics Toolbox and the Milkyway@home client repository.

We describe the design and first-light early science performance of the Shane Adaptive optics infraRed Camera-Spectrograph (ShARCS) on Lick Observatory's 3-m Shane telescope. Designed to work with the new ShaneAO adaptive optics system, ShARCS is capable of high-efficiency, diffraction-limited imaging and low-dispersion grism spectroscopy in J, H, and K-bands. ShARCS uses a HAWAII-2RG infrared detector, giving high quantum efficiency (>80%) and Nyquist sampling the diffraction limit in all three wavelength bands. The ShARCS instrument is also equipped for linear polarimetry and is sensitive down to 650 nm to support future visible-light adaptive optics capability. We report on the early science data taken during commissioning.

IRIS (InfraRed Imaging Spectrograph) is a first light near-infrared diffraction limited imager and integral field spectrograph being designed for the future Thirty Meter Telescope (TMT). IRIS is optimized to perform astronomical studies across a significant fraction of cosmic time, from our Solar System to distant newly formed galaxies (Barton et al. [1]). We present a selection of the innovative science cases that are unique to IRIS in the era of upcoming space and ground-based telescopes. We focus on integral field spectroscopy of directly imaged exoplanet atmospheres, probing fundamental physics in the Galactic Center, measuring 10^4 to 10^10 Msun supermassive black hole masses, resolved spectroscopy of young star-forming galaxies (1 < z < 5) and first light galaxies (6 < z < 12), and resolved spectroscopy of strong gravitational lensed sources to measure dark matter substructure. For each of these science cases we use the IRIS simulator (Wright et al. [2], Do et al. [3]) to explore IRIS capabilities. To highlight the unique IRIS capabilities, we also update the point and resolved source sensitivities for the integral field spectrograph (IFS) in all five broadband filters (Z, Y, J, H, K) for the finest spatial scale of 0.004" per spaxel. We briefly discuss future development plans for the data reduction pipeline and quicklook software for the IRIS instrument suite.

We present an overview of the design of IRIS, an infrared (0.84 - 2.4 micron) integral field spectrograph and imaging camera for the Thirty Meter Telescope (TMT). With extremely low wavefront error (<30 nm) and on-board wavefront sensors, IRIS will take advantage of the high angular resolution of the narrow field infrared adaptive optics system (NFIRAOS) to dissect the sky at the diffraction limit of the 30-meter aperture. With a primary spectral resolution of 4000 and spatial sampling starting at 4 milliarcseconds, the instrument will create an unparalleled ability to explore high redshift galaxies, the Galactic center, star forming regions and virtually any astrophysical object. This paper summarizes the entire design and basic capabilities. Among the design innovations is the combination of lenslet and slicer integral field units, new 4Kx4k detectors, extremely precise atmospheric dispersion correction, infrared wavefront sensors, and a very large vacuum cryogenic system.

The Gemini Planet Imager (GPI) is a new facility instrument for the Gemini Observatory designed to provide direct detection and characterization of planets and debris disks around stars in the solar neighborhood. In addition to its extreme adaptive optics and corona graphic systems which give access to high angular resolution and high-contrast imaging capabilities, GPI contains an integral field spectrograph providing low resolution spectroscopy across five bands between 0.95 and 2.5 $\mu$m. This paper describes the sequence of processing steps required for the spectro-photometric calibration of GPI science data, and the necessary calibration files. Based on calibration observations of the white dwarf HD 8049B we estimate that the systematic error in spectra extracted from GPI observations is less than 5%. The flux ratio of the occulted star and fiducial satellite spots within coronagraphic GPI observations, required to estimate the magnitude difference between a target and any resolved companions, was measured in the $H$-band to be $\Delta m = 9.23\pm0.06$ in laboratory measurements and $\Delta m = 9.39\pm 0.11$ using on-sky observations. Laboratory measurements for the $Y$, $J$, $K1$ and $K2$ filters are also presented. The total throughput of GPI, Gemini South and the atmosphere of the Earth was also measured in each photometric passband, with a typical throughput in $H$-band of 18% in the non-coronagraphic mode, with some variation observed over the six-month period for which observations were available. We also report ongoing development and improvement of the data cube extraction algorithm.

The Gemini Planet Imager is a newly commissioned facility instrument designed to measure the near-infrared spectra of young extrasolar planets in the solar neighborhood and obtain imaging polarimetry of circumstellar disks. GPI's science instrument is an integral field spectrograph that utilizes a HAWAII-2RG detector with a SIDECAR ASIC readout system. This paper describes the detector characterization and calibrations performed by the GPI Data Reduction Pipeline to compensate for effects including bad/hot/cold pixels, persistence, non-linearity, vibration induced microphonics and correlated read noise.

The Gemini Planet Imager (GPI) is a complex optical system designed to directly detect the self-emission of young planets within two arcseconds of their host stars. After suppressing the starlight with an advanced AO system and apodized coronagraph, the dominant residual contamination in the focal plane are speckles from the atmosphere and optical surfaces. Since speckles are diffractive in nature their positions in the field are strongly wavelength dependent, while an actual companion planet will remain at fixed separation. By comparing multiple images at different wavelengths taken simultaneously, we can freeze the speckle pattern and extract the planet light adding an order of magnitude of contrast. To achieve a bandpass of 20%, sufficient to perform speckle suppression, and to observe the entire two arcsecond field of view at diffraction limited sampling, we designed and built an integral field spectrograph with extremely low wavefront error and almost no chromatic aberration. The spectrograph is fully cryogenic and operates in the wavelength range 1 to 2.4 microns with five selectable filters. A prism is used to produce a spectral resolution of 45 in the primary detection band and maintain high throughput. Based on the OSIRIS spectrograph at Keck, we selected to use a lenslet-based spectrograph to achieve an rms wavefront error of approximately 25 nm. Over 36,000 spectra are taken simultaneously and reassembled into image cubes that have roughly 192x192 spatial elements and contain between 11 and 20 spectral channels. The primary dispersion prism can be replaced with a Wollaston prism for dual polarization measurements. The spectrograph also has a pupil-viewing mode for alignment and calibration.

The Gemini Planet Imager (GPI) is an "extreme" adaptive optics coronagraph system that is now on the Gemini South telescope in Chile. This instrument is composed of three different systems that historically have been separate instruments. These systems are the extreme Adaptive Optics system, with deformable mirrors, including a high-order 64x64 element MEMS system; the Science Instrument, which is a near-infrared integral field spectrograph; and the Calibration system, a precision IR wavefront sensor that also holds key coronagraph components. Each system coordinates actions that require precise timing. The observatory is responsible for starting these actions and has typically done this asynchronously across independent systems. Despite this complexity we strived to provide an interface that is as close to a one-button approach as possible. This paper will describe the sequencing of these systems both internally and externally through the observatory.

We present H- and K-band data from the inner arcsecond of the Seyfert 1.5 galaxy NGC 4151 obtained with the adaptive optics assisted near-infrared imaging field spectrograph OSIRIS at the Keck Observatory. We present the morphology and dynamics of most species detected but focus on the morphology and dynamics of the narrow line region (as traced by emission of [FeII] \lambda 1.644\mu m), the interplay between plasma ejected from the nucleus (as traced by 21 cm continuum radio data) and hot H2 gas and characterize the detected nuclear HeI$\lambda 2.058 \mu m absorption feature as a narrow absorption line (NAL) phenomenon.

The Iapetus -1:0 nodal bending wave, the first spiral wave ever described in Saturn's rings, has been seen again for the first time in 29 years. We demonstrate that it is in fact the nodal bending wave, not the 1:0 apsidal density wave as previously reported. We use wavelet analysis to determine the wavelength profile, thus deriving the surface density at every point in the region covered by the bending wave. This profile is consistent with surface densities measured from more localized spiral density waves in the outer Cassini Division and the inner and mid-A Ring, varying smoothly from the low values of the former to the higher values of the latter. Most remarkably, our analysis indicates that there is no significant change in surface density across the boundary between the outer Cassini Division and the inner-A ring, despite the very abrupt increase in optical depth and reflected brightness at this location. We consider anew the nature of the classically identified "inner edge of the A ring," given that it does not appear to be correlated with any abrupt increase in surface density. There is an abrupt increase in surface density at the Pandora 5:4 density wave, ~300 km outward of the A ring's inner edge. Further study is needed to robustly interpret our findings in terms of particle properties and abundances, much less to explain the origins of the implied structure.

We report the discovery of several large "propeller" moons in the outer part of Saturn's A ring, objects large enough to be followed over the 5-year duration of the Cassini mission. These are the first objects ever discovered that can be tracked as individual moons, but do not orbit in empty space. We infer sizes up to 1--2 km for the unseen moonlets at the center of the propeller-shaped structures, though many structural and photometric properties of propeller structures remain unclear. Finally, we demonstrate that some propellers undergo sustained non-keplerian orbit motion. (Note: This arXiv version of the paper contains supplementary tables that were left out of the ApJL version due to lack of space).

We have developed a semi-analytic method of parameterizing N-body simulations of self-gravity wakes in Saturn's rings, describing their properties by means of only 6 numbers: 3 optical depths and 3 weighting factors. These numbers are obtained using a density-estimation procedure that finds the frequencies of various values of local density within a simulated ring patch. Application of our parameterization to a suite of N-body simulations implies that the distribution of local optical depths is trimodal, rather than bimodal as previous authors have assumed. Rings dominated by self-gravity wakes appear to be mostly empty space. The implications of this result for the analysis of occultation data are more conceptual than practical. The only adjustment needed is that the model parameter $\tau_{gap}$ should be interpreted as representing the area-weighted average optical depth within the gaps (or inter-wake regions). The most significant consequence of our results applies to the question of why "propeller" structures observed in the mid-A ring are seen as relative-bright features, even though the most prominent features of simulated propellers are regions of relatively low density. We find preliminary quantitative support for the hypothesis that propellers would be bright if they involve a local and temporary disruption of self-gravity wakes, flooding the region with more "photometrically active" material (i.e., material that can contribute to the rings' local optical depth) even though the overall density is lower, and thus raising their apparent brightnesses in agreement with observations. We suggest that this mechanism be tested by future detailed numerical models.

We present 0.1" resolution near-infrared integral field spectroscopy of Halpha in a z=1.4781 star forming galaxy, Q2343-BM133. These observations were obtained with OSIRIS (OH Suppressing Infra-Red Imaging Spectrograph) using the W.M. Keck Observatory Laser Guide Star Adaptive Optics system. Halpha emission is resolved over a 0.8" (6.8 kpc) x 0.5" (4.3 kpc) region with a 0.1" spatial resolution. We find a global flux of 4.2+/-0.6x10^{-16} ergs s^{-1} cm^{-2}, and detect a spatially resolved velocity gradient of ~134 km s^{-1} across the galaxy and a global velocity dispersion of 73+/-9 km s^{-1}. An upper limit of NII/Halpha < 0.12 is inferred, which implies that this galaxy is not dominated by an active galactic nucleus and has a metallicity at or below 1/2 solar metallicity. We derive a star formation rate (SFR) of 47+/-6 Msun yr^{-1}, and a dereddened SFR of 66+/-9 Msun yr^{-1}. Two-dimensional kinematics for Q2343-BM133 fit well with an inclined-disk model, with which we estimate an enclosed mass of 4.3x10^{9} Msun within 5.5 kpc. A possible merger scenario is also presented, and can not be fully ruled out. We derive a virial mass of 1.1x10^{10} Msun for a disk geometry, using the observed velocity dispersion. We propose that Q2343-BM133 is currently at an early stage of disk formation at a look-back time of 9.3 Gyr.

We present commissioning data from the OSIRIS integral field spectrograph (IFS) on the Keck II 10 m telescope that demonstrate the utility of adaptive optics IFS spectroscopy in studying faint close-in sub-stellar companions in the haloes of bright stars. Our R~2000 J- and H-band spectra of the sub-stellar companion to the 1-10 Myr-old GQ Lup complement existing K-band spectra and photometry, and improve on the original estimate of its spectral type. We find that GQ Lup B is somewhat hotter (M6-L0) than reported in the discovery paper by Neuhauser and collaborators (M9-L4), mainly due to the surface-gravity sensitivity of the K-band spectral classification indices used by the discoverers. Spectroscopic features characteristic of low surface gravity objects, such as lack of alkali absorption and a triangular H-band continuum, are indeed prominent in our spectrum of GQ Lup B. The peculiar shape of the H-band continuum and the difference between the two spectral type estimates is well explained in the context of the diminishing strength of H2 collision induced absorption with decreasing surface gravity, as recently proposed for young ultra-cool dwarfs by Kirkpatrick and collaborators. Using our updated spectroscopic classification of GQ Lup B and a re-evaluation of the age and heliocentric distance of the primary, we perform a comparative analysis of the available sub-stellar evolutionary models to estimate the mass of the companion. We find that the mass of GQ Lup B is 0.010-0.040 Msun. Hence, it is unlikely to be a wide-orbit counterpart to the known radial-velocity extrasolar planets, whose masses are < 0.015 Msun. Instead, GQ Lup A/B is probably a member of a growing family of very low mass ratio widely separated binaries discovered through high-contrast imaging.

We present diffraction limited spectroscopic observations of an infrared flare associated with the radio source SgrA*. These are the first results obtained with OSIRIS, the new facility infrared imaging spectrograph for the Keck Observatory operated with the laser guide star adaptive optics system. After subtracting the spectrum of precursor emission at the location of Sgr A*, we find the flare has a spectral index of -2.6 +- 0.9. If we do not subtract the precursor light, then our spectral index is consistent with earlier observations by Ghez et al. (2005). All observations published so far suggest that the spectral index is a function of the flare's K-band flux.