Showing papers in "Monthly Notices of the Royal Astronomical Society in 2015"

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

Abstract: 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

The clustering of the SDSS DR7 Main Galaxy Sample I: a 4 per cent distance measure at z=0.15

Abstract: We create a sample of spectroscopically identified galaxies with z < 0.2 from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), covering 6813 deg(2). Galaxies are chosen to sample the highest mass haloes, with an effective bias of 1.5, allowing us to construct 1000 mock galaxy catalogues (described in Paper II), which we use to estimate statistical errors and test our methods. We use an estimate of the gravitational potential to 'reconstruct' the linear density fluctuations, enhancing the baryon acoustic oscillation (BAO) signal in the measured correlation function and power spectrum. Fitting to these measurements, we determine D-V(z(eff) = 0.15) = (664 +/- 25)(r(d)/r(d, fid)) Mpc; this is a better than 4 per cent distance measurement. This 'fills the gap' in BAO distance ladder between previously measured local and higher redshift measurements, and affords significant improvement in constraining the properties of dark energy. Combining our measurement with other BAO measurements from Baryon Oscillation Spectroscopic Survey and 6-degree Field Galaxy Redshift Survey galaxy samples provides a 15 per cent improvement in the determination of the equation of state of dark energy and the value of the Hubble parameter at z = 0 (H-0). Our measurement is fully consistent with the Planck results and the Lambda cold dark matter concordance cosmology, but increases the tension between Planck+BAO H-0 determinations and direct H-0 measurements.

The EAGLE simulations of galaxy formation: calibration of subgrid physics and model variations

Abstract: We present results from thirteen cosmological simulations that explore the parameter space of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

A new class of accurate, mesh-free hydrodynamic simulation methods

Abstract: We present two new Lagrangian methods for hydrodynamics, in a systematic comparison with moving-mesh, smoothed particle hydrodynamics (SPH), and stationary (non-moving) grid methods. The new methods are designed to simultaneously capture advantages of both SPH and grid-based/adaptive mesh refinement (AMR) schemes. They are based on a kernel discretization of the volume coupled to a high-order matrix gradient estimator and a Riemann solver acting over the volume ‘overlap’. We implement and test a parallel, second-order version of the method with self-gravity and cosmological integration, in the code gizmo:1 this maintains exact mass, energy and momentum conservation; exhibits superior angular momentum conservation compared to all other methods we study; does not require ‘artificial diffusion’ terms; and allows the fluid elements to move with the flow, so resolution is automatically adaptive. We consider a large suite of test problems, and find that on all problems the new methods appear competitive with moving-mesh schemes, with some advantages (particularly in angular momentum conservation), at the cost of enhanced noise. The new methods have many advantages versus SPH: proper convergence, good capturing of fluid-mixing instabilities, dramatically reduced ‘particle noise’ and numerical viscosity, more accurate sub-sonic flow evolution, and sharp shock-capturing. Advantages versus non-moving meshes include: automatic adaptivity, dramatically reduced advection errors and numerical overmixing, velocity-independent errors, accurate coupling to gravity, good angular momentum conservation and elimination of ‘grid alignment’ effects. We can, for example, follow hundreds of orbits of gaseous discs, while AMR and SPH methods break down in a few orbits. However, fixed meshes minimize ‘grid noise’. These differences are important for a range of astrophysical problems.

Rotation-invariant convolutional neural networks for galaxy morphology prediction

