@article{13466,
  abstract     = {Context. A majority of massive stars are part of binary systems, a large fraction of which will inevitably interact during their lives. Binary-interaction products (BiPs), that is, stars affected by such interaction, are expected to be commonly present in stellar populations. BiPs are thus a crucial ingredient in the understanding of stellar evolution.
Aims. We aim to identify and characterize a statistically significant sample of BiPs by studying clusters of 10 − 40 Myr, an age at which binary population models predict the abundance of BiPs to be highest. One example of such a cluster is NGC 330 in the Small Magellanic Cloud.
Methods. Using MUSE WFM-AO observations of NGC 330, we resolved the dense cluster core for the first time and were able to extract spectra of its entire massive star population. We developed an automated spectral classification scheme based on the equivalent widths of spectral lines in the red part of the spectrum.
Results. We characterize the massive star content of the core of NGC 330, which contains more than 200 B stars, 2 O stars, 6 A-type supergiants, and 11 red supergiants. We find a lower limit on the Be star fraction of 32 ± 3% in the whole sample. It increases to at least 46 ± 10% when we only consider stars brighter than V = 17 mag. We estimate an age of the cluster core between 35 and 40 Myr and a total cluster mass of 88−18+17 × 103 M⊙.
Conclusions. We find that the population in the cluster core is different than the population in the outskirts: while the stellar content in the core appears to be older than the stars in the outskirts, the Be star fraction and the observed binary fraction are significantly higher. Furthermore, we detect several BiP candidates that will be subject of future studies.},
  author       = {Bodensteiner, J. and Sana, H. and Mahy, L. and Patrick, L. R. and de Koter, A. and de Mink, S. E. and Evans, C. J. and Götberg, Ylva Louise Linsdotter and Langer, N. and Lennon, D. J. and Schneider, F. R. N. and Tramper, F.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {stars: massive / stars: emission-line / Be / binaries: spectroscopic / blue stragglers / Magellanic Clouds},
  publisher    = {EDP Sciences},
  title        = {{The young massive SMC cluster NGC 330 seen by MUSE}},
  doi          = {10.1051/0004-6361/201936743},
  volume       = {634},
  year         = {2020},
}

@article{13467,
  abstract     = {Massive stars are often found in binary systems, and it has been argued that binary products boost the ionizing radiation of stellar populations. Accurate predictions for binary products are needed to understand and quantify their contribution to cosmic reionization. We investigate the contribution of stars stripped in binaries because (1) they are, arguably, the best-understood products of binary evolution, (2) we recently produced the first radiative transfer calculations for the atmospheres of these stripped stars that predict their ionizing spectra, and (3) they are very promising sources because they boost the ionizing emission of stellar populations at late times. This allows stellar feedback to clear the surroundings such that a higher fraction of their photons can escape and ionize the intergalactic medium. Combining our detailed predictions for the ionizing spectra with a simple cosmic reionization model, we estimate that stripped stars contributed tens of percent of the photons that caused cosmic reionization of hydrogen, depending on the assumed escape fractions. More importantly, stripped stars harden the ionizing emission. We estimate that the spectral index for the ionizing part of the spectrum can increase to −1 compared to ≲ − 2 for single stars. At high redshift, stripped stars and massive single stars combined dominate the He II-ionizing emission, but we expect that active galactic nuclei drive cosmic helium reionization. Further observational consequences we expect are (1) high ionization states for the intergalactic gas surrounding stellar systems, such as C IV and Si IV, and (2) additional heating of the intergalactic medium of up to a few thousand Kelvin. Quantifying these warrants the inclusion of accurate models for stripped stars and other binary products in full cosmological simulations.},
  author       = {Götberg, Ylva Louise Linsdotter and de Mink, S. E. and McQuinn, M. and Zapartas, E. and Groh, J. H. and Norman, C.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Contribution from stars stripped in binaries to cosmic reionization of hydrogen and helium}},
  doi          = {10.1051/0004-6361/201936669},
  volume       = {634},
  year         = {2020},
}

