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

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

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

KiDS-450: cosmological parameter constraints from tomographic weak gravitational lensing

Abstract: We present cosmological parameter constraints from a tomographic weak gravitational lensing analysis of ~450deg$^2$ of imaging data from the Kilo Degree Survey (KiDS). For a flat $\Lambda$CDM cosmology with a prior on $H_0$ that encompasses the most recent direct measurements, we find $S_8\equiv\sigma_8\sqrt{\Omega_{\rm m}/0.3}=0.745\pm0.039$. This result is in good agreement with other low redshift probes of large scale structure, including recent cosmic shear results, along with pre-Planck cosmic microwave background constraints. A $2.3$-$\sigma$ tension in $S_8$ and `substantial discordance' in the full parameter space is found with respect to the Planck 2015 results. We use shear measurements for nearly 15 million galaxies, determined with a new improved `self-calibrating' version of $lens$fit validated using an extensive suite of image simulations. Four-band $ugri$ photometric redshifts are calibrated directly with deep spectroscopic surveys. The redshift calibration is confirmed using two independent techniques based on angular cross-correlations and the properties of the photometric redshift probability distributions. Our covariance matrix is determined using an analytical approach, verified numerically with large mock galaxy catalogues. We account for uncertainties in the modelling of intrinsic galaxy alignments and the impact of baryon feedback on the shape of the non-linear matter power spectrum, in addition to the small residual uncertainties in the shear and redshift calibration. The cosmology analysis was performed blind. Our high-level data products, including shear correlation functions, covariance matrices, redshift distributions, and Monte Carlo Markov Chains are available at http://kids.strw.leidenuniv.nl.

Improving the full spectrum fitting method: accurate convolution with Gauss-Hermite functions

Abstract: I start by providing an updated summary of the penalized pixel-fitting (pPXF) method, which is used to extract the stellar and gas kinematics, as well as the stellar population of galaxies, via full spectrum fitting. I then focus on the problem of extracting the kinematic when the velocity dispersion $\sigma$ is smaller than the velocity sampling $\Delta V$, which is generally, by design, close to the instrumental dispersion $\sigma_{\rm inst}$. The standard approach consists of convolving templates with a discretized kernel, while fitting for its parameters. This is obviously very inaccurate when $\sigma<\Delta V/2$, due to undersampling. Oversampling can prevent this, but it has drawbacks. Here I present a more accurate and efficient alternative. It avoids the evaluation of the under-sampled kernel, and instead directly computes its well-sampled analytic Fourier transform, for use with the convolution theorem. A simple analytic transform exists when the kernel is described by the popular Gauss-Hermite parametrization (which includes the Gaussian as special case) for the line-of-sight velocity distribution. I describe how this idea was implemented in a significant upgrade to the publicly available pPXF software. The key advantage of the new approach is that it provides accurate velocities regardless of $\sigma$. This is important e.g. for spectroscopic surveys targeting galaxies with $\sigma\ll\sigma_{\rm inst}$, for galaxy redshift determinations, or for measuring line-of-sight velocities of individual stars. The proposed method could also be used to fix Gaussian convolution algorithms used in today's popular software packages.

The mass distribution and gravitational potential of the Milky Way

Abstract: We present mass models of the Milky Way created to fit observational constraints and to be consistent with expectations from theoretical modelling. The method used to create these models is that demonstrated in our previous study, and we improve on those models by adding gas discs to the potential, considering the effects of allowing the inner slope of the halo density profile to vary, and including new observations of maser sources in the Milky Way amongst the new constraints. We provide a best-fitting model, as well as estimates of the properties of the Milky Way. Under the assumptions in our main model, we find that the Sun is R0 = 8.20 ± 0.09 kpc from the Galactic Centre, with the circular speed at the Sun being v0 = 232.8 ± 3.0 kms-1; and that the Galaxy has a total stellar mass of (54.3 ± 5.7) × 109 M⊙, a total virial mass of (1.30 ± 0.30) × 1012M⊙ and a local dark-matter density of 0.40 ± 0.04 GeV cm-3, where the quoted uncertainties are statistical. These values are sensitive to our choice of priors and constraints. We investigate systematic uncertainties, which in some cases may be larger. For example, if we weaken our prior on R0, we find it to be 7.97 ± 0.15 kpc and that v0 = 226.8 ± 4.2 kms-1.We find that most of these properties, including the local dark-matter density, are remarkably insensitive to the assumed power-law density slope at the centre of the dark-matter halo. We find that it is unlikely that the local standard of rest differs significantly from that found under assumptions of axisymmetry. We have made code to compute the force from our potential, and to integrate orbits within it, publicly available. (Less)

Simulating galaxy formation with black hole driven thermal and kinetic feedback

