Massive black hole (MBH) binary inspiral time scales are uncertain, and their spins are even more poorly constrained. Spin misalignment, along with unequal mass ratios and spin magnitudes, introduces asymmetry in the gravitational radiation, which imparts a recoil kick to the merged MBH. Understanding how MBH binary spins evolve is crucial for determining their recoil velocities, their gravitational wave (GW) waveforms detectable with LISA, as well as their post-merger retention rate in galaxies and thus their subsequent merger rate. Here we present a novel study that introduces a sub-resolution model for gas- and GW-driven MBH binary spin evolution using a population of accreting MBHs from the Illustris cosmological hydrodynamics simulations. We also model sub-resolution binary inspiral via dynamical friction, stellar scattering, viscous gas drag, and GW emission. Our model assumes differential accretion, which causes greater alignment of the secondary MBH spin in unequal-mass mergers. We find that 47% of the MBHs in our population merge by $z=0$. Of these, 19% have misaligned primaries and 10% have misaligned secondaries at the time of merger in our (conservative) fiducial model. The MBH misalignment fraction depends strongly on the accretion disc parameters, however. Reducing accretion rates by a factor of 100, in a thicker disc, yields 79% and 42% misalignment for primaries and secondaries, respectively. Even in the fiducial model, more than 12% of binaries experience recoils of $>500$ km s$^{-1}$, which could displace them at least temporarily from galactic nuclei. We additionally find that a significant number of systems experience strong precession.

We present a new analysis of the profile data from the 47 millisecond pulsars comprising the 12.5-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), which is presented in a parallel paper (Alam et al.; NG12.5). Our reprocessing is performed using "wideband" timing methods, which use frequency-dependent template profiles, simultaneous time-of-arrival (TOA) and dispersion measure (DM) measurements from broadband observations, and novel analysis techniques. In particular, the wideband DM measurements are used to constrain the DM portion of the timing model. We compare the ensemble timing results to NG12.5 by examining the timing residuals, timing models, and noise model components. There is a remarkable level of agreement across all metrics considered. Our best-timed pulsars produce encouragingly similar results to those from NG12.5. In certain cases, such as high-DM pulsars with profile broadening, or sources that are weak and scintillating, wideband timing techniques prove to be beneficial, leading to more precise timing model parameters by 10-15%. The high-precision multi-band measurements in several pulsars indicate frequency-dependent DMs. The TOA volume is reduced by a factor of 33, which may ultimately facilitate computational speed-ups for complex pulsar timing array analyses. This first wideband pulsar timing data set is a stepping stone, and its consistent results with NG12.5 assure us that such data sets are appropriate for gravitational wave analyses.

We present time-of-arrival measurements and timing models of 47 millisecond pulsars (MSPs) observed from 2004 to 2017 at the Arecibo Observatory and the Green Bank Telescope by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). The observing cadence was three to four weeks for most pulsars over most of this time span, with weekly observations of six sources. These data were collected for use in low-frequency gravitational wave searches and for other astrophysical purposes. We detail our observational methods and present a set of time-of-arrival (TOA) measurements, based on "narrowband" analysis, in which many TOAs are calculated within narrow radio-frequency bands for data collected simultaneously across a wide bandwidth. A separate set of "wideband" TOAs will be presented in a companion paper. We detail a number of methodological changes compared to our previous work which yield a cleaner and more uniformly processed data set. Our timing models include several new astrometric and binary pulsar measurements, including previously unpublished values for the parallaxes of PSRs J1832-0836 and J2322+2057, the secular derivatives of the projected semi-major orbital axes of PSRs J0613-0200 and J2229+2643, and the first detection of the Shapiro delay in PSR J2145-0750. We report detectable levels of red noise in the time series for fourteen pulsars. As a check on timing model reliability, we investigate the stability of astrometric parameters across data sets of different lengths. We report flux density measurements for all pulsars observed. Searches for stochastic and continuous gravitational waves using these data will be subjects of forthcoming publications.

Massive black hole binaries (MBHBs) form as a consequence of galaxy mergers. However, it is still unclear whether they typically merge within a Hubble time, and how accretion may affect their evolution. These questions will be addressed by pulsar timing arrays (PTAs), which aim to detect the GW background (GWB) emitted by MBHBs during the last Myrs of inspiral. Here we investigate the influence of differential accretion on MBHB merger rates, chirp masses and the resulting GWB spectrum. We evolve a MBHB sample from the Illustris hydrodynamic cosmological simulation using semi-analytic models and for the first time self-consistently evolve their masses with binary accretion models. In all models, MBHBs coalesce with median total masses up to $1.5 \times 10^8 M_{\odot}$, up to $3-4$ times larger than in models neglecting accretion. In our model with the largest plausible impact, the median mass ratio of coalescing MBHBs increases by a factor $3.6$, the coalescence rate by $52.3\%$, and the GWB amplitude by a factor $4.0$, yielding a dimensionless GWB strain $A_{yr^{-1}} = 1 \times 10^{-15}$. Our model that favours accretion onto the primary MBH reduces the median mass ratio of coalescing MBHBs by a factor of $2.9$, and yields a GWB amplitude $A_{yr^{-1}} = 3.1 \times 10^{-16}$. This is nearly indistinguishable from our model neglecting accretion, despite higher MBHB masses at coalescence. \textbf{We further predict binary separation and mass ratio distributions of stalled MBHBs in the low-redshift universe, and find that these depend sensitively on binary accretion models. This presents the potential for combined EM and GW observational constraints on merger rates and accretion models of MBHB populations.}

