We present VLT/SINFONI imaging spectroscopy of the warm ionized gas in 33 powerful radio galaxies at redshifts z>~2, which are excellent sites to study the interplay of rapidly accreting active galactic nuclei and the interstellar medium of the host galaxy in the very late formation stages of massive galaxies. Our targets span two orders of magnitude in radio size (2-400 kpc) and kinetic jet energy (a few 10^46 to almost 10^48 erg s^-1). All sources have complex gas kinematics with broad line widths up to ~1300 km s^-1. About half have bipolar velocity fields with offsets up to 1500 km s^-1 and are consistent with global back-to-back outflows. The others have complex velocity distributions, often with multiple abrupt velocity jumps far from the nucleus of the galaxy, and are not associated with a major merger in any obvious way. We present several empirical constraints that show why gas kinematics and radio jets seem to be physically related. The gas kinetic energy from large scale bulk and local outflow or turbulent motion corresponds to a few 10^-3 to 10^-2 of the kinetic energy of the jet. In galaxies with jet power >~10^47 erg s^-1, the bulk kinetic energy dominates the total energy budget of the gas, suggesting that the outflows encompasses the global interstellar medium, perhaps facilitated by the strong gas turbulence. We compare with recent hydrodynamic simulations, and discuss the potential consequences for the subsequent evolution of massive high-z galaxies. The gas-phase metallicities in our galaxies are lower than in most low-z AGN, but nonetheless solar or even super-solar, suggesting that the ISM in these galaxies is very similar to the gas from which massive low-redshift galaxies formed most of their gas. This further highlights that we are seeing these galaxies near the end of their active formation phase.

Most successful galaxy formation scenarios now postulate that the intense star formation in massive, high-redshift galaxies during their major growth period was truncated when powerful AGNs launched galaxy-wide outflows of gas that removed large parts of the interstellar medium. The most powerful radio galaxies at z~2 show clear signatures of such winds, but are too rare to be good representatives of a generic phase in the evolution of all massive galaxies at high redshift. Here we present SINFONI imaging spectroscopy of 12 radio galaxies at z~2 that are intermediate between the most powerful radio and vigorous starburst galaxies in radio power, and common enough to represent a generic phase in the early evolution of massive galaxies. The kinematic properties are diverse, with regular velocity gradients with amplitudes of Delta v=200-400 km s^-1 as in rotating disks as well as irregular kinematics with multiple velocity jumps of a few 100 km s^-1. Line widths are generally high, typically around FWHM=800 km s^-1, consistent with wind velocities in hydrodynamic models. A broad H-alpha line in one target implies a black hole mass of a few 10^9 M$_sun. The ratio of line widths, sigma, to bulk velocity, v, is so large that even the gas in galaxies with regular velocity fields is unlikely to be gravitationally bound. It is unclear, however, whether the large line widths are due to turbulence or unresolved, local outflows as are sometimes observed at low redshifts. Comparison of the kinetic energy with the energy supply from the AGN through jet and radiation pressure suggests that the radio source still plays a dominant role for feedback, consistent with low-redshift radio-loud quasars.

We report the detection of extended warm ionized gas in two powerful high-redshift radio galaxies, NVSS J210626-314003 at z=2.10 and TXS 2353-003 at z=1.49, that does not appear to be associated with the radio jets. This is contrary to what would be expected from the alignment effect, a characteristic feature of distant, powerful radio galaxies at z> 0.6. The gas also has smaller velocity gradients and line widths than most other high-z radio galaxies with similar data. Both galaxies are part of a systematic study of 50 high-redshift radio galaxies with SINFONI, and are the only two that are characterized by the presence of high surface-brightness gas not associated with the jet axis and by the absence of such gas aligned with the jet. Both galaxies are spatially resolved with ISAAC broadband imaging covering the rest-frame R band, and have extended wings that cannot be attributed to line contamination. We argue that the gas and stellar properties of these galaxies are more akin to gas-rich brightest cluster galaxies in cool-core clusters than the general population of high-redshift radio galaxies at z>2. In support of this interpretation, one of our sources, TXS 2353-003, for which we have H\alpha\ narrowband imaging, is associated with an overdensity of candidate H\alpha\ emitters by a factor of 8 relative to the field at z=1.5. We discuss possible scenarios of the evolutionary state of these galaxies and the nature of their emission line gas within the context of cyclical AGN feedback.

We present an algorithm for the detection of Low Surface Brightness (LSB) galaxies in images, called MARSIAA (MARkovian Software for Image Analysis in Astronomy), which is based on multi-scale Markovian modeling. MARSIAA can be applied simultaneously to different bands. It segments an image into a user-defined number of classes, according to their surface brightness and surroundings - typically, one or two classes contain the LSB structures. We have developed an algorithm, called DetectLSB, which allows the efficient identification of LSB galaxies from among the candidate sources selected by MARSIAA. To assess the robustness of our method, the method was applied to a set of 18 B and I band images (covering 1.3 square degrees in total) of the Virgo cluster. To further assess the completeness of the results of our method, both MARSIAA, SExtractor, and DetectLSB were applied to search for (i) mock Virgo LSB galaxies inserted into a set of deep Next Generation Virgo Survey (NGVS) gri-band subimages and (ii) Virgo LSB galaxies identified by eye in a full set of NGVS square degree gri images. MARSIAA/DetectLSB recovered ~20% more mock LSB galaxies and ~40% more LSB galaxies identified by eye than SExtractor/DetectLSB. With a 90% fraction of false positives from an entirely unsupervised pipeline, a completeness of 90% is reached for sources with r_e > 3" at a mean surface brightness level of mu_g=27.7 mag/arcsec^2 and a central surface brightness of mu^0 g=26.7 mag/arcsec^2. About 10% of the false positives are artifacts, the rest being background galaxies. We have found our method to be complementary to the application of matched filters and an optimized use of SExtractor, and to have the following advantages: it is scale-free, can be applied simultaneously to several bands, and is well adapted for crowded regions on the sky.

MUSE, the Multi Unit Spectroscopic Explorer, is a 2nd generation integral-field spectrograph under final assembly to see first light at the Very Large Telescope in 2013. By capturing ~ 90000 optical spectra in a single exposure, MUSE represents a challenge for data reduction and analysis. We summarise here the main features of the Data Reduction System, as well as some of the tools under development by the MUSE consortium and the DAHLIA team to handle the large MUSE datacubes (about 4x?10^8 pixels) to recover the original astrophysical signal.