Galactic outflows are known to consist of several gas phases, however, so far the connection between these multiple phases has been investigated little and only in a few objects. In this paper, we analyse MUSE/VLT data of 26 local (U)LIRGs and study their ionised and neutral atomic phases. We also include objects from the literature to obtain a total sample of 31 galaxies with spatially resolved multi-phase outflow information. We find that the ionized phase of the outflows has on average an electron density three times higher than the disc ($n_{\rm e, disc}$ $\sim$ 150 cm$^{-3}$ vs $n_{\rm e, outflow}$ $\sim$ 500 cm$^{-3}$), suggesting that cloud compression in the outflow is more important that cloud dissipation. We find that the difference in extinction between outflow and disc correlates with the outflow gas mass. Together with the analysis of the outflow velocities, this suggests that at least some of the outflows are associated with the ejection of dusty clouds from the disc. The presence of dust in outflows is relevant for potential formation of molecules inside them. We combine our data with millimetre data to investigate the molecular phase. We find that the molecular phase accounts for more than 60 $\%$ of the total mass outflow rate in most objects and this fraction is higher in AGN-dominated systems. The neutral atomic phase contributes of the order of 10 $\%$, while the ionized phase is negligible. The ionized-to-molecular mass outflow rate declines slightly with AGN luminosity, although with a large scatter.

Ionized outflows are ubiquitous in non radio-loud obscured quasars (QSO2s) at different redshifts. We search for large-scale ionized outflows associated with six optically selected QSO2 (five non-radio-loud and one radio-loud) at $z\sim$ 0.2-0.5, targeting objects with extended radio structures. We have obtained OSIRIS/GTC optical long slit spectroscopy data for these six QSO2 with the slit located along the radio axis. We traced the gas kinematics with the [OIII]$\lambda$4959,5007 lines to investigate ionized outflows and characterize the dynamical state of the host galaxies. This second study has been complemented with previously published FORS2/VLT spectroscopic data of 13 more QSO2 at similar z. We identify ionized outflows in four out of the six QSO2 observed with the GTC. The outflows are spatially unresolved in two QSO2 and compact in a third (radial size of R=0.8$\pm$0.3 kpc). Of particular interest is the radio-quiet QSO2 SDSS 0741+3020 at z=0.47, associated with a giant $\sim$ 112 kpc nebula. An ionized outflow probably induced by the radio structures has been detected along the axis defined by the central $\sim$1\arcsec\ radio structure, extending up to at least $\sim$4 kpc from the active galactic nucleus (AGN). Turbulent gas ($\sigma\sim$130 km s$^{-1}$) has also been detected across the giant gas nebula up to $\sim$40 kpc from the AGN. This turbulence may have been induced by outflows triggered by the interaction between a so-far undetected large-scale radio source and the nebula. Regarding the dynamical state of the host galaxies, we find that the majority of the QSO2 show v/$\sigma<$ 1, implying that they are dominated by random motions (so-called dispersion-dominated systems). Most (17 of 19) fall in the area of the E/S0 galaxies in the dynamical diagram v/$\sigma$ versus $\sigma$. None are consistent with spiral or disk galaxies.