The gravitational waves from Massive black-hole (MBH) binaries are expected to be detected by pulsar timing arrays in the next few years. While they are a promising source for multimessenger observations as binary AGN, few convincing candidates have been identified in electromagnetic surveys. One approach to identifying candidates has been through spectroscopic surveys searching for offsets or time-dependent offsets of broad emission lines (BLs), which may be characteristic of Doppler shifts from binary orbital motion. In this study, we predict the parameter space of MBH binaries that should be kinematically detectable. There is a delicate trade-off between requiring binary separations to be large enough for BL regions to remain attached to one of the AGN, but also small enough such that their orbital velocity is detectable. We find that kinematic signatures are only observable for the lower-mass secondary AGN, for binaries with total-masses above about $10^8 M_\odot$, and separations between 0.1 and 1 pc. We motivate our usage of a kinematic-offset sensitivity of $10^3$ km/s, and a sensitivity to changing offsets of $10^2$ km/s. With these parameters, and an Eddington ratio of 0.1, we find that 0.5% of binaries have detectable offsets, and only 0.03% have detectable velocity changes. Overall, kinematic binary signatures should be expected in fewer than one in $10^4$ AGN. Better characterizing the intrinsic variability of BLs is crucial to understanding and vetting MBH binary candidates. This requires multi-epoch spectroscopy of large populations of AGN over a variety of timescales.

We present a detailed analysis of two well-localized, highly-offset short gamma-ray bursts---GRB 070809 and GRB 090515---investigating the kinematic evolution of their progenitors from compact object formation until merger. Calibrating to observations of their most probable host galaxies, we construct semi-analytic galactic models that account for star formation history and galaxy growth over time. We pair detailed kinematic evolution with compact binary population modeling to infer viable post-supernova velocities and inspiral times. By populating binary tracers according to the star formation history of the host and kinematically evolving their post-supernova trajectories through the time-dependent galactic potential, we find that systems matching the observed offsets of the bursts require post-supernova systemic velocities of $\gtrsim 250~\mathrm{km\,s}^{-1}$. Marginalizing over uncertainties in the stellar mass-halo mass relation, we find that the second-born neutron star in the GRB 070809 and GRB 090515 progenitor systems received a natal kick of $\gtrsim 200~\mathrm{km\,s}^{-1}$ at the 77% and 95% credible level, respectively. Applying our analysis to the full catalog of localized short gamma-ray bursts will provide unique constraints on their progenitors and help unravel the selection effects inherent to observing transients that are highly offset with respect to their hosts.

