This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

The space environment around the Earth is populated by more than 130 million objects of 1 mm in size and larger, and future predictions shows that this amount is destined to increase, even if mitigation measures are implemented at a far better rate than today. These objects can hit and damage a spacecraft or its components. It is thus necessary to assess the risk level for a satellite during its mission lifetime. Few software packages perform this analysis, and most of them employ time-consuming ray-tracing methodology, where particles are randomly sampled from relevant distributions. In addition, they tend not to consider the risk associated with the secondary debris clouds. The paper presents the development of a vulnerability assessment model, which relies on a fully statistical procedure: the debris fluxes are directly used combining them with the concept of the vulnerable zone, avoiding the random sampling the debris fluxes. A novel methodology is presented to predict damage to internal components. It models the interaction between the components and the secondary debris cloud through basic geometric operations, considering mutual shielding and shadowing between internal components. The methodologies are tested against state-of-the-art software for relevant test cases, comparing results on external structures and internal components.

Many problems in contemporary astrophysics---from understanding the formation of black holes to untangling the chemical evolution of galaxies---rely on knowledge about binary stars. This, in turn, depends on discovery and characterization of binary companions for large numbers of different kinds of stars in different chemical and dynamical environments. Current stellar spectroscopic surveys observe hundreds of thousands to millions of stars with (typically) few observational epochs, which allows binary discovery but makes orbital characterization challenging. We use a custom Monte Carlo sampler (The Joker) to perform discovery and characterization of binary systems through radial-velocities, in the regime of sparse, noisy, and poorly sampled multi-epoch data. We use it to generate posterior samplings in Keplerian parameters for 232,531 sources released in APOGEE Data Release 16. Our final catalog contains 19,635 high-confidence close-binary (P < few years, a < few AU) systems that show interesting relationships between binary occurrence rate and location in the color-magnitude diagram. We find notable faint companions at high masses (black-hole candidates), at low masses (substellar candidates), and at very close separations (mass-transfer candidates). We also use the posterior samplings in a (toy) hierarchical inference to measure the long-period binary-star eccentricity distribution. We release the full set of posterior samplings for the entire parent sample of 232,531 stars. This set of samplings involves no heuristic "discovery" threshold and therefore can be used for myriad statistical purposes, including hierarchical inferences about binary-star populations and sub-threshold searches.

Among the mitigation measures introduced to cope with the space debris issue there is the de-orbiting of decommissioned satellites. Guidelines for re-entering objects call for a ground casualty risk no higher than 0.0001. To comply with this requirement, satellites can be designed through a design-for-demise philosophy. Still, a spacecraft designed to demise has to survive the debris-populated space environment for many years. The demisability and the survivability of a satellite can both be influenced by a set of common design choices such as the material selection, the geometry definition, and the position of the components. Within this context, two models have been developed to analyse the demise and the survivability of satellites. Given the competing nature of the demisability and the survivability, a multi-objective optimisation framework was developed, with the aim to identify trade-off solutions for the preliminary design of satellites. As the problem is nonlinear and involves the combination of continuous and discrete variables, classical derivative based approaches are unsuited and a genetic algorithm was selected instead. The genetic algorithm uses the developed demisability and survivability criteria as the fitness functions of the multi-objective algorithm. The paper presents a test case, which considers the preliminary optimisation of tanks in terms of material, geometry, location, and number of tanks for a representative Earth observation mission. The configuration of the external structure of the spacecraft is fixed. Tanks were selected because they are sensitive to both design requirements: they represent critical components in the demise process and impact damage can cause the loss of the mission because of leaking and ruptures. The results present the possible trade off solutions, constituting the Pareto front obtained from the multi-objective optimisation.

