Lack of high-precision long-term continuous photometric data for large samples of stars has prevented the large-scale exploration of pulsational variability in the OB star regime. As a result, the candidates for in-depth asteroseismic modelling remained limited to a few tens of dwarfs. The TESS nominal space mission has surveyed the southern sky, yielding continuous data of at least 27 d for hundreds of OB stars. We aim to couple TESS data in the southern sky with spectroscopy to study the variability over mass and evolution. We focus mainly on the presence of coherent pulsation modes that may or may not be present in the theoretical instability domains and unravel all frequency behaviour in the amplitude spectra of the TESS data. We compose a sample of 98 OB-type stars observed by TESS in Sectors 1-13 and with available high-resolution spectroscopy gathered by the IACOB and OWN surveys. We present the short-cadence 2-min light curves of dozens of OB-type stars, that have one or more spectra in the IACOB or OWN database. Based on these light curves and their Lomb-Scargle periodograms we perform variability classification and frequency analysis, and place the stars in the spectroscopic Hertzsprung-Russell diagram to interpret the variability in an evolutionary context. We deduce diverse origins of the variability found in all of the 98 OB stars in the TESS data. Among these we find several new variable stars, including three hybrid pulsators, three eclipsing binaries, high frequency modes in a Be star, and potential heat-driven pulsations in two Oe stars. We identify stars for which future asteroseismic modelling is possible, provided mode identification is achieved. By comparing the position of the variables to theoretical instability strips we discuss the current shortcomings in non-adiabatic pulsation theory, and the distribution of pulsators in the upper Hertzsprung-Russell diagram.

(Abridged) Eclipsing, spectroscopic double-lined binary star systems (SB2) are excellent laboratories for calibrating theories of stellar interior structure and evolution. We aim to investigate the mass discrepancy in binary stars. We study the effect of near-core mixing on the mass of the convective core of the stars and interpret the results in the context of the mass discrepancy. Two scenarios are considered, where individual stellar components of a binary system are treated independent of each other and where they are forced to have the same age and initial chemical composition. We find that the mass discrepancy is present in our sample and that it is anti-correlated with the surface gravity of the star. No correlations are found with other fundamental and atmospheric parameters, including the stellar mass. The mass discrepancy can be partially accounted for by increasing the amount of near-core mixing in stellar evolution models. We also find that ignoring the microturbulent velocity and turbulent pressure in stellar atmosphere models of hot evolved stars results in overestimation of their effective temperature by up to 8%. Together with enhanced near-core mixing, this can almost entirely account for the 30% mass discrepancy found for the evolved primary component of V380 Cyg. We find a strong link between the mass discrepancy and the convective core mass. The mass discrepancy can be solved by considering the combined effect of extra near-core boundary mixing and consistent treatment in the spectrum analysis of hot evolved stars. Our binary modelling results in convective core masses between 17 and 35% of the stellar mass, in excellent agreement with results from gravity-mode asteroseismology of single stars. This implies larger helium core masses near the end of the main sequence than anticipated so far.

Aims: We investigate from a theoretical perspective if space asteroseismology can be used to distinguish between different thermal structures and shapes of the near-core mixing profiles for different types of coherent oscillation modes in massive stars with convective cores, and if this capacity depends on the evolutionary stage of the models along the main sequence. Methods: We compute 1D stellar structure and evolution models for four different prescriptions of the mixing and temperature gradient in the near-core region. Their effect on the frequencies of dipole prograde gravity modes in both Slowly Pulsating B and $\beta$ Cep stars is investigated, as well as for pressure modes in $\beta$ Cep stars. Results: A comparison between the mode frequencies of the different models at various stages during the main sequence evolution reveals that they are more sensitive to a change in temperature gradient than to the exact shape of the mixing profile in the near-core region. Depending on the duration of the observed light curve, one can distinguish between either just the temperature gradient, or also between the shapes of the mixing coefficient. The relative frequency differences are in general larger for more evolved models, and are largest for the higher-frequency pressure modes in $\beta$ Cep stars. Conclusions:In order to unravel the core boundary mixing and thermal structure of the near-core region, one must have asteroseismic masses and radii with $\sim 1\%$ relative precision for hundreds of stars.

Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pertinent for the blue supergiant phase. Stellar oscillations represent a unique method of probing stellar interiors, yet inference for blue supergiants is hampered by a dearth of observed pulsation modes. Here we report the detection of diverse variability in blue supergiants using the K2 and TESS space missions. The discovery of pulsation modes or an entire spectrum of low-frequency gravity waves in these stars allow us to map the evolution of hot massive stars towards the ends of their lives. Future asteroseismic modelling will provide constraints on ages, core masses, interior mixing, rotation and angular momentum transport. The discovery of variability in blue supergiants is a step towards a data-driven empirical calibration of theoretical evolution models for the most massive stars in the Universe.

The evolution of a star is driven by the physical processes in its interior making the theory of stellar structure and evolution the most crucial ingredient for not only stellar evolution studies, but any field of astronomy which relies on the yields along stellar evolution. High-precision time-series photometric data assembled by recent space missions revealed that current models of stellar structure and evolution show major shortcomings already in the two earliest nuclear burning phases, impacting all subsequent phases prior to the formation of the end-of-life remnant. This white paper focuses specifically on the transport of chemical elements and of angular momentum in the stellar structure and evolution models of stars born with convective core, as revealed by their gravity-mode oscillations.

We propose a methodological framework to perform forward asteroseismic modeling of stars with a convective core, based on gravity-mode oscillations. These probe the near-core region in the deep stellar interior. The modeling relies on a set of observed high-precision oscillation frequencies of low-degree coherent gravity modes with long lifetimes and their observational uncertainties. Identification of the mode degree and azimuthal order is assumed to be achieved from rotational splitting and/or from period spacing patterns. This paper has two major outcomes. The first is a comprehensive list and discussion of the major uncertainties of theoretically predicted gravity-mode oscillation frequencies based on linear pulsation theory, caused by fixing choices of the input physics for evolutionary models. Guided by a hierarchy among these uncertainties of theoretical frequencies, we subsequently provide a global methodological scheme to achieve forward asteroseismic modeling. We properly take into account correlations amongst the free parameters included in stellar models. Aside from the stellar mass, metalicity and age, the major parameters to be estimated are the near-core rotation rate, the amount of convective core overshooting, and the level of chemical mixing in the radiative zones. This modeling scheme allows for maximum likelihood estimation of the stellar parameters for fixed input physics of the equilibrium models, followed by stellar model selection considering various choices of the input physics. Our approach uses the Mahalanobis distance instead of the often used $\chi^2$ statistic and includes heteroscedasticity. It provides estimation of the unknown variance of the theoretically predicted oscillation frequencies.