When galaxies merge, the supermassive black holes in their centers may form binaries and, during the process of merger, emit low-frequency gravitational radiation in the process. In this paper we consider the galaxy 3C66B, which was used as the target of the first multi-messenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C66B to less than $(1.65\pm0.02) \times 10^9~{M_\odot}$ using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits, and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data to `blind' pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

The regularity of pulsar emissions becomes apparent once we reference the pulses' times of arrivals to the inertial rest frame of the solar system. It follows that errors in the determination of Earth's position with respect to the solar-system barycenter can appear as a time-correlated bias in pulsar-timing residual time series, affecting the searches for low-frequency gravitational waves performed with pulsar timing arrays. Indeed, recent array datasets yield different gravitational-wave background upper limits and detection statistics when analyzed with different solar-system ephemerides. Crucially, the ephemerides do not generally provide usable error representations. In this article we describe the motivation, construction, and application of a physical model of solar-system ephemeris uncertainties, which focuses on the degrees of freedom (Jupiter's orbital elements) most relevant to gravitational-wave searches with pulsar timing arrays. This model, BayesEphem, was used to derive ephemeris-robust results in NANOGrav's 11-yr stochastic-background search, and it provides a foundation for future searches by NANOGrav and other consortia. The analysis and simulations reported here suggest that ephemeris modeling reduces the gravitational-wave sensitivity of the 11-yr dataset; and that this degeneracy will vanish with improved ephemerides and with the longer pulsar timing datasets that will become available in the near future.

The mergers of supermassive black hole binaries (SMBHBs) promise to be incredible sources of gravitational waves (GWs). While the oscillatory part of the merger gravitational waveform will be outside the frequency sensitivity range of pulsar timing arrays (PTAs), the non-oscillatory GW memory effect is detectable. Further, any burst of gravitational waves will produce GW memory, making memory a useful probe of unmodeled exotic sources and new physics. We searched the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GW memory. This dataset is sensitive to very low frequency GWs of $\sim3$ to $400$ nHz (periods of $\sim11$ yr $-$ $1$ mon). Finding no evidence for GWs, we placed limits on the strain amplitude of GW memory events during the observation period. We then used the strain upper limits to place limits on the rate of GW memory causing events. At a strain of $2.5\times10^{-14}$, corresponding to the median upper limit as a function of source sky position, we set a limit on the rate of GW memory events at $<0.4$ yr$^{-1}$. That strain corresponds to a SMBHB merger with reduced mass of $\eta M \sim 2\times10^{10}M_\odot$ and inclination of $\iota=\pi/3$ at a distance of 1 Gpc. As a test of our analysis, we analyzed the NANOGrav 9-year data set as well. This analysis found an anomolous signal, which does not appear in the 11-year data set. This signal is not a GW, and its origin remains unknown.

ESA and NASA are moving forward with plans to launch LISA around 2034. With data from the Illustris cosmological simulation, we provide analysis of LISA detection rates accompanied by characterization of the merging massive black hole population. Massive black holes of total mass $\sim10^5-10^{10} M_\odot$ are the focus of this study. We evolve Illustris massive black hole mergers, which form at separations on the order of the simulation resolution ($\sim$kpc scales), through coalescence with two different treatments for the binary massive black hole evolutionary process. The coalescence times of the population, as well as physical properties of the black holes, form a statistical basis for each evolutionary treatment. From these bases, we Monte Carlo synthesize many realizations of the merging massive black hole population to build mock LISA detection catalogs. We analyze how our massive black hole binary evolutionary models affect detection rates and the associated parameter distributions measured by LISA. With our models, we find massive black hole binary detection rates with LISA of $\sim0.5-1$ yr$^{-1}$ for massive black holes with masses greater than $10^5M_\odot$. This should be treated as a lower limit primarily because our massive black hole sample does not include masses below $10^5M_\odot$, which may significantly add to the observed rate. We suggest reasons why we predict lower detection rates compared to much of the literature.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a spectral index of -2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new "time slicing" analysis of the 11-year data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential "false positive" signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and show how independently modeling noise across 30 Fourier components, spanning NANOGrav's frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red noise process with low-spectral index noise, and use a tailored noise model to remove these effects.

Observations indicate that nearly all galaxies contain supermassive black holes (SMBHs) at their centers. When galaxies merge, their component black holes form SMBH binaries (SMBHBs), which emit low-frequency gravitational waves (GWs) that can be detected by pulsar timing arrays (PTAs). We have searched the recently-released North American Nanohertz Observatory for Gravitational Waves (NANOGrav) 11-year data set for GWs from individual SMBHBs in circular orbits. As we did not find strong evidence for GWs in our data, we placed 95\% upper limits on the strength of GWs from such sources as a function of GW frequency and sky location. We placed a sky-averaged upper limit on the GW strain of $h_0 < 7.3(3) \times 10^{-15}$ at $f_\mathrm{gw}= 8$ nHz. We also developed a technique to determine the significance of a particular signal in each pulsar using ``dropout' parameters as a way of identifying spurious signals in measurements from individual pulsars. We used our upper limits on the GW strain to place lower limits on the distances to individual SMBHBs. At the most-sensitive sky location, we ruled out SMBHBs emitting GWs with $f_\mathrm{gw}= 8$ nHz within 120 Mpc for $\mathcal{M} = 10^9 \, M_\odot$, and within 5.5 Gpc for $\mathcal{M} = 10^{10} \, M_\odot$. We also determined that there are no SMBHBs with $\mathcal{M} > 1.6 \times 10^9 \, M_\odot$ emitting GWs in the Virgo Cluster. Finally, we estimated the number of potentially detectable sources given our current strain upper limits based on galaxies in Two Micron All-Sky Survey (2MASS) and merger rates from the Illustris cosmological simulation project. Only 34 out of 75,000 realizations of the local Universe contained a detectable source, from which we concluded it was unsurprising that we did not detect any individual sources given our current sensitivity to GWs.

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.

We present the full public release of all data from the TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive blackholes from early time to the present day, z=0. Each of the flagship runs -- TNG50, TNG100, and TNG300 -- are accompanied by lower-resolution and dark-matter only counterparts, and we discuss scientific and numerical cautions and caveats relevant when using TNG. Full volume snapshots are available at 100 redshifts; halo and subhalo catalogs at each snapshot and merger trees are also released. The data volume now directly accessible online is ~750 TB, including 1200 full volume snapshots and ~80,000 high time-resolution subbox snapshots. This will increase to ~1.1 PB with the future release of TNG50. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform which supports user computation with fully local access to TNG data, alleviating the need to download large simulated datasets.

Precision timing of large arrays (>50) of millisecond pulsars will detect the nanohertz gravitational-wave emission from supermassive binary black holes within the next ~3-7 years. We review the scientific opportunities of these detections, the requirements for success, and the synergies with electromagnetic instruments operating in the 2020s.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Massive black-hole (MBH) binaries, which are expected to form following the merger of their parent galaxies, produce gravitational waves which will be detectable by Pulsar Timing Arrays at nanohertz frequencies (year periods). While no confirmed, compact MBH binary systems have been seen in electromagnetic observations, a large number of candidates have recently been identified in optical surveys of AGN variability. Using a combination of cosmological, hydrodynamic simulations; comprehensive, semi-analytic binary merger models; and analytic AGN spectra and variability prescriptions; we calculate the expected electromagnetic detection rates of MBH binaries as periodically variable AGN. In particular, we consider two independent variability models: (i) Doppler boosting due to large orbital velocities, and (ii) hydrodynamic variability in which the fueling of MBH accretion disks is periodically modulated by the companion. Our models predict that numerous MBH binaries should be present and distinguishable in the existing data. In particular, our fiducial models produce an expectation value of $0.2$ (Doppler) and $5$ (hydrodynamic) binaries to be identifiable in CRTS, while $20$ and $100$ are expected after five years of LSST observations. The brightness variations in most systems are too small to be distinguishable, but almost $1\%$ of AGN at redshifts $z \lesssim 0.6$ could be in massive binaries. We analyze the predicted binary parameters of observable systems and their selection biases, and include an extensive discussion of our model parameters and uncertainties.

The advent of time domain astronomy is revolutionizing our understanding of the Universe. Programs such as the Catalina Real-time Transient Survey (CRTS) or the Palomar Transient Factory (PTF) surveyed millions of objects for several years, allowing variability studies on large statistical samples. The inspection of $\approx$250k quasars in CRTS resulted in a catalogue of 111 potentially periodic sources, put forward as supermassive black hole binary (SMBHB) candidates. A similar investigation on PTF data yielded 33 candidates from a sample of $\approx$35k quasars. Working under the SMBHB hypothesis, we compute the implied SMBHB merger rate and we use it to construct the expected gravitational wave background (GWB) at nano-Hz frequencies, probed by pulsar timing arrays (PTAs). After correcting for incompleteness and assuming virial mass estimates, we find that the GWB implied by the CRTS sample exceeds the current most stringent PTA upper limits by almost an order of magnitude. After further correcting for the implicit bias in virial mass measurements, the implied GWB drops significantly but is still in tension with the most stringent PTA upper limits. Similar results hold for the PTF sample. Bayesian model selection shows that the null hypothesis (whereby the candidates are false positives) is preferred over the binary hypothesis at about $2.3\sigma$ and $3.6\sigma$ for the CRTS and PTF samples respectively. Although not decisive, our analysis highlights the potential of PTAs as astrophysical probes of individual SMBHB candidates and indicates that the CRTS and PTF samples are likely contaminated by several false positives.

We calculate the properties, occurrence rates and detection prospects of individually resolvable 'single sources' in the low frequency gravitational wave (GW) spectrum. Our simulations use the population of galaxies and massive black hole binaries from the Illustris cosmological hydrodynamic simulations, coupled to comprehensive semi-analytic models of the binary merger process. Using mock pulsar timing arrays (PTA) with, for the first time, varying red-noise models, we calculate plausible detection prospects for GW single sources and the stochastic GW background (GWB). Contrary to previous results, we find that single sources are at least as detectable as the GW background. Using mock PTA, we find that these 'foreground' sources (also 'deterministic'/'continuous') are likely to be detected with $\sim 20 \, \textrm{yr}$ total observing baselines. Detection prospects, and indeed the overall properties of single sources, are only moderately sensitive to binary evolution parameters---namely eccentricity & environmental coupling, which can lead to differences of $\sim 5 \,\textrm{yr}$ in times to detection. Red noise has a stronger effect, roughly doubling the time to detection of the foreground between a white-noise only model ($\sim 10$ - $15 \, \textrm{yr}$) and severe red noise ($\sim 20$ - $30 \, \textrm{yr}$). The effect of red noise on the GWB is even stronger, suggesting that single source detections may be more robust. We find that typical signal-to-noise ratios for the foreground peak near $f = 0.1 \, \textrm{yr}^{-1}$, and are much less sensitive to the continued addition of new pulsars to PTA.

Pulsar Timing Arrays (PTA) around the world are using the incredible consistency of millisecond pulsars to measure low frequency gravitational waves from (super)Massive Black Hole (MBH) binaries. We use comprehensive MBH merger models based on cosmological hydrodynamic simulations to predict the spectrum of the stochastic Gravitational-Wave Background (GWB). We use real Time-of-Arrival (TOA) specifications from the European, NANOGrav, Parkes, and International PTA (IPTA) to calculate realistic times to detection of the GWB across a wide range of model parameters. In addition to exploring the parameter space of environmental hardening processes (in particular: stellar scattering efficiencies), we have expanded our models to include eccentric binary evolution which can have a strong effect on the GWB spectrum. Our models show that strong stellar scattering and high characteristic eccentricities enhance the GWB strain amplitude near the PTA sensitive "sweet-spot" (near the frequency $f = 1 \, \mathrm{yr}^{-1}$), slightly improving detection prospects in these cases. While the GWB $amplitude$ is degenerate between cosmological and environmental parameters, the location of a spectral turnover at low frequencies ($f \lesssim 0.1 \, \mathrm{yr}^{-1}$) is strongly indicative of environmental coupling. At high frequencies ($f\gtrsim 1 \, \mathrm{yr}^{-1}$), the GWB spectral index can be used to infer the number density of sources and possibly their eccentricity distribution. Even with merger models that use pessimistic environmental and eccentricity parameters, if the current rate of PTA expansion continues, we find that the International PTA is highly likely to make a detection within about 10 years.

We present the Open Supernova Catalog, an online collection of observations and metadata for presently 36,000+ supernovae and related candidates. The catalog is freely available on the web (https://sne.space), with its main interface having been designed to be a user-friendly, rapidly-searchable table accessible on desktop and mobile devices. In addition to the primary catalog table containing supernova metadata, an individual page is generated for each supernova which displays its available metadata, light curves, and spectra spanning X-ray to radio frequencies. The data presented in the catalog is automatically rebuilt on a daily basis and is constructed by parsing several dozen sources, including the data presented in the supernova literature and from secondary sources such as other web-based catalogs. Individual supernova data is stored in the hierarchical, human- and machine-readable JSON format, with the entirety of each supernova's data being contained within a single JSON file bearing its name. The setup we present here, which is based upon open source software maintained via git repositories hosted on github, enables anyone to download the entirety of the supernova dataset to their home computer in minutes, and to make contributions of their own data back to the catalog via git. As the supernova dataset continues to grow, especially in the upcoming era of all-sky synoptic telescopes which will increase the total number of events by orders of magnitude, we hope that the catalog we have designed will be a valuable tool for the community to analyze both historical and contemporary supernovae.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

Gravitational Waves (GW) have now been detected from stellar-mass black hole binaries, and the first observations of GW from Massive Black Hole (MBH) Binaries are expected within the next decade. Pulsar Timing Arrays (PTA), which can measure the years long periods of GW from MBHB, have excluded many standard predictions for the amplitude of a stochastic GW Background (GWB). We use coevolved populations of MBH and galaxies from hydrodynamic, cosmological simulations ('Illustris') to calculate a predicted GWB. The most advanced predictions so far have included binary hardening mechanisms from individual environmental processes. We present the first calculation including all of the environmental mechanisms expected to be involved: dynamical friction, stellar 'loss-cone' scattering, and viscous drag from a circumbinary disk. We find that MBH binary lifetimes are generally multiple gigayears, and only a fraction coalesce by redshift zero. For a variety of parameters, we find all GWB amplitudes to be below the most stringent PTA upper limit of $A_{\textrm{yr}^{-1}} \approx 10^{-15}$. Our fairly conservative fiducial model predicts an amplitude of $A_{\textrm{yr}^{-1}} \approx 0.4\times 10^{-15}$---less than a factor of three below the current limit. At lower frequencies, we find $A_{0.1\,\textrm{yr}^{-1}} \approx 1.5\times 10^{-15}$ with spectral indices between $-0.4$ and $-0.6$---significantly flatter than the canonical value of $-2/3$ due to purely GW-driven evolution. Typical MBHB driving the GWB signal come from redshifts around $0.3$, with total masses of a few times $10^9\,M_\odot$, and in host galaxies with very large stellar masses. Even without GWB detections, our results can be connected to observations of dual AGN to constrain binary evolution.

Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas-richness of host galaxies. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly-oriented spins, less than about 10 spatially-offset AGN may be detectable in HST-COSMOS, and > 10^3 could be found with Pan-STARRS, LSST, Euclid, and WFIRST. Nearly a thousand velocity-offset AGN are predicted within the SDSS footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. Nonetheless, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN.

The transient Swift J1644+57 is believed to have been produced by an unlucky star wandering too close to a supermassive black hole (BH) leading to a tidal disruption event. This unusual flare displayed highly super-Eddington X-ray emission which likely originated in a relativistic, collimated jet. This presents challenges to modern accretion and jet theory as upper limits of prior BH activity, which we obtain from the radio afterglow of this event, imply that both the pre-disruption BH and stellar magnetic fluxes fall many orders of magnitude short of what is required to power the observed X-ray luminosity. We argue that a pre-existing, "fossil" accretion disc can contain a sufficient reservoir of magnetic flux and that the stellar debris stream is capable of dragging this flux into the BH. To demonstrate this, we perform local, 3D magnetohydrodynamic simulations of the disc--stream interaction and demonstrate that the interface between the two is unstable to mixing. This mixing entrains a sufficient amount of fossil disc magnetic flux into the infalling stellar debris to power the jet. We argue that the interaction with the fossil disc can have a pronounced effect on the structure and dynamics of mass fallback and likely the resulting transient. Finally, we describe possible ramifications of these interactions on unresolved problems in tidal disruption dynamics, in particular, the efficiency of debris circularization, and effects of the disruption on the preexisting black hole system. Animations online: http://goo.gl/T84tLs

The detection of an electromagnetic transient which may originate from a binary neutron star merger can increase the probability that a given segment of data from the LIGO-Virgo ground-based gravitational-wave detector network contains a signal from a binary coalescence. Additional information contained in the electromagnetic signal, such as the sky location or distance to the source, can help rule out false alarms, and thus lower the necessary threshold for a detection. Here, we develop a framework for determining how much sensitivity is added to a gravitational-wave search by triggering on an electromagnetic transient. We apply this framework to a variety of relevant electromagnetic transients, from short GRBs to signatures of r-process heating to optical and radio orphan afterglows. We compute the expected rates of multi-messenger observations in the Advanced detector era, and find that searches triggered on short GRBs --- with current high-energy instruments, such as Fermi --- and nucleosynthetic `kilonovae' --- with future optical surveys, like LSST --- can boost the number of multi-messenger detections by 15% and 40%, respectively, for a binary neutron star progenitor model. Short GRB triggers offer precise merger timing, but suffer from detection rates decreased by beaming and the high a priori probability that the source is outside the LIGO-Virgo sensitive volume. Isotropic kilonovae, on the other hand, could be commonly observed within the LIGO-Virgo sensitive volume with an instrument roughly an order of magnitude more sensitive than current optical surveys. We propose that the most productive strategy for making multi-messenger gravitational-wave observations is using triggers from future deep, optical all-sky surveys, with characteristics comparable to LSST, which could make as many as ten such coincident observations a year.

