We present the discoveries of WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5, two low-temperature (1200$-$1400 K), high proper motion T-type subdwarfs. Both objects were discovered via their high proper motion ($>$0.5 arcsec yr$^{-1}$); WISEA J181006.18-101000.5 as part of the NEOWISE proper motion survey and WISEA J041451.67-585456.7 as part of the citizen science project Backyard Worlds; Planet 9. We have confirmed both as brown dwarfs with follow-up near-infrared spectroscopy. Their spectra and near-infrared colors are unique amongst known brown dwarfs, with some colors consistent with L-type brown dwarfs and other colors resembling those of the latest-type T dwarfs. While no forward model consistently reproduces the features seen in their near-infrared spectra, the closest matches suggest very low metallicities ([Fe/H] $\leq$ -1), making these objects likely the first examples of extreme subdwarfs of the T spectral class (esdT). WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5 are found to be part of a small population of objects that occupy the "substellar transition zone," and have the lowest masses and effective temperatures of all objects in this group.

The radial spatial distribution of low-mass satellites around a Milky Way (MW)-like host is an important benchmark for simulations of small-scale structure. The distribution is sensitive to the disruption of subhalos by the central disk and can indicate whether the disruption observed in simulations of MW analogs is artificial (i.e., numeric) or physical in origin. We consider a sample of 12 well-surveyed satellite systems of MW-like hosts in the Local Volume that are complete to $M_V<-9$ and within 150 projected kpc. We investigate the radial distribution of satellites and compare with $\Lambda$CDM cosmological simulations, including big-box cosmological simulations and high resolution zoom in simulations of MW sized halos. We find that the observed satellites are significantly more centrally concentrated than the simulated systems. Several of the observed hosts, including the MW, are $\sim2\sigma$ outliers relative to the simulated hosts in being too concentrated, while none of the observed hosts are less centrally concentrated than the simulations. This result is robust to different ways of measuring the radial concentration. We find that this discrepancy is more significant for bright, $M_V<-12$ satellites, suggestive that this is not the result of observational incompleteness. We argue that the discrepancy is possibly due to artificial disruption in the simulations, but, if so, this has important ramifications for what stellar to halo mass relation is allowed in the low-mass regime by the observed abundance of satellites.

Low-mass satellites around Milky Way (MW)-like galaxies are important probes of small scale structure and galaxy formation. However, confirmation of satellite candidates with distance measurements remains a key barrier to fast progress in the Local Volume (LV). We measure the surface brightness fluctuation (SBF) distances to recently cataloged candidate dwarf satellites around 10 massive hosts within $D<12$ Mpc to confirm association. The satellite systems of these hosts are complete and mostly cleaned of contaminants down to $M_g{\sim}-9$ to $-10$, within the area of the search footprints. Joining this sample with hosts surveyed to comparable or better completeness in the literature, we explore how well cosmological simulations combined with common stellar to halo mass relations (SHMR) match observed satellite luminosity functions in the classical satellite luminosity regime. Adopting a SHMR that matches hydrodynamic simulations, the predicted overall satellite abundance agrees well with the observations. The MW is remarkably typical in its luminosity function amongst LV hosts. Contrary to recent results, we find that the host-to-host scatter predicted by the model is in close agreement with the scatter between the observed systems, once the different masses of the observed systems are taken into account. However, we find significant evidence that the observed systems have more bright and fewer faint satellites than the SHMR model predicts, necessitating a higher normalization of the SHMR around halo masses of $10^{11}$ \msun\ than present in common SHMRs. These results demonstrate the utility of nearby satellite systems in inferring the galaxy-subhalo connection in the low-mass regime.

We describe an empirical, self-contained method to isolate faint, large-scale emission in imaging data of low spatial resolution. Multi-resolution filtering (MRF) uses independent data of superior spatial resolution and point source depth to create a model for all compact and high surface brightness objects in the field. This model is convolved with an appropriate kernel and subtracted from the low resolution image. The halos of bright stars are removed in a separate step and artifacts are masked. The resulting image only contains extended emission fainter than a pre-defined surface brightness limit. The method was developed for the Dragonfly Telephoto Array, which produces images that have excellent low surface brightness sensitivity but poor spatial resolution. We demonstrate the MRF technique using Dragonfly images of a satellite of the spiral galaxy M101, the tidal debris surrounding M51, two ultra-diffuse galaxies in the Coma cluster, and the galaxy NGC5907. As part of the analysis we present a newly-identified very faint galaxy in the filtered Dragonfly image of the M101 field. We also discuss variations of the technique for cases when no low resolution data are available (self-MRF and cross-MRF). The method is implemented in mrf, an open-source MIT licensed Python package.

