2023ApJ...948...45I

A puzzling variety of superorbital modulations has been discovered in several supergiant high mass X-ray binaries (sgHMXBs). To investigate the mechanisms driving these superorbital modulations, we have analyzed long-term Neil Gehrels Swift Observatory (Swift) Burst Alert Telescope (BAT) observations of three sgHMXBs: 4U 1909 + 07, IGR J16418-4532, and IGR J16479-4514, and constructed their dynamic power spectra and superorbital intensity profiles. These Swift BAT observations are complemented by pointed Swift X-ray Telescope (XRT) and Nuclear Spectroscopic Telescope Array (NuSTAR) observations performed near the predicted maximum and minimum phase of a single superorbital cycle for each of these sources. The BAT dynamic power spectra show changes in the strength of the superorbital modulation on timescales of years, with either the peak at the fundamental frequency and/or the second harmonic present at different times for all three sources. The pointed Swift XRT and NuSTAR observations show no significant differences between the pulse profiles and spectral parameters at the superorbital maximum and minimum phase. This is likely due to the fact the superorbital modulation had weakened significantly during the times when the NuSTAR observations were carried out for all three sources. The results from the Swift XRT, BAT, and NuSTAR analysis indicate the possible presence of multiple corotating interaction regions (CIRs) in the stellar winds of the supergiant stars, although a structured stellar wind from the supergiant star due to tidal oscillations cannot be ruled out.