During most of their life, stars fuse hydrogen into helium in their cores. Mixing of chemical elements in the radiative envelope of stars with a convective

There are many Slowly Pulsating B (SPB) stars and γ Dor stars in the Kepler Mission data set. The light curves of these pulsating stars have been classified phenomenologically into stars with symmetric light curves and with asymmetric light curves. In the same effective temperature ranges as the γ Dor and SPB stars, there are variable stars with downward light curves that have been conjectured to be caused by spots. Among these phenomenological classes of stars, some show 'frequency groups' in their amplitude spectra that have not previously been understood. While it has been recognised that nonlinear pulsation gives rise to combination frequencies in a Fourier description of the light curves of these stars, such combination frequencies have been considered to be a only a minor constituent of the amplitude spectra. In this paper we unify the Fourier description of the light curves of these groups of stars, showing that many of them can be understood in terms of only a few base frequencies, which we attribute to g mode pulsations, and combination frequencies, where sometimes a very large number of combination frequencies dominate the amplitude spectra. The frequency groups seen in these stars are thus tremendously simplified. We show observationally that the combination frequencies can have amplitudes greater than the base frequency amplitudes, and we show theoretically how this arises. Thus for some γ Dor and SPB stars combination frequencies can have the highest observed amplitudes. Among the B stars are pulsating Be stars that show emission lines in their spectra from occasional ejection of material into a circumstellar disk. Our analysis gives strong support to the understanding of these pulsating Be stars as rapidly rotating SPB stars, explained entirely by g mode pulsations.

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.