We present a description of the Dragonfly Wide Field Survey (DWFS), a deep photometric survey of a wide area of sky. The DWFS covers 330 $\mathrm{deg}^2$ in the equatorial GAMA fields and the Stripe 82 fields in the SDSS $g$ and $r$ bands. It is carried out with the 48-lens Dragonfly Telephoto Array, a telescope that is optimized for the detection of low surface brightness emission. The main goal of the survey is to study the dwarf galaxy population beyond the Local Group. In this paper, we describe the survey design and show early results. We reach $1\sigma$ depths of $\mu_g\approx 31$ mag arcsec$^{-2}$ on arcminute scales and show that Milky Way satellites such as Sextans, Bootes, and Ursa Major should be detectable out to $D\gtrsim 10$ Mpc. We also provide an overview of the elements and operation of the 48-lens Dragonfly telescope and a detailed description of its data reduction pipeline. The pipeline is fully automated, with individual frames subjected to a rigorous series of quality tests. The sky subtraction is performed in two stages, ensuring that emission features with spatial scales up to $\sim 0.^{\circ}9 \times 0.^{\circ}6$ are preserved. The DWFS provides unparalleled sensitivity to low surface brightness features on arcminute scales.

We present an in-depth study of surface brightness fluctuations (SBFs) in low-luminosity stellar systems. Using the MIST models, we compute theoretical predictions for absolute SBF magnitudes in the LSST, HST, and proposed WFIRST filter systems. We compare our calculations to the observed SBF magnitudes of dwarf galaxies that have independent distance measurements from the tip of the red giant branch method. Consistent with previous studies, we find that single-age population models show excellent agreement with the observed SBF-color relation of low-mass galaxies with $0.5 \lesssim g - i \lesssim 0.9$. For bluer galaxies, the observed relation is better fit by models with composite stellar populations. To study SBF recovery from low-luminosity galaxies, we perform detailed image simulations in which we inject fully populated model galaxies into deep ground-based images from real observations. We demonstrate that measurements of SBF magnitudes from these simulated data correspond to the theoretical values with negligible bias ($\lesssim0.01$ mag). We then use the simulations to show that LSST will provide data of sufficient quality and depth to measure SBF distances with precisions of ${\sim}10$-20% to ultra-faint $\left(\mathrm{10^4 \leq M_\star/M_\odot \leq 10^5}\right)$ and low-mass classical ($\leq10^7$ M$_\odot$) dwarf galaxies out to ${\sim}4$ Mpc and ${\sim}25$ Mpc, respectively, within the first few years of its deep-wide-fast survey. Many systematic uncertainties remain, including an irreducible "sampling scatter" in the SBFs of ultra-faint dwarfs due to their undersampled stellar mass functions. We nonetheless conclude that SBFs in the new generation of wide-field imaging surveys have the potential to play a critical role in the efficient confirmation and characterization of dwarf galaxies in the nearby universe.

We present the first satellite system of the Large Binocular Telescope Satellites Of Nearby Galaxies Survey (LBT-SONG), a survey to characterize the close satellite populations of Large Magellanic Cloud to Milky Way-mass, star-forming galaxies in the Local Volume. In this paper, we describe our satellite finding and completeness measurement methodology and apply this framework to NGC 628, an isolated galaxy with ~1/4 the stellar mass of the Milky Way. We present two new dwarf satellite galaxy candidates: NGC 628 dwA, and dwB with $\text{M}_{\text{V}}$ = $-12.2$ and $-7.7$, respectively. NGC 628 dwA is a classical dwarf while NGC 628 dwB is a low-luminosity galaxy that appears to have been quenched after reionization. Completeness corrections indicate that the presence of these two satellites is consistent with CDM predictions. The satellite colors indicate that the galaxies are neither actively star-forming nor do they have the purely ancient stellar populations characteristic of ultrafaint dwarfs. Instead, and consistent with our previous work on the NGC 4214 system, they show signs of recent quenching, further indicating that environmental quenching can play a role in modifying satellite populations even for hosts smaller than the Milky Way.

Like massive galaxies, dwarf galaxies are expected to undergo major mergers with other dwarfs. However, the end state of these mergers and the role that merging plays in regulating dwarf star formation is uncertain. Using imaging from the Hyper Suprime-Cam Subaru Strategic program, we construct a sample of dwarf-dwarf mergers and examine the star formation and host properties of the merging systems. These galaxies are selected via an automated detection algorithm from a sample of 6875 spectroscopically selected isolated dwarf galaxies at $z<0.12$ and $\log(M_\star/M_\odot)<9.6$ from the Galaxy and Mass Assembly (GAMA) and Sloan Digital Sky Survey (SDSS) spectroscopic campaigns. We find a total tidal feature detection fraction of 3.29% (6.1% when considering only galaxies at $z<0.05$). The tidal feature detection fraction rises strongly as a function of star formation activity; 15%-20% of galaxies with extremely high H$\alpha$ equivalent width (H$\alpha$ EW > 250 Angstrom) show signs of tidal debris. Galaxies that host tidal debris are also systematically bluer than the average galaxy at fixed stellar mass. These findings extend the observed dwarf-dwarf merger sequence with a significant sample of dwarf galaxies, indicating that star formation triggered in mergers between dwarf galaxies continues after coalescence.

