We report the discovery of two ultra-faint stellar systems found in early data from the DECam Local Volume Exploration survey (DELVE). The first system, Centaurus I (DELVE J1238-4054), is identified as a resolved overdensity of old and metal-poor stars with a heliocentric distance of ${\rm D}_{\odot} = 116.3_{-0.6}^{+0.6}$ kpc, a half-light radius of $r_h = 2.3_{-0.3}^{+0.4}$ arcmin, an age of $\tau > 12.85$ Gyr, a metallicity of $Z = 0.0002_{-0.0002}^{+0.0001}$, and an absolute magnitude of $M_V = -5.55_{-0.11}^{+0.11}$ mag. This characterization is consistent with the population of ultra-faint satellites, and confirmation of this system would make Centaurus I one of the brightest recently discovered ultra-faint dwarf galaxies. Centaurus I is detected in Gaia DR2 with a clear and distinct proper motion signal, confirming that it is a real association of stars distinct from the Milky Way foreground; this is further supported by the clustering of blue horizontal branch stars near the centroid of the system. The second system, DELVE 1 (DELVE J1630-0058), is identified as a resolved overdensity of stars with a heliocentric distance of ${\rm D}_{\odot} = 19.0_{-0.6}^{+0.5} kpc$, a half-light radius of $r_h = 0.97_{-0.17}^{+0.24}$ arcmin, an age of $\tau = 12.5_{-0.7}^{+1.0}$ Gyr, a metallicity of $Z = 0.0005_{-0.0001}^{+0.0002}$, and an absolute magnitude of $M_V = -0.2_{-0.6}^{+0.8}$ mag, consistent with the known population of faint halo star clusters. Given the low number of probable member stars at magnitudes accessible with Gaia DR2, a proper motion signal for DELVE 1 is only marginally detected. We compare the spatial position and proper motion of both Centaurus I and DELVE 1 with simulations of the accreted satellite population of the Large Magellanic Cloud (LMC) and find that neither is likely to be associated with the LMC.

The population of Milky Way (MW) satellites contains the faintest known galaxies and thus provides essential insight into galaxy formation and dark matter microphysics. Here we combine a model of the galaxy--halo connection with newly derived observational selection functions based on searches for satellites in photometric surveys over nearly the entire high Galactic latitude sky. In particular, we use cosmological zoom-in simulations of MW-like halos that include realistic Large Magellanic Cloud (LMC) analogs to fit the position-dependent MW satellite luminosity function. We report decisive evidence for the statistical impact of the LMC on the MW satellite population due to an estimated $6\pm 2$ observed LMC-associated satellites, consistent with the number of LMC satellites inferred from Gaia proper-motion measurements, confirming the predictions of cold dark matter models for the existence of satellites within satellite halos. Moreover, we infer that the LMC fell into the MW within the last $2\ \rm{Gyr}$ at high confidence. Based on our detailed full-sky modeling, we find that the faintest observed satellites inhabit halos with peak virial masses below $3.2\times 10^{8}\ M_{\rm{\odot}}$ at $95\%$ confidence, and we place the first robust constraints on the fraction of halos that host galaxies in this regime. We predict that the faintest detectable satellites occupy halos with peak virial masses above $10^{6}\ M_{\rm{\odot}}$, highlighting the potential for powerful galaxy formation and dark matter constraints from future dwarf galaxy searches.

There is increasing evidence that a substantial fraction of Milky Way satellite galaxies align in a rotationally-supported plane of satellites, a rare configuration in cosmological simulations of galaxy formation. It has been suggested that other Milky Way substructures (namely young halo globular clusters and stellar/gaseous streams) similarly tend to align with this plane, accordingly dubbed the Vast Polar Structure (VPOS). Using systemic proper motions inferred from Gaia data, we find that globular cluster orbital poles are not clustered in the VPOS direction, though the population with the highest VPOS membership fraction is the young halo clusters (~30%). We additionally provide a current census of stellar streams, including new streams discovered using the Dark Energy Survey and Gaia datasets, and find that stellar stream normals are also not clustered in the direction of the VPOS normal. We also find that, based on orbit modeling, there is a likely association between NGC 3201 and the Gj\"{o}ll stellar stream and that, based on its orbital pole, NGC 4147 is likely not a Sagittarius globular cluster. That the Milky Way's accreted globular clusters and streams do not align in the same planar configuration as its satellites suggests that the plane of satellites is either a particularly stable orbital configuration or a population of recently accreted satellites. Neither of these explanations is particularly likely in light of other recent studies, leaving the plane of satellites problem as one of the more consequential open problems in galaxy formation and cosmology.

