The Wide Field InfraRed Survey Telescope (WFIRST) was the highest ranked large space-based mission of the 2010 New Worlds, New Horizons decadal survey. It is now a NASA mission in formulation with a planned launch in the mid-2020s. A primary mission objective is to precisely constrain the nature of dark energy through multiple probes, including Type Ia supernovae (SNe Ia). Here, we present the first realistic simulations of the WFIRST SN survey based on current hardware specifications and using open-source tools. We simulate SN light curves and spectra as viewed by the WFIRST wide-field channel (WFC) imager and integral-field channel (IFC) spectrometer, respectively. We examine 11 survey strategies with different time allocations between the WFC and IFC, two of which are based upon the strategy described by the WFIRST Science Definition Team, which measures SN distances exclusively from IFC data. We propagate statistical and, crucially, systematic uncertainties to predict the Dark Energy Task Force figure of merit (FoM) for each strategy. Of the strategies investigated, we find the most successful to be WFC-focused. However, further work in constraining systematics is required to fully optimize the use of the IFC. Even without improvements to other cosmological probes, the WFIRST SN survey has the potential to increase the FoM by more than an order of magnitude from the current values. Although the survey strategies presented here have not been fully optimized, these initial investigations are an important step in the development of the final hardware design and implementation of the WFIRST mission.

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

We present RAPID (Real-time Automated Photometric IDentification), a novel time-series classification tool capable of automatically identifying transients from within a day of the initial alert, to the full lifetime of a light curve. Using a deep recurrent neural network with Gated Recurrent Units (GRUs), we present the first method specifically designed to provide early classifications of astronomical time-series data, typing 12 different transient classes. Our classifier can process light curves with any phase coverage, and it does not rely on deriving computationally expensive features from the data, making RAPID well-suited for processing the millions of alerts that ongoing and upcoming wide-field surveys such as the Zwicky Transient Facility (ZTF), and the Large Synoptic Survey Telescope (LSST) will produce. The classification accuracy improves over the lifetime of the transient as more photometric data becomes available, and across the 12 transient classes, we obtain an average area under the receiver operating characteristic curve of 0.95 and 0.98 at early and late epochs, respectively. We demonstrate RAPID's ability to effectively provide early classifications of observed transients from the ZTF data stream. We have made RAPID available as an open-source software package (https://astrorapid.readthedocs.io) for machine learning-based alert-brokers to use for the autonomous and quick classification of several thousand light curves within a few seconds.

We describe the simulated data sample for the "Photometric LSST Astronomical Time Series Classification Challenge" (PLAsTiCC), a publicly available challenge to classify transient and variable events that will be observed by the Large Synoptic Survey Telescope (LSST), a new facility expected to start in the early 2020s. The challenge was hosted by Kaggle, ran from 2018 September 28 to 2018 December 17, and included 1,094 teams competing for prizes. Here we provide details of the 18 transient and variable source models, which were not revealed until after the challenge, and release the model libraries at https://doi.org/10.5281/zenodo.2612896. We describe the LSST Operations Simulator used to predict realistic observing conditions, and we describe the publicly available SNANA simulation code used to transform the models into observed fluxes and uncertainties in the LSST passbands (ugrizy). Although PLAsTiCC has finished, the publicly available models and simulation tools are being used within the astronomy community to further improve classification, and to study contamination in photometrically identified samples of type Ia supernova used to measure properties of dark energy. Our simulation framework will continue serving as a platform to improve the PLAsTiCC models, and to develop new models.

We present the analysis underpinning the measurement of cosmological parameters from 207 spectroscopically classified type Ia supernovae (SNe Ia) from the first three years of the Dark Energy Survey Supernova Program (DES-SN), spanning a redshift range of 0.017<$z$<0.849. We combine the DES-SN sample with an external sample of 122 low-redshift ($z$<0.1) SNe Ia, resulting in a "DES-SN3YR" sample of 329 SNe Ia. Our cosmological analyses are blinded: after combining our DES-SN3YR distances with constraints from the Cosmic Microwave Background (CMB; Planck Collaboration 2016), our uncertainties in the measurement of the dark energy equation-of-state parameter, $w$, are .042 (stat) and .059 (stat+syst) at 68% confidence. We provide a detailed systematic uncertainty budget, which has nearly equal contributions from photometric calibration, astrophysical bias corrections, and instrumental bias corrections. We also include several new sources of systematic uncertainty. While our sample is <1/3 the size of the Pantheon sample, our constraints on $w$ are only larger by 1.4$\times$, showing the impact of the DES SN Ia light curve quality. We find that the traditional stretch and color standardization parameters of the DES SNe Ia are in agreement with earlier SN Ia samples such as Pan-STARRS1 and the Supernova Legacy Survey. However, we find smaller intrinsic scatter about the Hubble diagram (0.077 mag). Interestingly, we find no evidence for a Hubble residual step ( 0.007 $\pm$ 0.018 mag) as a function of host galaxy mass for the DES subset, in 2.4$\sigma$ tension with previous measurements. We also present novel validation methods of our sample using simulated SNe Ia inserted in DECam images and using large catalog-level simulations to test for biases in our analysis pipelines.

We present the first cosmological parameter constraints using measurements of type Ia supernovae (SNe Ia) from the Dark Energy Survey Supernova Program (DES-SN). The analysis uses a subsample of 207 spectroscopically confirmed SNe Ia from the first three years of DES-SN, combined with a low-redshift sample of 122 SNe from the literature. Our "DES-SN3YR" result from these 329 SNe Ia is based on a series of companion analyses and improvements covering SN Ia discovery, spectroscopic selection, photometry, calibration, distance bias corrections, and evaluation of systematic uncertainties. For a flat LCDM model we find a matter density Omega_m = 0.331 +_ 0.038. For a flat wCDM model, and combining our SN Ia constraints with those from the cosmic microwave background (CMB), we find a dark energy equation of state w = -0.978 +_ 0.059, and Omega_m = 0.321 +_ 0.018. For a flat w0waCDM model, and combining probes from SN Ia, CMB and baryon acoustic oscillations, we find w0 = -0.885 +_ 0.114 and wa = -0.387 +_ 0.430. These results are in agreement with a cosmological constant and with previous constraints using SNe Ia (Pantheon, JLA).