We present extremely deep imaging of the well-known edge-on galaxy NGC 4565 (the "Needle Galaxy"), observed as part of the Dragonfly Edge-on Galaxies Survey. Our imaging traces starlight over scales of ~100 kpc to surface brightnesses of ~32 mag/arcsec^2 in azimuthally averaged radial profiles. In narrow slice profiles along the major axis (with bin sizes ranging from 1.7 x 0.5 kpc to 1.7 x 7.8 kpc), we achieve a depth of ~29-30 mag/arcsec^2. We confirm the previously observed asymmetric disc truncation. In addition, the sharp northwest truncation turns over to a shallower component that coincides with a fan-like feature seen to wrap around the northwest disc limb. We propose that the fan may be a tidal ribbon, and qualitatively replicate the fan with simple simulations. Alternative explanations of the fan and the disc's asymmetry are considered as well. We conclude that we are most likely seeing evidence for accretion-based outer disk growth in NGC 4565.

We present the results of an extensive search for dwarf satellite galaxies around 10 primary host galaxies in the Local Volume (D$<$12 Mpc) using archival CFHT/MegaCam imaging data. The hosts span a wide range in properties, with stellar masses ranging from that of the LMC to ${\sim}3$ times that of the Milky Way (MW). The surveyed hosts are: NGC 1023, NGC 1156, NGC 2903, NGC 4258, NGC 4565, NGC 4631, NGC 5023, M51, M64, and M104. We detect satellite candidates using a consistent semi-automated detection algorithm that is optimized for the detection of low surface brightness objects. Depending on the host, our completeness limit is $M_g{\sim}-8$ to $-10$ (assuming the distance of the host). We detect objects with surface brightness down to $\mu_{0,g}{\sim}26$ mag arcsec$^{-2}$ at $\gtrsim90\%$ completeness. The survey areas of the six best-surveyed hosts cover most of the inner projected $R<150$ kpc area, which roughly doubles the number of MW-mass hosts surveyed at this level of area and luminosity completeness. The number of detected candidates range from 1 around M64 to 33 around NGC 4258. In total, 153 candidates are found, of which 93 are new. While we defer an analysis of the satellite luminosity functions of the hosts until distance information is available for the candidates, we do show that the candidates are primarily red, spheroid systems with properties roughly consistent with known satellites in the Local Group.

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

We present low-resolution near-infrared spectra of discoveries from an all-sky proper motion search conducted using multi-epoch data from the Wide-field Infrared Survey Explorer. Using the data from NEOWISE, along with the AllWISE catalog, Schneider et al. (2016) conducted an all-sky proper motion survey to search for nearby objects with high proper motions. Here, we present a follow-up spectroscopic survey of 65 of their discoveries, which focused primarily on potentially nearby objects (d < 25 pc), candidate late-type brown dwarfs (>L7), and subdwarf candidates. We found 31 new M dwarfs, 18 new L dwarfs, and 11 new T dwarfs. Of these, 13 are subdwarfs, including one new sdL1 and two new sdL7s. Eleven of these discovereies, with spectral types ranging from M7 to T7 (including one subdwarf) are predicted to be within 25 pc, adding to the number of known objects in the solar neighborhood. We also discovered three new early-type T subdwarf candidates, one sdT1, one sdT2, and one sdT3, which would increase the number of known early-type T subdwarfs from two to five.

The advent of new deep+wide photometric lensing surveys will open up the possibility of direct measurements of the dark matter halos of dwarf galaxies. The HSC wide survey will be the first with the statistical capability of measuring the lensing signal with high signal-to-noise at log(M*)=8. At this same mass scale, LSST will have the most overall constraining power with a predicted signal-to-noise for the galaxy-galaxy lensing signal around dwarfs of SN=200. WFIRST and LSST will have the greatest potential to push below the log(M*) = 7 mass scale thanks to the depth of their imaging data. Studies of the dark matter halos of dwarf galaxies at z=0.1 with gravitational lensing are soon within reach. However, further work will be required to develop optimized strategies for extracting dwarfs samples from these surveys, determining redshifts, and accurately measuring lensing on small radial scales. Dwarf lensing will be a new and powerful tool to constrain the halo masses and inner density slopes of dwarf galaxies and to distinguish between baryonic feedback and modified dark matter scenarios.

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

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

We explore the use of ground-based surface brightness fluctuation (SBF) measurements to constrain distances to nearby dwarf galaxies. Using archival CFHT Megacam imaging data for a sample of 27 nearby dwarfs, we demonstrate that reliable SBF measurements and distances accurate to 15\% are possible even for very low surface brightness (LSB, $\mu_{i0}>24$ mag/arcsec$^2$) galaxies with modest, $\sim$hour-long exposures with CFHT. Combining our sample with a recent sample of 7 dwarfs with SBF measured with HST from the literature, we provide the most robust empirical SBF calibration to-date for the blue colors expected for these low mass systems. Our calibration is credible over the color range $0.3\lesssim g-i\lesssim0.8$ mag. It is also the first SBF calibration tied completely to TRGB distances as each galaxy in the sample has a literature TRGB distance. We find that even though the intrinsic scatter in SBF increases for blue galaxies, the rms scatter in the calibration is still $\lesssim0.3$ mag. We verify our measurements by comparing with HST SBF measurements and detailed image simulations. We argue that ground-based SBF is a very useful tool for characterizing dwarf satellite systems and field dwarfs in the nearby, D$\lesssim$20 Mpc universe.

