On September 14, 2015 at 09:50:45 UTC the two detectors of the Laser Interferometer Gravitational-Wave Observatory simultaneously observed a transient gravitational-wave signal. The signal sweeps upwards in frequency from 35 to 250 Hz with a peak gravitational-wave strain of 1.0×10(-21). It matches the waveform predicted by general relativity for the inspiral and merger of a pair of black holes and the ringdown of the resulting single black hole. The signal was observed with a matched-filter signal-to-noise ratio of 24 and a false alarm rate estimated to be less than 1 event per 203,000 years, equivalent to a significance greater than 5.1σ. The source lies at a luminosity distance of 410(-180)(+160)  Mpc corresponding to a redshift z=0.09(-0.04)(+0.03). In the source frame, the initial black hole masses are 36(-4)(+5)M⊙ and 29(-4)(+4)M⊙, and the final black hole mass is 62(-4)(+4)M⊙, with 3.0(-0.5)(+0.5)M⊙c(2) radiated in gravitational waves. All uncertainties define 90% credible intervals. These observations demonstrate the existence of binary stellar-mass black hole systems. This is the first direct detection of gravitational waves and the first observation of a binary black hole merger.

/pdf/the-observation-of-gravitational-waves-from-a-binary-black-58srpupgg9.pdf

The Observation of Gravitational Waves from a Binary Black Hole Merger

to be done in this area. Face recognition is a problem that scarcely clamoured for attention before the computer age but, having surfaced, has involved a wide range of techniques and has attracted the attention of some fine minds (David Mumford was a Fields Medallist in 1974). This singular application of mathematical modelling to a messy applied problem of obvious utility and importance but with no unique solution is a pretty one to share with students: perhaps, returning to the source of our opening quotation, we may invert Duncan's earlier observation, 'There is an art to find the mind's construction in the face!'.

Linear and nonlinear waves, by G. B. Whitham. Pp.636. £50. 1999. ISBN 0 471 35942 4 (Wiley).

Using a two-dimensional hydrodynamics code (PROMETHEUS), we explore the continued evolution of rotating helium stars, Mα 10 M☉, in which iron-core collapse does not produce a successful outgoing shock but instead forms a black hole of 2-3 M☉. The model explored in greatest detail is the 14 M☉ helium core of a 35 M☉ main-sequence star. The outcome is sensitive to the angular momentum. For j16 ≡ j/(1016 cm2 s-1) 3, material falls into the black hole almost uninhibited. No outflows are expected. For j16 20, the infalling matter is halted by centrifugal force outside 1000 km where neutrino losses are negligible. The equatorial accretion rate is very low, and explosive oxygen burning may power a weak equatorial explosion. For 3 j16 20, however, a reasonable value for such stars, a compact disk forms at a radius at which the gravitational binding energy can be efficiently radiated as neutrinos or converted to beamed outflow by magnetohydrodynamical (MHD) processes. These are the best candidates for producing gamma-ray bursts (GRBs). Here we study the formation of such a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter α ≈ 0.001 and 0.1. Infall along the rotational axis is initially uninhibited, and an evacuated channel opens during the first few seconds. Meanwhile the black hole is spun up by the accretion (to a ≈ 0.9), and energy is dissipated in the disk by MHD processes and radiated by neutrinos. For the α = 0.1 model, appreciable energetic outflows develop between polar angles of 30° and 45°. These outflows, powered by viscous dissipation in the disk, have an energy of up to a few times 1051 ergs and a mass ~1 M☉ and are rich in 56Ni. They constitute a supernova-like explosion by themselves. Meanwhile accretion through the disk is maintained for approximately 10-20 s but is time variable (±30%) because of hydrodynamical instabilities at the outer edge in a region where nuclei are experiencing photodisintegration. Because the efficiency of neutrino energy deposition is sensitive to the accretion rate, this instability leads to highly variable energy deposition in the polar regions. Some of this variability, which has significant power at 50 ms and overtones, may persist in the time structure of the burst. During the time followed, the average accretion rate for the standard α = 0.1 and j16 = 10 model is 0.07 M☉ s-1. The total energy deposited along the rotational axes by neutrino annihilation is (1-14) × 1051 ergs, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of energy in the polar regions, at a constant rate of 5 × 1050 ergs s-1 per pole, results in strong relativistic outflow jets beamed to about 1% of the sky. These jets may be additionally modulated by instabilities in the sides of the "nozzle" through which they flow. The jets blow aside the accreting material, remain highly focused, and are capable of penetrating the star in ~10 s. After the jet breaks through the surface of the star, highly relativistic flow can emerge. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible, ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that collapsars will always make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to produce in this model and may require merging neutron stars and black holes for their explanation.

https://iopscience.iop.org/article/10.1086/307790/pdf

Collapsars: Gamma-ray bursts and explosions in 'failed supernovae'

Monthly Notices is one of the three largest general primary astronomical research publications. It is an international journal, published by the Royal Astronomical Society. This article 1 describes its publication policy and practice.

Monthly Notices of the Royal Astronomical Society

How massive stars die-what sort of explosion and remnant each produces-depends chiefly on the masses of their helium cores and hydrogen envelopes at death. For single stars, stellar winds are the only means of mass loss, and these are a function of the metallicity of the star. We discuss how metallicity, and a simplified prescription for its effect on mass loss, affects the evolution and final fate of massive stars. We map, as a function of mass and metallicity, where black holes and neutron stars are likely to form and where different types of supernovae are produced. Integrating over an initial mass function, we derive the relative populations as a function of metallicity. Provided that single stars rotate rapidly enough at death, we speculate on stellar populations that might produce gamma-ray bursts and jet-driven supernovae.