We describe catalog-level simulations of Type Ia supernova (SN~Ia) light curves in the Dark Energy Survey Supernova Program (DES-SN), and in low-redshift samples from the Center for Astrophysics (CfA) and the Carnegie Supernova Project (CSP). These simulations are used to model biases from selection effects and light curve analysis, and to determine bias corrections for SN~Ia distance moduli that are used to measure cosmological parameters. To generate realistic light curves, the simulation uses a detailed SN~Ia model, incorporates information from observations (PSF, sky noise, zero point), and uses summary information (e.g., detection efficiency vs. signal to noise ratio) based on 10,000 fake SN light curves whose fluxes were overlaid on images and processed with our analysis pipelines. The quality of the simulation is illustrated by predicting distributions observed in the data. Averaging within redshift bins, we find distance modulus biases up to 0.05~mag over the redshift ranges of the low-z and DES-SN samples. For individual events, particularly those with extreme red or blue color, distance biases can reach 0.4~mag. Therefore, accurately determining bias corrections is critical for precision measurements of cosmological parameters. Files used to make these corrections are available at https://des.ncsa.illinois.edu/releases/sn.

The combination of multiple observational probes has long been advocated as a powerful technique to constrain cosmological parameters, in particular dark energy. The Dark Energy Survey has measured 207 spectroscopically--confirmed Type Ia supernova lightcurves; the baryon acoustic oscillation feature; weak gravitational lensing; and galaxy clustering. Here we present combined results from these probes, deriving constraints on the equation of state, $w$, of dark energy and its energy density in the Universe. Independently of other experiments, such as those that measure the cosmic microwave background, the probes from this single photometric survey rule out a Universe with no dark energy, finding $w=-0.80^{+0.09}_{-0.11}$. The geometry is shown to be consistent with a spatially flat Universe, and we obtain a constraint on the baryon density of $\Omega_b=0.069^{+0.009}_{-0.012}$ that is independent of early Universe measurements. These results demonstrate the potential power of large multi-probe photometric surveys and pave the way for order of magnitude advances in our constraints on properties of dark energy and cosmology over the next decade.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

Over the past century, major advances in astronomy and astrophysics have been largely driven by improvements in instrumentation and data collection. With the amassing of high quality data from new telescopes, and especially with the advent of deep and large astronomical surveys, it is becoming clear that future advances will also rely heavily on how those data are analyzed and interpreted. New methodologies derived from advances in statistics, computer science, and machine learning are beginning to be employed in sophisticated investigations that are not only bringing forth new discoveries, but are placing them on a solid footing. Progress in wide-field sky surveys, interferometric imaging, precision cosmology, exoplanet detection and characterization, and many subfields of stellar, Galactic and extragalactic astronomy, has resulted in complex data analysis challenges that must be solved to perform scientific inference. Research in astrostatistics and astroinformatics will be necessary to develop the state-of-the-art methodology needed in astronomy. Overcoming these challenges requires dedicated, interdisciplinary research. We recommend: (1) increasing funding for interdisciplinary projects in astrostatistics and astroinformatics; (2) dedicating space and time at conferences for interdisciplinary research and promotion; (3) developing sustainable funding for long-term astrostatisics appointments; and (4) funding infrastructure development for data archives and archive support, state-of-the-art algorithms, and efficient computing.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

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.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

The Wide Field Infrared Survey Telescope (WFIRST) is a 2.4m space telescope with a 0.281 deg^2 field of view for near-IR imaging and slitless spectroscopy and a coronagraph designed for > 10^8 starlight suppresion. As background information for Astro2020 white papers, this article summarizes the current design and anticipated performance of WFIRST. While WFIRST does not have the UV imaging/spectroscopic capabilities of the Hubble Space Telescope, for wide field near-IR surveys WFIRST is hundreds of times more efficient. Some of the most ambitious multi-cycle HST Treasury programs could be executed as routine General Observer (GO) programs on WFIRST. The large area and time-domain surveys planned for the cosmology and exoplanet microlensing programs will produce extraordinarily rich data sets that enable an enormous range of Archival Research (AR) investigations. Requirements for the coronagraph are defined based on its status as a technology demonstration, but its expected performance will enable unprecedented observations of nearby giant exoplanets and circumstellar disks. WFIRST is currently in the Preliminary Design and Technology Completion phase (Phase B), on schedule for launch in 2025, with several of its critical components already in production.

We analyze a set of 89 Type Ia supernovae (SN Ia) that have both optical and near-infrared (NIR) photometry to derive distances and construct low redshift ($z < 0.04$) Hubble diagrams. We construct mean light curve (LC) templates using a hierarchical Bayesian model. We explore both Gaussian process (GP) and template methods for fitting the LCs and estimating distances, while including peculiar velocity and photometric uncertainties. For the 56 SN Ia with both optical and NIR observations near maximum light, the GP method yields a NIR-only Hubble-diagram with a RMS of $0.117 \pm 0.014$ mag when referenced to the NIR maxima. For each NIR band, a comparable GP method RMS is obtained when referencing to NIR-max or B-max. Using NIR LC templates referenced to B-max yields a larger RMS value of $0.138 \pm 0.014$ mag. Fitting the corresponding optical data using standard LC fitters that use LC shape and color corrections yields larger RMS values of $0.179 \pm 0.018$ mag with SALT2 and $0.174 \pm 0.021$ mag with SNooPy. Applying our GP method to subsets of SN Ia NIR LCs at NIR maximum light, even without corrections for LC shape, color, or host-galaxy dust reddening, provides smaller RMS in the inferred distances, at the $\sim 2.3 - 4.1\sigma$ level, than standard optical methods that do correct for those effects. Our ongoing RAISIN program on the Hubble Space Telescope will exploit this promising infrared approach to limit systematic errors when measuring the expansion history of the universe to constrain dark energy.

We have established a network of 19 faint (16.5 mag $< V < $19 mag) northern and equatorial DA white dwarfs as spectrophotometric standards for present and future wide-field observatories. Our analysis infers SED models for the stars that are tied to the three CALSPEC primary standards. Our SED models are consistent with panchromatic Hubble Space Telescope ($HST$) photometry to better than 1%. The excellent agreement between observations and models validates the use of non-local-thermodynamic-equilibrium (NLTE) DA white dwarf atmospheres extinguished by interstellar dust as accurate spectrophotometric references. Our standards are accessible from both hemispheres and suitable for ground and space-based observatories covering the ultraviolet to the near infrared. The high-precision of these faint sources make our network of standards ideally suited for any experiment that has very stringent requirements on flux calibration, such as studies of dark energy using the Large Synoptic Survey Telescope (LSST) and the Wide-Field Infrared Survey Telescope ($WFIRST$).