The unusual transient Swift J1644+57 likely resulted from a collimated relativistic jet powered by accretion onto a massive black hole (BH) following the tidal disruption (TD) of a star. Several mysteries cloud the interpretation of this event: (1) extreme flaring and `plateau' shape of the X-ray/gamma-ray light curve during the first 10 days after the gamma-ray trigger; (2) unexpected rebrightening of the forward shock radio emission months after trigger; (3) no obvious evidence for jet precession, despite misalignment typically expected between the angular momentum of the accretion disk and BH; (4) recent abrupt shut-off in jet X-ray emission after 1.5 years. Here we show that all of these seemingly disparate mysteries are naturally resolved by one assumption: the presence of strong magnetic flux Phi threading the BH. Initially, Phi is weak relative to high fall-back mass accretion rate, Mdot, and the disk and jets precess about the BH axis = our line of sight. As Mdot drops, Phi becomes dynamically important and leads to a magnetically-arrested disk (MAD). MAD naturally aligns disk and jet axis along the BH spin axis, but only after a violent rearrangement phase (jet wobbling). This explains the erratic light curve at early times and the lack of precession at later times. We use our model for Swift J1644+57 to constrain BH and disrupted star properties, finding that a solar-mass main sequence star disrupted by a relatively low mass, M~10^5-10^6 Msun, BH is consistent with the data, while a WD disruption (though still possible) is disfavored. The magnetic flux required to power Swift J1644+57 is too large to be supplied by the star itself, but it could be collected from a quiescent `fossil' accretion disk present in the galactic nucleus prior to the TD. The presence (lack of) of such a fossil disk could be a deciding factor in what TD events are accompanied by powerful jets.[abridged]