Abstract: Measuring the morphological parameters of galaxies is a key requirement for studying their formation and evolution. Surveys such as the Sloan Digital Sky Survey have resulted in the availability of very large collections of images, which have permitted population-wide analyses of galaxy morphology. Morphological analysis has traditionally been carried out mostly via visual inspection by trained experts, which is time consuming and does not scale to large (≳104) numbers of images. Although attempts have been made to build automated classification systems, these have not been able to achieve the desired level of accuracy. The Galaxy Zoo project successfully applied a crowdsourcing strategy, inviting online users to classify images by answering a series of questions. Unfortunately, even this approach does not scale well enough to keep up with the increasing availability of galaxy images. We present a deep neural network model for galaxy morphology classification which exploits translational and rotational symmetry. It was developed in the context of the Galaxy Challenge, an international competition to build the best model for morphology classification based on annotated images from the Galaxy Zoo project. For images with high agreement among the Galaxy Zoo participants, our model is able to reproduce their consensus with near-perfect accuracy (>99 per cent) for most questions. Confident model predictions are highly accurate, which makes the model suitable for filtering large collections of images and forwarding challenging images to experts for manual annotation. This approach greatly reduces the experts’ workload without affecting accuracy. The application of these algorithms to larger sets of training data will be critical for analysing results from future surveys such as the Large Synoptic Survey Telescope.

The GALAH survey: scientific motivation

Abstract: The Galactic Archaeology with HERMES (GALAH) survey is a large high-resolution spectroscopic survey using the newly commissioned High Efficiency and Resolution Multi-Element Spectrograph (HERMES) on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~ 28 000) spectra in four passbands for 392 stars simultaneously over a 2 deg field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V ~ 14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, α-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

Comprehensive nucleosynthesis analysis for ejecta of compact binary mergers

Abstract: We present the first comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black hole (BH)–torus systems. The evolution of the BH–accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS–NS) and NS–BH merger models, producing ∼3–6 M_⊙ BHs with rotation parameters of A_BH ∼ 0.8 and tori of 0.03–0.3 M_⊙. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH–torus systems. While typically ∼20–25 per cent of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ∼1 per cent, but neutrino heating enhances the viscous ejecta significantly. Since BH–torus ejecta possess a wide distribution of electron fractions (0.1–0.6) and entropies, they produce heavy elements from A ∼ 80 up to the actinides, with relative contributions of A ≳ 130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH–torus phases can reproduce the solar abundances amazingly well for A ≳ 90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40 ≤ Z ≤ 56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS–BH mergers, might explain the composition of heavy-element deficient stars.

Galaxy formation in the Planck cosmology – I. Matching the observed evolution of star formation rates, colours and stellar masses

Abstract: We have updated the Munich galaxy formation model to the Planck first-year cosmology, while modifying the treatment of baryonic processes to reproduce recent data on the abundance and passive fractions of galaxies from z = 3 down to z = 0. Matching these more extensive and more precise observational results requires us to delay the reincorporation of wind ejecta, to lower the threshold for turning cold gas into stars, to eliminate ram-pressure stripping in halos less massive than 10 14 M⊙, and to modify our model for radio-mode feedback. These changes cure the most obvious failings of our previous models, namely the overly early formation of low-mass galaxies and the overly large fraction of them that are passive at late times. The new model reproduces the observed evolution both of the stellar mass function and of the distribution of star-formation rate at each stellar mass. Massive galaxies assemble most of their mass before z = 1 and are predominantly old and passive at z = 0, while lower mass galaxies assemble later and are predominantly blue and star-forming at z = 0. This phenomenological but physically based model allows the observations to be interpreted in terms of the eciency of the various processes t hat control the formation and evolution of galaxies as a function of their stellar mass, gas content, environment and time.

Gusty, gaseous flows of FIRE: Galactic winds in cosmological simulations with explicit stellar feedback