The Wide-Field Infrared Survey Telescope (WFIRST) is expected to launch in the mid-2020s. With its wide-field near-infrared (NIR) camera, it will survey the sky to unprecedented detail. As part of normal operations and as the result of multiple expected dedicated surveys, WFIRST will produce several relatively wide-field (tens of square degrees) deep (limiting magnitude of 28 or fainter) fields. In particular, a planned supernova survey is expected to image 3 deep fields in the LSST footprint roughly every 5 days over 2 years. Stacking all data, this survey will produce, over all WFIRST supernova fields in the LSST footprint, ~12-25 deg^2 and ~5-15 deg^2 regions to depths of ~28 mag and ~29 mag, respectively. We suggest LSST undertake mini-surveys that will match the WFIRST cadence and simultaneously observe the supernova survey fields during the 2-year WFIRST supernova survey, achieving a stacked depth similar to that of the WFIRST data. We also suggest additional observations of these same regions throughout the LSST survey to get deep images earlier, have long-term monitoring in the fields, and produce deeper images overall. These fields will provide a legacy for cosmology, extragalactic, and transient/variable science.

We present spectroscopy from the first three seasons of the Dark Energy Survey Supernova Program (DES-SN). We describe the supernova spectroscopic program in full: strategy, observations, data reduction, and classification. We have spectroscopically confirmed 307 supernovae, including 251 type Ia supernovae (SNe Ia) over a redshift range of $0.017 < z < 0.85$. We determine the effective spectroscopic selection function for our sample, and use it to investigate the redshift-dependent bias on the distance moduli of SNe Ia we have classified. We also provide a full overview of the strategy, observations, and data products of DES-SN, which has discovered 12,015 likely supernovae during these first three seasons. The data presented here are used for the first cosmology analysis by DES-SN ('DES-SN3YR'), the results of which are given in DES Collaboration (2018a).

Empirical models of supernova (SN) spectral energy distributions (SEDs) are widely used for SN survey simulations and photometric classifications. The existing library of SED models has excellent optical templates but limited, poorly constrained coverage of ultraviolet (UV) and infrared (IR) wavelengths. However, both regimes are critical for the design and operation of future SN surveys, particularly at IR wavelengths that will be accessible with the James Webb Space Telescope (JWST) and the Wide-Field Infrared Survey Telescope (WFIRST). We create a public repository of improved empirical SED templates using a sampling of Type Ia and core-collapse (CC) photometric light curves to extend the Type Ia parameterized SALT2 model and a set of SN Ib, SN Ic, and SN II SED templates into the UV and near-IR. We apply this new repository of extrapolated SN SED models to examine how future surveys can discriminate between CC and Type Ia SNe at UV and IR wavelengths, and present an open-source software package written in Python, SNSEDextend, that enables a user to generate their own extrapolated SEDs.

The Photometric LSST Astronomical Time Series Classification Challenge (PLAsTiCC) is an open data challenge to classify simulated astronomical time-series data in preparation for observations from the Large Synoptic Survey Telescope (LSST), which will achieve first light in 2019 and commence its 10-year main survey in 2022. LSST will revolutionize our understanding of the changing sky, discovering and measuring millions of time-varying objects. In this challenge, we pose the question: how well can we classify objects in the sky that vary in brightness from simulated LSST time-series data, with all its challenges of non-representativity? In this note we explain the need for a data challenge to help classify such astronomical sources and describe the PLAsTiCC data set and Kaggle data challenge, noting that while the references are provided for context, they are not needed to participate in the challenge.

Classification of transient and variable light curves is an essential step in using astronomical observations to develop an understanding of their underlying physical processes. However, upcoming deep photometric surveys, including the Large Synoptic Survey Telescope (LSST), will produce a deluge of low signal-to-noise data for which traditional labeling procedures are inappropriate. Probabilistic classification is more appropriate for the data but are incompatible with the traditional metrics used on deterministic classifications. Furthermore, large survey collaborations intend to use these classification probabilities for diverse science objectives, indicating a need for a metric that balances a variety of goals. We describe the process used to develop an optimal performance metric for an open classification challenge that seeks probabilistic classifications and must serve many scientific interests. The Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) is an open competition aiming to identify promising techniques for obtaining classification probabilities of transient and variable objects by engaging a broader community both within and outside astronomy. Using mock classification probability submissions emulating archetypes of those anticipated of PLAsTiCC, we compare the sensitivity of metrics of classification probabilities under various weighting schemes, finding that they yield qualitatively consistent results. We choose as a metric for PLAsTiCC a weighted modification of the cross-entropy because it can be meaningfully interpreted. Finally, we propose extensions of our methodology to ever more complex challenge goals and suggest some guiding principles for approaching the choice of a metric of probabilistic classifications.

High precision proper motion (PM) measurements are available for approximately 20% of all known dwarf satellite galaxies of the Milky Way (MW). Here we extend the Bayesian framework of Patel et al. (2017b) to include all MW satellites with measured 6D phase space information and apply it with the Illustris-Dark simulation to constrain the MW's mass. Using the properties of each MW satellite individually, we find that the scatter among mass estimates is reduced when the magnitude of specific orbital angular momentum (j) is adopted rather than their combined instantaneous positions and velocities. We also find that high j satellites (i.e. Leo II) constrain the upper limits for the MW's mass and low j satellites rather than the highest speed satellites (i.e. Leo I and LMC), set the lower mass limits. When j of all classical satellites is used to simultaneously estimate the MW's mass, we conclude the halo mass is $0.85^{+0.23}_{-0.26}\times 10^{12} \rm \, M_{\odot}$ (including Sagittarius dSph) and $0.96^{+0.29}_{-0.28} \times 10^{12} \rm \, M_{\odot}$ (excluding Sagittarius dSph), cautioning that low j satellites on decaying orbits like Sagittarius dSph may bias the distribution. These estimates markedly reduce the current factor of two spread in the mass range of the MW. We also find a well-defined relationship between host halo mass and satellite j distribution, which yields the prediction that upcoming PMs for ultra-faint dwarfs should reveal j within $ 5\times10^3 - 10^4$ kpc km s$^{-1}$. This is a promising method to significantly constrain the cosmologically expected mass range for the MW and eventually M31 as more satellite PMs become available.