We discuss two approaches to searches for gravitational-wave (GW) and electromagnetic (EM) counterparts of binary neutron star mergers. The first approach relies on triggering archival searches of GW detector data based on detections of EM transients. We introduce a quantitative approach to evaluate the improvement to GW detector reach due to the extra information gained from the EM transient and the increased confidence in the presence of a signal from a binary merger. We also advocate utilizing other transients in addition to short gamma ray bursts. The second approach involves following up GW candidates with targeted EM observations. We argue for the use of slower but optimal parameter-estimation techniques to localize the source on the sky, and for a more sophisticated use of astrophysical prior information, including galaxy catalogs, to find preferred followup locations.

Merging compact binaries are the most viable and best studied candidates for gravitational wave (GW) detection by the fully operational network of ground-based observatories. In anticipation of the first detections, the expected distribution of GW sources in the local universe is of considerable interest. Here we investigate the full phase space distribution of coalescing compact binaries at $z = 0$ using dark matter simulations of structure formation. The fact that these binary systems acquire large barycentric velocities at birth ("kicks") results in merger site distributions that are more diffusely distributed with respect to their putative hosts, with mergers occurring out to distances of a few Mpc from the host halo. Redshift estimates based solely on the nearest galaxy in projection can, as a result, be inaccurate. On the other hand, large offsets from the host galaxy could aid the detection of faint optical counterparts and should be considered when designing strategies for follow-up observations. The degree of isotropy in the projected sky distributions of GW sources is found to be augmented with increasing kick velocity and to be severely enhanced if progenitor systems possess large kicks as inferred from the known population of pulsars and double compact binaries. Even in the absence of observed electromagnetic counterparts, the differences in sky distributions of binaries produced by disparate kick-velocity models could be discerned by GW observatories, within the expected accuracies and detection rates of advanced LIGO--in particular with the addition of more interferometers.

