@article{11500, abstract = {We report the discovery of diffuse extended Lyα emission from redshift 3.1 to 4.5, tracing cosmic web filaments on scales of 2.5−4 cMpc. These structures have been observed in overdensities of Lyα emitters in the MUSE Extremely Deep Field, a 140 h deep MUSE observation located in the Hubble Ultra-Deep Field. Among the 22 overdense regions identified, five are likely to harbor very extended Lyα emission at high significance with an average surface brightness of 5 × 10−20 erg s−1 cm−2 arcsec−2. Remarkably, 70% of the total Lyα luminosity from these filaments comes from beyond the circumgalactic medium of any identified Lyα emitter. Fluorescent Lyα emission powered by the cosmic UV background can only account for less than 34% of this emission at z ≈ 3 and for not more than 10% at higher redshift. We find that the bulk of this diffuse emission can be reproduced by the unresolved Lyα emission of a large population of ultra low-luminosity Lyα emitters (< 1040 erg s−1), provided that the faint end of the Lyα luminosity function is steep (α ⪅ −1.8), it extends down to luminosities lower than 1038 − 1037 erg s−1, and the clustering of these Lyα emitters is significant (filling factor < 1/6). If these Lyα emitters are powered by star formation, then this implies their luminosity function needs to extend down to star formation rates < 10−4 M⊙ yr−1. These observations provide the first detection of the cosmic web in Lyα emission in typical filamentary environments and the first observational clue indicating the existence of a large population of ultra low-luminosity Lyα emitters at high redshift.}, author = {Bacon, R. and Mary, D. and Garel, T. and Blaizot, J. and Maseda, M. and Schaye, J. and Wisotzki, L. and Conseil, S. and Brinchmann, J. and Leclercq, F. and Abril-Melgarejo, V. and Boogaard, L. and Bouché, N. F. and Contini, T. and Feltre, A. and Guiderdoni, B. and Herenz, C. and Kollatschny, W. and Kusakabe, H. and Matthee, Jorryt J and Michel-Dansac, L. and Nanayakkara, T. and Richard, J. and Roth, M. and Schmidt, K. B. and Steinmetz, M. and Tresse, L. and Urrutia, T. and Verhamme, A. and Weilbacher, P. M. and Zabl, J. and Zoutendijk, S. L.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: high-redshift / galaxies: groups: general / cosmology: observations}, publisher = {EDP Sciences}, title = {{The MUSE Extremely Deep Field: The cosmic web in emission at high redshift}}, doi = {10.1051/0004-6361/202039887}, volume = {647}, year = {2021}, } @article{11605, abstract = {Context. The discovery of moderate differential rotation between the core and the envelope of evolved solar-like stars could be the signature of a strong magnetic field trapped inside the radiative interior. The population of intermediate-mass red giants presenting surprisingly low-amplitude mixed modes (i.e. oscillation modes that behave as acoustic modes in their external envelope and as gravity modes in their core) could also arise from the effect of an internal magnetic field. Indeed, stars more massive than about 1.1 solar masses are known to develop a convective core during their main sequence. The field generated by the dynamo triggered by this convection could be the progenitor of a strong fossil magnetic field trapped inside the core of the star for the remainder of its evolution. Aims. Observations of mixed modes can constitute an excellent probe of the deepest layers of evolved solar-like stars, and magnetic fields in those regions can impact their propagation. The magnetic perturbation on mixed modes may therefore be visible in asteroseismic data. To unravel which constraints can be obtained from observations, we theoretically investigate the effects of a plausible mixed axisymmetric magnetic field with various amplitudes on the mixed-mode frequencies of evolved solar-like stars. Methods. First-order frequency perturbations due to an axisymmetric magnetic field were computed for dipolar and quadrupolar mixed modes. These computations were carried out for a range of stellar ages, masses, and metallicities. Conclusions. We show that typical fossil-field strengths of 0.1 − 1 MG, consistent with the presence of a dynamo in the convective core during the main sequence, provoke significant asymmetries on mixed-mode frequency multiplets during the red giant branch. We provide constraints and methods for the detectability of such magnetic signatures. We show that these signatures may be detectable in asteroseismic data for field amplitudes small enough for the amplitude of the modes not to be affected by the conversion of gravity into Alfvén waves inside the magnetised interior. Finally, we infer an upper limit for the strength of the field and the associated lower limit for the timescale of its action in order to redistribute angular momentum in stellar interiors.}, author = {Bugnet, Lisa Annabelle and Prat, V. and Mathis, S. and Astoul, A. and Augustson, K. and García, R. A. and Mathur, S. and Amard, L. and Neiner, C.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, stars, oscillations / stars, magnetic field / stars, interiors / stars, evolution / stars, rotation}, publisher = {EDP Sciences}, title = {{Magnetic signatures on mixed-mode frequencies: I. An axisymmetric fossil field inside the core of red giants}}, doi = {10.1051/0004-6361/202039159}, volume = {650}, year = {2021}, } @article{11606, abstract = {Context. Our knowledge of the dynamics of stars has undergone a revolution through the simultaneous large amount of high-quality photometric observations collected by space-based asteroseismology and ground-based high-precision spectropolarimetry. They allowed us to probe the internal rotation of stars and their surface magnetism in the whole Hertzsprung-Russell diagram. However, new methods should still be developed to probe the deep magnetic fields in these stars. Aims. Our goal is to provide seismic diagnoses that allow us to probe the internal magnetism of stars. Methods. We focused on asymptotic low-frequency gravity modes and high-frequency acoustic modes. Using a first-order perturbative theory, we derived magnetic splittings of their frequencies as explicit functions of stellar parameters. Results. As in the case of rotation, we show that asymptotic gravity and acoustic modes can allow us to probe the different components of the magnetic field in the cavities in which they propagate. This again demonstrates the high potential of using mixed-modes when this is possible.}, author = {Mathis, S. and Bugnet, Lisa Annabelle and Prat, V. and Augustson, K. and Mathur, S. and Garcia, R. A.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, asteroseismology / waves / stars, magnetic field / stars, oscillations / methods, analytical}, publisher = {EDP Sciences}, title = {{Probing the internal magnetism of stars using asymptotic magneto-asteroseismology}}, doi = {10.1051/0004-6361/202039180}, volume = {647}, year = {2021}, } @article{11608, abstract = {In order to understand stellar evolution, it is crucial to efficiently determine stellar surface rotation periods. Indeed, while they are of great importance in stellar models, angular momentum transport processes inside stars are still poorly understood today. Surface rotation, which is linked to the age of the star, is one of the constraints needed to improve the way those processes are modelled. Statistics of the surface rotation periods for a large sample of stars of different spectral types are thus necessary. An efficient tool to automatically determine reliable rotation periods is needed when dealing with large samples of stellar photometric datasets. The objective of this work is to develop such a tool. For this purpose, machine learning classifiers constitute relevant bases to build our new methodology. Random forest learning abilities are exploited to automate the extraction of rotation periods in Kepler light curves. Rotation periods and complementary parameters are obtained via three different methods: a wavelet analysis, the autocorrelation function of the light curve, and the composite spectrum. We trained three different classifiers: one to detect if rotational modulations are present in the light curve, one to flag close binary or classical pulsators candidates that can bias our rotation period determination, and finally one classifier to provide the final rotation period. We tested our machine learning pipeline on 23 431 stars of the Kepler K and M dwarf reference rotation catalogue for which 60% of the stars have been visually inspected. For the sample of 21 707 stars where all the input parameters are provided to the algorithm, 94.2% of them are correctly classified (as rotating or not). Among the stars that have a rotation period in the reference catalogue, the machine learning provides a period that agrees within 10% of the reference value for 95.3% of the stars. Moreover, the yield of correct rotation periods is raised to 99.5% after visually inspecting 25.2% of the stars. Over the two main analysis steps, rotation classification and period selection, the pipeline yields a global agreement with the reference values of 92.1% and 96.9% before and after visual inspection. Random forest classifiers are efficient tools to determine reliable rotation periods in large samples of stars. The methodology presented here could be easily adapted to extract surface rotation periods for stars with different spectral types or observed by other instruments such as K2, TESS or by PLATO in the near future.}, author = {Breton, S. N. and Santos, A. R. G. and Bugnet, Lisa Annabelle and Mathur, S. and García, R. A. and Pallé, P. L.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, methods: data analysis / stars: solar-type / stars: activity / stars: rotation / starspots}, publisher = {EDP Sciences}, title = {{ROOSTER: A machine-learning analysis tool for Kepler stellar rotation periods}}, doi = {10.1051/0004-6361/202039947}, volume = {647}, year = {2021}, } @article{11609, abstract = {Context. Stellar interiors are the seat of efficient transport of angular momentum all along their evolution. In this context, understanding the dependence of the turbulent transport triggered by the instabilities of the vertical and horizontal shears of the differential rotation in stellar radiation zones as a function of their rotation, stratification, and thermal diffusivity is mandatory. Indeed, it constitutes one of the cornerstones of the rotational transport and mixing theory, which is implemented in stellar evolution codes to predict the rotational and chemical evolutions of stars. Aims. We investigate horizontal shear instabilities in rotating stellar radiation zones by considering the full Coriolis acceleration with both the dimensionless horizontal Coriolis component f̃ and the vertical component f. Methods. We performed a linear stability analysis using linearized equations derived from the Navier-Stokes and heat transport equations in the rotating nontraditional f-plane. We considered a horizontal shear flow with a hyperbolic tangent profile as the base flow. The linear stability was analyzed numerically in wide ranges of parameters, and we performed an asymptotic analysis for large vertical wavenumbers using the Wentzel-Kramers-Brillouin-Jeffreys (WKBJ) approximation for nondiffusive and highly-diffusive fluids. Results. As in the traditional f-plane approximation, we identify two types of instabilities: the inflectional and inertial instabilities. The inflectional instability is destabilized as f̃ increases and its maximum growth rate increases significantly, while the thermal diffusivity stabilizes the inflectional instability similarly to the traditional case. The inertial instability is also strongly affected; for instance, the inertially unstable regime is also extended in the nondiffusive limit as 0 < f < 1 + f̃ 2/N2, where N is the dimensionless Brunt-Väisälä frequency. More strikingly, in the high thermal diffusivity limit, it is always inertially unstable at any colatitude θ except at the poles (i.e., 0° < θ <  180°). We also derived the critical Reynolds numbers for the inertial instability using the asymptotic dispersion relations obtained from the WKBJ analysis. Using the asymptotic and numerical results, we propose a prescription for the effective turbulent viscosities induced by the inertial and inflectional instabilities that can be possibly used in stellar evolution models. The characteristic time of this turbulence is short enough so that it is efficient to redistribute angular momentum and to mix chemicals in stellar radiation zones.}, author = {Park, J. and Prat, V. and Mathis, S. and Bugnet, Lisa Annabelle}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, hydrodynamics / turbulence / stars, rotation / stars, evolution}, publisher = {EDP Sciences}, title = {{Horizontal shear instabilities in rotating stellar radiation zones: II. Effects of the full Coriolis acceleration}}, doi = {10.1051/0004-6361/202038654}, volume = {646}, year = {2021}, } @article{13455, abstract = {The majority of massive stars live in binary or multiple systems and will interact with a companion during their lifetimes, which helps to explain the observed diversity of core-collapse supernovae. Donor stars in binary systems can lose most of their hydrogen-rich envelopes through mass transfer. As a result, not only are the surface properties affected, but so is the core structure. However, most calculations of the core-collapse properties of massive stars rely on single-star models. We present a systematic study of the difference between the pre-supernova structures of single stars and stars of the same initial mass (11–21 M⊙) that have been stripped due to stable post-main-sequence mass transfer at solar metallicity. We present the pre-supernova core composition with novel diagrams that give an intuitive representation of the isotope distribution. As shown in previous studies, at the edge of the carbon-oxygen core, the binary-stripped star models contain an extended gradient of carbon, oxygen, and neon. This layer remains until core collapse and is more extended in mass for higher initial stellar masses. It originates from the receding of the convective helium core during core helium burning in binary-stripped stars, which does not occur in single-star models. We find that this same evolutionary phase leads to systematic differences in the final density and nuclear energy generation profiles. Binary-stripped star models have systematically higher total masses of carbon at the moment of core collapse compared to single-star models, which likely results in systematically different supernova yields. In about half of our models, the silicon-burning and oxygen-rich layers merge after core silicon burning. We discuss the implications of our findings for the “explodability”, supernova observations, and nucleosynthesis of these stars. Our models are publicly available and can be readily used as input for detailed supernova simulations.}, author = {Laplace, E. and Justham, S. and Renzo, M. and Götberg, Ylva Louise Linsdotter and Farmer, R. and Vartanyan, D. and de Mink, S. E.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{Different to the core: The pre-supernova structures of massive single and binary-stripped stars}}, doi = {10.1051/0004-6361/202140506}, volume = {656}, year = {2021}, } @article{13457, abstract = {Context. Observations of massive stars in open clusters younger than ∼8 Myr have shown that a majority of them are in binary systems, most of which will interact during their life. While these can be used as a proxy of the initial multiplicity properties, studying populations of massive stars older than ∼20 Myr allows us to probe the outcome of these interactions after a significant number of systems have experienced mass and angular momentum transfer and may even have merged. Aims. Using multi-epoch integral-field spectroscopy, we aim to investigate the multiplicity properties of the massive-star population in the dense core of the ∼40 Myr old cluster NGC 330 in the Small Magellanic Cloud in order to search for possible imprints of stellar evolution on the multiplicity properties. Methods. We obtained six epochs of VLT/MUSE observations operated in wide-field mode with the extended wavelength setup and supported by adaptive optics. We extracted spectra and measured radial velocities for stars brighter than mF814W = 19. We identified single-lined spectroscopic binaries through significant RV variability with a peak-to-peak amplitude larger than 20 km s−1. We also identified double-lined spectroscopic binaries, and quantified the observational biases for binary detection. In particular, we took into account that binary systems with similar line strengths are difficult to detect in our data set. Results. The observed spectroscopic binary fraction among stars brighter than mF814W = 19 (approximately 5.5 M⊙ on the main sequence) is fSBobs = 13.2 ± 2.0%. Considering period and mass ratio ranges from log(P) = 0.15−3.5 (about 1.4 to 3160 d), q = 0.1−1.0, and a representative set of orbital parameter distributions, we find a bias-corrected close binary fraction of fcl = 34−7+8%. This fraction seems to decline for the fainter stars, which indicates either that the close binary fraction drops in the B-type domain, or that the period distribution becomes more heavily weighted toward longer orbital periods. We further find that both fractions vary strongly in different regions of the color-magnitude diagram, which corresponds to different evolutionary stages. This probably reveals the imprint of the binary history of different groups of stars. In particular, we find that the observed spectroscopic binary fraction of Be stars (fSBobs = 2 ± 2%) is significantly lower than that of B-type stars (fSBobs = 9 ± 2%). Conclusions. We provide the first homogeneous radial velocity study of a large sample of B-type stars at a low metallicity ([Fe/H] ≲ −1.0). The overall bias-corrected close binary fraction (log(P) < 3.5 d) of the B-star population in NGC 330 is lower than the fraction reported for younger Galactic and Large Magellanic Cloud clusters in previous works. More data are needed, however, to establish whether the observed differences are caused by an age or a metallicity effect.}, author = {Bodensteiner, J. and Sana, H. and Wang, C. and Langer, N. and Mahy, L. and Banyard, G. and de Koter, A. and de Mink, S. E. and Evans, C. J. and Götberg, Ylva Louise Linsdotter and Patrick, L. R. and Schneider, F. R. N. and Tramper, F.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{The young massive SMC cluster NGC 330 seen by MUSE. II. Multiplicity properties of the massive-star population}}, doi = {10.1051/0004-6361/202140507}, volume = {652}, year = {2021}, } @article{11501, abstract = {We investigated the ultraviolet (UV) spectral properties of faint Lyman-α emitters (LAEs) in the redshift range 2.9 ≤ z ≤ 4.6, and we provide material to prepare future observations of the faint Universe. We used data from the MUSE Hubble Ultra Deep Survey to construct mean rest-frame spectra of continuum-faint (median MUV of −18 and down to MUV of −16), low stellar mass (median value of 108.4 M⊙ and down to 107 M⊙) LAEs at redshift z ≳ 3. We computed various averaged spectra of LAEs, subsampled on the basis of their observational (e.g., Lyα strength, UV magnitude and spectral slope) and physical (e.g., stellar mass and star-formation rate) properties. We searched for UV spectral features other than Lyα, such as higher ionization nebular emission lines and absorption features. We successfully observed the O III]λ1666 and [C III]λ1907+C III]λ1909 collisionally excited emission lines and the He IIλ1640 recombination feature, as well as the resonant C IVλλ1548,1551 doublet either in emission or P-Cygni. We compared the observed spectral properties of the different mean spectra and find the emission lines to vary with the observational and physical properties of the LAEs. In particular, the mean spectra of LAEs with larger Lyα equivalent widths, fainter UV magnitudes, bluer UV spectral slopes, and lower stellar masses show the strongest nebular emission. The line ratios of these lines are similar to those measured in the spectra of local metal-poor galaxies, while their equivalent widths are weaker compared to the handful of extreme values detected in individual spectra of z >  2 galaxies. This suggests that weak UV features are likely ubiquitous in high z, low-mass, and faint LAEs. We publicly released the stacked spectra, as they can serve as empirical templates for the design of future observations, such as those with the James Webb Space Telescope and the Extremely Large Telescope.}, author = {Feltre, Anna and Maseda, Michael V. and Bacon, Roland and Pradeep, Jayadev and Leclercq, Floriane and Kusakabe, Haruka and Wisotzki, Lutz and Hashimoto, Takuya and Schmidt, Kasper B. and Blaizot, Jeremy and Brinchmann, Jarle and Boogaard, Leindert and Cantalupo, Sebastiano and Carton, David and Inami, Hanae and Kollatschny, Wolfram and Marino, Raffaella A. and Matthee, Jorryt J and Nanayakkara, Themiya and Richard, Johan and Schaye, Joop and Tresse, Laurence and Urrutia, Tanya and Verhamme, Anne and Weilbacher, Peter M.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution / galaxies: high-redshift / ISM: lines and bands / ultraviolet: ISM / ultraviolet: galaxies}, publisher = {EDP Sciences}, title = {{The MUSE Hubble Ultra Deep Field Survey: XV. The mean rest-UV spectra of Lyα emitters at z > 3}}, doi = {10.1051/0004-6361/202038133}, volume = {641}, year = {2020}, } @article{11503, abstract = {Context. The Lyα emitter (LAE) fraction, XLAE, is a potentially powerful probe of the evolution of the intergalactic neutral hydrogen gas fraction. However, uncertainties in the measurement of XLAE are still under debate. Aims. Thanks to deep data obtained with the integral field spectrograph Multi Unit Spectroscopic Explorer (MUSE), we can measure the evolution of the LAE fraction homogeneously over a wide redshift range of z ≈ 3–6 for UV-faint galaxies (down to UV magnitudes of M1500 ≈ −17.75). This is a significantly fainter range than in former studies (M1500 ≤ −18.75) and it allows us to probe the bulk of the population of high-redshift star-forming galaxies. Methods. We constructed a UV-complete photometric-redshift sample following UV luminosity functions and measured the Lyα emission with MUSE using the latest (second) data release from the MUSE Hubble Ultra Deep Field Survey. Results. We derived the redshift evolution of XLAE for M1500 ∈ [ − 21.75; −17.75] for the first time with a equivalent width range EW(Lyα) ≥ 65 Å and found low values of XLAE ≲ 30% at z ≲ 6. The best-fit linear relation is XLAE = 0.07+0.06−0.03z − 0.22+0.12−0.24. For M1500 ∈ [ − 20.25; −18.75] and EW(Lyα) ≥ 25 Å, our XLAE values are consistent with those in the literature within 1σ at z ≲ 5, but our median values are systematically lower than reported values over the whole redshift range. In addition, we do not find a significant dependence of XLAE on M1500 for EW(Lyα) ≥ 50 Å at z ≈ 3–4, in contrast with previous work. The differences in XLAE mainly arise from selection biases for Lyman Break Galaxies (LBGs) in the literature: UV-faint LBGs are more easily selected if they have strong Lyα emission, hence XLAE is biased towards higher values when those samples are used. Conclusions. Our results suggest either a lower increase of XLAE towards z ≈ 6 than previously suggested, or even a turnover of XLAE at z ≈ 5.5, which may be the signature of a late or patchy reionization process. We compared our results with predictions from a cosmological galaxy evolution model. We find that a model with a bursty star formation (SF) can reproduce our observed LAE fractions much better than models where SF is a smooth function of time.}, author = {Kusakabe, Haruka and Blaizot, Jérémy and Garel, Thibault and Verhamme, Anne and Bacon, Roland and Richard, Johan and Hashimoto, Takuya and Inami, Hanae and Conseil, Simon and Guiderdoni, Bruno and Drake, Alyssa B. and Christian Herenz, Edmund and Schaye, Joop and Oesch, Pascal and Matthee, Jorryt J and Anna Marino, Raffaella and Borello Schmidt, Kasper and Pelló, Roser and Maseda, Michael and Leclercq, Floriane and Kerutt, Josephine and Mahler, Guillaume}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, dark ages / reionization / first stars / early Universe / cosmology: observations / galaxies: evolution / galaxies: high-redshift / intergalactic medium}, publisher = {EDP Sciences}, title = {{The MUSE Hubble Ultra Deep Field Survey: XIV. Evolution of the Lyα emitter fraction from z = 3 to z = 6}}, doi = {10.1051/0004-6361/201937340}, volume = {638}, year = {2020}, } @article{11504, abstract = {We present spatially resolved maps of six individually-detected Lyman α haloes (LAHs) as well as a first statistical analysis of the Lyman α (Lyα) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies (−17.5 >  MUV >  −21.5) using the Multi-Unit Spectroscopic Explorer. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lyα line profile. Using a three-dimensional two-component model for the Lyα emission, we measured the full width at half maximum (FWHM), the peak velocity shift, and the asymmetry of the Lyα line in the core and in the halo of 19 galaxies. We find that the Lyα line shape is statistically different in the halo compared to the core (in terms of width, peak wavelength, and asymmetry) for ≈40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lyα line both at the interstellar and circum-galactic scales. This trend has been predicted by radiative transfer simulations of galactic winds as a result of resonant scattering in outflows. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies (MUV <  −20) show broader, more redshifted, and less asymmetric Lyα lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lyα lines. The generally broad and red line shapes found in the halo component suggest that the Lyα haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lyα line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs.}, author = {Leclercq, Floriane and Bacon, Roland and Verhamme, Anne and Garel, Thibault and Blaizot, Jérémy and Brinchmann, Jarle and Cantalupo, Sebastiano and Claeyssens, Adélaïde and Conseil, Simon and Contini, Thierry and Hashimoto, Takuya and Herenz, Edmund Christian and Kusakabe, Haruka and Marino, Raffaella Anna and Maseda, Michael and Matthee, Jorryt J and Mitchell, Peter and Pezzulli, Gabriele and Richard, Johan and Schmidt, Kasper Borello and Wisotzki, Lutz}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics galaxies: high-redshift / galaxies: formation / galaxies: evolution / cosmology: observations}, publisher = {EDP Sciences}, title = {{The MUSE Hubble Ultra Deep field survey: XIII. Spatially resolved spectral properties of Lyman α haloes around star-forming galaxies at z > 3}}, doi = {10.1051/0004-6361/201937339}, volume = {635}, year = {2020}, } @article{13463, abstract = {Present and upcoming time-domain astronomy efforts, in part driven by gravitational-wave follow-up campaigns, will unveil a variety of rare explosive transients in the sky. Here, we focus on pulsational pair-instability evolution, which can result in signatures that are observable with electromagnetic and gravitational waves. We simulated grids of bare helium stars to characterize the resulting black hole (BH) masses together with the ejecta composition, velocity, and thermal state. We find that the stars do not react “elastically” to the thermonuclear ignition in the core: there is not a one-to-one correspondence between pair-instability driven ignition and mass ejections, which causes ambiguity as to what is an observable pulse. In agreement with previous studies, we find that for initial helium core masses of 37.5 M⊙ ≲ MHe, init ≲ 41 M⊙, corresponding to carbon-oxygen core masses 27.5 M⊙ ≲ MCO ≲ 30.1 M⊙, the explosions are not strong enough to affect the surface. With increasing initial helium core mass, they become progressively stronger causing first large radial expansion (41 M⊙ ≲ MHe, init ≲ 42 M⊙, corresponding to 30.1 M⊙ ≲ MCO ≲ 30.8 M⊙) and, finally, also mass ejection episodes (for MHe, init ≳ 42 M⊙, or MCO ≳ 30.8 M⊙). The lowest mass helium core to be fully disrupted in a pair-instability supernova is MHe, init ≃ 80 M⊙, corresponding to MCO ≃ 55 M⊙. Models with MHe, init ≳ 200 M⊙ (MCO ≳ 114 M⊙) reach the photodisintegration regime, resulting in BHs with masses of MBH ≳ 125 M⊙. Although this is currently considered unlikely, if BHs from these models form via (weak) explosions, the previously-ejected material might be hit by the blast wave and convert kinetic energy into observable electromagnetic radiation. We characterize the hydrogen-free circumstellar material from the pulsational pair-instability of helium cores by simply assuming that the ejecta maintain a constant velocity after ejection. We find that our models produce helium-rich ejecta with mass of 10−3 M⊙ ≲ MCSM ≲ 40 M⊙, the larger values corresponding to the more massive progenitor stars. These ejecta are typically launched at a few thousand km s−1 and reach distances of ∼1012 − 1015 cm before the core-collapse of the star. The delays between mass ejection events and the final collapse span a wide and mass-dependent range (from subhour to 104 years), and the shells ejected can also collide with each other, powering supernova impostor events before the final core-collapse. The range of properties we find suggests a possible connection with (some) type Ibn supernovae.}, author = {Renzo, M. and Farmer, R. and Justham, S. and Götberg, Ylva Louise Linsdotter and de Mink, S. E. and Zapartas, E. and Marchant, P. and Smith, N.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{Predictions for the hydrogen-free ejecta of pulsational pair-instability supernovae}}, doi = {10.1051/0004-6361/202037710}, volume = {640}, year = {2020}, } @article{13464, abstract = {Massive binaries that merge as compact objects are the progenitors of gravitational-wave sources. Most of these binaries experience one or more phases of mass transfer, during which one of the stars loses all or part of its outer envelope and becomes a stripped-envelope star. The evolution of the size of these stripped stars is crucial in determining whether they experience further interactions and understanding their ultimate fate. We present new calculations of stripped-envelope stars based on binary evolution models computed with MESA. We use these to investigate their radius evolution as a function of mass and metallicity. We further discuss their pre-supernova observable characteristics and potential consequences of their evolution on the properties of supernovae from stripped stars. At high metallicity, we find that practically all of the hydrogen-rich envelope is removed, which is in agreement with earlier findings. Only progenitors with initial masses below 10 M⊙ expand to large radii (up to 100 R⊙), while more massive progenitors remain compact. At low metallicity, a substantial amount of hydrogen remains and the progenitors can, in principle, expand to giant sizes (> 400 R⊙) for all masses we consider. This implies that they can fill their Roche lobe anew. We show that the prescriptions commonly used in population synthesis models underestimate the stellar radii by up to two orders of magnitude. We expect that this has consequences for the predictions for gravitational-wave sources from double neutron star mergers, particularly with regard to their metallicity dependence.}, author = {Laplace, E. and Götberg, Ylva Louise Linsdotter and de Mink, S. E. and Justham, S. and Farmer, R.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, publisher = {EDP Sciences}, title = {{The expansion of stripped-envelope stars: Consequences for supernovae and gravitational-wave progenitors}}, doi = {10.1051/0004-6361/201937300}, volume = {637}, year = {2020}, } @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}, }