The paper is concerned with examining the effects that design-for-demise solutions can have not only on the demisability of components, but also on their survivability that is their capability to withstand impacts from space debris. First two models are introduced. A demisability model to predict the behaviour of spacecraft components during the atmospheric re-entry and a survivability model to assess the vulnerability of spacecraft structures against space debris impacts. Two indices that evaluate the level of demisability and survivability are also proposed. The two models are then used to study the sensitivity of the demisability and of the survivability indices as a function of typical design-for-demise options. The demisability and the survivability can in fact be influenced by the same design parameters in a competing fashion that is while the demisability is improved, the survivability is worsened and vice versa. The analysis shows how the design-for-demise solutions influence the demisability and the survivability independently. In addition, the effect that a solution has simultaneously on the two criteria is assessed. Results shows which, among the design-for-demise parameters mostly influence the demisability and the survivability. For such design parameters maps are presented, describing their influence on the demisability and survivability indices. These maps represent a useful tool to quickly assess the level of demisability and survivability that can be expected from a component, when specific design parameters are changed.

Stars of a common origin are thought to have similar, if not nearly identical, chemistry. Chemical tagging seeks to exploit this fact to identify Milky Way subpopulations through their unique chemical fingerprints. In this work, we compare the chemical abundances of dwarf stars in wide binaries to test the abundance consistency of stars of a common origin. Our sample of 31 wide binaries is identified from a catalog produced by cross-matching APOGEE stars with UCAC5 astrometry, and we confirm the fidelity of this sample with precision parallaxes from Gaia DR2. For as many as 14 separate elements, we compare the abundances between components of our wide binaries, finding they have very similar chemistry (typically within 0.1 dex). This level of consistency is more similar than can be expected from stars with different origins (which show typical abundance differences of 0.3-0.4 dex within our sample). For the best measured elements, Fe, Si, K, Ca, Mn, and Ni, these differences are reduced to 0.05-0.08 dex when selecting pairs of dwarf stars with similar temperatures. Our results suggest that APOGEE dwarf stars may currently be used for chemical tagging at the level of $\sim$0.1 dex or at the level of $\sim$0.05 dex when restricting for the best-measured elements in stars of similar temperatures. Larger wide binary catalogs may provide calibration sets, in complement to open cluster samples, for on-going spectroscopic surveys.

We report the first detection of a credible progenitor system for a Type Ic supernova (SN Ic), SN 2017ein. We present spectra and photometry of the SN, finding it to be similar to carbon-rich, low-luminosity SNe Ic. Using a post-explosion Keck adaptive optics image, we precisely determine the position of SN 2017ein in pre-explosion \hst\ images, finding a single source coincident with the SN position. This source is marginally extended, and is consistent with being a stellar cluster. However, under the assumption that the emission of this source is dominated by a single point source, we perform point-spread function photometry, and correcting for line-of-sight reddening, we find it to have $M_{\rm F555W} = -7.5\pm0.2$ mag and $m_{\rm F555W}-m_{\rm F814W}$=$-0.67\pm0.14$ mag. This source is bluer than the main sequence and brighter than almost all Wolf-Rayet stars, however it is similar to some WC+O- and B-star binary systems. Under the assumption that the source is dominated by a single star, we find that it had an initial mass of $55\substack{+20-15} M_{\odot}$. We also examined binary star models to look for systems that match the overall photometry of the pre-explosion source and found that the best-fitting model is a $80$+$48 M_{\odot}$ close binary system in which the $80 M_{\odot}$ star is stripped and explodes as a lower mass star. Late-time photometry after the SN has faded will be necessary to cleanly separate the progenitor star emission from the additional coincident emission.

We use the multi-epoch radial velocities acquired by the APOGEE survey to perform a large scale statistical study of stellar multiplicity for field stars in the Milky Way, spanning the evolutionary phases between the main sequence and the red clump. We show that the distribution of maximum radial velocity shifts (\drvm) for APOGEE targets is a strong function of \logg, with main sequence stars showing \drvm\ as high as $\sim$300 \kms, and steadily dropping down to $\sim$30 \kms\ for \logg$\sim$0, as stars climb up the Red Giant Branch (RGB). Red clump stars show a distribution of \drvm\ values comparable to that of stars at the tip of the RGB, implying they have similar multiplicity characteristics. The observed attrition of high \drvm\ systems in the RGB is consistent with a lognormal period distribution in the main sequence and a multiplicity fraction of 0.35, which is truncated at an increasing period as stars become physically larger and undergo mass transfer after Roche Lobe Overflow during H shell burning. The \drvm\ distributions also show that the multiplicity characteristics of field stars are metallicity dependent, with metal-poor ([Fe/H]$\lesssim-0.5$) stars having a multiplicity fraction a factor 2-3 higher than metal-rich ([Fe/H]$\gtrsim0.0$) stars. This has profound implications for the formation rates of interacting binaries observed by astronomical transient surveys and gravitational wave detectors, as well as the habitability of circumbinary planets.

