On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ∼ 1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40 − 8 + 8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 M ⊙ . An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ∼ 40 Mpc ) less than 11 hours after the merger by the One-Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ∼10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ∼ 9 and ∼ 16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC 4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta.

We present a series of high-resolution cosmological simulations 1 of galaxy formation to z = 0, spanning halo masses ∼ 10 8 − 10 13 M , and stellar masses ∼ 10 4 − 10 11 M . Our simulations include fully explicit treatment of the multi-phase ISM & stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M halo ∼ 10 12 M . We predict weak redshift evolution in the M * − M halo relation, consistent with current constraints to z > 6. We find that the M * − M halo relation is insensitive to numerical details, but is sensitive to feedback physics. Simulations with only supernova feedback fail to reproduce observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to destroy GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars, and is determined by self-regulation via stellar feedback. Feedback reduces star formation rates and produces reservoirs of gas that lead to rising late-time star formation histories, significantly different from halo accretion histories. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. These properties are not captured by common "sub-grid" wind models.

This paper presents a theoretical framework for understanding plasma turbulence in astrophysical plasmas. It is motivated by observations of electromagnetic and density fluctuations in the solar wind, interstellar medium and galaxy clusters, as well as by models of particle heating in accretion disks. All of these plasmas and many others have turbulent motions at weakly collisional and collisionless scales. The paper focuses on turbulence in a strong mean magnetic field. The key assumptions are that the turbulent fluctuations are small compared to the mean field, spatially anisotropic with respect to it and that their frequency is low compared to the ion cyclotron frequency. The turbulence is assumed to be forced at some system-specific outer scale. The energy injected at this scale has to be dissipated into heat, which ultimately cannot be accomplished without collisions. A kinetic cascade develops that brings the energy to collisional scales both in space and velocity. The nature of the kinetic cascade in various scale ranges depends on the physics of plasma fluctuations that exist there. There are four special scales that separate physically distinct regimes: the electron and ion gyroscales, the mean free path and the electron diffusion scale. In each of the scale ranges separated by these scales, the fully kinetic problem is systematically reduced to a more physically transparent and computationally tractable system of equations, which are derived in a rigorous way. In the "inertial range" above the ion gyroscale, the kinetic cascade separates into two parts: a cascade of Alfvénic fluctuations and a passive cascade of density and magnetic-fieldstrength fluctuations. The former are governed by the Reduced Magnetohydrodynamic (RMHD) equations at both the collisional and collisionless scales; the latter obey a linear kinetic equation along the (moving) field lines associated with the Alfvénic component (in the collisional limit, these compressive fluctuations become the slow and entropy modes of the conventional MHD). In the "dissipation range" below ion gyroscale, there are again two cascades: the kinetic-Alfvén-wave (KAW) cascade governed by two fluid-like Electron Reduced Magnetohydrodynamic (ERMHD) equations and a passive cascade of ion entropy fluctuations both in space and velocity. The latter cascade brings the energy of the inertial-range fluctuations that was Landau-damped at the ion gyroscale to collisional scales in the phase space and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and converted into electron heat. Kolmogorov-style scaling relations are derived for all of these cascades. The relationship between the theoretical models proposed in this paper and astrophysical applications and observations is discussed in detail.

We investigate large-scale galactic winds driven by momentum deposition. Momentum injection is provided by (1) radiation pressure produced by the continuum absorption and scattering of photons on dust grains and (2) supernovae (momentum injection by supernovae is important even if the supernovae energy is radiated away). Radiation can be produced by a starburst or AGN activity.We argue that momentum-driven winds are an efficient mechanism for feedback during the formation of galaxies. We show that above a limiting luminosity, momentum deposition from star formation can expel a significant fraction of the gas in a galaxy. The limiting, Eddington-like luminosity is L M ≃ (4 f g c/G) σ 4 , where σ is the galaxy velocity dispersion and f g is the gas fraction; the subscript M refers to momentum driving. A starburst that attains L M moderates its star formation rate and its luminosity does not increase significantly further. We argue that ellipticals attain this limit during their growth at z 1 and that this is the origin of the Faber-Jackson relation. We show that Lyman break galaxies and ultra-luminous infrared galaxies have luminosities near L M . Since these starbursting galaxies account for a significant fraction of the star formation at z 1, this supports our hypothesis that much of the observed stellar mass in early type galaxies was formed during Eddington-limited star formation.Star formation is unlikely to efficiently remove gas from very small scales in galactic nuclei, i.e., scales much smaller than that of a nuclear starburst. This gas is available to fuel a central black hole (BH). We argue that a BH clears gas out of its galactic nucleus when the luminosity of the BH itself reaches ≈ L M . This shuts off the fuel supply to the BH and may also terminate star formation in the surrounding galaxy. As a result, the BH mass is fixed to be M BH ≃ ( f g κ es /πG 2 )σ 4 , where κ es is the electron scattering opacity. This limit is in accord with the observed M BH − σ relation.