We study how the use of several stellar subtraction methods and line fitting approaches can affect the derivation of the main kinematic parameters (velocity and velocity dispersion fields) of the ionized gas component. The target of this work is the nearby galaxy NGC 2906, observed with the MUSE instrument at Very Large Telescope. A sample of twelve spectra is selected from the inner (nucleus) and outer (spiral arms) regions, characterized by different ionization mechanisms. We compare three different methods to subtract the stellar continuum (FIT3D, STARLIGHT and pPXF), combined with one of the following stellar libraries: MILES, STELIB and GRANADA+MILES. The choice of the stellar subtraction method is the most important ingredient affecting the derivation of the gas kinematics, followed by the choice of the stellar library and by the line fitting approach. In our data, typical uncertainties in the observed wavelength and width of the H\alpha and [NII] lines are of the order of <\delta\lambda>_rms \sim 0.1\AA\ and <\delta\sigma>_rms \sim 0.2\AA\ (\sim 5 and 10km/s, respectively). The results obtained from the [NII] line seem to be slightly more robust, as it is less affected by stellar absorption than H\alpha. All methods considered yield statistically consistent measurements once a mean systemic contribution \Delta\bar\lambda=\Delta\bar\sigma=0.2xDelta_{MUSE} is added in quadrature to the line fitting errors, where \Delta_{MUSE} = 1.1\AA\ \sim 50 km/s denotes the instrumental resolution of the MUSE spectra. Although the subtraction of the stellar continuum is critical in order to recover line fluxes, any method (including none) can be used in order to measure the gas kinematics, as long as an additional component of 0.2 x Delta_MUSE is added to the error budget.

Recent observations have revealed massive galactic molecular outflows that may have physical conditions (high gas densities) required to form stars. Indeed, several recent models predict that such massive galactic outflows may ignite star formation within the outflow itself. This star-formation mode, in which stars form with high radial velocities, could contribute to the morphological evolution of galaxies, to the evolution in size and velocity dispersion of the spheroidal component of galaxies, and would contribute to the population of high-velocity stars, which could even escape the galaxy. Such star formation could provide in-situ chemical enrichment of the circumgalactic and intergalactic medium (through supernova explosions of young stars on large orbits), and some models also predict that it may contribute substantially to the global star formation rate observed in distant galaxies. Although there exists observational evidence for star formation triggered by outflows or jets into their host galaxy, as a consequence of gas compression, evidence for star formation occurring within galactic outflows is still missing. Here we report new spectroscopic observations that unambiguously reveal star formation occurring in a galactic outflow at a redshift of 0.0448. The inferred star formation rate in the outflow is larger than 15 Msun/yr. Star formation may also be occurring in other galactic outflows, but may have been missed by previous observations owing to the lack of adequate diagnostics.

This paper describes the Bayesian Technique for Multi-image Analysis (BaTMAn), a novel image-segmentation technique based on Bayesian statistics that characterizes any astronomical dataset containing spatial information and performs a tessellation based on the measurements and errors provided as input. The algorithm iteratively merges spatial elements as long as they are statistically consistent with carrying the same information (i.e. identical signal within the errors). We illustrate its operation and performance with a set of test cases including both synthetic and real Integral-Field Spectroscopic data. The output segmentations adapt to the underlying spatial structure, regardless of its morphology and/or the statistical properties of the noise. The quality of the recovered signal represents an improvement with respect to the input, especially in regions with low signal-to-noise ratio. However, the algorithm may be sensitive to small-scale random fluctuations, and its performance in presence of spatial gradients is limited. Due to these effects, errors may be underestimated by as much as a factor of two. Our analysis reveals that the algorithm prioritizes conservation of all the statistically-significant information over noise reduction, and that the precise choice of the input data has a crucial impact on the results. Hence, the philosophy of BaTMAn is not to be used as a `black box' to improve the signal-to-noise ratio, but as a new approach to characterize spatially-resolved data prior to its analysis. The source code is publicly available at http://astro.ft.uam.es/SELGIFS/BaTMAn .

The kinematic characterization of different galaxy populations is a key observational input to distinguish between different galaxy evolutionary scenarios, since it helps to determine the number ratio of rotating disks to mergers at different cosmic epochs. Local (U)LIRGs offer a unique opportunity to study at high linear resolution and S/N extreme star forming events and compare them with those observed at high-z. We obtained Very Large Telescope (VLT) VIMOS optical integral field spectroscopy (IFS) data of a sample of 38 (U)LIRGs. The `unweighted' and `weighted' {kinemetry}-based methods are used to kinematically classify our galaxies in `disk' and `merger'. We simulate our systems at z=3 to evaluate how a loss of angular resolution affects our results. From the kinemetry-based analysis we are able classify our local (U)LIRGs in three distinct kinematic groups according to their total kinematic asymmetry values (K$_{tot}$) as derived when using the weighted (unweighted) method: 1) 25 out of 50 galaxies are kinematically classified as `disk', with a K$_{tot}\leq$ 0.16 (0.14); 2) 9 out of 50 galaxies are kinematically classified as `merger', with a K$_{tot}\geq$ 0.94 (0.66); 3) 16 out of 50 galaxies lie in the `transition region', in which disks and mergers coexist, with 0.16 (0.14) $<$ K$_{tot}$ $<$ 0.94 (0.66). When we apply our criteria to the high-z simulated systems, a lower total kinematic asymmetry frontier value (K$_{tot}$ $\sim$ 0.16 ($\sim$ 0.14)) is derived with respect to that found locally. The loss of angular resolution smears out the kinematic features, thus making objects to appear more kinematically regular than actually they are.