Abstract: The inefficiency of star formation in massive elliptical galaxies is widely believed to be caused by the interactions of an active galactic nucleus (AGN) with the surrounding gas. Achieving a sufficiently rapid reddening of moderately massive galaxies without expelling too many baryons has however proven difficult for hydrodynamical simulations of galaxy formation, prompting us to explore a new model for the accretion and feedback effects of supermassive black holes. For high-accretion rates relative to the Eddington limit, we assume that a fraction of the accreted rest mass energy heats the surrounding gas thermally, similar to the ‘quasar mode’ in previous work. For low-accretion rates, we invoke a new, pure kinetic feedback model that imparts momentum to the surrounding gas in a stochastic manner. These two modes of feedback are motivated both by theoretical conjectures for the existence of different types of accretion flows as well as recent observational evidence for the importance of kinetic AGN winds in quenching galaxies. We find that a large fraction of the injected kinetic energy in this mode thermalizes via shocks in the surrounding gas, thereby providing a distributed heating channel. In cosmological simulations, the resulting model produces red, non-star-forming massive elliptical galaxies, and achieves realistic gas fractions, black hole growth histories and thermodynamic profiles in large haloes.

H0LiCOW - V. New COSMOGRAIL time delays of HE 0435-1223: H0 to 3.8 per cent precision from strong lensing in a flat ΛCDM model

Abstract: We present a new measurement of the Hubble Constant H-0 and other cosmological parameters based on the joint analysis of three multiply imaged quasar systems with measured gravitational time delays. First, we measure the time delay of HE 0435-1223 from 13-yr light curves obtained as part of the COSMOGRAIL project. Companion papers detail the modelling of the main deflectors and line-of-sight effects, and how these data are combined to determine the time-delay distance of HE 0435-1223. Crucially, the measurements are carried out blindly with respect to cosmological parameters in order to avoid confirmation bias. We then combine the time-delay distance of HE 0435-1223 with previous measurements from systems B1608+656 and RXJ1131-1231 to create a Time Delay Strong Lensing probe (IDSL). In flat A cold dark matter (ACDM) with free matter and energy density, we find H-0 = 71.9(-3.0)(+2.4) km s(-1) Mpc(-1) and Omega(Lambda) = 0.62(-0.35)(+0.24) This measurement is completely independent of, and in agreement with, the local distance ladder measurements of H-0. We explore more general cosmological models combining TDSL with other probes, illustrating its power to break degeneracies inherent to other methods. The joint constraints from IDSL and Planck are H-0 = 69.2(-2.2)(+1.4) km s(-1) Mpc(-1), Omega(Lambda) = 0.70(-0.01)(+0.01) and Omega(k) = 0.003(-0.006)(+0.004) in open ACDM and H-0 = 79.0(-4.2)(+4.4) km s(-1) Mpc(-1), Omega(de) = 0.77(-0.03)(+0.02) and w = -1.38(-0.16)(+0.14) in flat wCDM. In combination with Planck and baryon acoustic oscillation data, when relaxing the constraints on the numbers of relativistic species we find N-eff = 3.34(-0.21)(+0.21) in N-eff Lambda CDM and when relaxing the total mass of neutrinos we find Sigma rn(nu) <= 0.182 eV in m(nu) Lambda CDM. Finally, in an open wCDM in combination with Planck and cosmic microwave background lensing, we find H-0 = 77.9(-4.2)(+5.0) km s(-1) Mpc(-1), Omega(de) = 0.77(-0.03)(+0.03), Omega(k) = -0.003(-0.004)(+0.004) and w = -1.37(-0.23)(+0.18).

GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey - I. A low-frequency extragalactic catalogue

Abstract: Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilizing the first year of observations. The catalogue covers 24 831 square degrees, over declinations south of +30° and Galactic latitudes outside 10° of the Galactic plane, excluding some areas such as theMagellanic Clouds. It contains 307 455 radio sources with 20 separate flux density measurements across 72-231 MHz, selected from a time- and frequency-integrated image centred at 200 MHz, with a resolution of ≈2 arcmin. Over the catalogued region, we estimate that the catalogue is 90 per cent complete at 170 mJy, and 50 per cent complete at 55 mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97 per cent above the detection threshold of 5σ, which itself is typically 50 mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low-frequency flux density scale of the southern sky to better than 10 per cent. This paper presents details of the flagging, imaging, mosaicking and source extraction/characterization, as well as estimates of the completeness and reliability. All source measurements and images are available online. 1 This is the first in a series of publications describing the GLEAM survey results.

The Hubble Space Telescope UV Legacy Survey of Galactic globular clusters – IX. The Atlas of multiple stellar populations

Abstract: We use high-precision photometry of red-giant-branch (RGB) stars in 57 Galactic globular clusters (GCs), mostly from the 'Hubble Space Telescope (HST ) UV Legacy Survey of Galactic GCs', to identify and characterize their multiple stellar populations. For each cluster the pseudo-two-colour diagram (or 'chromosome map') is presented, built with a suitable combination of stellar magnitudes in the F275W, F336W, F438W, and F814W filters that maximizes the separation between multiple populations. In the chromosome map of most GCs (type-I clusters), stars separate in two distinct groups that we identify with the first (1G) and the second generation (2G). This identification is further supported by noticing that 1G stars have primordial (oxygen-rich, sodium-poor) chemical composition, whereas 2G stars are enhanced in sodium and depleted in oxygen. This 1G-2G separation is not possible for a few GCs where the two sequences have apparently merged into an extended, continuous sequence. In some GCs (type-II clusters) the 1G and/or the 2G sequences appear to be split, hence displaying more complex chromosome maps. These clusters exhibit multiple subgiant branches (SGBs) also in purely optical colour-magnitude diagrams, with the fainter SGB joining into a red RGB which is populated by stars with enhanced heavy-element abundance. We measure the RGB width by using appropriate colours and pseudo-colours. When the metallicity dependence is removed, the RGB width correlates with the cluster mass. The fraction of 1G stars ranges from ~8 per cent to ~67 per cent and anticorrelates with the cluster mass, indicating that incidence and complexity of the multiple population phenomenon both increase with cluster mass.