The diffuse gamma radiation arising from the interaction of cosmic ray particles with matter and radiation in the Galaxy is one of the few probes available to study the origin of the cosmic rays. Milagro is a water Cherenkov detector that continuously views the entire overhead sky. The large field-of-view combined with the long observation time makes Milagro the most sensitive instrument available for the study of large, low surface brightness sources such as the diffuse gamma radiation arising from interactions of cosmic radiation with interstellar matter. In this paper we present spatial and flux measurements of TeV gamma-ray emission from the Cygnus Region. The TeV image shows at least one new source MGRO J2019+37 as well as correlations with the matter density in the region as would be expected from cosmic-ray proton interactions. However, the TeV gamma-ray flux as measured at ~12 TeV from the Cygnus region (after excluding MGRO J2019+37) exceeds that predicted from a conventional model of cosmic ray production and propagation. This observation indicates the existence of either hard-spectrum cosmic-ray sources and/or other sources of TeV gamma rays in the region.

The Milagro gamma-ray observatory employs a water Cherenkov detector to observe extensive air showers produced by high energy particles interacting in the Earth's atmosphere. Milagro has a wide field of view and high duty cycle, monitoring the northern sky almost continuously in the 100 GeV to 100 TeV energy range. Milagro is, thus, uniquely capable of searching for very high-energy emission from gamma-ray bursts (GRBs) during the prompt emission phase. Detection of >100 GeV counterparts would place powerful constraints on GRB mechanisms. Twenty-five satellite-triggered GRBs occurred within the field of view of Milagro between January 2000 and December 2001. We have searched for counterparts to these GRBs and found no significant emission from any of the burst positions. Due to the absorption of high-energy gamma rays by the extragalactic background light, detections are only expected to be possible for redshifts less than ~0.5. Three of the GRBs studied have measured redshifts. GRB 010921 has a redshift low enough (0.45) to allow an upper limit on the fluence to place an observational constraint on potential GRB models.