/pdf/how-massive-single-stars-end-their-life-3v2i056c9w.pdf

How Massive Single Stars End Their Life

We present preliminary results of a time-dependent, three-dimensional hydrodynamic simulation of LMC-X4, an HMXB known to be undergoing RLOF. The simulation is initialized with the collapsed companion embedded in the undisturbed primary wind. The primary is in contact with the Roche surface, although no tidal stream or accretion disk is initialed; they are allowed to form independently. Several features of general interest to disk-fed HMXBs are apparent in the simulation. First, the primary immediately develops a compressed-wind disk in the orbital plane. This may be a natural result in most disk-fed HMXBs. Fora circularized system in an orbit close enough for RLOF to take place, we may expect the primary to be in corotation. The surface velocity may then be a significant fraction of the breakup velocity, leading to a compressed-wind disk.

/pdf/3-d-hydro-dynamic-simulation-of-accretion-disk-formation-in-4z5d2x10kc.pdf

3-D Hydro dynamic Simulation of Accretion Disk Formation in LMC-X4

The mechanism for core collapse supernova explosions remains undefined in detail and perhaps even in broad brush. Past multidimensional simulations point to the important role neutrino transport, fluid instabilities, rotation, and magnetic fields play, or may play, in generating core collapse supernova explosions, but the fundamental question as to whether or not these events are powered by neutrinos with the aid of some or all of these other phenomena or by magnetic fields or by a combination of both is unanswered. Here we present the results from two sets of simulations, in two and three spatial dimensions. In two dimensions, the simulations include multifrequency flux‐limited diffusion neutrino transport in the “ray‐by‐ray‐plus” approximation, two‐dimensional self gravity in the Newtonian limit, and nuclear burning through a 14‐isotope alpha network. The three‐dimensional simulations are model simulations constructed to reflect the post stellar core bounce conditions during neutrino shock reheating at the onset of explosion. They are hydrodynamics‐only models that focus on critical aspects of the shock stability and dynamics and their impact on the supernova mechanism and explosion. The two‐dimensional simulations demonstrate the important role nuclear burning may play despite the relatively small total energy deposition behind the shock. The three‐dimensional simulations demonstrate the need for three‐dimensional multi‐physics core collapse supernova models. In two dimensions, with the inclusion of nuclear burning, we obtain explosions (although in one case weak) for two progenitors (11 and 15 M⊙ models). Moreover, in both cases the explosion is initiated when the inner edge of the oxygen layer accretes through the shock. Thus, the shock is not revived while in the iron core, as previously discussed in the literature. The three‐dimensional studies of the development of the stationary accretion shock instability (SASI) demonstrate the fundamentally new dynamics allowed when simulations are performed in three spatial dimensions. The predominant l = 1 SASI mode gives way to a stable m = 1 mode, which in turn has significant ramifications for the distribution of angular momentum in the region between the shock and proto‐neutron star and, ultimately, for the spin of the remnant neutron star. Moreover, the three‐dimensional simulations make clear, given the increased number of degrees of freedom, that two‐dimensional models are severely limited by artificially imposed symmetries.

Ascertaining the Core Collapse Supernova Mechanism: An Emerging Picture?

Collimation of Astrophysical Jets: the Proto-Planetary Nebula He 3-1475

Supergiant High Mass X-ray Binary systems (sgHMXBs) consist of a massive, late type, star and a neutron star. The massive stars exhibits strong, radiatively driven, stellar winds. Wind accretion onto compact object triggers X-ray emission, which alters the stellar wind significantly. Hydrodynamic simulation has been used to study the neutron star - stellar wind interaction it two sgHMXBs: i) A heavily obscured sgHMXB (IGR J17252-3616) discovered by INTEGRAL. To account for observable quantities (i.e., absorbing column density) we have to assume a very slow wind terminal velocity of about 500 km/s and a rather massive neutron star. If confirmed in other obscured systems, this could provide a completely new stellar wind diagnostics. ii) A classical sgHMXB (Vela X-1) has been studied in depth to understand the origin of the off-states observed in this system. Among many models used to account for this observed behavior (clumpy wind, gating mechanism) we propose that self-organized criticality of the accretion stream is the likely reason for the observed behavior. In conclusion, the neutron star, in these two examples, acts very efficiently as a probe to study stellar winds.

/pdf/accretion-in-supergiant-high-mass-x-ray-binaries-1n7g7uncg4.pdf

Accretion in supergiant High Mass X-ray Binaries

We investigate core-collapse supernova (CCSN) nucleosynthesis in polar axisymmetric simulations using the multidimensional radiation hydrodynamics code CHIMERA. Computational costs have traditionally constrained the evolution of the nuclear composition in CCSN models to, at best, a 14-species $\alpha$-network. Such a simplified network limits the ability to accurately evolve detailed composition, neutronization and the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks in post-processing nucleosynthesis calculations. Limitations such as poor spatial resolution of the tracer particles, estimation of the expansion timescales, and determination of the "mass-cut" at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of these uncertainties on post-processing nucleosynthesis calculations and implications for future models.

John M. Blondin

Papers

3-D Hydro dynamic Simulation of Accretion Disk Formation in LMC-X4

Ascertaining the Core Collapse Supernova Mechanism: An Emerging Picture?

Collimation of Astrophysical Jets: the Proto-Planetary Nebula He 3-1475

Accretion in supergiant High Mass X-ray Binaries

Advancing Nucleosynthesis in Self-consistent, Multidimensional Models of Core-Collapse Supernovae