scispace - formally typeset
Search or ask a question

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


Journal ArticleDOI
TL;DR: In this paper, the authors used the suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies.
Abstract: Predicting structural properties of dark matter haloes is one of the fundamental goals of modern cosmology. We use the suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. We find that in order to understand the structure of dark matter haloes and to make 1–2 per cent accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the ratio of the virial radius to the core radius in the Navarro–Frenk–White (NFW) profile. For massive haloes, the average density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter α in the Einasto approximation. We show that haloes progress through three stages of evolution. They start as rare density peaks and experience fast and nearly radial infall that brings mass closer to the centre, producing a highly concentrated halo. Here, the halo concentration increases with increasing halo mass and the concentration is defined by the α parameter with a nearly constant core radius. Later haloes slide into the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. At this stage, the concentration does not depend on halo mass. Once the rate of accretion and merging slows down, haloes move into the domain of declining concentration–mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. Accurate analytical fits are provided.

808 citations


Journal ArticleDOI
T. M. C. Abbott, F. B. Abdalla1, Jelena Aleksić2, S. Allam3  +153 moreInstitutions (43)
TL;DR: In this paper, the authors presented the results of the Dark Energy Survey (DES) 2013, 2014, 2015, 2016, 2017, 2018, 2019 and 2019 at the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign.
Abstract: US Department of Energy; US National Science Foundation; Ministry of Science and Education of Spain; Science and Technology Facilities Council of the United Kingdom; Higher Education Funding Council for England; National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign; Kavli Institute of Cosmological Physics at the University of Chicago; Center for Cosmology and Astro-Particle Physics at the Ohio State University; Mitchell Institute for Fundamental Physics and Astronomy at Texas AM University; Financiadora de Estudos e Projetos; Fundacao Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro; Conselho Nacional de Desenvolvimento Cientifico e Tecnologico and the Ministerio da Ciencia; Tecnologia e Inovacao; Deutsche Forschungsgemeinschaft; Collaborating Institutions in the Dark Energy Survey; National Science Foundation [AST-1138766]; University of California at Santa Cruz; University of Cambridge, Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas-Madrid; University of Chicago, University College London; DES-Brazil Consortium; University of Edinburgh; Eidgenossische Technische Hochschule (ETH) Zurich, Fermi National Accelerator Laboratory; University of Illinois at Urbana-Champaign; Institut de Ciencies de l'Espai (IEEC/CSIC); Institut de Fisica d'Altes Energies, Lawrence Berkeley National Laboratory; Ludwig-Maximilians Universitat Munchen; European Research Council [FP7/291329]; MINECO [AYA2012-39559, ESP2013-48274, FPA2013-47986]; Centro de Excelencia Severo Ochoa [SEV-2012-0234]; European Research Council under the European Union [240672, 291329, 306478]

789 citations


Journal ArticleDOI
TL;DR: In this article, a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members is presented. And the results from APOSTLE (A Project Of Simulating The Local Environment) are compared to the observed abundance of LG galaxies.
Abstract: The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" ({\Lambda}CDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in {\Lambda}CDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

501 citations


Journal ArticleDOI
TL;DR: In this article, the authors explore a newly proposed channel to create binary black holes of stellar origin, where mixing induced by rotation and tides transports the products of hydrogen burning throughout the stellar envelopes.
Abstract: We explore a newly proposed channel to create binary black holes of stellar origin. This scenario applies to massive, tight binaries where mixing induced by rotation and tides transports the products of hydrogen burning throughout the stellar envelopes. This slowly enriches the entire star with helium, preventing the build-up of an internal chemical gradient. The stars remain compact as they evolve nearly chemically homogeneously, eventually forming two black holes, which we estimate typically merge 4–11 Gyr after formation. Like other proposed channels, this evolutionary pathway suffers from significant theoretical uncertainties, but could be constrained in the near future by data from advanced ground-based gravitational-wave detectors. We perform Monte Carlo simulations of the expected merger rate over cosmic time to explore the implications and uncertainties. Our default model for this channel yields a local binary black hole merger rate of about 10 Gpc−3 yr−1 at redshift z = 0, peaking at twice this rate at z = 0.5. This means that this channel is competitive, in terms of expected rates, with the conventional formation scenarios that involve a common-envelope phase during isolated binary evolution or dynamical interaction in a dense cluster. The events from this channel may be distinguished by the preference for nearly equal-mass components and high masses, with typical total masses between 50 and 110 M⊙. Unlike the conventional isolated binary evolution scenario that involves shrinkage of the orbit during a common-envelope phase, short time delays are unlikely for this channel, implying that we do not expect mergers at high redshift.

496 citations


Journal ArticleDOI
TL;DR: In this paper, the authors report on timing, flux density, and polarimetric observations of the transient magnetar and 5.54 s radio pulsar XTE J1810-197 using the GBT, Nancay, and Parkes radio telescopes beginning in early 2006, until its sudden disappearance as a radio source in late 2008.
Abstract: We report on timing, flux density, and polarimetric observations of the transient magnetar and 5.54 s radio pulsar XTE J1810-197 using the GBT, Nancay, and Parkes radio telescopes beginning in early 2006, until its sudden disappearance as a radio source in late 2008. Repeated observations through 2016 have not detected radio pulsations again. The torque on the neutron star, as inferred from its rotation frequency derivative f-dot, decreased in an unsteady manner by a factor of 3 in the first year of radio monitoring. In contrast, during its final year as a detectable radio source, the torque decreased steadily by only 9%. The period-averaged flux density, after decreasing by a factor of 20 during the first 10 months of radio monitoring, remained steady in the next 22 months, at an average of 0.7+/-0.3 mJy at 1.4 GHz, while still showing day-to-day fluctuations by factors of a few. There is evidence that during this last phase of radio activity the magnetar had a steep radio spectrum, in contrast to earlier behavior. There was no secular decrease that presaged its radio demise. During this time the pulse profile continued to display large variations, and polarimetry indicates that the magnetic geometry remained consistent with that of earlier times. We supplement these results with X-ray timing of the pulsar from its outburst in 2003 up to 2014. For the first 4 years, XTE J1810-197 experienced non-monotonic excursions in f-dot by at least a factor of 8. But since 2007, its f-dot has remained relatively stable near its minimum observed value. The only apparent event in the X-ray record that is possibly contemporaneous with the radio shut-down is a decrease of ~20% in the hot-spot flux in 2008-2009, to a stable, minimum value. However, the permanence of the high-amplitude, thermal X-ray pulse, even after the radio demise, implies continuing magnetar activity.