Assisted inspirals of stellar mass black holes embedded in AGN discs: solving the ‘final au problem’

Abstract: We explore the evolution of stellar mass black hole binaries (BHBs) which are formed in the self-gravitating disks of active galactic nuclei (AGN). Hardening due to three-body scattering and gaseous drag are effective mechanisms that reduce the semi-major axis of a BHB to radii where gravitational waves take over, on timescales shorter than the typical lifetime of the AGN disk. Taking observationally-motivated assumptions for the rate of star formation in AGN disks, we find a rate of disk-induced BHB mergers ($\mathcal{R} \sim 3~{\rm yr}^{-1}~{\rm Gpc}^{-3}$, but with large uncertainties) that is comparable with existing estimates of the field rate of BHB mergers, and the approximate BHB merger rate implied by the recent Advanced LIGO detection of GW150914. BHBs formed thorough this channel will frequently be associated with luminous AGN, which are relatively rare within the sky error regions of future gravitational wave detector arrays. This channel could also possess a (potentially transient) electromagnetic counterpart due to super-Eddington accretion onto the stellar mass black hole following the merger.

The Auriga Project: the properties and formation mechanisms of disc galaxies across cosmic time

Abstract: We introduce a suite of 30 cosmological magneto-hydrodynamical zoom simulations of the formation of galaxies in isolated Milky Way mass dark haloes. These were carried out with the moving mesh code arepo, together with a comprehensive model for galaxy formation physics, including active galactic nuclei (AGN) feedback and magnetic fields, which produces realistic galaxy populations in large cosmological simulations. We demonstrate that our simulations reproduce a wide range of present-day observables, in particular, two-component disc-dominated galaxies with appropriate stellar masses, sizes, rotation curves, star formation rates and metallicities. We investigate the driving mechanisms that set present-day disc sizes/scalelengths, and find that they are related to the angular momentum of halo material. We show that the largest discs are produced by quiescent mergers that inspiral into the galaxy and deposit high-angular momentum material into the pre-existing disc, simultaneously increasing the spin of dark matter and gas in the halo. More violent mergers and strong AGN feedback play roles in limiting disc size by destroying pre-existing discs and by suppressing gas accretion on to the outer disc, respectively. The most important factor that leads to compact discs, however, is simply a low angular momentum for the halo. In these cases, AGN feedback plays an important role in limiting central star formation and the formation of a massive bulge.

A deep ALMA image of the Hubble Ultra Deep Field

Abstract: European Research Council via the award of an Advanced Grant; EC [312725]; European Research Council via the award of a Consolidator Grant; UK Science and Technology Facilities Council; FWO Pegasus Marie Curie Fellowship; European Research Council through the Advanced Grant [321302, 669253, 670193]; JSPS KAKENHI [JP15K17604]; Chulalongkorn University's CUniverse (CUAASC); Royal Society

H0LiCOW – I. H0 Lenses in COSMOGRAIL's wellspring: program overview

