Impact of the primordial stellar initial mass function on the 21-cm signal

Properties of the first generation of stars [referred to as the Population III (Pop III) stars], such as their initial mass function (IMF), are poorly constrained by observations and have yet to converge between simulations. The cosmological 21-cm signal of neutral hydrogen is predicted to be sensitive to Lyman-band photons produced by these stars, thus providing a unique way to probe the first stellar population. In this paper, we investigate the impacts of the Pop III IMF on the cosmic-dawn 21-cm signal via the Wouthuysen-Field effect, Lyman-Werner feedback, Ly alpha heating, and cosmic microwave background heating. We calculate the emission spectra of star-forming haloes for different IMFs by integrating over individual metal-free stellar spectra, computed from a set of stellar evolution histories and stellar atmospheres, and taking into account variability of the spectra with stellar age. Through this study, we therefore relax two common assumptions: that the zero-age main-sequence emission rate of a Pop III star is representative of its lifetime mean emission rate, and that Pop III emission can be treated as instantaneous. Exploring bottom-heavy, top-heavy, and intermediate IMFs, we show that variations in the 21-cm signal are driven by stars lighter than 20 M-circle dot. For the explored models, we find maximum relative differences of 59 per cent in the cosmic-dawn global 21-cm signal, and 131 per cent between power spectra. Although this impact is modest, precise modelling of the first stars and their evolution is necessary for accurate prediction and interpretation of the 21-cm signal.UK Research & Innovation (UKRI)Science & Technology Facilities Council (STFC)Science and Technology Development Fund (STDF) ST/V506606/1AF's Royal Society University Research Fellowship 181073 180523UK Research & Innovation (UKRI)Science & Technology Facilities Council (STFC) ST/R000603/1Max-Planck-Gesellschaft via the fellowship of the International Max Planck Research School for Astronomy and Cosmic Physics at the University of Heidelberg (IMPRS-HD)UK Research & Innovation (UKRI)Science & Technology Facilities Council (STFC) ST/R000603/1 ST/L003910/1Royal Society of LondonIsrael Science Foundation 2359/20Ambrose Monell FoundationInstitute for Advanced StudyVera Rubin Presidential Chair in AstronomyThe David & Lucile Packard Foundatio

Discovery of starspots on Vega - First spectroscopic detection of surface structures on a normal A-type star