Much of the progress made in time-domain astronomy is accomplished by relating observational multi-wavelength time series data to models derived from our understanding of physical laws. This goal is typically accomplished by dividing the task in two: collecting data (observing), and constructing models to represent that data (theorizing). Owing to the natural tendency for specialization, a disconnect can develop between the best available theories and the best available data, potentially delaying advances in our understanding new classes of transients. We introduce MOSFiT: the Modular Open-Source Fitter for Transients, a Python-based package that downloads transient datasets from open online catalogs (e.g., the Open Supernova Catalog), generates Monte Carlo ensembles of semi-analytical light curve fits to those datasets and their associated Bayesian parameter posteriors, and optionally delivers the fitting results back to those same catalogs to make them available to the rest of the community. MOSFiT is designed to help bridge the gap between observations and theory in time-domain astronomy; in addition to making the application of existing models and creation of new models as simple as possible, MOSFiT yields statistically robust predictions for transient characteristics, with a standard output format that includes all the setup information necessary to reproduce a given result. As large-scale surveys such as LSST discover entirely new classes of transients, tools such as MOSFiT will be critical for enabling rapid comparison of models against data in statistically consistent, reproducible, and scientifically beneficial ways.

The rate of tidal disruption events (TDEs), $R_\text{TDE}$, is predicted to depend on stellar conditions near the super-massive black hole (SMBH), which are on difficult-to-measure sub-parsec scales. We test whether $R_\text{TDE}$ depends on kpc-scale global galaxy properties, which are observable. We concentrate on stellar surface mass density, $\Sigma_{M_\star}$, and velocity dispersion, $\sigma_v$, which correlate with the stellar density and velocity dispersion of the stars around the SMBH. We consider 35 TDE candidates, with and without known X-ray emission. The hosts range from star-forming to quiescent to quiescent with strong Balmer absorption lines. The last (often with post-starburst spectra) are overrepresented in our sample by a factor of $35^{+21}_{-17}$ or $18^{+8}_{-7}$, depending on the strength of the H$\delta$ absorption line. For a subsample of hosts with homogeneous measurements, $\Sigma_{M_\star}=10^9$-$10^{10}~{\rm M_\odot / kpc^2}$, higher on average than for a volume-weighted control sample of Sloan Digital Sky Survey galaxies with similar redshifts and stellar masses. This is because: (1) most of the TDE hosts here are quiescent galaxies, which tend to have higher $\Sigma_{M_\star}$ than the star-forming galaxies that dominate the control, and (2) the star-forming hosts have higher average $\Sigma_{M_\star}$ than the star-forming control. There is also a weak suggestion that TDE hosts have lower $\sigma_v$ than for the quiescent control. Assuming that $R_{\rm TDE}\propto \Sigma_{M_\star}^\alpha \times \sigma_v^\beta$, and applying a statistical model to the TDE hosts and control sample, we estimate $\hat{\alpha}=0.9 \pm 0.2$ and $\hat{\beta}=-1.0 \pm 0.6$. This is broadly consistent with $R_\text{TDE}$ being tied to the dynamical relaxation of stars surrounding the SMBH.

Space observatories like the Hubble Space Telescope and Gaia are providing unprecedented 6D phase space information of satellite galaxies. Such measurements can shed light on the structure and assembly history of the Local Group, but improved statistical methods are needed to use them efficiently. Here we illustrate such a method using analogs of the Local Group's two most massive satellite galaxies, the Large Magellanic Cloud (LMC) and Triangulum (M33), from the Illustris dark-matter-only cosmological simulation. We use a Bayesian inference scheme combining measurements of positions, velocities, and specific orbital angular momenta ( j ) of the LMC/M33 with importance sampling of their simulated analogs to compute posterior estimates of the Milky Way (MW) and Andromeda's (M31) halo masses. We conclude the resulting host halo mass is more susceptible to bias when using measurements of the current position and velocity of satellites, especially when satellites are at short-lived phases of their orbits (i.e. at pericentre). Instead, the j value of a satellite is well-conserved over time and provides a more reliable constraint on host mass. The inferred virial mass of the MW (M31) using j of the LMC (M33) is $\rm M_{vir, MW}=1.02^{+0.77}_{-0.55}\times10^{12}\; M_{\odot}$ ($\rm M_{vir, M31}=1.37^{+1.39}_{-0.75}\times10^{12}\; M_{\odot}$). Choosing simulated analogs whose j values are consistent with the conventional picture of a previous (< 3 Gyr ago), close encounter (< 100 kpc) of M33 about M31 results in a very low virial mass for M31 ($\rm\sim\!10^{12}\; M_{\odot}$). This supports the new scenario put forth in Patel et al. (2017), wherein M33 is on its first passage about M31 or on a long period orbit. We conclude that this Bayesian inference scheme, utilising satellite j, is a promising method to reduce the current factor of two spread in the mass range of the MW and M31 moving forward.

Measurements of radial velocity variations from the spectroscopic monitoring of stars and their companions are essential for a broad swath of astrophysics, providing access to the fundamental physical properties that dictate all phases of stellar evolution and facilitating the quantitative study of planetary systems. The conversion of those measurements into both constraints on the orbital architecture and individual component spectra can be a serious challenge, however, especially for extreme flux ratio systems and observations with relatively low sensitivity. Gaussian processes define sampling distributions of flexible, continuous functions that are well-motivated for modeling stellar spectra, enabling proficient search for companion lines in time-series spectra. We introduce a new technique for spectral disentangling, where the posterior distributions of the orbital parameters and intrinsic, rest-frame stellar spectra are explored simultaneously without needing to invoke cross-correlation templates. To demonstrate its potential, this technique is deployed on red-optical time-series spectra of the mid-M dwarf binary LP661-13. We report orbital parameters with improved precision compared to traditional radial velocity analysis and successfully reconstruct the primary and secondary spectra. We discuss potential applications for other stellar and exoplanet radial velocity techniques and extensions to time-variable spectra. The code used in this analysis is freely available as an open source Python package.