This work presents the first search for RR Lyrae stars (RRLs) in four of the ultra-faint systems imaged by the Dark Energy Survey (DES) using SOAR/Goodman and Blanco/DECam imagers. We have detected two RRLs in the field of Grus I, none in Kim 2, one in Phoenix II, and four in Grus II. With the detection of these stars, we accurately determine the distance moduli for these ultra-faint dwarf satellite galaxies; $\mu_0$=20.51$\pm$0.10 mag (D$_{\odot}$=127$\pm$6 kpc) for Grus I and $\mu_0$=20.01$\pm$0.10 mag (D$_{\odot}$=100$\pm$5 kpc) for Phoenix II. These measurements are larger than previous estimations by Koposov et al. 2015 and Bechtol et al. 2015, implying larger physical sizes; 5\% for Grus I and 33\% for Phoenix II. For Grus II, out of the four RRLs detected, one is consistent with being a member of the galactic halo (D$_\odot$=24$\pm$1 kpc, $\mu_0$=16.86$\pm$0.10 mag), another is at D$_\odot$=55$\pm$2 kpc ($\mu_0$=18.71$\pm$0.10 mag), which we associate with Grus II, and the two remaining at D$_\odot$=43$\pm$2 kpc ($\mu_0$=18.17$\pm$0.10 mag). Moreover, the appearance of a subtle red horizontal branch in the color-magnitude diagram of Grus II at the same brightness level of the latter two RRLs, which are at the same distance and in the same region, suggests that a more metal-rich system may be located in front of Grus II. The most plausible scenario is the association of these stars with the Chenab/Orphan Stream. Finally, we performed a comprehensive and updated analysis of the number of RRLs in dwarf galaxies. This allows us to predict that the method of finding new ultra-faint dwarf galaxies by using two or more clumped RRLs will work only for systems brighter than M$_V\sim-6$ mag.

We use the Fermi-LAT gamma-ray observatory to search for gamma-ray emission from four nearby, debris disk-hosting main sequence stars: $\tau$ Ceti, $\epsilon$ Eridani, Fomalhaut, and Vega. For three stars ($\tau$ Ceti, Fomalhaut, and Vega), we establish upper limits that are consistent with theoretical expectations. For $\epsilon$ Eridani, we find a possible spatially coincident source with a soft energy spectrum of $dN/dE \sim E^{-3.6}$. However, at this stage we are unable to rule out that this emission is due to a more extended feature in the diffuse background. In the interpretation that the emission is due to $\epsilon$ Eridani, the $> 100$ MeV gamma-ray luminosity is $\sim 10^{27}$ erg/s $\simeq 3\times 10^{-7}$ L$_\odot$, which is $\sim 10^{10}$ times the gamma-ray luminosity from the disk of the quiet Sun. We find $\lesssim 2 \sigma$ evidence of source variability over a $\sim 7$ year timescale. In the interpretation that the gamma-ray emission from $\epsilon$ Eridani itself, we consider two possible models: 1) cosmic-ray collisions with solid bodies in the debris disk which extends out $\sim$60 AU from the host star, and 2) emission from the stellar activity. For the former model, assuming a total disk mass consistent with infrared measurements, we find that the size distribution of bodies is steeper than expected for a collisional cascade state. If confirmed as being associated with $\epsilon$ Eridani, this would be the first indication of gamma-ray emission from the vicinity of a main sequence star other than the Sun.

We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $\beta$, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from $\beta\sim -2$ at $r\sim20$ kpc to $\beta\sim 0.5$ at $r\sim200$ kpc, indicating that satellites closer to the Galactic centre have tangentially-biased motions while those farther out have radially-biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the $\beta(r)$ profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the $\beta(r)$ profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the Milky Way system.

Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag) in a single exposure; over 5000 deg$^{2}$) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multiband light curves from the DES wide-field survey (median 4 observations in each of grizY over the first three years) pose a challenge for traditional techniques used to detect RRL. We present an empirically motivated and computationally efficient template fitting method to identify these variable stars using three years of DES data. When tested on DES light curves of previously classified objects in SDSS stripe 82, our algorithm recovers 89% of RRL periods to within 1% of their true value with 85% purity and 76% completeness. Using this method, we identify 5783 RRL candidates, ~31% of which are previously undiscovered. This method will be useful for identifying RRL in other sparse multiband data sets.