Abstract: Strong gravitational lens systems with time delays between the multiple images allow measurements of time-delay distances, which are primarily sensitive to the Hubble constant that is key to probing dark energy, neutrino physics and the spatial curvature of the Universe, as well as discovering new physics. We present H0LiCOW (H-0 Lenses in COSMOGRAIL's Wellspring), a program that aims to measure H-0 with <3.5 per cent uncertainty from five lens systems (B1608+ 656, RXJ1131-1231, HE 0435-1223, WFI2033-4723 and HE 1104-1805). We have been acquiring (1) time delays through COSMOGRAIL and Very Large Array monitoring, (2) high-resolution Hubble Space Telescope imaging for the lens mass modelling, (3) wide-field imaging and spectroscopy to characterize the lens environment and (4) moderate-resolution spectroscopy to obtain the stellar velocity dispersion of the lenses for mass modelling. In cosmological models with one-parameter extension to flat Lambda cold dark matter, we expect to measure H-0 to <3.5 per cent in most models, spatial curvature Omega(k) to 0.004, w to 0.14 and the effective number of neutrino species to 0.2 (1s uncertainties) when combined with current cosmic microwave background (CMB) experiments. These are, respectively, a factor of similar to 15, similar to 2 and similar to 1.5 tighter than CMB alone. Our data set will further enable us to study the stellar initial mass function of the lens galaxies, and the co-evolution of supermassive black holes and their host galaxies. This program will provide a foundation for extracting cosmological distances from the hundreds of time-delay lenses that are expected to be discovered in current and future surveys.

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: RSD measurement from the power spectrum and bispectrum of the DR12 BOSS galaxies

Abstract: We measure and analyse the bispectrum of the final, Data Release 12, galaxy sample provided by the Baryon Oscillation Spectroscopic Survey, splitting by selection algorithm into LOWZ and CMASS galaxies. The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. Combining the power spectrum, measured relative to the line-of-sight, with the spherically averaged bispectrum, we are able to constrain the product of the growth of structure parameter, $f$, and the amplitude of dark matter density fluctuations, $\sigma_8$, along with the geometric Alcock-Paczynski parameters, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch, $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon, $D_A(z)/r_s(z_d)$. After combining pre-reconstruction RSD analyses of the power spectrum monopole, quadrupole and bispectrum monopole; with post-reconstruction analysis of the BAO power spectrum monopole and quadrupole, we find $f(z_{\rm LOWZ})\sigma_8(z_{\rm LOWZ})=0.427\pm 0.056$, $D_A(z_{\rm LOWZ})/r_s(z_d)=6.60 \pm 0.13$, $H(z_{\rm LOWZ})r_s(z_d)=(11.55\pm 0.38)10^3\,{\rm kms}^{-1}$ for the LOWZ sample, and $f(z_{\rm CMASS})\sigma_8(z_{\rm CMASS})=0.426\pm 0.029$, $D_A(z_{\rm CMASS})/r_s(z_d)=9.39 \pm 0.10$, $H(z_{\rm CMASS})r_s(z_d)=(14.02\pm 0.22)10^3\,{\rm kms}^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 and DR12 measurements. Combining our dataset with {\it Planck15} we perform a null test of General Relativity (GR) through the $\gamma$-parametrisation finding $\gamma=0.733^{+0.068}_{-0.069}$, which is $\sim2.7\sigma$ away from the GR predictions.

Lyα and C iii] emission in z = 7–9 Galaxies: accelerated reionization around luminous star-forming systems?

Abstract: We discuss new Keck/MOSFIRE spectroscopic observations of four luminous galaxies at z ≃ 7–9 selected to have intense optical line emission by Roberts-Borsani et al. Previous follow-up has revealed Lyα in two of the four galaxies. Our new MOSFIRE observations confirm that Lyα is present in the entire sample. We detect Lyα emission in the galaxy COS-zs7-1, confirming its redshift as zLyα = 7.154, and we detect Lyα in EGS-zs8-2 at zLyα = 7.477, verifying an earlier tentative detection. The ubiquity of Lyα in this sample is puzzling given that the IGM is expected to be significantly neutral over 7 7 is expected to be strongly luminosity-dependent, with transmission accelerated in systems with intense star formation.

The BAHAMAS project: Calibrated hydrodynamical simulations for large-scale structure cosmology

Abstract: The evolution of the large-scale distribution of matter is sensitive to a variety of fundamental parameters that characterise the dark matter, dark energy, and other aspects of our cosmological framework. Since the majority of the mass density is in the form of dark matter that cannot be directly observed, to do cosmology with large-scale structure one must use observable (baryonic) quantities that trace the underlying matter distribution in a (hopefully) predictable way. However, recent numerical studies have demonstrated that the mapping between observable and total mass, as well as the total mass itself, are sensitive to unresolved feedback processes associated with galaxy formation, motivating explicit calibration of the feedback efficiencies. Here we construct a new suite of large-volume cosmological hydrodynamical simulations (called BAHAMAS, for BAryons and HAloes of MAssive Systems) where subgrid models of stellar and Active Galactic Nucleus (AGN) feedback have been calibrated to reproduce the present-day galaxy stellar mass function and the hot gas mass fractions of groups and clusters in order to ensure the effects of feedback on the overall matter distribution are broadly correct. We show that the calibrated simulations reproduce an unprecedentedly wide range of properties of massive systems, including the various observed mappings between galaxies, hot gas, total mass, and black holes, and represent a significant advance in our ability to mitigate the primary systematic uncertainty in most present large-scale structure tests.

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: anisotropic galaxy clustering in Fourier space

Abstract: We investigate the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 sample, which consists of 1198 006 galaxies in the redshift range 0.2 < z < 0.75 and a sky coverage of 10 252 deg(2). We analyse this data set in Fourier space, using the power-spectrum multipoles to measure redshift-space distortions simultaneously with the Alcock-Paczynski effect and the baryon acoustic oscillation scale. We include the power-spectrum monopole, quadrupole and hexadecapole in our analysis and compare our measurements with a perturbation-theory-based model, while properly accounting for the survey window function. To evaluate the reliability of our analysis pipeline, we participate in a mock challenge, which results in systematic uncertainties significantly smaller than the statistical uncertainties. While the high-redshift constraint on f sigma(8) at z(eff) = 0.61 indicates a small (similar to 1.4 sigma) deviation from the prediction of the Planck Lambda CDM (Lambda cold dark matter) model, the low-redshift constraint is in good agreement with Planck Lambda CDM. This paper is part of a set that analyses the final galaxy clustering data set from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. to produce the final cosmological constraints from BOSS.

