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.

We explore the relative role of small-scale fragmentation and global collapse in low-metallicity clouds, pointing out that in such clouds the cooling time may be longer than the dynamical time, allowing the cloud to collapse globally before it can fragment. This, we suggest, may help to explain the formation of the low-metallicity globular cluster population, since such dense stellar systems need a large amount of gas to be collected in a small region (without significant feedback during the collapse). To explore this further, we carry out numerical simulations of low-metallicity Bonner-Ebert stable gas clouds, demonstrating that there exists a critical metallicity (between 0.001 and 0.01 $Z_\odot$) below which the cloud collapses globally without fragmentation. We also run simulations including a background radiative heating source, showing that this can also produce clouds that do not fragment, and that the critical metallicity -- which can exceed the no-radiation case -- increases with the heating rate.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

We present Sedonu, a new open source, steady-state, special relativistic Monte Carlo (MC) neutrino transport code, available at bitbucket.org/srichers/sedonu. The code calculates the energy- and angle-dependent neutrino distribution function on fluid backgrounds of any number of spatial dimensions, calculates the rates of change of fluid internal energy and electron fraction, and solves for the equilibrium fluid temperature and electron fraction. We apply this method to snapshots from two-dimensional simulations of accretion disks left behind by binary neutron star mergers, varying the input physics and comparing to the results obtained with a leakage scheme for the case of a central black hole and a central hypermassive neutron star. Neutrinos are guided away from the densest regions of the disk and escape preferentially around 45 degrees from the equatorial plane. Neutrino heating is strengthened by MC transport a few scale heights above the disk midplane near the innermost stable circular orbit, potentially leading to a stronger neutrino-driven wind. Neutrino cooling in the dense midplane of the disk is stronger when using MC transport, leading to a globally higher cooling rate by a factor of a few and a larger leptonization rate by an order of magnitude. We calculate neutrino pair annihilation rates and estimate that an energy of 2.8e46 erg is deposited within 45 degrees of the symmetry axis over 300 ms when a central BH is present. Similarly, 1.9e48 erg is deposited over 3 s when an HMNS sits at the center, but neither estimate is likely to be sufficient to drive a GRB jet.