We study the physical and kinematic properties of the narrow line region (NLR) of the nearest obscured quasar MRK 477 (z=0.037), using optical and near-infrared spectroscopy. We explore a diversity of aspects that provide a more complete understanding of the nature of this object, example of a type 2 quasar in the nearby Universe, as well as a starburst-AGN hybrid system [abridged].

We present the analysis of the integrated spectral energy distribution (SED) from the ultraviolet (UV) to the far-infrared and H$\alpha$ of a sample of 29 local systems and individual galaxies with infrared (IR) luminosities between 10^11 Lsun and 10^11.8 Lsun. We have combined new narrow-band H$\alpha$+[NII] and broad-band g, r optical imaging taken with the Nordic Optical Telescope (NOT), with archival GALEX, 2MASS, Spitzer, and Herschel data. The SEDs (photometry and integrated H$\alpha$ flux) have been fitted with a modified version of the MAGPHYS code using stellar population synthesis models for the UV-near-IR range and thermal emission models for the IR emission taking into account the energy balance between the absorbed and re-emitted radiation. From the SED fits we derive the star-formation histories (SFH) of these galaxies. For nearly half of them the star-formation rate appears to be approximately constant during the last few Gyrs. In the other half, the current star-formation rate seems to be enhanced by a factor of 3-20 with respect to that occured ~1 Gyr ago. Objects with constant SFH tend to be more massive than starbursts and they are compatible with the expected properties of a main-sequence (M-S) galaxy. Likewise, the derived SFHs show that all our objects were M-S galaxies ~1 Gyr ago with stellar masses between 10^10.1 and 10^11.5 Msun. We also derived from our fits the average extinction (A_v=0.6-3 mag) and the polycyclic aromatic hydrocarbons (PAH) luminosity to L(IR) ratio (0.03-0.16). We combined the A_v with the total IR and H$\alpha$ luminosities into a diagram which can be used to identify objects with rapidly changing (increasing or decreasing) SFR during the last 100 Myr.

We study the kinematic properties of the ambient ionized ISM and ionized gas outflows in a large and representative sample of local luminous and ultraluminous infrared galaxies (U/LIRGs) (58 systems, 75 galaxies), on the basis of integral field spectroscopy (IFS)-based high S/N integrated spectra at galactic and sub-galactic, i.e. star forming (SF) clumps, scales. Ambient ionized gas. The velocity dispersion of the ionized ISM in U/LIRGs (<sigma> ~ 70 kms-1) is larger than in lower luminosity local star forming galaxies (<sigma> ~ 25 kms-1). While for isolated disc LIRGs star formation appears to sustain turbulence, gravitational energy release associated to interactions and mergers plays an important role driving sigma in the U/LIRG range. We also find that the impact of an AGN in ULIRGs is strong, increasing sigma by a factor 1.5 on average. The observed weak dependency of sigma with SFR surface density for local U/LIRGs is in very good agreement with that measured in some high-z samples. Ionized outflows. The presence of ionized gas outflows in U/LIRGs seems universal based on the detection of a broad, usually blueshifted, Halpha line. AGNs in U/LIRGs are able to generate faster (x2) and more massive (x1.4) ionized gas outflows than pure starbursts. The derived ionized mass loading factors are in general below one, with only a few AGNs above this limit. Only a small fraction of the ionized material from low mass LIRGs (log(Mdyn/Msun) < 10.4) could reach the intergalactic medium, with more massive galaxies retaining the gas. The observed average outflow properties in U/LIRGs are similar to high-z galaxies of comparable SFR. In the bright SF clumps found in LIRGs, ionized gas outflows appear to be very common. For a given SFR surface density, outflows in LIRG clumps would be about one to two orders of magnitude less energetic than those launched by clumps in high-z SF 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.

