Dynamical Architectures of S-type Transiting Planets in Binaries. II. A Dichotomy in Orbital Alignment of Small Planets in Close Binary Systems

Abstract 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 <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mover> <mml:mrow> <mml:mo>.</mml:mo> </mml:mrow> <mml:mrow> <mml:mtext>″</mml:mtext> </mml:mrow> </mml:mover> </mml:math> 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 &lt; 5 R ⊕ ). 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 ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mi>log</mml:mi> <mml:mi>Z</mml:mi> <mml:mo>≳</mml:mo> <mml:mn>13.9</mml:mn> </mml:math> , or ≈5 σ ), indicating the observed planet-binary Δ I los values likely originate from two underlying populations: one nearly aligned ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msub> <mml:mrow> <mml:mi>σ</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>1</mml:mn> </mml:mrow> </mml:msub> <mml:mo>=</mml:mo> <mml:mn>2</mml:mn> <mml:mover> <mml:mrow> <mml:mo>.</mml:mo> </mml:mrow> <mml:mrow> <mml:mtext>°</mml:mtext> </mml:mrow> </mml:mover> <mml:msubsup> <mml:mrow> <mml:mn>4</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>0.9</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.7</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> ) and one with more scattered mutual inclinations ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:msub> <mml:mrow> <mml:mi>σ</mml:mi> </mml:mrow> <mml:mrow> <mml:mn>2</mml:mn> </mml:mrow> </mml:msub> <mml:mo>=</mml:mo> <mml:mn>23</mml:mn> <mml:mover> <mml:mrow> <mml:mo>.</mml:mo> </mml:mrow> <mml:mrow> <mml:mtext>°</mml:mtext> </mml:mrow> </mml:mover> <mml:msubsup> <mml:mrow> <mml:mn>6</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>−</mml:mo> <mml:mn>7.1</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>8.8</mml:mn> </mml:mrow> </mml:msubsup> </mml:math> ). 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 &gt;40 au, while systems with close-in or eccentric stellar companions (periastron distances &lt;40 au) preserve planet−binary alignment.