Abstract: We present an analysis of the galaxy-scale gaseous outflows from the Feedback in Realistic Environments (FIRE) simulations. This suite of hydrodynamic cosmological zoom simulations resolves formation of star-forming giant molecular clouds to z = 0, and features an explicit stellar feedback model on small scales. Our simulations reveal that high-redshift galaxies undergo bursts of star formation followed by powerful gusts of galactic outflows that eject much of the interstellar medium and temporarily suppress star formation. At low redshift, however, sufficiently massive galaxies corresponding to L* progenitors develop stable discs and switch into a continuous and quiescent mode of star formation that does not drive outflows far into the halo. Mass-loading factors for winds in L* progenitors are η ≈ 10 at high redshift, but decrease to η ≪ 1 at low redshift. Although lower values of η are expected as haloes grow in mass over time, we show that the strong suppression of outflows with decreasing redshift cannot be explained by mass evolution alone. Circumgalactic outflow velocities are variable and broadly distributed, but typically range between one and three times the circular velocity of the halo. Much of the ejected material builds a reservoir of enriched gas within the circumgalactic medium, some of which could be later recycled to fuel further star formation. However, a fraction of the gas that leaves the virial radius through galactic winds is never regained, causing most haloes with mass M_h ≤ 10^(12) M_⊙ to be deficient in baryons compared to the cosmic mean by z = 0.

European Pulsar Timing Array Limits On An Isotropic Stochastic Gravitational-Wave Background

Abstract: We present new limits on an isotropic stochastic gravitational-wave background (GWB) using a six pulsar dataset spanning 18 yr of observations from the 2015 European Pulsar Timing Array data release Performing a Bayesian analysis, we fit simultaneously for the intrinsic noise parameters for each pulsar, along with common correlated signals including clock, and Solar System ephemeris errors, obtaining a robust 95$\%$ upper limit on the dimensionless strain amplitude $A$ of the background of $A<30\times 10^{-15}$ at a reference frequency of $1\mathrm{yr^{-1}}$ and a spectral index of $13/3$, corresponding to a background from inspiralling super-massive black hole binaries, constraining the GW energy density to $\Omega_\mathrm{gw}(f)h^2 < 11\times10^{-9}$ at 28 nHz We also present limits on the correlated power spectrum at a series of discrete frequencies, and show that our sensitivity to a fiducial isotropic GWB is highest at a frequency of $\sim 5\times10^{-9}$~Hz Finally we discuss the implications of our analysis for the astrophysics of supermassive black hole binaries, and present 95$\%$ upper limits on the string tension, $G\mu/c^2$, characterising a background produced by a cosmic string network for a set of possible scenarios, and for a stochastic relic GWB For a Nambu-Goto field theory cosmic string network, we set a limit $G\mu/c^2<13\times10^{-7}$, identical to that set by the {\it Planck} Collaboration, when combining {\it Planck} and high-$\ell$ Cosmic Microwave Background data from other experiments For a stochastic relic background we set a limit of $\Omega^\mathrm{relic}_\mathrm{gw}(f)h^2<12 \times10^{-9}$, a factor of 9 improvement over the most stringent limits previously set by a pulsar timing array

The Illustris simulation: the evolving population of black holes across cosmic time

Abstract: We study the properties of black holes and their host galaxies across cosmic time in the Illustris simulation. Illustris is a large-scale cosmological hydrodynamical simulation which resolves a (106.5 Mpc)^3 volume with more than 12 billion resolution elements and includes state-of-the-art physical models relevant for galaxy formation. We find that the black hole mass density for redshifts z = 0–5 and the black hole mass function at z = 0 predicted by Illustris are in very good agreement with the most recent observational constraints. We show that the bolometric and hard X-ray luminosity functions of active galactic nuclei (AGN) at z = 0 and 1 reproduce observational data very well over the full dynamic range probed. Unless the bolometric corrections are largely underestimated, this requires radiative efficiencies to be on average low, ϵ_r ≲ 0.1, noting however that in our model radiative efficiencies are degenerate with black hole feedback efficiencies. Cosmic downsizing of the AGN population is in broad agreement with the findings from X-ray surveys, but we predict a larger number density of faint AGN at high redshifts than currently inferred. We also study black hole–host galaxy scaling relations as a function of galaxy morphology, colour and specific star formation rate. We find that black holes and galaxies co-evolve at the massive end, but for low mass, blue and star-forming galaxies there is no tight relation with either their central black hole masses or the nuclear AGN activity.

The merger rate of galaxies in the Illustris simulation: a comparison with observations and semi-empirical models