The Milagro telescope monitors the northern sky for 100 GeV to 100 TeV transient emission through continuous very high energy wide-field observations. The large effective area and ~100 GeV energy threshold of Milagro allow it to detect very high energy (VHE) gamma-ray burst emission with much higher sensitivity than previous instruments and a fluence sensitivity at VHE energies comparable to that of dedicated gamma-ray burst satellites at keV to MeV energies. Even in the absence of a positive detection, VHE observations can place important constraints on gamma-ray burst (GRB) progenitor and emission models. We present limits on the VHE flux of 40 s -- 3 h duration transients nearby to earth, as well as sensitivity distributions which have been corrected for gamma-ray absorption by extragalactic background light and cosmological effects. The sensitivity distributions suggest that the typical intrinsic VHE fluence of GRBs is similar or weaker than the keV -- MeV emission, and we demonstrate how these sensitivity distributions may be used to place observational constraints on the absolute VHE luminosity of gamma-ray bursts for any GRB emission and progenitor model.

Evidence of TeV emission from GRB970417a has been previously reported using data from the Milagrito detector. Constraints on the TeV fluence and the energy spectrum are now derived using additional data from a scaler system that recorded the rate of signals from the Milagrito photomultipliers. This analysis shows that if emission from GRB970417a has been observed, it must contain photons with energies above 650 GeV. Some consequences of this observation are discussed.

Milagrito, a detector sensitive to very high energy gamma rays, monitored the northern sky from February 1997 through May 1998. With a large field of view and a high duty cycle, this instrument was well suited to perform a search for TeV gamma-ray bursts (GRBs). We report on a search made for TeV counterparts to GRBs observed by BATSE. BATSE detected 54 GRBs within the field of view of Milagrito during this period. An excess of events coincident in time and space with one of these bursts, GRB 970417a, was observed by Milagrito. The excess has a chance probability of $2.8 \times 10^{-5}$ of being a fluctuation of the background. The probability for observing an excess at least this large from any of the 54 bursts is $1.5 \times 10^{-3}$. No significant correlations were detected from the other bursts.

Milagrito, a detector sensitive to gamma-rays at TeV energies, monitored the northern sky during the period February 1997 through May 1998. With a large field of view and high duty cycle, this instrument was used to perform a search for TeV counterparts to gamma-ray bursts. Within the Milagrito field of view 54 gamma-ray bursts at keV energies were observed by the Burst And Transient Satellite Experiment (BATSE) aboard the Compton Gamma-Ray Observatory. This paper describes the results of a preliminary analysis to search for TeV emission correlated with BATSE detected bursts. Milagrito detected an excess of events coincident both spatially and temporally with GRB 970417a, with chance probability $2.8 \times 10^{-5}$ within the BATSE error radius. No other significant correlations were detected. Since 54 bursts were examined the chance probability of observing an excess with this significance in any of these bursts is $1.5 \times 10^{-3}$. The statistical aspects and physical implications of this result are discussed.

Previous works have shown that water Cherenkov detectors have superior sensitivity to those of scintillation counters as applied to detecting extensive air showers (EAS). This is in large part due to their much higher sensitivity to EAS photons which are more than five times more numerous than EAS electrons. Large area water Cherenkov detectors can be constructed relatively cheaply and operated reliably. A sparse detector array has been designed which uses these types of detectors to substantially increase the area over which the Milagro Gamma Ray Observatory collects EAS information. Improvements to the Milagro detector's performance characteristics and sensitivity derived from this array and preliminary results from a prototype array currently installed near the Milagro detector will be presented.