429 citations


Journal ArticleDOI
Joris P. W. Verbiest1, Joris P. W. Verbiest2, L. Lentati, George Hobbs3, R. van Haasteren4, Paul Demorest5, Gemma H. Janssen, J. B. Wang6, Gregory Desvignes1, R. N. Caballero1, Michael Keith, D. J. Champion1, Zaven Arzoumanian7, Stanislav Babak8, C. G. Bassa9, N. D. R. Bhat10, A. Brazier11, P. Brem8, M. Burgay12, Sarah Burke-Spolaor5, S. J. Chamberlin13, Sourav Chatterjee11, B. Christy14, Ismaël Cognard15, Ismaël Cognard16, James M. Cordes11, Shi Dai3, Shi Dai17, Timothy Dolch18, Timothy Dolch11, Justin A. Ellis4, Robert D. Ferdman, Emmanuel Fonseca19, Jonathan R. Gair20, N. Garver-Daniels21, Peter A. Gentile21, Marjorie Gonzalez22, E. Graikou1, Lucas Guillemot15, Lucas Guillemot16, Jason W. T. Hessels9, Jason W. T. Hessels23, Glenn Jones24, Ramesh Karuppusamy, Matthew Kerr3, Michael Kramer, Michael T. Lam11, Paul D. Lasky25, A. Lassus1, P. Lazarus1, T. J. W. Lazio4, Kejia Lee17, Lina Levin26, Lina Levin21, Kang Liu1, R. S. Lynch5, Andrew Lyne, J. W. McKee26, Maura McLaughlin21, Sean T. McWilliams21, D. R. Madison5, Richard N. Manchester3, Chiara M. F. Mingarelli1, Chiara M. F. Mingarelli4, David J. Nice27, Stefan Oslowski2, Stefan Oslowski1, Nipuni Palliyaguru28, Timothy T. Pennucci29, Benetge Perera, Delphine Perrodin12, A. Possenti12, Antoine Petiteau30, Scott M. Ransom5, Daniel J. Reardon3, Daniel J. Reardon25, Pablo Rosado31, S. A. Sanidas23, Alberto Sesana32, G. Shaifullah1, G. Shaifullah2, Ryan Shannon10, Ryan Shannon3, X. Siemens33, Joseph Simon33, R. Smits, Renée Spiewak33, Ingrid H. Stairs19, Benjamin Stappers, Daniel R. Stinebring34, Kevin Stovall35, J. K. Swiggum21, Stephen Taylor4, Gilles Theureau16, Gilles Theureau30, Gilles Theureau15, Caterina Tiburzi2, Caterina Tiburzi1, L. Toomey3, Michele Vallisneri4, W. van Straten31, Alberto Vecchio32, Yue-Fei Wang36, Linqing Wen37, X. P. You38, Weiwei Zhu1, Xing-Jiang Zhu37 
TL;DR: In this article, the first joint analysis of the data from the three regional pulsar timing arrays (IPTA) is presented, i.e. of the first IPTA data set, and the approach presently followed for its combination and suggest improvements for future PTA research.
Abstract: The highly stable spin of neutron stars can be exploited for a variety of (astro)physical investigations. In particular, arrays of pulsars with rotational periods of the order of milliseconds can be used to detect correlated signals such as those caused by gravitational waves. Three such 'pulsar timing arrays' (PTAs) have been set up around the world over the past decades and collectively form the 'International' PTA (IPTA). In this paper, we describe the first joint analysis of the data from the three regional PTAs, i.e. of the first IPTA data set. We describe the available PTA data, the approach presently followed for its combination and suggest improvements for future PTA research. Particular attention is paid to subtle details (such as underestimation of measurement uncertainty and long-period noise) that have often been ignored but which become important in this unprecedentedly large and inhomogeneous data set. We identify and describe in detail several factors that complicate IPTA research and provide recommendations for future pulsar timing efforts. The first IPTA data release presented here (and available on-line) is used to demonstrate the IPTA's potential of improving upon gravitational-wave limits

412 citations


