@article{11557,
  abstract     = {Deep narrow-band surveys have revealed a large population of faint Ly α emitters (LAEs) in the distant Universe, but relatively little is known about the most luminous sources (⁠LLyα≳1042.7 erg s−1; LLyα≳L∗Lyα⁠). Here we present the spectroscopic follow-up of 21 luminous LAEs at z ∼ 2–3 found with panoramic narrow-band surveys over five independent extragalactic fields (≈4 × 106 Mpc3 surveyed at z ∼ 2.2 and z ∼ 3.1). We use WHT/ISIS, Keck/DEIMOS, and VLT/X-SHOOTER to study these sources using high ionization UV lines. Luminous LAEs at z ∼ 2–3 have blue UV slopes (⁠β=−2.0+0.3−0.1⁠) and high Ly α escape fractions (⁠50+20−15 per cent) and span five orders of magnitude in UV luminosity (MUV ≈ −19 to −24). Many (70 per cent) show at least one high ionization rest-frame UV line such as C IV, N V, C III], He II or O III], typically blue-shifted by ≈100–200 km s−1 relative to Ly α. Their Ly α profiles reveal a wide variety of shapes, including significant blue-shifted components and widths from 200 to 4000 km s−1. Overall, 60 ± 11  per cent appear to be active galactic nucleus (AGN) dominated, and at LLyα > 1043.3 erg s−1 and/or MUV < −21.5 virtually all LAEs are AGNs with high ionization parameters (log U = 0.6 ± 0.5) and with metallicities of ≈0.5 − 1 Z⊙. Those lacking signatures of AGNs (40 ± 11  per cent) have lower ionization parameters (⁠logU=−3.0+1.6−0.9 and log ξion = 25.4 ± 0.2) and are apparently metal-poor sources likely powered by young, dust-poor ‘maximal’ starbursts. Our results show that luminous LAEs at z ∼ 2–3 are a diverse population and that 2×L∗Lyα and 2×M∗UV mark a sharp transition in the nature of LAEs, from star formation dominated to AGN dominated.},
  author       = {Sobral, David and Matthee, Jorryt J and Darvish, Behnam and Smail, Ian and Best, Philip N and Alegre, Lara and Röttgering, Huub and Mobasher, Bahram and Paulino-Afonso, Ana and Stroe, Andra and Oteo, Iván},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: active, galaxies: evolution, galaxies: high-redshift, galaxies: ISM, galaxies: starburst, cosmology: observations},
  number       = {2},
  pages        = {2817--2840},
  publisher    = {Oxford University Press},
  title        = {{The nature of luminous Ly α emitters at z ∼ 2–3: Maximal dust-poor starbursts and highly ionizing AGN}},
  doi          = {10.1093/mnras/sty782},
  volume       = {477},
  year         = {2018},
}

@article{11558,
  abstract     = {We present and explore deep narrow- and medium-band data obtained with the Subaru and the Isaac Newton Telescopes in the ∼2 deg2 COSMOS field. We use these data as an extremely wide, low-resolution (R ∼ 20–80) Integral Field Unit survey to slice through the COSMOS field and obtain a large sample of ∼4000 Ly α emitters (LAEs) from z ∼ 2 to 6 in 16 redshift slices (SC4K). We present new Ly α luminosity functions (LFs) covering a comoving volume of ∼108 Mpc3. SC4K extensively complements ultradeep surveys, jointly covering over 4 dex in Ly α luminosity and revealing a global (2.5 < z < 6) synergy LF with α=−1.93+0.12−0.12⁠, log10Φ∗Lyα=−3.45+0.22−0.29 Mpc−3, and log10L∗Lyα=42.93+0.15−0.11 erg s−1. The Schechter component of the Ly α LF reveals a factor ∼5 rise in L∗Lyα and a ∼7 × decline in Φ∗Lyα from z ∼ 2 to 6. The data reveal an extra power-law (or Schechter) component above LLy α ≈ 1043.3 erg s−1 at z ∼ 2.2–3.5 and we show that it is partially driven by X-ray and radio active galactic nucleus (AGN), as their Ly α LF resembles the excess. The power-law component vanishes and/or is below our detection limits above z > 3.5, likely linked with the evolution of the AGN population. The Ly α luminosity density rises by a factor ∼2 from z ∼ 2 to 3 but is then found to be roughly constant (⁠1.1+0.2−0.2×1040 erg s−1 Mpc−3) to z ∼ 6, despite the ∼0.7 dex drop in ultraviolet (UV) luminosity density. The Ly α/UV luminosity density ratio rises from 4 ± 1 per cent to 30 ± 6 per cent from z ∼ 2.2 to 6. Our results imply a rise of a factor of ≈2 in the global ionization efficiency (ξion) and a factor ≈4 ± 1 in the Ly α escape fraction from z ∼ 2 to 6, hinting for evolution in both the typical burstiness/stellar populations and even more so in the typical interstellar medium conditions allowing Ly α photons to escape.},
  author       = {Sobral, David and Santos, Sérgio and Matthee, Jorryt J and Paulino-Afonso, Ana and Ribeiro, Bruno and Calhau, João and Khostovan, Ali A},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: formation, galaxies: high-redshift, galaxies: luminosity function, mass function, galaxies: statistics},
  number       = {4},
  pages        = {4725--4752},
  publisher    = {Oxford University Press},
  title        = {{Slicing COSMOS with SC4K: The evolution of typical Ly α emitters and the Ly α escape fraction from z ∼ 2 to 6}},
  doi          = {10.1093/mnras/sty378},
  volume       = {476},
  year         = {2018},
}

@article{11584,
  abstract     = {Observations show that star-forming galaxies reside on a tight 3D plane between mass, gas-phase metallicity, and star formation rate (SFR), which can be explained by the interplay between metal-poor gas inflows, SFR and outflows. However, different metals are released on different time-scales, which may affect the slope of this relation. Here, we use central, star-forming galaxies with Mstar = 109.0–10.5 M⊙ from the EAGLE hydrodynamical simulation to examine 3D relations between mass, SFR, and chemical enrichment using absolute and relative C, N, O, and Fe abundances. We show that the scatter is smaller when gas-phase α-enhancement is used rather than metallicity. A similar plane also exists for stellar α-enhancement, implying that present-day specific SFRs are correlated with long time-scale star formation histories. Between z = 0 and 1, the α-enhancement plane is even more insensitive to redshift than the plane using metallicity. However, it evolves at z > 1 due to lagging iron yields. At fixed mass, galaxies with higher SFRs have star formation histories shifted towards late times, are more α-enhanced, and this α-enhancement increases with redshift as observed. These findings suggest that relations between physical properties inferred from observations may be affected by systematic variations in α-enhancements.},
  author       = {Matthee, Jorryt J and Schaye, Joop},
  issn         = {1745-3933},
  journal      = {Monthly Notices of the Royal Astronomical Society: Letters},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: abundances, galaxies: evolution, galaxies: formation, galaxies: star formation},
  number       = {1},
  pages        = {L34 -- L39},
  publisher    = {Oxford University Press},
  title        = {{Star-forming galaxies are predicted to lie on a fundamental plane of mass, star formation rate, and α-enhancement}},
  doi          = {10.1093/mnrasl/sly093},
  volume       = {479},
  year         = {2018},
}

@article{11618,
  abstract     = {Asteroseismology provides global stellar parameters such as masses, radii, or surface gravities using mean global seismic parameters and effective temperature for thousands of low-mass stars (0.8 M⊙ < M < 3 M⊙). This methodology has been successfully applied to stars in which acoustic modes excited by turbulent convection are measured. Other methods such as the Flicker technique can also be used to determine stellar surface gravities, but only works for log g above 2.5 dex. In this work, we present a new metric called FliPer (Flicker in spectral power density, in opposition to the standard Flicker measurement which is computed in the time domain); it is able to extend the range for which reliable surface gravities can be obtained (0.1 < log g < 4.6 dex) without performing any seismic analysis for stars brighter than Kp < 14. FliPer takes into account the average variability of a star measured in the power density spectrum in a given range of frequencies. However, FliPer values calculated on several ranges of frequency are required to better characterize a star. Using a large set of asteroseismic targets it is possible to calibrate the behavior of surface gravity with FliPer through machine learning. This calibration made with a random forest regressor covers a wide range of surface gravities from main-sequence stars to subgiants and red giants, with very small uncertainties from 0.04 to 0.1 dex. FliPer values can be inserted in automatic global seismic pipelines to either give an estimation of the stellar surface gravity or to assess the quality of the seismic results by detecting any outliers in the obtained νmax values. FliPer also constrains the surface gravities of main-sequence dwarfs using only long-cadence data for which the Nyquist frequency is too low to measure the acoustic-mode properties.},
  author       = {Bugnet, Lisa Annabelle and García, R. A. and Davies, G. R. and Mathur, S. and Corsaro, E. and Hall, O. J. and Rendle, B. M.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, asteroseismology / methods, data analysis / stars, oscillations},
  publisher    = {EDP Sciences},
  title        = {{FliPer: A global measure of power density to estimate surface gravities of main-sequence solar-like stars and red giants}},
  doi          = {10.1051/0004-6361/201833106},
  volume       = {620},
  year         = {2018},
}

@article{11619,
  abstract     = {We report on the confirmation and mass determination of π Men c, the first transiting planet discovered by NASA’s TESS space mission. π Men is a naked-eye (V = 5.65 mag), quiet G0 V star that was previously known to host a sub-stellar companion (π Men b) on a longperiod (Porb = 2091 days), eccentric (e = 0.64) orbit. Using TESS time-series photometry, combined with Gaia data, published UCLES at AAT Doppler measurements, and archival HARPS at ESO-3.6m radial velocities, we found that π Men c is a close-in planet with an orbital period of Porb = 6.27 days, a mass of Mc = 4.52 ± 0.81 M⊕, and a radius of Rc = 2.06 ± 0.03 R⊕. Based on the planet’s orbital period and size, π Men c is a super-Earth located at, or close to, the radius gap, while its mass and bulk density suggest it may have held on to a significant atmosphere. Because of the brightness of the host star, this system is highly suitable for a wide range of further studies to characterize the planetary atmosphere and dynamical properties. We also performed an asteroseismic analysis of the TESS data and detected a hint of power excess consistent with the seismic values expected for this star, although this result depends on the photometric aperture used to extract the light curve. This marginal detection is expected from pre-launch simulations hinting at the asteroseismic potential of the TESS mission for longer, multi-sector observations and/or for more evolved bright stars.},
  author       = {Gandolfi, D. and Barragán, O. and Livingston, J. H. and Fridlund, M. and Justesen, A. B. and Redfield, S. and Fossati, L. and Mathur, S. and Grziwa, S. and Cabrera, J. and García, R. A. and Persson, C. M. and Van Eylen, V. and Hatzes, A. P. and Hidalgo, D. and Albrecht, S. and Bugnet, Lisa Annabelle and Cochran, W. D. and Csizmadia, Sz. and Deeg, H. and Eigmüller, Ph. and Endl, M. and Erikson, A. and Esposito, M. and Guenther, E. and Korth, J. and Luque, R. and Montañes Rodríguez, P. and Nespral, D. and Nowak, G. and Pätzold, M. and Prieto-Arranz, J.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, planetary systems / planets and satellites, detection / planets and satellites, fundamental parameters / planets and satellites, terrestrial planets / stars, fundamental parameters},
  publisher    = {EDP Sciences},
  title        = {{TESS’s first planet: A super-Earth transiting the naked-eye star π Mensae}},
  doi          = {10.1051/0004-6361/201834289},
  volume       = {619},
  year         = {2018},
}

@article{11620,
  abstract     = {We report the discovery and characterization of HD 89345b (K2-234b; EPIC 248777106b), a Saturn-sized planet orbiting a slightly evolved star. HD 89345 is a bright star (V = 9.3 mag) observed by the K2 mission with 1 min time sampling. It exhibits solar-like oscillations. We conducted asteroseismology to determine the parameters of the star, finding the mass and radius to be 1.12+0.04−0.01M⊙ and 1.657+0.020−0.004R⊙⁠, respectively. The star appears to have recently left the main sequence, based on the inferred age, 9.4+0.4−1.3Gyr⁠, and the non-detection of mixed modes. The star hosts a ‘warm Saturn’ (P = 11.8 d, Rp = 6.86 ± 0.14 R⊕). Radial-velocity follow-up observations performed with the FIbre-fed Echelle Spectrograph, HARPS, and HARPS-N spectrographs show that the planet has a mass of 35.7 ± 3.3 M⊕. The data also show that the planet’s orbit is eccentric (e ≈ 0.2). An investigation of the rotational splitting of the oscillation frequencies of the star yields no conclusive evidence on the stellar inclination angle. We further obtained Rossiter–McLaughlin observations, which result in a broad posterior of the stellar obliquity. The planet seems to confirm to the same patterns that have been observed for other sub-Saturns regarding planet mass and multiplicity, orbital eccentricity, and stellar metallicity.},
  author       = {Van Eylen, V and Dai, F and Mathur, S and Gandolfi, D and Albrecht, S and Fridlund, M and García, R A and Guenther, E and Hjorth, M and Justesen, A B and Livingston, J and Lund, M N and Pérez Hernández, F and Prieto-Arranz, J and Regulo, C and Bugnet, Lisa Annabelle and Everett, M E and Hirano, T and Nespral, D and Nowak, G and Palle, E and Silva Aguirre, V and Trifonov, T and Winn, J N and Barragán, O and Beck, P G and Chaplin, W J and Cochran, W D and Csizmadia, S and Deeg, H and Endl, M and Heeren, P and Grziwa, S and Hatzes, A P and Hidalgo, D and Korth, J and Mathis, S and Montañes Rodriguez, P and Narita, N and Patzold, M and Persson, C M and Rodler, F and Smith, A M S},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, asteroseismology, planets and satellites: composition, planets and satellites: formation, planets and satellites: fundamental parameters},
  number       = {4},
  pages        = {4866--4880},
  publisher    = {Oxford University Press},
  title        = {{HD 89345: A bright oscillating star hosting a transiting warm Saturn-sized planet observed by K2}},
  doi          = {10.1093/mnras/sty1390},
  volume       = {478},
  year         = {2018},
}

@article{13473,
  abstract     = {Stripped-envelope stars form in binary systems after losing mass through Roche-lobe overflow. They bear astrophysical significance as sources of UV and ionizing radiation in older stellar populations and, if sufficiently massive, as stripped supernova progenitors. Binary evolutionary models predict that they are common, but only a handful of subdwarfs with B-type companions are known. The question is whether a large population of such systems has evaded detection as a result of biases, or whether the model predictions are wrong. We reanalyze the well-studied post-interaction binary φ Persei. Recently, new data have improved the orbital solution of the system, which contains an ~1.2M⊙ stripped-envelope star and a rapidly rotating ~9.6M⊙ Be star. We compare with an extensive grid of evolutionary models using a Bayesian approach and constrain the initial masses of the progenitor to 7.2 ± 0.4M⊙ for the stripped star and 3.8 ± 0.4M⊙ for the Be star. The system must have evolved through near-conservative mass transfer. These findings are consistent with earlier studies. The age we obtain, 57 ± 9 Myr, is in excellent agreement with the age of the α Persei cluster. We note that neither star was initially massive enough to produce a core-collapse supernova, but mass exchange pushed the Be star above the mass threshold. We find that the subdwarf is overluminous for its mass by almost an order of magnitude, compared to the expectations for a helium core burning star. We can only reconcile this if the subdwarf resides in a late phase of helium shell burning, which lasts only 2–3% of the total lifetime as a subdwarf. Assuming continuous star formation implies that up to ~50 less evolved, dimmer subdwarfs exist for each system similar to φ Persei, but have evaded detection so far. Our findings can be interpreted as a strong indication that a substantial population of stripped-envelope stars indeed exists, but has so far evaded detection because of observational biases and lack of large-scale systematic searches.},
  author       = {Schootemeijer, A. and Götberg, Ylva Louise Linsdotter and de Mink, S. E. and Gies, D. and Zapartas, E.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Clues about the scarcity of stripped-envelope stars from the evolutionary state of the sdO+Be binary system φ Persei}},
  doi          = {10.1051/0004-6361/201731194},
  volume       = {615},
  year         = {2018},
}

@article{13474,
  abstract     = {Recent surveys of the Magellanic Clouds have revealed a subtype of Wolf–Rayet (WR) star with peculiar properties. WN3/O3 spectra exhibit both WR-like emission and O3 V-like absorption – but at lower luminosity than O3 V or WN stars. We examine the projected spatial distribution of WN3/O3 stars in the Large Magellanic Cloud as compared to O-type stars. Surprisingly, WN3/O3 stars are among the most isolated of all classes of massive stars; they have a distribution similar to red supergiants dominated by initial masses of 10–15 M⊙, and are far more dispersed than classical WR stars or luminous blue variables. Their lack of association with clusters of O-type stars suggests strongly that WN3/O3 stars are not the descendants of single massive stars (30 M⊙ or above). Instead, they are likely products of interacting binaries at lower initial mass (10–18 M⊙). Comparison with binary models suggests a probable origin with primaries in this mass range that were stripped of their H envelopes through non-conservative mass transfer by a low-mass secondary. We show that model spectra and positions on the Hertzsprung–Russell diagram for binary-stripped stars are consistent with WN3/O3 stars. Monitoring radial velocities with high-resolution spectra can test for low-mass companions or runaway velocities. With lower initial mass and environments that avoid very massive stars, the WN3/O3 stars fit expectations for progenitors of Type Ib and possibly Type Ibn supernovae.},
  author       = {Smith, Nathan and Götberg, Ylva Louise Linsdotter and de Mink, Selma E},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {1},
  pages        = {772--782},
  publisher    = {Oxford University Press},
  title        = {{Extreme isolation of WN3/O3 stars and implications for their evolutionary origin as the elusive stripped binaries}},
  doi          = {10.1093/mnras/stx3181},
  volume       = {475},
  year         = {2018},
}

@article{13475,
  abstract     = {Stars stripped of their hydrogen-rich envelope through interaction with a binary companion are generally not considered when accounting for ionizing radiation from stellar populations, despite the expectation that stripped stars emit hard ionizing radiation, form frequently, and live 10–100 times longer than single massive stars. We compute the first grid of evolutionary and spectral models specially made for stars stripped in binaries for a range of progenitor masses (2–20 M⊙) and metallicities ranging from solar to values representative for pop II stars. For stripped stars with masses in the range 0.3–7 M⊙, we find consistently high effective temperatures (20 000–100 000 K, increasing with mass), small radii (0.2–1 R⊙), and high bolometric luminosities, comparable to that of their progenitor before stripping. The spectra show a continuous sequence that naturally bridges subdwarf-type stars at the low-mass end and Wolf-Rayet-like spectra at the high-mass end. For intermediate masses we find hybrid spectral classes showing a mixture of absorption and emission lines. These appear for stars with mass-loss rates of 10−8−10−6 M⊙ yr−1, which have semi-transparent atmospheres. At low metallicity, substantial hydrogen-rich layers are left at the surface and we predict spectra that resemble O-type stars instead. We obtain spectra undistinguishable from subdwarfs for stripped stars with masses up to 1.7 M⊙, which questions whether the widely adopted canonical value of 0.47 M⊙ is uniformly valid. Only a handful of stripped stars of intermediate mass have currently been identified observationally. Increasing this sample will provide necessary tests for the physics of interaction, internal mixing, and stellar winds. We use our model spectra to investigate the feasibility to detect stripped stars next to an optically bright companion and recommend systematic searches for their UV excess and possible emission lines, most notably HeII λ4686 in the optical and HeII λ1640 in the UV. Our models are publicly available for further investigations or inclusion in spectral synthesis simulations.},
  author       = {Götberg, Ylva Louise Linsdotter and de Mink, S. E. and Groh, J. H. and Kupfer, T. and Crowther, P. A. and Zapartas, E. and Renzo, M.},
  issn         = {1432-0746},
  journal      = {Astronomy &amp; Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Spectral models for binary products: Unifying subdwarfs and Wolf-Rayet stars as a sequence of stripped-envelope stars}},
  doi          = {10.1051/0004-6361/201732274},
  volume       = {615},
  year         = {2018},
}

@article{11518,
  abstract     = {We present spectroscopic follow-up observations of CR7 with ALMA, targeted at constraining the infrared (IR) continuum and [C II]158 mm line-emission at high spatial resolution matched to the HST/WFC3 imaging. CR7 is a luminous Lyα emitting galaxy at z = 6.6 that consists of three separated UV-continuum components. Our observations reveal several well-separated components of [C II] emission. The two most luminous components in [C II] coincide with the brightest UV components (A and B), blueshifted by »150 km s−1 with respect to the
peak of Lyα emission. Other [C II] components are observed close to UV clumps B and C and are blueshifted by »300 and ≈80 km s−1 with respect to the systemic redshift. We do not detect FIR continuum emission due to dust with a 3σ limiting luminosity LIR T L d 35 K 3.1 10 = <´ 10 ( ) . This allows us to mitigate uncertainties in the dust-corrected SFR and derive SFRs for the three UV clumps A, B, and C of 28, 5, and 7 M yr−1. All clumps have [C II] luminosities consistent within the scatter observed in the local relation between SFR and L[ ] C II , implying that strong Lyα emission does not necessarily anti-correlate with [C II] luminosity. Combining
our measurements with the literature, we show that galaxies with blue UV slopes have weaker [C II] emission at fixed SFR, potentially due to their lower metallicities and/or higher photoionization. Comparison with hydrodynamical simulations suggests that CR7ʼs clumps have metallicities of 0.1 Z Z 0.2 < < . The observed ISM structure of CR7 indicates that we are likely witnessing the build up of a central galaxy in the early universe through complex accretion of satellites.},
  author       = {Matthee, Jorryt J and Sobral, D. and Boone, F. and Röttgering, H. and Schaerer, D. and Girard, M. and Pallottini, A. and Vallini, L. and Ferrara, A. and Darvish, B. and Mobasher, B.},
  issn         = {1538-4357},
  journal      = {The Astrophysical Journal},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, dark ages, reionization, first stars – galaxies: formation – galaxies: high-redshift – galaxies: ISM – galaxies: kinematics and dynamics},
  number       = {2},
  publisher    = {IOP Publishing},
  title        = {{ALMA reveals metals yet no dust within multiple components in CR7}},
  doi          = {10.3847/1538-4357/aa9931},
  volume       = {851},
  year         = {2017},
}

@article{11561,
  abstract     = {We present a sample of ∼1000 emission-line galaxies at z = 0.4–4.7 from the ∼0.7deg2 High-z Emission-Line Survey in the Boötes field identified with a suite of six narrow-band filters at ≈0.4–2.1 μm. These galaxies have been selected on their Ly α (73), [O II] (285), H β/[O III] (387) or H α (362) emission line, and have been classified with optical to near-infrared colours. A subsample of 98 sources have reliable redshifts from multiple narrow-band (e.g. [O II]–H α) detections and/or spectroscopy. In this survey paper, we present the observations, selection and catalogues of emitters. We measure number densities of Ly α, [O II], H β/[O III] and H α and confirm strong luminosity evolution in star-forming galaxies from z ∼ 0.4 to ∼5, in agreement with previous results. To demonstrate the usefulness of dual-line emitters, we use the sample of dual [O II]–H α emitters to measure the observed [O II]/H α ratio at z = 1.47. The observed [O II]/H α ratio increases significantly from 0.40 ± 0.01 at z = 0.1 to 0.52 ± 0.05 at z = 1.47, which we attribute to either decreasing dust attenuation with redshift, or due to a bias in the (typically) fibre measurements in the local Universe that only measure the central kpc regions. At the bright end, we find that both the H α and Ly α number densities at z ≈ 2.2 deviate significantly from a Schechter form, following a power law. We show that this is driven entirely by an increasing X-ray/active galactic nucleus fraction with line luminosity, which reaches ≈100 per cent at line luminosities L ≳ 3 × 1044 erg s−1.},
  author       = {Matthee, Jorryt J and Sobral, David and Best, Philip and Smail, Ian and Bian, Fuyan and Darvish, Behnam and Röttgering, Huub and Fan, Xiaohui},
  issn         = {0035-8711},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics galaxies, active, galaxies, evolution, galaxies, high-redshift, galaxies, luminosity function, mass function, galaxies: star formation},
  number       = {1},
  pages        = {629--649},
  publisher    = {Oxford University Press},
  title        = {{Boötes-HiZELS: An optical to near-infrared survey of emission-line galaxies at z = 0.4–4.7}},
  doi          = {10.1093/mnras/stx1569},
  volume       = {471},
  year         = {2017},
}

@article{11562,
  abstract     = {We present the CAlibrating LYMan-α with Hα (CALYMHA) pilot survey and new results on Lyman α (Lyα) selected galaxies at z ∼ 2. We use a custom-built Lyα narrow-band filter at the Isaac Newton Telescope, designed to provide a matched volume coverage to the z = 2.23 Hα HiZELS survey. Here, we present the first results for the COSMOS and UDS fields. Our survey currently reaches a 3σ line flux limit of ∼4 × 10−17 erg s−1 cm−2, and a Lyα luminosity limit of ∼1042.3 erg s−1. We find 188 Lyα emitters over 7.3 × 105 Mpc3, but also find significant numbers of other line-emitting sources corresponding to He II, C III] and C IV emission lines. These sources are important contaminants, and we carefully remove them, unlike most previous studies. We find that the Lyα luminosity function at z = 2.23 is very well described by a Schechter function up to LLy α ≈ 1043 erg s−1 with L∗=1042.59+0.16−0.08 erg s−1, ϕ∗=10−3.09+0.14−0.34 Mpc−3 and α = −1.75 ± 0.25. Above LLy α ≈ 1043 erg s−1, the Lyα luminosity function becomes power-law like, driven by X-ray AGN. We find that Lyα-selected emitters have a high escape fraction of 37 ± 7 per cent, anticorrelated with Lyα luminosity and correlated with Lyα equivalent width. Lyα emitters have ubiquitous large (≈40 kpc) Lyα haloes, ∼2 times larger than their Hα extents. By directly comparing our Lyα and Hα luminosity functions, we find that the global/overall escape fraction of Lyα photons (within a 13 kpc radius) from the full population of star-forming galaxies is 5.1 ± 0.2 per cent at the peak of the star formation history. An extra 3.3 ± 0.3 per cent of Lyα photons likely still escape, but at larger radii.},
  author       = {Sobral, David and Matthee, Jorryt J and Best, Philip and Stroe, Andra and Röttgering, Huub and Oteo, Iván and Smail, Ian and Morabito, Leah and Paulino-Afonso, Ana},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: haloes, galaxies: high-redshift, galaxies: luminosity function, mass function, galaxies: statistics, cosmology: observations},
  number       = {1},
  pages        = {1242--1258},
  publisher    = {Oxford University Press},
  title        = {{The CALYMHA survey: Lyα luminosity function and global escape fraction of Lyα photons at z = 2.23}},
  doi          = {10.1093/mnras/stw3090},
  volume       = {466},
  year         = {2017},
}

@article{11564,
  abstract     = {We study the production rate of ionizing photons of a sample of 588 Hα emitters (HAEs) and 160 Lyman-α emitters (LAEs) at z = 2.2 in the COSMOS field in order to assess the implied emissivity from galaxies, based on their ultraviolet (UV) luminosity. By exploring the rest-frame Lyman Continuum (LyC) with GALEX/NUV data, we find fesc < 2.8 (6.4) per cent through median (mean) stacking. By combining the Hα luminosity density with intergalactic medium emissivity measurements from absorption studies, we find a globally averaged 〈fesc〉 of 5.9+14.5−4.2 per cent at z = 2.2 if we assume HAEs are the only source of ionizing photons. We find similarly low values of the global 〈fesc〉 at z ≈ 3–5, also ruling out a high 〈fesc〉 at z < 5. These low escape fractions allow us to measure ξion, the number of produced ionizing photons per unit UV luminosity, and investigate how this depends on galaxy properties. We find a typical ξion ≈ 1024.77 ± 0.04 Hz erg−1 for HAEs and ξion ≈ 1025.14 ± 0.09 Hz erg−1 for LAEs. LAEs and low-mass HAEs at z = 2.2 show similar values of ξion as typically assumed in the reionization era, while the typical HAE is three times less ionizing. Due to an increasing ξion with increasing EW(Hα), ξion likely increases with redshift. This evolution alone is fully in line with the observed evolution of ξion between z ≈ 2 and 5, indicating a typical value of ξion ≈ 1025.4 Hz erg−1 in the reionization era.},
  author       = {Matthee, Jorryt J and Sobral, David and Best, Philip and Khostovan, Ali Ahmad and Oteo, Iván and Bouwens, Rychard and Röttgering, Huub},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: high-redshift, cosmology: observations, dark ages, reionization, first stars},
  number       = {3},
  pages        = {3637--3655},
  publisher    = {Oxford University Press},
  title        = {{The production and escape of Lyman-Continuum radiation from star-forming galaxies at z ∼ 2 and their redshift evolution}},
  doi          = {10.1093/mnras/stw2973},
  volume       = {465},
  year         = {2017},
}

@article{11565,
  abstract     = {We use the hydrodynamical EAGLE simulation to study the magnitude and origin of the scatter in the stellar mass–halo mass relation for central galaxies. We separate cause and effect by correlating stellar masses in the baryonic simulation with halo properties in a matched dark matter only (DMO) simulation. The scatter in stellar mass increases with redshift and decreases with halo mass. At z = 0.1, it declines from 0.25 dex at M200, DMO ≈ 1011 M⊙ to 0.12 dex at M200, DMO ≈ 1013 M⊙, but the trend is weak above 1012 M⊙. For M200, DMO < 1012.5 M⊙ up to 0.04 dex of the scatter is due to scatter in the halo concentration. At fixed halo mass, a larger stellar mass corresponds to a more concentrated halo. This is likely because higher concentrations imply earlier formation times and hence more time for accretion and star formation, and/or because feedback is less efficient in haloes with higher binding energies. The maximum circular velocity, Vmax, DMO, and binding energy are therefore more fundamental properties than halo mass, meaning that they are more accurate predictors of stellar mass, and we provide fitting formulae for their relations with stellar mass. However, concentration alone cannot explain the total scatter in the Mstar−M200,DMO relation, and it does not explain the scatter in Mstar–Vmax, DMO. Halo spin, sphericity, triaxiality, substructure and environment are also not responsible for the remaining scatter, which thus could be due to more complex halo properties or non-linear/stochastic baryonic effects.},
  author       = {Matthee, Jorryt J and Schaye, Joop and Crain, Robert A. and Schaller, Matthieu and Bower, Richard and Theuns, Tom},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: formation, galaxies: haloes, cosmology: theory},
  number       = {2},
  pages        = {2381--2396},
  publisher    = {Oxford University Press},
  title        = {{The origin of scatter in the stellar mass–halo mass relation of central galaxies in the EAGLE simulation}},
  doi          = {10.1093/mnras/stw2884},
  volume       = {465},
  year         = {2017},
}

@article{11566,
  abstract     = {While traditionally associated with active galactic nuclei (AGN), the properties of the C II] (λ = 2326 Å), C III] (λ, λ = 1907, 1909 Å) and C IV (λ, λ = 1549, 1551 Å) emission lines are still uncertain as large, unbiased samples of sources are scarce. We present the first blind, statistical study of C II], C III] and C IV emitters at z ∼ 0.68, 1.05, 1.53, respectively, uniformly selected down to a flux limit of ∼4 × 10−17 erg s−1 cm−1 through a narrow-band survey covering an area of ∼1.4 deg2 over COSMOS and UDS. We detect 16 C II], 35 C III] and 17 C IV emitters, whose nature we investigate using optical colours as well as Hubble Space Telescope (HST), X-ray, radio and far-infrared data. We find that z ∼ 0.7 C II] emitters are consistent with a mixture of blue (UV slope β = −2.0 ± 0.4) star-forming (SF) galaxies with discy HST structure and AGN with Seyfert-like morphologies. Bright C II] emitters have individual X-ray detections as well as high average black hole accretion rates (BHARs) of ∼0.1 M⊙ yr−1. C III] emitters at z ∼ 1.05 trace a general population of SF galaxies, with β = −0.8 ± 1.1, a variety of optical morphologies, including isolated and interacting galaxies and low BHAR (<0.02 M⊙ yr−1). Our C IV emitters at z ∼ 1.5 are consistent with young, blue quasars (β ∼ −1.9) with point-like optical morphologies, bright X-ray counterparts and large BHAR (0.8  M⊙ yr−1). We also find some surprising C II], C III] and C IV emitters with rest-frame equivalent widths (EWs) that could be as large as 50–100 Å. AGN or spatial offsets between the UV continuum stellar disc and the line-emitting regions may explain the large EW. These bright C II], C III] and C IV emitters are ideal candidates for spectroscopic follow-up to fully unveil their nature.},
  author       = {Stroe, Andra and Sobral, David and Matthee, Jorryt J and Calhau, João and Oteo, Ivan},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: active, galaxies: high-redshift, quasars: emission lines, galaxies: star formation, cosmology: observations},
  number       = {3},
  pages        = {2558--2574},
  publisher    = {Oxford University Press},
  title        = {{A 1.4 deg2 blind survey for C II], C III] and C IV at z ∼ 0.7–1.5 – I. Nature, morphologies and equivalent widths }},
  doi          = {10.1093/mnras/stx1712},
  volume       = {471},
  year         = {2017},
}

@article{11567,
  abstract     = {Recently, the C III] and C IV emission lines have been observed in galaxies in the early Universe (z > 5), providing new ways to measure their redshift and study their stellar populations and active galactic nuclei (AGN). We explore the first blind C II], C III] and C IV survey (z ∼ 0.68, 1.05, 1.53, respectively) presented in Stroe et al. (2017). We derive luminosity functions (LF) and study properties of C II], C III] and C IV line emitters through comparisons to the LFs of H α and Ly α emitters, UV selected star-forming (SF) galaxies and quasars at similar redshifts. The C II] LF at z ∼ 0.68 is equally well described by a Schechter or a power-law LF, characteristic of a mixture of SF and AGN activity. The C III] LF (z ∼ 1.05) is consistent to a scaled down version of the Schechter H α and Ly α LF at their redshift, indicating a SF origin. In stark contrast, the C IV LF at z ∼ 1.53 is well fit by a power-law, quasar-like LF. We find that the brightest UV sources (MUV < −22) will universally have C III] and C IV emission. However, on average, C III] and C IV are not as abundant as H α or Ly α emitters at the same redshift, with cosmic average ratios of ∼0.02–0.06 to H α and ∼0.01–0.1 to intrinsic Ly α. We predict that the C III] and C IV lines can only be truly competitive in confirming high-redshift candidates when the hosts are intrinsically bright and the effective Ly α escape fraction is below 1 per cent. While C III] and C IV were proposed as good tracers of young, relatively low-metallicity galaxies typical of the early Universe, we find that, at least at z ∼ 1.5, C IV is exclusively hosted by AGN/quasars, especially at large line equivalent widths.},
  author       = {Stroe, Andra and Sobral, David and Matthee, Jorryt J and Calhau, João and Oteo, Ivan},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: active, galaxies: high redshift, galaxies: luminosity function, mass function, quasars: emission lines, star formation, cosmology: observations},
  number       = {3},
  pages        = {2575--2586},
  publisher    = {Oxford University Press},
  title        = {{A 1.4 deg2 blind survey for C II], C III] and C IV at z ∼ 0.7–1.5 – II. Luminosity functions and cosmic average line ratios}},
  doi          = {10.1093/mnras/stx1713},
  volume       = {471},
  year         = {2017},
}

@article{11572,
  abstract     = {We present spectroscopic follow-up of candidate luminous Ly α emitters (LAEs) at z = 5.7–6.6 in the SA22 field with VLT/X-SHOOTER. We confirm two new luminous LAEs at z = 5.676 (SR6) and z = 6.532 (VR7), and also present HST follow-up of both sources. These sources have luminosities LLy α ≈ 3 × 1043 erg s−1, very high rest-frame equivalent widths of EW0 ≳ 200 Å and narrow Ly α lines (200–340 km s−1). VR7 is the most UV-luminous LAE at z > 6.5, with M1500 = −22.5, even brighter in the UV than CR7. Besides Ly α, we do not detect any other rest-frame UV lines in the spectra of SR6 and VR7, and argue that rest-frame UV lines are easier to observe in bright galaxies with low Ly α equivalent widths. We confirm that Ly α line widths increase with Ly α luminosity at z = 5.7, while there are indications that Ly α lines of faint LAEs become broader at z = 6.6, potentially due to reionization. We find a large spread of up to 3 dex in UV luminosity for >L⋆ LAEs, but find that the Ly α luminosity of the brightest LAEs is strongly related to UV luminosity at z = 6.6. Under basic assumptions, we find that several LAEs at z ≈ 6–7 have Ly α escape fractions ≳ 100  per cent, indicating bursty star formation histories, alternative Ly α production mechanisms, or dust attenuating Ly α emission differently than UV emission. Finally, we present a method to compute ξion, the production efficiency of ionizing photons, and find that LAEs at z ≈ 6–7 have high values of log10(ξion/Hz erg−1) ≈ 25.51 ± 0.09 that may alleviate the need for high Lyman-Continuum escape fractions required for reionization.},
  author       = {Matthee, Jorryt J and Sobral, David and Darvish, Behnam and Santos, Sérgio and Mobasher, Bahram and Paulino-Afonso, Ana and Röttgering, Huub and Alegre, Lara},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution – galaxies: high-redshift, dark ages, reionization, first stars, cosmology: observations},
  number       = {1},
  pages        = {772--787},
  publisher    = {Oxford University Press},
  title        = {{Spectroscopic properties of luminous Ly α emitters at z ≈ 6–7 and comparison to the Lyman-break population}},
  doi          = {10.1093/mnras/stx2061},
  volume       = {472},
  year         = {2017},
}

@article{11573,
  abstract     = {We present dynamical measurements from the KMOS (K-band multi-object spectrograph) Deep Survey (KDS), which comprises 77 typical star-forming galaxies at z ≃ 3.5 in the mass range 9.0 < log (M⋆/M⊙) < 10.5. These measurements constrain the internal dynamics, the intrinsic velocity dispersions (σint) and rotation velocities (VC) of galaxies in the high-redshift Universe. The mean velocity dispersion of the galaxies in our sample is σint=70.8+3.3−3.1kms−1⁠, revealing that the increasing average σint with increasing redshift, reported for z ≲ 2, continues out to z ≃ 3.5. Only 36 ± 8 per cent of our galaxies are rotation-dominated (VC/σint > 1), with the sample average VC/σint value much smaller than at lower redshift. After carefully selecting comparable star-forming samples at multiple epochs, we find that the rotation-dominated fraction evolves with redshift with a z−0.2 dependence. The rotation-dominated KDS galaxies show no clear offset from the local rotation velocity–stellar mass (i.e. VC–M⋆) relation, although a smaller fraction of the galaxies are on the relation due to the increase in the dispersion-dominated fraction. These observations are consistent with a simple equilibrium model picture, in which random motions are boosted in high-redshift galaxies by a combination of the increasing gas fractions, accretion efficiency, specific star formation rate and stellar feedback and which may provide significant pressure support against gravity on the galactic disc scale.},
  author       = {Turner, O. J. and Cirasuolo, M. and Harrison, C. M. and McLure, R. J. and Dunlop, J. S. and Swinbank, A. M. and Johnson, H. L. and Sobral, D. and Matthee, Jorryt J and Sharples, R. M.},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: high-redshift, galaxies: kinematics and dynamics},
  number       = {2},
  pages        = {1280--1320},
  publisher    = {Oxford University Press},
  title        = {{The KMOS Deep Survey (KDS) – I. Dynamical measurements of typical star-forming galaxies at z ≃ 3.5}},
  doi          = {10.1093/mnras/stx1366},
  volume       = {471},
  year         = {2017},
}

@article{13476,
  abstract     = {Understanding ionizing fluxes of stellar populations is crucial for various astrophysical problems including the epoch of reionization. Short-lived massive stars are generally considered as the main stellar sources. We examine the potential role of less massive stars that lose their envelope through interaction with a binary companion. Here, we focus on the role of metallicity (Z). For this purpose we used the evolutionary code MESA and created tailored atmosphere models with the radiative transfer code CMFGEN. We show that typical progenitors, with initial masses of 12 M⊙, produce hot and compact stars (~ 4 M⊙, 60–80 kK, ~1 R⊙). These stripped stars copiously produce ionizing photons, emitting 60–85% and 30–60% of their energy as HI and HeI ionizing radiation, for Z = 0.0001–0.02, respectively. Their output is comparable to what massive stars emit during their Wolf-Rayet phase, if we account for their longer lifetimes and the favorable slope of the initial mass function. Their relative importance for reionization may be further favored since they emit their photons with a time delay (~ 20 Myr after birth in our fiducial model). This allows time for the dispersal of the birth clouds, allowing the ionizing photons to escape into the intergalactic medium. At low Z, we find that Roche stripping fails to fully remove the H-rich envelope, because of the reduced opacity in the subsurface layers. This is in sharp contrast with the assumption of complete stripping that is made in rapid population synthesis simulations, which are widely used to simulate the binary progenitors of supernovae and gravitational waves. Finally, we discuss the urgency to increase the observed sample of stripped stars to test these models and we discuss how our predictions can help to design efficient observational campaigns.},
  author       = {Götberg, Ylva Louise Linsdotter and de Mink, S. E. and Groh, J. H.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Ionizing spectra of stars that lose their envelope through interaction with a binary companion: Role of metallicity}},
  doi          = {10.1051/0004-6361/201730472},
  volume       = {608},
  year         = {2017},
}

@article{13477,
  abstract     = {Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, 15+9-8%, of core-collapse supernovae are “late”, that is, they occur 50–200 Myr after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass (4–8 M⊙). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by 14+15-14% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that φ Persei is a likely progenitor and that eccentric neutron star – white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.},
  author       = {Zapartas, E. and de Mink, S. E. and Izzard, R. G. and Yoon, S.-C. and Badenes, C. and Götberg, Ylva Louise Linsdotter and de Koter, A. and Neijssel, C. J. and Renzo, M. and Schootemeijer, A. and Shrotriya, T. S.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {A&A},
  publisher    = {EDP Sciences},
  title        = {{Delay-time distribution of core-collapse supernovae with late events resulting from binary interaction}},
  doi          = {10.1051/0004-6361/201629685},
  volume       = {601},
  year         = {2017},
}