The Milagro detector is an air shower array which uses the water Cherenkov technique and is capable of continuously monitoring the sky at energies near 1 TeV. The detector consists of 20000 metric tons of pure water instrumented with 723 photo-multiplier tubes (PMTs). The PMTs are arranged in a two-layer structure on a lattice of 3 m spacing covering 5000 $m^2$ area. The direction of the shower is determined from the relative timing of the PMT signals, necessitating a common time reference and amplitude slewing corrections to improve the time resolution. The calibration system to provide these consists of a pulsed laser driving 30 diffusing light sources deployed in the pond to allow cross-calibration of the PMTs. The system is capable of calibrating times and the pulse-heights from the PMTs using the time-over-threshold technique. The absolute energy scale is provided using single muons passing through the detector. The description of the calibration system of the Milagro detector and its prototype Milagrito will be presented.

Observing gamma ray bursts (GRBs) in the TeV energy range can be extremely valuable in providing insight to GRB radiation mechanisms and in constraining source distances. The Milagrito detector was an air shower array which used the water Cherenkov technique to search for TeV sources. Data from this detector was analyzed to look for a TeV component of GRBs coincident with low energy -rays detected by the BATSE instrument on the Compton Gamma Ray Observatory. A sample of 54 BATSE GRBs which were in the field of view of the Milagrito detector during its lifetime (February 1997 to May 1998) was used.

The Milagro Gamma Ray Observatory is the world's first large-area water Cherenkov detector capable of continuously monitoring the sky at TeV energies. Located in northern New Mexico, Milagro will perform an all sky survey of the Northern Hemisphere at energies between ~250 GeV and 50 TeV. With a high duty cycle, large detector area (~5000 square meters), and a wide field-of-view (~1 sr), Milagro is uniquely capable of searching for transient and DC sources of high-energy gamma-ray emission. Milagro has been operating since February, 1999. The current status of the Milagro Observatory and initial results will be discussed.

The Milagro water Cherenkov detector near Los Alamos, New Mexico, has been operated as a sky monitor at energies of a few TeV between February 1997 and April 1998. Serving as a test run for the full Milagro detector, Milagrito has taken data during the strong and long-lasting 1997 flare of Mrk 501. We present results from the analysis of Mrk 501 and compare the excess and background rates with expectations from the detector simulations.

Milagrito, a prototype for the Milagro detector, operated for 15 months in 1997-8 and collected 8.9 billion events. It was the first extensive air shower (EAS) array sensitive to showers intiated by primaries with energy below 1 TeV. The shadows of the sun and moon observed with cosmic rays can be used to study systematic pointing shifts and measure the angular resolution of EAS arrays. Below a few TeV, the paths of cosmic rays coming toward the earth are bent by the helio- and geo-magnetic fields. This is expected to distort and displace the shadows of the sun and the moon. The moon shadow, offset from the nominal (undeflected) position, has been observed with high statistical significance in Milagrito. This can be used to establish energy calibrations, as well as to search for the anti-matter content of the VHE cosmic ray flux. The shadow of the sun has also been observed with high significance.

The Milagrito water Cherenkov telescope operated for over a year. The most probable gamma-ray energy was ~1 TeV and the trigger rate was as high as 400 Hz. We have developed an efficient technique for searching the entire sky for short duration bursts of TeV photons. Such bursts may result from "traditional" gamma-ray bursts that were not in the field-of-view of any other instruments, the evaporation of primordial black holes, or some as yet undiscovered phenomenon. We have begun to search the Milagrito data set for bursts of duration 10 seconds. Here we will present the technique and the expected results. Final results will be presented at the conference.

Solar Energetic Particles from the 6 November 1997 solar flare/CME(coronal mass ejection) with energies exceeding 10 GeV have been detected by Milagrito, a prototype of the Milagro Gamma Ray Observatory. While particle acceleration beyond 1 GeV at the Sun is well established, few data exist for protons or ions beyond 10 GeV. The Milagro observatory, a ground based water Cherenkov detector designed for observing very high energy gamma ray sources, can also be used to study the Sun. Milagrito, which operated for approximately one year in 1997/98, was sensitive to solar proton and neutron fluxes above ~5- 10 GeV. Milagrito operated in a scaler mode, which was primarily sensitive to muons, low energy photons, and electrons, and the detector operated in a mode sensitive to showers and high zenith angle muons. In its scaler mode, Milagrito registered a rate increase coincident with the 6 November 1997 ground level event observed by Climax and other neutron monitors. A preliminary analysis suggests the presence of >10 GeV particles.

The Milagrito water Cherenkov detector in the Jemez Mountains near Los Alamos, New Mexico took data from February 1997 to April 1998. Milagrito served as a prototype for the larger Milagro detector, which has just begun operations. Milagrito was the first large-aperture gamma-ray detector with sensitivity to gamma rays below 1 TeV. We report here on a search for steady emission from point sources over most of the northern sky using data from Milagrito.