The Lunar Gateway is expected to be positioned on-orbit around the Moon or in a Halo orbit at the L2 Lagrange point. The proposed Lunar Gateway is a game-changer for enabling new, high-priority lunar science utilizing Cu-beSats and presents a refreshing new opportunity for utilization of these small spacecraft as explorers. In context, CubeSats are being stretched to their limits as interplanetary explorers. The main technological hurdles include high-bandwidth communications and reliable high delta-v propulsion. Advances in deep-space attitude determination and control has been made possible from the recent NASA JPL MarCO missions. Due to these limitations, CubeSats are primarily designed to be dropped-off from a larger mission. The limited mass and volume have required compromises of the onboard science instruments, longer wait times to send back science data to Earth, shorter mission durations or higher accepted risk. With the Lunar Gateway being planned to be closer to the Moon, it will provide significant savings for a propulsion system and provide a primary relay for communication apart from the DSN and enable tele-operated command/control. These three factors can simplify the mission enabling routine deployment of CubeSats into lunar orbit and enable surface missions. In this paper, we present preliminary designs of 2 CubeSat lunar landers that will explore the lunar pits, Mare Tranquilitatis and the remnant magnetic fields Reiner Gamma.

The recent discovery of gravitational waves (GWs) from stellar-mass binary black holes (BBHs) provided direct evidence of the existence of these systems. BBH lenses would have gravitational microlensing signatures that are distinct from single-lens signals. We apply Bayesian statistics to examine the distinguishability of BBH microlensing events from single-lens events under ideal observing conditions, using the photometric capabilities of the Korean Microlensing Telescope Network. Given one year of observations, a source star at the Galactic Centre, a GW150914-like BBH lens (total mass 65M$_\odot$, mass ratio 0.8) at half that distance, and an impact parameter of 0.4 Einstein radii, we find that binarity is detectable for BBHs with separations down to 0.0250 Einstein radii, which is nearly 3.5 times greater than the maximum separation for which such BBHs would merge within the age of the Universe. Microlensing searches are thus sensitive to more widely separated BBHs than GW searches, perhaps allowing the discovery of BBH populations produced in different channels of binary formation.

We use the distance probability density function (DPDF) formalism of Ellsworth-Bowers et al. (2013, 2015) to derive physical properties for the collection of 1,710 Bolocam Galactic Plane Survey (BGPS) version 2 sources with well-constrained distance estimates. To account for Malmquist bias, we estimate that the present sample of BGPS sources is 90% complete above 400 $M_\odot$ and 50% complete above 70 $M_\odot$. The mass distributions for the entire sample and astrophysically motivated subsets are generally fitted well by a lognormal function, with approximately power-law distributions at high mass. Power-law behavior emerges more clearly when the sample population is narrowed in heliocentric distance (power-law index $\alpha = 2.0\pm0.1$ for sources nearer than 6.5 kpc and $\alpha = 1.9\pm0.1$ for objects between 2 kpc and 10 kpc). The high-mass power-law indices are generally $1.85 \leq \alpha \leq 2.05$ for various subsamples of sources, intermediate between that of giant molecular clouds and the stellar initial mass function. The fit to the entire sample yields a high-mass power-law $\hat{\alpha} = 1.94_{-0.10}^{+0.34}$. Physical properties of BGPS sources are consistent with large molecular cloud clumps or small molecular clouds, but the fractal nature of the dense interstellar medium makes difficult the mapping of observational categories to the dominant physical processes driving the observed structure. The face-on map of the Galactic disk's mass surface density based on BGPS dense molecular cloud structures reveals the high-mass star-forming regions W43, W49, and W51 as prominent mass concentrations in the first quadrant. Furthermore, we present a 0.25-kpc resolution map of the dense gas mass fraction across the Galactic disk that peaks around 5%.