Context. Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and rather high artificial viscosity and thermal diffusivity are needed for numerical stability. Direct comparison with observations is thus limited. Aims. We explore the capabilities of our compressible multidimensional hydrodynamics code SLH to simulate stellar oscillations. Methods. We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic $3\mathrm{M}_\odot$ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results. Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with observed spectra of massive stars. Conclusion. The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows the coupling of gravity and pressure waves to be studied too.

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.

We present the first three-dimensional (3D), hydrodynamic simulations of the core convection zone (CZ) and extended radiative zone spanning from 1% to 90% of the stellar radius of an intermediate mass (3 $\mathrm{M}_\odot$) star. This allows us to self-consistently follow the generation of internal gravity waves (IGWs) at the convective boundary and their propagation to the surface. We find that convection in the core is dominated by plumes. The frequency spectrum in the CZ and that of IGW generation is a double power law as seen in previous two-dimensional (2D) simulations. The spectrum is significantly flatter than theoretical predictions using excitation through Reynolds stresses induced by convective eddies alone. It is compatible with excitation through plume penetration. An empirically determined distribution of plume frequencies generally matches the one necessary to explain a large part of the observed spectrum. We observe waves propagating in the radiation zone and excited standing modes, which can be identified as gravity and fundamental modes. They show similar frequencies and node patterns to those predicted by the stellar oscillation code GYRE. The continuous part of the spectrum fulfills the IGW dispersion relation. A spectrum of tangential velocity and temperature fluctuations close to the surface is extracted, which are directly related to observable brightness variations in stars. Unlike 2D simulations we do not see the high frequencies associated with wave breaking, likely because these 3D simulations are more heavily damped.

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.

Internal gravity waves (IGW) propagate in the radiation zones of all stars. During propagation, their amplitudes are affected by two main features: radiative diffusion and density stratification. We have studied the implications of these two features on waves traveling within the radiative zones of non-rotating stars with stellar parameters obtained from the one dimensional stellar evolution code, MESA. As a simple measure of induced wave dynamics, we define a criterion to see if waves can become nonlinear and if so, under what conditions. This was done to understand the role IGW may play in angular momentum transport and mixing within stellar interiors. We find that the IGW generation spectrum, convective velocities and the strength of density stratification all play major roles in whether waves become nonlinear. With increasing stellar mass, there is an increasing trend in nonlinear wave energies. The trends with different metallicities and ages depend on the generation spectrum.

Stochastically-triggered photospheric light variations reaching $\sim$$40$ mmag peak-to-valley amplitudes have been detected in the O8Iaf supergiant V973 Scorpii as the outcome of two months of high-precision time-resolved photometric observations with the BRIght Target Explorer (BRITE) nanosatellites. The amplitude spectrum of the time series photometry exhibits a pronounced broad bump in the low-frequency regime ($\lesssim$$0.9$ d$^{-1}$) where several prominent frequencies are detected. A time-frequency analysis of the observations reveals typical mode lifetimes of the order of $5-10$ days. The overall features of the observed brightness amplitude spectrum of V973 Sco match well with those extrapolated from two-dimensional hydrodynamical simulations of convectively-driven internal gravity waves randomly excited from deep in the convective cores of massive stars. An alternative or additional possible source of excitation from a subsurface convection zone needs to be explored in future theoretical investigations.