Context. Massive stars are predicted to excite internal gravity waves (IGWs) by turbulent core convection and from turbulent pressure fluctuations in their near-surface layers. These IGWs are extremely efficient at transporting angular momentum and chemical species within stellar interiors, but they remain largely unconstrained observationally. Aims. We aim to characterise the photometric detection of IGWs across a large number of O and early-B stars in the Hertzsprung-Russell diagram, and explain the ubiquitous detection of stochastic variability in the photospheres of massive stars. Methods. We combined high-precision time-series photometry from the NASA Transiting Exoplanet Survey Satellite with high-resolution ground-based spectroscopy of 70 stars with spectral types O and B to probe the relationship between the photometric signatures of IGWs and parameters such as spectroscopic mass, luminosity, and macroturbulence. Results. A relationship is found between the location of a star in the spectroscopic Hertzsprung-Russell diagram and the amplitudes and frequencies of stochastic photometric variability in the light curves of massive stars. Furthermore, the properties of the stochastic variability are statistically correlated with macroturbulent velocity broadening in the spectral lines of massive stars. Conclusions. The common ensemble morphology for the stochastic low-frequency variability detected in space photometry and its relationship to macroturbulence is strong evidence for IGWs in massive stars, since these types of waves are unique in providing the dominant tangential velocity field required to explain the observed spectroscopy.

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.