Journal ArticleDOI
TL;DR: In this article, the authors used the Illustris simulation to study the relative contributions of in situ star formation and stellar accretion to the build-up of galaxies over an unprecedentedly wide range of masses.
Abstract: We use the Illustris simulation to study the relative contributions of in situ star formation and stellar accretion to the build-up of galaxies over an unprecedentedly wide range of masses (M_* = 10^9-10^(12) M_⊙), galaxy types, environments, and assembly histories. We find that the ‘two-phase’ picture of galaxy formation predicted by some models is a good approximation only for the most massive galaxies in our simulation – namely, the stellar mass growth of galaxies below a few times 10^(11) M_⊙ is dominated by in situ star formation at all redshifts. The fraction of the total stellar mass of galaxies at z = 0 contributed by accreted stars shows a strong dependence on galaxy stellar mass, ranging from about 10 per cent for Milky Way-sized galaxies to over 80 per cent for M_* ≈ 10^(12) M_⊙ objects, yet with a large galaxy-to-galaxy variation. At a fixed stellar mass, elliptical galaxies and those formed at the centres of younger haloes exhibit larger fractions of ex situ stars than disc-like galaxies and those formed in older haloes. On average, ∼50 per cent of the ex situ stellar mass comes from major mergers (stellar mass ratio μ > 1/4), ∼20 per cent from minor mergers (1/10 < μ < 1/4), ∼20 per cent from very minor mergers (μ < 1/10), and ∼10 per cent from stars that were stripped from surviving galaxies (e.g. flybys or ongoing mergers). These components are spatially segregated, with in situ stars dominating the innermost regions of galaxies, and ex situ stars being deposited at larger galactocentric distances in order of decreasing merger mass ratio.

384 citations


Journal ArticleDOI
TL;DR: The Baryon Oscillation Spectroscopic Survey (BOSS) as discussed by the authors provides the largest survey of galaxy redshifts available to date, in terms of both the number of galaxies measured by a single survey, and the effective cosmological volume covered.
Abstract: The Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey (SDSS) III project, has provided the largest survey of galaxy redshifts available to date, in terms of both the number of galaxy redshifts measured by a single survey, and the effective cosmological volume covered. Key to analysing the clustering of these data to provide cosmological measurements is understanding the detailed properties of this sample. Potential issues include variations in the target catalogue caused by changes either in the targeting algorithm or properties of the data used, the pattern of spectroscopic observations, the spatial distribution of targets for which redshifts were not obtained, and variations in the target sky density due to observational systematics. We document here the target selection algorithms used to create the galaxy samples that comprise BOSS. We also present the algorithms used to create large-scale structure catalogues for the final Data Release (DR12) samples and the associated random catalogues that quantify the survey mask. The algorithms are an evolution of those used by the BOSS team to construct catalogues from earlier data, and have been designed to accurately quantify the galaxy sample. The code used, designated MKSAMPLE, is released with this paper.

373 citations


Journal ArticleDOI
TL;DR: In this article, a general-relativistic simulation of binary neutron star mergers with a temperature and composition dependent nuclear equation of state is presented, where the outflow is composed of a combination of tidally and shock-driven ejecta, mostly distributed over a broad ∼60∘ angle from the orbital plane and, to a lesser extent, by thermally driven winds at high latitudes.
Abstract: We present fully general-relativistic simulations of binary neutron star mergers with a temperature and composition dependent nuclear equation of state. We study the dynamical mass ejection from both quasi-circular and dynamical-capture eccentric mergers. We systematically vary the level of our treatment of the microphysics to isolate the effects of neutrino cooling and heating and we compute the nucleosynthetic yields of the ejecta. We find that eccentric binaries can eject significantly more material than quasi-circular binaries and generate bright infrared and radio emission. In all our simulations the outflow is composed of a combination of tidally- and shock-driven ejecta, mostly distributed over a broad ∼60∘ angle from the orbital plane, and, to a lesser extent, by thermally driven winds at high latitudes. Ejecta from eccentric mergers are typically more neutron rich than those of quasi-circular mergers. We find neutrino cooling and heating to affect, quantitatively and qualitatively, composition, morphology, and total mass of the outflows. This is also reflected in the infrared and radio signatures of the binary. The final nucleosynthetic yields of the ejecta are robust and insensitive to input physics or merger type in the regions of the second and third r-process peaks. The yields for elements on the first peak vary between our simulations, but none of our models is able to explain the Solar abundances of first-peak elements without invoking additional first-peak contributions from either neutrino and viscously-driven winds operating on longer timescales after the mergers, or from core-collapse supernovae.

372 citations


Journal ArticleDOI
TL;DR: In this paper, the authors explore the predictions for detectable gravitational-wave signals from merging binary black holes formed through chemically homogeneous evolution in massive short-period stellar binaries and conclude that ∼500 events per year could be detected with advanced ground-based detectors operating at full sensitivity.
Abstract: We explore the predictions for detectable gravitational-wave signals from merging binary black holes formed through chemically homogeneous evolution in massive short-period stellar binaries. We find that ∼500 events per year could be detected with advanced ground-based detectors operating at full sensitivity. We analyse the distribution of detectable events, and conclude that there is a very strong preference for detecting events with nearly equal components (mass ratio >0.66 at 90 per cent confidence in our default model) and high masses (total source-frame mass between 57 and 103 M⊙ at 90 per cent confidence). We consider multiple alternative variations to analyse the sensitivity to uncertainties in the evolutionary physics and cosmological parameters, and conclude that while the rates are sensitive to assumed variations, the mass distributions are robust predictions. Finally, we consider the recently reported results of the analysis of the first 16 double-coincident days of the O1 LIGO (Laser Interferometer Gravitational-wave Observatory) observing run, and find that this formation channel is fully consistent with the inferred parameters of the GW150914 binary black hole detection and the inferred merger rate.

371 citations


