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

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to ≲0.2 dex over the full mass range, 108<M∗/M⊙≲1011, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for M∗≳109M⊙, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

/pdf/the-eagle-project-simulating-the-evolution-and-assembly-of-2fhhyaoie1.pdf

The EAGLE project: Simulating the evolution and assembly of galaxies and their environments

Supermassive black holes (BHs) have been found in 85 galaxies by dynamical modeling of spatially resolved kinematics. The Hubble Space Telescope revolutionized BH research by advancing the subject from its proof-of-concept phase into quantitative studies of BH demographics. Most influential was the discovery of a tight correlation between BH mass and the velocity dispersion σ of the bulge component of the host galaxy. Together with similar correlations with bulge luminosity and mass, this led to the widespread belief that BHs and bulges coevolve by regulating each other's growth. Conclusions based on one set of correlations from in brightest cluster ellipticals to in the smallest galaxies dominated BH work for more than a decade. New results are now replacing this simple story with a richer and more plausible picture in which BHs correlate differently with different galaxy components. A reasonable aim is to use this progress to refine our understanding of BH-galaxy coevolution. BHs with masses of 105−106M...

https://www.annualreviews.org/doi/suppl/10.1146/annurev-astro-082708-101811/suppl_file/aa51_kormendy_supmat.pdf

Coevolution (Or Not) of Supermassive Black Holes and Host Galaxies

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51 and 0.61. We measure the angular diameter distance and Hubble parameter H from the baryon acoustic oscillation (BAO) method, in combination with a cosmic microwave background prior on the sound horizon scale, after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product D_MH from the Alcock–Paczynski (AP) effect and the growth of structure, quantified by fσ_8(z), from redshift-space distortions (RSD). We combine individual measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method; in particular, the AP measurement from sub-BAO scales sharpens constraints from post-reconstruction BAOs by breaking degeneracy between D_M and H. Combined with Planck 2016 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature Ω_K = 0.0003 ± 0.0026 and a dark energy equation-of-state parameter w = −1.01 ± 0.06, in strong affirmation of the spatially flat cold dark matter (CDM) model with a cosmological constant (ΛCDM). Our RSD measurements of fσ_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H_0 = 67.3 ± 1.0 km s^−1 Mpc^−1 even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H_0 = 67.8 ± 1.2 km s^−1 Mpc^−1. Assuming flat ΛCDM, we find Ω_m = 0.310 ± 0.005 and H_0 = 67.6 ± 0.5 km s^−1 Mpc^−1, and we find a 95 per cent upper limit of 0.16 eV c^−2 on the neutrino mass sum.

/pdf/the-clustering-of-galaxies-in-the-completed-sdss-iii-baryon-9u22inggxj.pdf

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological analysis of the DR12 galaxy sample

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

We present a novel implementation of Smoothed Particle Hydrodynamics (SPHS) that uses the spatial derivative of the velocity divergence as a higher order dissipation switch. Our switch -- which is second order accurate -- detects flow convergence before it occurs. If particle trajectories are going to cross, we switch on the usual SPH artificial viscosity, as well as conservative dissipation in all advected fluid quantities (for example, the entropy). The viscosity and dissipation terms (that are numerical errors) are designed to ensure that all fluid quantities remain single-valued as particles approach one another, to respect conservation laws, and to vanish on a given physical scale as the resolution is increased. SPHS alleviates a number of known problems with `classic' SPH, successfully resolving mixing, and recovering numerical convergence with increasing resolution. An additional key advantage is that -- treating the particle mass similarly to the entropy -- we are able to use multimass particles, giving significantly improved control over the refinement strategy. We present a wide range of code tests including the Sod shock tube, Sedov-Taylor blast wave, Kelvin-Helmholtz Instability, the `blob test', and some convergence tests. Our method performs well on all tests, giving good agreement with analytic expectations.

SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch

The origin of the gas in between the Magellanic Clouds (MCs) -- known as the `Magellanic Bridge' (MB) -- is puzzling. Numerical simulations suggest that the MB formed from tidally stripped gas and stars in a recent interaction between the MCs. However, the apparent lack of stripped intermediate- or old-age stars associated with the MB is at odds with this picture. In this paper, we present the first results from the MAGellanic Inter-Cloud program (MAGIC) aimed at probing the stellar populations in the inter-cloud region. We present observations of the stellar populations in two large fields located in between the Large and Small Magellanic Clouds (LMC/SMC), secured using the WFI camera on the 2.2 m telescope in La Silla. Using a synthetic color-magnitude diagram (CMD) technique, we present the first quantitative evidence for the presence of intermediate-age and old stars in the inter-Cloud region. The intermediate-age stars -- that make up ~28% of all stars in the region -- are not present in fields at a similar distance from the SMC in a direction pointing away from the LMC. This provides potential evidence that these intermediate-age stars could have been tidally stripped from the SMC. However, spectroscopic studies will be needed to confirm or rule out the tidal origin for the inter-Cloud gas and stars.

The MAGellanic Inter-Cloud project (MAGIC) I: Evidence for intermediate-age stellar populations in between the Magellanic Clouds

The origin of the gas in between the Magellanic Clouds (MCs), known as the Magellanic Bridge (MB), has always been the subject of controversy. To shed light into this, we present the results from the MAGIC II project aimed at probing the stellar populations in ten large fields located perpendicular to the main ridge-line of HI in the Inter-Cloud region. We secured these observations of the stellar populations in between the MCs using the WFI camera on the 2.2 m telescope in La Silla. Using colour-magnitude diagrams (CMDs), we trace stellar populations across the Inter-Cloud region. In good agreement with MAGIC I, we find significant intermediate-age stars in the Inter-Cloud region as well as young stars of a similar age to the last pericentre passage in between the MCs (~200 Myr ago). We show here that the young, intermediate-age and old stars have distinct spatial distributions. The young stars correlate well with the HI gas suggesting that they were either recently stripped from the SMC or formed in-situ. The bulk of intermediate-age stars are located mainly in the bridge region where the HI column density is higher, but they are more spread out than the young stars. They have very similar properties to stars located ~2 Kpc from the SMC centre, suggesting that they were tidally stripped from this region. Finally, the old stars extend to some 8 Kpc from the SMC supporting the idea that all galaxies have a large extended metal poor stellar halo.

The Magellanic Inter-Cloud Project (MAGIC) II: Slicing up the Bridge

We revisit systematics in determining the local dark matter density (rho_dm) from the vertical motion of stars in the Solar Neighbourhood. Using a simulation of a Milky Way-like galaxy, we determine the data-quality required to detect the dark matter density at its expected local value. We introduce a new method for recovering rho_dm that uses moments of the Jeans equations, combined with a Monte Carlo Markov Chain technique to marginalise over the unknown parameters. Given sufficiently good data, we show that our method can recover the correct local dark matter density even in the face of disc inhomogeneities, non-isothermal tracers and a non-separable distribution function. We illustrate the power of our technique by applying it to Hipparcos data [Holmberg & Flynn 2000,2004]. We first make the assumption that the A and F star tracer populations are isothermal. This recovers rho_dm=0.003^{+0.009}_{-0.007}Msun/pc^3 (with 90 per cent confidence), consistent with previous determinations. However, the vertical dispersion profile of these tracers is poorly known. If we assume instead a non-isothermal profile similar to the blue disc stars from SDSS DR-7 [Abazajian et al. 2009] measured by Bond et al. (2009), we obtain a fit with a very similar chi^2 value, but with rho_dm=0.033^{+0.008}_{-0.009}Msun/pc^3 (with 90 per cent confidence). This highlights that it is vital to measure the vertical dispersion profile of the tracers to recover an unbiased estimate of the local dark matter density.

/pdf/limits-on-the-local-dark-matter-density-4e75rk78ya.pdf

Limits on the local dark matter density

We examine systematic problems in determining the local matter density from the vertical motion of stars, i.e. the 'Oort limit'. Using collisionless simulations and a Monte Carlo Markov Chain technique, we determine the data quality required to detect local dark matter at its expected density. We find that systematic errors are more important than observational errors and apply our technique to Hipparcos data to reassign realistic error bars to the local dark matter density.

/pdf/measuring-the-local-dark-matter-density-48dsgxxpex.pdf

Justin I. Read

Papers

SPHS: Smoothed Particle Hydrodynamics with a higher order dissipation switch

The MAGellanic Inter-Cloud project (MAGIC) I: Evidence for intermediate-age stellar populations in between the Magellanic Clouds

The Magellanic Inter-Cloud Project (MAGIC) II: Slicing up the Bridge

Limits on the local dark matter density

Measuring the local dark matter density