@article{11499,
  abstract     = {Deep optical spectroscopic surveys of galaxies provide a unique opportunity to investigate rest-frame ultra-violet (UV) emission line properties of galaxies at z ∼ 2 − 4.5. Here we combine VLT/MUSE Guaranteed Time Observations of the Hubble Deep Field South, Ultra Deep Field, COSMOS, and several quasar fields with other publicly available data from VLT/VIMOS and VLT/FORS2 to construct a catalogue of He II λ1640 emitters at z ≳ 2. The deepest areas of our MUSE pointings reach a 3σ line flux limit of 3.1 × 10−19 erg s−1 cm−2. After discarding broad-line active galactic nuclei, we find 13 He II λ1640 detections from MUSE with a median MUV = −20.1 and 21 tentative He II λ1640 detections from other public surveys. Excluding Lyα, all except two galaxies in our sample show at least one other rest-UV emission line, with C III] λ1907, λ1909 being the most prominent. We use multi-wavelength data available in the Hubble legacy fields to derive basic galaxy properties of our sample through spectral energy distribution fitting techniques. Taking advantage of the high-quality spectra obtained by MUSE (∼10 − 30 h of exposure time per pointing), we use photo-ionisation models to study the rest-UV emission line diagnostics of the He II λ1640 emitters. Line ratios of our sample can be reproduced by moderately sub-solar photo-ionisation models, however, we find that including effects of binary stars lead to degeneracies in most free parameters. Even after considering extra ionising photons produced by extreme sub-solar metallicity binary stellar models, photo-ionisation models are unable to reproduce rest-frame He II λ1640 equivalent widths (∼0.2 − 10 Å), thus additional mechanisms are necessary in models to match the observed He II λ1640 properties.},
  author       = {Nanayakkara, Themiya and Brinchmann, Jarle and Boogaard, Leindert and Bouwens, Rychard and Cantalupo, Sebastiano and Feltre, Anna and Kollatschny, Wolfram and Marino, Raffaella Anna and Maseda, Michael and Matthee, Jorryt J and Paalvast, Mieke and Richard, Johan and Verhamme, Anne},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: ISM / galaxies: star formation / galaxies: evolution / galaxies: high-redshift},
  publisher    = {EDP Sciences},
  title        = {{Exploring He II λ1640 emission line properties at z ∼2−4}},
  doi          = {10.1051/0004-6361/201834565},
  volume       = {648},
  year         = {2019},
}

@article{11505,
  abstract     = {Contact. This paper presents the results obtained with the Multi-Unit Spectroscopic Explorer (MUSE) at the ESO Very Large Telescope on the faint end of the Lyman-alpha luminosity function (LF) based on deep observations of four lensing clusters. The goal of our project is to set strong constraints on the relative contribution of the Lyman-alpha emitter (LAE) population to cosmic reionization.

Aims. The precise aim of the present study is to further constrain the abundance of LAEs by taking advantage of the magnification provided by lensing clusters to build a blindly selected sample of galaxies which is less biased than current blank field samples in redshift and luminosity. By construction, this sample of LAEs is complementary to those built from deep blank fields, whether observed by MUSE or by other facilities, and makes it possible to determine the shape of the LF at fainter levels, as well as its evolution with redshift.

Methods. We selected a sample of 156 LAEs with redshifts between 2.9 ≤ z ≤ 6.7 and magnification-corrected luminosities in the range 39 ≲ log LLyα [erg s−1] ≲43. To properly take into account the individual differences in detection conditions between the LAEs when computing the LF, including lensing configurations, and spatial and spectral morphologies, the non-parametric 1/Vmax method was adopted. The price to pay to benefit from magnification is a reduction of the effective volume of the survey, together with a more complex analysis procedure to properly determine the effective volume Vmax for each galaxy. In this paper we present a complete procedure for the determination of the LF based on IFU detections in lensing clusters. This procedure, including some new methods for masking, effective volume integration and (individual) completeness determinations, has been fully automated when possible, and it can be easily generalized to the analysis of IFU observations in blank fields.

Results. As a result of this analysis, the Lyman-alpha LF has been obtained in four different redshift bins: 2.9 <  z <  6, 7, 2.9 <  z <  4.0, 4.0 <  z <  5.0, and 5.0 <  z <  6.7 with constraints down to log LLyα = 40.5. From our data only, no significant evolution of LF mean slope can be found. When performing a Schechter analysis also including data from the literature to complete the present sample towards the brightest luminosities, a steep faint end slope was measured varying from α = −1.69−0.08+0.08 to α = −1.87−0.12+0.12 between the lowest and the highest redshift bins.

Conclusions. The contribution of the LAE population to the star formation rate density at z ∼ 6 is ≲50% depending on the luminosity limit considered, which is of the same order as the Lyman-break galaxy (LBG) contribution. The evolution of the LAE contribution with redshift depends on the assumed escape fraction of Lyman-alpha photons, and appears to slightly increase with increasing redshift when this fraction is conservatively set to one. Depending on the intersection between the LAE/LBG populations, the contribution of the observed galaxies to the ionizing flux may suffice to keep the universe ionized at z ∼ 6.},
  author       = {de La Vieuville, G. and Bina, D. and Pello, R. and Mahler, G. and Richard, J. and Drake, A. B. and Herenz, E. C. and Bauer, F. E. and Clément, B. and Lagattuta, D. and Laporte, N. and Martinez, J. and Patrício, V. and Wisotzki, L. and Zabl, J. and Bouwens, R. J. and Contini, T. and Garel, T. and Guiderdoni, B. and Marino, R. A. and Maseda, M. V. and Matthee, Jorryt J and Schaye, J. and Soucail, G.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, gravitational lensing: strong / galaxies: high-redshift / dark ages, reionization, first stars / galaxies: clusters: general / galaxies: luminosity function, mass function},
  publisher    = {EDP Sciences},
  title        = {{Faint end of the z ∼ 3–7 luminosity function of Lyman-alpha emitters behind lensing clusters observed with MUSE}},
  doi          = {10.1051/0004-6361/201834471},
  volume       = {628},
  year         = {2019},
}

@article{11507,
  abstract     = {Lyman-α (Lyα) is intrinsically the brightest line emitted from active galaxies. While it originates from many physical processes, for star-forming galaxies the intrinsic Lyα luminosity is a direct tracer of the Lyman-continuum (LyC) radiation produced by the most massive O- and early-type B-stars (M⋆ ≳ 10 M⊙) with lifetimes of a few Myrs. As such, Lyα luminosity should be an excellent instantaneous star formation rate (SFR) indicator. However, its resonant nature and susceptibility to dust as a rest-frame UV photon makes Lyα very hard to interpret due to the uncertain Lyα escape fraction, fesc, Lyα. Here we explore results from the CAlibrating LYMan-α with Hα (CALYMHA) survey at z = 2.2, follow-up of Lyα emitters (LAEs) at z = 2.2 − 2.6 and a z ∼ 0−0.3 compilation of LAEs to directly measure fesc, Lyα with Hα. We derive a simple empirical relation that robustly retrieves fesc, Lyα as a function of Lyα rest-frame EW (EW0): fesc,Lyα = 0.0048 EW0[Å] ± 0.05 and we show that it constrains a well-defined anti-correlation between ionisation efficiency (ξion) and dust extinction in LAEs. Observed Lyα luminosities and EW0 are easy measurable quantities at high redshift, thus making our relation a practical tool to estimate intrinsic Lyα and LyC luminosities under well controlled and simple assumptions. Our results allow observed Lyα luminosities to be used to compute SFRs for LAEs at z ∼ 0−2.6 within ±0.2 dex of the Hα dust corrected SFRs. We apply our empirical SFR(Lyα,EW0) calibration to several sources at z ≥ 2.6 to find that star-forming LAEs have SFRs typically ranging from 0.1 to 20 M⊙ yr−1 and that our calibration might be even applicable for the most luminous LAEs within the epoch of re-ionisation. Our results imply high ionisation efficiencies (log10[ξion/Hz erg−1] = 25.4−25.6) and low dust content in LAEs across cosmic time, and will be easily tested with future observations with JWST which can obtain Hα and Hβ measurements for high-redshift LAEs.},
  author       = {Sobral, David and Matthee, Jorryt J},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / galaxies: star formation / galaxies: statistics / galaxies: evolution / galaxies: formation / galaxies: ISM},
  publisher    = {EDP Sciences},
  title        = {{Predicting Lyα escape fractions with a simple observable: Lyα in emission as an empirically calibrated star formation rate indicator}},
  doi          = {10.1051/0004-6361/201833075},
  volume       = {623},
  year         = {2019},
}

@article{11614,
  abstract     = {The NASA Transiting Exoplanet Survey Satellite (TESS) is about to provide full-frame images of almost the entire sky. The amount of stellar data to be analysed represents hundreds of millions stars, which is several orders of magnitude more than the number of stars observed by the Convection, Rotation and planetary Transits satellite (CoRoT), and NASA Kepler and K2 missions. We aim at automatically classifying the newly observed stars with near real-time algorithms to better guide the subsequent detailed studies. In this paper, we present a classification algorithm built to recognise solar-like pulsators among classical pulsators. This algorithm relies on the global amount of power contained in the power spectral density (PSD), also known as the flicker in spectral power density (FliPer). Because each type of pulsating star has a characteristic background or pulsation pattern, the shape of the PSD at different frequencies can be used to characterise the type of pulsating star. The FliPer classifier (FliPerClass) uses different FliPer parameters along with the effective temperature as input parameters to feed a ML algorithm in order to automatically classify the pulsating stars observed by TESS. Using noisy TESS-simulated data from the TESS Asteroseismic Science Consortium (TASC), we classify pulsators with a 98% accuracy. Among them, solar-like pulsating stars are recognised with a 99% accuracy, which is of great interest for a further seismic analysis of these stars, which are like our Sun. Similar results are obtained when we trained our classifier and applied it to 27-day subsets of real Kepler data. FliPerClass is part of the large TASC classification pipeline developed by the TESS Data for Asteroseismology (T’DA) classification working group.},
  author       = {Bugnet, Lisa Annabelle and García, R. A. and Mathur, S. and Davies, G. R. and Hall, O. J. and Lund, M. N. and Rendle, B. M.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Science},
  title        = {{FliPerClass: In search of solar-like pulsators among TESS targets}},
  doi          = {10.1051/0004-6361/201834780},
  volume       = {624},
  year         = {2019},
}

@article{13468,
  abstract     = {Hydrogen-rich supernovae, known as Type II (SNe II), are the most common class of explosions observed following the collapse of the core of massive stars. We used analytical estimates and population synthesis simulations to assess the fraction of SNe II progenitors that are expected to have exchanged mass with a companion prior to explosion. We estimate that 1/3 to 1/2 of SN II progenitors have a history of mass exchange with a binary companion before exploding. The dominant binary channels leading to SN II progenitors involve the merger of binary stars. Mergers are expected to produce a diversity of SN II progenitor characteristics, depending on the evolutionary timing and properties of the merger. Alternatively, SN II progenitors from interacting binaries may have accreted mass from their companion, and subsequently been ejected from the binary system after their companion exploded. We show that the overall fraction of SN II progenitors that are predicted to have experienced binary interaction is robust against the main physical uncertainties in our models. However, the relative importance of different binary evolutionary channels is affected by changing physical assumptions. We further discuss ways in which binarity might contribute to the observed diversity of SNe II by considering potential observational signatures arising from each binary channel. For supernovae which have a substantial H-rich envelope at explosion (i.e., excluding Type IIb SNe), a surviving non-compact companion would typically indicate that the supernova progenitor star was in a wide, non-interacting binary. We argue that a significant fraction of even Type II-P SNe are expected to have gained mass from a companion prior to explosion.},
  author       = {Zapartas, Emmanouil and de Mink, Selma E. and Justham, Stephen and Smith, Nathan and de Koter, Alex and Renzo, Mathieu and Arcavi, Iair and Farmer, Rob and Götberg, Ylva Louise Linsdotter and Toonen, Silvia},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{The diverse lives of progenitors of hydrogen-rich core-collapse supernovae: The role of binary interaction}},
  doi          = {10.1051/0004-6361/201935854},
  volume       = {631},
  year         = {2019},
}

@article{13469,
  abstract     = {Stars stripped of their envelopes from interaction with a binary companion emit a significant fraction of their radiation as ionizing photons. They are potentially important stellar sources of ionizing radiation, however, they are still often neglected in spectral synthesis simulations or simulations of stellar feedback. In anticipating the large datasets of galaxy spectra from the upcoming James Webb Space Telescope, we modeled the radiative contribution from stripped stars by using detailed evolutionary and spectral models. We estimated their impact on the integrated spectra and specifically on the emission rates of H I-, He I-, and He II-ionizing photons from stellar populations. We find that stripped stars have the largest impact on the ionizing spectrum of a population in which star formation halted several Myr ago. In such stellar populations, stripped stars dominate the emission of ionizing photons, mimicking a younger stellar population in which massive stars are still present. Our models also suggest that stripped stars have harder ionizing spectra than massive stars. The additional ionizing radiation, with which stripped stars contribute affects observable properties that are related to the emission of ionizing photons from stellar populations. In co-eval stellar populations, the ionizing radiation from stripped stars increases the ionization parameter and the production efficiency of hydrogen ionizing photons. They also cause high values for these parameters for about ten times longer than what is predicted for massive stars. The effect on properties related to non-ionizing wavelengths is less pronounced, such as on the ultraviolet continuum slope or stellar contribution to emission lines. However, the hard ionizing radiation from stripped stars likely introduces a characteristic ionization structure of the nebula, which leads to the emission of highly ionized elements such as O2+ and C3+. We, therefore, expect that the presence of stripped stars affects the location in the BPT diagram and the diagnostic ratio of O III to O II nebular emission lines. Our models are publicly available through CDS database and on the STARBURST99 website.},
  author       = {Götberg, Ylva Louise Linsdotter and de Mink, S. E. and Groh, J. H. and Leitherer, C. and Norman, C.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{The impact of stars stripped in binaries on the integrated spectra of stellar populations}},
  doi          = {10.1051/0004-6361/201834525},
  volume       = {629},
  year         = {2019},
}

@article{13470,
  abstract     = {Context. Massive Wolf–Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core collapse. It is not known whether core He-burning WR stars (classical WR; cWR) form predominantly through wind stripping (w-WR) or binary stripping (b-WR). Whereas spectroscopy of WR binaries has so-far largely been avoided because of its complexity, our study focuses on the 44 WR binaries and binary candidates of the Large Magellanic Cloud (LMC; metallicity Z ≈ 0.5 Z⊙), which were identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities.

Aims. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at subsolar metallicity and constraining the impact of binary interaction in forming these stars.

Methods. Spectroscopy was performed using the Potsdam Wolf–Rayet (PoWR) code and cross-correlation techniques. Disentanglement was performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status was interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically homogeneous evolution.

Results. Among our sample, 28/44 objects show composite spectra and are analyzed as such. An additional five targets show periodically moving WR primaries but no detected companions (SB1); two (BAT99 99 and 112) are potential WR + compact-object candidates owing to their high X-ray luminosities. We cannot confirm the binary nature of the remaining 11 candidates. About two-thirds of the WN components in binaries are identified as cWR, and one-third as hydrogen-burning WR stars. We establish metallicity-dependent mass-loss recipes, which broadly agree with those recently derived for single WN stars, and in which so-called WN3/O3 stars are clear outliers. We estimate that 45  ±  30% of the cWR stars in our sample have interacted with a companion via mass transfer. However, only ≈12  ±  7% of the cWR stars in our sample naively appear to have formed purely owing to stripping via a companion (12% b-WR). Assuming that apparently single WR stars truly formed as single stars, this comprises ≈4% of the whole LMC WN population, which is about ten times less than expected. No obvious differences in the properties of single and binary WN stars, whose luminosities extend down to log L ≈ 5.2 [L⊙], are apparent. With the exception of a few systems (BAT99 19, 49, and 103), the equatorial rotational velocities of the OB-type companions are moderate (veq ≲ 250 km s−1) and challenge standard formalisms of angular-momentum accretion. For most objects, chemically homogeneous evolution can be rejected for the secondary, but not for the WR progenitor.

Conclusions. No obvious dichotomy in the locations of apparently single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.},
  author       = {Shenar, T. and Sablowski, D. P. and Hainich, R. and Todt, H. and Moffat, A. F. J. and Oskinova, L. M. and Ramachandran, V. and Sana, H. and Sander, A. A. C. and Schnurr, O. and St-Louis, N. and Vanbeveren, D. and Götberg, Ylva Louise Linsdotter and Hamann, W.-R.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{The Wolf–Rayet binaries of the nitrogen sequence in the Large Magellanic Cloud}},
  doi          = {10.1051/0004-6361/201935684},
  volume       = {627},
  year         = {2019},
}

@article{13471,
  abstract     = {We perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (i) identify which predictions are robust against model uncertainties and assess their implications, (ii) investigate which physical processes leave a clear imprint and may therefore be constrained observationally, and (iii) provide a suite of publicly available model predictions to allow for the use of kinematic constraints from the Gaia mission. We find that 22+26−8% of all massive binary systems merge prior to the first core-collapse in the system. Of the remainder, 86+11−9% become unbound because of the core-collapse. Remarkably, this rarely produces runaway stars (observationally defined as stars with velocities above 30 km s−1). These are outnumbered by more than an order of magnitude by slower unbound companions, or “walkaway stars”. This is a robust outcome of our simulations and is due to the reversal of the mass ratio prior to the explosion and widening of the orbit, as we show analytically and numerically. For stars more massive than 15 M⊙, we estimate that 10+5−8% are walkaways and only 0.5+1.0−0.4% are runaways, nearly all of which have accreted mass from their companion. Our findings are consistent with earlier studies; however, the low runaway fraction we find is in tension with observed fractions of about 10%. Thus, astrometric data on presently single massive stars can potentially constrain the physics of massive binary evolution. Finally, we show that the high end of the mass distributions of runaway stars is very sensitive to the assumed black hole natal kicks, and we propose this as a potentially stringent test for the explosion mechanism. We also discuss companions remaining bound that can evolve into X-ray and gravitational wave sources.},
  author       = {Renzo, M. and Zapartas, E. and de Mink, S. E. and Götberg, Ylva Louise Linsdotter and Justham, S. and Farmer, R. J. and Izzard, R. G. and Toonen, S. and Sana, H.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Massive runaway and walkaway stars}},
  doi          = {10.1051/0004-6361/201833297},
  volume       = {624},
  year         = {2019},
}

@article{13472,
  abstract     = {Massive stars in binaries can give rise to extreme phenomena such as X-ray binaries and gravitational wave sources after one or both stars end their lives as core-collapse supernovae. Stars in close orbit around a stellar or compact companion are expected to explode as “stripped-envelope supernovae”, showing no (Type Ib/c) or little (Type IIb) signs of hydrogen in the spectra, because hydrogen-rich progenitors are too large to fit. The physical processes responsible for the stripping process and the fate of the companion are still very poorly understood. Aiming to find new clues, we investigate Cas A, which is a very young (∼340 yr) and near (∼3.4 kpc) remnant of a core-collapse supernova. Cas A has been subject to several searches for possible companions, all unsuccessfully. We present new measurements of the proper motions and photometry of stars in the vicinity based on deep HST ACS/WFC and WFC3-IR data. We identify stellar sources that are close enough in projection but using their proper motions we show that none are compatible with being at the location of center at the time of explosion, in agreement with earlier findings. Our photometric measurements allow us to place much deeper (order-of-magnitude) upper limits on the brightness of possible undetected companions. We systematically compare them with model predictions for a wide variety of scenarios. We can confidently rule out the presence of any stellar companion of any reasonable mass and age (main sequence, pre main sequence or stripped) ruling out what many considered to be likely evolutionary scenarios for Type IIb supernova (SN IIb). More exotic scenarios that predict the presence of a compact companion (white dwarf, neutron star or black hole) are still possible as well as scenarios where the progenitor of Cas A was single at the moment of explosion (either because it was truly single, or resulted from a binary that was disrupted, or from a binary merger). The presence of a compact companion would imply that Cas A is of interest to study exotic outcomes of binary evolution. The single-at-death solution would still require fine-tuning of the process that removed most of the envelope through a mass-loss mechanism yet to be identified. We discuss how future constraints from Gaia and even deeper photometric studies may help to place further constraints.},
  author       = {Kerzendorf, Wolfgang E. and Do, Tuan and de Mink, Selma E. and Götberg, Ylva Louise Linsdotter and Milisavljevic, Dan and Zapartas, Emmanouil and Renzo, Mathieu and Justham, Stephen and Podsiadlowski, Philipp and Fesen, Robert A.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{No surviving non-compact stellar companion to Cassiopeia A}},
  doi          = {10.1051/0004-6361/201732206},
  volume       = {623},
  year         = {2019},
}

@article{11508,
  abstract     = {Distant luminous Lyman-α emitters (LAEs) are excellent targets for spectroscopic observations of galaxies in the epoch of reionisation (EoR). We present deep high-resolution (R = 5000) VLT/X-shooter observations, along with an extensive collection of photometric data of COLA1, a proposed double peaked LAE at z = 6.6. We rule out the possibility that COLA1’s emission line is an [OII] doublet at z = 1.475 on the basis of i) the asymmetric red line-profile and flux ratio of the peaks (blue/red=0.31 ± 0.03) and ii) an unphysical [OII]/Hα ratio ([OII]/Hα >  22). We show that COLA1’s observed B-band flux is explained by a faint extended foreground LAE, for which we detect Lyα and [OIII] at z = 2.142. We thus conclude that COLA1 is a real double-peaked LAE at z = 6.593, the first discovered at z >  6. COLA1 is UV luminous (M1500 = −21.6 ± 0.3), has a high equivalent width (EW0,Lyα = 120−40+50 Å) and very compact Lyα emission (r50,Lyα = 0.33−0.04+0.07 kpc). Relatively weak inferred Hβ+[OIII] line-emission from Spitzer/IRAC indicates an extremely low metallicity of Z <  1/20 Z⊙ or reduced strength of nebular lines due to high escape of ionising photons. The small Lyα peak separation of 220 ± 20 km s−1 implies a low HI column density and an ionising photon escape fraction of ≈15 − 30%, providing the first direct evidence that such galaxies contribute actively to the reionisation of the Universe at z >  6. Based on simple estimates, we find that COLA1 could have provided just enough photons to reionise its own ≈0.3 pMpc (2.3 cMpc) bubble, allowing the blue Lyα line to be observed. However, we also discuss alternative scenarios explaining the detected double peaked nature of COLA1. Our results show that future high-resolution observations of statistical samples of double peaked LAEs at z >  5 are a promising probe of the occurrence of ionised regions around galaxies in the EoR.},
  author       = {Matthee, Jorryt J and Sobral, David and Gronke, Max and Paulino-Afonso, Ana and Stefanon, Mauro and Röttgering, Huub},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / galaxies: formation / dark ages / reionization / first stars / techniques: spectroscopic / intergalactic medium},
  publisher    = {EDP Sciences},
  title        = {{Confirmation of double peaked Lyα emission at z = 6.593: Witnessing a galaxy directly contributing to the reionisation of the universe}},
  doi          = {10.1051/0004-6361/201833528},
  volume       = {619},
  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{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{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{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{13478,
  abstract     = {HR 8799 is a star accompanied by four massive planets on wide orbits. The observed planetary configuration has been shown to be unstable on a timescale much shorter than the estimated age of the system (~30 Myr) unless the planets are locked into mean motion resonances. This condition is characterised by small-amplitude libration of one or more resonant angles that stabilise the system by preventing close encounters. We simulate planetary systems similar to the HR 8799 planetary system, exploring the parameter space in separation between the orbits, planetary masses and distance from the Sun to the star. We find systems that look like HR 8799 and remain stable for longer than the estimated age of HR 8799. None of our systems are forced into resonances. We find, with nominal masses (Mb = 5 MJup and Mc,d,e = 7 MJup) and in a narrow range of orbit separations, that 5 of 100 systems match the observations and lifetime. Considering a broad range of orbit separations, we find 12 of 900 similar systems. The systems survive significantly longer because of their slightly increased initial orbit separations compared to assuming circular orbits from the observed positions. A small increase in separation leads to a significant increase in survival time. The low eccentricity the orbits develop from gravitational interaction is enough for the planets to match the observations. With lower masses, but still comfortably within the estimated planet mass uncertainty, we find 18 of 100 matching and long-lived systems in a narrow orbital separation range. In the broad separation range, we find 82 of 900 matching systems. Our results imply that the planets in the HR 8799 system do not have to be in strong mean motion resonances. We also investigate the future of wide-orbit planetary systems using our HR 8799 analogues. We find that 80% of the systems have two planets left after strong planet-planet scattering and these are on eccentric orbits with semi-major axes of a1 ~ 10 AU and a2 ~ 30−1000 AU. We speculate that other wide-orbit planetary systems, such as AB Pic and HD 106906, are the remnants of HR 8799 analogues that underwent close encounters and dynamical instability.},
  author       = {Götberg, Ylva Louise Linsdotter and Davies, Melvyn B. and Mustill, Alexander J. and Johansen, Anders and Church, Ross P.},
  issn         = {1432-0746},
  journal      = {Astronomy & Astrophysics},
  keywords     = {Space and Planetary Science, Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Long-term stability of the HR 8799 planetary system without resonant lock}},
  doi          = {10.1051/0004-6361/201526309},
  volume       = {592},
  year         = {2016},
}