@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{11498, abstract = {Rest-frame ultraviolet (UV) emission lines probe electron densities, gas-phase abundances, metallicities, and ionization parameters of the emitting star-forming galaxies and their environments. The strongest main UV emission line, Lyα, has been instrumental in advancing the general knowledge of galaxy formation in the early universe. However, observing Lyα emission becomes increasingly challenging at z ≳ 6 when the neutral hydrogen fraction of the circumgalactic and intergalactic media increases. Secondary weaker UV emission lines provide important alternative methods for studying galaxy properties at high redshift. We present a large sample of rest-frame UV emission line sources at intermediate redshift for calibrating and exploring the connection between secondary UV lines and the emitting galaxies’ physical properties and their Lyα emission. The sample of 2052 emission line sources with 1.5 < z < 6.4 was collected from integral field data from the MUSE-Wide and MUSE-Deep surveys taken as part of Guaranteed Time Observations. The objects were selected through untargeted source detection (i.e., no preselection of sources as in dedicated spectroscopic campaigns) in the three-dimensional MUSE data cubes. We searched optimally extracted one-dimensional spectra of the full sample for UV emission features via emission line template matching, resulting in a sample of more than 100 rest-frame UV emission line detections. We show that the detection efficiency of (non-Lyα) UV emission lines increases with survey depth, and that the emission line strength of He IIλ1640 Å, [O III] λ1661 + O III] λ1666, and [Si III] λ1883 + Si III] λ1892 correlate with the strength of [C III] λ1907 + C III] λ1909. The rest-frame equivalent width (EW0) of [C III] λ1907 + C III] λ1909 is found to be roughly 0.22 ± 0.18 of EW0(Lyα). We measured the velocity offsets of resonant emission lines with respect to systemic tracers. For C IVλ1548 + C IVλ1551 we find that ΔvC IV ≲ 250 km s−1, whereas ΔvLyα falls in the range of 250−500 km s−1 which is in agreement with previous results from the literature. The electron density ne measured from [Si III] λ1883 + Si III] λ1892 and [C III] λ1907 + C III] λ1909 line flux ratios is generally < 105 cm−3 and the gas-phase abundance is below solar at 12 + log10(O/H)≈8. Lastly, we used “PhotoIonization Model Probability Density Functions” to infer physical parameters of the full sample and individual systems based on photoionization model parameter grids and observational constraints from our UV emission line searches. This reveals that the UV line emitters generally have ionization parameter log10(U) ≈ −2.5 and metal mass fractions that scatter around Z ≈ 10−2, that is Z ≈ 0.66 Z⊙. Value-added catalogs of the full sample of MUSE objects studied in this work and a collection of UV line emitters from the literature are provided with this paper.}, author = {Schmidt, K. B. and Kerutt, J. and Wisotzki, L. and Urrutia, T. and Feltre, A. and Maseda, M. V. and Nanayakkara, T. and Bacon, R. and Boogaard, L. A. and Conseil, S. and Contini, T. and Herenz, E. C. and Kollatschny, W. and Krumpe, M. and Leclercq, F. and Mahler, G. and Matthee, Jorryt J and Mauerhofer, V. and Richard, J. and Schaye, J.}, issn = {1432-0746}, journal = {Astronomy & Astrophysics}, keywords = {Space and Planetary Science, Astronomy and Astrophysics, ultraviolet: galaxies / galaxies: high-redshift / galaxies: ISM / ISM: lines and bands / methods: observational / techniques: imaging spectroscopy}, publisher = {EDP Sciences}, title = {{Recovery and analysis of rest-frame UV emission lines in 2052 galaxies observed with MUSE at 1.5 < z < 6.4}}, doi = {10.1051/0004-6361/202140876}, volume = {654}, year = {2021}, } @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{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{13460, abstract = {Binary interaction can cause stellar envelopes to be stripped, which significantly reduces the radius of the star. The orbit of a binary composed of a stripped star and a compact object can therefore be so tight that the gravitational radiation the system produces reaches frequencies accessible to the Laser Interferometer Space Antenna (LISA). Two such stripped stars in tight orbits with white dwarfs are known so far (ZTF J2130+4420 and CD−30°11223), but many more are expected to exist. These binaries provide important constraints for binary evolution models and may be used as LISA verification sources. We develop a Monte Carlo code that uses detailed evolutionary models to simulate the Galactic population of stripped stars in tight orbits with either neutron star or white dwarf companions. We predict 0–100 stripped star + white dwarf binaries and 0–4 stripped star + neutron star binaries with a signal-to-noise ratio >5 after 10 yr of observations with LISA. More than 90% of these binaries are expected to show large radial velocity shifts of ≳200 $\,\mathrm{km}\,{{\rm{s}}}^{-1}$, which are spectroscopically detectable. Photometric variability due to tidal deformation of the stripped star is also expected and has been observed in ZTF J2130+4420 and CD−30°11223. In addition, the stripped star + neutron star binaries are expected to be X-ray bright with LX ≳ 1033–1036 $\,\mathrm{erg}\,{{\rm{s}}}^{-1}$. Our results show that stripped star binaries are promising multimessenger sources for the upcoming electromagnetic and gravitational wave facilities.}, author = {Götberg, Ylva Louise Linsdotter and Korol, V. and Lamberts, A. and Kupfer, T. and Breivik, K. and Ludwig, B. and Drout, M. R.}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {1}, publisher = {American Astronomical Society}, title = {{Stars stripped in binaries: The living gravitational-wave sources}}, doi = {10.3847/1538-4357/abbda5}, volume = {904}, year = {2020}, } @article{13461, abstract = {High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (fesc) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is ≳70%. Moreover, simulations so far have yielded values of fesc falling only on the lower end of the ∼10%–20% range, the amount presumed necessary to reionize the universe. Analyzing a high-resolution (4 pc) cosmological radiation-hydrodynamic simulation, we study how fesc changes when we include two different products of binary stellar evolution—stars stripped of their hydrogen envelopes and massive blue stragglers. Both produce significant amounts of ionizing photons 10–200 Myr after each starburst. We find the relative importance of these photons to be amplified with respect to escaped ionizing photons, because peaks in star formation rates (SFRs) and fesc are often out of phase by this 10–200 Myr. Additionally, low-mass, bursty galaxies emit Lyman continuum radiation primarily from binary products when SFRs are low. Observations of these galaxies by the James Webb Space Telescope could provide crucial information on the evolution of binary stars as a function of redshift. Overall, including stripped stars and massive blue stragglers increases our photon-weighted mean escape fraction ($\langle {f}_{\mathrm{esc}}\rangle $) by ∼13% and ∼10%, respectively, resulting in $\langle {f}_{\mathrm{esc}}\rangle =17 \% $. Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization.}, author = {Secunda, Amy and Cen, Renyue and Kimm, Taysun and Götberg, Ylva Louise Linsdotter and de Mink, Selma E.}, issn = {1538-4357}, journal = {The Astrophysical Journal}, keywords = {Space and Planetary Science, Astronomy and Astrophysics}, number = {1}, publisher = {American Astronomical Society}, title = {{Delayed photons from binary evolution help reionize the universe}}, doi = {10.3847/1538-4357/abaefa}, volume = {901}, year = {2020}, }