Not so lumpy after all: modelling the depletion of dark matter subhaloes by Milky Way-like galaxies

Abstract: Among the most important goals in cosmology is detecting and quantifying small (M_(halo)≃10^(6−9) M⊙) dark matter (DM) subhaloes. Current probes around the Milky Way (MW) are most sensitive to such substructure within ∼20 kpc of the halo centre, where the galaxy contributes significantly to the potential. We explore the effects of baryons on subhalo populations in ΛCDM using cosmological zoom-in baryonic simulations of MW-mass haloes from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two haloes run using (1) DM-only (DMO), (2) full baryonic physics and (3) DM with an embedded disc potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain ∼2 × as many subhaloes within 100 kpc of the halo centre; this excess is ≳5 × within 25 kpc. At z = 0, the baryonic simulations are completely devoid of subhaloes down to 3×10^6M⊙ within 15 kpc of the MW-mass galaxy, and fewer than 20 surviving subhaloes have orbital pericentres <20 kpc. Despite the complexities of baryonic physics, the simple addition of an embedded central disc potential to DMO simulations reproduces this subhalo depletion, including trends with radius, remarkably well. Thus, the additional tidal field from the central galaxy is the primary cause of subhalo depletion. Subhaloes on radial orbits that pass close to the central galaxy are preferentially destroyed, causing the surviving population to have tangentially biased orbits compared to DMO predictions. Our method of embedding a potential in DMO simulations provides a fast and accurate alternative to full baryonic simulations, thus enabling suites of cosmological simulations that can provide accurate and statistical predictions of substructure populations.

The cosmic baryon cycle and galaxy mass assembly in the FIRE simulations

Abstract: We use cosmological simulations from the FIRE (Feedback In Realistic Environments) project to study the baryon cycle and galaxy mass assembly for central galaxies in the halo mass range Mhalo ∼ 10^(10)–10^(13) M_⊙. By tracing cosmic inflows, galactic outflows, gas recycling and merger histories, we quantify the contribution of physically distinct sources of material to galaxy growth. We show that in situ star formation fuelled by fresh accretion dominates the early growth of galaxies of all masses, while the re-accretion of gas previously ejected in galactic winds often dominates the gas supply for a large portion of every galaxy’s evolution. Externally processed material contributes increasingly to the growth of central galaxies at lower redshifts. This includes stars formed ex situ and gas delivered by mergers, as well as smooth intergalactic transfer of gas from other galaxies, an important but previously underappreciated growth mode. By z = 0, wind transfer, i.e. the exchange of gas between galaxies via winds, can dominate gas accretion on to ∼L* galaxies over fresh accretion and standard wind recycling. Galaxies of all masses re-accrete ≳50 per cent of the gas ejected in winds and recurrent recycling is common. The total mass deposited in the intergalactic medium per unit stellar mass formed increases in lower mass galaxies. Re-accretion of wind ejecta occurs over a broad range of time-scales, with median recycling times (∼100–350 Myr) shorter than previously found. Wind recycling typically occurs at the scale radius of the halo, independent of halo mass and redshift, suggesting a characteristic recycling zone around galaxies that scales with the size of the inner halo and the galaxy’s stellar component.

Dynamical modelling of the galactic bulge and bar: the Milky Way's pattern speed, stellar and dark matter mass distribution

Abstract: We construct a large set of dynamical models of the galactic bulge, bar and inner disk using the Made-to-Measure method. Our models are constrained to match the red clump giant density from a combination of the VVV, UKIDSS and 2MASS infrared surveys together with stellar kinematics in the bulge from the BRAVA and OGLE surveys, and in the entire bar region from the ARGOS survey. We are able to recover the bar pattern speed and the stellar and dark matter mass distributions in the bar region, thus recovering the entire galactic effective potential. We find a bar pattern speed of $39.0 \pm 3.5 \,\rm{km\,s^{-1}\,kpc^{-1}}$, placing the bar corotation radius at $6.1 \pm 0.5 \rm{kpc}$ and making the Milky Way bar a typical fast rotator. We evaluate the stellar mass of the long bar and bulge structure to be $M_{\rm{bar/bulge}} = 1.88 \pm 0.12 \times 10^{10} \, \rm{M}_{\odot}$, larger than the mass of disk in the bar region, $M_{\rm{inner\ disk}} = 1.29\pm0.12 \times 10^{10} \, \rm{M}_{\odot}$. The total dynamical mass in the bulge volume is $1.85\pm0.05\times 10^{10} \, \rm{M}_{\odot}$. Thanks to more extended kinematic data sets and recent measurement of the bulge IMF our models have a low dark matter fraction in the bulge of $17\%\pm2\%$. We find a dark matter density profile which flattens to a shallow cusp or core in the bulge region. Finally, we find dynamical evidence for an extra central mass of $\sim0.2\times10^{10} \,\rm{M}_{\odot}$, probably in a nuclear disk or disky pseudobulge.

The SCUBA-2 Cosmology Legacy Survey: 850 μm maps, catalogues and number counts