Motivated by the ever increasing pursuit of science with the transient sky (dubbed Time Domain Astronomy or TDA), we are fabricating and will commission a new deployable tertiary mirror for the Keck I telescope (K1DM3) at the W.M. Keck Observatory. This paper presents the detailed design of K1DM3 with emphasis on the opto-mechanics. This project has presented several design challenges. Foremost are the competing requirements to avoid vignetting the light path when retracted against a sufficiently rigid system for high-precision and repeatable pointing. The design utilizes an actuated swing arm to retract the mirror or deploy it into a kinematic coupling. The K1DM3 project has also required the design and development of custom connections to provide power, communications, and compressed air to the system. This NSF-MRI funded project is planned to be commissioned in Spring 2017.

ESA and NASA maintain asteroid hazard lists that contain all known asteroids with a non zero chance of colliding with the Earth in the future. Some software tools exist that are, either, capable of calculating the impact points of those asteroids, or that can estimate the impact effects of a given impact incident. However, no single tool is available that combines both aspects and enables a comprehensive risk analysis. The question is, thus, whether tools that can calculate impact location may be used to obtain a qualitative understanding of the asteroid impact risk distribution. To answer this question, two impact risk distributions that control for impact effect modelling were generated and compared. The Asteroid Risk Mitigation Optimization and Research (ARMOR) tool, in conjunction with the freely available software OrbFit, was used to project the impact probabilities of listed asteroids with a minimum diameter of 30 m onto the surface of the Earth representing a random sample (15% of all objects) of the hazard list. The resulting 261 impact corridors were visualized on a global map. Furthermore, the impact corridors were combined with Earth population data to estimate the simplified risk (without impact effects) and advanced risk (with impact effects) associated with the direct asteroid impacts that each nation faces from present to 2100 based on this sample. The relationship between risk and population size was examined for the 40 most populous countries and it was apparent that population size is a good proxy for relative risk. The advanced and simplified risk distributions were compared and the alteration of the results based on the introduction of physical impact effects was discussed. Population remained a valid proxy for relative impact risk, but the inclusion of impact effects resulted in significantly different risks, especially when considered at the national level.

Asteroids that could collide with the Earth are listed on the publicly available Near-Earth object (NEO) hazard web sites maintained by the US and European space agencies (NASA and ESA). The impact probability distribution of 69 potentially threatening NEOs from these lists that produce 261 dynamically distinct impact instances, or Virtual Impactors (VIs), were calculated using the Asteroid Risk Mitigation and Optimization Research (ARMOR) tool in conjunction with OrbFit. ARMOR projected the impact probability of each VI onto the surface of the Earth as a spatial probability distribution. The projection considers orbit solution accuracy and the global impact probability. The method of ARMOR is introduced and the tool is validated against two asteroid-Earth collision cases with objects 2008 TC3 and 2014 AA. In the analysis, the natural distribution of impact corridors is contrasted against the impact probability distribution to evaluate the distributions conformity with the uniform impact distribution assumption. The distribution of impact corridors is based on the NEO population and orbital mechanics. The analysis shows that the distribution of impact corridors matches the common assumption of uniform impact distribution and the result extends the evidence base for the uniform assumption from qualitative analysis of historic impact events into the future in a quantitative way. This finding is confirmed in a parallel analysis of impact points belonging to a synthetic population of 10006 VIs. Taking into account the impact probabilities introduced significant variation into the results and the impact probability distribution, consequently, deviates markedly from uniformity. The concept of impact probabilities is a product of the asteroid observation and orbit determination technique and, thus, represents a man-made component that is largely disconnected from natural processes.