Discovery of a Jupiter Analog Misaligned to the Inner Planetary System in HD 73344

We present the discovery of a Jupiter-like planet, HD 73344 d ($m_{d}=2.55^{+0.56}_{-0.46}\ \mathrm{M_{J}}$, $a_{d}=6.70^{+0.25}_{-0.26}$ AU, $e_{d}=0.18^{+0.14}_{-0.12}$) based on 27-year radial velocity observations. HD 73344 also hosts a compact inner planetary system, including a transiting sub-Neptune HD 73344 b ($P_{b}=15.61\ \mathrm{days}$, $r_{b}=2.88^{+0.08}_{-0.07}\ \mathrm{R_{\oplus}}$) and a non-transiting Saturn-mass planet ($P_{c}=65.936\ \mathrm{days}$, $m_{c}\sin{i_c}=0.36^{+0.02}_{-0.02}\ \mathrm{M_{J}}$). By analyzing TESS light curves, we identified a stellar rotation period of $9.03\pm{1.3}$ days. Combining this with $v\sin{i_*}$ measurements from stellar spectra, we derived a stellar inclination of $63.6^{+17.4}_{-16.5}\ \rm{deg}$. Furthermore, by combining radial velocities and Hipparcos-Gaia astrometric acceleration, we characterized the three-dimensional orbit of the outer giant planet and constrained its mutual inclination relative to the innermost transiting planet to be $46 <\Delta I_{bd}< 134\ \rm{deg}\ (1\sigma)$ and $20 <\Delta I_{bd}< 160\ \rm{deg}\ (2\sigma)$, strongly disfavoring coplanar architectures. Our analytical calculations and N-body simulation reveal that the two inner planets are strongly coupled with each other and undergo nodal precession together around the orbital axis of the giant planet. During nodal precession, the orbital inclination of inner planets oscillate with time and therefore become misaligned relative to the stellar spin axis. The formation of such systems suggests a history of planet-planet scattering or misalignment between the inner and outer components of protoplanetary disks. The upcoming release of Gaia DR4 will uncover more systems similar to HD 73344 and enable the study of the flatness of exoplanet systems with a mixture of inner and outer planetary systems on a statistical level.