Conventional Type Ia supernova (SN Ia) cosmology analyses currently use a simplistic linear regression of magnitude versus color and light curve shape, which does not model intrinsic SN Ia variations and host galaxy dust as physically distinct effects, resulting in low color-magnitude slopes. We construct a probabilistic generative model for the dusty distribution of extinguished absolute magnitudes and apparent colors as the convolution of a intrinsic SN Ia color-magnitude distribution and a host galaxy dust reddening-extinction distribution. If the intrinsic color-magnitude ($M_B$ vs. $B-V$) slope $\beta_{int}$ differs from the host galaxy dust law $R_B$, this convolution results in a specific curve of mean extinguished absolute magnitude vs. apparent color. The derivative of this curve smoothly transitions from $\beta_{int}$ in the blue tail to $R_B$ in the red tail of the apparent color distribution. The conventional linear fit approximates this effective curve near the average apparent color, resulting in an apparent slope $\beta_{app}$ between $\beta_{int}$ and $R_B$. We incorporate these effects into a hierarchical Bayesian statistical model for SN Ia light curve measurements, and analyze a dataset of SALT2 optical light curve fits of 248 nearby SN Ia at z < 0.10. The conventional linear fit obtains $\beta_{app} \approx 3$. Our model finds a $\beta_{int} = 2.3 \pm 0.3$ and a distinct dust law of $R_B = 3.8 \pm 0.3$, consistent with the average for Milky Way dust, while correcting a systematic distance bias of $\sim 0.10$ mag in the tails of the apparent color distribution. Finally, we extend our model to examine the SN Ia luminosity-host mass dependence in terms of intrinsic and dust components.

The gravitational field of a galaxy can act as a lens and deflect the light emitted by a more distant object such as a quasar. Strong gravitational lensing causes multiple images of the same quasar to appear in the sky. Since the light in each gravitationally lensed image traverses a different path length from the quasar to the Earth, fluctuations in the source brightness are observed in the several images at different times. The time delay between these fluctuations can be used to constrain cosmological parameters and can be inferred from the time series of brightness data or light curves of each image. To estimate the time delay, we construct a model based on a state-space representation for irregularly observed time series generated by a latent continuous-time Ornstein-Uhlenbeck process. We account for microlensing, an additional source of independent long-term extrinsic variability, via a polynomial regression. Our Bayesian strategy adopts a Metropolis-Hastings within Gibbs sampler. We improve the sampler by using an ancillarity-sufficiency interweaving strategy and adaptive Markov chain Monte Carlo. We introduce a profile likelihood of the time delay as an approximation of its marginal posterior distribution. The Bayesian and profile likelihood approaches complement each other, producing almost identical results; the Bayesian method is more principled but the profile likelihood is simpler to implement. We demonstrate our estimation strategy using simulated data of doubly- and quadruply-lensed quasars, and observed data from quasars Q0957+561 and J1029+2623.

We report evidence for excess blue light from the Type Ia supernova SN 2012cg at fifteen and sixteen days before maximum B-band brightness. The emission is consistent with predictions for the impact of the supernova on a non-degenerate binary companion. This is the first evidence for emission from a companion to a SN Ia. Sixteen days before maximum light, the B-V color of SN 2012cg is 0.2 mag bluer than for other normal SN~Ia. At later times, this supernova has a typical SN Ia light curve, with extinction-corrected M_B = -19.62 +/- 0.02 mag and Delta m_{15}(B) = 0.86 +/- 0.02. Our data set is extensive, with photometry in 7 filters from 5 independent sources. Early spectra also show the effects of blue light, and high-velocity features are observed at early times. Near maximum, the spectra are normal with a silicon velocity v_{Si} = -10,500$ km s^{-1}. Comparing the early data with models by Kasen (2010) favors a main-sequence companion of about 6 solar masses. It is possible that many other SN Ia have main-sequence companions that have eluded detection because the emission from the impact is fleeting and faint.

We present a modular, extensible likelihood framework for spectroscopic inference based on synthetic model spectra. The subtraction of an imperfect model from a continuously sampled spectrum introduces covariance between adjacent datapoints (pixels) into the residual spectrum. For the high signal-to-noise data with large spectral range that is commonly employed in stellar astrophysics, that covariant structure can lead to dramatically underestimated parameter uncertainties (and, in some cases, biases). We construct a likelihood function that accounts for the structure of the covariance matrix, utilizing the machinery of Gaussian process kernels. This framework specifically address the common problem of mismatches in model spectral line strengths (with respect to data) due to intrinsic model imperfections (e.g., in the atomic/molecular databases or opacity prescriptions) by developing a novel local covariance kernel formalism that identifies and self-consistently downweights pathological spectral line "outliers." By fitting many spectra in a hierarchical manner, these local kernels provide a mechanism to learn about and build data-driven corrections to synthetic spectral libraries. An open-source software implementation of this approach is available at http://iancze.github.io/Starfish, including a sophisticated probabilistic scheme for spectral interpolation when using model libraries that are sparsely sampled in the stellar parameters. We demonstrate some salient features of the framework by fitting the high resolution $V$-band spectrum of WASP-14, an F5 dwarf with a transiting exoplanet, and the moderate resolution $K$-band spectrum of Gliese 51, an M5 field dwarf.

CfAIR2 is a large homogeneously reduced set of near-infrared (NIR) light curves for Type Ia supernovae (SN Ia) obtained with the 1.3m Peters Automated InfraRed Imaging TELescope (PAIRITEL). This data set includes 4607 measurements of 94 SN Ia and 4 additional SN Iax observed from 2005-2011 at the Fred Lawrence Whipple Observatory on Mount Hopkins, Arizona. CfAIR2 includes JHKs photometric measurements for 88 normal and 6 spectroscopically peculiar SN Ia in the nearby universe, with a median redshift of z~0.021 for the normal SN Ia. CfAIR2 data span the range from -13 days to +127 days from B-band maximum. More than half of the light curves begin before the time of maximum and the coverage typically contains ~13-18 epochs of observation, depending on the filter. We present extensive tests that verify the fidelity of the CfAIR2 data pipeline, including comparison to the excellent data of the Carnegie Supernova Project. CfAIR2 contributes to a firm local anchor for supernova cosmology studies in the NIR. Because SN Ia are more nearly standard candles in the NIR and are less vulnerable to the vexing problems of extinction by dust, CfAIR2 will help the supernova cosmology community develop more precise and accurate extragalactic distance probes to improve our knowledge of cosmological parameters, including dark energy and its potential time variation.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

We investigate the statistical dependence of the peak intrinsic colors of Type Ia supernovae (SN Ia) on their expansion velocities at maximum light, measured from the Si II 6355 spectral feature. We construct a new hierarchical Bayesian regression model, accounting for the random effects of intrinsic scatter, measurement error, and reddening by host galaxy dust, and implement a Gibbs sampler and deviance information criteria to estimate the correlation. The method is applied to the apparent colors from BVRI light curves and Si II velocity data for 79 nearby SNe Ia. The apparent color distributions of high (HV) and normal velocity (NV) supernovae exhibit significant discrepancies for B-V and B-R, but not other colors. Hence, they are likely due to intrinsic color differences originating in the B-band, rather than dust reddening. The mean intrinsic B-V and B-R color differences between HV and NV groups are 0.06 +/- 0.02 and 0.09 +/- 0.02 mag, respectively. A linear model finds significant slopes of -0.021 +/- 0.006 and -0.030 +/- 0.009 mag/(1000 km/s) for intrinsic B-V and B-R colors versus velocity, respectively. Since the ejecta velocity distribution is skewed towards high velocities, these effects imply non-Gaussian intrinsic color distributions with skewness up to +0.3. Accounting for the intrinsic color-velocity correlation results in corrections to A_V extinction estimates as large as -0.12 mag for HV SNe Ia and +0.06 mag for NV events. Velocity measurements from SN Ia spectra have potential to diminish systematic errors from the confounding of intrinsic colors and dust reddening affecting supernova distances.

