Spectral Variability of VHS J1256-1257b from 1 to 5 μm

Multiwavelength time-resolved observations of rotationally modulated variability from brown dwarfs and giant exoplanets are the most effective method for constraining their heterogeneous atmospheric structures. In a companion paper, we reported the discovery of strong near-infrared variability in HST/WFC3/G141 light curves of the very red L-dwarf companion VHS J1256-1257b. In this paper, we present a follow-up 36 hr Spitzer/IRAC Channel 2 light curve together with an in-depth analysis of the Hubble space telescope (HST) and the Spitzer data. The combined data set provides time-resolved light curves of VHS1256b sampling 1.1 to 4.5 μm. The Spitzer light curve is best fit with a single sine wave with a period of 22.04 ± 0.05 hr and a peak-to-peak amplitude of 5.76 ± 0.04%. Combining the period with a previously measured projected rotational velocity ( $v\,\sin \,i$ ), we find that VHS1256b is most consistent with equatorial viewing geometry. The HST/G141+Spitzer spectral energy distribution favors a model with a ${T}_{\mathrm{eff}}$ of 1000 K and low surface gravity with disequilibrium chemistry. The spectral variability of VHS1256b is consistent with predictions from partly cloudy models, suggesting that heterogeneous clouds are the dominant source of the observed modulations. We find evidence at the 3.3σ level for amplitude variations within the 1.67 μm CH₄ band, which is the first such detection for a variable L-dwarf. We compare the HST/G141 time-resolved spectra of three red L-dwarfs with high-amplitude near-infrared rotational modulations and find that although their time-averaged spectra are similar, their spectroscopic variabilities exhibit notable differences. This diversity reinforces the advantage of time-resolved spectroscopic observations for understanding the atmospheres of brown dwarfs and directly imaged exoplanets.