Abstract: We present a catalogue of similar to 3000 submillimetre sources detected (>= 3.5 sigma) at 850 mu m over similar to 5 deg(2) surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850 mu m, increasing the sample size of 850 mu m selected submillimetre galaxies by an order of magnitude. The wide 850 mu m survey component of S2CLS covers the extragalactic fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1s depth of S2CLS is 1.2 mJy beam(-1), approaching the SCUBA-2 850 mu m confusion limit, which we determine to be sigma(c) approximate to 0.8 mJy beam(-1). We measure the 850 mu m number counts, reducing the Poisson errors on the differential counts to approximately 4 per cent at S-850 approximate to 3 mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5 degrees-1 degrees scales are generally within 50 per cent of the S2CLS mean for S-850 > 3 mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2 sigma) density enhancement in GOODS-North. The observed counts are in reasonable agreement with recent phenomenological and semi-analytic models, although determining the shape of the faint-end slope (S-850 10 mJy there are approximately 10 sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850 mu m emission

Age-dating luminous red galaxies observed with the Southern African Large Telescope

Abstract: We measure a value for the cosmic expansion of $H(z) = 89 \pm 23$(stat) $\pm$ 44(syst) km s$^{-1}$ Mpc$^{-1}$ at a redshift of $z \simeq 0.47$ based on the differential age technique. This technique, also known as cosmic chronometers, uses the age difference between two redshifts for a passively evolving population of galaxies to calculate the expansion rate of the Universe. Our measurement is based on analysis of high quality spectra of Luminous Red Galaxies (LRGs) obtained with the Southern African Large Telescope (SALT) in two narrow redshift ranges of $z \simeq 0.40$ and $z \simeq 0.55$ as part of an initial pilot study. Ages were estimated by fitting single stellar population models to the observed spectra. This measurement presents one of the best estimates of $H(z)$ via this method at $z\sim0.5$ to date.

Breaking the chains: hot super-Earth systems from migration and disruption of compact resonant chains

Abstract: 'Hot super-Earths' (or 'mini-Neptunes') between one and four times Earth's size with period shorter than 100 d orbit 30-50 per cent of Sun-like stars Their orbital configuration - measured as the period ratio distribution of adjacent planets in multiplanet systems - is a strong constraint for formation models Here, we use N-body simulations with synthetic forces from an underlying evolving gaseous disc to model the formation and long-term dynamical evolution of super-Earth systems While the gas disc is present, planetary embryos grow and migrate inward to form a resonant chain anchored at the inner edge of the disc These resonant chains are far more compact than the observed super-Earth systems Once the gas dissipates, resonant chains may become dynamically unstable They undergo a phase of giant impacts that spreads the systems out Disc turbulence has no measurable effect on the outcome Our simulations match observations if a small fraction of resonant chains remain stable, while most super- Earths undergo a late dynamical instability Our statistical analysis restricts the contribution of stable systems to less than 25 per cent Our results also suggest that the large fraction of observed single-planet systems does not necessarily imply any dichotomy in the architecture of planetary systems Finally, we use the low abundance of resonances in Kepler data to argue that, in reality, the survival of resonant chains happens likely only in ~5 per cent of the cases This leads to a mystery: in our simulations only 50-60 per cent of resonant chains became unstable, whereas at least 75 per cent (and probably 90-95 per cent) must be unstable to match observations

CFHTLenS revisited: assessing concordance with Planck including astrophysical systematics

Abstract: We investigate the impact of astrophysical systematics on cosmic shear cosmological parameter constraints from the Canada–France–Hawaii Telescope Lensing Survey (CFHTLenS) and the concordance with cosmic microwave background measurements by Planck. We present updated CFHTLenS cosmic shear tomography measurements extended to degree scales using a covariance calibrated by a new suite of N-body simulations. We analyse these measurements with a new model fitting pipeline, accounting for key systematic uncertainties arising from intrinsic galaxy alignments, baryonic effects in the non-linear matter power spectrum, and photometric redshift uncertainties. We examine the impact of the systematic degrees of freedom on the cosmological parameter constraints, both independently and jointly. When the systematic uncertainties are considered independently, the intrinsic alignment amplitude is the only degree of freedom that is substantially preferred by the data. When the systematic uncertainties are considered jointly, there is no consistently strong preference in favour of the more complex models. We quantify the level of concordance between the CFHTLenS and Planck data sets by employing two distinct data concordance tests, grounded in Bayesian evidence and information theory. We find that the two data concordance tests largely agree with one another and that the level of concordance between the CFHTLenS and Planck data sets is sensitive to the exact details of the systematic uncertainties included in our analysis, ranging from decisive discordance to substantial concordance as the treatment of the systematic uncertainties becomes more conservative. The least conservative scenario is the one most favoured by the cosmic shear data, but it is also the one that shows the greatest degree of discordance with Planck. The data and analysis code are publicly available at https://github.com/sjoudaki/cfhtlens_revisited.

The Fermi blazar sequence

Abstract: We revisit the blazar sequence exploiting the complete, flux limited sample of blazars with known redshift detected by the Fermi satellite after 4 years of operations (the 3LAC sample). We divide the sources into gamma-ray luminosity bins, collect all the archival data for all blazars, and construct their spectral energy distribution (SED). We describe the average SED of blazars in the same luminosity bin through a very simple, completely phenomenological function consisting of two broken power laws connecting with a power law of fixed slope describing the radio emission. We do that separately for BL Lacs and for flat spectrum radio quasars (FSRQs) and also for all blazars together. The main results are: i) FSRQs display approximately the same SED as the luminosity increases, except for the fact that the relative importance of the high energy peak increases; ii) as a consequence, X-ray spectra of FSRQs become harder for larger luminosities; iii) BL Lacs form indeed a sequence: they become redder (i.e. the peak frequencies becomes smaller) for increasing luminosities, with a steeper gamma-ray slope and a larger dominance of the high energy peak; iv) for all blazars (BL Lacs+FSRQs) these properties becomes more prominent, as the highest luminosity bin is populated mostly by FSRQs and the lowest luminosity bin mostly by BL Lacs. This agrees with the original blazar sequence, although BL Lacs never have an average gamma-ray slope as hard as found in the original sequence. v) At high luminosities, a large fraction of FSRQs shows signs of thermal emission from the accretion disc, contributing in the optical-UV.