Abstract: We have constructed merger trees for galaxies in the Illustris simulation by directly tracking the baryonic content of subhaloes. These merger trees are used to calculate the galaxy–galaxy merger rate as a function of descendant stellar mass, progenitor stellar mass ratio, and redshift. We demonstrate that the most appropriate definition for the mass ratio of a galaxy–galaxy merger consists in taking both progenitor masses at the time when the secondary progenitor reaches its maximum stellar mass. Additionally, we avoid effects from ‘orphaned’ galaxies by allowing some objects to ‘skip’ a snapshot when finding a descendant, and by only considering mergers which show a well-defined ‘infall’ moment. Adopting these definitions, we obtain well-converged predictions for the galaxy–galaxy merger rate with the following main features, which are qualitatively similar to the halo–halo merger rate except for the last one: a strong correlation with redshift that evolves as ∼(1 + z)^(2.4–2.8), a power law with respect to mass ratio, and an increasing dependence on descendant stellar mass, which steepens significantly for descendant stellar masses greater than ∼2 × 10^(11)M_⊙. These trends are consistent with observational constraints for medium-sized galaxies (M* ≳ 10^(10) M_⊙), but in tension with some recent observations of the close pair fraction for massive galaxies (M* ≳ 10^(11) M_⊙), which report a nearly constant or decreasing evolution with redshift. Finally, we provide a fitting function for the galaxy–galaxy merger rate which is accurate over a wide range of stellar masses, progenitor mass ratios, and redshifts.

A self-consistent, absolute isochronal age scale for young moving groups in the solar neighbourhood

Abstract: We present a self-consistent, absolute isochronal age scale for young (< 200 Myr), nearby (< 100 pc) moving groups in the solar neighbourhood based on homogeneous fitting of semi-empirical pre-main-sequence model isochrones using the tau^2 maximum-likelihood fitting statistic of Naylor & Jeffries in the M_V, V-J colour-magnitude diagram. The final adopted ages for the groups are: 149+51-19 Myr for the AB Dor moving group, 24+/-3 Myr for the {\beta} Pic moving group (BPMG), 45+11-7 Myr for the Carina association, 42+6-4 Myr for the Columba association, 11+/-3 Myr for the {\eta} Cha cluster, 45+/-4 Myr for the Tucana-Horologium moving group (Tuc-Hor), 10+/-3 Myr for the TW Hya association, and 22+4-3 Myr for the 32 Ori group. At this stage we are uncomfortable assigning a final, unambiguous age to the Argus association as our membership list for the association appears to suffer from a high level of contamination, and therefore it remains unclear whether these stars represent a single population of coeval stars. Our isochronal ages for both the BPMG and Tuc-Hor are consistent with recent lithium depletion boundary (LDB) ages, which unlike isochronal ages, are relatively insensitive to the choice of low-mass evolutionary models. This consistency between the isochronal and LDB ages instills confidence that our self-consistent, absolute age scale for young, nearby moving groups is robust, and hence we suggest that these ages be adopted for future studies of these groups. Software implementing the methods described in this study is available from http: //www.astro.ex.ac.uk/people/timn/tau-squared/.

Galaxy and mass assembly (GAMA): End of survey report and data release 2

Abstract: The Galaxy And Mass Assembly (GAMA) survey is one of the largest contemporary spectroscopic surveys of low redshift galaxies. Covering an area of ∼286 deg2 (split among five survey regions) down to a limiting magnitude of r < 19.8 mag, we have collected spectra and reliable redshifts for 238 000 objects using the AAOmega spectrograph on the Anglo-Australian Telescope. In addition, we have assembled imaging data from a number of independent surveys in order to generate photometry spanning the wavelength range 1 nm–1 m. Here, we report on the recently completed spectroscopic survey and present a series of diagnostics to assess its final state and the quality of the redshift data. We also describe a number of survey aspects and procedures, or updates thereof, including changes to the input catalogue, redshifting and re-redshifting, and the derivation of ultraviolet, optical and near-infrared photometry. Finally, we present the second public release of GAMA data. In this release, we provide input catalogue and targeting information, spectra, redshifts, ultraviolet, optical and near-infrared photometry, single-component Sersic fits, stellar masses, Hα-derived star formation rates, environment information, and group properties for all galaxies with r < 19.0 mag in two of our survey regions, and for all galaxies with r < 19.4 mag in a third region (72 225 objects in total). The data base serving these data is available at http://www.gama-survey.org/.