Journal ArticleDOI
TL;DR: In this article, the authors consider radio bursts that originate from extragalactic neutron stars (NSs) by addressing three questions about source distances: What are the physical limitations on coherent radiation at GHz frequencies? Do they permit detection at cosmological distances? How many bursts per NS are needed to produce the inferred burst rate 10 3 -10 4 sky 1 day 1?
Abstract: We consider radio bursts that originate from extragalactic neutron stars (NSs) by addressing three questions about source distances. What are the physical limitations on coherent radiation at GHz frequencies? Do they permit detection at cosmological distances? How many bursts per NS are needed to produce the inferred burst rate 10 3 -10 4 sky 1 day 1 ? The burst rate is comparable to the NS formation rate in a Hubble volume, requiring only one per NS if they are bright enough. However, radiation physics causes us to favor a closer population. More bursts per NS are then required but repeats in 10 to 100 yr could still be negligible. Bursts are modeled as sub-ns, coherent shot pulses superposed incoherently to produce msduration 1 Jy amplitudes; each shot-pulse can be much weaker than the burst amplitude, placing less restrictive requirements on the emission process. Nonetheless, single shot pulses are similar to the extreme, unresolved (< 0:4 ns) MJy shot pulse seen from the Crab pulsar, which is consistent with coherent curvature radiation emitted near the light cylinder by an almost neutral clump with net charge 10 21 e and total energy & 10 23 ergs. Bursts from Gpc distances require incoherent superposition of 10 12 d 2 shot pulses or a total energy & 10 35 d 2 erg. The energy reservoir near the light cylinder limits the detection distance to . few 100 Mpc for a fluence 1 Jy ms unless conditions are more extreme than for the Crab pulsar. Similarly, extreme single pulses from ordinary pulsars and magnetars could be detectable from throughout the Local Group and perhaps farther. Contributions to dispersion measures from galaxy clusters will be significant for some of the bursts. We discuss tests for the signatures of bursts associated with extragalactic NSs.

Journal ArticleDOI
TL;DR: In this article, the authors presented results of a spectroscopic survey for new K and M-type members of Scorpius-Centaurus (Sco-Cen), the nearest OB Association (~100-200 pc) using an X-ray, proper motion and color-magnitude selected sample.
Abstract: We present results of a spectroscopic survey for new K- and M-type members of Scorpius-Centaurus (Sco-Cen), the nearest OB Association (~100-200 pc). Using an X-ray, proper motion and color-magnitude selected sample, we obtained spectra for 361 stars, for which we report spectral classifications and Li and Halpha equivalent widths. We identified 156 new members of Sco-Cen, and recovered 51 previously published members. We have combined these with previously known members to form a sample of 493 solar-mass (~0.7-1.3 Msun) members of Sco-Cen. We investigated the star-formation history of this sample, and re-assessed the ages of the massive main-sequence turn-off and the G-type members in all three subgroups. We performed a census for circumstellar disks in our sample using WISE infrared data and find a protoplanetary disk fraction for K-type stars of 4.4$^{+1.6}_{-0.9}$% for Upper Centaurus-Lupus and Lower Centaurus-Crux at ~16 Myr and 9.0$^{+4.0}_{-2.2}$% for Upper Scorpius at ~10 Myr. These data are consistent with a protoplanetary disk e-folding timescale of ~4-5 Myr for ~1 Msun stars, twice that previously quoted (Mamajek 2009), but consistent with the Bell et al. revised age scale of young clusters. Finally, we construct an age map of Scorpius-Centaurus which clearly reveals substructure consisting of concentrations of younger and older stars. We find evidence for strong age gradients within all three subgroups. None of the subgroups are consistent with being simple, coeval populations which formed in single bursts, but likely represents a multitude of smaller star formation episodes of hundreds to tens of stars each.

Journal ArticleDOI
TL;DR: In this paper, the authors used high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass-metallicity relations (MZR) from z=0-6.
Abstract: We use high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environment (FIRE) project to study the galaxy mass–metallicity relations (MZR) from z=0–6. These simulations include explicit models of the multi-phase ISM, star formation, and stellar feedback. The simulations cover halo masses M_(halo) = 10^9–10^(13) M_☉ and stellar masses M_* = 10^4–10^(11) M_☉ at z = 0 and have been shown to produce many observed galaxy properties from z = 0–6. For the first time, our simulations agree reasonably well with the observed mass–metallicity relations at z = 0–3 for a broad range of galaxy masses. We predict the evolution of the MZR from z = 0–6, as log(Z_(gas)/Z_☉) = 12+log(O/H)-9.0 = 0.35 [log(M_*/M_☉) - 10] + 0.93exp(-0.43z) - 1.05 and log(Z_*/Z_☉) = [Fe=H] + 0.2 = 0.40 [log(M_*/M_☉)-10]+0.67exp(-0.50z)-1.04, for gas-phase and stellar metallicity, respectively. Our simulations suggest that the evolution of MZR is associated with the evolution of stellar/gas mass fractions at different redshifts, indicating the existence of a universal metallicity relation between stellar mass, gas mass, and metallicities. In our simulations, galaxies above M_* = 10^6 M_☉ are able to retain a large fraction of their metals inside the halo, because metal-rich winds fail to escape completely and are recycled into the galaxy. This resolves a long-standing discrepancy between “sub-grid” wind models (and semi-analytic models) and observations, where common sub-grid models cannot simultaneously reproduce the MZR and the stellar mass functions.

Journal ArticleDOI
TL;DR: The relationship between stellar populations and the ionizing flux with which they irradiate their surroundings has profound implications for the evolution of the intergalactic medium (IGM).
Abstract: The relationship between stellar populations and the ionizing flux with which they irradiate their surroundings has profound implications for the evolution of the intergalactic medium (IGM). We quantify the ionizing flux arising from synthetic stellar populations which incorporate the evolution of interacting binary stars. We determine that these show ionizing flux boosted by 60 per cent at 0.05 ≤ Z ≤ 0.3 Z⊙ and a more modest 10–20 per cent at near-solar metallicities relative to star-forming populations in which stars evolve in isolation. The relation of ionizing flux to observables such as 1500 A continuum and ultraviolet spectral slope is sensitive to attributes of the stellar population including age, star formation history and initial mass function (IMF). For a galaxy forming 1 M⊙ yr−1, observed at >100 Myr after the onset of star formation, we predict a production rate of photons capable of ionizing hydrogen, Nion = 1.4 × 1053 s−1 at Z = Z⊙ and 3.5 × 1053 s−1 at 0.1 Z⊙, assuming a Salpeter-like IMF. We evaluate the impact of these issues on the ionization of the IGM, finding that the known galaxy populations can maintain the ionization state of the Universe back to z ∼ 9, assuming that their luminosity functions continue to MUV = −10, and that constraints on the IGM at z ∼ 2–5 can be satisfied with modest Lyman-continuum photon escape fractions of 4–24 per cent depending on assumed metallicity.

Journal ArticleDOI
TL;DR: In this article, a new version of the GALFORM semi-analytical model of galaxy formation is presented, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos and a new empirical star formation law in galaxy disks based on their molecular gas content.
Abstract: We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850μm. This was not previously achieved by a single physical model within the ΛCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z = 0 to z = 6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z = 0, as well as far-IR number counts, and far-UV luminosity functions at z ∼ 3 − 6. Discrepancies are however found in galaxy sizes and metallicities at low luminosities, and in the abundance of low mass galaxies at high-z, suggesting the need for a more sophisticated model of supernova feedback.

Journal ArticleDOI
TL;DR: In this article, ground-based and Swift photometric and spectroscopic observations of the candidate tidal disruption event (TDE) ASASSN-14li, found at the center of PGC043234 (d'90 Mpc) by the All-Sky Automated Survey for SuperNovae (ASAS-SN).
Abstract: We present ground-based and Swift photometric and spectroscopic observations of the candidate tidal disruption event (TDE) ASASSN-14li, found at the center of PGC043234 (d ' 90 Mpc) by the All-Sky Automated Survey for SuperNovae (ASAS-SN). The source had a peak bolometric luminosity of L ' 1044 ergs s

Journal ArticleDOI
TL;DR: The detection of five new fast radio bursts (FRBs) found in the 1.4 GHz High Time Resolution Universe high-latitude survey at Parkes, is presented in this article.
Abstract: The detection of five new fast radio bursts (FRBs) found in the 1.4-GHz High Time Resolution Universe high-latitude survey at Parkes, is presented. The rate implied is 7(-3)(+5) x 10(3) (95 per cent) FRBs sky(-1) d(-1) above a fluence of 0.13 Jy ms for an FRB of 0.128 ms duration to 1.5 Jy ms for 16 ms duration. One of these FRBs has a two-component profile, in which each component is similar to the known population of single component FRBs and the two components are separated by 2.4 +/- 0.4 ms. All the FRB components appear to be unresolved following deconvolution with a scattering tail and accounting for intrachannel smearing. The two-component burst, FRB 121002, also has the highest dispersion measure (1629 pc cm(-3)) of any FRB to-date. Many of the proposed models to explain FRBs use a single high-energy event involving compact objects (such as neutron-star mergers) and therefore cannot easily explain a two-component FRB. Models that are based on extreme versions of flaring, pulsing, or orbital events, however, could produce multiple component profiles. The compatibility of these models and the FRB rate implied by these detections is discussed.

Journal ArticleDOI
TL;DR: In this article, the authors presented the detection of Lyman continuum (LyC) radiation of four other compact star-forming galaxies observed with the Cosmic Origins Spectrograph (COS) onboard the HST.
Abstract: Following our first detection reported in Izotov et al. (2016), we present the detection of Lyman continuum (LyC) radiation of four other compact star-forming galaxies observed with the Cosmic Origins Spectrograph (COS) onboard the Hubble Space Telescope (HST). These galaxies, at redshifts of z~0.3, are characterized by high emission-line flux ratios [OIII]5007/[OII]3727 > 5. The escape fractions of the LyC radiation fesc(LyC) in these galaxies are in the range of ~6%-13%, the highest values found so far in low-redshift star-forming galaxies. Narrow double-peaked Lyalpha emission lines are detected in the spectra of all four galaxies, compatible with predictions for Lyman continuum leakers. We find escape fractions of Lyalpha, fesc(Lyalpha) ~20%-40%, among the highest known for Lyalpha emitters (LAEs). Surface brightness profiles produced from the COS acquisition images reveal bright star-forming regions in the center and exponential discs in the outskirts with disc scale lengths alpha in the range ~0.6-1.4 kpc. Our galaxies are characterized by low metallicity, ~1/8-1/5 solar, low stellar mass ~(0.2 - 4)e9 Msun, high star formation rates SFR~14-36 Msun/yr, and high SFR densities Sigma~2-35 Msun/yr/kpc^2. These properties are comparable to those of high-redshift star-forming galaxies. Finally, our observations, combined with our first detection reported in Izotov et al. (2016), reveal that a selection for compact star-forming galaxies showing high [OIII]5007/[OII]3727 ratios appears to pick up very efficiently sources with escaping Lyman continuum radiation: all five of our selected galaxies are LyC leakers.

Journal ArticleDOI
TL;DR: In this article, the bolometric light curve of 38 stripped-envelope core-collapse supernovae (SE SNe) is recovered and template light curves provided.
Abstract: Literature data are collated for 38 stripped-envelope core-collapse supernovae (SE SNe; i.e. SNe IIb, Ib, Ic and Ic-BL) that have good light curve coverage in more than one optical band. Using bolometric corrections derived in previous work, the bolometric light curve of each SN is recovered and template bolometric light curves provided. Peak light distributions and decay rates are investigated; SNe subtypes are not cleanly distinguished in this parameter space, although some grouping of types does occur and there is a suggestion of a Phillips-like relation for most SNe Ic-BL. The bolometric light curves are modelled with a simple analytical prescription and compared to results from more detailed modelling. Distributions of the explosion parameters shows the extreme nature of SNe Ic-BL in terms of their 56Ni mass and the kinetic energy, however ejected masses are similar to other subtypes. SNe Ib and Ic have very similar distributions of explosion parameters, indicating a similarity in progenitors. SNe~IIb are the most homogeneous subtype and have the lowest average values for 56Ni mass, ejected mass, and kinetic energy. Ejecta masses for each subtype and SE SNe as a whole are inconsistent with those expected from very massive stars. The majority of the ejecta mass distribution is well described by more moderately massive progenitors in binaries, indicating these are the dominant progenitor channel for SE SNe.

Journal ArticleDOI
TL;DR: In this paper, a large galaxy sample was used to explore the sensitivity of the TDE rates to observational uncertainties, such as the parametrization of galaxy light profiles and the stellar mass function.
Abstract: Rates of stellar tidal disruption events (TDEs) by supermassive black holes (SMBHs) due to two-body relaxation are calculated using a large galaxy sample (N=146) in order to explore the sensitivity of the TDE rates to observational uncertainties, such as the parametrization of galaxy light profiles and the stellar mass function. The largest uncertainty arises due to the poorly constrained occupation fraction of SMBHs in low-mass galaxies, which otherwise dominate the total TDE rate. The detection rate of TDE flares by optical surveys is calculated as a function of SMBH mass and other observables for several physically-motivated models of TDE emission. We also quantify the fraction of galaxies that produce deeply penetrating disruption events. If the majority of the detected events are characterized by super-Eddington luminosities (such as disk winds, or synchrotron radiation from an off-axis relativistic jet), then the measured SMBH mass distribution will tightly constrain the low-end SMBH occupation fraction. If Eddington-limited emission channels dominate, however, then the occupation fraction sensitivity is much less pronounced in a flux-limited survey (although still present in a volume-complete event sample). The SMBH mass distribution of the current sample of TDEs, though highly inhomogeneous and encumbered by selection effects, already suggests that Eddington-limited emission channels dominate. Even our most conservative rate estimates appear to be in tension with much lower observationally inferred TDE rates, and we discuss several possible resolutions to this discrepancy.

Journal ArticleDOI
TL;DR: In this paper, the authors studied the spatially resolved excitation properties of the ionised gas in a sample of 646 galaxies using integral field spectroscopy data from SDSS-IV MaNGA.
Abstract: We study the spatially resolved excitation properties of the ionised gas in a sample of 646 galaxies using integral field spectroscopy data from SDSS-IV MaNGA. Making use of Baldwin-Philips-Terlevich diagnostic diagrams we demonstrate the ubiquitous presence of extended (kpc scale) low ionisation emission-line regions (LIERs) in both star forming and quiescent galaxies. In star forming galaxies LIER emission can be associated with diffuse ionised gas, most evident as extra-planar emission in edge-on systems. In addition, we identify two main classes of galaxies displaying LIER emission: `central LIER' (cLIER) galaxies, where central LIER emission is spatially extended, but accompanied by star formation at larger galactocentric distances, and `extended LIER' (eLIER) galaxies, where LIER emission is extended throughout the whole galaxy. In eLIER and cLIER galaxies, LIER emission is associated with radially flat, low H$\alpha$ equivalent width of line emission ($

Journal ArticleDOI
TL;DR: In this article, the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass were analyzed.
Abstract: Using cosmological simulations, we address the properties of high-redshift star-forming galaxies (SFGs) across their main sequence (MS) in the plane of star-formation rate (SFR) versus stellar mass. We relate them to the evolution of galaxies through phases of gas compaction, depletion, possible replenishment, and eventual quenching. We find that the high-SFR galaxies in the upper envelope of the MS are compact, with high gas fractions and short depletion times ("blue nuggets"), while the lower-SFR galaxies in the lower envelope have lower central gas densities, lower gas fractions and longer depletion times, consistent with observed gradients across the MS. Stellar-structure gradients are negligible. The SFGs oscillate about the MS ridge on timescales $\sim0.4~t_{\mathrm{Hubble}}$ ($\sim1$ Gyr at $z\sim3$). The propagation upwards is due to gas compaction, triggered, e.g., by mergers, counter-rotating streams, and/or violent disc instabilities. The downturn at the upper envelope is due to central gas depletion by peak star formation and outflows while inflow from the shrunken gas disc is suppressed. An upturn at the lower envelope can occur once the extended disc has been replenished by fresh gas and a new compaction can be triggered, namely as long as the replenishment time is shorter than the depletion time. The mechanisms of gas compaction, depletion and replenishment confine the SFGs to the narrow ($\pm0.3$ dex) MS. Full quenching occurs in massive haloes ($M_{\mathrm{vir}}>10^{11.5}~M_\odot$) and/or at low redshifts ($z<3$), where the replenishment time is long compared to the depletion time, explaining the observed bending down of the MS at the massive end.

Journal ArticleDOI
TL;DR: In this paper, the authors presented a UV-extended E-MILES stellar population synthesis model covering the spectral range of 1680-50 000 angstrom at moderately high resolution.
Abstract: We present UV-extended E-MILES stellar population synthesis models covering the spectral range.. 1680-50 000 angstrom at moderately high resolution. We employ the NGSL space-based stellar library to compute spectra of single-age, single-metallicity stellar populations in the wavelength range from 1680 to 3540 angstrom. These models represent a significant improvement in resolution and age/metallicity coverage over previous studies based on earlier space-based libraries. These model spectra were joined with those we computed in the visible using MILES, and other empirical libraries for redder wavelengths. The models span the metallicity range -1.79 <= [M/H] <= +0.26 and ages above 30 Myr, for a suite of initial mass function types with varying slopes. We focus on the behaviour of colours, spectra and line-strength indices in the UV range as a function of relevant stellar population parameters. Whereas some indices strengthen with increasing age and metallicity, as most metallicity indicators in the visible, other indices peak around 3 Gyr for metal-rich stellar populations, such as Mg at 2800 angstrom. Our models provide reasonably good fits to the integrated colours and most line strengths of the stellar clusters of the Milky Way and Large Magellanic Cloud. Our full spectrum fits in the UV range for a representative set of early-type galaxies (ETGs) of varying mass yield age and metallicity estimates in very good agreement with those obtained in the optical range. The comparison of UV colours and line strengths of massive ETGs with our models reveals the presence of young stellar components, with ages in the range 0.1-0.5 Gyr and mass fractions 0.1-0.5 per cent, on the top of an old stellar population.

Journal ArticleDOI
TL;DR: In this article, the Parkes Pulsar Timing Array was used to estimate the parallax of 20 millisecond pulsars, including the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars.
Abstract: We present timing models for 20 millisecond pulsars in the Parkes Pulsar Timing Array. The precision of the parameter measurements in these models has been improved over earlier results by using longer data sets and modelling the non-stationary noise. We describe a new noise modelling procedure and demonstrate its effectiveness using simulated data. Our methodology includes the addition of annual dispersion measure (DM) variations to the timing models of some pulsars. We present the first significant parallax measurements for PSRs J1024-0719, J1045-4509, J1600-3053, J1603-7202, and J1730-2304, as well as the first significant measurements of some post-Keplerian orbital parameters in six binary pulsars, caused by kinematic effects. Improved Shapiro delay measurements have resulted in much improved pulsar mass measurements, particularly for PSRs J0437-4715 and J1909-3744 with M-p = 1.44 +/- 0.07 and 1.47 +/- 0.03 M-circle dot, respectively. The improved orbital period-derivative measurement for PSR J0437-4715 results in a derived distance measurement at the 0.16 per cent level of precision, D = 156.79 +/- 0.25 pc, one of the most fractionally precise distance measurements of any star to date.

Journal ArticleDOI
TL;DR: In this article, the role of pair-enriched bubbles produced by young magnetars, rapidly-rotating neutron stars, and magnetized white dwarfs in the escape of high-energy gamma-rays and radio waves was studied.
Abstract: Tenuous wind bubbles, which are formed by the spin-down activity of central compact remnants, are relevant in some models of fast radio bursts (FRBs) and super-luminous supernovae. We study their high-energy signatures, focusing on the role of pair-enriched bubbles produced by young magnetars, rapidly-rotating neutron stars, and magnetized white dwarfs. (i) First, we study the nebular properties and the conditions allowing for escape of high-energy gamma-rays and radio waves, showing that their escape is possible for nebulae with ages of >10-100 yr. In the rapidly-rotating neutron star scenario, we find that radio emission from the quasi-steady nebula itself may be bright enough to be detected especially at sub-mm frequencies, which is relevant as a possible counterpart of pulsar-driven SNe and FRBs. (ii) Second, we consider the fate of bursting emission in the nebulae. We suggest that an impulsive burst may lead to a highly relativistic flow, which would interact with the nebula. If the shocked nebula is still relativistic, pre-existing non-thermal particles in the nebula can be significantly boosted by the forward shock, leading to short-duration (maybe millisecond or longer) high-energy gamma-ray flashes. Possible dissipation at the reverse shock may also lead to gamma-ray emission. (iii) After such flares, interactions with the baryonic ejecta may lead to afterglow emission with a duration of days to weeks. In the magnetar scenario, this burst-in-bubble model leads to the expectation that nearby (<10-100 Mpc) high-energy gamma-ray flashes may be detected by HAWC and CTA, and the subsequent afterglow emission may be seen by radio telescopes such as VLA. (iv) Finally, we discuss several implications specific to FRBs, including constraints on the emission regions and limits on soft gamma-ray counterparts.

Journal ArticleDOI
TL;DR: The mufasa suite of cosmological hydrodynamic simulations as mentioned in this paper employs the gizmo meshless finite mass (MFM) code including H_2-based star formation, nine-element chemical evolution, two-phase kinetic outflows following scalings from the Feedback in Realistic Environments zoom simulations, and evolving halo mass-based quenching.
Abstract: We present the mufasa suite of cosmological hydrodynamic simulations, which employs the gizmo meshless finite mass (MFM) code including H_2-based star formation, nine-element chemical evolution, two-phase kinetic outflows following scalings from the Feedback in Realistic Environments zoom simulations, and evolving halo mass-based quenching. Our fiducial (50 h^(−1) Mpc)^3 volume is evolved to z = 0 with a quarter billion elements. The predicted galaxy stellar mass functions (GSMFs) reproduces observations from z = 4 → 0 to ≲ 1.2σ in cosmic variance, providing an unprecedented match to this key diagnostic. The cosmic star formation history and stellar mass growth show general agreement with data, with a strong archaeological downsizing trend such that dwarf galaxies form the majority of their stars after z ∼ 1. We run 25 and 12.5 h^(−1) Mpc volumes to z = 2 with identical feedback prescriptions, the latter resolving all hydrogen-cooling haloes, and the three runs display fair resolution convergence. The specific star formation rates broadly agree with data at z = 0, but are underpredicted at z ∼ 2 by a factor of 3, re-emphasizing a longstanding puzzle in galaxy evolution models. We compare runs using MFM and two flavours of smoothed particle hydrodynamics, and show that the GSMF is sensitive to hydrodynamics methodology at the ∼×2 level, which is sub-dominant to choices for parametrizing feedback.

Journal ArticleDOI
TL;DR: In this paper, the critical role of active galactic nuclei (AGN) in setting up the correct galaxy morphology for the massive end of the population was identified by comparing the state-of-the-art hydrodynamical cosmological simulations Horizon-AGN and Horizon-noAGN.
Abstract: The interplay between cosmic gas accretion on to galaxies and galaxy mergers drives the observed morphological diversity of galaxies. By comparing the state-of-the-art hydrodynamical cosmological simulations Horizon-AGN and Horizon-noAGN, we unambiguously identify the critical role of active galactic nuclei (AGN) in setting up the correct galaxy morphology for the massive end of the population. With AGN feedback, typical kinematic and morpho-metric properties of galaxy populations as well as the galaxy-halo mass relation are in much better agreement with observations. Only AGN feedback allows massive galaxies at the centre of groups and clusters to become ellipticals, while without AGN feedback those galaxies reform discs. It is the merger-enhanced AGN activity that is able to freeze the morphological type of the post-merger remnant by durably quenching its quiescent star formation. Hence morphology is shown to be driven not only by mass but also by the nature of cosmic accretion: at constant galaxy mass, ellipticals are galaxies that are mainly assembled through mergers, while discs are preferentially built from the in situ star formation fed by smooth cosmic gas infall.

Journal ArticleDOI
TL;DR: In this paper, the authors acknowledge support from the STFC and the Ernest Rutherford Fellowship scheme for the SDSS and SDSSE-II, and the Higher Education Funding Council for England.
Abstract: CLM acknowledges support from the STFC Consolidated Grant (Ref. St/M001229/1). NPR acknowledges support from the STFC and the Ernest Rutherford Fellowship scheme. KH acknowledges support from STFC grant ST/M001296/1. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the US Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England.

Journal ArticleDOI
TL;DR: In this article, the authors calculated the rates, time-scales and mass distributions for binary black hole (BH) mergers as a function of metallicity using the Binary Population and Spectral Synthesis code, BPASS.
Abstract: Using the Binary Population and Spectral Synthesis code, BPASS, we have calculated the rates, time-scales and mass distributions for binary black hole (BH) mergers as a function of metallicity. We consider these in the context of the recently reported first Laser Interferometer Gravitational-Wave Observatory (LIGO) event detection. We find that the event has a very low probability of arising from a stellar population with initial metallicity mass fraction above Z = 0.010 (Z ≳ 0.5 Z⊙). Binary BH merger events with the reported masses are most likely in populations below 0.008 (Z ≲ 0.4 Z⊙). Events of this kind can occur at all stellar population ages from 3 Myr up to the age of the Universe, but constitute only 0.1–0.4 per cent of binary BH mergers between metallicities of Z = 0.001 and 0.008. However at metallicity Z = 10−4, 26 per cent of binary BH mergers would be expected to have the reported masses. At this metallicity, the progenitor merger times can be close to ≈10 Gyr and rotationally mixed stars evolving through quasi-homogeneous evolution, due to mass transfer in a binary, dominate the rate. The masses inferred for the BHs in the binary progenitor of GW 150914 are amongst the most massive expected at anything but the lowest metallicities in our models. We discuss the implications of our analysis for the electromagnetic follow-up of future LIGO event detections.

Journal ArticleDOI
TL;DR: In this paper, the authors made use of data from the first public release of the WASP data (Butters et al. 2010) as provided by the NASA Exoplanet Archive, which is operated by the California Institute of Technology under contract with the National Aeronautics and Space Administration under the ERC grant number 279973.
Abstract: TSB acknowledges support provided through NASA grant ADAP12-0172 and ADAP14-0245. MCW and GMK acknowledge the support of the European Union through ERC grant number 279973. The authors acknowledge support from the Hungarian Research Grants OTKA K-109276, OTKA K-113117, the Lendulet-2009 and Lendulet-2012 Program (LP2012-31) of the Hungarian Academy of Sciences, the Hungarian National Research, Development and Innovation Office – NKFIH K-115709, and the ESA PECS Contract No. 4000110889/14/NL/NDe. This work was supported by the Momentum grant of the MTA CSFK Lendulet Disc Research Group. GH acknowledges support by the Polish NCN grant 2011/01/B/ST9/05448. Based on observations made with the NOT, operated by the Nordic Optical Telescope Scientific Association at the Observatorio del Roque de los Muchachos, La Palma, Spain, of the Instituto de Astrofisica de Canarias. This research made use of The DASCH project; we are also grateful for partial support from NSF grants AST-0407380, AST-0909073, and AST-1313370. The research leading to these results has received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under grant agreements no. 269194 (IRSES/ASK) and no. 312844 (SPACEINN). We thank Scott Dahm, Julie Rivera, and the Keck Observatory staff for their assistance with these observations. This research was supported in part by NSF grant AST-0909222 awarded to M. Liu. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. KS gratefully acknowledges support from Swiss National Science Foundation Grant PP00P2_138979/1. HJD and DN acknowledge support by grant AYA2012-39346-C02-02 of the Spanish Secretary of State for R&D&i (MINECO). This paper makes use of data from the first public release of the WASP data (Butters et al. 2010) as provided by the WASP consortium and services at the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, and NEOWISE, which is a project of the Jet Propulsion Laboratory/California Institute of Technology. WISE and NEOWISE are funded by the National Aeronautics and Space Administration. This research made use of the SIMBAD and VIZIER Astronomical Databases, operated at CDS, Strasbourg, France (http://cdsweb.u-strasbg.fr/), and of NASA's Astrophysics Data System.