Type Ia Supernovae (SNe Ia) are powerful standardizable candles for constraining cosmological models and provided the first evidence of the accelerated expansion of the universe. Their precision derives from empirical correlations, now measured from $>1000$ SNe Ia, between their luminosities, light-curve shapes, colors and most recently with the stellar mass of their host galaxy. As mass correlates with other galaxy properties, alternative parameters have been investigated to improve SN Ia standardization though none have been shown to significantly alter the determination of cosmological parameters. We re-examine a recent claim, based on 34 SN Ia in nearby passive host galaxies, of a 0.05 mag/Gyr dependence of standardized SN Ia luminosity on host age which if extrapolated to higher redshifts, would be a bias up to 0.25 mag, challenging the inference of dark energy. We reanalyze this sample of hosts using both the original method and a Bayesian hierarchical model and find after a fuller accounting of the uncertainties the significance of a dependence on age to be $\leq2\sigma$ and $\sim1\sigma$ after the removal of a single poorly-sampled SN Ia. To test the claim that a trend seen in old stellar populations can be applied to younger ages, we extend our analysis to a larger sample which includes young hosts. We find the residual dependence of host age (after all standardization typically employed for cosmological measurements) to be consistent with zero for 254 SNe Ia from the Pantheon sample, ruling out the large but low significance trend seen in passive hosts.

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).

To determine if the SuperNova Empirical Model (SNEMO) can improve Type Ia supernova (SN Ia) standardization of several currently available photometric data sets, we perform an initial test, comparing results with the much-used SALT2 approach. We fit the SNEMO light-curve parameters and pass them to the Bayesian hierarchical model UNITY1.2 to estimate the Tripp-like standardization coefficients, including a host mass term as a proxy for redshift dependent astrophysical systematics. We find that, among the existing large data sets, only the Carnegie Supernova Project data set consistently provides the signal-to-noise and time sampling necessary to constrain the additional five parameters that SNEMO7 incorporates beyond SALT2. This is an important consideration for future SN Ia surveys like LSST and WFIRST. Although the SNEMO7 parameters are poorly constrained by most of the other available data sets of light curves, we find that the SNEMO2 parameters are just as well-constrained as the SALT2 parameters. In addition, SNEMO2 and SALT2 have comparable unexplained intrinsic scatter when fitting the same data. When looking at the total scatter, SNEMO7 reduces the Hubble-Lemaitre diagram RMS from 0.148~mag to 0.141~mag. It is not then, the SNEMO methodology, but the interplay of data quality and the increased number of degrees of freedom that is behind these reduced constraints. With this in mind, we recommend further investigation into the data required to use SNEMO7 and the possibility of fitting the poorer photometry data with intermediate SNEMO-like models with three to six components.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

Astrophysical transients have been observed for millennia and have shaped our most basic assumptions about the Universe. In the last century, systematic searches have grown from detecting handfuls of transients per year to over 7000 in 2018 alone. As these searches have matured, we have discovered both large samples of "normal" classes and new, rare classes. Recently, a transient was the first object observed in both gravitational waves and light. Ground-based observatories, including LSST, will discover thousands of transients in the optical, but these facilities will not provide the high-fidelity near-infrared (NIR) photometry and high-resolution imaging of a space-based observatory. WFIRST can fill this gap. With its survey designed to measure the expansion history of the Universe with Type Ia supernovae, WFIRST will also discover and monitor thousands of other transients in the NIR, revealing the physics for these high-energy events. Small-scale GO programs, either as a supplement to the planned survey or as specific target-of-opportunity observations, would significantly expand the scope of transient science that can be studied with WFIRST.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Opening up a new window of millimeter-wave observations that span frequency bands in the range of 30 to 500 GHz, survey half the sky, and are both an order of magnitude deeper (about 0.5 uK-arcmin) and of higher-resolution (about 10 arcseconds) than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. In particular, such a survey would allow for major advances in measuring the distribution of dark matter and gas on small-scales, and yield needed insight on 1.) dark matter particle properties, 2.) the evolution of gas and galaxies, 3.) new light particle species, 4.) the epoch of inflation, and 5.) the census of bodies orbiting in the outer Solar System.