ias15: a fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits

Abstract: We present IAS15, a 15th-order integrator to simulate gravitational dynamics. The integrator is based on a Gaus-Radau quadrature and can handle conservative as well as non-conservative forces. We develop a step-size control that can automatically choose an optimal timestep. The algorithm can handle close encounters and high-eccentricity orbits. The systematic errors are kept well below machine precision and long-term orbit integrations over 10 9 orbits show that IAS15 is optimal in the sense that it follows Brouwer’s law, i.e. the energy error behaves like a random walk. Our tests show that IAS15 is superior to a mixed-variable symplectic integrator (MVS) and other high-order integrators in both speed and accuracy. In fact, IAS15 preserves the symplecticity of Hamiltonian systems better than the commonly-used nominally symplectic integrators to which we compared it. We provide an open-source implementation of IAS15. The package comes with several easy-to-extend examples involving resonant planetary systems, Kozai-Lidov cycles, close encounters, radiation pressure, quadrupole moment, and generic damping functions that can, among other things, be used to simulate planet-disc interactions. Other non-conservative forces can be added easily.

Compaction and quenching of high-z galaxies in cosmological simulations: blue and red nuggets

Abstract: We use cosmological simulations to study a characteristic evolution pattern of high-redshift galaxies. Early, stream-fed, highly perturbed, gas-rich discs undergo phases of dissipative contraction into compact, star-forming systems (‘blue’ nuggets) at z ∼ 4–2. The peak of gas compaction marks the onset of central gas depletion and inside-out quenching into compact ellipticals (red nuggets) by z ∼ 2. These are sometimes surrounded by gas rings or grow extended dry stellar envelopes. The compaction occurs at a roughly constant specific star formation rate (SFR), and the quenching occurs at a constant stellar surface density within the inner kpc (Σ 1 ). Massive galaxies quench earlier, faster, and at a higher Σ 1 than lower mass galaxies, which compactify and attempt to quench more than once. This evolution pattern is consistent with the way galaxies populate the SFR-size–mass space, and with gradients and scatter across the main sequence. The compaction is triggered by an intense inflow episode, involving (mostly minor) mergers, counter-rotating streams or recycled gas, and is commonly associated with violent disc instability. The contraction is dissipative, with the inflow rate >SFR, and the maximum Σ 1 anticorrelated with the initial spin parameter. The central quenching is triggered by the high SFR and stellar/supernova feedback (maybe also active galactic nucleus feedback) due to the high central gas density, while the central inflow weakens as the disc vanishes. Suppression of fresh gas supply by a hot halo allows the long-term maintenance of quenching once above a threshold halo mass, inducing the quenching downsizing.

PARSEC evolutionary tracks of massive stars up to 350 M ☉ at metallicities 0.0001 ≤ Z ≤ 0.04