We present a new method for probabilistically classifying supernovae (SNe) without using SN spectral or photometric data. Unlike all previous studies to classify SNe without spectra, this technique does not use any SN photometry. Instead, the method relies on host-galaxy data. We build upon the well-known correlations between SN classes and host-galaxy properties, specifically that core-collapse SNe rarely occur in red, luminous, or early-type galaxies. Using the nearly spectroscopically complete Lick Observatory Supernova Search sample of SNe, we determine SN fractions as a function of host-galaxy properties. Using these data as inputs, we construct a Bayesian method for determining the probability that a SN is of a particular class. This method improves a common classification figure of merit by a factor of >2, comparable to the best light-curve classification techniques. Of the galaxy properties examined, morphology provides the most discriminating information. We further validate this method using SN samples from the Sloan Digital Sky Survey and the Palomar Transient Factory. We demonstrate that this method has wide-ranging applications, including separating different subclasses of SNe and determining the probability that a SN is of a particular class before photometry or even spectra can. Since this method uses completely independent data from light-curve techniques, there is potential to further improve the overall purity and completeness of SN samples and to test systematic biases of the light-curve techniques. Further enhancements to the host-galaxy method, including additional host-galaxy properties, combination with light-curve methods, and hybrid methods should further improve the quality of SN samples from past, current, and future transient surveys.

We investigate the shape of the extinction law in two 1-degree square fields of the Perseus Molecular Cloud complex. We combine deep red-optical (r, i, and z-band) observations obtained using Megacam on the MMT with UKIDSS near-infrared (J, H, and K-band) data to measure the colours of background stars. We develop a new hierarchical Bayesian statistical model, including measurement error, intrinsic colour variation, spectral type, and dust reddening, to simultaneously infer parameters for individual stars and characteristics of the population. We implement an efficient MCMC algorithm utilising generalised Gibbs sampling to compute coherent probabilistic inferences. We find a strong correlation between the extinction (Av) and the slope of the extinction law (parameterized by Rv). Because the majority of the extinction toward our stars comes from the Perseus molecular cloud, we interpret this correlation as evidence of grain growth at moderate optical depths. The extinction law changes from the diffuse value of Rv = 3 to the dense cloud value of Rv = 5 as the column density rises from Av = 2 mags to Av = 10 mags. This relationship is similar for the two regions in our study, despite their different physical conditions, suggesting that dust grain growth is a fairly universal process.

In this first paper in a series we present 1298 low-redshift (z\leq0.2) optical spectra of 582 Type Ia supernovae (SNe Ia) observed from 1989 through 2008 as part of the Berkeley SN Ia Program (BSNIP). 584 spectra of 199 SNe Ia have well-calibrated light curves with measured distance moduli, and many of the spectra have been corrected for host-galaxy contamination. Most of the data were obtained using the Kast double spectrograph mounted on the Shane 3 m telescope at Lick Observatory and have a typical wavelength range of 3300-10,400 Ang., roughly twice as wide as spectra from most previously published datasets. We present our observing and reduction procedures, and we describe the resulting SN Database (SNDB), which will be an online, public, searchable database containing all of our fully reduced spectra and companion photometry. In addition, we discuss our spectral classification scheme (using the SuperNova IDentification code, SNID; Blondin & Tonry 2007), utilising our newly constructed set of SNID spectral templates. These templates allow us to accurately classify our entire dataset, and by doing so we are able to reclassify a handful of objects as bona fide SNe Ia and a few other objects as members of some of the peculiar SN Ia subtypes. In fact, our dataset includes spectra of nearly 90 spectroscopically peculiar SNe Ia. We also present spectroscopic host-galaxy redshifts of some SNe Ia where these values were previously unknown. [Abridged]

We present 2603 spectra of 462 nearby Type Ia supernovae (SN Ia) obtained during 1993-2008 through the Center for Astrophysics Supernova Program. Most of the spectra were obtained with the FAST spectrograph at the FLWO 1.5m telescope and reduced in a consistent manner, making data set well suited for studies of SN Ia spectroscopic diversity. We study the spectroscopic and photometric properties of SN Ia as a function of spectroscopic class using the classification schemes of Branch et al. and Wang et al. The width-luminosity relation appears to be steeper for SN Ia with broader lines. Based on the evolution of the characteristic Si II 6355 line, we propose improved methods for measuring velocity gradients, revealing a larger range than previously suspected, from ~0 to ~400 km/s/day considering the instantaneous velocity decline rate at maximum light. We find a weaker and less significant correlation between Si II velocity and intrinsic B-V color at maximum light than reported by Foley et al., owing to a more comprehensive treatment of uncertainties and host galaxy dust. We study the extent of nuclear burning and report new detections of C II 6580 in 23 early-time spectra. The frequency of C II detections is not higher in SN Ia with bluer colors or narrower light curves, in conflict with the recent results of Thomas et al. Based on nebular spectra of 27 SN Ia, we find no relation between the FWHM of the iron emission feature at ~4700 A and Dm15(B) after removing the two low-luminosity SN 1986G and SN 1991bg, suggesting that the peak luminosity is not strongly dependent on the kinetic energy of the explosion for most SN Ia. Finally, we confirm the correlation of velocity shifts in some nebular lines with the intrinsic B-V color of SN Ia at maximum light, although several outliers suggest a possible non-monotonic behavior for the largest blueshifts.