The reliability of Type Ia supernovae (SNIa) may be limited by the imprint of their galactic origins. To investigate the connection between supernovae and their host characteristics, we developed an improved method to estimate the stellar population age of the host as well as the local environment around the site of the supernova. We use a Bayesian method to estimate the star formation history and mass weighted age of a supernova's environment by matching observed spectral energy distributions to a synthesized stellar population. Applying this age estimator to both the photometrically and spectroscopically classified Sloan Digital Sky Survey II supernovae (N=103) we find a $0.114 \pm 0.039~{\rm mag}$ `step' in the average Hubble residual at a stellar age of $\sim 8~\text{Gyr}$; it is nearly twice the size of the currently popular mass step. We then apply a principal component analysis on the SALT2 parameters, host stellar mass, and local environment age. We find that a new parameter, PC$_1$, consisting of a linear combination of stretch, host stellar mass, and local age, shows a very significant ($4.7\sigma$) correlation with Hubble residuals. There is a much broader range of PC$_1$ values found in the Hubble flow sample when compared with the Cepheid calibration galaxies. These samples have mildly statistically different average PC$_1$ values, at $\sim 2.5\sigma$, resulting in at most a 1.3% reduction in the evaluation of H$_0$. Despite accounting for the highly significant trend in SNIa Hubble residuals, there remains a 9% discrepancy between the most recent precision estimates of H$_0$ using SNIa and the CMB.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We re-analyze the detectability of large scale dark flow (or local bulk flow) with respect to the CMB background based upon the redshift-distance relation for Type Ia supernovae (SN Ia). We made two independent analyses: one based upon identifying the three Cartesian velocity components; and the other based upon the cosine dependence of the deviation from Hubble flow on the sky. We apply these analyses to the Union2.1 SN Ia data and to the SDSS-II supernova survey. For both methods, results for low redshift, $z < 0.05$, are consistent with previous searches. We find a local bulk flow of $v_{\rm bf} \sim 300$ km s$^{-1}$ in the direction of $(l,b) \sim (270, 35)^{\circ}$. However, the search for a dark flow at $z>0.05$ is inconclusive. Based upon simulated data sets, we deduce that the difficulty in detecting a dark flow at high redshifts arises mostly from the observational error in the distance modulus. Thus, even if it exists, a dark flow is not detectable at large redshift with current SN Ia data sets. We estimate that a detection would require both significant sky coverage of SN Ia out to $z = 0.3$ and a reduction in the effective distance modulus error from 0.2 mag to $\lesssim 0.02$ mag. We estimate that a greatly expanded data sample of $\sim 10^4$ SN Ia might detect a dark flow as small as 300 km s$^{-1}$ out to $z = 0.3$ even with a distance modulus error of $0.2$ mag. This may be achievable in a next generation large survey like LSST.

We present optical photometry and spectroscopy of the new eclipsing Cataclysmic Variable MASTER OTJ192328.22+612413.5, discovered by the MASTER team. We find the orbital period to be P=0.16764612(5) day /4.023507(1) hour. The depth of the eclipse (2.9$\pm$0.1 mag) suggests that the system is nearly edge on, and modeling of the system confirms the inclination to be between 81.3-83.6 degree. The brightness outside of eclipse varies between observations, with a change of 1.6$\pm$0.1 mag. Spectroscopy reveals double-peaked Balmer emission lines. By using spectral features matching a late M-type companion, we bound the distance to be 750$\pm$250 pc, depending on the companion spectral type. The source displays 2 mag brightness changes on timescales of days. The amplitude of these changes, along with the spectrum at the faint state, suggest the system is possibly a dwarf nova. The lack of any high excitation HeII lines suggests this system is not magnetically dominated. The light curve in both quiescence and outburst resembles that of Lanning 386, implying MASTER OTJ192328.22+612413.5 is a possible cross between a dwarf nova and a SW Sextantis star.

We summarize three recent efforts to constrain the first few moments of cosmic creation before and during the epoch of inflation. We consider two means to explain a slight dip in the power spectrum of the cosmic microwave background for multipoles in the range of $\ell= 10-30$ from both the {\it Planck} and {\it WMAP} data. We show that such a dip could be the result of resonant creation of a massive particle that couples to the inflaton field. For best-fit models, the epoch of resonant particle creation reenters the horizon at wave numbers of $k_* \sim 0.00011 \pm 0.0004 $ ($h$ Mpc$^{-1}$). The amplitude and location of these features correspond to the creation of a number of degenerate fermion species of mass $\sim 15/\lambda^{3/2} $ $m_{pl}$ during inflation where $\lambda$ is the coupling constant between the inflaton field and the created fermion species. Alternatively, one can explain the existence of such a dip as due to a jump in the inflation generating potential. We show that such a jump can also resolve the excessively large dark flow predicted from the M-theory landscape. Finally, we summarize our efforts to quantify constraints on the cosmic dark flow from a new analysis of the Type Ia supernova distance-redshift relation.

We present photometric and spectroscopic observations of the cataclysmic variable MASTER OT J132104.04+560957.8 which strongly indicate that it is a polar with an orbital period of 91 minutes. The optical light curve shows two maxima and two minima per orbital cycle, with considerable variation in the strength of the secondary maximum and in the morphology and depth of the minima.