We present the detection and morphological characterization of hot molecular gas outflows in nearby ultra-luminous infrared galaxies, using the near-IR integral-field spectrograph SINFONI on the VLT. We detect outflows observed in the 2.12 micron H$_{2}$ 1-0 S(1) line for three out of four ULIRGs analyzed; IRAS 12112+0305, 14348-1447, and 22491-1808. The outflows are mapped on scales of 0.7-1.6 kpc, show typical outflow velocities of 300-500 km/s, and appear to originate from the nuclear region. The outflows comprise hot molecular gas masses of ~6-8x10$^3$ M(sun). Assuming a hot-to-cold molecular gas mass ratio of 6x10$^{-5}$, as found in nearby luminous IR galaxies, the total (hot+cold) molecular gas mass in these outflows is expected to be ~1x10$^{8}$ M(sun). This translates into molecular mass outflow rates of ~30-85 M(sun)/yr, which is a factor of a few lower than the star formation rate in these ULIRGs. In addition, most of the outflowing molecular gas does not reach the escape velocity of these merger systems, which implies that the bulk of the outflowing molecular gas is re-distributed within the system and thus remains available for future star formation. The fastest H$_{2}$ outflow is seen in the Compton-thick AGN of IRAS 14348-1447, reaching a maximum outflow velocity of ~900 km/s. Another ULIRG, IRAS 17208-0014, shows asymmetric H$_{2}$ line profiles different from the outflows seen in the other three ULIRGs. We discuss several alternative explanations for its line asymmetries, including a very gentle galactic wind, internal gas dynamics, low-velocity gas outside the disk, or two superposed gas disks. We do not detect the hot molecular counterpart to the outflow previously detected in CO(2-1) in IRAS 17208-0014, but we note that our SINFONI data are not sensitive enough to detect this outflow if it has a small hot-to-cold molecular gas mass ratio of < 9x10$^{-6}$.

The nuclei of merging galaxies are often deeply buried in dense layers of gas and dust. In these regions, gas outflows driven by starburst and AGN activity are believed to play a crucial role in the evolution of these galaxies. However, to fully understand this process it is essential to resolve the morphology and kinematics of such outflows. Using near-IR integral-field spectroscopy obtained with VLT/SINFONI, we detect a kpc-scale structure of high-velocity molecular hydrogen (H2) gas associated with the deeply buried secondary nucleus of the IR-luminous merger NGC3256. We show that this structure is likely the hot component of a molecular outflow, which is detected also in the cold molecular gas by Sakamoto et al. This outflow, with a molecular gas mass of M(H2)~2x10^7 Msun, is among the first to be spatially resolved in both the hot H2 gas with VLT/SINFONI and the cold CO-emitting gas with ALMA. The hot and cold components share a similar morphology and kinematics, with a hot-to-cold molecular gas mass ratio of ~6x10^-5. The high (~100 pc) resolution at which we map the geometry and velocity structure of the hot outflow reveals a biconical morphology with opening angle ~40 deg and gas spread across a FWZI~1200 km/s. Because this collimated outflow is oriented close to the plane of the sky, the molecular gas may reach maximum intrinsic outflow velocities of ~1800 km/s, with an average mass outflow rate of at least ~20 Msun/yr. By modeling the line-ratios of various near-IR H2 transitions, we show that the H2 gas in the outflow is heated through shocks or X-rays to a temperature of ~1900K. The energy needed to drive the outflow is likely provided by a hidden Compton-thick AGN or by the nuclear starburst. We show that the global kinematics of the molecular outflow in NGC3256 mimic those of CO-outflows that have been observed at low spatial resolution in starburst- and active galaxies.

We present a multi-wavelength integral field spectroscopic study of the low-z LIRG IRAS F11506-3851, on the basis of VIMOS and SINFONI (ESO-VLT) observations. The morphology and the 2D kinematics of the gaseous (neutral and ionized) and stellar components have been mapped using the NaD doublet, the H$\alpha$ line, and the near-IR CO(2-0) and CO(3-1) bands. The kinematics of the ionized gas and the stars are dominated by rotation, with large observed velocity amplitudes and centrally peaked velocity dispersion maps. The stars lag behind the warm gas and represent a dynamically hotter system, as indicated by the observed dynamical ratios. Thanks to these IFS data we have disentangled the contribution of the stars and the ISM to the NaD feature, finding that it is dominated by the absorption of neutral gas clouds in the ISM. The neutral gas 2D kinematics shows a complex structure dominated by two components. On the one hand, the thick slowly rotating disk lags significantly compared to the ionized gas and the stars, with an irregular and off-center velocity dispersion map. On the other hand, a kpc-scale neutral gas outflow is observed along the semi-minor axis of the galaxy, as revealed by large blueshifted velocities (30-154 km/s). We derive an outflowing mass rate in neutral gas of about 48 $\dot{M_{\rm w}}$/yr. Although this implies a global mass loading factor of 1.4, the 2D distribution of the ongoing SF suggests a much larger value of mass loading factor associated with the inner regions (R$<$200 pc), where the current SF represents only 3 percent of the total. All together these results strongly suggest that we are witnessing (nuclear) quenching due to SF feedback in IRAS F11506-3851. However, the relatively large mass of molecular gas detected in the nuclear region via the H2 1-0 S(1) line suggests that further episodes of SF may take place again.