Eta Carinae (Eta Car) is one of the most massive binary stars in the Milky Way. It became the second-brightest star in the sky during its mid-19th century "Great Eruption," but then faded from view (with only naked-eye estimates of brightness). Its eruption is unique among known astronomical transients in that it exceeded the Eddington luminosity limit for 10 years. Because it is only 2.3 kpc away, spatially resolved studies of the nebula have constrained the ejected mass and velocity, indicating that in its 19th century eruption, Eta Car ejected more than 10 M_solar in an event that had 10% of the energy of a typical core-collapse supernova without destroying the star. Here we report the discovery of light echoes of Eta Carinae which appear to be from the 1838-1858 Great Eruption. Spectra of these light echoes show only absorption lines, which are blueshifted by -210 km/s, in good agreement with predicted expansion speeds. The light-echo spectra correlate best with those of G2-G5 supergiant spectra, which have effective temperatures of ~5000 K. In contrast to the class of extragalactic outbursts assumed to be analogs of Eta Car's Great Eruption, the effective temperature of its outburst is significantly cooler than allowed by standard opaque wind models. This indicates that other physical mechanisms like an energetic blast wave may have triggered and influenced the eruption.

We compare models for Type Ia supernova (SN Ia) light curves and spectra with an extensive set of observations. The models come from a recent survey of 44 two-dimensional delayed-detonation models computed by Kasen, Roepke & Woosley (2009), each viewed from multiple directions. The data include optical light curves of 251 SNe Ia and 2231 low-dispersion spectra from the Center for Astrophysics, plus data from the literature. The analysis uses standard techniques employed by observers, including MLCS2k2, SALT2, and SNooPy for light-curve analysis, and the Supernova Identification (SNID) code of Blondin & Tonry for spectroscopic comparisons to assess how well the models match the data. We show that the models that match observed spectra best lie systematically on the observed width-luminosity relation. Conversely, we reject six models with highly asymmetric ignition conditions and a large amount (>1 M_sun) of synthesized 56Ni that yield poor matches to observed SN Ia spectra. More subtle features of the comparison include the general difficulty of the models to match the U-band flux at early times, caused by a hot ionized ejecta that affect the subsequent redistribution of flux at longer wavelengths. We examine ways in which the asymptotic kinetic energy of the explosion affects both the predicted velocity and velocity gradient in the Si II and Ca II lines. Models with an asymmetric distribution of 56Ni are found to result in a larger variation of photometric and spectroscopic properties with viewing angle, regardless of the initial ignition setup. We discuss more generally whether highly anisotropic ignition conditions are ruled out by observations, and how detailed comparisons between models and observations involving both light curves and spectra can lead to a better understanding of SN Ia explosion mechanisms.

We have constructed a comprehensive statistical model for Type Ia supernova (SN Ia) light curves spanning optical through near infrared (NIR) data. A hierarchical framework coherently models multiple random and uncertain effects, including intrinsic supernova light curve covariances, dust extinction and reddening, and distances. An improved BayeSN MCMC code computes probabilistic inferences for the hierarchical model by sampling the global probability density of parameters describing individual supernovae and the population. We have applied this hierarchical model to optical and NIR data of 127 SN Ia from PAIRITEL, CfA3, CSP, and the literature. We find an apparent population correlation between the host galaxy extinction A_V and the the ratio of total-to-selective dust absorption R_V. For SN with low dust extinction, A_V < 0.4, we find R_V = 2.5 - 2.9, while at high extinctions, A_V > 1, low values of R_V < 2 are favored. The NIR luminosities are excellent standard candles and are less sensitive to dust extinction. They exhibit low correlation with optical peak luminosities, and thus provide independent information on distances. The combination of NIR and optical data constrains the dust extinction and improves the predictive precision of individual SN Ia distances by about 60%. Using cross-validation, we estimate an rms distance modulus prediction error of 0.11 mag for SN with optical and NIR data versus 0.15 mag for SN with optical data alone. Continued study of SN Ia in the NIR is important for improving their utility as precise and accurate cosmological distance indicators.

The light echo systems of historical supernovae in the Milky Way and local group galaxies provide an unprecedented opportunity to reveal the effects of asymmetry on observables, particularly optical spectra. Scattering dust at different locations on the light echo ellipsoid witnesses the supernova from different perspectives and the light consequently scattered towards Earth preserves the shape of line profile variations introduced by asymmetries in the supernova photosphere. However, the interpretation of supernova light echo spectra to date has not involved a detailed consideration of the effects of outburst duration and geometrical scattering modifications due to finite scattering dust filament dimension, inclination, and image point-spread function and spectrograph slit width. In this paper, we explore the implications of these factors and present a framework for future resolved supernova light echo spectra interpretation, and test it against Cas A and SN 1987A light echo spectra. We conclude that the full modeling of the dimensions and orientation of the scattering dust using the observed light echoes at two or more epochs is critical for the correct interpretation of light echo spectra. Indeed, without doing so one might falsely conclude that differences exist when none are actually present.

[Abridged] We investigate the use of a wide variety of spectroscopic measurements to determine distances to low-redshift Type Ia supernovae (SN Ia). We consider linear models for predicting distances to SN Ia using light-curve width and color parameters (determined using the SALT2 light-curve fitter) and a spectroscopic indicator, and evaluate the resulting Hubble diagram scatter using a cross-validation procedure. We confirm the ability of spectral flux ratios alone at maximum light to reduce the scatter of Hubble residuals by ~10% with respect to the standard combination of light-curve width and color. When used in combination with the SALT2 color parameter, the color-corrected flux ratio R^c(6420/5290) at maximum light leads to an even lower scatter, although the improvement has low statistical significance (<2 sigma) given the size of our sample (26 SN Ia). We highlight the importance of an accurate relative flux calibration and the failure of this method for highly-reddened objects. Comparison with synthetic spectra from 2D delayed-detonation explosion models shows that the correlation of R(6630/4400) with SN Ia absolute magnitudes can be largely attributed to intrinsic color variations and not to reddening by dust in the host galaxy. We consider flux ratios at other ages, as well as the use of pairs of flux ratios, revealing the presence of small-scale intrinsic spectroscopic variations in the iron-group dominated absorption features around ~4300 A and ~4800 A. The best flux ratio overall is the color-corrected R^c(4610/4260) at t=-2.5d from maximum light, which leads to ~30% lower scatter with respect to the standard combination of light-curve width and color. We examine other spectroscopic indicators related to line-profile morphology, but none appear to lead to a significant improvement over the standard light-curve width and color parameters.