We use surface brightness fluctuation (SBF) measurements to constrain the distance to low surface brightness (LSB) dwarfs in the vicinity of M101. Recent work has discovered many LSB candidate satellite companions of M101. However, without accurate distances, it is problematic to identify these dwarfs as physical satellites of M101. We use CFHT Legacy Survey (CFHTLS) data to measure the SBF signal for 43 candidate dwarfs. The data is deep enough that we constrain 29 of these to be unassociated background galaxies by their lack of SBF. We measure high S/N SBF signals for two of the candidate dwarfs, which are consistent with being at the distance of M101. The remaining candidates are too LSB and/or small for their distances to be constrained. Still, by comparison with Local Group dwarfs, we argue that the M101 satellite system is likely now complete down to stellar masses of $\sim5\times10^5$ M$_\odot$. We also provide a new SBF distance for the nearby dwarf UGC 8882, which suggests that it is significantly outside of the virial radius of M101 and is thus not a physical satellite. By constraining the distances to a majority of the candidates using only archival data, our work demonstrates the usefulness of SBF for nearby LSB galaxies and for studying the satellite systems of nearby massive galaxies.

Shock revival in core-collapse supernovae (CCSNe) may be due to the neutrino mechanism. While it is known that in a neutrino-powered CCSN, explosion begins when the neutrino luminosity of the proto-neutron star exceeds a critical value, the physics of this condition in time-dependent, multidimensional simulations are not fully understood. \citet{Pejcha2012} found that an `antesonic condition' exists for time-steady spherically symmetric models, potentially giving a physical explanation for the critical curve observed in simulations. In this paper, we extend that analysis to time-dependent, spherically symmetric polytropic models. We verify the critical antesonic condition in our simulations, showing that models exceeding it drive transonic winds whereas models below it exhibit steady accretion. In addition, we find that (1) high spatial resolution is needed for accurate determination of the antesonic ratio and shock radius at the critical curve, and that low resolution simulations systematically underpredict these quantities, making explosion more difficult at lower resolution; (2) there is an important physical connection between the critical mass accretion rate at explosion and the mass loss rate of the post-explosion wind: the two are directly proportional at criticality, implying that, at criticality, the wind kinetic power is tied directly to the accretion power; (3) the value of the post-shock adiabatic index $\Gamma$ has a large effect on the length and time scales of the post-bounce evolution of the explosion larger values of $\Gamma$ result in a longer transition from the accretion to wind phases.

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.

We present optical long-slit spectroscopy and far-ultraviolet to near-infrared spectral energy distribution fitting of two diffuse dwarf galaxies, LSBG-285 and LSBG-750, which were recently discovered by the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure redshifts using H$\alpha$ line emission, and find that these galaxies are at comoving distances of ${\approx}25$ and ${\approx}41$ Mpc, respectively, after correcting for the local velocity field. They have effective radii of $r_\mathrm{eff}=1.2$ and 1.8 kpc and stellar masses of $M_\star\approx2$-$3\times10^{7}~M_\odot$. There are no massive galaxies ($M_\star>10^{10} M_\odot$) within a comoving separation of at least 1.5 Mpc from LSBG-285 and 2 Mpc from LSBG-750. These sources are similar in size and surface brightness to ultra-diffuse galaxies, except they are isolated, star-forming objects that were optically selected in an environmentally blind survey. Both galaxies likely have low stellar metallicities $[Z_\star/Z_\odot] < -1.0$ and are consistent with the stellar mass-metallicity relation for dwarf galaxies. We set an upper limit on LSBG-750's rotational velocity of ${\sim}50$ km s$^{-1}$, which is comparable to dwarf galaxies of similar stellar mass with estimated halo masses $<10^{11}~M_\odot$. We find tentative evidence that the gas-phase metallicities in both of these diffuse systems are high for their stellar mass, though a statistically complete, optically-selected galaxy sample at very low surface brightness will be necessary to place these results into context with the higher-surface-brightness galaxy population.

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.

We present a catalog of extended low-surface-brightness galaxies (LSBGs) identified in the Wide layer of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). Using the first ${\sim}$200 deg$^2$ of the survey, we have uncovered 781 LSBGs, spanning red ($g-i\geq0.64$) and blue ($g-i<0.64$) colors and a wide range of morphologies. Since we focus on extended galaxies ($r_\mathrm{eff}=2.5$-$14^{\prime\prime}$), our sample is likely dominated by low-redshift objects. We define LSBGs to have mean surface brightnesses $\bar{\mu}_\mathrm{eff}(g)>24.3$ mag arcsec$^{-2}$, which allows nucleated galaxies into our sample. As a result, the central surface brightness distribution spans a wide range of $\mu_0(g)=18$-$27.4$ mag arcsec$^{-2}$, with 50% and 95% of galaxies fainter than 24.3 and 22 mag arcsec$^{-2}$, respectively. Furthermore, the surface brightness distribution is a strong function of color, with the red distribution being much broader and generally fainter than that of the blue LSBGs, and this trend shows a clear correlation with galaxy morphology. Red LSBGs typically have smooth light profiles that are well-characterized by single-component S\'{e}rsic functions. In contrast, blue LSBGs tend to have irregular morphologies and show evidence for ongoing star formation. We crossmatch our sample with existing optical, HI, and ultraviolet catalogs to gain insight into the physical nature of the LSBGs. We find that our sample is diverse, ranging from dwarf spheroidals and ultra-diffuse galaxies in nearby groups to gas-rich irregulars to giant LSB spirals, demonstrating the potential of the HSC-SSP to provide a truly unprecedented view of the LSBG population.

The relation between the mass and integrated electron pressure of galaxy group and cluster halos can be probed by stacking maps of the thermal Sunyaev-Zel'dovich (tSZ) effect. Perhaps surprisingly, recent observational results have indicated that the scaling relation between integrated pressure and mass follows the prediction of simple, self-similar models down to halo masses as low as $10^{12.5} \, M_{\odot}$. Hydrodynamical simulations that incorporate energetic feedback processes suggest that gas should be depleted from such low-mass halos, thus decreasing their tSZ signal relative to self-similar predictions. Here, we build on the modeling of Vikram, Lidz, and Jain (2017) to evaluate the bias in the interpretation of stacked tSZ measurements due to the signal from correlated halos (the "two-halo" term), which has generally been neglected in the literature. We fit theoretical models to a measurement of the tSZ -- galaxy group cross-correlation function, accounting explicitly for the one- and two- halo contributions. We find moderate evidence of a deviation from self-similarity in the pressure -- mass relation, even after marginalizing over conservative miscentering effects. We explore pressure -- mass models with a break at $10^{14} \, M_{\odot}$, as well as other variants. We discuss and test for sources of uncertainty in our analysis, in particular a possible bias in the halo mass estimates and the coarse resolution of the Planck beam. We compare our findings with earlier analyses by exploring the extent to which halo isolation criteria can reduce the two-halo contribution. Finally, we show that ongoing third-generation CMB experiments will explicitly resolve the one-halo term in low-mass groups; our methodology can be applied to these upcoming data sets to obtain a clear answer to the question of self-similarity and an improved understanding of hot gas in low-mass halos.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

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

Collisions and interactions between gas-rich galaxies are thought to be pivotal stages in their formation and evolution, causing the rapid production of new stars, and possibly serving as a mechanism for fueling supermassive black holes (BH). Harnessing the exquisite spatial resolution (~0.5 arcsec) afforded by the first ~170 deg^2 of the Hyper Suprime-Cam (HSC) Survey, we present our new constraints on the importance of galaxy-galaxy major mergers (1:4) in growing BHs throughout the last ~8 Gyrs. Utilizing mid-infrared observations in the WISE All-Sky survey, we robustly select active galactic nuclei (AGN) and mass-matched control galaxy samples, totaling ~140,000 spectroscopically confirmed systems at i<22 mag. We identify galaxy interaction signatures using a novel machine-learning random forest decision tree technique allowing us to select statistically significant samples of major-mergers, minor-mergers/irregular-systems, and non-interacting galaxies. We use these samples to show that galaxies undergoing mergers are a factor ~2-7 more likely to contain luminous obscured AGN than non-interacting galaxies, and this is independent of both stellar mass and redshift to z < 0.9. Furthermore, based on our comparison of AGN fractions in mass-matched samples, we determine that the most luminous AGN population (L_AGN > 10^45 erg/s) systematically reside in merging systems over non-interacting galaxies. Our findings show that galaxy-galaxy interactions do, on average, trigger luminous AGN activity substantially more often than in secularly evolving non-interacting galaxies, and we further suggest that the BH growth rate may be closely tied to the dynamical time of the merger system.

We study 379 central and 159 satellite early-type galaxies with two-dimensional kinematics from the integral-field survey Mapping Nearby Galaxies at APO (MaNGA) to determine how their angular momentum content depends on stellar and halo mass. Using the Yang et. al. (2007) group catalog, we identify central and satellite galaxies in groups with halo masses in the range 10^12.5 h^-1 M_sun < M_200b < 10^15 h^-1 M_sun. As in previous work, we see a sharp dependence on stellar mass, in the sense that ~ 70% of galaxies with stellar mass M_* > 10^11 h^-2 M_sun tend to have very little rotation, while nearly all galaxies at lower mass show some net rotation. The ~ 30% of high-mass galaxies that have significant rotation do not stand out in other galaxy properties except for a higher incidence of ionized gas emission. Our data are consistent with recent simulation results suggesting that major merging and gas accretion have more impact on the rotational support of lower-mass galaxies. When carefully matching the stellar mass distributions, we find no residual differences in angular momentum content between satellite and central galaxies at the 20\% level. Similarly, at fixed mass, galaxies have consistent rotation properties across a wide range of halo mass. However, we find that errors in classification of centrals and satellites with group finders systematically lowers differences between satellite and central galaxies at a level that is comparable to current measurement uncertainties. To improve constraints, the impact of group finding methods will have to be forward modeled via mock catalogs.

We present the data reduction pipeline for CHARIS, a high-contrast integral-field spectrograph for the Subaru Telescope. The pipeline constructs a ramp from the raw reads using the measured nonlinear pixel response, and reconstructs the data cube using one of three extraction algorithms: aperture photometry, optimal extraction, or $\chi^2$ fitting. We measure and apply both a detector flatfield and a lenslet flatfield and reconstruct the wavelength- and position-dependent lenslet point-spread function (PSF) from images taken with a tunable laser. We use these measured PSFs to implement a $\chi^2$-based extraction of the data cube, with typical residuals of ~5% due to imperfect models of the undersampled lenslet PSFs. The full two-dimensional residual of the $\chi^2$ extraction allows us to model and remove correlated read noise, dramatically improving CHARIS' performance. The $\chi^2$ extraction produces a data cube that has been deconvolved with the line-spread function, and never performs any interpolations of either the data or the individual lenslet spectra. The extracted data cube also includes uncertainties for each spatial and spectral measurement. CHARIS' software is parallelized, written in Python and Cython, and freely available on github with a separate documentation page. Astrometric and spectrophotometric calibrations of the data cubes and PSF subtraction will be treated in a forthcoming paper.

The "kinematic" morphology-density relation for early-type galaxies posits that those galaxies with low angular momentum are preferentially found in the highest-density regions of the universe. We use a large sample of galaxy groups with halo masses 10^12.5 < M_halo < 10^14.5 M_sun/h observed with the Mapping Nearby Galaxies at APO (MaNGA) survey to examine whether there is a correlation between local environment and rotational support that is independent of stellar mass. We find no compelling evidence for a relationship between the angular momentum content of early-type galaxies and either local overdensity or radial position within the group at fixed stellar mass.

We report the discovery of a diffuse stellar cloud with an angular extent $\gtrsim30^{\prime\prime}$, which we term "Sumo Puff", in data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). While we do not have a redshift for this object, it is in close angular proximity to a post-merger galaxy at redshift $z=0.0431$ and is projected within a few virial radii (assuming similar redshifts) of two other ${\sim}L_\star$ galaxies, which we use to bracket a potential redshift range of $0.0055 < z < 0.0431$. The object's light distribution is flat, as characterized by a low Sersic index ($n\sim0.3$). It has a low central $g$-band surface brightness of ${\sim}26.4$ mag arcsec$^{-2}$, large effective radius of ${\sim}13^{\prime\prime}$ (${\sim}11$ kpc at $z=0.0431$ and ${\sim}1.5$ kpc at $z=0.0055$), and an elongated morphology ($b/a\sim0.4$). Its red color ($g-i\sim1$) is consistent with a passively evolving stellar population and similar to the nearby post-merger galaxy, and we may see tidal material connecting Sumo Puff with this galaxy. We offer two possible interpretations for the nature of this object: (1) it is an extreme, galaxy-size tidal feature associated with a recent merger event, or (2) it is a foreground dwarf galaxy with properties consistent with a quenched, disturbed ultra-diffuse galaxy. We present a qualitative comparison with simulations that demonstrates the feasibility of forming a structure similar to this object in a merger event. Follow-up spectroscopy and/or deeper imaging to confirm the presence of the bridge of tidal material will be necessary to reveal the true nature of this object.

{\iota} Ori is a well studied massive binary consisting of an O9 III + B1 III/IV star. Due to its high eccentricity (e = 0.764) and short orbital period (P orb = 29.13376 d), it has been considered to be a good candidate to show evidence of tidal effects; however, none have previously been identified. Using photometry from the BRITE-Constellation space photometry mission we have confirmed the existence of tidal distortions through the presence of a heartbeat signal at periastron. We combine spectroscopic and light curve analyses to measure the masses and radii of the components, revealing {\iota} Ori to be the most massive heartbeat system known to date. In addition, using a thorough frequency analysis, we also report the unprecedented discovery of multiple tidally induced oscillations in an O star. The amplitudes of the pulsations allow us to empirically estimate the tidal circularization rate, yielding an effective tidal quality factor Q $\approx 4 \times 10^{4}$ .

The late-type Be star $\beta$ CMi is remarkably stable compared to other Be stars that have been studied. This has led to a realistic model of the outflowing Be disk by Klement et al. These results showed that the disk is likely truncated at a finite radius from the star, which Klement et al.~suggest is evidence for an unseen binary companion in orbit. Here we report on an analysis of the Ritter Observatory spectroscopic archive of $\beta$ CMi to search for evidence of the elusive companion. We detect periodic Doppler shifts in the wings of the H$\alpha$ line with a period of 170 d and an amplitude of 2.25 km s$^{-1}$, consistent with a low-mass binary companion ($M\approx 0.42 M_\odot$). We then compared the small changes in the violet-to-red peak height changes ($V/R$) with the orbital motion. We find weak evidence that it does follow the orbital motion, as suggested by recent Be binary models by Panoglou et al. Our results, which are similar to those for several other Be stars, suggest that $\beta$ CMi may be a product of binary evolution where Roche lobe overflow has spun up the current Be star, likely leaving a hot subdwarf or white dwarf in orbit around the star. Unfortunately, no direct sign of this companion star is found in the very limited archive of {\it International Ultraviolet Explorer} spectra.

The recovery of an exoplanet's atmospheric parameters from its spectrum requires accurate knowledge of the spectral errors and covariances. Unfortunately, the complex image processing used in high-contrast integral-field spectrograph (IFS) observations generally produces spectral covariances that are poorly understood and often ignored. In this work, we show how to measure the spectral errors and covariances and include them self-consistently in parameter retrievals. By combining model exoplanet spectra with a realistic noise model generated from GPI early science data, we show that ignoring spectral covariance in high-contrast IFS data can both bias inferred parameters and lead to unreliable confidence regions on those parameters. This problem is made worse by the common practice of scaling the $\chi^2$ per degree of freedom to unity; the input parameters then fall outside the $95\%$ confidence regions in as many as ${\sim}80\%$ of noise realizations. The biases we observe can approach the typical levels of precision achieved in high-contrast spectroscopy. Accounting for realistic priors in fully Bayesian retrievals can also have a significant impact on the inferred parameters. Plausible priors on effective temperature and surface gravity can vary by an order of magnitude across the confidence regions appropriate for objects with weak age constraints; priors for objects with good age constraints are dominated by modeling uncertainties. Our methods are directly applicable to existing high-contrast IFSs including GPI and SPHERE, as well as upcoming instruments like CHARIS and, ultimately, WFIRST-AFTA.

We use the AllWISE Data Release to continue our search for WISE-detected motions. In this paper, we publish another 27,846 motion objects, bringing the total number to 48,000 when objects found during our original AllWISE motion survey are included. We use this list, along with the lists of confirmed WISE-based motion objects from the recent papers by Luhman and by Schneider et al. and candidate motion objects from the recent paper by Gagne et al. to search for widely separated, common-proper-motion systems. We identify 1,039 such candidate systems. All 48,000 objects are further analyzed using color-color and color-mag plots to provide possible characterizations prior to spectroscopic follow-up. We present spectra of 172 of these, supplemented with new spectra of 23 comparison objects from the literature, and provide classifications and physical interpretations of interesting sources. Highlights include: (1) the identification of three G/K dwarfs that can be used as standard candles to study clumpiness and grain size in nearby molecular clouds because these objects are currently moving behind the clouds, (2) the confirmation/discovery of several M, L, and T dwarfs and one white dwarf whose spectrophotometric distance estimates place them 5-20 pc from the Sun, (3) the suggestion that the Na 'D' line be used as a diagnostic tool for interpreting and classifying metal-poor late-M and L dwarfs, (4) the recognition of a triple system including a carbon dwarf and late-M subdwarf, for which model fits of the late-M subdwarf (giving [Fe/H] ~ -1.0) provide a measured metallicity for the carbon star, and (5) a possible 24-pc-distant K5 dwarf + peculiar red L5 system with an apparent physical separation of 0.1 pc.

The Wide-field Infrared Survey Explorer (WISE) was reactivated in December of 2013 (NEOWISE) to search for potentially hazardous near-Earth objects. We have conducted a survey using the first sky pass of NEOWISE data and the AllWISE catalog to identify nearby stars and brown dwarfs with large proper motions ($\mu_{\rm total}$ $\gtrsim$ 250 mas yr$^{-1}$). A total of 20,548 high proper motion objects were identified, 1,006 of which are new discoveries. This survey has uncovered a significantly larger sample of fainter objects (W2 $\gtrsim$13 mag) than the previous WISE motion surveys of Luhman (2014a) and Kirkpatrick et al. (2014). Many of these objects are predicted to be new L and T dwarfs based on near- and mid-infrared colors. Using estimated spectral types along with distance estimates, we have identified several objects likely belonging to the nearby Solar neighborhood (d $<$ 25 pc). We have followed up 19 of these new discoveries with near-infrared or optical spectroscopy, focusing on potentially nearby objects, objects with the latest predicted spectral types, and potential late-type subdwarfs. This subset includes 6 M dwarfs, 5 of which are likely subdwarfs, as well as 8 L dwarfs and 5 T dwarfs, many of which have blue near-infrared colors. As an additional supplement, we provide 2MASS and AllWISE positions and photometry for every object found in our search, as well as 2MASS/AllWISE calculated proper motions.

We use the Planck full mission temperature maps to examine the stacked thermal Sunyaev-Zel'dovich (tSZ) signal of 188042 "locally brightest galaxies'" (LBGs) selected from the Sloan Digital Sky Survey Data Release 7. Our LBG sample closely matches that of Planck Collaboration XI (2013, PCXI) but our analysis differs in several ways. We work directly in terms of physically observable quantities, requiring minimal assumptions about the gas pressure profile. We explicitly model the dust emission from each LBG and simultaneously measure both the stacked tSZ and dust signals as a function of stellar mass $M_*$. There is a small residual bias in stacked tSZ measurements; we measure this bias and subtract it from our results, finding that the effects are non-negligible at the lowest masses in the LBG sample. Finally, we compare our measurements with two pressure profile models, finding that the profile from Battaglia et al. (2012b) provides a better fit to the results than the Arnaud et al. (2010) "universal pressure profile." However, within the uncertainties, we find that the data are consistent with a self-similar scaling with mass --- more precise measurements are needed to detect the relatively small deviations from self-similarity predicted by these models. Consistent with PCXI, we measure the stacked tSZ signal from LBGs with stellar masses down to $\log_{10}(M_*/M_{\odot}) \sim 11.1-11.3$. For lower stellar masses, however, we do not see evidence for a stacked tSZ signal. We note that the stacked dust emission is comparable to, or larger than, the stacked tSZ signal for $\log_{10}(M_*/M_{\odot}) \lesssim 11.3$. Future tSZ analyses with larger samples and lower noise levels should be able to probe deviations from self-similarity and thus provide constraints on models of feedback and the evolution of hot halo gas over cosmic time.

Motivated by the possibility that a coronagraph will be put on WFIRST/AFTA, we explore the direct detectability of extrasolar giant planets (EGPs) in the optical. We quantify a planet's detectability by the fraction of its orbit for which it is in an observable configuration ($f_\mathrm{obs}$). Using a suite of Monte Carlo experiments, we study the dependence of $f_\mathrm{obs}$ upon the inner working angle (IWA) and minimum achievable contrast ($C_\mathrm{min}$) of the direct-imaging observatory; the planet's phase function, geometric albedo, single-scattering albedo, radius, and distance from Earth; and the semi-major axis distribution of EGPs. We calculate phase functions for a given geometric or single-scattering albedo, assuming various scattering mechanisms. We find that the Lambertian phase function can predict significantly larger $f_\mathrm{obs}$'s with respect to the more realistic Rayleigh phase function. For observations made with WFIRST/AFTA's baseline capabilities ($C_\mathrm{min}\sim10^{-9}$, $\mathrm{IWA}\sim0.2''$), Jupiter-like planets orbiting stars within 10, 30, and 50 parsecs of Earth have volume-averaged observability fractions of ${\sim}$12%, 3%, and 0.5%, respectively. At 10 parsecs, such observations yield $f_\mathrm{obs}>1\%$ for low- to modest-eccentricity planets with semi-major axes in the range ${\sim}2 - 10$ AU. If $C_\mathrm{min}=10^{-10}$, this range extends to ${\sim}35$ AU. We find that, in all but the most optimistic configurations, the probability for detection in a blind search is low (${<}\,5\%$). However, with orbital parameter constraints from long-term radial-velocity campaigns and Gaia astrometry, the tools we develop in this work can be used to determine both the most promising systems to target and when to observe them.

We present a K-band spectroscopic study of the stellar and gas kinematics, mass distribution, and stellar populations of the archetypical starburst galaxy M82. Our results are based on a single spectrum at a position angle of 67.5 degrees through the K-band nucleus. We used the CO stellar absorption band head at 2.29 {\mu}m (CO_2.29) to measure the rotation curve out to nearly 4 kpc radius on both the eastern and western sides of the galaxy. Our data show that the rotation curve is flat from 1 - 4 kpc. This stands in sharp contrast to some previous studies, which have interpreted H I and CO emission-line position-velocity diagrams as evidence for a declining rotation curve. The kinematics of the Br\gamma, H_2, and He I emission lines are consistent with, although characterized by slightly higher velocities than, the stellar kinematics. We derived M82's mass distribution from our stellar kinematic measurements and estimate its total dynamical mass is ~10^10 Msun. We measured the equivalent width of CO_2.29 (W_2.29) as a function of distance from the center of the galaxy to investigate the spatial extent of the red supergiant (RSG) population. The variation in W_2.29 with radius clearly shows that RSGs dominate the light inside 500 pc radius. M82's superwind is likely launched from this region, where we estimate the enclosed mass is <= 2 x 10^9 Msun.