Fast radio burst source properties and curvature radiation model

Abstract: We use the observed properties of fast radio bursts (FRBs) and a number of general physical considerations to provide a broad-brush model for the physical properties of FRB sources and the radiation mechanism. We show that the magnetic field in the source region should be at least 10^{14} Gauss. This strong field is required to ensure that the electrons have sufficiently high ground state Landau energy so that particle collisions, instabilities, and strong electric and magnetic fields associated with the FRB radiation do not perturb electrons' motion in the direction transverse to the magnetic field and destroy their coherent motion; coherence is required by the high observed brightness temperature of FRB radiation. The electric field in the source region required to sustain particle motion for a wave period is estimated to be of order 10^{11} esu. These requirements suggest that FRBs are produced near the surface of magnetars perhaps via forced reconnection of magnetic fields to produce episodic, repeated, outbursts. The beaming-corrected energy release in these bursts is estimated to be ~10^{36} ergs, whereas the total energy in the magnetic field is at least ~10^{45} ergs. We provide a number of predictions for this model which can be tested by future observations. One of which is that short duration FRB-like bursts should exist at much higher frequencies, possibly up to optical.

Very massive stars, pair-instability supernovae and intermediate-mass black holes with the sevn code

Abstract: Understanding the link between massive ($\gtrsim 30$ M$_{\odot{}}$) stellar black holes (BHs) and their progenitor stars is a crucial step to interpret observations of gravitational-wave events. In this paper, we discuss the final fate of very massive stars (VMSs), with zero-age main sequence (ZAMS) mass $>150$ M$_{\odot}$, accounting for pulsational pair-instability supernovae (PPISNe) and for pair-instability supernovae (PISNe). We describe an updated version of our population synthesis code SEVN, in which we added stellar evolution tracks for VMSs with ZAMS mass up to $350 $M$_{\odot{}}$ and we included analytical prescriptions for PPISNe and PISNe. We use the new version of SEVN to study the BH mass spectrum at different metallicity $Z$, ranging from $Z=2.0\times 10^{-4}$ to $Z=2.0\times 10^{-2}$. The main effect of PPISNe and PISNe is to favour the formation of BHs in the mass range of the first gravitational-wave event (GW150914), while they prevent the formation of remnants with mass 60 - 120 M$_{\odot{}}$. In particular, we find that PPISNe significantly enhance mass loss of metal-poor ($Z\leq 2.0\times 10^{-3}$) stars with ZAMS mass $60\leq $M$_{\mathrm{ZAMS}}/$M$_{\odot{}}\leq 125$. In contrast, PISNe become effective only for moderately metal-poor ($Z<8.0\times 10^{-3}$) VMSs. VMSs with M$_{\rm ZAMS}\gtrsim{}220$ M$_\odot$ and $Z<10^{-3}$ do not undergo PISNe and form intermediate-mass BHs (IMBHs, with mass $\gtrsim 200 $M$_{\odot{}}$) via direct collapse.

Water loss from terrestrial planets orbiting ultracool dwarfs: implications for the planets of TRAPPIST-1

Abstract: Ultracool dwarfs (UCDs) encompass the population of extremely low mass stars (later than M6-type) and brown dwarfs.Because UCDs cool monotonically, their habitable zone (HZ) sweeps inward in time.Assuming they possess water, planets found in the HZ of UCDs have experienced a runaway greenhouse phase too hot for liquid water prior to entering the HZ.It has been proposed that such planets are desiccated by this hot early phase and enter the HZ as dry, inhospitable worlds.Here we model the water loss during this pre-HZ hot phase taking into account recent upper limits on the XUV emission of UCDs and using 1D radiation-hydrodynamic simulations.We address the whole range of UCDs but also focus on the planets b, c and d recently found around the 0.08 M⊙ dwarf TRAPPIST-1.Despite assumptions maximizing the FUV-photolysis of water and the XUV-driven escape of hydrogen, we find that planets can retain significant amounts of water in the HZ of UCDs, with a sweet spot in the 0.04-0.06 M⊙ range.We also studied the TRAPPIST-1 system using observed constraints on the XUV-flux.We found that TRAPPIST-1b and c can lose as much as 15 Earth Ocean and planet d -- which may be inside the HZ depending on its actual period -- may have lost less than 1 Earth Ocean.Depending on its initial content, they could have enough water to remain habitable.TRAPPIST-1 planets are key targets for atmospheric characterization and could provide strong constraints on the water erosion around UCDs.

