2025ApJ...994..133M

NASA's TESS mission has unveiled a plethora of eclipsing binaries (EBs), among them hundreds of triples and higher-order, hierarchical systems. These complex targets require follow-up observations to enable full characterization of system architectures and identify the most compact multiples expected to undergo the most dramatic dynamical evolution. We report first results from a long-term effort to perform such follow-up, focusing here on multiband speckle imaging of a majority (57) of the sample of 97 quadruple- and higher-order eclipsing binaries (Q+EBs) identified via TESS light curves by V. B. Kostov et al. Diffraction-limited imaging with the Differential Speckle Survey Instrument on the Astrophysical Research Consortium 3.5 m telescope and HRCam on the Southern Astrophysical Research 4.1 m telescope reveals nearly 60% of the 57 to resolve into two sources separated by ≥0$\mathop{.}\limits^{^{\prime\prime} }$03. For these partly resolved systems, we report derived characteristics (e.g., relative position angle, angular separation, and magnitude differences in multiple passbands) from the speckle imaging. We find those Q+EBs partly resolved with 4 m class telescopes to have significantly inflated Gaia parallax errors and large Gaia renormalized unit weight errors, particularly for systems with separations comparable to Gaia's resolution limit (∼0$\mathop{.}\limits^{^{\prime\prime} }$6). For unresolved systems we report upper limits on angular and linear projected separations. We find two partly resolved Q+EBs with wide linear separations having eclipse timing variations that are therefore candidates of higher-than-quadruple multiplicity. Finally, we demonstrate how speckle imaging of resolved Q+EBs during an eclipse can clarify which speckle-resolved Q+EB subsystem is associated with a particular set of TESS eclipses.