Abstract: We complement the PARSEC data base of stellar evolutionary tracks with new models of massive stars, from the pre-main sequence phase to the central carbon ignition. We consider a broad range of metallicities, 0.0001$\leq Z \leq$0.04 and initial masses up to $M_{\rm ini}=350\,M_\odot$. The main difference with respect to our previous models of massive stars is the adoption of a recent formalism accounting for the mass-loss enhancement when the ratio of the stellar to the Eddington luminosity, $\Gamma_e$, approaches unity. With this new formalism, the models are able to reproduce the Humphreys-Davidson limit observed in the Galactic and Large Magellanic Cloud colour-magnitude diagrams, without an ad hoc mass-loss enhancement. We also follow the predictions of recent wind models indicating that the metallicity dependence of the mass-loss rates becomes shallower when $\Gamma_e$ approaches unity. We thus find that the more massive stars may suffer from substantial mass-loss even at low metallicity. We also predict that the Humphreys-Davidson limit should become brighter at decreasing metallicity. We supplement the evolutionary tracks with new tables of theoretical bolometric corrections, useful to compare tracks and isochrones with the observations. For this purpose, we homogenize existing stellar atmosphere libraries of hot and cool stars (PoWR, ATLAS9 and Phoenix) and we add, where needed, new atmosphere models computed with WM-basic. The mass, age and metallicity grids are fully adequate to perform detailed investigations of the properties of very young stellar systems, both in local and distant galaxies. The new tracks supersede the previous old Padova models of massive stars.

The SAMI Galaxy Survey: Instrument specification and target selection

Abstract: The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi- object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume-limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar-mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes rpet < 19.4, stellar masses 107– 1012M⊙, and environments from isolated field galaxies through groups to clusters of _ 1015M⊙.

The ATLAS3D Project - XXX. Star formation histories and stellar population scaling relations of early-type galaxies

Abstract: We present the stellar population content of early-type galaxies from the ATLAS^(3D) survey. Using spectra integrated within apertures covering up to one effective radius, we apply two methods: one based on measuring line-strength indices and applying single stellar population (SSP) models to derive SSP-equivalent values of stellar age, metallicity, and alpha enhancement; and one based on spectral fitting to derive non-parametric star formation histories, mass-weighted average values of age, metallicity, and half-mass formation time-scales. Using homogeneously derived effective radii and dynamically determined galaxy masses, we present the distribution of stellar population parameters on the Mass Plane (M_(JAM), σ_e, R^(maj)_e), showing that at fixed mass, compact early-type galaxies are on average older, more metal-rich, and more alpha-enhanced than their larger counterparts. From non-parametric star formation histories, we find that the duration of star formation is systematically more extended in lower mass objects. Assuming that our sample represents most of the stellar content of today's local Universe, approximately 50 per cent of all stars formed within the first 2 Gyr following the big bang. Most of these stars reside today in the most massive galaxies (>10^(10.5) M⊙), which themselves formed 90 per cent of their stars by z ∼ 2. The lower mass objects, in contrast, have formed barely half their stars in this time interval. Stellar population properties are independent of environment over two orders of magnitude in local density, varying only with galaxy mass. In the highest density regions of our volume (dominated by the Virgo cluster), galaxies are older, alpha-enhanced, and have shorter star formation histories with respect to lower density regions.

Model-independent evidence in favour of an end to reionization by z ≈ 6

Abstract: We present new upper limits on the volume-weighted neutral hydrogen fraction, , at z~5-6 derived from spectroscopy of bright quasars. The fraction of the Lyman-alpha and Lyman-beta forests that is "dark" (with zero flux) provides the only model-independent upper limit on , requiring no assumptions about the physical conditions in the intergalactic medium or the quasar's unabsorbed UV continuum. In this work we update our previous results using a larger sample (22 objects) of medium-depth (~ few hours) spectra of high-redshift quasars obtained with the Magellan, MMT, and VLT. This significantly improves the upper bound on derived from dark pixel analysis to <= 0.04 + 0.05 (1{\sigma}) at z=5.6. These results provide robust constraints for theoretical models of reionization, and provide the strongest available evidence that reionization has completed (or is very nearly complete) by z~6.

The Canadian Cluster Comparison Project: detailed study of systematics and updated weak lensing masses