From Sloan Digital Sky Survey u'g'r'i'z' imaging, we estimate the stellar masses of the host galaxies of 70 low redshift SN Ia (0.015 < z < 0.08) from the hosts' absolute luminosities and mass-to-light ratios. These nearby SN were discovered largely by searches targeting luminous galaxies, and we find that their host galaxies are substantially more massive than the hosts of SN discovered by the flux-limited Supernova Legacy Survey. Testing four separate light curve fitters, we detect ~2.5{\sigma} correlations of Hubble residuals with both host galaxy size and stellar mass, such that SN Ia occurring in physically larger, more massive hosts are ~10% brighter after light curve correction. The Hubble residual is the deviation of the inferred distance modulus to the SN, calculated from its apparent luminosity and light curve properties, away from the expected value at the SN redshift. Marginalizing over linear trends in Hubble residuals with light curve parameters shows that the correlations cannot be attributed to a light curve-dependent calibration error. Combining 180 higher-redshift ESSENCE, SNLS, and HigherZ SN with 30 nearby SN whose host masses are less than 10^10.8 solar masses in a cosmology fit yields 1+w=0.22 +0.152/-0.143, while a combination where the 30 nearby SN instead have host masses greater than 10^10.8 solar masses yields 1+w=-0.03 +0.217/-0.108. Progenitor metallicity, stellar population age, and dust extinction correlate with galaxy mass and may be responsible for these systematic effects. Host galaxy measurements will yield improved distances to SN Ia.

We present a comprehensive statistical analysis of the properties of Type Ia SN light curves in the near infrared using recent data from PAIRITEL and the literature. We construct a hierarchical Bayesian framework, incorporating several uncertainties including photometric error, peculiar velocities, dust extinction and intrinsic variations, for coherent statistical inference. SN Ia light curve inferences are drawn from the global posterior probability of parameters describing both individual supernovae and the population conditioned on the entire SN Ia NIR dataset. The logical structure of the hierarchical model is represented by a directed acyclic graph. Fully Bayesian analysis of the model and data is enabled by an efficient MCMC algorithm exploiting the conditional structure using Gibbs sampling. We apply this framework to the JHK_s SN Ia light curve data. A new light curve model captures the observed J-band light curve shape variations. The intrinsic variances in peak absolute magnitudes are: sigma(M_J) = 0.17 +/- 0.03, sigma(M_H) = 0.11 +/- 0.03, and sigma(M_Ks) = 0.19 +/- 0.04. We describe the first quantitative evidence for correlations between the NIR absolute magnitudes and J-band light curve shapes, and demonstrate their utility for distance estimation. The average residual in the Hubble diagram for the training set SN at cz > 2000 km/s is 0.10 mag. The new application of bootstrap cross-validation to SN Ia light curve inference tests the sensitivity of the model fit to the finite sample and estimates the prediction error at 0.15 mag. These results demonstrate that SN Ia NIR light curves are as effective as optical light curves, and, because they are less vulnerable to dust absorption, they have great potential as precise and accurate cosmological distance indicators.

We have obtained 1087 NIR (JHKs) measurements of 21 SNe Ia using PAIRITEL, nearly doubling the number of well-sampled NIR SN Ia light curves. These data strengthen the evidence that SNe Ia are excellent standard candles in the NIR, even without correction for optical light-curve shape. We construct fiducial NIR templates for normal SNe Ia from our sample, excluding only the three known peculiar SNe Ia: SN 2005bl, SN 2005hk, and SN 2005ke. The H-band absolute magnitudes in this sample of 18 SNe Ia have an intrinsic rms of only 0.15 mag with no correction for light-curve shape. We found a relationship between the H-band extinction and optical color excess of AH=0.2E(B-V). This variation is as small as the scatter in distance modulus measurements currently used for cosmology based on optical light curves after corrections for light-curve shape. Combining the homogeneous PAIRITEL measurements with 23 SNe Ia from the literature, these 41 SNe Ia have standard H-band magnitudes with an rms scatter of 0.16 mag. The good match of our sample with the literature sample suggests there are few systematic problems with the photometry. We present a nearby NIR Hubble diagram that shows no correlation of the residuals from the Hubble line with light-curve properties. Future samples that account for optical and NIR light-curve shapes, absorption, spectroscopic variation, or host-galaxy properties may reveal effective ways to improve the use of SNe Ia as distance indicators. Since systematic errors due to dust absorption in optical bands remain the leading difficulty in the cosmological use of supernovae, the good behavior of SN Ia NIR light curves and their relative insensitivity to reddening make these objects attractive candidates for future cosmological work.

We describe the effects of weak gravitational lensing by cosmological large scale structure on the diffuse emission of 21 centimeter radiation from neutral hydrogen at high redshifts during the era of reionization. The ability to observe radial information through the frequency, and thus three-dimensional regions of the background radiation at different redshifts, suggests that 21 cm studies may provide a useful context for studying weak lensing effects. We focus on the gravitational lensing effects on both the angular power spectra and the intrinsic, three-dimensional power spectra. We present a new approach for calculating the weak lensing signature based on integrating differential Fourier-space shells of the deflection field and approximating the magnification matrix. This method is applied to reionization models of the 21 cm spectra up to small angular scales over a range in redshift. The effect on the angular power spectrum is typically < 1% on small angular scales, and very small on scales corresponding to the feature imprinted by reionization bubbles, due to the near-scale invariance of the angular power spectrum of the 21 cm signal on these scales. We describe the expected effect of weak lensing on three-dimensional 21 cm power spectra, and show that lensing creates aspherical perturbations to the intrinsic power spectrum which depend on the polar angle of the wavevector. The effect on the 3D power spectrum is < 1% on scales k < 0.1 h/Mpc, but can be > 1% for highly inclined modes for k > 1 h/Mpc. The angular variation of the lensing effect on these scales is well described by a quartic polynomial in the cosine of the polar angle. The detection of the gravitational lensing effects on 21 cm power spectra will require very sensitive, high resolution observations by future low-frequency radio arrays.

Our fundamental lack of understanding of the acceleration of the Universe suggests that we consider a kinematic description. The simplest formalism involves the third derivative of the scale factor through a jerk parameter. A new approach is presented for describing the results of astronomical observations in terms of the contemporary jerk parameter and this is related to the equation of state approach. Simple perturbative expansions about $\Lambda$CDM are given.

We present exact analytic formulae for the eclipse of a star described by quadratic or nonlinear limb darkening. In the limit that the planet radius is less than a tenth of the stellar radius, we show that the exact lightcurve can be well approximated by assuming the region of the star blocked by the planet has constant surface brightness. We apply these results to the HST observations of HD 209458, showing that the ratio of the planetary to stellar radii is 0.1207+-0.0003. These formulae give a fast and accurate means of computing lightcurves using limb-darkening coefficients from model atmospheres which should aid in the detection, simulation, and parameter fitting of planetary transits.