@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{6756,
  abstract     = {We study the topology generated by the temperature fluctuations of the cosmic microwave background (CMB) radiation, as quantified by the number of components and holes, formally given by the Betti numbers, in the growing excursion sets. We compare CMB maps observed by the Planck satellite with a thousand simulated maps generated according to the ΛCDM paradigm with Gaussian distributed fluctuations. The comparison is multi-scale, being performed on a sequence of degraded maps with mean pixel separation ranging from 0.05 to 7.33°. The survey of the CMB over 𝕊2 is incomplete due to obfuscation effects by bright point sources and other extended foreground objects like our own galaxy. To deal with such situations, where analysis in the presence of “masks” is of importance, we introduce the concept of relative homology. The parametric χ2-test shows differences between observations and simulations, yielding p-values at percent to less than permil levels roughly between 2 and 7°, with the difference in the number of components and holes peaking at more than 3σ sporadically at these scales. The highest observed deviation between the observations and simulations for b0 and b1 is approximately between 3σ and 4σ at scales of 3–7°. There are reports of mildly unusual behaviour of the Euler characteristic at 3.66° in the literature, computed from independent measurements of the CMB temperature fluctuations by Planck’s predecessor, the Wilkinson Microwave Anisotropy Probe (WMAP) satellite. The mildly anomalous behaviour of the Euler characteristic is phenomenologically related to the strongly anomalous behaviour of components and holes, or the zeroth and first Betti numbers, respectively. Further, since these topological descriptors show consistent anomalous behaviour over independent measurements of Planck and WMAP, instrumental and systematic errors may be an unlikely source. These are also the scales at which the observed maps exhibit low variance compared to the simulations, and approximately the range of scales at which the power spectrum exhibits a dip with respect to the theoretical model. Non-parametric tests show even stronger differences at almost all scales. Crucially, Gaussian simulations based on power-spectrum matching the characteristics of the observed dipped power spectrum are not able to resolve the anomaly. Understanding the origin of the anomalies in the CMB, whether cosmological in nature or arising due to late-time effects, is an extremely challenging task. Regardless, beyond the trivial possibility that this may still be a manifestation of an extreme Gaussian case, these observations, along with the super-horizon scales involved, may motivate the study of primordial non-Gaussianity. Alternative scenarios worth exploring may be models with non-trivial topology, including topological defect models.},
  author       = {Pranav, Pratyush and Adler, Robert J. and Buchert, Thomas and Edelsbrunner, Herbert and Jones, Bernard J.T. and Schwartzman, Armin and Wagner, Hubert and Van De Weygaert, Rien},
  issn         = {1432-0746},
  journal      = {Astronomy and Astrophysics},
  publisher    = {EDP Sciences},
  title        = {{Unexpected topology of the temperature fluctuations in the cosmic microwave background}},
  doi          = {10.1051/0004-6361/201834916},
  volume       = {627},
  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},
}