Abstract: Masses of clusters of galaxies from weak gravitational lensing analyses of ever larger samples are increasingly used as the reference to which baryonic scaling relations are compared. In this paper we revisit the analysis of a sample of 50 clusters studied as part of the Canadian Cluster Comparison Project. We examine the key sources of systematic error in cluster masses. We quantify the robustness of our shape measurements and calibrate our algorithm empirically using extensive image simulations. The source redshift distribution is revised using the latest state-of-the-art photometric redshift catalogues that include new deep near-infrared observations. None the less we find that the uncertainty in the determination of photometric redshifts is the largest source of systematic error for our mass estimates. We use our updated masses to determine b, the bias in the hydrostatic mass, for the clusters detected by Planck. Our results suggest 1 − b = 0.76 ± 0.05 (stat) ± 0.06 (syst), which does not resolve the tension with the measurements from the primary cosmic microwave background.

Ages and fundamental properties of Kepler exoplanet host stars from asteroseismology

Abstract: We present a study of 33 {\it Kepler} planet-candidate host stars for which asteroseismic observations have sufficiently high signal-to-noise ratio to allow extraction of individual pulsation frequencies. We implement a new Bayesian scheme that is flexible in its input to process individual oscillation frequencies, combinations of them, and average asteroseismic parameters, and derive robust fundamental properties for these targets. Applying this scheme to grids of evolutionary models yields stellar properties with median statistical uncertainties of 1.2\% (radius), 1.7\% (density), 3.3\% (mass), 4.4\% (distance), and 14\% (age), making this the exoplanet host-star sample with the most precise and uniformly determined fundamental parameters to date. We assess the systematics from changes in the solar abundances and mixing-length parameter, showing that they are smaller than the statistical errors. We also determine the stellar properties with three other fitting algorithms and explore the systematics arising from using different evolution and pulsation codes, resulting in 1\% in density and radius, and 2\% and 7\% in mass and age, respectively. We confirm previous findings of the initial helium abundance being a source of systematics comparable to our statistical uncertainties, and discuss future prospects for constraining this parameter by combining asteroseismology and data from space missions. Finally we compare our derived properties with those obtained using the global average asteroseismic observables along with effective temperature and metallicity, finding an excellent level of agreement. Owing to selection effects, our results show that the majority of the high signal-to-noise ratio asteroseismic {\it Kepler} host stars are older than the Sun.

Evidence of patchy hydrogen reionization from an extreme Lyα trough below redshift six

Abstract: We report the discovery of an extremely long (∼110 Mpc/h) and dark (τeff≳7) Lyα trough extending down to z≃5.5 towards the zem≃6.0 quasar ULAS J0148+0600. We use these new data in combination with Lyα forest measurements from 42 quasars at 4.5≤zem≤6.4 to conduct an updated analysis of the line-of-sight variance in the intergalactic Lyα opacity over 4≤z≤6. We find that the scatter in transmission among lines of sight near z∼6 significantly exceeds theoretical expectations for either a uniform ultraviolet background (UVB) or simple fluctuating UVB models in which the mean free path to ionizing photons is spatially invariant. The data, particularly near z≃5.6-5.8, instead require fluctuations in the volume-weighted hydrogen neutral fraction that are a factor of 3 or more beyond those expected from density variations alone. We argue that these fluctuations are most likely driven by large-scale variations in the mean free path, consistent with expectations for the final stages of inhomogeneous hydrogen reionization. Even by z≃5.6, however, a large fraction of the data are consistent with a uniform UVB, and by z∼5 the data are fully consistent with opacity fluctuations arising solely from the density field. This suggests that while reionization may be ongoing at z∼6, it has fully completed by z∼5.

Evolution of galaxy stellar masses and star formation rates in the EAGLE simulations

Abstract: We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation.

Gravitational lensing analysis of the Kilo-Degree Survey

