Interstellar complex organic molecules towards outflows from the G351.16+0.70 (NGC 6334 V) massive protostellar system

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 sources. 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 analysed the emission spectra in the 1mm band within 8GHz bandwidth coverage, finding that the observed emission spectra could be explained by 11 iCOMs and 5 molecular isotopologues. Additionally, we analysed the unusually 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 CH3OH and CH3OCHO present evidence of extended emission that does not fit with spherical morphology and that follows the extension of the low-velocity (CO)-C-13 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 an extremely high-velocity outflow (|V-Max - V-LSR| similar to 60kms(-1)). Finally, we provide the complete catalogue of iCOMs produced by the protostellar core within the 8GHz bandwidth of our observations, along with a compilation of iCOM abundances studied in interferometric works from the literature and a comparison with single-dish results.