Title: Interstellar Complex Organic Molecules towards outflows from the G351.16+0.70 (NGC 6334 V) massive protostellar system

Abstract: G351.16+0.70 is a relatively well-studied High Mass Star Forming Region with at least two main bipolar outflow structures originating from an OB embedded star and multiple IR cores. Using high-resolution and large bandwidth SMA observations, we studied its molecular content to probe the emission of iCOMs which could be related to the bipolar outflows or their jets. We analyzed the emission spectra in the 1mm band within 8 GHz bandwidth coverage, from 216.75 to 220.75 GHz, and from 228.75 to 232.75 GHz. Employing the LTE approximation using the XCLASS software, we identified 260 emission lines arising from iCOMs. The emission lines in the synthetic emission spectra could be explained by 11 iCOMs and 5 molecular isotopologues. Additionally, we analyzed the outstanding broad iCOM emission lines by using integrated and velocity field maps, searching for extended emission and velocity gradients related to molecular outflows. Ro-vibrational transitions of CH$_3$OH, C$_2$H$_3$CN, CH$_3$OCHO, CH$_3$COCH$_3$, and aGg'-(CH$_2$OH)$_2$ present evidence of extended emission that does not fit with spherical morphology and that follows the path of the low-velocity $^{13}$CO outflow. The multiple outflows in the system are revealed also by the CO (2-1) and SiO (5-4) emission, but in particular we have discovered a so-called extremely high velocity outflow ($|V_{Max}-V_{LSR}|\sim 60$ km s$^{-1}$). Additionally, we provide the full line catalog of iCOMs along the 8 GHz bandwidth produced by the main protostellar core.