We present and discuss the 2D kinematic properties of the ionized gas (Halpha) in a sample of 38 local (ultra) luminous infrared galaxies [(U)LIRGs] (31 LIRGs and 7 ULIRGs) observed with VIMOS at the VLT using integral field spectroscopy. This sample covers well the less studied LIRG luminosity range and includes isolated disks, interacting systems, and mergers. The majority of the galaxies have two main kinematically distinct components. One component (i.e., narrow or systemic) extends over the whole line-emitting region and is characterized by small to intermediate velocity dispersions (i.e., sigma from 30 to 160 km s^-1). It traces the overall velocity field. The second component (broad) has in general a larger velocity dispersion (up to 320 km s^-1), mainly found in the inner regions and generally blueshifted with respect to the systemic component. Most of the objects (76%) are dominated by rotation, more relevant in LIRGs than in ULIRGs. Isolated disks, interacting galaxies, and merging systems define a sequence of increasing mean velocity dispersion, and decreasing velocity field amplitude.The LIRGs classified as isolated disks have similar velocity amplitudes but larger mean velocity dispersions (44 vs. 24 km s^-1) than local spirals, implying a larger turbulence and thicker disks. Interacting systems and mergers have values closer to those of low velocity dispersion ellipticals/lenticular galaxies (E/SOs). The (U)LIRGs classified as mergers have kinematic properties similar to those shown by the Lyman break analogs (LBAs). The dynamical masses range from \sim 0.04 m* to 1.4 m* (i.e., m* = 1.4x10^{11} Msun), with ULIRGs (M{dyn} sim 0.5 +/- 0.2 m*) being more massive than LIRGs by, on average, a factor of about 2. The mass ratio of individual pre-coalescence galaxies is <2.5 for most of the systems, confirming that most (U)LIRG mergers involve sub-m* galaxies of similar mass.

