2026AJ....171...77Z

Stellar multiplicity plays a crucial role in shaping planet formation and dynamical evolution. We present a survey of 54 TESS objects of interest (TOIs) within 300 pc that exhibit significant Hipparcos-Gaia astrometric accelerations. We identified 35 TOIs with stellar companions at projected separations between 0$\mathop{.}\limits^{\unicode{x02033}}$1 and 2″ (or 10-200 au). We also identified 12 TOIs that could host planetary-mass or brown dwarf companions, including six that are newly discovered. Furthermore, we perform 3D orbital characterization for 12 binaries hosting confirmed planets or planet candidates, allowing us to constrain the line-of-sight mutual inclination, ΔI_los, between the planetary and binary orbits. Combining our sample with previous measurements, we apply Bayesian hierarchical analysis to a total of 26 binary systems with S-type transiting planets (r_p < 5R_⊕). Specifically, we fit the ΔI_los distribution with both single (Rayleigh) and mixture models (two-component Rayleigh and Rayleigh-isotropic mixture). We find the mixture models are strongly favored ($\mathrm{log}Z\gtrsim 13.9$, or ≈5σ), indicating the observed planet-binary ΔI_los values likely originate from two underlying populations: one nearly aligned (${\sigma }_{1}=2\mathop{.}\limits^{\unicode{x000b0}}{4}_{-0.9}^{+0.7}$) and one with more scattered mutual inclinations (${\sigma }_{2}=23\mathop{.}\limits^{\unicode{x000b0}}{6}_{-7.1}^{+8.8}$). Alternatively, the misaligned systems can be equally well described by an isotropic distribution of inclinations. This observed dichotomy likely reflects different dynamical histories. Notably, the misaligned population only emerges in systems with stellar periastron distances >40 au, while systems with close-in or eccentric stellar companions (periastron distances <40 au) preserve planet-binary alignment.