We develop analytic and numerical models of the properties of super-Eddington stellar winds, motivated by phases in stellar evolution when super-Eddington energy deposition (via, e.g., unstable fusion, wave heating, or a binary companion) heats a region near the stellar surface. This appears to occur in luminous blue variables (LBVs), Type IIn supernovae progenitors, classical novae, and X-ray bursts. We show that when the wind kinetic power exceeds Eddington, the photons are trapped and behave like a fluid. Convection does not play a significant role in the wind energy transport. The wind properties depend on the ratio of a characteristic speed in the problem vc ~ (Edot G)^{1/5} (where Edot is the heating rate) to the stellar escape speed near the heating region vesc(r_h). For vc > vesc(r_h) the wind kinetic power at large radii Edot_w ~ Edot. For vc < vesc(r_h), most of the energy is used to unbind the wind material and thus Edot_w < Edot. Multidimensional hydrodynamic simulations without radiation diffusion using FLASH and one-dimensional hydrodynamic simulations with radiation diffusion using MESA are in good agreement with the analytic predictions. The photon luminosity from the wind is itself super-Eddington but in many cases the photon luminosity is likely dominated by `internal shocks' in the wind. We discuss the application of our models to eruptive mass loss from massive stars and argue that the wind models described here can account for the broad properties of LBV outflows and the enhanced mass loss in the years prior to Type IIn core-collapse supernovae.

We study the radioactively-powered transients produced by accretion disk winds following a compact object merger. Starting with the outflows generated in two-dimensional hydrodynamical disk models, we use wavelength-dependent radiative transfer calculations to generate synthetic light curves and spectra. We show that the brightness and color of the resulting kilonova transients carry information about the merger physics. In the regions of the wind where neutrino irradiation raises the electron fraction to Ye > 0.25, r-process nucleosynthesis halts before producing high-opacity, complex ions (the lanthanides). The kilonova light curves thus show two distinct components: a brief (~2 day) blue optical transient produced in the outer lanthanide-free ejecta, and a longer (~10 day) infrared transient produced in the inner, lanthanide line-blanketed region. Mergers producing a longer-lived neutron star, or a more rapidly spinning black hole, have stronger neutrino irradiation, generate more lanthanide-free ejecta, and are optically brighter and bluer. At least some optical emission is produced in all disk wind models, which should enhance the detectability of electromagnetic counterparts to gravitational wave sources. However, the presence of even a small amount (10^{-4} Msun) of overlying, neutron-rich dynamical ejecta will act as a "lanthanide-curtain", obscuring the optical wind emission from certain viewing angles. Because the disk outflows have moderate velocities (~10,000 km/s), numerous resolved line features are discernible in the spectra, distinguishing disk winds from fast-moving dynamical ejecta, and offering a potential diagnostic of the detailed composition of freshly produced r-process material.

We resolve spatially the star formation history of 300 nearby galaxies from the CALIFA integral field survey to investigate: a) the radial structure and gradients of the present stellar populations properties as a function of the Hubble type; and b) the role that plays the galaxy stellar mass and stellar mass surface density in governing the star formation history and metallicity enrichment of spheroids and the disks of galaxies. We apply the fossil record method based on spectral synthesis techniques to recover spatially and temporally resolved maps of stellar population properties of spheroids and spirals with galaxy mass from 10$^9$ to 7$\times$10$^{11}$ M$_{\odot}$. The individual radial profiles of the stellar mass surface density ($\mu_{*}$), stellar extinction (A$_{V}$), luminosity weighted ages ($\langle$ log age $\rangle_{L}$), and mass weighted metallicity ($\langle$ log Z/Z$_{\odot}$$\rangle_{M}$) are stacked in seven bins of galaxy morphology (E, S0, Sa, Sb, Sbc, Sc and Sd). All these properties show negative gradients as a sight of the inside-out growth of massive galaxies. However, the gradients depend on the Hubble type in different ways. For the same galaxy mass, E and S0 galaxies show the largest inner gradients in $\mu_{*}$; and Andromeda-like galaxies (Sb with log M$_{*}$(M$_{\odot}$) $\sim$ 11) show the largest inner age and metallicity gradients. In average, spiral galaxies have a stellar metallicity gradient $\sim$ -0.1 dex per half-light radius, in agreement with the value estimated for the ionized gas oxygen abundance gradient by CALIFA. A global (M$_{*}$-driven) and local ($\mu_{*}$- driven) stellar metallicity relation are derived. We find that in disks, the stellar mass surface density regulates the stellar metallicity; in spheroids, the galaxy stellar mass dominates the physics of star formation and chemical enrichment.

The presence of quasars at redshifts z > 6 indicates the existence of supermassive black holes (SMBHs) as massive as a few times 10^9 Msun, challenging models for SMBH formation. One pathway is through the direct collapse of gas in T_{vir} > 10^4 K halos; however, this requires the suppression of H_2 cooling to prevent fragmentation. In this paper, we examine a proposed new mechanism for this suppression which relies on cold-mode accretion flows leading to shocks at high densities (n > 10^4 cm^{-3}) and temperatures (T > 10^4 K). In such gas, H_2 is efficiently collisionally dissociated. We use high-resolution numerical simulations to test this idea, demonstrating that such halos typically have lower temperature progenitors, in which cooling is efficient. Those halos do show filamentary flows; however, the gas shocks at or near the virial radius (at low densities), thus preventing the proposed collisional mechanism from operating. We do find that, if we artificially suppress H_2 formation with a high UV background, so as to allow gas in the halo center to enter the high-temperature, high-density "zone of no return", it will remain there even if the UV flux is turned off, collapsing to high density at high temperature. Due to computational limitations, we simulated only three halos. However, we demonstrate, using Monte Carlo calculations of 10^6 halo merger histories, that a few rare halos could assemble rapidly enough to avoid efficient H_2 cooling in all of their progenitor halos, provided that the UV background exceeds J_{21} ~ few at redshifts as high as z ~ 20.

The thermal emission detected from the millisecond pulsar J0437-4715 is not explained by standard cooling models of neutron stars without a heating mechanism. We investigated three heating mechanisms controlled by the rotational braking of the pulsar: breaking of the solid crust, superfluid vortex creep, and non-equilibrium reactions ('rotochemical heating'). We find that the crust cracking mechanism does not produce detectable heating. Given the dependence of the heating mechanisms on spin-down parameters, which leads to different temperatures for different pulsars, we study the thermal evolution for two types of pulsars: young, slowly rotating 'classical' pulsars and old, fast rotating millisecond pulsars (MSPs). We find that the rotochemical heating and vortex creep mechanism can be important both for classical pulsars and MSPs.

Although the surface of a magnetar is a source of bright thermal X-rays, its spectrum contains substantial non-thermal components. The X-ray emission is pulsed, with pulsed fractions that can be as high as ~ 70%. Several properties of magnetars indicate the presence of persistent, static currents flowing across the stellar surface and closing within the magnetosphere. The charges supporting these currents supply a significant optical depth to resonant cyclotron scattering in the 1-100 keV band. Here we describe a Monte Carlo approach to calculating the redistribution of thermal seed photons in frequency and angle by multiple resonant scattering in the magnetosphere. The calculation includes the full angular dependence of the cyclotron scattering cross section, the relativistic Doppler effect due to the motion of the charges, and allows for an arbitrary particle velocity distribution and magnetic field geometry. We construct synthetic spectra and pulse profiles for arbitrary orientations of the spin axis, magnetic axis, and line of sight, using a self-similar, twisted dipole field geometry, and assuming that the seed photons are supplied by single-temperature black body emission from the stellar surface. Pulse profiles and 1-10 keV spectra typical of AXPs are easily produced by this model, with pulsed fractions of ~ 50%. However, this model cannot reproduce the hard, rising energy spectra that are observed from SGRs during periods of activity, without overproducing the thermal emission peak. This suggests that the 1-100 keV emission of SGRs has a common origin with the hard X-ray emission detected from some AXPs above ~20 keV.

The equilibrium composition of neutron star matter is achieved through weak interactions (direct and inverse beta decays), which proceed on relatively long time scales. If the density of a matter element is perturbed, it will relax to the new chemical equilibrium through non-equilibrium reactions, which produce entropy that is partly released through neutrino emission, while a similar fraction heats the matter and is eventually radiated as thermal photons. We examined two possible mechanisms causing such density perturbations: 1) the reduction in centrifugal force caused by spin-down (particularly in millisecond pulsars), leading to "rotochemical heating", and 2) a hypothetical time-variation of the gravitational constant, as predicted by some theories of gravity and current cosmological models, leading to "gravitochemical heating". If only slow weak interactions are allowed in the neutron star (modified Urca reactions, with or without Cooper pairing), rotochemical heating can account for the observed ultraviolet emission from the closest millisecond pulsar, PSR J0437-4715, which also provides a constraint on |dG/dt| of the same order as the best available in the literature.

A hypothetical time-variation of the gravitational constant $G$ would cause neutron star matter to depart from beta equilibrium, due to the changing hydrostatic equilibrium. This forces non-equilibrium beta processes to occur, which release energy that is invested partly in neutrino emission and partly in heating the stellar interior. Eventually, the star arrives at a stationary state in which the temperature remains nearly constant, as the forcing through the change of $G$ is balanced by the ongoing reactions. Comparing the surface temperature of the nearest millisecond pulsar, PSR J0437-4715, inferred from ultraviolet observations, with our predicted stationary temperature, we estimate two upper limits for this variation: (1) $|\dot G/G| < 2 \times 10^{-10}$ yr$^{-1}$, if we allow direct Urca reactions operating in the neutron star core, and (2) $|\dot G/G| < 4 \times 10^{-12}$ yr$^{-1}$, considering only modified Urca reactions. Both results are competitive with those obtained by other methods, with (2) being among the most restrictive.

The electrostatic potential that keeps approximate charge neutrality in neutron star matter is self-consistently introduced into the formalism for rotochemical heating presented in a previous paper by Fernandez and Reisenegger. Although the new formalism is more rigorous, we show that its observable consequences are indistinguishable from those of the previous one, leaving the conclusions of the previous paper unchanged.

We have measured the expansion velocities and proper motion of the ansae in NGC7009 using high dispersion echelle spectra and archive narrow band HST images. Assuming that the ansae are moving at equal and opposite velocities from the central star we obtain an average system radial velocity of -54 +-2 km/s, the eastern ansa approaching and the western ansa receding at Vr = 5.5 +-1 km/s relative to this value. Only the proper motion of the eastern ansa could be measured, leading to 2.8 +- 0.8 arcsec/century, or Vt=130 +- 40.d km/s, where d is the distance to the nebula in kpc. Additionally, the electron temperature and density for each ansa was measured using line intensity ratios. The results are Te = 9000 K and ne = 2000 /cm3 for both ansae within the errors.