We have obtained VIMOS/VLT optical integral field spectroscopy (IFS) data for a sample of 4 LIRGs which have been selected at a similar distance ($\sim$ 70 Mpc) to avoid relative resolution effects. They have been classified in two groups (isolated disk and post-coalescence mergers) according to their morphology. The $kinemetry$ method (developed by Krajnovic and coworkers) is used to characterize the kinematic properties of these galaxies and to discuss new criteria for distinguishing their status. We present and discuss new kinematic maps (i.e., velocity field and velocity dispersion) for these four galaxies. The morphological and kinematic classifications of these systems are consistent, with disks having lower kinematic asymmetries than post-coalescence mergers. We then propose and discuss a new kinematic criterion to differentiate these two groups. This criterion distinguishes better these two categories and has the advantage of being less sensitive to angular resolution effects. According to the previous criteria,the present post-coalescence systems would have been classified as disks. This indicates that the separation of disks from mergers is subjective to the definition of `merger'. It also suggests that previous estimates of the merger/disk ratio could have been underestimated, but larger samples are necessary to establish a firmer conclusion.

The "Extreme starbursts in the local universe" workshop was held at the Insituto de Astrofisica de Andalucia in Granada, Spain on 21-25 June 2010. Bearing in mind the advent of a new generation of facilities such as JWST, Herschel, ALMA, eVLA and eMerlin, the aim of the workshop was to bring together observers and theorists to review the latest results. The purpose of the workshop was to address the following issues: what are the main modes of triggering extreme starbursts in the local Universe? How efficiently are stars formed in extreme starbursts? What are the star formation histories of local starburst galaxies? How well do the theoretical simulations model the observations? What can we learn about starbursts in the distant Universe through studies of their local counterparts? How important is the role of extreme starbursts in the hierarchical assembly of galaxies? How are extreme starbursts related to the triggering of AGN in the nuclei of galaxies? Overall, 41 talks and 4 posters with their corresponding 10 minutes short talks were presented during the workshop. In addition, the workshop was designed with emphasis on discussions, and therefore, there were 6 discussion sessions of up to one hour during the workshop. Here is presented a summary of the purposes of the workshop as well as a compilation of the abstracts corresponding to each of the presentations. The summary and conclusions of the workshop along with a description of the future prospects by Sylvain Veilleux can be found in the last section of this document. A photo of the assistants is included.

We combine near-to-mid-IR Spitzer data with shorter wavelength observations (optical to X-rays) to get insights on the properties of a sample of luminous, obscured Active Galactic Nuclei (AGN). We aim at modeling their broad-band Spectral Energy Distributions (SEDs) in order to estimate the main parameters related to the dusty torus. The sample comprises 16 obscured high-redshift (0.9<z<2.1) xray luminous quasars (L_2-10 ~ 10^44 erg s-1) selected from the HELLAS2XMM survey. The SEDs are described by a multi-component model including a stellar component, an AGN component and a starburst. The majority (~80%) of the sources show moderate optical depth (tau_9.7um<3) and the derived column densities N_H are consistent with the xray inferred values (10^22 <N_H< 3x10^23 cm-2) for most of the objects, confirming that the sources are moderately obscured Compton-thin AGN. Accretion luminosities in the range 5x10^44 < Lbol < 4x10^46 erg s-1 are inferred. We compare model luminosities with those obtained by integrating the observed SED, finding that the latter are lower by a factor of ~2 in the median. The discrepancy can be as high as an order of magnitude for models with high optical depth (tau_9.7um=10). The ratio between the luminosities obtained by the fitting procedure and from the observed SED suggest that, at least for Type~2 AGN, observed bolometric luminosities are likely to underestimate intrinsic ones and the effect is more severe for highly obscured sources. Bolometric corrections from the hard X-ray band are computed and have a median value of k_2-10kev ~ 20. The obscured AGN in our sample are characterized by relatively low Eddington ratios (median lambda_Edd~0.08). On average, they are consistent with the Eddington ratio increasing at increasing bolometric correction (e.g. Vasudevan & Fabiam 2009).

We present multi-wavelength observations (from optical to sub-millimeter, including Spitzer and SCUBA) of H2XMMJ 003357.2-120038 (also GD158_19), an X-ray selected, luminous narrow-line (Type 2) quasar at z=1.957 selected from the HELLAS2XMM survey. Its broad-band properties can be reasonably well modeled assuming three components: a stellar component to account for the optical and near-IR emission, an AGN component (i.e., dust heated by an accreting active nucleus), dominant in the mid-IR, with an optical depth at 9.7 micron along the line of sight (close to the equatorial plane of the obscuring matter) of tau(9.7)=1 and a full covering angle of the reprocessing matter (torus) of 140 degrees, and a far-IR starburst component (i.e., dust heated by star formation) to reproduce the wide bump observed longward of 70 micron. The derived star-formation rate is about 1500 solar masses per year. The overall modeling indicates that GD158_19 is a high-redshift X-ray luminous, obscured quasar with coeval powerful AGN activity and intense star formation. It is probably caught before the process of expelling the obscuring gas has started, thus quenching the star formation.