@article{11512, abstract = {We study the molecular gas content of 24 star-forming galaxies at z = 3–4, with a median stellar mass of 109.1 M⊙, from the MUSE Hubble Ultra Deep Field (HUDF) Survey. Selected by their Lyα λ1216 emission and HF160W-band magnitude, the galaxies show an average $\langle {\mathrm{EW}}_{\mathrm{Ly}\alpha }^{0}\rangle \approx 20$ Å, below the typical selection threshold for Lyα emitters (${\mathrm{EW}}_{\mathrm{Ly}\alpha }^{0}\gt 25$ Å), and a rest-frame UV spectrum similar to Lyman-break galaxies. We use rest-frame optical spectroscopy from KMOS and MOSFIRE, and the UV features observed with MUSE, to determine the systemic redshifts, which are offset from Lyα by 〈Δv(Lyα)〉 = 346 km s−1, with a 100 to 600 km s−1 range. Stacking 12CO J = 4 → 3 and [C i]3P1 → 3P0 (and higher-J CO lines) from the ALMA Spectroscopic Survey of the HUDF, we determine 3σ upper limits on the line luminosities of 4.0 × 108 K km s−1pc2 and 5.6 × 108 K km s−1pc2, respectively (for a 300 km s−1 line width). Stacking the 1.2 mm and 3 mm dust-continuum flux densities, we find a 3σ upper limits of 9 μJy and 1.2 μJy, respectively. The inferred gas fractions, under the assumption of a "Galactic" CO-to-H2 conversion factor and gas-to-dust ratio, are in tension with previously determined scaling relations. This implies a substantially higher αCO ≥ 10 and δGDR ≥ 1200, consistent with the subsolar metallicity estimated for these galaxies ($12+\mathrm{log}({\rm{O}}/{\rm{H}})\approx 7.8\pm 0.2$). The low metallicity of z ≥ 3 star-forming galaxies may thus make it very challenging to unveil their cold gas through CO or dust emission, warranting further exploration of alternative tracers, such as [C ii].}, author = {Boogaard, Leindert A. and Bouwens, Rychard J. and Riechers, Dominik and van der Werf, Paul and Bacon, Roland and Matthee, Jorryt J and Stefanon, Mauro and Feltre, Anna and Maseda, Michael and Inami, Hanae and Aravena, Manuel and Brinchmann, Jarle and Carilli, Chris and Contini, Thierry and Decarli, Roberto and González-López, Jorge and Nanayakkara, Themiya and Walter, Fabian}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {1}, publisher = {IOP Publishing}, title = {{Measuring the average molecular gas content of star-forming galaxies at z = 3–4}}, doi = {10.3847/1538-4357/ac01d7}, volume = {916}, year = {2021}, } @article{11523, abstract = {We present the first results from the X-SHOOTER Lyman α survey at z = 2 (XLS-z2). XLS-z2 is a deep spectroscopic survey of 35 Lyman α emitters (LAEs) utilizing ≈90 h of exposure time with Very Large Telescope/X-SHOOTER and covers rest-frame Ly α to H α emission with R ≈ 4000. We present the sample selection, the observations, and the data reduction. Systemic redshifts are measured from rest-frame optical lines for 33/35 sources. In the stacked spectrum, our LAEs are characterized by an interstellar medium with little dust, a low metallicity, and a high ionization state. The ionizing sources are young hot stars that power strong emission lines in the optical and high-ionization lines in the ultraviolet (UV). The LAEs exhibit clumpy UV morphologies and have outflowing kinematics with blueshifted Si II absorption, a broad [O III] component, and a red-skewed Ly α line. Typically, 30 per cent of the Ly α photons escape, of which one quarter on the blue side of the systemic velocity. A fraction of Ly α photons escape directly at the systemic suggesting clear channels enabling an ≈10 per cent escape of ionizing photons, consistent with an inference based on Mg II. A combination of a low effective H I column density, a low dust content, and young starburst determines whether a star-forming galaxy is observed as an LAE. The first is possibly related to outflows and/or a fortunate viewing angle, while we find that the latter two in LAEs are typical for their stellar mass of 109 M⊙.}, author = {Matthee, Jorryt J and Sobral, David and Hayes, Matthew and Pezzulli, Gabriele and Gronke, Max and Schaerer, Daniel and Naidu, Rohan P and Röttgering, Huub and Calhau, João and Paulino-Afonso, Ana and Santos, Sérgio and Amorín, Ricardo}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: formation, galaxies: ISM, galaxies: starburst, dark ages, reionization, first stars}, number = {1}, pages = {1382--1412}, publisher = {Oxford University Press}, title = {{The X-SHOOTER Lyman α survey at z = 2 (XLS-z2) I: What makes a galaxy a Lyman α emitter?}}, doi = {10.1093/mnras/stab1304}, volume = {505}, year = {2021}, } @article{11524, abstract = {We measure the evolution of the rest-frame UV luminosity function (LF) and the stellar mass function (SMF) of Lyman-α (Ly α) emitters (LAEs) from z ∼ 2 to z ∼ 6 by exploring ∼4000 LAEs from the SC4K sample. We find a correlation between Ly α luminosity (LLy α) and rest-frame UV (MUV), with best fit MUV=−1.6+0.2−0.3log10(LLyα/ergs−1)+47+12−11 and a shallower relation between LLy α and stellar mass (M⋆), with best fit log10(M⋆/M⊙)=0.9+0.1−0.1log10(LLyα/ergs−1)−28+4.0−3.8⁠. An increasing LLy α cut predominantly lowers the number density of faint MUV and low M⋆ LAEs. We estimate a proxy for the full UV LFs and SMFs of LAEs with simple assumptions of the faint end slope. For the UV LF, we find a brightening of the characteristic UV luminosity (M∗UV⁠) with increasing redshift and a decrease of the characteristic number density (Φ*). For the SMF, we measure a characteristic stellar mass (⁠M∗⋆/M⊙⁠) increase with increasing redshift, and a Φ* decline. However, if we apply a uniform luminosity cut of log10(LLyα/ergs−1)≥43.0⁠, we find much milder to no evolution in the UV and SMF of LAEs. The UV luminosity density (ρUV) of the full sample of LAEs shows moderate evolution and the stellar mass density (ρM) decreases, with both being always lower than the total ρUV and ρM of more typical galaxies but slowly approaching them with increasing redshift. Overall, our results indicate that both ρUV and ρM of LAEs slowly approach the measurements of continuum-selected galaxies at z > 6, which suggests a key role of LAEs in the epoch of reionization.}, author = {Santos, S and Sobral, D and Butterworth, J and Paulino-Afonso, A and Ribeiro, B and da Cunha, E and Calhau, J and Khostovan, A A and Matthee, Jorryt J and Arrabal Haro, P}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, galaxies: evolution, galaxies: high-redshift, galaxies: luminosity function, mass function}, number = {1}, pages = {1117--1134}, publisher = {Oxford University Press}, title = {{The evolution of the UV luminosity and stellar mass functions of Lyman-α emitters from z ∼ 2 to z ∼ 6}}, doi = {10.1093/mnras/stab1218}, volume = {505}, year = {2021}, } @article{11525, abstract = {The intensity of the Cosmic UV background (UVB), coming from all sources of ionizing photons such as star-forming galaxies and quasars, determines the thermal evolution and ionization state of the intergalactic medium (IGM) and is, therefore, a critical ingredient for models of cosmic structure formation. Most of the previous estimates are based on the comparison between observed and simulated Lyman-α forest. We present the results of an independent method to constrain the product of the UVB photoionization rate and the covering fraction of Lyman limit systems (LLSs) by searching for the fluorescent Lyman-α emission produced by self-shielded clouds. Because the expected surface brightness is well below current sensitivity limits for direct imaging, we developed a new method based on 3D stacking of the IGM around Lyman-α emitting galaxies (LAEs) between 2.9 < z < 6.6 using deep MUSE observations. Combining our results with covering fractions of LLSs obtained from mock cubes extracted from the EAGLE simulation, we obtain new and independent constraints on the UVB at z > 3 that are consistent with previous measurements, with a preference for relatively low UVB intensities at z = 3, and which suggest a non-monotonic decrease of ΓH I with increasing redshift between 3 < z < 5. This could suggest a possible tension between some UVB models and current observations which however require deeper and wider observations in Lyman-α emission and absorption to be confirmed. Assuming instead a value of UVB from current models, our results constrain the covering fraction of LLSs at 3 < z < 4.5 to be less than 25 per cent within 150 kpc from LAEs.}, author = {Gallego, Sofia G and Cantalupo, Sebastiano and Sarpas, Saeed and Duboeuf, Bastien and Lilly, Simon and Pezzulli, Gabriele and Marino, Raffaella Anna and Matthee, Jorryt J and Wisotzki, Lutz and Schaye, Joop and Richard, Johan and Kusakabe, Haruka and Mauerhofer, Valentin}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {1}, pages = {16--32}, publisher = {Oxford University Press}, title = {{Constraining the cosmic UV background at z > 3 with MUSE Lyman-α emission observations}}, doi = {10.1093/mnras/stab796}, volume = {504}, year = {2021}, } @article{11526, abstract = {We present the results from a MUSE survey of twelve z ≃ 3.15 quasars, which were selected to be much fainter (20 < iSDSS < 23) than in previous studies of giant Ly α nebulae around the brightest quasars (16.6 < iAB < 18.7). We detect H I Ly α nebulae around 100 per cent of our target quasars, with emission extending to scales of at least 60 physical kpc, and up to 190 pkpc. We explore correlations between properties of the nebulae and their host quasars, with the goal of connecting variations in the properties of the illuminating QSO to the response in nebular emission. We show that the surface brightness profiles of the nebulae are similar to those of nebulae around bright quasars, but with a lower normalization. Our targeted quasars are on average 3.7 mag (≃30 times) fainter in UV continuum than our bright reference sample, and yet the nebulae around them are only 4.3 times fainter in mean Ly α surface brightness, measured between 20 and 50 pkpc. We find significant correlations between the surface brightness of the nebula and the luminosity of the quasar in both UV continuum and Ly α. The latter can be interpreted as evidence for a substantial contribution from unresolved inner parts of the nebulae to the narrow components seen in the Ly α lines of some of our faint quasars, possibly from the inner circumgalactic medium or from the host galaxy’s interstellar medium.}, author = {Mackenzie, Ruari and Pezzulli, Gabriele and Cantalupo, Sebastiano and Marino, Raffaella A and Lilly, Simon and Muzahid, Sowgat and Matthee, Jorryt J and Schaye, Joop and Wisotzki, Lutz}, issn = {1365-2966}, journal = {Monthly Notices of the Royal Astronomical Society}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, techniques: imaging spectroscopy, intergalactic medium, quasars: emission lines, quasars: general}, number = {1}, pages = {494--509}, publisher = {Oxford University Press}, title = {{Revealing the impact of quasar luminosity on giant Lyα nebulae}}, doi = {10.1093/mnras/staa3277}, volume = {502}, year = {2021}, } @article{11604, abstract = {The NASA Transiting Exoplanet Survey Satellite (TESS) is observing tens of millions of stars with time spans ranging from ∼27 days to about 1 yr of continuous observations. This vast amount of data contains a wealth of information for variability, exoplanet, and stellar astrophysics studies but requires a number of processing steps before it can be fully utilized. In order to efficiently process all the TESS data and make it available to the wider scientific community, the TESS Data for Asteroseismology working group, as part of the TESS Asteroseismic Science Consortium, has created an automated open-source processing pipeline to produce light curves corrected for systematics from the short- and long-cadence raw photometry data and to classify these according to stellar variability type. We will process all stars down to a TESS magnitude of 15. This paper is the next in a series detailing how the pipeline works. Here, we present our methodology for the automatic variability classification of TESS photometry using an ensemble of supervised learners that are combined into a metaclassifier. We successfully validate our method using a carefully constructed labeled sample of Kepler Q9 light curves with a 27.4 days time span mimicking single-sector TESS observations, on which we obtain an overall accuracy of 94.9%. We demonstrate that our methodology can successfully classify stars outside of our labeled sample by applying it to all ∼167,000 stars observed in Q9 of the Kepler space mission.}, author = {Audenaert, J. and Kuszlewicz, J. S. and Handberg, R. and Tkachenko, A. and Armstrong, D. J. and Hon, M. and Kgoadi, R. and Lund, M. N. and Bell, K. J. and Bugnet, Lisa Annabelle and Bowman, D. M. and Johnston, C. and García, R. A. and Stello, D. and Molnár, L. and Plachy, E. and Buzasi, D. and Aerts, C.}, issn = {1538-3881}, journal = {The Astronomical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {5}, publisher = {IOP Publishing}, title = {{TESS Data for Asteroseismology (T’DA) stellar variability classification pipeline: Setup and application to the Kepler Q9 data}}, doi = {10.3847/1538-3881/ac166a}, volume = {162}, 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{13453, abstract = {Most massive stars are born in binaries close enough for mass transfer episodes. These modify the appearance, structure, and future evolution of both stars. We compute the evolution of a 100-day-period binary, consisting initially of a 25 M⊙ star and a 17 M⊙ star, which experiences stable mass transfer. We focus on the impact of mass accretion on the surface composition, internal rotation, and structure of the accretor. To anchor our models, we show that our accretor broadly reproduces the properties of ζ Ophiuchi, which has long been proposed to have accreted mass before being ejected as a runaway star when the companion exploded. We compare our accretor to models of single rotating stars and find that the later and stronger spin-up provided by mass accretion produces significant differences. Specifically, the core of the accretor retains higher spin at the end of the main sequence, and a convective layer develops that changes its density profile. Moreover, the surface of the accretor star is polluted by CNO-processed material donated by the companion. Our models show effects of mass accretion in binaries that are not captured in single rotating stellar models. This possibly impacts the further evolution (either in a binary or as single stars), the final collapse, and the resulting spin of the compact object.}, author = {Renzo, M. and Götberg, Ylva Louise Linsdotter}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {2}, publisher = {American Astronomical Society}, title = {{Evolution of accretor stars in massive binaries: Broader implications from modeling ζ Ophiuchi}}, doi = {10.3847/1538-4357/ac29c5}, volume = {923}, year = {2021}, } @article{13454, abstract = {Helium star–carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code MESA, we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the time of explosion. We convert these model properties to near-UV/optical magnitudes assuming a blackbody spectrum and support this approach using a matched stellar atmosphere model. These models will be valuable to compare with pre-explosion imaging for future supernovae, though we emphasize the observational difficulty of detecting extremely blue companions. The pre-explosion source detected in association with SN 2012Z has been interpreted as a helium star binary containing an initially ultra-massive WD in a multiday orbit. However, extending our binary models to initial CO WD masses of up to 1.2 M⊙, we find that these systems undergo off-center carbon ignitions and thus are not expected to produce thermonuclear supernovae. This tension suggests that, if SN 2012Z is associated with a helium star–WD binary, then the pre-explosion optical light from the system must be significantly modified by the binary environment and/or the WD does not have a carbon-rich interior composition.}, author = {Wong, Tin Long Sunny and Schwab, Josiah and Götberg, Ylva Louise Linsdotter}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {2}, publisher = {American Astronomical Society}, title = {{Pre-explosion properties of Helium star donors to thermonuclear supernovae}}, doi = {10.3847/1538-4357/ac27ae}, volume = {922}, 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{13456, abstract = {While most simulations of the epoch of reionization have focused on single-stellar populations in star-forming dwarf galaxies, products of binary evolution are expected to significantly contribute to emissions of hydrogen-ionizing photons. Among these products are stripped stars (or helium stars), which have their envelopes stripped from interactions with binary companions, leaving an exposed helium core. Previous work has suggested these stripped stars can dominate the Lyman Continuum (LyC) photon output of high-redshift, low-luminosity galaxies post-starburst. Other sources of hard radiation in the early universe include zero-metallicity Population iii stars, which may have similar spectral energy distribution (SED) properties to galaxies with radiation dominated by stripped-star emissions. Here, we use four metrics (the power-law exponent over wavelength intervals 240–500 Å, 600–900 Å, and 1200–2000 Å, and the ratio of total luminosity in FUV wavelengths to LyC wavelengths) to compare the SEDs of simulated galaxies with only single-stellar evolution, galaxies containing stripped stars, and galaxies containing Population iii stars, with four different initial mass functions (IMFs). We find that stripped stars significantly alter SEDs in the LyC range of galaxies at the epoch of reionization. SEDs in galaxies with stripped stars have lower power-law indices in the LyC range and lower FUV to LyC luminosity ratios. These differences in SEDs are present at all considered luminosities (${M}_{\mathrm{UV}}\gt -15$, AB system), and are most pronounced for lower-luminosity galaxies. Intrinsic SEDs as well as those with interstellar medium absorption of galaxies with stripped stars and Population iii stars are found to be distinct for all tested Population iii IMFs.}, author = {Berzin, Elizabeth and Secunda, Amy and Cen, Renyue and Menegas, Alexander and Götberg, Ylva Louise Linsdotter}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {1}, publisher = {American Astronomical Society}, title = {{Spectral signatures of population III and envelope-stripped stars in galaxies at the epoch of reionization}}, doi = {10.3847/1538-4357/ac0af6}, volume = {918}, 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{13458, abstract = {Most massive stars experience binary interactions in their lifetimes that can alter both the surface and core structure of the stripped star with significant effects on their ultimate fate as core-collapse supernovae. However, core-collapse supernovae simulations to date have focused almost exclusively on the evolution of single stars. We present a systematic simulation study of single and binary-stripped stars with the same initial mass as candidates for core-collapse supernovae (11–21 M⊙). Generally, we find that binary-stripped stars core tend to have a smaller compactness parameter, with a more prominent, deeper silicon/oxygen interface, and explode preferentially to the corresponding single stars of the same initial mass. Such a dichotomy of behavior between these two modes of evolution would have important implications for supernovae statistics, including the final neutron star masses, explosion energies, and nucleosynthetic yields. Binary-stripped remnants are also well poised to populate the possible mass gap between the heaviest neutron stars and the lightest black holes. Our work presents an improvement along two fronts, as we self-consistently account for the pre-collapse stellar evolution and the subsequent explosion outcome. Even so, our results emphasize the need for more detailed stellar evolutionary models to capture the sensitive nature of explosion outcome.}, author = {Vartanyan, David and Laplace, Eva and Renzo, Mathieu and Götberg, Ylva Louise Linsdotter and Burrows, Adam and de Mink, Selma E.}, issn = {2041-8213}, journal = {The Astrophysical Journal Letters}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {1}, publisher = {American Astronomical Society}, title = {{Binary-stripped stars as core-collapse supernovae progenitors}}, doi = {10.3847/2041-8213/ac0b42}, volume = {916}, year = {2021}, } @article{13459, abstract = {The B emission-line stars are rapid rotators that were probably spun up by mass and angular momentum accretion through mass transfer in an interacting binary. Mass transfer will strip the donor star of its envelope to create a small and hot subdwarf remnant. Here we report on Hubble Space Telescope/STIS far-ultraviolet spectroscopy of a sample of Be stars that reveals the presence of the hot sdO companion through the calculation of cross-correlation functions of the observed and model spectra. We clearly detect the spectral signature of the sdO star in 10 of the 13 stars in the sample, and the spectral signals indicate that the sdO stars are hot, relatively faint, and slowly rotating as predicted by models. A comparison of their temperatures and radii with evolutionary tracks indicates that the sdO stars occupy the relatively long-lived, He-core burning stage. Only 1 of the 10 detections was a known binary prior to this investigation, which emphasizes the difficulty of finding such Be+sdO binaries through optical spectroscopy. However, these results and others indicate that many Be stars probably host hot subdwarf companions.}, author = {Wang, Luqian and Gies, Douglas R. and Peters, Geraldine J. and Götberg, Ylva Louise Linsdotter and Chojnowski, S. Drew and Lester, Kathryn V. and Howell, Steve B.}, issn = {1538-3881}, journal = {The Astronomical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {5}, publisher = {American Astronomical Society}, title = {{The detection and characterization of Be+sdO binaries from HST/STIS FUV spectroscopy}}, doi = {10.3847/1538-3881/abf144}, volume = {161}, year = {2021}, } @article{12583, abstract = {Peruvian glaciers are important contributors to dry season runoff for agriculture and hydropower, but they are at risk of disappearing due to climate change. We applied a physically based, energy balance melt model at five on-glacier sites within the Peruvian Cordilleras Blanca and Vilcanota. Net shortwave radiation dominates the energy balance, and despite this flux being higher in the dry season, melt rates are lower due to losses from net longwave radiation and the latent heat flux. The sensible heat flux is a relatively small contributor to melt energy. At three of the sites the wet season snowpack was discontinuous, forming and melting within a daily to weekly timescale, and resulting in highly variable melt rates closely related to precipitation dynamics. Cold air temperatures due to a strong La Niña year at Shallap Glacier (Cordillera Blanca) resulted in a continuous wet season snowpack, significantly reducing wet season ablation. Sublimation was most important at the highest site in the accumulation zone of the Quelccaya Ice Cap (Cordillera Vilcanota), accounting for 81% of ablation, compared to 2%–4% for the other sites. Air temperature and precipitation inputs were perturbed to investigate the climate sensitivity of the five glaciers. At the lower sites warmer air temperatures resulted in a switch from snowfall to rain, so that ablation was increased via the decrease in albedo and increase in net shortwave radiation. At the top of Quelccaya Ice Cap warming caused melting to replace sublimation so that ablation increased nonlinearly with air temperature.}, author = {Fyffe, Catriona L. and Potter, Emily and Fugger, Stefan and Orr, Andrew and Fatichi, Simone and Loarte, Edwin and Medina, Katy and Hellström, Robert Å. and Bernat, Maud and Aubry‐Wake, Caroline and Gurgiser, Wolfgang and Perry, L. Baker and Suarez, Wilson and Quincey, Duncan J. and Pellicciotti, Francesca}, issn = {2169-8996}, journal = {Journal of Geophysical Research: Atmospheres}, keywords = {Space and Planetary Science, Earth and Planetary Sciences (miscellaneous), Atmospheric Science, Geophysics}, number = {23}, publisher = {American Geophysical Union}, title = {{The energy and mass balance of Peruvian Glaciers}}, doi = {10.1029/2021jd034911}, volume = {126}, 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}, }