The 100 most-cited articles in Monthly Notices of the Royal Astronomical Society, 2017

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 D M 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 Hubble Fellow.

We introduce an updated physical model to simulate the formation and evolution of galaxies in cosmological, large-scale gravity+magnetohydrodynamical simulations with the moving mesh code AREPO. The overall framework builds upon the successes of the Illustris galaxy formation model, and includes prescriptions for star formation, stellar evolution, chemical enrichment, primordial and metal-line cooling of the gas, stellar feedback with galactic outflows, and black hole formation, growth and multi-mode feedback. In this paper we give a comprehensive description of the physical and numerical advances which form the core of the IllustrisTNG (The Next Generation) framework. We focus on the revised implementation of the galactic winds, of which we modify the directionality, velocity, thermal content, and energy scalings, and explore its effects on the galaxy population. As described in earlier works, the model also includes a new black hole driven kinetic feedback at low accretion rates, magnetohydrodynamics, and improvements to the numerical scheme. Using a suite of (25 Mpc $h^{-1}$)$^3$…

Hydrodynamical simulations of galaxy formation have now reached sufficient volume to make precision predictions for clustering on cosmologically relevant scales. Here, we use our new IllustrisTNG simulations to study the non-linear correlation functions and power spectra of baryons, dark matter, galaxies, and haloes over an exceptionally large range of scales. We find that baryonic effects increase the clustering of dark matter on small scales and damp the total matter power spectrum on scales up to k ∼ 10 h Mpc−1 by 20 per cent. The non-linear two-point correlation function of the stellar mass is close to a power-law over a wide range of scales and approximately invariant in time from very high redshift to the present. The two-point correlation function of the simulated galaxies agrees well with Sloan Digital Sky Survey at its mean redshift z ≃ 0.1, both as a function of stellar mass and when split according to galaxy colour, apart from a mild excess in the clustering of red galaxies in the stellar mass range of109–1010 h−2 M⊙. Given this agreement, the TNG simulations can make…

The IllustrisTNG project is a new suite of cosmological magnetohydrodynamical simulations of galaxy formation performed with the arepo code and updated models for feedback physics. Here, we introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters (1013 ≤ M200c/M⊙ ≤ 1015) at recent times (z ≤ 1). The richest clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intracluster light. Haloes more massive than about 5 × 1014 M⊙ have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial mass density distribution as shallow as the dark matter's ( − 3.5 ≲ α3D ≲ −3). Total halo mass is a very good predictor of stellar mass, and vice versa: at z = 0, the 3D stellar mass measured within 30 kpc scales as ∝(M500c)0.49 with a ∼0.12 dex scatter. This is possibly too steep in comparison to the available observational constraints, even though the abundance of…

We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colours of galaxies. We explore TNG100, a rerun of the original Illustris box, and TNG300, which includes 2 × 25003 resolution elements in a volume 20 times larger. Here, we present first results on the galaxy colour bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g − r) colours of 109 < M⋆/M⊙ < 1012.5 galaxies to the observed distribution from the Sloan Digital Sky Survey. We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement with the observations, with a sharp transition in median colour from blue to red at a characteristic M⋆ ∼ 1010.5 M⊙. Investigating the build-up of the colour–mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy colour transition is supermassive black hole feedback in its low accretion state. Across the entire population the…

We measure a value for the cosmic expansion of H(z) = 89 ± 23(stat) ± 44(syst) km s−1 Mpc−1 at a redshift of z ≃ 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 the analysis of high-quality spectra of luminous red galaxies obtained with the Southern African Large Telescope in two narrow redshift ranges of z ≃ 0.40 and 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 ∼ 0.5 to date.

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…

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 ∼ 1010–1013 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…

We present measurements of the Baryon Acoustic Oscillation (BAO) scale in redshift-space using the clustering of quasars. We consider a sample of 147,000 quasars from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) distributed over 2044 square degrees with redshifts $0.8 < z < 2.2$ and measure their spherically-averaged clustering in both configuration and Fourier space. Our observational dataset and the 1400 simulated realizations of the dataset allow us to detect a preference for BAO that is greater than 2.8$σ$. We determine the spherically averaged BAO distance to $z = 1.52$ to 3.8 per cent precision: $D_V(z=1.52)=3843\pm147 \left(r_{\rm d}/r_{\rm d, fid}\right)\ $Mpc. This is the first time the location of the BAO feature has been measured between redshifts 1 and 2. Our result is fully consistent with the prediction obtained by extrapolating the Planck flat $Λ$CDM best-fit cosmology. All of our results are consistent with basic large-scale structure (LSS) theory, confirming quasars to be a reliable tracer of LSS, and provide a starting point for numerous…

‘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…

The cosmic web is one of the most striking features of the distribution of galaxies and dark matter on the largest scales in the Universe. It is composed of dense regions packed full of galaxies, long filamentary bridges, flattened sheets and vast low-density voids. The study of the cosmic web has focused primarily on the identification of such features, and on understanding the environmental effects on galaxy formation and halo assembly. As such, a variety of different methods have been devised to classify the cosmic web – depending on the data at hand, be it numerical simulations, large sky surveys or other. In this paper, we bring 12 of these methods together and apply them to the same data set in order to understand how they compare. In general, these cosmic-web classifiers have been designed with different cosmological goals in mind, and to study different questions. Therefore, one would not a priori expect agreement between different techniques; however, many of these methods do converge on the identification of specific features. In this paper, we study the agreements and…

The APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) is an unbiased 870 µm submillimetre survey of the inner Galactic plane (|ℓ| < 60 • with |b| < 1.• 5).It is the largest and most sensitive ground-based submillimetre wavelength Galactic survey to date and has provided a large and systematic inventory of all massive, dense clumps in the Galaxy (≥1000 M at a heliocentric distance of 20 kpc) and includes representative samples of all of the earliest embedded stages of high-mass star formation.Here, we present the first detailed census of the properties (velocities, distances, luminosities and masses) and spatial distribution of a complete sample of ∼8000 dense clumps located in the Galactic disc (5 • < |ℓ| < 60 • ).We derive highly reliable velocities and distances to ∼97 per cent of the sample and use midand far-infrared survey data to develop an evolutionary classification scheme that we apply to the whole sample.Comparing the evolutionary subsamples reveals trends for increasing dust temperatures, luminosities and linewidths as a function of evolution indicating that the…

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 (H0 Lenses in COSMOGRAIL's Wellspring), a program that aims to measure H0 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 Λ cold dark matter, we expect to measure H0 to <3.5 per cent in most models, spatial curvature Ωk to 0.004, w to 0.14 and the effective number of…

Among the most important goals in cosmology is detecting and quantifying small ($M_{\rm halo}\simeq10^{6-9}~\mathrm{M}_\odot$) dark matter (DM) subhalos. Current probes around the Milky Way (MW) are most sensitive to such substructure within $\sim20$ kpc of the halo center, 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 halos from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two halos run using (1) DM-only (DMO), (2) full baryonic physics, and (3) DM with an embedded disk potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain $\sim2\times$ as many subhalos within 100 kpc of the halo center; this excess is $\gtrsim5\times$ within 25 kpc. At $z=0$, the baryonic simulations are completely devoid of subhalos down to $3\times10^6~\mathrm{M}_\odot$ within $15$ kpc of the MW-mass galaxy, and fewer than 20…

Understanding the link between massive (≳30 M⊙) 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⊙, 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⊙ 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 × 10−4 to 2.0 × 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⊙. In particular, we find that PPISNe significantly enhance mass-loss of metal-poor (Z ≤ 2.0 × 10−3) stars with ZAMS mass 60 ≤ MZAMS/ M⊙ ≤ 125. In…

We present a novel implementation of supermassive black hole (SMBH) formation, dynamics and accretion in the massively parallel tree+SPH code, ChaNGa. This approach improves the modelling of SMBHs in fully cosmological simulations, allowing for a more detailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, dynamics and sinking timescales within galaxies and SMBH accretion of rotationally supported gas. The sub-grid parameters that regulate star formation (SF) and feedback from SMBHs and SNe are optimized against a comprehensive set of z = 0 galaxy scaling relations using a novel, multidimensional parameter search. We have incorporated our new SMBH implementation and parameter optimization into a new set of high-resolution, large-scale cosmological simulations called Romulus. We present initial results from our flagship simulation, Romulus25, showing that our SMBH model results in SF efficiency, SMBH masses and global SF and SMBH accretion histories at high redshift that are consistent with…

We present a theoretical analysis of some unexplored aspects of relaxed Bose–Einstein condensate dark matter (BECDM) haloes. This type of ultralight bosonic scalar field dark matter is a viable alternative to the standard cold dark matter (CDM) paradigm, as it makes the same large-scale predictions as CDM and potentially overcomes CDM's small-scale problems via a galaxy-scale de Broglie wavelength. We simulate BECDM halo formation through mergers, evolved under the Schrödinger–Poisson equations. The formed haloes consist of a soliton core supported against gravitational collapse by the quantum pressure tensor and an asymptotic r−3 NFW-like profile. We find a fundamental relation of the core-to-halo mass with the dimensionless invariant Ξ ≡ |E|/M3/(Gm/ℏ)2 or Mc/M ≃ 2.6Ξ1/3, linking the soliton to global halo properties. For r ≥ 3.5 rc core radii, we find equipartition between potential, classical kinetic and quantum gradient energies. The haloes also exhibit a conspicuous turbulent behaviour driven by the continuous reconnection of vortex lines due to wave interference. We analyse…

Accurately predicting the demographics of dark matter (DM) substructure is of paramount importance for many fields of astrophysics, including gravitational lensing, galaxy evolution, halo occupation modelling, and constraining the nature of DM.Because of its strongly nonlinear nature, DM substructure is typically modelled using N-body simulations, which reveal that large fractions of DM subhaloes undergo complete disruption.In this paper, we use both analytical estimates and idealized numerical simulations to investigate whether this disruption is mainly physical, due to tidal heating and stripping, or numerical (i.e.artificial).We show that, contrary to naive expectation, subhaloes that experience a tidal shock E that exceeds the subhalo's binding energy, |E b |, do not undergo disruption, even when E/|E b | is as large as 100.Along the same line, and contrary to existing claims in the literature, instantaneously stripping matter from the outskirts of a DM subhalo also does not result in its complete disruption, even when the instantaneous remnant has positive binding energy.In…

We describe a new multiscale deconvolution algorithm that can also be used in a multifrequency mode. The algorithm only affects the minor clean loop. In single-frequency mode, the minor loop of our improved multiscale algorithm is over an order of magnitude faster than the casa multiscale algorithm, and produces results of similar quality. For multifrequency deconvolution, a technique named joined-channel cleaning is used. In this mode, the minor loop of our algorithm is two to three orders of magnitude faster than casa msmfs. We extend the multiscale mode with automated scale-dependent masking, which allows structures to be cleaned below the noise. We describe a new scale-bias function for use in multiscale cleaning. We test a second deconvolution method that is a variant of the moresane deconvolution technique, and uses a convex optimization technique with isotropic undecimated wavelets as dictionary. On simple well-calibrated data, the convex optimization algorithm produces visually more representative models. On complex or imperfect data, the convex optimization algorithm has…

We study the gas phase metallicity (O/H) and nitrogen abundance gradients traced by star-forming regions in a representative sample of 550 nearby galaxies in the stellar mass range 109–1011.5 M⊙ with resolved spectroscopic data from the Sloan Digital Sky Survey IV Mapping Nearby Galaxies at Apache Point Observatory survey. Using strong-line ratio diagnostics (R23 and O3N2 for metallicity and N2O2 for N/O) and referencing to the effective (half-light) radius (Re), we find that the metallicity gradient steepens with stellar mass, lying roughly flat among galaxies with log (M⋆/M⊙) = 9.0 but exhibiting slopes as steep as −0.14 dex |$R_{\rm e}^{-1}$| at log (M⋆/M⊙) = 10.5 (using R23, but equivalent results are obtained using O3N2). At higher masses, these slopes remain typical in the outer regions of our sample (R > 1.5Re), but a flattening is observed in the central regions (R < 1Re). In the outer regions (R > 2.0Re), we detect a mild flattening of the metallicity gradient in stacked profiles, although with low significance. The N/O ratio gradient provides complementary constraints on…

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 1014 G. This strong field is required to ensure that the electrons have sufficiently high ground state Landau energy so that particle collisions, instabilities and strong electromagnetic 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 the order of 1011 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…

We present Hubble Space Telescope imaging confirming the optical disappearance of the failed supernova (SN) candidate identified by Gerke et al. (2015). This $\sim 25~M_{\odot}$ red supergiant experienced a weak $\sim 10^{6}~L_{\odot}$ optical outburst in 2009 and is now at least 5 magnitudes fainter than the progenitor in the optical. The mid-IR flux has slowly decreased to the lowest levels since the first measurements in 2004. There is faint ($2000-3000~L_{\odot}$) near-IR emission likely associated with the source. We find the late-time evolution of the source to be inconsistent with obscuration from an ejected, dusty shell. Models of the spectral energy distribution indicate that the remaining bolometric luminosity is $>6$ times fainter than that of the progenitor and is decreasing as $\sim t^{-4/3}$. We conclude that the transient is unlikely to be a SN impostor or stellar merger. The event is consistent with the ejection of the envelope of a red supergiant in a failed SN and the late-time emission could be powered by fallback accretion onto a newly-formed black hole. Future…

The first four gravitational wave events detected by LIGO were all interpreted as merging black hole binaries (BHBs), opening a new perspective on the study of such systems. Here we use our new population-synthesis code mobse, an upgraded version of bse, to investigate the demography of merging BHBs. mobse includes metallicity-dependent prescriptions for mass-loss of massive hot stars. It also accounts for the impact of the electron-scattering Eddington factor on mass-loss. We perform >108 simulations of isolated massive binaries, with 12 different metallicities, to study the impact of mass-loss, core-collapse supernovae and common envelope on merging BHBs. Accounting for the dependence of stellar winds on the Eddington factor leads to the formation of black holes (BHs) with mass up to 65 M⊙ at metallicity Z ∼ 0.0002. However, most BHs in merging BHBs have masses ≲ 40 M⊙. We find merging BHBs with mass ratios in the 0.1–1.0 range, even if mass ratios >0.6 are more likely. We predict that systems like GW150914, GW170814 and GW170104 can form only from progenitors with metallicity Z ≤…

We analyze scaling relations and evolution histories of galaxy sizes in TNG100, part of the IllustrisTNG simulation suite. Observational qualitative trends of size with stellar mass, star-formation rate and redshift are reproduced, and a quantitative comparison of projected r-band sizes at 0~~10^{9.5}Msun, the evolution of the median main progenitor differs, with quenched galaxies hardly growing in median size before quenching, whereas main-sequence galaxies grow their median size continuously, thus opening a gap from the progenitors of quenched galaxies. This is partly because the main-sequence high-redshift progenitors of quenched z=0 galaxies are drawn from the lower end of the size distribution of the overall population of main-sequence high-redshift galaxies. (ii) Quenched galaxies with M_{*,z=0}>~10^{9.5}Msun experience a steep size growth on the size-mass plane after their quenching time, but with the exception of galaxies with M_{*,z=0}>~10^{11}Msun, the size growth after quenching is small in absolute terms, such that most of the size (and mass) growth of quenched galaxies…

We perform a combined analysis of cosmic shear tomography, galaxy-galaxy lensing tomography, and redshift-space multipole power spectra (monopole and quadrupole) using 450 deg 2 of imaging data by the Kilo Degree Survey (KiDS-450) overlapping with two spectroscopic surveys: the 2-degree Field Lensing Survey (2dFLenS) and the Baryon Oscillation Spectroscopic Survey (BOSS).We restrict the galaxy-galaxy lensing and multipole power spectrum measurements to the overlapping regions with KiDS, and self-consistently compute the full covariance between the different observables using a large suite of N-body simulations.We methodically analyse different combinations of the observables, finding that the galaxy-galaxy lensing measurements are particularly useful in improving the constraint on the intrinsic alignment amplitude, while the multipole power spectra are useful in tightening the constraints along the lensing degeneracy direction.The fully combined constraint on S 8 ≡ σ 8 √ m /0.3 = 0.742 ± 0.035, which is an improvement by 20 per cent compared to KiDS alone, corresponds to a 2.6σ…

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 ma 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 ma ≳ 10 − 22 eV. We derive a scaling law between ma 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…

Interpretations of exoplanetary transmission spectra have been undermined by apparent obscuration due to clouds/hazes. Debate rages on whether weak H2O features seen in exoplanet spectra are due to clouds or inherently depleted oxygen. Assertions of solar H2O abundances have relied on making a priori model assumptions, for example, chemical/radiative equilibrium. In this work, we attempt to address this problem with a new retrieval paradigm for transmission spectra. We introduce poseidon, a two-dimensional atmospheric retrieval algorithm including generalized inhomogeneous clouds. We demonstrate that this prescription allows one to break vital degeneracies between clouds and prominent molecular abundances. We apply poseidon to the best transmission spectrum presently available, for the hot Jupiter HD 209458b, uncovering new insights into its atmosphere at the day–night terminator. We extensively explore the parameter space with an unprecedented 108 models, spanning the continuum from fully cloudy to cloud-free atmospheres, in a fully Bayesian retrieval framework. We report the first…

We revisit the blazar sequence exploiting the complete, flux-limited sample of blazars with known redshift detected by the Fermi satellite after 4 yr of operations (the 3LAC sample). We divide the sources into γ-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 simple phenomenological function consisting of two broken power laws connecting with a power law 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, but the relative importance of the high-energy peak increases; (ii) as a consequence, the X-ray spectra of FSRQs become harder for larger luminosities; (iii) BL Lacs indeed form a sequence: they become redder (i.e. smaller peak frequencies) with increasing luminosities, with a softer γ-ray slope and a larger dominance of the high-energy peak; (iv)…

We describe the Next Generation Transit Survey (NGTS), which is a ground-based project searching for transiting exoplanets orbiting bright stars. NGTS builds on the legacy of previous surveys, most notably WASP, and is designed to achieve higher photometric precision and hence find smaller planets than have previously been detected from the ground. It also operates in red light, maximizing sensitivity to late K and early M dwarf stars. The survey specifications call for photometric precision of 0.1 per cent in red light over an instantaneous field of view of 100 deg2, enabling the detection of Neptune-sized exoplanets around Sun-like stars and super-Earths around M dwarfs. The survey is carried out with a purpose-built facility at Cerro Paranal, Chile, which is the premier site of the European Southern Observatory (ESO). An array of twelve 20 cm f/2.8 telescopes fitted with back-illuminated deep-depletion CCD cameras is used to survey fields intensively at intermediate Galactic latitudes. The instrument is also ideally suited to ground-based photometric follow-up of exoplanet…

We study the impact of large-scale perturbations from convective shell burning on the core-collapse supernova explosion mechanism using 3D multigroup neutrino hydrodynamics simulations of an 18M⊙ progenitor. Seed asphericities in the O shell, obtained from a recent 3D model of O shell burning, help trigger a neutrino-driven explosion 330 ms after bounce whereas the shock is not revived in a model based on a spherically symmetric progenitor for at least another 300 ms. We tentatively infer a reduction of the critical luminosity for shock revival by |$\mathord {\sim }20 \, {\rm per \, cent}$| 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 18M⊙ model into the explosion phase for more than 2 s and find that the cycle of accretion and mass ejection is still ongoing at this stage. With a preliminary value of 7.7 × 1050 erg for the diagnostic explosion energy, a baryonic neutron star mass of 1.85M⊙, a neutron star kick of |$\mathord {\sim }600 \,…

The methods in the Gasoline2 smoothed particle hydrodynamics (SPH) code are described and tested. Gasoline2 is the most recent version of the Gasoline code for parallel hydrodynamics and gravity with identical hydrodynamics to the Changa code. As with other Modern SPH codes, we prevent sharp jumps in time-steps, use upgraded kernels and larger neighbour numbers and employ local viscosity limiters. Unique features in Gasoline2 include its Geometric Density Average Force expression, explicit Turbulent Diffusion terms and Gradient-Based shock detection to limit artificial viscosity. This last feature allows Gasoline2 to completely avoid artificial viscosity in non-shocking compressive flows. We present a suite of tests demonstrating the value of these features with the same code configuration and parameter choices used for production simulations.

We explore strategies to extract cosmological constraints from a joint analysis of cosmic shear, galaxy–galaxy lensing, galaxy clustering, cluster number counts and cluster weak lensing. We utilize the cosmolike software to simulate results from a Large Synoptic Survey Telescope (LSST) like data set, specifically, we (1) compare individual and joint analyses of the different probes, (2) vary the selection criteria for lens and source galaxies, (3) investigate the impact of blending, (4) investigate the impact of the assumed cosmological model in multiprobe covariances, (6) quantify information content as a function of scales and (7) explore the impact of intrinsic galaxy alignment in a multiprobe context. Our analyses account for all cross-correlations within and across probes and include the higher-order (non-Gaussian) terms in the multiprobe covariance matrix. We simultaneously model cosmological parameters and a variety of systematics, e.g. uncertainties arising from shear and photo-z calibration, cluster mass-observable relation, galaxy intrinsic alignment and galaxy bias (up to…

We present the Pristine survey, a new narrow-band photometric survey focused on the metallicity-sensitive Ca H&K lines and conducted in the Northern hemisphere with the wide-field imager MegaCam on the Canada–France–Hawaii Telescope. This paper reviews our overall survey strategy and discusses the data processing and metallicity calibration. Additionally we review the application of these data to the main aims of the survey, which are to gather a large sample of the most metal-poor stars in the Galaxy, to further characterize the faintest Milky Way satellites, and to map the (metal-poor) substructure in the Galactic halo. The current Pristine footprint comprises over 1000 deg2 in the Galactic halo ranging from b ∼ 30° to ∼78° and covers many known stellar substructures. We demonstrate that, for Sloan Digital Sky Survey (SDSS) stellar objects, we can calibrate the photometry at the 0.02-mag level. The comparison with existing spectroscopic metallicities from SDSS/Sloan Extension for Galactic Understanding and Exploration (SEGUE) and Large Sky Area Multi-Object Fiber Spectroscopic…

We find that clouds of optically thin, pressure-confined gas are prone to fragmentation as they cool below ∼106 K. This fragmentation follows the lengthscale ∼cstcool, ultimately reaching very small scales (∼0.1 pc/n), as they reach the temperature ∼104 K at which hydrogen recombines. While this lengthscale depends on the ambient pressure confining the clouds, we find that the column density through an individual fragment Ncloudlet ∼ 1017 cm−2 is essentially independent of environment; this column density represents a characteristic scale for atomic gas at 104 K. We therefore suggest that ‘clouds’ of cold, atomic gas may, in fact, have the structure of a mist or a fog, composed of tiny fragments dispersed throughout the ambient medium. We show that this scale emerges in hydrodynamic simulations, and that the corresponding increase in the surface area may imply rapid entrainment of cold gas. We also apply it to a number of observational puzzles, including the large covering fraction of diffuse gas in galaxy haloes, the broad-line widths seen in quasar and AGN spectra and the…

Cometary ices are rich in CO2, CO and organic volatile compounds, but the carbon content of cometary dust was only measured for the Oort Cloud comet 1P/Halley, during its flyby in 1986. The COmetary Secondary Ion Mass Analyzer (COSIMA)/Rosetta mass spectrometer analysed dust particles with sizes ranging from 50 to 1000 mu m, collected over 2 yr, from 67P/Churyumov-Gerasimenko (67P), a Jupiter family comet. Here, we report 67P dust composition focusing on the elements C and O. It has a high carbon content (atomic C/Si = 5.5(-1.2)(+1.4) on average) close to the solar value and comparable to the 1P/Halley data. From COSIMA measurements, we conclude that 67P particles are made of nearly 50 per cent organic matter in mass, mixed with mineral phases that are mostly anhydrous. The whole composition, rich in carbon and non-hydrated minerals, points to a primitive matter that likely preserved its initial characteristics since the comet accretion in the outer regions of the protoplanetary disc.

We describe data release 3 (DR3) of the Galaxy And Mass Assembly (GAMA) survey. The GAMA survey is a spectroscopic redshift and multiwavelength photometric survey in three equatorial regions each of 60.0 deg2 (G09, G12, and G15), and two southern regions of 55.7 deg2 (G02) and 50.6 deg2 (G23). DR3 consists of: the first release of data covering the G02 region and of data on H-ATLAS (Herschel – Astrophysical Terahertz Large Area Survey) sources in the equatorial regions; and updates to data on sources released in DR2. DR3 includes 154 809 sources with secure redshifts across four regions. A subset of the G02 region is 95.5 per cent redshift complete to r < 19.8 mag over an area of 19.5 deg2, with 20 086 galaxy redshifts, that overlaps substantially with the XXL survey (X-ray) and VIPERS (redshift survey). In the equatorial regions, the main survey has even higher completeness (98.5 per cent), and spectra for about 75 per cent of H-ATLAS filler targets were also obtained. This filler sample extends spectroscopic redshifts, for probable optical counterparts to H-ATLAS submillimetre…

We report the detection of the Lyman continuum (LyC) radiation of the compact star-forming galaxy (SFG) J1154+2443 observed with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope. This galaxy, at a redshift of z = 0.3690, is characterized by a high emission-line flux ratio O32 = [O iii] λ5007/[O ii] λ3727 = 11.5. The escape fraction of the LyC radiation fesc(LyC) in this galaxy is 46 per cent, the highest value found so far in low-redshift SFGs and one of the highest values found in galaxies at any redshift. The narrow double-peaked Ly α emission line is detected in the spectrum of J1154+2443 with a separation between the peaks Vsep of 199 km s−1, one of the lowest known for Ly α-emitting galaxies, implying a high fesc(Ly α). Comparing the extinction-corrected Ly α/H β flux ratio with the case B value, we find fesc(Ly α) = 98 per cent. Our observations, combined with previous detections in the literature, reveal an increase of O32 with increasing fesc(LyC). We also find a tight anticorrelation between fesc(LyC) and Vsep. The surface brightness profile derived…

We study the dust content of galaxies from $z=0$ to $z=9$ in semi-analytic models of galaxy formation that include new recipes to track the production and destruction of dust. We include condensation of dust in stellar ejecta, the growth of dust in the interstellar medium (ISM), the destruction of dust by supernovae and in the hot halo, and dusty winds and inflows. The rate of dust growth in the ISM depends on the metallicity and density of molecular clouds. Our fiducial model reproduces the relation between dust mass and stellar mass from $z=0$ to $z=7$, the number density of galaxies with dust masses less than $10^{8.3}\,\rm{M}_\odot$, and the cosmic density of dust at $z=0$. The model accounts for the double power-law trend between dust-to-gas (DTG) ratio and gas-phase metallicity of local galaxies and the relation between DTG ratio and stellar mass. The dominant mode of dust formation is dust growth in the ISM, except for galaxies with $M_*<10^7\,\rm{M}_\odot$, where condensation of dust in supernova ejecta dominates. The dust-to-metal ratio of galaxies depends on the gas-phase…

We present a numerical code to simulate maps of Galactic emission in intensity and polarization at microwave frequencies, aiding in the design of cosmic microwave background experiments. This python code builds on existing efforts to simulate the sky by providing an easy-to-use interface and is based on publicly available data from the WMAP (Wilkinson Microwave Anisotropy Probe) and Planck satellite missions. We simulate synchrotron, thermal dust, free–free and anomalous microwave emission over the whole sky, in addition to the cosmic microwave background, and include a set of alternative prescriptions for the frequency dependence of each component, for example, polarized dust with multiple temperatures and a decorrelation of the signals with frequency, which introduce complexity that is consistent with current data. We also present a new prescription for adding small-scale realizations of these components at resolutions greater than current all-sky measurements. The usefulness of the code is demonstrated by forecasting the impact of varying foreground complexity on the recovered…

We present a suite of 15 cosmological zoom-in simulations of isolated dark matter halos, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_\odot$ at $z=0$, in order to understand the relationship between halo assembly, galaxy formation, and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) project and are performed at extremely high resolution. The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way ($M_\star/{\rm M}_\odot\approx 10^5-10^7$). This non-negligible spread in stellar mass at $z=0$ in halos within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity $V_{\rm max}$. Much of this dependence of $M_\star$ on a second parameter (beyond $M_{\rm halo}$) is a direct consequence of the $M_{\rm halo}\sim10^{10}\,{\rm M}_\odot$ mass scale coinciding with the threshold for strong…

We present new determinations of the stellar-to-halo mass relation (SHMR) at $z=0-10$ that match the evolution of the galaxy stellar mass function, the SFR$-M_*$ relation,and the cosmic star formation rate. We utilize a compilation of 40 observational studies from the literature and correct them for potential biases. Using our robust determinations of halo mass assembly and the SHMR, we infer star formation histories, merger rates, and structural properties for average galaxies, combining star-forming and quenched galaxies. Our main findings: (1) The halo mass $M_{50}$ above which 50\% of galaxies are quenched coincides with sSFR/sMAR$\sim1$, where sMAR is the specific halo mass accretion rate. (2) $M_{50}$ increases with redshift, presumably due to cold streams being more efficient at high redshift while virial shocks and AGN feedback become more relevant at lower redshifts. (3) The ratio sSFR/sMAR has a peak value, which occurs around $M_{\rm vir}\sim2\times10^{11}M_{\odot}$. (4) The stellar mass density within 1 kpc, $Σ_1$, is a good indicator of the galactic global sSFR. (5)…

We introduce the MOdelling Star cluster population Assembly In Cosmological Simulations within EAGLE (E-MOSAICS) project. E-MOSAICS incorporates models describing the formation, evolution, and disruption of star clusters into the EAGLE galaxy formation simulations, enabling the examination of the co-evolution of star clusters and their host galaxies in a fully cosmological context. A fraction of the star formation rate of dense gas is assumed to yield a cluster population; this fraction and the population's initial properties are governed by the physical properties of the natal gas. The subsequent evolution and disruption of the entire cluster population are followed accounting for two-body relaxation, stellar evolution, and gravitational shocks induced by the local tidal field. This introductory paper presents a detailed description of the model and initial results from a suite of 10 simulations of ∼L⋆ galaxies with disc-like morphologies at z = 0. The simulations broadly reproduce key observed characteristics of young star clusters and globular clusters (GCs), without invoking…

Abstract We introduce the Cluster-EAGLE (c-eagle) simulation project, a set of cosmological hydrodynamical zoom simulations of the formation of 30 galaxy clusters in the mass range of 1014 < M200/M⊙ < 1015.4 that incorporates the Hydrangea sample of Bahé et al. (2017). The simulations adopt the state-of-the-art eagle galaxy formation model, with a gas particle mass of 1.8 × 106 M⊙ and physical softening length of 0.7 kpc. In this paper, we introduce the sample and present the low-redshift global properties of the clusters. We calculate the X-ray properties in a manner consistent with observational techniques, demonstrating the bias and scatter introduced by using estimated masses. We find the total stellar content and black hole masses of the clusters to be in good agreement with the observed relations. However, the clusters are too gas rich, suggesting that the active galactic nucleus (AGN) feedback model is not efficient enough at expelling gas from the high-redshift progenitors of the clusters. The X-ray properties, such as the spectroscopic temperature and the soft-band…

Massive black holes (BHs) inhabit local galaxies, including the Milky Way and some dwarf galaxies. BH formation, occurring at early cosmic times, must account for the properties of BHs in today’s galaxies, notably why some galaxies host a BH, and others do not. We investigate the formation, distribution and growth of BH ‘seeds’ by using the adaptive mesh refinement code ramses. We develop an implementation of BH formation in dense, low-metallicity environments, as advocated by models invoking the collapse of the first generation of stars, or of dense nuclear star clusters. The seed masses are computed one-by-one on-the-fly, based on the star formation rate and the stellar initial mass function. This self-consistent method to seed BHs allows us to study the distribution of BHs in a cosmological context and their evolution over cosmic time. We find that all high-mass galaxies tend to host a BH, whereas low-mass counterparts have a lower probability of hosting a BH. After the end of the epoch of BH formation, this probability is modulated by the growth of the galaxy. The simulated BHs…

Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Several research groups have developed techniques independently to achieve this goal. While these methods have all provided remarkably high-precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we present for the first time a detailed comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modelling community was provided simulated cluster images that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that several properties of the lensing…

Chemical compositions of giant planets provide a means to constrain how and where they form. Traditionally, super-stellar elemental abundances in giant planets were thought to be possible due to accretion of metal-rich solids. Such enrichments are accompanied by oxygen-rich compositions (i.e. C/O below the disc's value, assumed to be solar, C/O = 0.54). Without solid accretion, the planets are expected to have sub-solar metallicity, but high C/O ratios. This arises because the solids are dominated by oxygen-rich species, e.g. H2O and CO2, which freeze out in the disc earlier than CO, leaving the gas metal poor but carbon rich. Here we demonstrate that super-solar metallicities can be achieved by gas accretion alone when growth and radial drift of pebbles are considered in protoplanetary discs. Through this mechanism, planets may simultaneously acquire super-solar metallicities and super-solar C/O ratios. This happens because the pebbles transport volatile species inwards as they migrate through the disc, enriching the gas at snow lines where the volatiles sublimate. Furthermore, the…

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not.The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

We compare the predictions of Horizon-AGN, a hydrodynamical cosmological simulation that uses an adaptive mesh refinement code, to observational data in the redshift range 0 < z < 6. We study the reproduction, by the simulation, of quantities that trace the aggregate stellar-mass growth of galaxies over cosmic time: luminosity and stellar-mass functions, the star formation main sequence, rest-frame UV–optical–near-infrared colours and the cosmic star formation history. We show that Horizon-AGN, which is not tuned to reproduce the local Universe, produces good overall agreement with these quantities, from the present day to the epoch when the Universe was 5 per cent of its current age. By comparison to Horizon-noAGN, a twin simulation without active galactic nuclei feedback, we quantify how feedback from black holes is likely to help shape galaxy stellar-mass growth in the redshift range 0 < z < 6, particularly in the most massive galaxies. Our results demonstrate that Horizon-AGN successfully captures the evolutionary trends of observed galaxies over the lifetime of the Universe,…

We introduce the Hydrangea simulations, a suite of 24 cosmological hydrodynamic zoom-in simulations of massive galaxy clusters (M200c = 1014–1015.4 M⊙) with baryon particle masses of ∼106 M⊙. Designed to study the impact of the cluster environment on galaxy formation, they are a key part of the ‘Cluster–EAGLE’ project. They use a galaxy formation model developed for the EAGLE project, which has been shown to yield both realistic field galaxies and hot gas fractions of galaxy groups consistent with observations. The total stellar mass content of the simulated clusters agrees with observations, but central cluster galaxies are too massive, by up to 0.6 dex. Passive satellite fractions are higher than in the field, and at stellar masses Mstar > 1010 M⊙, this environmental effect is quantitatively consistent with observations. The predicted satellite stellar mass function matches data from local cluster surveys. Normalized to total mass, there are fewer low-mass (Mstar ≲ 1010 M⊙) galaxies within the virial radius of clusters than in the field, primarily due to star formation quenching.…

Early observations of supernovae (SNe) indicate that enhanced mass-loss and pre-SN outbursts may occur in progenitors of many types of SNe. We investigate the role of energy transport via waves driven by vigorous convection during late-stage nuclear burning of otherwise typical 15 M⊙ red supergiant SN progenitors. Using mesa stellar evolution models including 1D hydrodynamics, we find that waves carry ∼107 L⊙ of power from the core to the envelope during core neon/oxygen burning in the final years before core collapse. The waves damp via shocks and radiative diffusion at the base of the hydrogen envelope, which heats up fast enough to launch a pressure wave into the overlying envelope that steepens into a weak shock near the stellar surface, causing a mild stellar outburst and ejecting a small (≲1 M⊙) amount of mass at low speed (≲50 km s−1) roughly one year before the SN. The wave heating inflates the stellar envelope but does not completely unbind it, producing a non-hydrostatic pre-SN envelope density structure different from prior expectations. In our models, wave heating is…

This is the first of a series of papers presenting the results from our survey of 25 Galactic globular clusters with the MUSE integral-field spectrograph.In combination with our dedicated algorithm for source deblending, MUSE provides unique multiplex capabilities in crowded stellar fields and allows us to acquire samples of up to 20 000 stars within the half-light radius of each cluster.The present paper focuses on the analysis of the internal dynamics of 22 out of the 25 clusters, using about 500 000 spectra of 200 000 individual stars.Thanks to the large stellar samples per cluster, we are able to perform a detailed analysis of the central rotation and dispersion fields using both radial profiles and two-dimensional maps.The velocity dispersion profiles we derive show a good general agreement with existing radial velocity studies but typically reach closer to the cluster centres.By comparison with proper motion data, we derive or update the dynamical distance estimates to 14 clusters.Compared to previous dynamical distance estimates for 47 Tuc, our value is in much better…

Variability in the light curves of spotted, rotating stars is often non-sinusoidal and quasi-periodic – spots move on the stellar surface and have finite lifetimes, causing stellar flux variations to slowly shift in phase. A strictly periodic sinusoid therefore cannot accurately model a rotationally modulated stellar light curve. Physical models of stellar surfaces have many drawbacks preventing effective inference, such as highly degenerate or high-dimensional parameter spaces. In this work, we test an appropriate effective model: a Gaussian Process with a quasi-periodic covariance kernel function. This highly flexible model allows sampling of the posterior probability density function of the periodic parameter, marginalizing over the other kernel hyperparameters using a Markov Chain Monte Carlo approach. To test the effectiveness of this method, we infer rotation periods from 333 simulated stellar light curves, demonstrating that the Gaussian process method produces periods that are more accurate than both a sine-fitting periodogram and an autocorrelation function method. We also…

We present a study of the spectral properties of 441 pulsars observed with the Parkes radio telescope near the centre frequencies of 728, 1382 and 3100 MHz. The observations at 728 and 3100 MHz were conducted simultaneously using the dual-band 10-50cm receiver. These high-sensitivity, multi-frequency observations provide a systematic and uniform sample of pulsar flux densities. We combine our measurements with spectral data from the literature in order to derive the spectral properties of these pulsars. Using techniques from robust regression and information theory we classify the observed spectra in an objective, robust and unbiased way into five morphological classes: simple or broken power law, power law with either low or high-frequency cut-off and log-parabolic spectrum. While about $79 \%$ of the pulsars that could be classified have simple power law spectra, we find significant deviations in 73 pulsars, 35 of which have curved spectra, 25 with a spectral break and 10 with a low-frequency turn-over. We identify 11 gigahertz-peaked spectrum (GPS) pulsars, with 3 newly…

Double white dwarf (DWD) binaries are expected to be very common in the Milky Way, but their intrinsic faintness challenges the detection of these systems. Currently, only a few tens of detached DWDs are know. Such systems offer the best chance of extracting the physical properties that would allow us to address a wealth of outstanding questions ranging from the nature of white dwarfs and thermonuclear supernovae, over stellar and binary evolution to mapping the Galaxy. In this paper, we explore the prospects for detections of ultra-compact (with binary separations of a few solar radii or less) detached DWDs in (1) optical radiation with Gaia and the Large Synoptic Survey Telescope (LSST) and (2) gravitational wave radiation with Large Synoptic Survey Telescope (LISA). We show that Gaia, LSST and LISA have the potential to detect, respectively, around a few hundred, a thousand and 25 thousand DWD systems. Moreover, Gaia and LSST data will extend by, respectively, a factor of 2 and 7 the guaranteed sample of LISA verification sources, binaries detectable in electromagnetic and…

Sulphur-bearing species are often used to probe the physical structure of star forming regions of the interstellar medium, but the chemistry of sulphur in these regions is still poorly understood. In dark clouds, sulphur is supposed to be depleted under a form which is still unknown despite numerous observations and chemical modeling studies that have been performed. In order to improve the modeling of sulphur chemistry, we propose an enhancement of the sulphur chemical network using experimental and theoretical literature. We test the effect of the updated network on the outputs of a three phases gas-grain chemical model for dark cloud conditions using different elemental sulphur abundances. More particularly, we focus our study on the main sulphur reservoirs as well as on the agreement between model predictions and the abundances observed in the dark cloud TMC-1 (CP). Our results show that depending on the age of the observed cloud, the reservoir of sulphur could either be atomic sulphur in the gas-phase or HS and H2S in icy grain bulks. We also report the first chemical model…

We show that black hole high-mass X-ray binaries (HMXBs) with O- or B-type donor stars and relatively short orbital periods, of order one week to several months may survive spiral-in, to then form Wolf–Rayet (WR) X-ray binaries with orbital periods of order a day to a few days; while in systems where the compact star is a neutron star, HMXBs with these orbital periods never survive spiral-in. We therefore predict that WR X-ray binaries can only harbour black holes. The reason why black hole HMXBs with these orbital periods may survive spiral-in is: the combination of a radiative envelope of the donor star and a high mass of the compact star. In this case, when the donor begins to overflow its Roche lobe, the systems are able to spiral in slowly with stable Roche lobe overflow, as is shown by the system SS433. In this case, the transferred mass is ejected from the vicinity of the compact star (so-called isotropic re-emission mass-loss mode, or SS433-like mass-loss), leading to gradual spiral-in. If the mass ratio of donor and black hole is ≳3.5, these systems will go into…

We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star formation rates for over 200 000 GAMA galaxies, 170 000 G10-COSMOS galaxies, and 200 000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous data set spanning a broad range in stellar-mass (10 8 -10 12 M ), dust-mass (10 6 -10 9 M ), and star formation rates (0.01-100 M yr -1 ), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous data set using consistent mass and star formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our CSFH we precisely reproduce the SMD with an interstellar medium replenishment factor of…

Motivated by recent interest in the phenomenon of waves transport in massive stars, we examine whether the heat-driven gravity (g) modes excited in slowly pulsating B (SPB) stars can significantly modify the stars’ internal rotation. We develop a formalism for the differential torque exerted by g modes, and implement this formalism using the gyre oscillation code and the mesastar stellar evolution code. Focusing first on a 4.21M⊙ model, we simulate 1 000 yr of stellar evolution under the combined effects of the torque due to a single unstable prograde g mode (with an amplitude chosen on the basis of observational constraints), and diffusive angular momentum transport due to convection, overshooting, and rotational instabilities. We find that the g mode rapidly extracts angular momentum from the surface layers, depositing it deeper in the stellar interior. The angular momentum transport is so efficient that by the end of the simulation, the initially non-rotating surface layers are spun in the retrograde direction to ≈ 30 per cent of the critical rate. However, the additional…

We present high signal-to-noise galaxy-galaxy lensing measurements of the BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS and CS82. We compare this signal with predictions from mock catalogs trained to match observables including the stellar mass function and the projected and two dimensional clustering of CMASS. We show that the clustering of CMASS, together with standard models of the galaxy-halo connection, robustly predicts a lensing signal that is 20-40% larger than observed. Detailed tests show that our results are robust to a variety of systematic effects. Lowering the value of $S_{\rm 8}=σ_{\rm 8} \sqrt{Ω_{\rm m}/0.3}$ compared to Planck2015 reconciles the lensing with clustering. However, given the scale of our measurement ($r<10$ $h^{-1}$ Mpc), other effects may also be at play and need to be taken into consideration. We explore the impact of baryon physics, assembly bias, massive neutrinos, and modifications to general relativity on $ΔΣ$ and show that several of these effects may be non-negligible given the precision of our measurement.…

Magnetic reconnection in relativistic collisionless plasmas can accelerate particles and power high-energy emission in various astrophysical systems. Whereas most previous studies focused on relativistic reconnection in pair plasmas, less attention has been paid to electron–ion plasma reconnection, expected in black hole accretion flows and relativistic jets. We report a comprehensive particle-in-cell numerical investigation of reconnection in an electron–ion plasma, spanning a wide range of ambient ion magnetizations σi, from the semirelativistic regime (ultrarelativistic electrons but non-relativistic ions, 10−3 ≪ σi ≪ 1) to the fully relativistic regime (both species are ultrarelativistic, σi ≫ 1). We investigate how the reconnection rate, electron and ion plasma flows, electric and magnetic field structures, electron/ion energy partitioning, and non-thermal particle acceleration depend on σi. Our key findings are: (1) the reconnection rate is about 0.1 of the Alfvénic rate across all regimes; (2) electrons can form concentrated moderately relativistic outflows even in the…

We fit the rotation curves of isolated dwarf galaxies to directly measure the stellar mass–halo mass relation (M⋆–M200) over the mass range 5 × 105 ≲ M⋆/M⊙ ≲ 108. By accounting for cusp-core transformations due to stellar feedback, we find a monotonic relation with little scatter. Such monotonicity implies that abundance matching should yield a similar M⋆–M200 if the cosmological model is correct. Using the ‘field galaxy’ stellar mass function from the Sloan Digital Sky Survey (SDSS) and the halo mass function from the Λ cold dark matter Bolshoi simulation, we find remarkable agreement between the two. This holds down to M200 ∼ 5 × 109 M⊙, and to M200 ∼ 5 × 108 M⊙ if we assume a power-law extrapolation of the SDSS stellar mass function below M⋆ ∼ 107 M⊙. However, if instead of SDSS we use the stellar mass function of nearby galaxy groups, then the agreement is poor. This occurs because the group stellar mass function is shallower than that of the field below M⋆ ∼ 109 M⊙, recovering the familiar ‘missing satellites’ and ‘too big to fail’ problems. Our result demonstrates that both…

The role of the cosmic web in shaping galaxy properties is investigated in the GAMA spectroscopic survey in the redshift range $0.03 \leq z \leq 0.25$. The stellar mass, $u - r$ dust corrected colour and specific star formation rate (sSFR) of galaxies are analysed as a function of their distances to the 3D cosmic web features, such as nodes, filaments and walls, as reconstructed by DisPerSE. Significant mass and type/colour gradients are found for the whole population, with more massive and/or passive galaxies being located closer to the filament and wall than their less massive and/or star-forming counterparts. Mass segregation persists among the star-forming population alone. The red fraction of galaxies increases when closing in on nodes, and on filaments regardless of the distance to nodes. Similarly, the star-forming population reddens (or lowers its sSFR) at fixed mass when closing in on filament, implying that some quenching takes place. Comparable trends are also found in the state-of-the-art hydrodynamical simulation Horizon-AGN. These results suggest that on top of stellar…

Binary black holes (BBHs) may form both through isolated binary evolution and through dynamical interactions in dense stellar environments. The formation channel leaves an imprint on the alignment between the BH spins and the orbital angular momentum. Gravitational waves (GW) from these systems directly encode information about the spin–orbit misalignment angles, allowing them to be (weakly) constrained. Identifying subpopulations of spinning BBHs will inform us about compact binary formation and evolution. We simulate a mixed population of BBHs with spin–orbit misalignments modelled under a range of assumptions. We then develop a hierarchical analysis and apply it to mock GW observations of these populations. Assuming a population with dimensionless spin magnitudes of χ = 0.7, we show that tens of observations will make it possible to distinguish the presence of subpopulations of coalescing binary black holes based on their spin orientations. With 100 observations, it will be possible to infer the relative fraction of coalescing BBHs with isotropic spin directions (corresponding to…

Mergers of galaxies are thought to cause significant gas inflows to the inner parsecs, which can activate rapid accretion on to supermassive black holes (SMBHs), giving rise to active galactic nuclei (AGN). During a significant fraction of this process, SMBHs are predicted to be enshrouded by gas and dust. Studying 52 galactic nuclei in infrared-selected local luminous and ultraluminous infrared galaxies in different merger stages in the hard X-ray band, where radiation is less affected by absorption, we find that the amount of material around SMBHs increases during the last phases of the merger. We find that the fraction of Compton-thick (CT, N H ≥ 1024 cm− 2) AGN in late-merger galaxies is higher (⁠|$f_{\rm \,CT}=65^{+12}_{-13}{\rm per\, cent}$|⁠) than in local hard X-ray selected AGN (f CT = 27 ± 4 per cent), and that obscuration reaches its maximum when the nuclei of the two merging galaxies are at a projected distance of D12 ≃ 0.4–10.8 kpc (⁠|$f_{\rm \,CT}=77_{-17}^{+13}{\rm per\, cent}$|⁠). We also find that all AGN of our sample in late-merger galaxies have N H > 1023 cm− 2,…

We identify six ultra-diffuse galaxies (UDGs) outside clusters in three nearby isolated groups (0.014 < z < 0.026) using very deep imaging in three different Sloan Digital Sky Survey filters (g, r and i bands) from the IAC Stripe 82 Legacy Project. By comparing with the abundance of UDGs in rich galaxy clusters, we find that the density of UDGs (i.e. the number per unit mass of the host structure where they are located) decreases towards the most massive systems. This is compatible with a scenario where UDGs are formed preferentially outside clusters. In the periphery (D > 250 kpc) of our three groups, we identify a population of potential UDG progenitors (two of them confirmed spectroscopically). These progenitors have similar masses, shapes and sizes but are bluer, g − i ∼ 0.45 [and for this reason brighter, μg(0) < 24 mag arcsec−2] than traditional UDGs (g − i ∼ 0.76). Passive evolution of these progenitors will transform them into regular [i.e. μg(0) > 24 mag arcsec−2] UDGs after ∼6 Gyr. If confirmed, our observations support a scenario where UDGs are old, extended, low surface…

We revisit double neutron star (DNS) formation in the classical binary evolution scenario in light of the recent Laser Interferometer Gravitational-wave Observatory (LIGO)/Virgo DNS detection (GW170817). The observationally estimated Galactic DNS merger rate of |$R_{\rm MW} = 21^{+28}_{-14}$| Myr−1, based on three Galactic DNS systems, fully supports our standard input physics model with RMW = 24 Myr−1. This estimate for the Galaxy translates in a non-trivial way (due to cosmological evolution of progenitor stars in chemically evolving Universe) into a local (z ≈ 0) DNS merger rate density of Rlocal = 48 Gpc−3 yr−1, which is not consistent with the current LIGO/Virgo DNS merger rate estimate (⁠|$1540^{+3200}_{-1220}$| Gpc−3 yr−1). Within our study of the parameter space, we find solutions that allow for DNS merger rates as high as |$R_{\rm local} \approx 600^{+600}_{-300}$| Gpc−3 yr−1 which are thus consistent with the LIGO/Virgo estimate. However, our corresponding BH–BH merger rates for the models with high DNS merger rates exceed the current LIGO/Virgo estimate of local BH–BH…

The ESA Rosetta spacecraft followed comet 67P at a close distance for more than 2 yr. In addition, it deployed the lander Philae on to the surface of the comet. The (surface) composition of the comet is of great interest to understand the origin and evolution of comets. By combining measurements made on the comet itself and in the coma, we probe the nature of this surface material and compare it to remote sensing observations. We compare data from the double focusing mass spectrometer (DFMS) of the ROSINA experiment on ESA's Rosetta mission and previously published data from the two mass spectrometers COSAC (COmetary Sampling And Composition) and Ptolemy on the lander. The mass spectra of all three instruments show very similar patterns of mainly CHO-bearing molecules that sublimate at temperatures of 275 K. The DFMS data also show a great variety of CH-, CHN-, CHS-, CHO2- and CHNO-bearing saturated and unsaturated species. Methyl isocyanate, propanal and glycol aldehyde suggested by the earlier analysis of the measured COSAC spectrum could not be confirmed. The presence of…

We present measurements of the weak gravitational lensing shear power spectrum based on |$450 \deg ^2$| of imaging data from the Kilo Degree Survey. We employ a quadratic estimator in two and three redshift bins and extract band powers of redshift autocorrelation and cross-correlation spectra in the multipole range 76 ≤ ℓ ≤ 1310. The cosmological interpretation of the measured shear power spectra is performed in a Bayesian framework assuming a ΛCDM model with spatially flat geometry, while accounting for small residual uncertainties in the shear calibration and redshift distributions as well as marginalizing over intrinsic alignments, baryon feedback and an excess-noise power model. Moreover, massive neutrinos are included in the modelling. The cosmological main result is expressed in terms of the parameter combination |$S_8 \equiv \sigma _8 \sqrt{\Omega _{\rm m}/0.3}$| yielding S8 = 0.651 ± 0.058 (three z-bins), confirming the recently reported tension in this parameter with constraints from Planck at 3.2σ (three z-bins). We cross-check the results of the three z-bin analysis with…

We present 15 GHz stacked VLBA images of 373 jets associated with active galactic nuclei (AGNs) having at least five observing epochs within a 20 yr time interval 1994–2015 from the Monitoring Of Jets in Active galactic nuclei with VLBA Experiments (MOJAVE) programme and/or its precursor, the 2-cm VLBA Survey. These data are supplemented by 1.4 GHz single-epoch VLBA observations of 135 MOJAVE AGNs to probe larger scale jet structures. The typical jet geometry is found to be close to conical on scales from hundreds to thousands of parsecs, while a number of galaxies show quasi-parabolic streamlines on smaller scales. A true jet geometry in a considerable fraction of AGNs appears only after stacking epochs over several years. The jets with significant radial accelerated motion undergo more active collimation. We have analysed total intensity jet profiles transverse to the local jet ridgeline and derived both apparent and intrinsic opening angles of the flows, with medians of 21|$_{.}^{\circ}$|5 and 1|$_{.}^{\circ}$|3, respectively. The Fermi LAT-detected gamma-ray AGNs in our sample…

We test extensions to the standard cosmological model with weak gravitational lensing tomography using 450 deg2 of imaging data from the Kilo Degree Survey (KiDS). In these extended cosmologies, which include massive neutrinos, non-zero curvature, evolving dark energy, modified gravity and running of the scalar spectral index, we also examine the discordance between KiDS and cosmic microwave background (CMB) measurements from Planck. The discordance between the two data sets is largely unaffected by a more conservative treatment of the lensing systematics and the removal of angular scales most sensitive to non-linear physics. The only extended cosmology that simultaneously alleviates the discordance with Planck and is at least moderately favoured by the data includes evolving dark energy with a time-dependent equation of state (in the form of the w0 − wa parametrization). In this model, the respective |$S_8=\sigma _8\sqrt{\Omega _{\mathrm{m}}/0.3}$| constraints agree at the 1σ level, and there is ‘substantial concordance’ between the KiDS and Planck data sets when accounting for the…

Mergers and the spin of the dark matter halo are factors traditionally believed to determine the morphology of galaxies within a Λ cold dark matter (ΛCDM) cosmology. We study this hypothesis by considering approximately 18 000 central galaxies at z = 0 with stellar masses M* = 109–1012 M⊙ selected from the Illustris cosmological hydrodynamic simulation. The fraction of accreted stars – which measures the importance of massive, recent and dry mergers – increases steeply with galaxy stellar mass, from less than 5 per cent in dwarfs to 80 per cent in the most massive objects, and the impact of mergers on galaxy morphology increases accordingly. For galaxies with M* ≳ 1011 M⊙, mergers have the expected effect: If gas-poor, they promote the formation of spheroidal galaxies, whereas gas-rich mergers favour the formation and survivability of massive discs. This trend, however, breaks at lower masses. For objects with M* ≲ 1011 M⊙, mergers do not seem to play any significant role in determining the morphology, with accreted stellar fractions and mean merger gas fractions that are…

We study the z=0 gas kinematics, morphology, and angular momentum content of isolated galaxies in a suite of cosmological zoom-in simulations from the FIRE project spanning $M_{\star}=10^{6-11}M_{\odot}$. Gas becomes increasingly rotationally supported with increasing galaxy mass. In the lowest-mass galaxies ($M_{\star}<10^{8}M_{\odot}$), gas fails to form a morphological disk and is primarily dispersion and pressure supported. At intermediate masses ($M_{\star}=10^{8-10}M_{\odot}$), galaxies display a wide range of gas kinematics and morphologies, from thin, rotating disks, to irregular spheroids with negligible net rotation. All the high-mass ($M_{\star}=10^{10-11}M_{\odot}$) galaxies form rotationally supported gas disks. Many of the halos whose galaxies fail to form disks harbor high angular momentum gas in their circumgalactic medium. The ratio of the specific angular momentum of gas in the central galaxy to that of the dark-matter halo increases significantly with galaxy mass, from $j_{\rm gas}/j_{\rm DM}\sim0.1$ at $M_{\star}=10^{6-7}M_{\odot}$ to $j_{\rm gas}/j_{\rm…

Author: Elia, D. et al.; Genre: Journal Article; Published in Print: 2017-10; Keywords: catalogues, ISM: clouds, dust, extinction, local interstellar matter, infrared: ISM, submillimetre: ISM; Title: The Hi-GAL compact source catalogue - I. The physical properties of the clumps in the inner Galaxy (-71.0° < ℓ < 67.0°)

We investigate the nucleosynthesis of heavy elements in the winds ejected by accretion discs formed in neutron star mergers. We compute the element formation in disc outflows from hypermassive neutron star (HMNS) remnants of variable lifetime, including the effect of angular momentum transport in the disc evolution. We employ long-term axisymmetric hydrodynamic disc simulations to model the ejecta, and compute r-process nucleosynthesis with tracer particles using a nuclear reaction network containing ∼8000 species. We find that the previously known strong correlation between HMNS lifetime, ejected mass and average electron fraction in the outflow is directly related to the amount of neutrino irradiation on the disc, which dominates mass ejection at early times in the form of a neutrino-driven wind. Production of lanthanides and actinides saturates at short HMNS lifetimes (≲10 ms), with additional ejecta contributing to a blue optical kilonova component for longer-lived HMNSs. We find good agreement between the abundances from the disc outflow alone and the solar r-process…

The early star-forming Universe is still poorly constrained, with the properties of high-redshift stars, the first heating sources, and reionization highly uncertain. This leaves observers planning 21-cm experiments with little theoretical guidance. In this work we explore the possible range of high-redshift parameters including the star formation efficiency and the minimal mass of star-forming halos; the efficiency, spectral energy distribution, and redshift evolution of the first X-ray sources; and the history of reionization. These parameters are only weakly constrained by available observations, mainly the optical depth to the cosmic microwave background. We use realistic semi-numerical simulations to produce the global 21-cm signal over the redshift range $z = 6-40$ for each of 193 different combinations of the astrophysical parameters spanning the allowed range. We show that the expected signal fills a large parameter space, but with a fixed general shape for the global 21-cm curve. Even with our wide selection of models we still find clear correlations between the key…

The processes that led to the formation of the planetary bodies in the Solar system are still not fully understood. Using the results obtained with the comprehensive suite of instruments onboard the European Space Agency's Rosetta mission, we present evidence that comet 67P/Churyumov-Gerasimenko likely formed through the gentle gravitational collapse of a bound clump of mm-sized dust aggregates (‘pebbles’), intermixed with microscopic ice particles. This formation scenario leads to a cometary make-up that is simultaneously compatible with the global porosity, homogeneity, tensile strength, thermal inertia, vertical temperature profiles, sizes and porosities of emitted dust and the steep increase in water-vapour production rate with decreasing heliocentric distance, measured by the instruments onboard the Rosetta spacecraft and the Philae lander. Our findings suggest that the pebbles observed to be abundant in protoplanetary discs around young stars provide the building material for comets and other minor bodies.

We present mock optical images, broad-band and H α fluxes, and D4000 spectral indices for 30 145 galaxies from the eagle hydrodynamical simulation at redshift z = 0.1, modelling dust with the skirt Monte Carlo radiative transfer code. The modelling includes a subgrid prescription for dusty star-forming regions, with both the subgrid obscuration of these regions and the fraction of metals in diffuse interstellar dust calibrated against far-infrared fluxes of local galaxies. The predicted optical colours as a function of stellar mass agree well with observation, with the skirt model showing marked improvement over a simple dust-screen model. The orientation dependence of attenuation is weaker than observed because eagle galaxies are generally puffier than real galaxies, due to the pressure floor imposed on the interstellar medium (ISM). The mock H α luminosity function agrees reasonably well with the data, and we quantify the extent to which dust obscuration affects observed H α fluxes. The distribution of D4000 break values is bimodal, as observed. In the simulation, 20 per cent of…

We perform a tomographic baryon acoustic oscillations (BAOs) analysis using the two-point galaxy correlation function measured from the combined sample of Baryon Oscillation Spectroscopic Survey Data Release 12 (BOSS DR12), which covers the redshift range of 0.2 < z < 0.75. Upon splitting the sample into multiple overlapping redshift slices to extract the redshift information of galaxy clustering, we obtain a measurement of DA(z)/rd and H(z)rd at nine effective redshifts with the full covariance matrix calibrated using MultiDark-Patchy mock catalogues. Using the reconstructed galaxy catalogues, we obtain the precision of 1.3–2.2 per cent for DA(z)/rd and 2.1–6.0 per cent for H(z)rd. To quantify the gain from the tomographic information, we compare the constraints on the cosmological parameters using our nine-bin BAO measurements, the consensus three-bin BAO and redshift space distortion measurements at three effective redshifts in Alam et al., and the non-tomographic (one-bin) BAO measurement at a single effective redshift. Comparing the nine-bin with one-bin constraint result, it…

We present calculations of the wide angle emission of short-duration gamma-ray bursts from compact binary merger progenitors. Such events are expected to be localized by their gravitational wave emission, fairly irrespective of the orientation of the angular momentum vector of the system, along which the gamma-ray burst outflow is expected to propagate. We show that both the prompt and afterglow emission are dim and challenging to detect for observers lying outside the cone within which the relativistic outflow is propagating. If the jet initially propagates through a baryon contaminated region surrounding the merger site, however, a hot cocoon forms around it. The cocoon subsequently expands quasi-isotropically producing its own prompt emission and external shock powered afterglow. We show that the cocoon prompt emission is detectable by Swift BAT and Fermi GBM. We also show that the cocoon afterglow peaks a few hours to a few days after the burst and is detectable for up to a few weeks at all wavelengths. The timing and brightness of the transient are however uncertain due to…

Short gamma-ray bursts are believed to arise from compact binary mergers (either neutron star-neutron star or black hole-neutron star). If so, their jets must penetrate outflows that are ejected during the merger. As a jet crosses the ejecta, it dissipates its energy, producing a hot cocoon that surrounds it. We present here 3D numerical simulations of jet propagation in mergers’ outflows, and we calculate the resulting emission. This emission consists of two components: the cooling emission, the leakage of the thermal energy of the hot cocoon, and the cocoon macronova that arises from the radioactive decay of the cocoon's material. This emission gives a brief (∼1 h) blue, wide angle signal. While the parameters of the outflow and jet are uncertain, for the configurations we have considered, the signal is bright (∼−14 to −15 absolute magnitude) and outshines all other predicted ultraviolet–optical signals. The signal is brighter when the jet breakout time is longer, and its peak brightness does not depend strongly on the highly uncertain opacity. A rapid search for such a signal is…

The inner few hundred parsecs of the Milky Way harbours gas densities, pressures, velocity dispersions, an interstellar radiation field and a cosmic ray ionisation rate orders of magnitude higher than the disc; akin to the environment found in star-forming galaxies at high-redshift. Previous studies have shown that this region is forming stars at a rate per unit mass of dense gas which is at least an order of magnitude lower than in the disc, potentially violating theoretical predictions. We show that all observational star formation rate diagnostics - both direct counting of young stellar objects and integrated light measurements - are in agreement within a factor two, hence the low star formation rate is not the result of the systematic uncertainties that affect any one method. As these methods trace the star formation over different timescales, from $0.1 - 5$ Myr, we conclude that the star formation rate has been constant to within a factor of a few within this time period. We investigate the progression of star formation within gravitationally bound clouds on $\sim$ parsec…

We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102.Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst.No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities.None of the FRBs were seen to repeat.Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time.FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = 2596.1±0.3 pc cm -3 ) detected to date.Three of the FRBs have high dispersion measures (DM >1500 pc cm -3 ), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium.The slope of the Parkes FRB source counts distribution with fluences > 2 Jy…

Reionization in the early Universe is likely driven by dwarf galaxies. Using cosmological radiation-hydrodynamic simulations, we study star formation and the escape of Lyman continuum (LyC) photons from mini-haloes with |${M_{\rm halo}}\lesssim 10^8\,{\rm M_{\odot }}$|⁠. Our simulations include a new thermo-turbulent star formation model, non-equilibrium chemistry and relevant stellar feedback processes (photoionization by young massive stars, radiation pressure and mechanical supernova explosions). We find that feedback reduces star formation very efficiently in mini-haloes, resulting in the stellar mass consistent with the slope and normalization reported in Kimm & Cen and the empirical stellar mass-to-halo mass relation derived in the local Universe. Because star formation is stochastic and dominated by a few gas clumps, the escape fraction in mini-haloes is generally determined by radiation feedback (heating due to photoionization), rather than supernova explosions. We also find that the photon number-weighted mean escape fraction in mini-haloes is higher (∼20–40 per cent) than…

Galaxy-scale strong gravitational lensing can not only provide a valuable probe of the dark matter distribution of massive galaxies, but also provide valuable cosmological constraints, either by studying the population of strong lenses or by measuring time delays in lensed quasars. Due to the rarity of galaxy-scale strongly lensed systems, fast and reliable automated lens finding methods will be essential in the era of large surveys such as Large Synoptic Survey Telescope, Euclid and Wide-Field Infrared Survey Telescope. To tackle this challenge, we introduce CMU DeepLens, a new fully automated galaxy–galaxy lens finding method based on deep learning. This supervised machine learning approach does not require any tuning after the training step which only requires realistic image simulations of strongly lensed systems. We train and validate our model on a set of 20 000 LSST-like mock observations including a range of lensed systems of various sizes and signal-to-noise ratios (S/N). We find on our simulated data set that for a rejection rate of non-lenses of 99 per cent, a…

We release the AllWISE counterparts and Gaia matches to 106 573 and 17 665 X-ray sources detected in the ROSAT 2RXS and XMMSL2 surveys with |b| > 15°. These are the brightest X-ray sources in the sky, but their position uncertainties and the sparse multi-wavelength coverage until now rendered the identification of their counterparts a demanding task with uncertain results. New all-sky multi-wavelength surveys of sufficient depth, like AllWISE and Gaia, and a new Bayesian statistics based algorithm, Nway, allow us, for the first time, to provide reliable counterpart associations. Nway extends previous distance and sky density based association methods and, using one or more priors (e.g. colours, magnitudes), weights the probability that sources from two or more catalogues are simultaneously associated on the basis of their observable characteristics. Here, counterparts have been determined using a Wide-field Infrared Survey Explorer (WISE) colour-magnitude prior. A reference sample of 4524 XMM/Chandra and Swift X-ray sources demonstrates a reliability of ∼94.7 per cent (2RXS) and…

The volume of data that will be produced by new-generation surveys requires automatic classification methods to select and analyze sources. Indeed, this is the case for the search for strong gravitational lenses, where the population of the detectable lensed sources is only a very small fraction of the full source population. We apply for the first time a morphological classification method based on a Convolutional Neural Network (CNN) for recognizing strong gravitational lenses in $255$ square degrees of the Kilo Degree Survey (KiDS), one of the current-generation optical wide surveys. The CNN is currently optimized to recognize lenses with Einstein radii $\gtrsim 1.4$ arcsec, about twice the $r$-band seeing in KiDS. In a sample of $21789$ colour-magnitude selected Luminous Red Galaxies (LRG), of which three are known lenses, the CNN retrieves 761 strong-lens candidates and correctly classifies two out of three of the known lenses. The misclassified lens has an Einstein radius below the range on which the algorithm is trained. We down-select the most reliable 56 candidates by a…

The cosmic merger rate density of black hole binaries (BHBs) can give us an essential clue to constraining the formation channels of BHBs, in light of current and forthcoming gravitational wave detections. Following a Monte Carlo approach, we couple new population-synthesis models of BHBs with the Illustris cosmological simulation, to study the cosmic history of BHB mergers. We explore six population-synthesis models, varying the prescriptions for supernovae, common envelope and natal kicks. In most considered models, the cosmic BHB merger rate follows the same trend as the cosmic star formation rate. The normalization of the cosmic BHB merger rate strongly depends on the treatment of common envelope and on the distribution of natal kicks. We find that most BHBs merging within LIGO's instrumental horizon come from relatively metal-poor progenitors (<0.2 Z⊙). The total masses of merging BHBs span a large range of values, from ∼6 to ∼82 M⊙. In our fiducial model, merging BHBs consistent with GW150914, GW151226 and GW170104 represent ∼6, 3 and 12 per cent of all BHBs merging within the…

Nearby dwarf galaxies provide a unique laboratory in which to test stellar population models below Z⊙/2. Such tests are particularly important for interpreting the surprising high-ionization ultraviolet (UV) line emission detected at z > 6 in recent years. We present HST/COS UV spectra of 10 nearby metal-poor star-forming galaxies selected to show He ii emission in SDSS optical spectra. The targets span nearly a dex in gas-phase oxygen abundance (7.8 < 12 + log O/H < 8.5) and present uniformly large specific star formation rates (sSFR ∼102 Gyr−1). The UV spectra confirm that metal-poor stellar populations can power extreme nebular emission in high-ionization UV lines, reaching C iii] equivalent widths comparable to those seen in systems at z ∼ 6–7. Our data reveal a marked transition in UV spectral properties with decreasing metallicity, with systems below 12 + log O/H ≲ 8.0 (Z/Z⊙ ≲ 1/5) presenting minimal stellar wind features and prominent nebular emission in He ii and C iv. This is consistent with nearly an order of magnitude increase in ionizing photon production beyond the…

The measurement of the structure of stellar populations in the Milky Way disc places fundamental constraints on models of galaxy formation and evolution. Previously, the disc's structure has been studied in terms of populations defined geometrically and/or chemically, but a decomposition based on stellar ages provides a more direct connection to the history of the disc, and stronger constraint on theory. Here, we use positions, abundances and ages for 31 244 red giant branch stars from the Sloan Digital Sky Survey (SDSS)-APOGEE survey, spanning 3 < Rgc < 15 kpc, to dissect the disc into mono-age and mono-[Fe/H] populations at low and high |$\mathrm{[ \alpha \mathrm{/Fe]}}$|⁠. For each population, with Δage < 2 Gyr and Δ[Fe/H] < 0.1 dex, we measure the structure and surface-mass density contribution. We find that low |$\mathrm{[ \alpha \mathrm{/Fe]}}$| mono-age populations are fit well by a broken exponential, which increases to a peak radius and decreases thereafter. We show that this profile becomes broader with age, interpreted here as a new signal of disc heating and radial…

Abstract The magnetic fields observed in the Milky Way and nearby galaxies appear to be in equipartition with the turbulent, thermal and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky Way like galaxies, carried out with a moving-mesh magnetohydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly 100 Myr in the centre of haloes until saturation occurs around z = 2–3, when the magnetic energy density reaches about 10 per cent of the turbulent energy density with a typical strength of $10\text{--}50\,\rm {\mu G}$. In the galactic centres, the ratio between magnetic and…

We present a new spectral fitting code, firefly, for deriving the stellar population properties of stellar systems. firefly is a chi-squared minimization fitting code that fits combinations of single-burst stellar population models to spectroscopic data, following an iterative best-fitting process controlled by the Bayesian information criterion. No priors are applied, rather all solutions within a statistical cut are retained with their weight. Moreover, no additive or multiplicative polynomials are employed to adjust the spectral shape. This fitting freedom is envisaged in order to map out the effect of intrinsic spectral energy distribution degeneracies, such as age, metallicity, dust reddening on galaxy properties, and to quantify the effect of varying input model components on such properties. Dust attenuation is included using a new procedure, which was tested on Integral Field Spectroscopic data in a previous paper. The fitting method is extensively tested with a comprehensive suite of mock galaxies, real galaxies from the Sloan Digital Sky Survey and Milky Way globular…

The variations of galaxy stellar masses and colour-types with the distance to projected cosmic filaments are quantified using the precise photometric redshifts of the COSMOS2015 catalogue extracted from COSMOS field (2 deg 2 ).Realistic mock catalogues are also extracted from the lightcone of the cosmological hydrodynamical simulation HORIZON-AGN.They show that the photometric redshift accuracy of the observed catalogue (σ z < 0.015 at M * > 10 10 M and z < 0.9) is sufficient to provide 2D filaments that closely match their projected 3D counterparts.Transverse stellar mass gradients are measured in projected slices of thickness 75 Mpc between 0.5 < z < 0.9, showing that the most massive galaxies are statistically closer to their neighbouring filament.At fixed stellar mass, passive galaxies are also found closer to their filament while active star-forming galaxies statistically lie further away.The contributions of nodes and local density are removed from these gradients to highlight the specific role played by the geometry of the filaments.We find that the measured signal does…

We present gravitational wave (GW) signal predictions from four 3D multigroup neutrino hydrodynamics simulations of core-collapse supernovae of progenitors with 11.2, 20 and 27 M⊙. GW emission in the pre-explosion phase strongly depends on whether the post-shock flow is dominated by the standing accretion shock instability (SASI) or convection and differs considerably from 2D models. SASI activity produces a strong signal component below 250 Hz through asymmetric mass motions in the gain layer and a non-resonant coupling to the proto-neutron star (PNS). Both convection- and SASI-dominated models show GW emission above 250 Hz, but with considerably lower amplitudes than in 2D. This is due to a different excitation mechanism for high-frequency l = 2 motions in the PNS surface, which are predominantly excited by PNS convection in 3D. Resonant excitation of high-frequency surface g modes in 3D by mass motions in the gain layer is suppressed compared to 2D because of smaller downflow velocities and a lack of high-frequency variability in the downflows. In the exploding 20 M⊙ model, shock…

We use deep Chandra X-ray imaging to measure the distribution of specific black hole accretion rates (LX relative to the stellar mass of the galaxy) and thus trace active galactic nucleus (AGN) activity within star-forming and quiescent galaxies, as a function of stellar mass (from 108.5 to 1011.5 M⊙) and redshift (to z ∼ 4). We adopt near-infrared-selected samples of galaxies from the CANDELS and UltraVISTA surveys, extract X-ray data for every galaxy, and use a flexible Bayesian method to combine these data and to measure the probability distribution function of specific black hole accretion rates, λsBHAR. We identify a broad distribution of λsBHAR in both star-forming and quiescent galaxies – likely reflecting the stochastic nature of AGN fuelling – with a roughly power-law shape that rises towards lower λsBHAR, a steep cut-off at λsBHAR ≳ 0.1–1 (in Eddington equivalent units), and a turnover or flattening at |$\lambda _\mathrm{sBHAR} \lesssim 10^{-3}\hbox{ {to} }10^{-2}$|⁠. We find that the probability of a star-forming galaxy hosting a moderate λsBHAR AGN depends on stellar…

We report on five compact, extremely young (<10 Myr) and blue (βUV < −2.5, Fλ = λβ) objects observed with VLT/Multi Unit Spectroscopic Explorer at redshifts 3.1169 and 3.235, in addition to three objects at z = 6.145. These sources are strongly magnified (3–40 times) by the Hubble Frontier Field galaxy clusters MACS J0416 and AS1063. Their delensed half-light radii (Re) are between 16 and 140 pc, the stellar masses are ≃1–20 × 106 M⊙, the magnitudes are mUV = 28.8–31.4 (−17 < MUV < −15) and specific star formation rates can be as large as ∼800 Gyr−1. Remarkably, the inferred physical properties of two objects are similar to those expected in some globular cluster formation scenarios, representing the best candidate proto-GCs discovered so far. Rest-frame optical high-dispersion spectroscopy of one of them at z = 3.1169 yields a velocity dispersion σv ≃ 20 km s−1, implying a dynamical mass dominated by the stellar mass. Another object at z = 6.145, with delensed MUV ≃ −15.3 (mUV ≃ 31.4), shows a stellar mass and a star formation rate surface density consistent with the values…

We describe ipole, a new public ray-tracing code for covariant, polarized radiative transport. The code extends the ibothros scheme for covariant, unpolarized transport using two representations of the polarized radiation field: In the coordinate frame, it parallel transports the coherency tensor; in the frame of the plasma it evolves the Stokes parameters under emission, absorption, and Faraday conversion. The transport step is implemented to be as spacetime- and coordinate- independent as possible. The emission, absorption, and Faraday conversion step is implemented using an analytic solution to the polarized transport equation with constant coefficients. As a result, ipole is stable, efficient, and produces a physically reasonable solution even for a step with high optical depth and Faraday depth. We show that the code matches analytic results in flat space, and that it produces results that converge to those produced by Dexter's grtrans polarized transport code on a complicated model problem. We expect ipole will mainly find applications in modelling Event Horizon Telescope…

Atmospheric chemical abundances of giant planets lead to important constraints on planetary formation and migration. Studies have shown that giant planets that migrate through the protoplanetary disc can accrete substantial amounts of oxygen-rich planetesimals, leading to supersolar metallicities in the envelope and solar or subsolar C/O ratios. Pebble accretion has been demonstrated to play an important role in core accretion and to have growth rates that are consistent with planetary migration. The high pebble accretion rates allow planetary cores to start their growth beyond 10 au and subsequently migrate to cold (≳1 au), warm (∼0.1–1 au) or hot (≲0.1 au) orbits. In this work we investigate how the formation of giant planets via pebble accretion influences their atmospheric chemical compositions. We find that under the standard pebble accretion scenario, where the core is isolated from the envelope, the resulting metallicities (O/H and C/H ratios) are subsolar, while the C/O ratios are supersolar. Planets that migrate through the disc to become hot Jupiters accrete substantial…

Stellar-remnant black holes (BH) in dense stellar clusters have always drawn attention due to their potential in a number of phenomena, especially the dynamical formation of binary black holes (BBH), which potentially coalesce via gravitational-wave radiation. This study presents a preliminary set of evolutionary models of compact stellar clusters with initial masses ranging over 1.0 × 104–5.0 × 104 M⊙, and half-mass radius of 2 or 1 pc, which is typical for young massive and starburst clusters. They have metallicities between 0.05 Z⊙ and Z⊙. Including contemporary schemes for stellar wind and remnant formation, such model clusters are evolved, for the first time, using the state-of-the-art direct N-body evolution program nbody7, until their dissolution or at least for 10 Gyr. That way, a self-regulatory behaviour in the effects of dynamical interactions among the BHs is demonstrated. In contrast to earlier studies, the BBH coalescences obtained in these models show a prominence in triple-mediated coalescences while being bound to the clusters, compared to those occurring among the…

I use distance- and mass-limited subsamples of the Spitzer Survey of Stellar Structure in Galaxies (S4G) to investigate how the presence of bars in spiral galaxies depends on mass, colour, and gas content and whether large, Sloan Digital Sky Survey (SDSS)-based investigations of bar frequencies agree with local data. Bar frequency reaches a maximum of fbar ≈ 0.70 at M⋆ ∼ 109.7M⊙, declining to both lower and higher masses. It is roughly constant over a wide range of colours (g − r ≈ 0.1–0.8) and atomic gas fractions (⁠|$\log (M_{\mathrm{H} \small {I}}/ M_{\star }) \approx -2.5$| to 1). Bars are thus as common in blue, gas-rich galaxies are they are in red, gas-poor galaxies. This is in sharp contrast to many SDSS-based studies of z ∼ 0.01–0.1 galaxies, which report fbar increasing strongly to higher masses (from M⋆ ∼ 1010 to 1011M⊙), redder colours, and lower gas fractions. The contradiction can be explained if SDSS-based studies preferentially miss bars in, and underestimate the bar fraction for, lower mass (bluer, gas-rich) galaxies due to poor spatial resolution and the…

The tight correlation between total galaxy stellar mass and star formation rate (SFR) has become known as the star-forming main sequence. Using ∼487 000 spaxels from galaxies observed as part of the Sloan Digital Sky Survey Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we confirm previous results that a correlation also exists between the surface densities of star formation (ΣSFR) and stellar mass (Σ⋆) on kpc scales, representing a ‘resolved’ main sequence. Using a new metric (ΔΣSFR), which measures the relative enhancement or deficit of star formation on a spaxel-by-spaxel basis relative to the resolved main sequence, we investigate the SFR profiles of 864 galaxies as a function of their position relative to the global star-forming main sequence (ΔSFR). For galaxies above the global main sequence (positive ΔSFR) ΔΣSFR is elevated throughout the galaxy, but the greatest enhancement in star formation occurs at small radii (<3 kpc, or 0.5Re). Moreover, galaxies that are at least a factor of 3 above the main sequence show diluted gas phase metallicities out to…