The theoretically studied impact of rapid rotation on stellar evolution needs to be confronted with the results of high resolution spectroscopy-velocimetry observations. A weak surface magnetic field had recently been detected in the A0 prototype star Vega, potentially leading to a (yet undetected) structured surface. The goal of this article is to present a thorough analysis of the line profile variations and associated estimators in the early-type standard star Vega (A0) in order reveal potential activity tracers, exoplanet companions and stellar oscillations. Vega was monitored in high-resolution spectroscopy with the velocimeter Sophie/OHP. A total of 2588 high S/N spectra was obtained during 5 nights (August 2012) at R = 75000 and covering the visible domain. For each reduced spectrum, Least Square Deconvolved (LSD) equivalent photospheric profiles were calculated with a Teff = 9500 and logg = 4.0 spectral line mask. Several methods were applied to study the dynamic behavior of the profile variations (evolution of radial velocity, bisectors, vspan, 2D profiles, amongst others). We present the discovery of a starspotted stellar surface in an A-type standard star with faint spot amplitudes Delta F/Fc ~5 10^{-4}. A rotational modulation of spectral lines with a period of rotation P = 0.68 d has clearly been exhibited, confirming the results of previous spectropolarimetric studies. Either a very thin convective layer can be responsible for magnetic field generation at small amplitudes, or a new mechanism has to be invoked in order to explain the existence of activity tracing starspots. This first strong evidence that standard A-type stars can show surface structures opens a new field of research and asks the question about a potential link with the recently discovered weak magnetic field discoveries in this category of stars.Comment: accepted for publication by Astronomy & Astrophysics (23rd of March 2015

A new model for mixing by double-diffusive convection (semi-convection): I. The conditions for layer formation

The process referred to as "semi-convection" in astrophysics and "double-diffusive convection in the diffusive regime" in Earth and planetary sciences, occurs in stellar and planetary interiors in regions which are stable according to the Ledoux criterion but unstable according to the Schwarzschild criterion. In this series of papers, we analyze the results of an extensive suite of 3D numerical simulations of the process, and ultimately propose a new 1D prescription for heat and compositional transport in this regime which can be used in stellar or planetary structure and evolution models. In a preliminary study of the phenomenon, Rosenblum et al. (2011) showed that, after saturation of the primary instability, a system can evolve in one of two possible ways: the induced turbulence either remains homogeneous, with very weak transport properties, or transitions into a thermo-compositional staircase where the transport rate is much larger (albeit still smaller than in standard convection). In this paper, we show that this dichotomous behavior is a robust property of semi-convection across a wide region of parameter space. We propose a simple semi-analytical criterion to determine whether layer formation is expected or not, and at what rate it proceeds, as a function of the background stratification and of the diffusion parameters (viscosity, thermal diffusivity and compositional diffusivity) only. The theoretical criterion matches the outcome of our numerical simulations very adequately in the numerically accessible "planetary" parameter regime, and can easily be extrapolated to the stellar parameter regime. Subsequent papers will address more specifically the question of quantifying transport in the layered case and in the non-layered case.Comment: Submitted to Ap

Detailed equilibrium and dynamical tides: impact on circularization and synchronization in open clusters

Binary stars evolve into chemically-peculiar objects and are a major driver of the Galactic enrichment of heavy elements. During their evolution they undergo interactions, including tides, that circularize orbits and synchronize stellar spins, impacting both individual systems and stellar populations. Using Zahn's tidal theory and MESA main-sequence model grids, we derive the governing parameters $\lambda_{lm}$ and $E_2$, and implement them in the new MINT library of the stellar population code BINARY_C. Our MINT equilibrium tides are 2 to 5 times more efficient than the ubiquitous BSE prescriptions while the radiative-tide efficiency drops sharply with increasing age. We also implement precise initial distributions based on bias-corrected observations. We assess the impact of tides and initial orbital-parameter distributions on circularization and synchronization in eight open clusters, comparing synthetic populations and observations through a bootstrapping method. We find that changing the tidal prescription yields no statistically-significant improvement as both calculations typically lie within 0.5$\sigma$. The initial distribution, especially the primordial concentration of systems at $\log_{10}(P/{\rm d}) \approx 0.8, e\approx 0.05$ dominates the statistics even when artificially increasing tidal strength. This confirms the inefficiency of tides on the main sequence and shows that constraining tidal-efficiency parameters using the $e-\log_{10}(P/{\rm d})$ distribution alone is difficult or impossible. Orbital synchronization carries a more striking age-dependent signature of tidal interactions. In M35 we find twice as many synchronized rotators in our MINT calculation as with BSE. This measure of tidal efficiency is verifiable with combined measurements of orbital parameters and stellar spins.Comment: 24 pages, 29 figures includings appendices. Accepted for publication in MNRA

Spectroscopic detection of Altair's non-radial pulsations

Context: Rapid rotation is a common feature of early-type stars but which remains a challenge for the models. The understanding of its effect on stellar evolution is however imperative to interpret the observed properties of numerous stars. Aims: We wish to bring more observational constraints on the properties of fast rotating stars, especially on their oscillation modes. Methods: We focus on the nearby star Altair which is known as a very rapidly rotating star with an equatorial velocity estimated recently at 313 km/s. We observed this star with the high-resolution spectropolarimeter Neo-Narval during six nights, with one night of interruption, in September 2020. Results: We detect significant line profile variations on the mean line profile of the spectra. Their time-frequency analysis shows that these variations are induced by gravito-inertial waves propagating at Altair's surface with azimuthal wavenumbers of order $m=10-15$. With a preliminary computation of the eigenspectrum using the most recent concordance model of Altair we can give a first modelling of the observed waves. Conclusions: Altair was known as the brightest $\delta$ Scuti star. We now see that it is the brightest hybrid oscillating star with excited gravito-inertial waves and acoustic waves. Clearly, more observations and more advanced models are needed to explain the observations in greater detailsComment: 8 pages, 8 figures, submitted to Astronomy and Astrophysic

Dating young open clusters using delta Scuti stars. Results for Trumpler 10 and Praesepe

Aims. The main goal of this work is to date young open clusters using $\delta$ Sct stars. Seismic indices such as the large separation and the frequency at maximum power can help to constrain the models to better characterise the stars. We propose a reliable method to identify some radial modes, which gives us greater confidence in the constrained models. Methods. We extract the frequency content of a sample of $\delta$ Sct stars belonging to the same open cluster. We estimate the low-order large separation by means of different techniques and the frequency at maximum power for each member of the sample. We use a grid of models built with the typical parameters of $\delta$ Sct stars, including mass, metallicity and rotation as independent variables, and determine the oscillation modes. We select the observed frequencies whose ratios match those of the models. Once we find a range of radial modes matching the observed frequencies, mainly the fundamental mode, we add it to the other seismic parameters to derive the stellar age. Assuming star groups have similar chemistry and age, we estimate their mean age by computing a weighted probability density function fit to the age distribution of the seismically constrained models. Results. We estimate the age of Trumpler 10 to be $30_{-20}{+30}$ Myr, and that of Praesepe to be $580 \pm 230$ Myr. In this latter case, we find two apparent populations of $\delta$ Sct stars in the same cluster, one at $510 \pm 140$ Myr and another at $890 \pm 140$ Myr. This may be due to two different formation events, different rotational velocities of the members in our sample of stars (as rapid rotation may modify the observed large separation), or to membership of unresolved binary systems.Comment: 20 pages, 25 figure, uses the open-source code MultiModes (see https://github.com/davidpamos/MultiModes). It will be published in A&

The PL diagram for δ\delta Sct stars: back in business as distance estimators

In this work, we focus on the period-luminosity relation (PLR) of $\delta$ Sct stars, in which mode excitation and selection mechanisms are still poorly constrained, and whose structure and oscillations are affected by rotation. We review the PLRs in the recent literature, and add a new inference from a large sample of $\delta$ Sct. We highlight the difficulty in identifying the fundamental mode and show that rotation-induced surface effects can impact the measured luminosities, explaining the broadening of the PLR. We derive a tight relation between the low-order large separation and the fundamental radial mode frequency (F0) that holds for rotating stars, thus paving the way towards mode identification. We show that the PLRs we obtain for different samples are compatible with each other and with the recent literature, and with most observed $\delta$ Sct stars when taking rotation effects into account. We also find that the highest-amplitude peak in the frequency spectrum corresponds to the fundamental mode in most $\delta$ Sct, thus shedding some light on their elusive mode selection mechanism.Comment: 10 pages, 5 figures, 1 table, IAU conference proceedin

Population III X-ray Binaries and their Impact on the Early Universe

The first population of X-ray binaries (XRBs) is expected to affect the thermal and ionization states of the gas in the early Universe. Although these X-ray sources are predicted to have important implications for high-redshift observable signals, such as the hydrogen 21-cm signal from cosmic dawn and the cosmic X-ray background, their properties are poorly explored, leaving theoretical models largely uninformed. In this paper we model a population of X-ray binaries arising from zero metallicity stars. We explore how their properties depend on the adopted initial mass function (IMF) of primordial stars, finding a strong effect on their number and X-ray production efficiency. We also present scaling relations between XRBs and their X-ray emission with the local star formation rate, which can be used in sub-grid models in numerical simulations to improve the X-ray feedback prescriptions. Specifically, we find that the uniformity and strength of the X-ray feedback in the intergalactic medium is strongly dependant on the IMF. Bottom-heavy IMFs result in a smoother distribution of XRBs, but have a luminosity orders of magnitude lower than more top-heavy IMFs. Top-heavy IMFs lead to more spatially uneven, albeit strong, X-ray emission. An intermediate IMF has a strong X-ray feedback while sustaining an even emission across the intergalactic medium. These differences in X-ray feedback could be probed in the future with measurements of the cosmic dawn 21-cm line of neutral hydrogen, which offers us a new way of constraining population III IMF.Comment: Accepted for publication in MNRAS, 17 pages, 9 figure

Asteroseismology of massive stars with the TESS mission: the runaway Beta Cep pulsator PHL 346 = HN Aqr

We report an analysis of the first known Beta Cep pulsator observed by the TESS mission, the runaway star PHL 346 = HN Aqr. The star, previously known as a singly-periodic pulsator, has at least 34 oscillation modes excited, 12 of those in the g-mode domain and 22 p modes. Analysis of archival data implies that the amplitude and frequency of the dominant mode and the stellar radial velocity were variable over time. A binary nature would be inconsistent with the inferred ejection velocity from the Galactic disc of 420 km/s, which is too large to be survivable by a runaway binary system. A kinematic analysis of the star results in an age constraint (23 +- 1 Myr) that can be imposed on asteroseismic modelling and that can be used to remove degeneracies in the modelling process. Our attempts to match the excitation of the observed frequency spectrum resulted in pulsation models that were too young. Hence, asteroseismic studies of runaway pulsators can become vital not only in tracing the evolutionary history of such objects, but to understand the interior structure of massive stars in general. TESS is now opening up these stars for detailed asteroseismic investigation.Comment: accepted for ApJ

Layered convection as the origin of Saturn's luminosity anomaly

As they keep cooling and contracting, Solar System giant planets radiate more energy than they receive from the Sun. Applying the first and second principles of thermodynamics, one can determine their cooling rate, luminosity, and temperature at a given age. Measurements of Saturn's infrared intrinsic luminosity, however, reveal that this planet is significantly brighter than predicted for its age. This excess luminosity is usually attributed to the immiscibility of helium in the hydrogen-rich envelope, leading to "rains" of helium-rich droplets. Existing evolution calculations, however, suggest that the energy released by this sedimentation process may not be sufficient to resolve the puzzle. Here, we demonstrate using planetary evolution models that the presence of layered convection in Saturn's interior, generated, like in some parts of Earth oceans, by the presence of a compositional gradient, significantly reduces its cooling. It can explain the planet's present luminosity for a wide range of configurations without invoking any additional source of energy. This suggests a revision of the conventional homogeneous adiabatic interior paradigm for giant planets, and questions our ability to assess their heavy element content. This reinforces the possibility for layered convection to help explaining the anomalously large observed radii of extrasolar giant planets.Comment: Published in Nature Geoscience. Online publication date: April 21st, 2013. Accepted version before journal editing and with Supplementary Informatio