Abstract: The Kilo-Degree Survey (KiDS) is a multi-band imaging survey designed for cosmological studies from weak lensing and photometric redshifts. It uses the European Southern Observatory VLT Survey Telescope with its wide-field camera OmegaCAM. KiDS images are taken in four filters similar to the Sloan Digital Sky Survey ugri bands. The best seeing time is reserved for deep r-band observations. The median 5σ limiting AB magnitude is 24.9 and the median seeing is below 0.7 arcsec. Initial KiDS observations have concentrated on the Galaxy and Mass Assembly (GAMA) regions near the celestial equator, where extensive, highly complete redshift catalogues are available. A total of 109 survey tiles, 1 square degree each, form the basis of the first set of lensing analyses of halo properties of GAMA galaxies. Nine galaxies per square arcminute enter the lensing analysis, for an effective inverse shear variance of 69 arcmin-2. Accounting for the shape measurement weight, the median redshift of the sources is 0.53. KiDS data processing follows two parallel tracks, one optimized for weak lensing measurement and one for accurate matched-aperture photometry (for photometric redshifts). This technical paper describes the lensing and photometric redshift measurements (including a detailed description of the Gaussian aperture and photometry pipeline), summarizes the data quality and presents extensive tests for systematic errors that might affect the lensing analyses. We also provide first demonstrations of the suitability of the data for cosmological measurements, and describe our blinding procedure for preventing confirmation bias in the scientific analyses. The KiDS catalogues presented in this paper are released to the community through http://kids.strw.leidenuniv.nl.

The unexpected diversity of dwarf galaxy rotation curves

Abstract: We examine the circular velocity profiles of galaxies in A cold dark matter (CDM) cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark-matter-dominated systems, reflecting the expected similarity of the underlying CDM haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for CDM haloes and include many galaxies where previous work claims the presence of a constant density 'core'. The 'cusp versus core' issue is thus better characterized as an 'inner mass deficit' problem than as a density slope mismatch. For several galaxies, the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the diversity in the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of 'inner mass deficit' galaxies inferred from kinematic data are incorrect.

Forged in FIRE: cusps, cores and baryons in low-mass dwarf galaxies

Abstract: We present multiple ultrahigh resolution cosmological hydrodynamic simulations of M_★ ≃ 10^(4–6.3) M_⊙ dwarf galaxies that form within two M_(vir) = 10^(9.5–10) M_⊙ dark matter halo initial conditions. Our simulations rely on the Feedback in Realistic Environments (FIRE) implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve structure on the ∼200 pc scales. The resultant galaxies sit on the M_★ versus M_(vir) relation required to match the Local Group stellar mass function via abundance matching. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (∼1 kpc) constant-density dark matter core in a cosmological simulation of an M_★ ≃ 10^(6.3) M_⊙ dwarf galaxy within a typical M_(vir) = 10^(10) M_⊙ halo – precisely the scale of interest for resolving the ‘too big to fail’ problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late (z ≲ 2), after the early epoch of cusp building mergers has ended. Our M_★ ≃ 10^4 M_⊙ dwarf retains a cuspy dark matter halo density profile that matches that of a dark-matter-only run of the same system. Though ancient, most of the stars in our ultrafaint form after reionization; the ultraviolet field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the protodwarf.

Kilonova light curves from the disc wind outflows of compact object mergers

Abstract: 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.

Active galactic nuclei flicker: an observational estimate of the duration of black hole growth phases of ∼105 yr

Abstract: We present an observational constraint for the typical active galactic nucleus (AGN) phase lifetime. The argument is based on the time lag between an AGN central engine switching on and becoming visible in X-rays, and the time the AGN then requires to photoionize a large fraction of the host galaxy. Based on the typical light travel time across massive galaxies, and the observed fraction of X-ray selected AGN without AGN-photoionized narrow lines, we estimate that the AGN phase typically lasts 10 5 years. This lifetime is short compared to the total growth time of 10 7 10 9 years estimated from e.g. the Soltan argument and implies

Baryon effects on the internal structure of ΛCDM haloes in the EAGLE simulations

Abstract: We investigate the internal structure and density profiles of haloes of mass 1010–1014 M⊙ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a Λ cold dark matter Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of haloes in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, haloes are well described by the Navarro–Frenk–White density profile at radii larger than ∼5 per cent of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in haloes of mass 1012–1013 M⊙, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass M* < 5 × 1010 M⊙. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass.