@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{15234,
  abstract     = {Using parallaxes from Gaia Data Release 2 (Gaia DR2), we estimate the distance to the globular clusters 47 Tuc and NGC 362, taking advantage of the background stars in the Small Magellanic Cloud and quasars to account for various parallax systematics. We found the parallax to be dependent on the Gaia DR2 G-band apparent magnitude for stars with 13 < G < 18, where brighter stars have a lower parallax zero point than fainter stars. The distance to 47 Tuc was found to be 4.45 ± 0.01 ± 0.12 kpc, and for NGC 362 8.54 ± 0.20 ± 0.44 kpc, with random and systematic errors listed, respectively. This is the first time a precise distance measurement directly using parallaxes has been determined for either of these two globular clusters.},
  author       = {Chen, Seery and Richer, Harvey and Caiazzo, Ilaria and Heyl, Jeremy},
  issn         = {1538-4357},
  journal      = {The Astrophysical Journal},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {2},
  publisher    = {American Astronomical Society},
  title        = {{Distances to the globular clusters 47 Tucanae and NGC 362 using Gaia DR2 parallaxes}},
  doi          = {10.3847/1538-4357/aae089},
  volume       = {867},
  year         = {2018},
}

@article{15236,
  abstract     = {Radio pulsars found in binary systems with short orbital periods are usually fast spinning as a consequence of recycling via mass transfer from their companion stars; this process is also thought to decrease the magnetic field of the neutron star being recycled. Here, we report on timing observations of the recently discovered binary PSR J1755−2550 and find that this pulsar is an exception: with a characteristic age of 2.1 Myr, it is relatively young; furthermore, with a spin period of 315 ms and a surface magnetic field strength at its poles of 0.88 × 1012 G, the pulsar shows no sign of having been recycled. Based on its timing and orbital characteristics, the pulsar either has a massive white dwarf (WD) or a neutron star (NS) companion. To distinguish between these two cases, we searched radio observations for a potential recycled pulsar companion and analysed archival optical data for a potential WD companion. Neither work returned conclusive detections. We apply population synthesis modelling and find that both solutions are roughly equally probable. Our population synthesis also predicts a minimum mass of 0.90 M⊙ for the companion star to PSR J1755−2550 and we simulate the systemic runaway velocities for the resulting WDNS systems which may merge and possibly produce Ca-rich supernovae. Whether PSR J1755−2550 hosts a WD or a NS companion star, it is certainly a member of a rare subpopulation of binary radio pulsars.},
  author       = {Ng, C and Kruckow, M U and Tauris, T M and Lyne, A G and Freire, P C C and Ridolfi, A and Caiazzo, Ilaria and Heyl, J and Kramer, M and Cameron, A D and Champion, D J and Stappers, B},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {4},
  pages        = {4315--4326},
  publisher    = {Oxford University Press},
  title        = {{PSR J1755−2550: A young radio pulsar with a massive, compact companion}},
  doi          = {10.1093/mnras/sty482},
  volume       = {476},
  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},
}

@article{15239,
  abstract     = {Using images from the Hubble Space Telescope Advanced Camera for Surveys, we measure the rate of cooling of white dwarfs in the globular cluster 47 Tucanae and compare it to modelled cooling curves. We examine the effects of the outer convective envelope reaching the nearly isothermal degenerate core and the release of latent heat during core crystallization on the white dwarf cooling rates. For white dwarfs typical of 47 Tuc, the onset of these effects occur at similar times. The latent heat released during crystallization is a small heat source. In contrast, the heat reservoir of the degenerate core is substantially larger. When the convective envelope reaches the nearly isothermal interior of the white dwarf, the star becomes brighter than it would be in the absence of this effect. Our modelled cooling curves that include this convective coupling closely match the observed luminosity function of the white dwarfs in 47 Tuc.},
  author       = {Obertas, Alysa and Caiazzo, Ilaria and Heyl, Jeremy and Richer, Harvey and Kalirai, Jason and Tremblay, Pier-Emmanuel},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {1},
  pages        = {677--682},
  publisher    = {Oxford University Press},
  title        = {{The onset of convective coupling and freezing in the white dwarfs of 47 Tucanae}},
  doi          = {10.1093/mnras/stx2759},
  volume       = {474},
  year         = {2017},
}

@article{15240,
  abstract     = {Multi-epoch observations with the Advanced Camera Survey and WFC3 on the Hubble Space Telescope provide a unique and comprehensive probe of stellar dynamics within 47 Tucanae. We confront analytic models of the globular cluster with the observed stellar proper motions that probe along the main sequence from just above 0.8–0.1M⊙ as well as white dwarfs younger than 1 Gyr. One field lies just beyond the half-light radius where dynamical models (e.g., lowered Maxwellian distributions) make robust predictions for the stellar proper motions. The observed proper motions in this outer field show evidence for anisotropy in the velocity distribution as well as skewness; the latter is evidence of rotation. The measured velocity dispersions and surface brightness distributions agree in detail with a rotating anisotropic model of the stellar distribution function with mild dependence of the proper-motion dispersion on mass. However, the best-fitting models underpredict the rotation and skewness of the stellar velocities. In the second field, centered on the core of the cluster, the mass segregation in proper motion is much stronger. Nevertheless the model developed in the outer field can be extended inward by taking this mass segregation into account in a heuristic fashion. The proper motions of the main-sequence stars yield a mass estimate of the cluster of 
 at a distance of 4.7 kpc. By comparing the proper motions of a sample of giant and subgiant stars with the observed radial velocities we estimate the distance to the cluster kinematically to be 4.29 ± 0.47 kpc.},
  author       = {Heyl, J. and Caiazzo, Ilaria and Richer, H. and Anderson, J. and Kalirai, J. and Parada, J.},
  issn         = {1538-4357},
  journal      = {The Astrophysical Journal},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {2},
  publisher    = {American Astronomical Society},
  title        = {{Deep HST imaging in 47 Tucanae: A global dynamical model}},
  doi          = {10.3847/1538-4357/aa974f},
  volume       = {850},
  year         = {2017},
}

@article{15241,
  abstract     = {We present a model to account for the observed debris discs around young white dwarfs and the presence of metal lines in their spectra. Stellar evolution models predict that the mass-loss on the AGB will be pulsed; furthermore, observations indicate that the bulk of the mass-loss occurs on the AGB. In this case, if the progenitors of the white dwarfs had remnants of planetary formation like the Sun’s Oort cloud or the Kuiper Belt and a planet lying within that cloud or nearby, we find that up to 2 per cent of the planetesimals will fall either into planet-crossing orbits or into chaotic regions after the mass-loss, depending on the location and mass of the planet (from Mars to Neptune). This yields a sufficient mass of comets that can be scattered towards the star, form a debris disc and pollute the atmosphere.},
  author       = {Caiazzo, Ilaria and Heyl, Jeremy S.},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {3},
  pages        = {2750--2759},
  publisher    = {Oxford University Press},
  title        = {{Polluting white dwarfs with perturbed exo-comets}},
  doi          = {10.1093/mnras/stx1036},
  volume       = {469},
  year         = {2017},
}

@article{15243,
  abstract     = {High-energy (>250 keV) emission has been detected persisting for several tens of seconds after the initial spike of magnetar giant flares (GFs). It has been conjectured that this emission might arise via inverse Compton scattering in a highly extended corona generated by super-Eddington outflows high up in the magnetosphere. In this paper, we undertake a detailed examination of this model. We investigate the properties of the required scatterers, and whether the mechanism is consistent with the degree of pulsed emission observed in the tail of the GF. We conclude that the mechanism is consistent with current data, although the origin of the scattering population remains an open question. We propose an alternative picture in which the emission is closer to that star and is dominated by synchrotron radiation. The Reuven Ramaty High Energy Solar Spectroscopic Imager observations of the 2004 December flare modestly favour this latter picture. We assess the prospects for the Fermi Gamma-ray Space Telescope to detect and characterize a similar high-energy component in a future GF. Such a detection should help to resolve some of the outstanding issues.},
  author       = {Elenbaas, C. and Huppenkothen, D. and Omand, C. and Watts, A. L. and Bissaldi, E. and Caiazzo, Ilaria and Heyl, J.},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  number       = {2},
  pages        = {1856--1872},
  publisher    = {Oxford University Press},
  title        = {{Magnetar giant flare high-energy emission}},
  doi          = {10.1093/mnras/stx1727},
  volume       = {471},
  year         = {2017},
}

@article{11574,
  abstract     = {We present new results from the widest narrow-band survey search for Lyα emitters at z = 5.7, just after reionization. We survey a total of 7 deg2 spread over the COSMOS, UDS and SA22 fields. We find over 11 000 line emitters, out of which 514 are robust Lyα candidates at z = 5.7 within a volume of 6.3 × 106 Mpc3. Our Lyα emitters span a wide range in Lyα luminosities, from faint to bright (LLyα ∼ 1042.5–44 erg s−1) and rest-frame equivalent widths (EW0 ∼ 25–1000 Å) in a single, homogeneous data set. By combining all our fields, we find that the faint end slope of the z = 5.7 Lyα luminosity function is very steep, with α=−2.3+0.4−0.3⁠. We also present an updated z = 6.6 Lyα luminosity function, based on comparable volumes and obtained with the same methods, which we directly compare with that at z = 5.7. We find a significant decline of the number density of faint Lyα emitters from z = 5.7 to 6.6 (by 0.5 ± 0.1 dex), but no evolution at the bright end/no evolution in L*. Faint Lyα emitters at z = 6.6 show much more extended haloes than those at z = 5.7, suggesting that neutral Hydrogen plays an important role, increasing the scattering and leading to observations missing faint Lyα emission within the epoch of reionization. Altogether, our results suggest that we are observing patchy reionization which happens first around the brightest Lyα emitters, allowing the number densities of those sources to remain unaffected by the increase of neutral Hydrogen fraction from z ∼ 5 to 7.},
  author       = {Santos, Sérgio and Sobral, David and Matthee, Jorryt J},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift, galaxies: luminosity function, mass function, cosmology: observations, dark ages, reionization, first stars},
  number       = {2},
  pages        = {1678--1691},
  publisher    = {Oxford University Press},
  title        = {{The Lyα luminosity function at z= 5.7–6.6 and the steep drop of the faint end: Implications for reionization}},
  doi          = {10.1093/mnras/stw2076},
  volume       = {463},
  year         = {2016},
}

@article{11575,
  abstract     = {We investigate correlations between different physical properties of star-forming galaxies in the ‘Evolution and Assembly of GaLaxies and their Environments’ (EAGLE) cosmological hydrodynamical simulation suite over the redshift range 0 ≤ z ≤ 4.5. A principal component analysis reveals that neutral gas fraction (fgas,neutral), stellar mass (Mstellar) and star formation rate (SFR) account for most of the variance seen in the population, with galaxies tracing a two-dimensional, nearly flat, surface in the three-dimensional space of fgas, neutral–Mstellar–SFR with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their fgas, neutral and SFR drop with redshift. The positions of galaxies along the plane are highly correlated with gas metallicity. The metallicity can therefore be robustly predicted from fgas, neutral, or from the Mstellar and SFR. We argue that the appearance of this ‘Fundamental Plane of star formation’ is a consequence of self-regulation, with the plane's curvature set by the dependence of the SFR on gas density and metallicity. We analyse a large compilation of observations spanning the redshift range 0 ≲ z ≲ 3, and find that such a plane is also present in the data. The properties of the observed Fundamental Plane of star formation are in good agreement with EAGLE's predictions.},
  author       = {Lagos, Claudia del P. and Theuns, Tom and Schaye, Joop and Furlong, Michelle and Bower, Richard G. and Schaller, Matthieu and Crain, Robert A. and Trayford, James W. and Matthee, Jorryt J},
  issn         = {1365-2966},
  journal      = {Monthly Notices of the Royal Astronomical Society},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics  stars: formation, ISM: evolution, galaxies: evolution, galaxies: formation, galaxies: ISM},
  number       = {3},
  pages        = {2632--2650},
  publisher    = {Oxford University Press},
  title        = {{The Fundamental Plane of star formation in galaxies revealed by the EAGLE hydrodynamical simulations}},
  doi          = {10.1093/mnras/stw717},
  volume       = {459},
  year         = {2016},
}