Constraining the mass of light bosonic dark matter using SDSS Lyman-α forest

Abstract: If a significant fraction of the dark matter in the Universe is made of an ultra-light scalar field, named fuzzy dark matter (FDM) with a mass of the order of 10^−22 − 10^−21 eV, then its de Broglie wavelength is large enough to impact the physics of large-scale structure formation. In particular, the associated cut-off in the linear matter power spectrum modifies the structure of the intergalactic medium (IGM) at the scales probed by the Lyman-α forest of distant quasars. We study this effect by making use of dedicated cosmological simulations which take into account the hydrodynamics of the IGM. We explore heuristically the amplitude of quantum pressure for the FDM masses considered here and conclude that quantum effects should not modify significantly the non-linear evolution of matter density at the scales relevant to the measured Lyman-α flux power, and for m_a ≳ 10^ − 22 eV. We derive a scaling law between m_a and the mass of the well-studied thermal warm dark matter model that is best adapted to the Lyman-α forest data, and differs significantly from the one inferred by a simple linear extrapolation. By comparing FDM simulations with the Lyman-α flux power spectra determined from the BOSS survey, and marginalizing over relevant nuisance parameters, we exclude FDM masses in the range 10^ − 22 ≲ m_a < 2.3 × 10^ − 21 eV at 95 per cent CL. Adding higher resolution Lyman-α spectra extends the exclusion range up to 2.9 × 10^−21 eV. This provides a significant constraint on FDM models tailored to solve the ‘small-scale problems’ of ΛCDM.

The clustering of galaxies in the completed SDSS-III Baryon Oscillation Spectroscopic Survey: baryon acoustic oscillations in the Fourier space

Abstract: Here, we analyse the Baryon Acoustic Oscillation (BAO) signal of the final Baryon Oscillation Spectroscopic Survey (BOSS) data release (DR12). Our analysis is performed in Fourier-space, using the power spectrum monopole and quadrupole. The dataset includes 1 198 006 galaxies over the redshift range 0.2 < z < 0.75. We divide this dataset into three (overlapping) redshift bins with the effective redshifts zeff = 0.38, 0.51 and 0.61. We demonstrate the reliability of our analysis pipeline using N-body simulations as well as 1000 MultiDark-Patchy mock catalogues, which mimic the BOSS-DR12 target selection. We apply density eld reconstruction to enhance the BAO signal-to-noise ratio. By including the power spectrum quadrupole we can sep-arate the line-of-sight and angular modes, which allows us to constrain the angular diameter distance DA(z) and the Hubble parameter H ( z ) separately. We obtain two independent 1 : 6% and 1 : 5% constraints on DA(z) and 2.9% and 2.3% constraints on H(z) for the low (zeff = 0.38) and high (zeff = 0.61) redshift bin, respectively. We obtain two independent 1% and 0.9% constraints on the angular averaged distance DV(z), when ignoring the Alcock-Paczynski e ect. The detection significance of the BAO signal is ofmore » the order of 8σ (post-reconstruction) for each of the three redshift bins. Our results are in good agreement with the Planck prediction within CDM. This paper is part of a set that analyses the final galaxy clustering dataset from BOSS. The measurements and likelihoods presented here are combined with others in Alam et al. (2016) to produce the final cosmological constraints from BOSS.« less

Supernova simulations from a 3D progenitor model - Impact of perturbations and evolution of explosion properties

Abstract: We study the impact of large-scale perturbations from convective shell burning on the core-collapse supernova explosion mechanism using three-dimensional (3D) multi-group neutrino hydrodynamics simulations of an 18 solar mass progenitor. Seed asphericities in the O shell, obtained from a recent 3D model of O shell burning, help trigger a neutrino-driven explosion 330ms after bounce whereas the shock is not revived in a model based on a spherically symmetric progenitor for at least another 300ms. We tentatively infer a reduction of the critical luminosity for shock revival by ~20% due to pre-collapse perturbations. This indicates that convective seed perturbations play an important role in the explosion mechanism in some progenitors. We follow the evolution of the 18 solar mass model into the explosion phase for more than 2s and find that the cycle of accretion and mass ejection is still ongoing at this stage. With a preliminary value of 0.77 Bethe for the diagnostic explosion energy, a baryonic neutron star mass of 1.85 solar masses, a neutron star kick of ~600km/s and a neutron star spin period of ~20ms at the end of the simulation, the explosion and remnant properties are slightly atypical, but still lie comfortably within the observed distribution. Although more refined simulations and a larger survey of progenitors are still called for, this suggests that a solution to the problem of shock revival and explosion energies in the ballpark of observations are within reach for neutrino-driven explosions in 3D.