The most promising astrophysical sources of kHz gravitational waves (GWs) are the inspiral and merger of binary neutron star(NS)/black hole systems. Maximizing the scientific return of a GW detection will require identifying a coincident electromagnetic (EM) counterpart. One of the most likely sources of isotropic EM emission from compact object mergers is a supernova‐like transient powered by the radioactive decay of heavy elements synthesized in ejecta from the merger. We present the first calculations of the optical transients from compact object mergers that self‐consistently determine the radioactive heating by means of a nuclear reaction network; using this heating rate, we model the light curve with a one‐dimensional Monte Carlo radiation transfer calculation. For an ejecta mass ∼10−2 M⊙ (10−3 M⊙) the resulting light‐curve peaks on a time‐scale ∼1 d at a V‐band luminosity νLν∼ 3 × 1041 (1041) erg s−1[MV=−15(−14)]; this corresponds to an effective ‘f’ parameter ∼3 × 10−6 in the Li–Paczynski toy model. We argue that these results are relatively insensitive to uncertainties in the relevant nuclear physics and to the precise early‐time dynamics and ejecta composition. Since NS merger transients peak at a luminosity that is a factor of ∼103 higher than a typical nova, we propose naming these events ‘kilo‐novae’. Because of the rapid evolution and low luminosity of NS merger transients, EM counterpart searches triggered by GW detections will require close collaboration between the GW and astronomical communities. NS merger transients may also be detectable following a short‐duration gamma‐ray burst or ‘blindly’ with present or upcoming optical transient surveys. Because the emission produced by NS merger ejecta is powered by the formation of rare r‐process elements, current optical transient surveys can directly constrain the unknown origin of the heaviest elements in the Universe.

The cosmic origin of elements heavier than iron has long been uncertain. Theoretical modelling shows that the matter that is expelled in the violent merger of two neutron stars can assemble into heavy elements such as gold and platinum in a process known as rapid neutron capture (r-process) nucleosynthesis. The radioactive decay of isotopes of the heavy elements is predicted to power a distinctive thermal glow (a 'kilonova'). The discovery of an electromagnetic counterpart to the gravitational-wave source GW170817 represents the first opportunity to detect and scrutinize a sample of freshly synthesized r-process elements. Here we report models that predict the electromagnetic emission of kilonovae in detail and enable the mass, velocity and composition of ejecta to be derived from observations. We compare the models to the optical and infrared radiation associated with the GW170817 event to argue that the observed source is a kilonova. We infer the presence of two distinct components of ejecta, one composed primarily of light (atomic mass number less than 140) and one of heavy (atomic mass number greater than 140) r-process elements. The ejected mass and a merger rate inferred from GW170817 imply that such mergers are a dominant mode of r-process production in the Universe.

We investigate radiatively inefficient accretion flow models for Sgr A*, the supermassive black hole in our Galactic center, in light of new observational constraints. Confirmation of linear polarization in the submillimeter emission argues for accretion rates much less than the canonical Bondi rate. We consider models with low accretion rates and calculate the spectra produced by a hybrid electron population consisting of both thermal and nonthermal particles. The thermal electrons produce the submillimeter emission and can account for its linear polarization properties. As noted in previous work, the observed low-frequency radio spectrum can be explained if a small fraction (%1.5%) of the electron thermal energy resides in a soft powerlaw tail. In the innermost region of the accretion flow, turbulence and/or magnetic reconnection events may occasionally accelerate a fraction of the electrons into a harder power-law tail. We show that the synchrotron emission from these electrons, or the Compton upscattering of synchrotron photons by the same electrons, may account for the X-ray flares observed by Chandra.

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code ("FIRE-1") for consistency. Motivated by the development of more accurate numerics -including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms -and exploration of new physics (e.g. magnetic fields), we introduce "FIRE-2", an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (∼kpc) mass concentrations in massive (> L * ) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.