Processing of raw data from modern astronomical instruments is nowadays often carried out using dedicated software, so-called "pipelines" which are largely run in automated operation. In this paper we describe the data reduction pipeline of the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph operated at ESO's Paranal observatory. This spectrograph is a complex machine: it records data of 1152 separate spatial elements on detectors in its 24 integral field units. Efficiently handling such data requires sophisticated software, a high degree of automation and parallelization. We describe the algorithms of all processing steps that operate on calibrations and science data in detail, and explain how the raw science data gets transformed into calibrated datacubes. We finally check the quality of selected procedures and output data products, and demonstrate that the pipeline provides datacubes ready for scientific analysis.

The merging of supermassive black holes (SMBHs) is a direct consequence of our hierarchical picture of galaxy evolution. It is difficult to track the merging process of SMBHs during mergers of galaxies as SMBHs are naturally difficult to observe. We want to characterise and confirm the presence of two independent active galactic nuclei (AGN) separated by a few kiloparsec in seven strongly interacting galaxies previously selected from the Sloan Digital Sky Survey (SDSS) as Seyfert-Seyfert pairs based on emission-line ratio diagnostics. Optical slit spectra taken with MODS at the Large Binocular Telescope (LBT) are presented to infer the detailed spatial distribution of optical emission lines, and their line ratios and AGN signatures with respect to the host galaxies, thereby quantifying the impact of beam smearing and large fibre apertures on the spectra captured by the SDSS. We find that at most two of the seven targets actually retain a Seyfert-Seyfert dual AGN, whereas the others may be more likely powered by post-AGB stars in retired galaxies or through shocks in the ISM based on spatially resolved optical line diagnostics. The major cause of this discrepancy is a bias caused by the spillover of flux from the primary source in the secondary SDSS fibre which can be more than an order of magnitude at <3" separations. Previously reported extremely low X-ray-to-[\ion{O}{iii}] luminosity ratios may be explained by this misclassification, as can heavily obscured AGN for the primaries. We also find that the nuclei with younger stellar ages host the primary AGN. Studies of close dual AGN selected solely from fibre-based spectroscopy can create severe biases in the sample selection and interpretation of the results. Spatially resolved spectroscopy should ideally be used in the future to characterise such compact systems together with multi-wavelength follow-up observations.

We study the global dust and (molecular) gas content in the Lyman Alpha Reference Sample (LARS), i.e. 14 local star-forming galaxies. We characterize their interstellar medium and relate newly derived properties to quantities relevant for Ly$\alpha$ escape. We observed LARS galaxies with Herschel/PACS, SOFIA/FIFI-LS, the IRAM 30m telescope and APEX, targeting far-infrared (FIR) continuum and emission lines of [C II]158$\mu$m, [O I]63$\mu$m, [O III]88$\mu$m and low-J CO lines. Using Bayesian methods we derive dust model parameters and estimate total gas masses for all LARS galaxies, taking into account a metallicity-dependent gas-to-dust ratio. Star formation rates were estimated from FIR, [C II]158$\mu$m and [O I]63$\mu$m luminosities. LARS covers a wide dynamic range in the derived properties, with FIR-based star formation rates from $\sim$0.5-100 $M_{\odot}\ yr^{-1}$, gas fractions between $\sim$15-80% and gas depletion times ranging from a few hundred Myr up to more than 10 Gyr. The distribution of LARS galaxies in the $\Sigma_{gas}$ vs. $\Sigma_{SFR}$ (Kennicutt-Schmidt plane) is thus quite heterogeneous. However, we find that LARS galaxies with the longest gas depletion times, i.e. relatively high gas surface densities ($\Sigma_{gas}$) and low star formation rate densities ($\Sigma_{SFR}$), have by far the highest Ly$\alpha$ escape fraction. A strong $\sim$linear relation is found between Ly$\alpha$ escape fraction and the total gas (HI+H$_2$) depletion time. We argue that the Ly$\alpha$ escape in those galaxies is driven by turbulence in the star-forming gas that shifts the Ly$\alpha$ photons out of resonance close to the places where they originate. We further report on an extreme [C II]158$\mu$m excess in LARS 5, corresponding to $\sim$14$\pm$3% of the FIR luminosity, i.e. the most extreme [C II]-to-FIR ratio observed in a non-AGN galaxy to date.

Stellar populations in barred galaxies save an imprint of the influence of the bar on the host galaxy's evolution. We present a detailed analysis of star formation histories (SFHs) and chemical enrichment of stellar populations in nine nearby barred galaxies from the TIMER project. We use integral field observations with the MUSE instrument to derive unprecedented spatially resolved maps of stellar ages, metallicities, [Mg/Fe] abundances and SFHs, as well as H$\alpha$ as a tracer of ongoing star formation. We find a characteristic V-shaped signature in the SFH perpendicular to the bar major axis which supports the scenario where intermediate age stars ($\sim 2$-$6\ \mathrm{Gyr}$) are trapped on more elongated orbits shaping a thinner part of the bar, while older stars ($> 8\ \mathrm{Gyr}$) are trapped on less elongated orbits shaping a rounder and thicker part of the bar. We compare our data to state-of-the-art cosmological magneto-hydrodynamical simulations of barred galaxies and show that such V-shaped SFHs arise naturally due to the dynamical influence of the bar on stellar populations with different ages and kinematic properties. Additionally, we find an excess of very young stars ($< 2\ \mathrm{Gyr}$) on the edges of the bars, predominantly on the leading side, confirming typical star formation patterns in bars. Furthermore, mass-weighted age and metallicity gradients are slightly shallower along the bar than in the disc likely due to orbital mixing in the bar. Finally, we find that bars are mostly more metal-rich and less [Mg/Fe]-enhanced than the surrounding discs. We interpret this as a signature that the bar quenches star formation in the inner region of discs, usually referred to as star formation deserts. We discuss these results and their implications on two different scenarios of bar formation and evolution.

The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field.

We present an empirical relation between the cold gas surface density ($\Sigma_{\rm gas}$) and the optical extinction (${\rm A_V}$) in a sample of 103 galaxies from the Extragalactic Database for Galaxy Evolution (EDGE) survey. This survey provides CARMA interferometric CO observations for 126 galaxies included in the Calar Alto Legacy Integral Field Area (CALIFA) survey. The matched, spatially resolved nature of these data sets allows us to derive the $\Sigma_{\rm gas}$-${\rm A_V}$ relation on global, radial, and kpc (spaxel) scales. We determine ${\rm A_V}$ from the Balmer decrement (H$\alpha$/H$\beta$). We find that the best fit for this relation is $\Sigma_{\rm gas} ({\rm M_\odot pc^{-2}})\sim~26~\times~ {\rm A_V}({\rm mag})$, and that it does not depend on the spatial scale used for the fit. However, the scatter in the fits increases as we probe smaller spatial scales, reflecting the complex relative spatial distributions of stars, gas, and dust. We investigate the $\Sigma_{\rm gas}$/ ${\rm A_V}$ ratio on radial and spaxel scales as a function of ${\rm EW(H\alpha)}$. We find that at larger values of ${\rm EW(H\alpha)}$ (i.e., actively star-forming regions) this ratio tend to converge to the value expected for dust-star mixed geometries ($\sim$ 30 $\mathrm{M_{\odot} \,pc^{-2}\,mag^{-1}}$). On radial scales, we do not find a significant relation between the $\Sigma_{\rm gas}$/${\rm A_V}$ ratio and the ionized gas metallicity. We contrast our estimates of $\Sigma_{\rm gas}$ using ${\rm A_V}$ with compilations in the literature of the gas fraction on global and radial scales as well as with well known scaling relations such as the radial star-formation law and the $\Sigma_{\rm gas}$-$\Sigma_*$ relation. These tests show that optical extinction is a reliable proxy for estimating $\Sigma_{\rm gas}$ in the absence of direct sub/millimeter observations of the cold gas.

[Abridged] We present the apparent stellar angular momentum of 300 galaxies across the Hubble sequence, using integral-field spectroscopic data from the CALIFA survey. Adopting the same $\lambda_\mathrm{R}$ parameter previously used to distinguish between slow and fast rotating early-type (elliptical and lenticular) galaxies, we show that spiral galaxies as expected are almost all fast rotators. Given the extent of our data, we provide relations for $\lambda_\mathrm{R}$ measured in different apertures, including conversions to long-slit 1D apertures. Our sample displays a wide range of $\lambda_\mathrm{Re}$ values, consistent with previous IFS studies. The fastest rotators are dominated by relatively massive and highly star-forming Sb galaxies, which preferentially reside in the main star-forming sequence. These galaxies reach $\lambda_\mathrm{Re}$ values of $\sim$0.85, are the largest galaxies at a given mass, and display some of the strongest stellar population gradients. Compared to the population of S0 galaxies, our findings suggest that fading may not be the dominant mechanism transforming spirals into lenticulars. Interestingly, we find that $\lambda_\mathrm{Re}$ decreases for late-type Sc and Sd spiral galaxies, with values than in occasions puts them in the slow-rotator regime. While for some of them this can be explained by their irregular morphologies and/or face-on configurations, others are edge-on systems with no signs of significant dust obscuration. The latter are typically at the low-mass end, but this does not explain their location in the classical ($V/\sigma$,$\varepsilon$) and ($\lambda_\mathrm{Re}$,$\varepsilon$) diagrams. Our initial investigations, based on dynamical models, suggest that these are dynamically hot disks, probably influenced by the observed important fraction of dark matter within R$_\mathrm{e}$.

Stellar feedback plays a significant role in modulating star formation, redistributing metals, and shaping the baryonic and dark structure of galaxies -- however, the efficiency of its energy deposition to the interstellar medium is challenging to constrain observationally. Here we leverage HST and ALMA imaging of a molecular gas and dust shell ($M_{H2} \sim 2\times 10^{5} ~{\rm M}_{\odot}$) in an outflow from the nuclear star forming ring of the galaxy NGC 3351, to serve as a boundary condition for a dynamical and energetic analysis of the outflowing ionised gas seen in our MUSE TIMER survey. We use \texttt{STARBURST99} models and prescriptions for feedback from simulations to demonstrate that the observed star formation energetics can reproduce the ionised and molecular gas dynamics -- provided a dominant component of the momentum injection comes from direct photon pressure from young stars, on top of supernovae, photoionisation heating and stellar winds. The mechanical energy budget from these sources is comparable to low luminosity AGN, suggesting that stellar feedback can be a relevant driver of bulk gas motions in galaxy centres - although here $\lesssim 10^{-3}$ of the ionized gas mass is escaping the galaxy. We test several scenarios for the survival/formation of the cold gas in the outflow, including in-situ condensation and cooling. Interestingly, the geometry of the molecular gas shell, observed magnetic field strengths and emission line diagnostics are consistent with a scenario where magnetic field lines aided survival of the dusty ISM as it was initially launched (with mass loading factor $\lesssim 1$) from the ring by stellar feedback. This system's unique feedback driven morphology can hopefully serve as a useful litmus test for feedback prescriptions in magnetohydrodynamical galaxy simulations.

We present a convenient, all-in-one framework for the scientific analysis of fully reduced, (integral-field) spectroscopic data. The GIST pipeline (Galaxy IFU Spectroscopy Tool) is entirely written in Python3 and conducts all steps from the preparation of input data, over the scientific analysis to the production of publication-quality plots. In its basic setup, it extracts stellar kinematics, performs an emission-line analysis and derives stellar population properties from full spectral fitting as well as via the measurement of absorption line-strength indices by exploiting the well-known pPXF and GandALF routines, where the latter has now been implemented in Python. The pipeline is not specific to any instrument or analysis technique and provides easy means of modification and further development, as of its modular code architecture. An elaborate, Python-native parallelisation is implemented and tested on various machines. The software further features a dedicated visualization routine with a sophisticated graphical user interface. This allows an easy, fully-interactive plotting of all measurements, spectra, fits, and residuals, as well as star formation histories and the weight distribution of the models. The pipeline has successfully been applied to both low and high-redshift data from MUSE, PPAK (CALIFA), and SINFONI, as well as to simulated data for HARMONI@ELT and WEAVE and is currently being used by the TIMER, Fornax3D, and PHANGS collaborations. We demonstrate its capabilities by applying it to MUSE TIMER observations of NGC 1433.

The [CII]158$\mu$m line is one of the strongest far-infrared (FIR) lines and an important coolant in the interstellar medium of galaxies that is accessible out to high redshifts. The excitation of [CII] is complex and can best be studied in detail at low redshifts. Here we report the discovery of the highest global [CII] excess with respect to the FIR luminosity in the nearby AGN host galaxy HE1353-1917. This galaxy is exceptional among a sample of five targets because the AGN ionization cone and radio jet directly intercept the cold galactic disk. As a consequence, a massive multiphase gas outflow on kiloparsec scales is embedded in an extended narrow-line region. Because HE1353-1917 is distinguished by these special properties from our four bright AGN, we propose that a global [CII] excess in AGN host galaxies could be a direct signature of a multiphase AGN-driven outflow with a high mass-loading factor.

[Abridged] We combine extensive spatially-resolved multi-wavelength observations, taken as part of the Close AGN Reference Survey (CARS), for the edge-on disc galaxy HE1353-1917 to characterize the impact of the AGN on its host galaxy via outflows and radiation. Multi-color broad-band photometry is combined with spatially-resolved optical, NIR and sub-mm and radio observations taken with VLT/MUSE, Gemini-N/NIFS, ALMA and the VLA to map the physical properties and kinematics of the multi-phase inter-stellar medium (ISM). We detect a biconical extended narrow-line region (ENLR) ionized by the luminous AGN oriented nearly parallel to the galaxy disc, extending out to at least 25kpc. The extra-planar gas originates from galactic fountains initiated by star formation processes in the disc, rather than an AGN outflow, as shown by the kinematics and the metallicity of the gas. Nevertheless, a fast multi-phase AGN-driven outflow with speeds up to 1000km/s is detected close to the nucleus at 1kpc distance. A radio jet, in connection with the AGN radiation field, is likely responsible for driving the outflow as confirmed by the energetics and the spatial alignment of the jet and multi-phase outflow. Evidence for negative AGN feedback suppressing the star formation rate (SFR) is mild and restricted to the central kpc. But while any SFR suppression must have happened recently, the outflow has the potential to greatly impact the future evolution of the galaxy disc due to its geometrical orientation. Our observations reveal that low-power radio jets can play a major role in driving fast multi-phase galaxy-scale outflows even in radio-quiet AGN. Since the outflow energetics for HE1353-1917 are consistent with literature scaling relations of AGN-driven outflows the contribution of radio jets as the driving mechanisms still needs to be systematically explored.

We present an unprecedented view on the morphology and kinematics of the extended narrow-line region (ENLR) and molecular gas around the prototypical hyper-luminous quasar 3C273 ($L\sim10^{47}$ erg/s at z=0.158) based on VLT-MUSE optical 3D spectroscopy and ALMA observations. We find that: 1) The ENLR size of 12.1$\pm$0.2kpc implies a smooth continuation of the size-luminosity relation out to large radii or a much larger break radius as previously proposed. 2) The kinematically disturbed ionized gas with line splits reaching 1000km/s out to 6.1$\pm$1.5kpc is aligned along the jet axis. 3) The extreme line broadening on kpc scales is caused by spatial and spectral blending of many distinct gas clouds separated on sub-arcsecond scales with different line-of-sight velocities. The ENLR velocity field combined with the known jet orientation rule out a simple scenario of a radiatively-driven radial expansion of the outflow. Instead we propose that a pressurized expanding hot gas cocoon created by the radio jet is impacting on an inclined gas disk leading to transverse and/or backflow motion with respect to our line-of-sight. The molecular gas morphology may either be explained by a density wave at the front of the outflow expanding along the jet direction as predicted by positive feedback scenario or the cold gas may be trapped in a stellar over-density caused by a recent merger event. Using 3C273 as a template for observations of high-redshift hyper-luminous AGN reveals that large-scale ENLRs and kpc scale outflows may often be missed due to the brightness of the nuclei and the limited sensitivity of current near-IR instrumentation.

The absence of star formation in the bar region that has been reported for some galaxies can theoretically be explained by shear. However, it is not clear how star-forming (SF) bars fit into this picture and how the dynamical state of the bar is related to other properties of the host galaxy. We used integral-field spectroscopy from VLT/MUSE to investigate how star formation within bars is connected to structural properties of the bar and the host galaxy. We derived spatially resolved H$\alpha$ fluxes from MUSE observations from the CARS survey to estimate star formation rates in the bars of 16 nearby ($0.01 < z < 0.06$) disc galaxies with stellar masses between $10^{10} M_\odot$ and $10^{11} M_\odot$. We further performed a detailed multicomponent photometric decomposition on images derived from the data cubes. We find that bars clearly divide into SF and non-star-forming (non-SF) types, of which eight are SF and eight are non-SF. Whatever the responsible quenching mechanism is, it is a quick process compared to the lifetime of the bar. The star formation of the bar appears to be linked to the flatness of the surface brightness profile in the sense that only the flattest bars $\left(n_\mathrm{bar} \leq 0.4\right)$ are actively SF $\left(\mathrm{SFR_{b}} > 0.5 M_\odot \mathrm{yr^{-1}}\right)$. Both parameters are uncorrelated with Hubble type. We find that star formation is 1.75 times stronger on the leading than on the trailing edge and is radially decreasing. The conditions to host non-SF bars might be connected to the presence of inner rings. Additionally, from testing an AGN feeding scenario, we report that the star formation rate of the bar is uncorrelated with AGN bolometric luminosity. The results of this study may only apply to type-1 AGN hosts and need to be confirmed for the full population of barred galaxies.

Research over the past decade has shown diminishing empirical evidence for major galaxy mergers being a dominating or even important mechanism for the growth of supermassive black holes in galaxies and the triggering of optically or X-ray selected active galactic nuclei (AGN). We here for the first time test whether such a connection exists at least in the most plausible part of parameter space for this mechanism: the highest specific accretion rate broad-line AGNs at the peak epoch of black hole activity around z = 2. To that end we examine 21 galaxies hosting a high accreting black hole (L/Ledd > 0.7) observed with HST/WFC3 and 92 stellar mass- and redshift- matched inactive galaxies taken from the CANDELS survey. We removed the AGN point sources from their host galaxies and avoided bias in visual classification by adding and then subtracting mock point sources to and from the comparison galaxies, producing matched residual structures for both sets. The resulting samples were joined, randomized, and subsequently visually ranked with respect to perceived strength of structural distortions by 10 experts. The ensuing individual rankings were combined into a consensus sequence and from this we derived merger fractions for both samples. With the merger fractions f$_{m,agn}$ = 0.24 $\pm$ 0.09 for the AGN host galaxy sample and f$_{m,ina}$ = 0.19 $\pm$ 0.04 for the inactive galaxies, we find no significant difference between the AGN host galaxies and inactive galaxies. Also, both samples display comparable fractions of disk-dominated galaxies. These findings are consistent with previous studies for different AGN populations, and we conclude that even black hole growth at the highest specific accretion rates and at the peak of cosmic AGN activity is not predominantly caused by major mergers. (abriged)

The formation of two stellar bars within a galaxy has proved challenging for numerical studies. It is yet not clear whether the inner bar is born via a star formation process promoted by gas inflow along the outer bar, or whether it is dynamically assembled from instabilities in a small-scale stellar disc. Observational constraints to these scenarios are scarce. We present a thorough study of the stellar content of two double-barred galaxies observed by the MUSE TIMER project, NGC 1291 and NGC 5850, combined with a two-dimensional multi-component photometric decomposition performed on the 3.6{\mu}m images from S4G. Our analysis confirms the presence of {\sigma}-hollows appearing in the stellar velocity dispersion distribution at the ends of the inner bars. Both galaxies host inner discs matching in size with the inner bars, suggestive of a dynamical formation for the inner bars from small-scale discs. The analysis of the star formation histories for the structural components shaping the galaxies provides constraints on the epoch of dynamical assembly of the inner bars, which took place >6.5 Gyr ago for NGC 1291 and >4.5 Gyr ago for NGC 5850. This implies that inner bars are long-lived structures.

Aims: We probe the radiatively-efficient, hot wind feedback mode in two nearby luminous unobscured (type 1) AGN from the Close AGN Reference Survey (CARS), which show intriguing kpc-scale arc-like features of extended [OIII] ionized gas as mapped with VLT-MUSE. We aimed to detect hot gas bubbles that would indicate the existence of powerful, galaxy-scale outflows in our targets, HE 0227-0931 and HE 0351+0240, from deep (200 ks) Chandra observations. Methods: By measuring the spatial and spectral properties of the extended X-ray emission and comparing with the sub kpc-scale IFU data, we are able to constrain feedback scenarios and directly test if the ionized gas is due to a shocked wind. Results: No extended hot gas emission on kpc-scales was detected. Unless the ambient medium density is low ($n_{H}\sim~1$ cm$^{-3}$ at 100 pc), the inferred upper limits on the extended X-ray luminosities are well below what is expected from theoretical models at matching AGN luminosities. Conclusions: We conclude that the highly-ionized gas structures on kpc scales are not inflated by a hot outflow in either target, and instead are likely caused by photo-ionization of pre-existing gas streams of different origins. Our non-detections suggest that extended X-ray emission from an AGN-driven wind is not universal, and may lead to conflicts with current theoretical predictions.

The Time Inference with MUSE in Extragalactic Rings (TIMER) project is a survey with the VLT-MUSE integral-field spectrograph of 24 nearby barred galaxies with prominent central structures (e.g., nuclear rings or inner discs). The main goals of the project are: (i) estimating the cosmic epoch when discs of galaxies settle, leading to the formation of bars; (ii) testing the hypothesis whereby discs in more massive galaxies are assembled first; and (iii) characterising the history of external gas accretion in disc galaxies. We present details on the sample selection, observations, data reduction, and derivation of high-level data products, including stellar kinematics, ages and metallicities. We also derive star formation histories and physical properties and kinematics of ionised gas. We illustrate how this dataset can be used for a plethora of scientific applications, e.g., stellar feedback, outflows, nuclear and primary bars, stellar migration and chemical enrichment, and the gaseous and stellar dynamics of nuclear spiral arms, barlenses, box/peanuts and bulges. Amongst our first results - based on a few selected galaxies -, we show that the dynamics of nuclear rings and inner discs is consistent with the picture in which they are formed by bars, that the central few hundred parsecs in massive disc galaxies tend to show a pronounced peak in stellar metallicity, and that nuclear rings can efficiently prevent star formation in this region. Finally, we present evidence that star-bursting nuclear rings can be fed with low-metallicity gas from low-mass companions.

We report spatially-resolved [CII]$\lambda 158$ $\mu$m observations of HE 0433-1028, which is the first detection of a nearby luminous AGN (redshift 0.0355) with FIFI-LS onboard the airborne observatory SOFIA. We compare the spatially-resolved star formation tracers [CII], as provided by our SOFIA observations, and H$\alpha$ from MUSE optical integral-field spectroscopy. We find that the [CII] emission is mainly matching the extended star formation as traced by the extinction-corrected H$\alpha$ line emission but some additional flux is present. While a larger sample is needed to statistically confirm our findings and investigate possible dependencies on AGN luminosity and star formation rate, our study underlines the necessity of collecting a spatially-resolved optical-FIR dataset for nearby AGNs, and shows that it is technically feasible to collect such datasets with FIFI-LS onboard SOFIA.

Theoretical models of galaxy formation suggest that the presence of an active galactic nucleus (AGN) is required to regulate the growth of its host galaxy through feedback mechanisms, produced by e.g. AGN-driven outflows. Although such outflows are common both at low and high redshift, a comprehensive picture is still missing. The peak epoch of galaxy assembly (1<z<3) has been poorly explored so far, and current observations in this redshift range are mostly limited to targets with high chances to be in an outflowing phase. This paper introduces SUPER (a SINFONI Survey for Unveiling the Physics and Effect of Radiative feedback), an ongoing ESO's VLT/SINFONI Large Programme. SUPER will perform the first systematic investigation of ionized outflows in a sizeable and blindly-selected sample of 39 X-ray AGN at z~2, which reaches high spatial resolutions (~2 kpc) thanks to the adaptive optics-assisted IFU observations. The outflow morphology and star formation in the host galaxy will be mapped through the broad component of [OIII] and the narrow component of Ha emission lines. The main aim of our survey is to infer the impact of outflows on the on-going star formation and to link the outflow properties to a number of AGN and host galaxy properties. We describe here the survey characteristics and goals, as well as the selection of the target sample. Moreover, we present a full characterization of its multi-wavelength properties: we measure, via spectral energy distribution fitting of UV-to-FIR photometry, stellar masses (4x10^9-2x10^11 Msun), star formation rates (25-680 Msun yr^-1) and AGN bolometric luminosities (2x10^44-8x10^47 erg s^-1), along with obscuring column densities (<2x10^24 cm^-2) and 2-10 keV luminosities (2x10^43-6x10^45 erg s^-1) derived through X-ray spectral analysis. Finally, we classify our AGN as jetted or non-jetted according to their radio and FIR emission.

The James Webb Space Telescope (JWST) will enable the detection of optical emission lines in galaxies spanning a broad range of luminosities out to redshifts z>10. Measurements of key galaxy properties, such as star formation rate and metallicity, through these observations will provide unique insight into, e.g., the role of feedback from stars and active galactic nuclei (AGNs) in regulating galaxy evolution, the co-evolution of AGNs and host galaxies, the physical origin of the 'main sequence' of star-forming galaxies and the contribution by star-forming galaxies to cosmic reionization. We present an original framework to simulate and analyse observations performed with the Near Infrared Spectrograph (NIRSpec) on board JWST. We use the BEAGLE tool (BayEsian Analysis of GaLaxy sEds) to build a semi-empirical catalogue of galaxy spectra based on photometric spectral energy distributions (SEDs) of dropout galaxies in the Hubble Ultra Deep Field (HUDF). We demonstrate that the resulting catalogue of galaxy spectra satisfies different types of observational constraints on high redshift galaxies, and use it as input to simulate NIRSpec/prism (R~100) observations. We show that a single 'deep' (~100 ks) NIRSpec/prism pointing in the HUDF will enable S/N>3 detections of multiple optical emission lines in ~30 (~60) galaxies at z>6 (z~4-6) down to m_F160W<30 AB mag. Such observations will allow measurements of galaxy star formation rates, ionization parameters and gas-phase metallicities within factors of 1.5, mass-to-light ratios within a factor of 2, galaxy ages within a factor of 3 and V-band attenuation optical depths with a precision of 0.3.

We present ALMA and MUSE observations of the Brightest Cluster Galaxy in Abell 2597, a nearby (z=0.0821) cool core cluster of galaxies. The data map the kinematics of a three billion solar mass filamentary nebula that spans the innermost 30 kpc of the galaxy's core. Its warm ionized and cold molecular components are both cospatial and comoving, consistent with the hypothesis that the optical nebula traces the warm envelopes of many cold molecular clouds that drift in the velocity field of the hot X-ray atmosphere. The clouds are not in dynamical equilibrium, and instead show evidence for inflow toward the central supermassive black hole, outflow along the jets it launches, and uplift by the buoyant hot bubbles those jets inflate. The entire scenario is therefore consistent with a galaxy-spanning "fountain", wherein cold gas clouds drain into the black hole accretion reservoir, powering jets and bubbles that uplift a cooling plume of low-entropy multiphase gas, which may stimulate additional cooling and accretion as part of a self-regulating feedback loop. All velocities are below the escape speed from the galaxy, and so these clouds should rain back toward the galaxy center from which they came, keeping the fountain long-lived. The data are consistent with major predictions of chaotic cold accretion, precipitation, and stimulated feedback models, and may trace processes fundamental to galaxy evolution at effectively all mass scales.

NGC7252, which is one of the nearest major-merger galaxy remnants, is an ideal laboratory to study the processes inherent to the transformation of disc galaxies to ellipticals as observed about ~1Gyr after the collision. We obtained wide-field IFU spectroscopy with the VLT-VIMOS integral-field spectrograph covering the central 50"x50" of NGC7252 to map the stellar and ionised gas kinematics, and the distribution and conditions of the ionised gas, revealing the extent of ongoing star formation and recent star formation history. We find that the inner gas disc is not counter-rotating with respect to the stars and that the stellar kinematics appear complex with a clear indication of a prolate-like rotation component suggesting a polar merger configuration. The ongoing star formation rate is 2.2+-0.6 M_sun/yr and implies a typical depletion time of 2Gyr given the molecular gas content. Furthermore, the spatially-resolved star formation history suggests a slight radial dependence, moving outwards at later times. We confirm a large AGN-ionised gas cloud previously discovered 5kpc south of the nucleus, and find higher ionisation also at the galaxy centre relative to the surrounding gas disc. Although the higher ionisation towards the centre is potentially degenerate within the central star forming ring, it may be associated with a low-luminosity AGN. Although NGC7252 has been classified as post-starburst galaxy at the centre, the elliptical-like major-merger remnant still appears very active. A central kpc-scale gas disc has presumably re-formed quickly within the last 100Myr after final coalescence. The disc features ongoing star formation, implying Gyr long timescale to reach the red sequence through gas consumption alone. While NGC7252 is useful to probe the transformation from discs to ellipticals, it is not well-suited to study the transformation from blue to red at this point.

We report on an optical (SDSS) and X-ray (XMM) study of an optically selected sample of four dual AGN systems at projected separations of 30--60~kpc. All sources are detected in the X-ray band (0.3-10 keV); seven objects are optically identified as Seyfert, while one source, optically classified as a LINER, is likely powered by accretion in virtue of its relatively high X-ray luminosity (1.2$\times10^{41}$ erg/s). Six of the eight objects are obscured in X-rays with N$_{\rm H} \geq$ 10$^{23}$ cm$^{-2}$; three of these, whose X-ray spectrum is dominated by a reflection component, are likely Compton-thick (N$_{\rm H} \geq$ 10$^{24}$ cm$^{-2}$). This finding is in agreement with the hypothesis that galaxy encounters are effective in driving gas inflow toward the nuclear region, thus increasing the obscuration. We compare the absorption properties in our dual AGN with those in larger samples observed in X-rays but selected in different ways (optical, IR and hard X-rays). We find that the obscured (N$_{\rm H} \geq$ 10$^{22}$ cm$^{-2}$) AGN fraction within the larger sample is 84$\pm$4 per cent (taking into account the 90 per cent error on the N$_{\rm H}$ measure) up to large pair separations ($\sim$100~kpc). This is statistically higher than the fraction of obscured AGN in isolated galaxies found in X-ray surveys.

We recently discovered that the luminous radio-quiet QSO LBQS0302-0019 at z=3.286 is likely accompanied by an obscured AGN at 20kpc projected distance, which we dubbed Jil. It represents the tightest candidate obscured/unobscured dual AGN system at z>3. To verify the dual AGN scenario we obtained deep $K_s$ band (rest-frame $V$ band) imaging with the VLT/HAWK-I+GRAAL instrument at 0.4" resolution during science verification in January 2018. Indeed, we detect the individual host galaxies of the QSO and Jil with estimated stellar masses of $\log(M_\star/M_{\odot})=11.4\pm0.5$ and $\log(M_\star/M_{\odot})=10.9\pm0.5$, respectively. Near-IR spectra obtained with VLT-KMOS reveal a clear [OIII] $\lambda$5007 line detection at the location of Jil which does not contribute significantly to the $K_s$ band flux. Both observations therefore corroborate the dual AGN scenario. A comparison to Illustris cosmological simulations suggests a parent halo mass of $\log(M_\mathrm{halo}/M_{\odot})=13.2\pm0.5$ for this interacting galaxy system, corresponding to a very massive dark matter halo at that epoch.

We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved $^{12}$CO(J=1-0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and H$\alpha$ rotation curves using the same geometric parameters out to $\gtrsim$1 $R_e$. We find that $\sim$75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence supported disk. Using observations of the H$\gamma$ transition also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and H$\alpha$ rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.

We present an original phenomenological model to describe the evolution of galaxy number counts, morphologies, and spectral energy distributions across a wide range of redshifts (0.2<z<15) and stellar masses [Log10 M/Msun >6]. Our model follows observed mass and luminosity functions of both star-forming and quiescent galaxies, and reproduces the redshift evolution of colors, sizes, star-formation and chemical properties of the observed galaxy population. Unlike other existing approaches, our model includes a self-consistent treatment of stellar and photoionized gas emission and dust attenuation based on the BEAGLE tool. The mock galaxy catalogs generated with our new model can be used to simulate and optimize extragalactic surveys with future facilities such as the James Webb Space Telescope (JWST), and to enable critical assessments of analysis procedures, interpretation tools, and measurement systematics for both photometric and spectroscopic data. As a first application of this work, we make predictions for the upcoming JWST Advanced Deep Extragalactic Survey (JADES), a joint program of the JWST/NIRCam and NIRSpec Guaranteed Time Observations teams. We show that JADES will detect, with NIRCam imaging, thousands of galaxies at z>6, and tens at z>10 at m_AB<30 (5-sigma) within the 236 arcmin^2 of the survey. The JADES data will enable accurate constraints on the evolution of the UV luminosity function at z>8, and resolve the current debate about the rate of evolution of galaxies at z>8. Ready to use mock catalogs and software to generate new realizations are publicly available as the JAdes extraGalactic Ultradeep Artificial Realizations (JAGUAR) package.

A prediction of the current paradigm of the hierarchical assembly of galaxies is the presence of supermassive dual black holes at separations of a few kpc or less. In this context, we report the detection of a narrow-line emitter within the extended LyA nebula (~120kpc diameter) of the luminous radio-quiet quasi-stellar object (QSO) LBQS0302-0019 at z=3.286. We identify several high-ionization narrow emission lines (HeII, CIV, CIII]) associated with this point-like source, which we have named "Jil", which is only ~20kpc (2.9") away from the QSO in projection. Emission-line diagnostics confirm that the source is likely powered by photoionization of an obscured active galactic nucleus (AGN) three orders of magnitude fainter than the QSO. The system represents the tightest unobscured/obscured dual AGN currently known at z>3, highlighting the power of MUSE to detect these elusive systems.

Observations and simulations show that outflows in active galactic nuclei (AGN) contain gas in different phases. To understand their true impact on galaxy evolution, we advocate consistent and unbiased investigation of these multiphase winds in large AGN samples.

Feedback from active galactic nuclei (AGN) remains controversial despite its wide acceptance as necessary to regulate massive galaxy growth. A dedicated workshop was held on 16-20 October 2017 at the Lorentz Center in Leiden to distinguish between the reality and myths of AGN feedback from the observational side. Here, we summarize briefly all the sessions and outcome of the stimulating workshop. More details on the outcome of the discussions are provided in a series of articles.

We present a kpc-scale analysis of the relationship between the molecular depletion time ($\tau_\mathrm{dep}^\mathrm{mol}$) and the orbital time ($\tau_\mathrm{orb}$) across the field of 39 face-on local galaxies, selected from the EDGE-CALIFA sample. We find that, on average, 5% of the available molecular gas is converted into stars per orbital time, or $\tau_\mathrm{dep}^\mathrm{mol}\sim20\tau_\mathrm{orb}$. The resolved relation shows a scatter of $\sim0.5$ dex. The scatter is ascribable to galaxies of different morphologies that follow different $\tau_\mathrm{dep}^\mathrm{mol}-\tau_\mathrm{orb}$ relations which decrease in steepness from early- to late-types. The morphologies appear to be linked with the star formation rate surface density, the molecular depletion time, and the orbital time, but they do not correlate with the molecular gas content of the galaxies in our sample. We speculate that in our molecular gas rich, early-type galaxies, the morphological quenching (in particular the disc stabilization via shear), rather than the absence of molecular gas, is the main factor responsible for their current inefficient star formation.

After changing optical AGN type from 1.9 to 1 in 1984, the AGN Mrk 1018 recently reverted back to its type 1.9 state. Our ongoing monitoring now reveals that the AGN has halted its dramatic dimming, reaching a minimum around October 2016. The minimum was followed by an outburst rising with $\sim$0.25 U-band mag/month. The rebrightening lasted at least until February 2017, as confirmed by joint Chandra and Hubble observations. Monitoring was resumed in July 2017 after the source emerged from sunblock, at which point the AGN was found only $\sim$0.4 mag brighter than its minimum. The intermittent outburst was accompanied by the appearance of a red wing asymmetry in broad-line shape, indicative of an inhomogeneous broad-line region. The current flickering brightness of Mrk 1018 following its rapid fading either suggests that the source has reignited, remains variable at a low level, or may continue dimming over the next few years. Distinguishing between these possibilities requires continuous multiwavelength monitoring.

We present MUSE observations of the field of the quasar Q0152$-$020 whose spectrum shows a Lyman limit system (LLS) at redshift $z_{\rm abs} = 0.38$, with a metallicity Z $\gtrsim 0.06$ Z$_\odot$. The low ionization metal lines associated with the LLS present two narrow distinct absorption components with a velocity separation of 26 km ${\rm s}^{-1}$. We detect six galaxies within 600 km ${\rm s}^{-1}$ from the absorption redshift; their projected distances from the quasar sightline range from 60 to 200 kpc. The optical spectra of five of these galaxies exhibit prominent nebular emission lines, from which we deduce extinction-corrected star formation rates in the range SFR = 0.06-1.3 M$_\odot$~yr$^{-1}$, and metallicities between 0.2 Z$_\odot$ and Z$_\odot$. The sixth galaxy is only detected in the stellar continuum. By combining our data with archival Keck/HIRES spectroscopy of the quasar and HST/WFPC2 imaging of the field, we can relate absorption line and galaxy kinematics; we conclude that the LLS is most likely associated with the galaxy closest to the quasar sight-line (galaxy "a"). Our morphokinematic analysis of galaxy "a" combined with the absorption line kinematics supports the interpretation that one of the absorption components originates from an extension of the stellar disk of galaxy "a", while the other component may arise in accreting gas in a warped disk with specific angular momentum $\sim 3$ times larger than the specific angular momentum of the galaxy halo. Such warped disks are common features in hydrodynamical simulations of cold-flow accretion onto galaxies; the data presented here provide observational evidence in favour of this scenario.

We present results from the EDGE survey, a spatially resolved CO(1-0) follow-up to CALIFA, an optical Integral Field Unit (IFU) survey of local galaxies. By combining the data products of EDGE and CALIFA, we study the variation in molecular gas depletion time ($\tau_{\rm dep}$) on kiloparsec scales in 52 galaxies. We divide each galaxy into two parts: the center, defined as the region within $0.1 \ R_{25}$, and the disk, defined as the region between $0.1$ and $0.7 \ R_{25}$. We find that 14 galaxies show a shorter $\tau_{\rm dep}$ ($\sim 1$ Gyr) in the center relative to that in the disk ($\tau_{\rm dep} \sim 2.4$ Gyrs), which means the central region in those galaxies is more efficient at forming stars per unit molecular gas mass. This finding implies that the centers with shorter $\tau_{\rm dep}$ resemble the intermediate regime between galactic disks and starburst galaxies. Furthermore, the central drop in $\tau_{\rm dep}$ is correlated with a central increase in the stellar surface density, suggesting that a shorter $\tau_{\rm dep}$ is associated with molecular gas compression by the stellar gravitational potential. We argue that varying the CO-to-H$_2$ conversion factor only exaggerates the central drop of $\tau_{\rm dep}$.

We present optical integral field spectroscopy for five $z<0.062$ narrow-line Seyfert 1 galaxies (NLS1s) host galaxies, probing their host galaxies at $\gtrsim 2-3$ kpc scales. Emission lines in the nuclear AGN spectra and the large-scale host galaxy are analyzed separately, based on an AGN-host decomposition technique. The host galaxy gas kinematics indicates large-scale gas rotation in all five sources. At the probed scales of $\gtrsim 2-3$ kpc, the host galaxy gas is found to be predominantly ionized by star formation without any evidence of a strong AGN contribution. None of the five objects shows specific star formation rates exceeding the main sequence of low-redshift star forming galaxies. The specific star formation rates for MCG-05-01-013 and WPVS 007 are roughly consistent with the main sequence, while ESO 399-IG20, MS 22549-3712, and TON S180 show lower specific star formation rates, intermediate to the main sequence and red quiescent galaxies. The host galaxy metallicities, derived for the two sources with sufficient data quality (ESO 399-IG20 and MCG-05-01-013), indicate central oxygen abundances just below the low-redshift mass-metallicity relation. Based on this initial case study, we outline a comparison of AGN and host galaxy parameters as a starting point for future extended NLS1 studies with similar methods.

We estimate the current extinction-corrected H$\alpha$ star formation rate (SFR) of the different morphological components that shape galaxies (bulges, bars, and disks). We use a multi-component photometric decomposition based on SDSS imaging to CALIFA Integral Field Spectroscopy datacubes for a sample of 219 galaxies. This analysis reveals an enhancement of the central SFR and specific SFR (sSFR $=$ SFR/$M_{\star}$) in barred galaxies. Along the Main Sequence, we find more massive galaxies in total have undergone efficient suppression (quenching) of their star formation, in agreement with many studies. We discover that more massive disks have had their star formation quenched as well. We evaluate which mechanisms might be responsible for this quenching process. The presence of type-2 AGNs plays a role at damping the sSFR in bulges and less efficiently in disks. Also, the decrease in the sSFR of the disk component becomes more noticeable for stellar masses around 10$^{10.5}$ M$_{\odot}$; for bulges, it is already present at $\sim$10$^{9.5}$ M$_{\odot}$. The analysis of the line-of-sight stellar velocity dispersions ($\sigma$) for the bulge component and of the corresponding Faber-Jackson relation shows that AGNs tend to have slightly higher $\sigma$ values than star-forming galaxies for the same mass. Finally, the impact of environment is evaluated by means of the projected galaxy density, $\Sigma$$_{5}$. We find that the SFR of both bulges and disks decreases in intermediate-to-high density environments. This work reflects the potential of combining IFS data with 2D multi-component decompositions to shed light on the processes that regulate the SFR.

We present interferometric CO observations made with the Combined Array for Millimeter-wave Astronomy (CARMA) of galaxies from the Extragalactic Database for Galaxy Evolution survey (EDGE). These galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) sample, mapped with optical integral field spectroscopy. EDGE provides good quality CO data (3$\sigma$ sensitivity $\Sigma_{\rm mol}$ $\sim$ 11 M$_\odot$ pc$^{-2}$ before inclination correction, resolution $\sim1.4$ kpc) for 126 galaxies, constituting the largest interferometric CO survey of galaxies in the nearby universe. We describe the survey, the data characteristics, the data products, and present initial science results. We find that the exponential scale-lengths of the molecular, stellar, and star-forming disks are approximately equal, and galaxies that are more compact in molecular gas than in stars tend to show signs of interaction. We characterize the molecular to stellar ratio as a function of Hubble type and stellar mass, present preliminary results on the resolved relations between the molecular gas, stars, and star formation rate, and discuss the dependence of the resolved molecular depletion time on stellar surface density, nebular extinction, and gas metallicity. EDGE provides a key dataset to address outstanding topics regarding gas and its role in star formation and galaxy evolution, which will be publicly available on completion of the quality assessment.

We present single-dish CO(1-0) and CO(2-1) observations for 14 low-redshift quasi-stellar objects (QSOs). In combination with optical integral field spectroscopy we study how the cold gas content relates to the star formation rate (SFR) and black hole accretion rate. CO(1-0) is detected in 8 of 14 targets and CO(2-1) is detected in 7 out of 11 cases. The majority of disc-dominated QSOs reveal gas fractions and depletion times well matching normal star forming systems. Two gas-rich major mergers show clear starburst signatures with higher than average gas fractions and shorter depletion times. Bulge-dominated QSO hosts are mainly undetected in CO(1-0) which corresponds, on average, to lower gas fractions than in disc-dominated counterparts. Their SFRs however imply shorter than average depletion times and higher star formation efficiencies. Negative QSO feedback through removal of cold gas seems to play a negligible role in our sample. We find a trend between black hole accretion rate and total molecular gas content for disc-dominated QSOs when combined with literature samples. We interpret this as an upper envelope for nuclear activity which is well represented by a scaling relation between the total and circum-nuclear gas reservoir accessible for accretion. Bulge-dominated QSOs significantly differ from that scaling relation and appear uncorrelated with the total molecular gas content. This could be explained either by a more compact gas reservoir, blow out of the gas envelope through outflows, or a different ISM phase composition.

Context: According to numerical simulations, stars are not always kept at their birth galactocentric distances but migrate. The importance of this radial migration in shaping galactic light distributions is still unclear. However, if it is indeed important, galaxies with different surface brightness (SB) profiles must display differences in their stellar population properties. Aims: We investigate the role of radial migration on the light distribution and the radial stellar content by comparing the inner colour, age and metallicity gradients for galaxies with different SB profiles. We define these inner parts avoiding the bulge and bar regions and up to around three disc scale-lengths (type I, pure exponential) or the break radius (type II, downbending; type III, upbending). Methods: We analyse 214 spiral galaxies from the CALIFA survey covering different SB profiles. We make use of GASP2D and SDSS data to characterise their light distribution and obtain colour profiles. The stellar age and metallicity profiles are computed using a methodology based on full-spectrum fitting techniques (pPXF, GANDALF, and STECKMAP) to the IFS CALIFA data. Results: The distributions of the colour, stellar age and stellar metallicity gradients in the inner parts for galaxies displaying different SB profiles are unalike as suggested by Kolmogorov-Smirnov and Anderson-Darling tests. We find a trend in which type II galaxies show the steepest profiles of all and type III the shallowest, with type I galaxies displaying an intermediate behaviour. Conclusions: These results are consistent with a scenario in which radial migration is more efficient for type III galaxies than for type I systems with type II galaxies presenting the lowest radial migration efficiency. In such scenario, radial migration mixes the stellar content flattening the radial stellar properties and... [abriged]

Understanding the nature of bulges in disc galaxies can provide important insights into the formation and evolution of galaxies. For instance, the presence of a classical bulge suggests a relatively violent history, in contrast, the presence of simply an inner disc (also referred to as a "pseudobulge") indicates the occurrence of secular evolution processes in the main disc. However, we still lack criteria to effectively categorise bulges, limiting our ability to study their impact on the evolution of the host galaxies. Here we present a recipe to separate inner discs from classical bulges by combining four different parameters from photometric and kinematic analyses: The bulge S\'ersic index $n_\mathrm{b}$, the concentration index $C_{20,50}$, the Kormendy (1977) relation and the inner slope of the radial velocity dispersion profile $\nabla\sigma$. With that recipe we provide a detailed bulge classification for a sample of 45 galaxies from the integral-field spectroscopic survey CALIFA. To aid in categorising bulges within these galaxies, we perform 2D image decomposition to determine bulge S\'ersic index, bulge-to-total light ratio, surface brightness and effective radius of the bulge and use growth curve analysis to derive a new concentration index, $C_{20,50}$. We further extract the stellar kinematics from CALIFA data cubes and analyse the radial velocity dispersion profile. The results of the different approaches are in good agreement and allow a safe classification for approximately $95\%$ of the galaxies. In particular, we show that our new "inner" concentration index performs considerably better than the traditionally used $C_{50,90}$ when yielding the nature of bulges. We also found that a combined use of this index and the Kormendy (1977) relation gives a very robust indication of the physical nature of the bulge.

Similarly to the cosmic star formation history, the black hole accretion rate density of the Universe peaked at 1<z<3. This cosmic epoch is hence best suited for investigating the effects of radiative feedback from AGN. Observational efforts are underway to quantify the impact of AGN feedback, if any, on their host galaxies. Here we present a study of the molecular gas content of AGN hosts at z~1.5 using CO[2-1] line emission observed with ALMA for a sample of 10 AGNs. We compare this with a sample of galaxies without an AGN matched in redshift, stellar mass, and star formation rate. We detect CO in 3 AGNs with $\mathrm{L_{CO} \sim 6.3-25.1\times 10^{9} L_{\odot}}$ which translates to a molecular hydrogen gas mass of $\mathrm{2.5-10\times 10^{10} M_{\odot}}$ assuming conventional conversion factor of $\mathrm{\alpha_{CO}}\sim3.6$. Our results indicate a >99% probability of lower depletion time scales and lower molecular gas fractions in AGN hosts with respect to the non-AGN comparison sample. We discuss the implications of these observations on the impact that AGN feedback may have on star formation efficiency of z>1 galaxies.

Dedicated photometric and spectroscopic surveys have provided unambiguous evidence for a strong stellar mass-size evolution of galaxies within the last 10 Gyr. The likely progenitors of today's most massive galaxies are remarkably small, disky, passive and have already assembled much of their stellar mass at redshift z=2. An in-depth analysis of these objects, however, is currently not feasible due to the lack of high-quality, spatially-resolved photometric and spectroscopic data. In this paper, we present a sample of nearby compact elliptical galaxies (CEGs), which bear resemblance to the massive and quiescent galaxy population at earlier times. Hubble Space Telescope (HST) and wide-field integral field unit (IFU) data have been obtained, and are used to constrain orbit-based dynamical models and stellar population synthesis (SPS) fits, to unravel their structural and dynamical properties. We first show that our galaxies are outliers in the present-day stellar mass-size relation. They are, however, consistent with the mass-size relation of compact, massive and quiescent galaxies at redshift z=2. The compact sizes of our nearby galaxies imply high central stellar mass surface densities, which are also in agreement with the massive galaxy population at higher redshift, hinting at strong dissipational processes during their formation. Corroborating evidence for a largely passive evolution within the last 10 Gyr is provided by their orbital distribution as well as their stellar populations, which are difficult to reconcile with a very active (major) merging history. This all supports that we can use nearby CEGs as local analogues of the high-redshift, massive and quiescent galaxy population, thus providing additional constraints for models of galaxy formation and evolution.

The light of the merger remnant NGC 1316 is dominated by old and intermediate-age stars. The only sign of current star formation in this big galaxy is the HII region SH2, an isolated star cluster complex with a ring-like morphology and an estimated age of 0.1 Gyr at a galactocentric distance of about 35 kpc. A nearby intermediate-age globular cluster, surrounded by weak line emission and a few more young star clusters, is kinematically associated. The origin of this complex is enigmatic. The nebular emission lines permit a metallicity determination which can discriminate between a dwarf galaxy or other possible precursors. We used the Integrated Field Unit of the VIMOS instrument at the Very Large Telescope of the European Southern Observatory to study the morphology, kinematics, and metallicity employing line maps, velocity maps, and line diagnostics of a few characteristic spectra. The line ratios of different spectra vary, indicating highly structured HII regions, but define a locus of uniform metallicity. The strong-line diagnostic diagrams and empirical calibrations point to a nearly solar or even super-solar oxygen abundance. The velocity dispersion of the gas is highest in the region offset from the bright clusters. Star formation may be active on a low level. There is evidence for a large-scale disk-like structure in the region of SH2, which would make the similar radial velocity of the nearby globular cluster easier to understand. The high metallicity does not fit to a dwarf galaxy as progenitor. We favour the scenario of a free-floating gaseous complex having its origin in the merger 2 Gyr ago. Over a long period the densities increased secularly until finally the threshold for star formation was reached. SH2 illustrates how massive star clusters can form outside starbursts and without a considerable field population.

Mrk 1216 is a nearby, early-type galaxy with a small effective radius of 2.8 kpc and a large stellar velocity dispersion of 308 km/s for its K-band luminosity of 1.4x10^11 L_sun. Using integral-field spectroscopy assisted by adaptive optics from Gemini North, we measure spatially resolved stellar kinematics within ~450 pc of the galaxy nucleus. The galaxy exhibits regular rotation with velocities of \pm 180 km/s and a sharply peaked velocity dispersion profile that reaches 425 km/s at the center. We fit axisymmetric, orbit-based dynamical models to the combination of these high angular resolution kinematics, large-scale kinematics extending to roughly three effective radii, and Hubble Space Telescope imaging, resulting in a constraint of the mass of the central black hole in Mrk 1216. After exploring several possible sources of systematics that commonly affect stellar-dynamical black hole mass measurements, we find a black hole mass of (4.9\pm1.7)x10^9 M_sun and a H-band stellar mass-to-light ratio of 1.3\pm0.4 M_sun/L_sun (1-sigma uncertainties). Mrk 1216 is consistent with the local black hole mass - stellar velocity dispersion relation, but is a factor of ~5-10 larger than expectations from the black hole mass - bulge luminosity and black hole mass - bulge mass correlations when conservatively using the galaxy's total luminosity or stellar mass. This behavior is quite similar to the extensively studied compact galaxy NGC 1277. Resembling the z~2 quiescent galaxies, Mrk 1216 may be a passively evolved descendant, and perhaps reflects a previous era when galaxies contained over-massive black holes relative to their bulge luminosities/masses, and the growth of host galaxies had yet to catch up.

We present the results of VLT-MUSE integral field spectroscopy of SECCO1, a faint, star-forming stellar system recently discovered as the stellar counterpart of an Ultra Compact High Velocity Cloud (HVC274.68+74.0), very likely residing within a substructure of the Virgo cluster of galaxies. We have obtained the radial velocity of a total of 38 individual compact sources identified as HII regions in the main and secondary body of the system, and derived the metallicity for 18 of them. We provide the first direct demonstration that the two stellar bodies of SECCO1 are physically associated and that their velocities match the HI velocities. The metallicity is quite uniform over the whole system, with a dispersion sigma_12+log(O/H/)=0.08, lower than the uncertainty on individual metallicity estimates. The mean abundance, 12+log(O/H)=8.44, is much higher than the typical values for local dwarf galaxies of similar stellar mass. This strongly suggests that the SECCO~1 stars were born from a pre-enriched gas cloud, possibly stripped from a larger galaxy. Using archival HST images we derive a total stellar mass of ~1.6 X 10^5 M_sun for SECCO1, confirming that it has a very high HI to stellar mass ratio for a dwarf galaxy, M_HI/M_*~ 100. The star formation rate, derived from the H_alpha flux is a factor of more than 10 higher than in typical dwarf galaxies of similar luminosity.

In this paper we present an overview of the MAGNA (Multiple AGN Activity) project aiming at a comprehensive study of multiple supemassive black hole systems. With the main goal to characterize the sources in merging systems at different stages of evolution, we selected a sample of objects optically classified as multiple systems on the basis of emission line diagnostics and started a massive multiband observational campaign. Here we report on the discovery of the exceptionally high AGN density compact group SDSS~J0959+1259. A multiband study suggests that strong interactions are taking place among its galaxies through tidal forces, therefore this system represents a case study for physical mechanisms that trigger nuclear activity and star formation. We also present a preliminary analysis of the multiple AGN system SDSS~J1038+3921.}

We present new MUSE observations of quasar field Q2131-1207 with a log N(HI)=19.50+/-0.15 sub-DLA at z_abs=0.42980. We detect four galaxies at a redshift consistent with that of the absorber where only one was known before this study. Two of these are star forming galaxies, while the ones further away from the quasar (>140 kpc) are passive galaxies. We report the metallicities of the HII regions of the closest objects (12+log(O/H)=8.98+/-0.02 and 8.32+/-0.16) to be higher or equivalent within the errors to the metallicity measured in absorption in the neutral phase of the gas (8.15+/-0.20). For the closest object, a detailed morpho-kinematic analysis indicates that it is an inclined large rotating disk with V_max=200+/-3 km/s. We measure the masses to be M_dyn=7.4+/-0.4 x 10^10 M_sun and M_halo=2.9+/-0.2 x 10^12 M_sun. Some of the gas seen in absorption is likely to be co-rotating with the halo of that object, possibly due to a warped disk. The azimuthal angle between the quasar line of sight and the projected major axis of the galaxy on the sky is 12+/-1 degrees which indicates that some other fraction of the absorbing gas might be associated with accreting gas. This is further supported by the galaxy to gas metallicity difference. Based on the same arguments, we exclude outflows as a possibility to explain the gas in absorption. The four galaxies form a large structure (at least 200 kpc wide) consistent with a filament or a galaxy group so that a fraction of the absorption could be related to intra-group gas.

We present the stellar kinematic maps of a large sample of galaxies from the integral-field spectroscopic survey CALIFA. The sample comprises 300 galaxies displaying a wide range of morphologies across the Hubble sequence, from ellipticals to late-type spirals. This dataset allows us to homogeneously extract stellar kinematics up to several effective radii. In this paper, we describe the level of completeness of this subset of galaxies with respect to the full CALIFA sample, as well as the virtues and limitations of the kinematic extraction compared to other well-known integral-field surveys. In addition, we provide averaged integrated velocity dispersion radial profiles for different galaxy types, which are particularly useful to apply aperture corrections for single aperture measurements or poorly resolved stellar kinematics of high-redshift sources. The work presented in this paper sets the basis for the study of more general properties of galaxies that will be explored in subsequent papers of the survey.

We recently discovered that the active galactic nucleus (AGN) of Mrk 1018 has changed optical type again after 30 years as a type 1 AGN. Here we combine Chandra, NuStar, Swift, Hubble Space Telescope and ground-based observations to explore the cause of this change. The 2-10keV flux declines by a factor of ~8 between 2010 and 2016. We show with our X-ray observation that this is not caused by varying neutral hydrogen absorption along the line-of-sight up to the Compton-thick level. The optical-UV spectral energy distributions are well fit with a standard geometrically thin optically thick accretion disc model that seems to obey the expected $L\sim T^4$ relation. It confirms that a decline in accretion disc luminosity is the primary origin for the type change. We detect a new narrow-line absorber in Lya blue-shifted by ~700km/s with respect to the systemic velocity of the galaxy. This new Lya absorber could be evidence for the onset of an outflow or a companion black hole with associated gas that could be related to the accretion rate change. However, the low column density of the absorber means that it is not the direct cause for Mrk 1018's changing-look nature.

We report the discovery that the known `changing look' AGN Mrk 1018 has changed spectral type for a second time. New VLT-MUSE data taken in 2015 as part of the Close AGN Reference Survey (CARS) shows that the AGN has returned to its original Seyfert 1.9 classification. The CARS sample is selected to contain only bright type 1 AGN, but Mrk 1018's broad emission lines and continuum, typical of type 1 AGN, have almost entirely disappeared. We use spectral fitting of the MUSE spectrum and previously available spectra to determine the drop in broad line flux and the Balmer decrement. We find that the broad line flux has decreased by a factor of 4.75+/- 0.5 in H{\alpha} since an SDSS spectrum was taken in 2000. The Balmer decrement has not changed significantly implying no enhanced reddening with time, but the remaining broad lines are more asymmetric than those present in the type 1 phase. We posit that the change is due to an intrinsic drop in flux from the accretion disk rather than variable extinction or a tidal disruption event.

We analysed the optical spectra of HII regions extracted from a sample of 350 galaxies of the CALIFA survey. We calculated total O/H abundances and, for the first time, N/O ratios using the semi-empirical routine HII-CHI-mistry, which, according to P\'erez-Montero (2014), is consistent with the direct method and reduces the uncertainty in the O/H derivation using [NII] lines owing to the dispersion in the O/H-N/O relation. Then we performed linear fittings to the abundances as a function of the de-projected galactocentric distances. The analysis of the radial distribution both for O/H and N/O in the non-interacting galaxies reveals that both average slopes are negative, but a non-negligible fraction of objects have a flat or even a positive gradient (at least 10\% for O/H and 4\% for N/O). The slopes normalised to the effective radius appear to have a slight dependence on the total stellar mass and the morphological type, as late low-mass objects tend to have flatter slopes. No clear relation is found, however, to explain the presence of inverted gradients in this sample, and there is no dependence between the average slopes and the presence of a bar. The relation between the resulting O/H and N/O linear fittings at the effective radius is much tighter (correlation coefficient $\rho_s$ = 0.80) than between O/H and N/O slopes ($\rho_s$ = 0.39) or for O/H and N/O in the individual \hii\ regions ($\rho_s$ = 0.37). These O/H and N/O values at the effective radius also correlate very tightly (less than 0.03 dex of dispersion) with total luminosity and stellar mass. The relation with other integrated properties, such as star formation rate, colour, or morphology, can be understood only in light of the found relation with mass.

The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population (60 < v_circ < 170 km/s), we use the HIPASS velocity function. For the massive galaxy population (170 < v_circ < 320 km/s), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume- corrected galaxy distribution functions. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations.

In this letter we describe how we use stellar dynamics information to constrain the shape of the stellar IMF in a sample of 27 early-type galaxies from the CALIFA survey. We obtain dynamical and stellar mass-to-light ratios, $\Upsilon_\mathrm{dyn}$ and $\Upsilon_{\ast}$, over a homogenous aperture of 0.5~$R_{e}$. We use the constraint $\Upsilon_\mathrm{dyn} \ge \Upsilon_{\ast}$ to test two IMF shapes within the framework of the extended MILES stellar population models. We rule out a single power law IMF shape for 75% of the galaxies in our sample. Conversely, we find that a double power law IMF shape with a varying high-mass end slope is compatible (within 1$\sigma$) with 95% of the galaxies. We also show that dynamical and stellar IMF mismatch factors give consistent results for the systematic variation of the IMF in these galaxies.

This paper describes the Third Public Data Release (DR3) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. Science-grade quality data for 667 galaxies are made public, including the 200 galaxies of the Second Public Data Release (DR2). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto Observatory. Three different spectral setups are available, i) a low-resolution V500 setup covering the wavelength range 3749-7500 AA (4240-7140 AA unvignetted) with a spectral resolution of 6.0 AA (FWHM), for 646 galaxies, ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 AA (3650-4620 AA unvignetted) with a spectral resolution of 2.3 AA (FWHM), for 484 galaxies, and iii) the combination of the cubes from both setups (called COMBO), with a spectral resolution of 6.0 AA and a wavelength range between 3700-7500 AA (3700-7140 AA unvignetted), for 446 galaxies. The Main Sample, selected and observed according to the CALIFA survey strategy covers a redshift range between 0.005 and 0.03, spans the color-magnitude diagram and probes a wide range of stellar mass, ionization conditions, and morphological types. The Extension Sample covers several types of galaxies that are rare in the overall galaxy population and therefore not numerous or absent in the CALIFA Main Sample. All the cubes in the data release were processed using the latest pipeline, which includes improved versions of the calibration frames and an even further improved im- age reconstruction quality. In total, the third data release contains 1576 datacubes, including ~1.5 million independent spectra. It is available at http://califa.caha.es/DR3.

Photoionization models of HII regions require as input a description of the ionizing SED and of the gas distribution, in terms of ionization parameter U and chemical abundances (e.g. O/H and N/O). A strong degeneracy exists between the hardness of the SED and U, which in turn leads to high uncertainties in the determination of the other parameters, including abundances. One way to resolve the degeneracy is to fix one of the parameters using additional information. For each of the ~ 20000 sources of the CALIFA HII regions catalog, a grid of photoionization models is computed assuming the ionizing SED being described by the underlying stellar population obtained from spectral synthesis modeling. The ionizing SED is then defined as the sum of various stellar bursts of different ages and metallicities. This solves the degeneracy between the shape of the ionizing SED and U. The nebular metallicity (associated to O/H) is defined using the classical strong line method O3N2 (which gives to our models the status of "hybrids"). The remaining free parameters are the abundance ratio N/O and the ionization parameter U, which are determined by looking for the model fitting [NII]/Ha and [OIII]/Hb. The models are also selected to fit [OII]/Hb. This process leads to a set of ~ 3200 models that reproduce simultaneously the three observations. We find that the regions associated to young stellar bursts suffer leaking of the ionizing photons, the proportion of escaping photons having a median of 80\%. The set of photoionization models satisfactorily reproduces the electron temperature derived from the [OIII]4363/5007 line ratio. We determine new relations between the ionization parameter U and the [OII]/[OIII] or [SII]/[SIII] line ratios. New relations between N/O and O/H and between U and O/H are also determined. All the models are publicly available on the 3MdB database.

[Abridged] Enormous observational efforts have been made to constrain the energetics of AGN feedback by mapping the kinematics of the ionized gas on kpc scale. We study how the observed kinematics and inferred energetics are affected by beam smearing of a bright unresolved NLR due to seeing. We analyse IFU spectroscopy of a sample of twelve QSOs initially presented by Liu et al. (2014). The PSF for the observations is directly obtained from the light distribution of the broad Hb line component. We are able to compare the ionized gas kinematics and derived energetics of the total, spatially extended, and unresolved [OIII] emission. We find that the spatially resolved [OIII] line width on kpc scales is significantly narrower than the one before PSF deblending. The ENLRs appear offset from the QSO position or more elongated which can be interpreted in favour of a conical outflow on large scales while a spherical geometry cannot be excluded for the unresolved NLR. We find that the kinetic power at 5kpc distance from the spherical model by Liu et al. (2013) is reduced by two magnitudes for a conical outflow and one magnitude for the unresolved NLR after PSF deblending. This reduced kinetic power corresponds to 0.01-0.1% of the bolometric AGN luminosity. This is smaller than the 5-10% feedback efficiency required by some simulations to reproduce the massive galaxy population. The injected momentum fluxes are close or below the simple radiation-pressure limit Lbol/c for the conical outflow model for the NLR and ENLR when beam smearing is considered. IFU spectroscopy is a powerful tool to investigate the energetics of AGN outflows, but the impact of beam smearing has to be taken into account. For the majority of observations in the literature, this has not been addressed carefully so that the incidence and energetics of presumed kpc-scale AGN-driven outflows still remain an unsolved issue.

AGN driven outflows are invoked in numerical simulations to reproduce several observed properties of local galaxies. The z > 1 epoch is of particular interest as it was during this time that the volume averaged star formation and the accretion rate of black holes were maximum. Radiatively driven outflows are therefore believed to be common during this epoch. We aim to trace and characterize outflows in AGN hosts with high mass accretion rates at z > 1 using integral field spectroscopy. We obtain spatially-resolved kinematics of the [OIII]5007 line in two targets which reveal the morphology and spatial extension of the outflows. We present J and H+K band SINFONI observations of 5 AGNs at 1.2 < z < 2.2. To maximize the chance of observing radiatively driven outflows, our sample was pre-selected based on peculiar values of the Eddington ratio and the hydrogen column density of the surrounding interstellar medium. We observe high velocity (~600-1900 km/s) and kiloparsec scale extended ionized outflows in at least 3 of our targets, using [OIII]5007 line kinematics tracing the AGN narrow line region. We estimate the total mass of the outflow, the mass outflow rate, and the kinetic power of the outflows based on theoretical models and report on the uncertainties associated with them. We find mass outflow rates of ~1-10 M_sun/yr for the sample presented in this paper. Based on the high star formation rates of the host galaxies, the observed outflow kinetic power and the expected power due to the AGN, we infer that both star formation and AGN radiation could be the dominant source for the outflows. The outflow models suffer from large uncertainties, hence we call for further detailed observations for an accurate determination of the outflow properties to confirm the exact source of these outflows.

This paper aims at providing aperture corrections for emission lines in a sample of spiral galaxies from the Calar Alto Legacy Integral Field Area Survey (CALIFA) database. In particular, we explore the behavior of the log([OIII]5007/Hbeta)/([NII]6583/Halpha) (O3N2) and log[NII]6583/Halpha (N2) flux ratios since they are closely connected to different empirical calibrations of the oxygen abundances in star forming galaxies. We compute median growth curves of Halpha, Halpha/Hbeta, O3N2 and N2 up to 2.5R_50 and 1.5 disk R_eff. The growth curves simulate the effect of observing galaxies through apertures of varying radii. The median growth curve of the Halpha/Hbeta ratio monotonically decreases from the center towards larger radii, showing for small apertures a maximum value of ~10% larger than the integrated one. The median growth curve of N2 shows a similar behavior, decreasing from the center towards larger radii. No strong dependence is seen with the inclination, morphological type and stellar mass for these growth curves. Finally, the median growth curve of O3N2 increases monotonically with radius. However, at small radii it shows systematically higher values for galaxies of earlier morphological types and for high stellar mass galaxies. Applying our aperture corrections to a sample of galaxies from the SDSS survey at 0.02<=z<=0.3 shows that the average difference between fiber-based and aperture corrected oxygen abundances, for different galaxy stellar mass and redshift ranges, reaches typically to ~11%, depending on the abundance calibration used. This average difference is found to be systematically biased, though still within the typical uncertainties of oxygen abundances derived from empirical calibrations. Caution must be exercised when using observations of galaxies for small radii (e.g. below 0.5R_eff) given the high dispersion shown around the median growth curves.

We present MUSE observations of the debated ultra faint stellar system Crater. We spectroscopically confirm 26 member stars of this system via radial velocity measurements. We derive the systematic instrumental velocity uncertainty of MUSE spectra to be 2.27$\rm \,km\,s^{-1}$. This new dataset increases the confirmed member stars of Crater by a factor of 3. One out of three bright blue stars and a fainter blue star just above the main-sequence-turn-off are also found to be likely members of the system. The observations reveal that Crater has a systemic radial velocity of $v_{\rm sys}=148.18^{\rm +1.08}_{\rm -1.15}\rm \,km\,s^{-1}$, whereas the most likely velocity dispersion of this system is $\sigma_{\rm v}=2.04^{\rm +2.19}_{\rm -1.06} \rm \,km\,s^{-1}$. The total dynamical mass of the system, assuming dynamical equilibrium is then $M_{\rm tot}=1.50^{+4.9}_{-1.2}\cdot 10^{\rm 5}M_{\odot}$ implying a mass-to-light ratio of M/L$_{\rm V}$=8.52$^{+28.0}_{-6.5}\, M_{\odot}/L_{\odot}$, which is consistent with a purely baryonic stellar population within its errors and no significant evidence for the presence dark matter was found. We also find evidence for a velocity gradient in the radial velocity distribution. We conclude that our findings strongly support that Crater is a faint intermediate-age outer halo globular cluster and not a dwarf galaxy.

The "main sequence of galaxies" $-$ defined in terms of the total star formation rate $\psi$ vs. the total stellar mass $M_*$ $-$ is a well-studied tight relation that has been observed at several wavelengths and at different redshifts. All earlier studies have derived this relation from integrated properties of galaxies. We recover the same relation from an analysis of spatially-resolved properties, with integral field spectroscopic (IFS) observations of 306 galaxies from the CALIFA survey. We consider the SFR surface density in units of log(M$_{\odot}$ yr$^{-1}$ Kpc$^{-2}$) and the stellar mass surface density in units of log(M$_{\odot}$ Kpc$^{-2}$) in individual spaxels which probe spatial scales of 0.5-1.5 Kpc. This local relation exhibits a high degree of correlation with small scatter ($\sigma = 0.23$ dex), irrespective of the dominant ionisation source of the host galaxy or its integrated stellar mass. We highlight: $(i)$ the integrated star formation main sequence formed by galaxies whose dominant ionisation process is related to star formation, for which we find a slope of 0.81 $\pm 0.02$; (ii) the spatially-resolved relation obtained with the spaxel analysis, we find a slope of 0.72 $\pm 0.04$; (iii) for the integrated main sequence we identified also a sequence formed by galaxies that are dominated by an old stellar population, which we have called the retired galaxies sequence.

We measured the gas abundance profiles in a sample of 122 face-on spiral galaxies observed by the CALIFA survey and included all spaxels whose line emission was consistent with star formation. This type of analysis allowed us to improve the statistics with respect to previous studies, and to properly estimate the oxygen distribution across the entire disc to a distance of up to 3-4 disc effective radii (r$_e$). We confirm the results obtained from classical HII region analysis. In addition to the general negative gradient, an outer flattening can be observed in the oxygen abundance radial profile. An inner drop is also found in some cases. There is a common abundance gradient between 0.5 and 2.0 r$_e$ of $\alpha_{O/H} = -\,0.075\,\rm{dex}/r_e$ with a scatter of $\sigma = 0.016\,\rm{dex}/r_e$ when normalising the distances to the disc effective radius. By performing a set of Kolmogorov-Smirnov tests, we determined that this slope is independent of other galaxy properties, such as morphology, absolute magnitude, and the presence or absence of bars. In particular, barred galaxies do not seem to display shallower gradients, as predicted by numerical simulations. Interestingly, we find that most of the galaxies in the sample with reliable oxygen abundance values beyond $\sim 2$ effective radii (57 galaxies) present a flattening of the abundance gradient in these outer regions. This flattening is not associated with any morphological feature, which suggests that it is a common property of disc galaxies. Finally, we detect a drop or truncation of the abundance in the inner regions of 27 galaxies in the sample; this is only visible for the most massive galaxies.

This is the fourth paper in a series that reports on our investigation of the clustering properties of active galactic nuclei (AGN) identified in the ROSAT All-Sky Survey (RASS) and Sloan Digital Sky Survey (SDSS). In this paper we investigate the cause of the X-ray luminosity dependence of the clustering of broad-line, luminous AGN at 0.16<z<0.36. We fit the H-alpha line profile in the SDSS spectra for all X-ray and optically-selected broad-line AGN, determine the mass of the super-massive black hole (SMBH), M_BH, and infer the accretion rate relative to Eddington (L/L_EDD). Since M_BH and L/L_EDD are correlated, we create AGN subsamples in one parameter while maintaining the same distribution in the other parameter. In both the X-ray and optically-selected AGN samples we detect a weak clustering dependence with M_BH and no statistically significant dependence on L/L_EDD. We find a difference of up to 2.7sigma when comparing the objects that belong to the 30% least and 30% most massive M_BH subsamples, in that luminous broad-line AGN with more massive black holes reside in more massive parent dark matter halos at these redshifts. These results provide evidence that higher accretion rates in AGN do not necessarily require dense galaxy environments in which more galaxy mergers and interactions are expected to channel large amounts of gas onto the SMBH. We also present semi-analytic models which predict a positive M_DMH dependence on M_BH, which is most prominent at M_BH ~ 10^{8-9} M_SUN.

The nearby lenticular galaxy NGC 1277 is thought to host one of the largest black holes known, however the black hole mass measurement is based on low spatial resolution spectroscopy. In this paper, we present Gemini Near-infrared Integral Field Spectrometer observations assisted by adaptive optics. We map out the galaxy's stellar kinematics within ~440 pc of the nucleus with an angular resolution that allows us to probe well within the region where the potential from the black hole dominates. We find that the stellar velocity dispersion rises dramatically, reaching ~550 km/s at the center. Through orbit-based, stellar-dynamical models we obtain a black hole mass of (4.9 \pm 1.6) x 10^9 Msun (1-sigma uncertainties). Although the black hole mass measurement is smaller by a factor of ~3 compared to previous claims based on large-scale kinematics, NGC 1277 does indeed contain one of the most massive black holes detected to date, and the black hole mass is an order of magnitude larger than expectations from the empirical relation between black hole mass and galaxy luminosity. Given the galaxy's similarities to the higher redshift (z~2) massive quiescent galaxies, NGC 1277 could be a relic, passively evolving since that period. A population of local analogs to the higher redshift quiescent galaxies that also contain over-massive black holes may suggest that black hole growth precedes that of the host galaxy.

We study the stellar properties of 44 face-on spiral galaxies from the Calar Alto Legacy Integral Field Area survey via full spectrum fitting techniques. We compare the age profiles with the surface brightness distribution in order to highlight differences between profile types (type I, exponential profile; and II, down-bending profile). We observe an upturn ("U-shape") in the age profiles for 17 out of these 44 galaxies with reliable stellar information up to their outer parts. This "U-shape" is not a unique feature for type II galaxies but can be observed in type I as well. These findings suggest that the mechanisms shaping the surface brightness and stellar population distributions are not directly coupled. This upturn in age is only observable in the light-weighted profiles while it flattens out in the mass-weighted profiles. Given recent results on the outer parts of nearby systems and the results presented in this Letter, one of the most plausible explanations for the age upturn is an early formation of the entire disc ($\sim$~10~Gyr ago) followed by an inside-out quenching of the star formation.

We investigate the compact, early-type galaxy NGC 1281 with integral field unit observations to map the stellar LOSVD out to 5 effective radii and construct orbit-based dynamical models to constrain its dark and luminous matter content. Under the assumption of mass-follows-light, the H-band stellar mass-to-light ratio (M/L) is {\Upsilon} = 2.7(+-0.1) {\Upsilon}_{sun}, higher than expected from our stellar population synthesis fits with either a canonical Kroupa ({\Upsilon} = 1.3 {\Upsilon}_{sun}) or Salpeter ({\Upsilon} = 1.7 {\Upsilon}_{sun}) stellar initial mass function. Such models also cannot reproduce the details of the LOSVD. Models with a dark halo recover the kinematics well and indicate that NGC 1281 is dark matter dominated, making up ~ 90 per cent of the total enclosed mass within the kinematic bounds. Parameterised as a spherical NFW profile, the dark halo mass is 11.5 < log(M_{DM}/M_{sun}) < 11.8 and the stellar M/L is 0.6 < {\Upsilon} < 1.1. However, this stellar M/L is lower than predicted by its old stellar population. Moreover, the halo mass within the kinematic extent is ten times larger than expected based on {\Lambda}CDM predictions, and an extrapolation yields cluster sized dark halo masses. Adopting {\Upsilon} = 1.7 {\Upsilon}_{sun} yields more moderate dark halo virial masses, but these models fit the kinematics worse. A non-NFW model might solve the discrepancy between the unphysical consequences of the best-fitting dynamical models and models based on more reasonable assumptions for the dark halo and stellar mass-to-light ratio, which are disfavoured according to our parameter estimation.

We report deep optical integral-field spectroscopy with the MUSE of the luminous radio-quiet QSO PG1307+085 (z=0.154) obtained during the commissioning of the instrument. Given the high sensitivity and spatial resolution delivered by MUSE, we are able to resolve the compact ($r_e$~1.3") elliptical host galaxy. After careful spectroscopic deblending of the QSO and host galaxy emission, we infer a stellar velocity dispersion of $155\pm19$km/s. This places PG1307+085 local $M_{BH}-\sigma_*$ relation within the intrinsic scatter but offset towards a higher black hole mass with respect to the mean relation. The observations with MUSE also reveal a large extended ENLR around PG1307+085 reaching out to 30kpc. In addition, we detect a faint bridge of ionized gas towards the most massive galaxy of the galaxy group being just 20" (50kpc) away. Previous long-slit spectroscopic observations missed most of these extended features due to a miss-aligned slit. The ionized gas kinematics does not show any evidence for gas outflows on kpc scales despite the high QSO luminosity of $L_\mathrm{bol}>10^{46}$ erg/s. Based on the ionized gas distribution, kinematics and metallicity we discuss the origin of the ENLR with respect to its group environments including minor mergers, ram-pressure stripping or filamentary gas accretion as the most likely scenarios. We conclude that PG1307+085 is a normal elliptical host in terms of the scaling relations, but that the gas is most likely affected by the environment through gravity or ambient pressure. It is possible that the ongoing interaction with the environment, mainly seen in the ionized gas, is also be responsible for driving sufficient gas to feed the black hole at the centre of the galaxy.

We use data from the Multi-Unit Spectroscopic Explorer (MUSE), recently commissioned at the Very Large Telescope (VLT), to study the kinematics and stellar population content of NGC 4371, an early-type massive barred galaxy in the core of the Virgo cluster. We integrate this study with a detailed structural analysis using imaging data from the Hubble and Spitzer space telescopes, which allows us to perform a thorough investigation of the physical properties of the galaxy. We show that the rotationally supported inner components in NGC 4371, an inner disc and a nuclear ring - which, according to the predominant scenario, are built with stars formed from gas brought to the inner region by the bar - are vastly dominated by stars older than 10 Gyr. Our results thus indicate that the formation of the bar occurred at a redshift of about $z=1.8^{+0.5}_{-0.4}$ (error bars are derived from 100 Monte Carlo realisations). NGC 4371 thus testifies to the robustness of bars. In addition, the mean stellar age of the fraction of the major disc of the galaxy covered by our MUSE data is above 7 Gyr, with a small contribution from younger stars. This suggests that the quenching of star formation in NGC 4371, likely due to environmental effects, was already effective at a redshift of about $z=0.8^{+0.2}_{-0.1}$. Our results point out that bar-driven secular evolution processes may have an extended impact in the evolution of galaxies, and thus on the properties of galaxies as observed today, not necessarily restricted to more recent cosmic epochs.

We present a multi-wavelength study of a newly discovered compact group (CG), SDSS J0959+1259, based data from XMM-Newton, SDSS and the Calar Alto optical imager BUSCA. With a maximum velocity offset of 500 km s$^{-1}$, a mean redshift of 0.035, and a mean spatial extension of 480 kpc, this CG is exceptional in having the highest concentration of nuclear activity in the local Universe, established with a sensitivity limit L$_{X}>4\times $10$^{40}$ erg s$^{-1}$ in 2--10 keV band and R-band magnitude $M_R < -19$. The group is composed of two type-2 Seyferts, one type-1 Seyfert, two LINERs and three star forming galaxies. Given the high X-ray luminosity of LINERs which reaches $\sim 10^{41}$ erg s$^{-1}$, it is likely that they are also accretion driven, bringing the number of active nuclei in this group to to 5 out of 8 (AGN fraction of 60\%). The distorted shape of one member of the CG suggests that strong interactions are taking place among its galaxies through tidal forces. Therefore, this system represents a case study for physical mechanisms that trigger nuclear activity and star formation in CGs.

The Star Formation Rate (SFR) is one of the main parameters used to analyze the evolution of galaxies through time. The need for recovering the light reprocessed by dust commonly requires the use of low spatial resolution far-infrared data. Recombination-line luminosities provide an alternative, although uncertain dust-extinction corrections based on narrow-band imaging or long-slit spectroscopy have traditionally posed a limit to their applicability. Integral Field Spectroscopy (IFS) is clearly the way to overcome such limitation. We obtain integrated H{\alpha}, ultraviolet (UV) and infrared (IR)-based SFR measurements for 272 galaxies from the CALIFA survey at 0.005 < z < 0.03 using single-band and hybrid tracers. We provide updated calibrations, both global and split by properties (including stellar mass and morphological type), referred to H{\alpha}. The extinction-corrected H{\alpha} luminosity agrees with the updated hybrid SFR estimators based on either UV or H{\alpha} plus IR luminosity over the full range of SFRs (0.03-20 M$_{\odot}$ yr$^{-1}$). The coefficient that weights the amount of energy produced by newly-born stars that is reprocessed by dust on the hybrid tracers, a$_{IR}$, shows a large dispersion. However, it does not became increasingly small at high attenuations, as expected if significant highly-obscured H$\alpha$ emission would be missed. Lenticulars, early-type spirals and type-2 AGN host galaxies show smaller coefficients due to the contribution of optical photons and AGN to dust heating. In the Local Universe the H{\alpha} luminosity derived from IFS observations can be used to measure SFR, at least in statistically-significant, optically-selected galaxy samples. The analysis of the SFR calibrations by galaxies properties could be potentially used by other works to study the impact of different selection criteria in the SFR values derived.

We present a dynamical analysis to infer the structural parameters and properties of the two nearby, compact, high velocity dispersion galaxies MRK1216 & NGC1277. Combining deep HST imaging, wide-field IFU stellar kinematics, and complementary long-slit spectroscopic data out to 3 R_e, we construct orbit-based models to constrain their black hole masses, dark matter content and stellar mass-to-light ratios. We obtain a black hole mass of log(Mbh/Msun) = 10.1(+0.1/-0.2) for NGC1277 and an upper limit of log(Mbh/Msun) = 10.0 for MRK1216, within 99.7 per cent confidence. The stellar mass-to-light ratios span a range of Upsilon_V = 6.5(+1.5/-1.5) in NGC1277 and Upsilon_H = 1.8(+0.5/-0.8) in MRK1216 and are in good agreement with SSP models of a single power-law Salpeter IMF. Even though our models do not place strong constraints on the dark halo parameters, they suggest that dark matter is a necessary ingredient in MRK1216, with a dark matter contribution of 22(+30/-20) per cent to the total mass budget within 1 R_e. NGC1277, on the other hand, can be reproduced without the need for a dark halo, and a maximal dark matter fraction of 13 per cent within the same radial extent. In addition, we investigate the orbital structures of both galaxies, which are rotationally supported and consistent with photometric multi-S\'ersic decompositions, indicating that these compact objects do not host classical, non-rotating bulges formed during recent (z <= 2) dissipative events or through violent relaxation. Finally, both MRK 1216 and NGC 1277 are anisotropic, with a global anisotropy parameter delta of 0.33 and 0.58, respectively. While MRK 1216 follows the trend of fast-rotating, oblate galaxies with a flattened velocity dispersion tensor in the meridional plane of the order of beta_z = delta, NGC 1277 is highly tangentially anisotropic and seems to belong kinematically to a distinct class of objects.

Located in the Perseus cluster, NGC 1271 is an early-type galaxy with a small effective radius of 2.2 kpc and a large stellar velocity dispersion of 276 km/s for its K-band luminosity of 8.9x10^{10} L_sun. We present a mass measurement for the black hole in this compact, high-dispersion galaxy using observations from the integral field spectrograph NIFS on the Gemini North telescope assisted by laser guide star adaptive optics, large-scale integral field unit observations with PPAK at the Calar Alto Observatory, and Hubble Space Telescope WFC3 imaging observations. We are able to map out the stellar kinematics on small spatial scales, within the black hole sphere of influence, and on large scales that extend out to four times the galaxy's effective radius. We find that the galaxy is rapidly rotating and exhibits a sharp rise in the velocity dispersion. Through the use of orbit-based stellar dynamical models, we determine that the black hole has a mass of (3.0^{+1.0}_{-1.1}) x 10^9 M_sun and the H-band stellar mass-to-light ratio is 1.40^{+0.13}_{-0.11} M_sun/L_sun (1-sigma uncertainties). NGC 1271 occupies the sparsely-populated upper end of the black hole mass distribution, but is very different from the Brightest Cluster Galaxies (BCGs) and giant elliptical galaxies that are expected to host the most massive black holes. Interestingly, the black hole mass is an order of magnitude larger than expectations based on the galaxy's bulge luminosity, but is consistent with the mass predicted using the galaxy's bulge stellar velocity dispersion. More compact, high-dispersion galaxies need to be studied using high spatial resolution observations to securely determine black hole masses, as there could be systematic differences in the black hole scaling relations between these types of galaxies and the BCGs/giant ellipticals, thereby implying different pathways for black hole and galaxy growth.

This paper characterizes the radial structure of stellar population properties of galaxies in the nearby universe, based on 300 galaxies from the CALIFA survey. The sample covers a wide range of Hubble types, and galaxy stellar mass. We apply the spectral synthesis techniques to recover the stellar mass surface density, stellar extinction, light and mass-weighted ages, and mass-weighted metallicity, for each spatial resolution element in our target galaxies. To study mean trends with overall galaxy properties, the individual radial profiles are stacked in seven bins of galaxy morphology. We confirm that more massive galaxies are more compact, older, more metal rich, and less reddened by dust. Additionally, we find that these trends are preserved spatially with the radial distance to the nucleus. Deviations from these relations appear correlated with Hubble type: earlier types are more compact, older, and more metal rich for a given mass, which evidences that quenching is related to morphology, but not driven by mass. Negative gradients of ages are consistent with an inside-out growth of galaxies, with the largest ages gradients in Sb-Sbc galaxies. Further, the mean stellar ages of disks and bulges are correlated, with disks covering a wider range of ages, and late type spirals hosting younger disks. The gradients in stellar mass surface density depend mostly on stellar mass, in the sense that more massive galaxies are more centrally concentrated. There is a secondary correlation in the sense that at the same mass early type galaxies have steeper gradients. We find mildly negative metallicity gradients, shallower than predicted from models of galaxy evolution in isolation. The largest gradients occur in Sb galaxies. Overall we conclude that quenching processes act in manners that are independent of mass, while metallicity and galaxy structure are influenced by mass-dependent processes.

We present spatially resolved stellar and/or ionized gas kinematic properties for a sample of 103 interacting galaxies, tracing all merger stages: close companions, pairs with morphological signatures of interaction, and coalesced merger remnants. We compare our sample with 80 non-interacting galaxies. We measure for the stellar and the ionized gas components the major (projected) kinematic position angles (PA$_{\mathrm{kin}}$, approaching and receding) directly from the velocity fields with no assumptions on the internal motions. This method allow us to derive the deviations of the kinematic PAs from a straight line ($\delta$PA$_{\mathrm{kin}}$). Around half of the interacting objects show morpho-kinematic PA misalignments that cannot be found in the control sample. Those misalignments are present mostly in galaxies with morphological signatures of interaction. Alignment between the kinematic sides for both samples is similar, with most of the galaxies displaying small misalignments. Radial deviations of the kinematic PA from a straight line in the stellar component measured by $\delta$PA$_{\mathrm{kin}}$ are large for both samples. However, for a large fraction of interacting galaxies the ionized gas $\delta$PA$_{\mathrm{kin}}$ is larger than typical values derived from isolated galaxies (48%), making this parameter a good indicator to trace the impact of interaction and mergers in the internal motions of galaxies. By comparing the stellar and ionized gas kinematic PA, we find that 42% (28/66) of the interacting galaxies have misalignments larger than 16 degrees, compared to 10% from the control sample. Our results show the impact of interactions in the internal structure of galaxies as well as the wide variety of their velocity distributions. This study also provides a local Universe benchmark for kinematic studies in merging galaxies at high redshift.

We use optical integral-field spectroscopic (IFS) data from 103 nearby galaxies at different stages of the merging event, from close pairs to merger remnants provided by the CALIFA survey, to study the impact of the interaction in the specific star formation and oxygen abundance on different galactic scales. To disentangle the effect of the interaction and merger from internal processes, we compared our results with a control sample of 80 non-interacting galaxies. We confirm the moderate enhancement (2-3 times) of specific star formation for interacting galaxies in central regions as reported by previous studies; however, the specific star formation is comparable when observed in extended regions. We find that control and interacting star-forming galaxies have similar oxygen abundances in their central regions, when normalized to their stellar masses. Oxygen abundances of these interacting galaxies seem to decrease compared to the control objects at the large aperture sizes measured in effective radius. Although the enhancement in central star formation and lower metallicities for interacting galaxies have been attributed to tidally induced inflows, our results suggest that other processes such as stellar feedback can contribute to the metal enrichment in interacting galaxies.

While studies of gas-phase metallicity gradients in disc galaxies are common, very little has been done in the acquisition of stellar abundance gradients in the same regions. We present here a comparative study of the stellar metallicity and age distributions in a sample of 62 nearly face-on, spiral galaxies with and without bars, using data from the CALIFA survey. We measure the slopes of the gradients and study their relation with other properties of the galaxies. We find that the mean stellar age and metallicity gradients in the disc are shallow and negative. Furthermore, when normalized to the effective radius of the disc, the slope of the stellar population gradients does not correlate with the mass or with the morphological type of the galaxies. Contrary to this, the values of both age and metallicity at $\sim$2.5 scale-lengths correlate with the central velocity dispersion in a similar manner to the central values of the bulges, although bulges show, on average, older ages and higher metallicities than the discs. One of the goals of the present paper is to test the theoretical prediction that non-linear coupling between the bar and the spiral arms is an efficient mechanism for producing radial migrations across significant distances within discs. The process of radial migration should flatten the stellar metallicity gradient with time and, therefore, we would expect flatter stellar metallicity gradients in barred galaxies. However, we do not find any difference in the metallicity or age gradients in galaxies with without bars. We discuss possible scenarios that can lead to this absence of difference.

Supernova remnants are among the most spectacular examples of astrophysical pistons in our cosmic neighborhood. The gas expelled by the supernova explosion is launched with velocities ~ 1000 km/s into the ambient, tenuous interstellar medium, producing shocks that excite hydrogen lines. We have used an optical integral-field spectrograph to obtain high-resolution spatial-spectral maps that allow us to study in detail the shocks in the northwestern rim of supernova 1006. The two-component H\alpha\ line is detected at 133 sky locations. Variations in the broad line widths and the broad-to-narrow line intensity ratios across tens of atomic mean free paths suggest the presence of suprathermal protons, the potential seed particles for generating high-energy cosmic-rays.

The bar pattern speed ($\Omega_{\rm b}$) is defined as the rotational frequency of the bar, and it determines the bar dynamics. Several methods have been proposed for measuring $\Omega_{\rm b}$. The non-parametric method proposed by Tremaine \& Weinberg (1984; TW) and based on stellar kinematics is the most accurate. This method has been applied so far to 17 galaxies, most of them SB0 and SBa types. We have applied the TW method to a new sample of 15 strong and bright barred galaxies, spanning a wide range of morphological types from SB0 to SBbc. Combining our analysis with previous studies, we investigate 32 barred galaxies with their pattern speed measured by the TW method. The resulting total sample of barred galaxies allows us to study the dependence of $\Omega_{\rm b}$ on galaxy properties, such as the Hubble type. We measured $\Omega_{\rm b}$ using the TW method on the stellar velocity maps provided by the integral-field spectroscopy data from the CALIFA survey. Integral-field data solve the problems that long-slit data present when applying the TW method, resulting in the determination of more accurate $\Omega_{\rm b}$. In addition, we have also derived the ratio $\cal{R}$ of the corotation radius to the bar length of the galaxies. According to this parameter, bars can be classified as fast ($\cal{R}$ $< 1.4$) and slow ($\cal{R}$>1.4). For all the galaxies, $\cal{R}$ is compatible within the errors with fast bars. We cannot rule out (at 95$\%$ level) the fast bar solution for any galaxy. We have not observed any significant trend between $\cal{R}$ and the galaxy morphological type. Our results indicate that independent of the Hubble type, bars have been formed and then evolve as fast rotators. This observational result will constrain the scenarios of formation and evolution of bars proposed by numerical simulations.

Circumnuclear star formation and AGN feedback is believed to play a critical role in the context of galaxy evolution. The low-luminosity QSO (LLQSO) sample that contains 99 of the closest AGN with redshift z<=0.06 fills the gap between the local AGN population and high-redshift QSOs that is essential to understand the AGN evolution with redshift. In this paper, we present the results of near-infrared H+K-integral field spectroscopy of the inner kiloparsecs of the LLQSO HE 1029-1831 with SINFONI. Line maps show that ionized hydrogen gas is located in spiral arms within the stellar bar and in a circumnuclear ring. Line fluxes and diagnostic line ratios indicate recent or ongoing star formation in the circumnuclear region and the presence of young and intermediate-age stellar populations in the bulge. In particular, we find traces of an intense starburst in the circumnuclear region that has begun around 100 Myr ago but has declined to a fraction of the maximum intensity now. We estimate the dynamical bulge mass and find that the galaxy follows published M_BH-M_bulge relations. However, bulge-disk decomposition of the K-band image with BUDDA reveals that HE 1029-1831 does not follow the M_BH-L_bulge relations of inactive galaxies. We conclude that the deviation from M_BH-L_bulge relations of inactive galaxies in this source is rather caused by young stellar populations and not by an undermassive black hole.

This paper describes the Second Public Data Release (DR2) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. The data for 200 objects are made public, including the 100 galaxies of the First Public Data Release (DR1). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5 m telescope at the Calar Alto observatory. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the wavelength range 3745-7500 \AA\ with a spectral resolution of 6.0 \AA\ (FWHM), and (ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 \AA\ with a spectral resolution of 2.3 \AA\ (FWHM). The sample covers a redshift range between 0.005 and 0.03, with a wide range of properties in the Color-Magnitude diagram, stellar mass, ionization conditions, and morphological types. All released cubes were reduced with the latest pipeline, including improved spectrophotometric calibration, spatial registration and spatial resolution. The spectrophotometric calibration is better than 6% and the median spatial resolution is 2.5". Altogether the second data release contains over 1.5 million spectra. It is available at http://califa.caha.es/DR2.

This work provides an overall characterization of the kinematic behavior of the ionized gas of the galaxies included in the Calar Alto Legacy Integral field Area (CALIFA), offering kinematic clues to potential users of this survey for including kinematical criteria for specific studies. From the first 200 galaxies observed by CALIFA, we present the 2D kinematic view of the 177 galaxies satisfying a gas detection threshold. After removing the stellar contribution, we used the cross-correlation technique to obtain the radial velocity of the dominant gaseous component. The main kinematic parameters were directly derived from the radial velocities with no assumptions on the internal motions. Evidence of the presence of several gaseous components with different kinematics were detected by using [OIII] profiles. Most objects in the sample show regular velocity fields, although the ionized-gas kinematics are rarely consistent with simple coplanar circular motions. 35% of the objects present evidence of a displacement between the photometric and kinematic centers larger than the original spaxel radii. Only 17% of the objects in the sample exhibit kinematic lopsidedness when comparing receding and approaching sides of the velocity fields, but most of them are interacting galaxies exhibiting nuclear activity. Early-type galaxies in the sample present clear photometric-kinematic misaligments. There is evidence of asymmetries in the emission line profiles suggesting the presence of kinematically distinct gaseous components at different distances from the nucleus. This work constitutes the first determination of the ionized gas kinematics of the galaxies observed in the CALIFA survey. The derived velocity fields, the reported kinematic peculiarities and the identification of the presence of several gaseous components might be used as additional criteria for selecting galaxies for specific studies.

We describe and discuss the selection procedure and statistical properties of the galaxy sample used by the Calar Alto Legacy Integral Field Area Survey (CALIFA), a public legacy survey of 600 galaxies using integral field spectroscopy. The CALIFA "mother sample" was selected from the Sloan Digital Sky Survey (SDSS) DR7 photometric catalogue to include all galaxies with an r-band isophotal major axis between 45" and 79.2" and with a redshift 0.005 < z < 0.03. The mother sample contains 939 objects, 600 of which will be observed in the course of the CALIFA survey. The selection of targets for observations is based solely on visibility and thus keeps the statistical properties of the mother sample. By comparison with a large set of SDSS galaxies, we find that the CALIFA sample is representative of galaxies over a luminosity range of -19 > Mr > -23.1 and over a stellar mass range between 10^9.7 and 10^11.4Msun. In particular, within these ranges, the diameter selection does not lead to any significant bias against - or in favour of - intrinsically large or small galaxies. Only below luminosities of Mr = -19 (or stellar masses < 10^9.7Msun) is there a prevalence of galaxies with larger isophotal sizes, especially of nearly edge-on late-type galaxies, but such galaxies form < 10% of the full sample. We estimate volume-corrected distribution functions in luminosities and sizes and show that these are statistically fully compatible with estimates from the full SDSS when accounting for large-scale structure. We also present a number of value-added quantities determined for the galaxies in the CALIFA sample. We explore different ways of characterizing the environments of CALIFA galaxies, finding that the sample covers environmental conditions from the field to genuine clusters. We finally consider the expected incidence of active galactic nuclei among CALIFA galaxies.

We use spatially and temporally resolved maps of stellar population properties of 300 galaxies from the CALIFA integral field survey to investigate how the stellar metallicity (Z*) relates to the total stellar mass (M*) and the local mass surface density ($\mu$*) in both spheroidal and disk dominated galaxies. The galaxies are shown to follow a clear stellar mass-metallicity relation (MZR) over the whole 10$^9$ to 10$^{12}$ M$_{\odot}$ range. This relation is steeper than the one derived from nebular abundances, which is similar to the flatter stellar MZR derived when we consider only young stars. We also find a strong relation between the local values of $\mu$* and Z* (the $\mu$ZR), betraying the influence of local factors in determining Z*. This shows that both local ($\mu$*-driven) and global (M*-driven) processes are important in determining the metallicity in galaxies. We find that the overall balance between local and global effects varies with the location within a galaxy. In disks, $\mu$* regulates Z*, producing a strong $\mu$ZR whose amplitude is modulated by M*. In spheroids it is M* who dominates the physics of star formation and chemical enrichment, with $\mu$* playing a minor, secondary role. These findings agree with our previous analysis of the star formation histories of CALIFA galaxies, which showed that mean stellar ages are mainly governed by surface density in galaxy disks and by total mass in spheroids.

[abridged] We present optical integral field spectroscopy for a flux-limited sample of 19 QSOs at z<0.2 and spatially resolve their ionized gas properties at a physical resolution of 2-5kpc. The extended narrow line regions (ENLRs), photoionized by the radiation of AGN, have sizes of up to several kpc and correlate more strongly with the QSO continuum luminosity than with the integrated [OIII] luminosity. We find a relation of the form log(r)~(0.46+-0.04)log(L_5100), reinforcing the picture of an approximately constant ionization parameter for the ionized clouds across the ENLR. Besides the ENLR, we also find gas ionized by young massive stars in more than 50 per cent of the galaxies on kpc scales. In more than half of the sample, the specific star formation rates based on the extinction-corrected Ha luminosity are consistent with those of inactive disc-dominated galaxies, even for some bulge-dominated QSO hosts. Enhanced SFRs of up to 70Msun/yr are rare and always associated with signatures of major mergers. Comparison with the SFR based on the 60+100micron FIR luminosity suggests that the FIR luminosity is systematically contaminated by AGN emission and Ha appears to be a more robust and sensitive tracer for the star formation rate. Evidence for efficient AGN feedback is scarce in our sample, but some of our QSO hosts lack signatures of ongoing star formation leading to a reduced specific SFR with respect to the main sequence of galaxies. Based on 12 QSOs where we can make measurements, we find that on average bulge-dominated QSO host galaxies tend to fall below the mass-metallicity relation compared to their disc-dominated counterparts. While not yet statistically significant for our small sample, this may provide a useful diagnostic for future large surveys if this metal dilution can be shown to be linked to recent or ongoing galaxy interactions.

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (>10^11Msol) gas-rich spirals NGC4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties and stellar population diagnostics, over the full optical extent of both galaxies. The Mice provide a perfect case study highlighting the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars likely induced in both galaxies. The barred spiral NGC4676B exhibits a strong twist in both its stellar and ionised gas disk. On the other hand, the impact of the merger on the stellar populations has been minimal thus far: star formation induced by the recent close passage has not contributed significantly to the global star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks extend to ~6.6kpc above the disk plane. The measured ram pressure and mass outflow rate (~8-20Msol/yr) are similar to superwinds from local ULIRGs, although NGC4676A has only a moderate infrared luminosity of 3x10^10Lsol. Energy beyond that provided by the mechanical energy of the starburst appears to be required to drive the outflow. We compare the observations to mock kinematic and stellar population maps from a merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

We present 80 stellar and ionised gas velocity maps from the Calar Alto Legacy Integral Field Area (CALIFA) survey in order to characterize the kinematic orientation of non-interacting galaxies. The study of galaxies in isolation is a key step towards understanding how fast-external processes, such as major mergers, affect kinematic properties in galaxies. We derived the global and individual (projected approaching and receding sides) kinematic position angles (PAs) for both the stellar and ionised gas line-of-sight velocity distributions. When compared to the photometric PA, we find that morpho-kinematic differences are smaller than 22 degrees in 90% of the sample for both components; internal kinematic misalignments are generally smaller than 16 degrees. We find a tight relation between the global stellar and ionised gas kinematic PA consistent with circular-flow pattern motions in both components. This relation also holds generally in barred galaxies across the bar and galaxy disk scales. Our findings suggest that even in the presence of strong bars, both the stellar and the gaseous components tend to follow the gravitational potential of the disk. As a result, kinematic orientation can be used to assess the degree of external distortions in interacting galaxies.

We present the kinematic properties of a type-2 QSO, SDSS J132323.33-015941.9 at z~0.35, based on the analysis of Very Large Telescope integral field spectroscopy and Hubble Space Telescope (HST) imaging, which suggest that the target is a binary active galactic nucleus (AGN) with double narrow-line regions. The QSO features double-peaked emission lines ([OIII] and Hb) which can be decomposed into two kinematic components. The flux-weighted centroids of the blue and red components are separated by ~0.2" (0.8 kpc in projection) and coincide with the location of the two stellar cores detected in the HST broad-band images, implying that both stellar cores host an active black hole. The line-of-sight velocity of the blue component is comparable to the luminosity-weighted velocity of stars in the host galaxy while the red component is redshifted by ~240 km/s, consistent with typical velocity offsets of two cores in a late stage of a galaxy merger. If confirmed, the target is one of the rare cases of sub-kpc scale binary AGNs, providing a test-bed for understanding the binary AGN population.

We present the largest and most homogeneous catalog of HII regions and associations compiled so far. The catalog comprises more than 7000 ionized regions, extracted from 306 galaxies observed by the CALIFA survey. We describe the procedures used to detect, select, and analyse the spectroscopic properties of these ionized regions. In the current study we focus on the characterization of the radial gradient of the oxygen abundance in the ionized gas, based on the study of the deprojected distribution of HII regions. We found that all galaxies without clear evidence of an interaction present a common gradient in the oxygen abundance, with a characteristic slope of alpha = -0.1 dex/re between 0.3 and 2 disk effective radii, and a scatter compatible with random fluctuations around this value, when the gradient is normalized to the disk effective radius. The slope is independent of morphology, incidence of bars, absolute magnitude or mass. Only those galaxies with evidence of interactions and/or clear merging systems present a significant shallower gradient, consistent with previous results. The majority of the 94 galaxies with H ii regions detected beyond 2 disk effective radii present a flattening in the oxygen abundance. The flattening is statistically significant. We cannot provide with a conclusive answer regarding the origin of this flattening. However, our results indicate that its origin is most probably related to the secular evolution of galaxies. Finally, we find a drop/truncation of the oxygen abundance in the inner regions for 26 of the galaxies. All of them are non-interacting, mostly unbarred, Sb/Sbc galaxies. This feature is associated with a central star-forming ring, which suggests that both features are produced by radial gas flows induced by resonance processes.

Over the past decade, 3D optical spectroscopy has become the preferred tool for understanding the properties of galaxies and is now increasingly used to carry out galaxy surveys. Low redshift surveys include SAURON, DiskMass, ATLAS3D, PINGS and VENGA. At redshifts above 0.7, surveys such as MASSIV, SINS, GLACE, and IMAGES have targeted the most luminous galaxies to study mainly their kinematic properties. The on-going CALIFA survey ($z\sim0.02$) is the first of a series of upcoming Integral Field Spectroscopy (IFS) surveys with large samples representative of the entire population of galaxies. Others include SAMI and MaNGA at lower redshift and the upcoming KMOS surveys at higher redshift. Given the importance of spatial scales in IFS surveys, the study of the effects of spatial resolution on the recovered parameters becomes important. We explore the capability of the CALIFA survey and a hypothetical higher redshift survey to reproduce the properties of a sample of objects observed with better spatial resolution at lower redshift. Using a sample of PINGS galaxies, we simulate observations at different redshifts. We then study the behaviour of different parameters as the spatial resolution degrades with increasing redshift.

We study the radial structure of the stellar mass surface density ($\mu$) and stellar population age as a function of the total stellar mass and morphology for a sample of 107 galaxies from the CALIFA survey. We use the fossil record to recover the star formation history (SFH) in spheroidal and disk dominated galaxies with masses from 10$^9$ to 10$^{12}$ M$_\odot$. We derive the half mass radius, and we find that galaxies are on average 15% more compact in mass than in light. HMR/HLR decreases with increasing mass for disk galaxies, but is almost constant in spheroidal galaxies. We find that the galaxy-averaged stellar population age, stellar extinction, and $\mu$ are well represented by their values at 1 HLR. Negative radial gradients of the stellar population ages support an inside-out formation. The larger inner age gradients occur in the most massive disk galaxies that have the most prominent bulges; shallower age gradients are obtained in spheroids of similar mass. Disk and spheroidal galaxies show negative $\mu$ gradients that steepen with stellar mass. In spheroidal galaxies $\mu$ saturates at a critical value that is independent of the galaxy mass. Thus, all the massive spheroidal galaxies have similar local $\mu$ at the same radius (in HLR units). The SFH of the regions beyond 1 HLR are well correlated with their local $\mu$, and follow the same relation as the galaxy-averaged age and $\mu$; suggesting that local stellar mass surface density preserves the SFH of disks. The SFH of bulges are, however, more fundamentally related to the total stellar mass, since the radial structure of the stellar age changes with galaxy mass even though all the spheroid dominated galaxies have similar radial structure in $\mu$. Thus, galaxy mass is a more fundamental property in spheroidal systems while the local stellar mass surface density is more important in disks.

Galaxies, which often contain ionised gas, sometimes also exhibit a so-called low-ionisation nuclear emission line region (LINER). For 30 years this was attributed to a central mass-accreting supermassive black hole (AGN) of low luminosity, making LINER galaxies the largest AGN-sub-population, dominating in numbers over higher luminosity Seyfert galaxies and quasars. This, however, poses a serious problem. While the inferred energy balance is plausible, many LINERs clearly do not contain any other independent signatures of an AGN. Using integral field spectroscopic data from the CALIFA survey, we aim at comparing the observed radial surface brightness profiles with what is expected from illumination by an AGN. Essential for this analysis is a proper extraction of emission-lines, especially weak lines such as the Balmer Hb line which is superposed on an absorption trough. To accomplish this, we use the GANDALF code which simultaneously fits the underlying stellar continuum and emission lines. We show for 48 galaxies with LINER-like emission, that the radial emission-line surface brightness profiles are inconsistent with ionisation by a central point-source and hence cannot be due to an AGN alone. The most probable explanation for the excess LINER-like emission is ionisation by evolved stars during the short but very hot and energetic phase known as post-AGB. This leads us to an entirely new interpretation. Post-AGB stars are ubiquitous and their ionising effect should be potentially observable in every galaxy with gas present and stars older than ~1 Gyr, unless a stronger radiation field from young hot stars or an AGN outshines them. This means that galaxies with LINER-like emission are in fact not a class defined by a property, but rather by the absence of a property. It also explains why LINER emission is observed mostly in massive galaxies with old stars and little star formation.

The use of IFS is since recently allowing to measure the emission line fluxes of an increasingly large number of star-forming galaxies both locally and at high redshift. The main goal of this study is to review the most widely used empirical oxygen calibrations, O3N2 and N2, by using new direct abundance measurements. We pay special attention to the expected uncertainty of these calibrations as a function of the index value or abundance derived and the presence of possible systematic offsets. This is possible thanks to the analysis of the most ambitious compilation of Te-based HII regions to date. This new dataset compiles the Te-based abundances of 603 HII regions extracted from the literature but also includes new measurements from the CALIFA survey. Besides providing new and improved empirical calibrations for the gas abundance, we also present here a comparison between our revisited calibrations with a total of 3423 additional CALIFA HII complexes with abundances derived using the ONS calibration by Pilyugin et al. (2010). The combined analysis of Te-based and ONS abundances allows us to derive their most accurate calibration to date for both the O3N2 and N2 single-ratio indicators, in terms of all statistical significance, quality and coverage of the space of parameters. In particular, we infer that these indicators show shallower abundance dependencies and statistically-significant offsets compared to those of Pettini and Pagel (2004), Nagao et al. (2006) and P\'erez-Montero and Contini (2009). The O3N2 and N2 indicators can be empirically applied to derive oxygen abundances calibrations from either direct abundance determinations with random errors of 0.18 and 0.16, respectively, or from indirect ones (but based on a large amount of data) reaching an average precision of 0.08 and 0.09 dex (random) and 0.02 and 0.08 dex (systematic; compared to the direct estimations),respectively.

In a companion paper we have presented many products derived from the application of the spectral synthesis code STARLIGHT to datacubes from the CALIFA survey, including 2D maps of stellar population properties and 1D averages in the temporal and spatial dimensions. Here we evaluate the uncertainties in these products. Uncertainties due to noise and spectral shape calibration errors and to the synthesis method are investigated by means of a suite of simulations based on 1638 CALIFA spectra for NGC 2916, with perturbations amplitudes gauged in terms of the expected errors. A separate study was conducted to assess uncertainties related to the choice of evolutionary synthesis models. We compare results obtained with the Bruzual & Charlot models, a preliminary update of them, and a combination of spectra derived from the Granada and MILES models. About 100k CALIFA spectra are used in this comparison. Noise and shape-related errors at the level expected for CALIFA propagate to 0.10-0.15 dex uncertainties in stellar masses, mean ages and metallicities. Uncertainties in A_V increase from 0.06 mag in the case of random noise to 0.16 mag for shape errors. Higher order products such as SFHs are more uncertain, but still relatively stable. Due to the large number statistics of datacubes, spatial averaging reduces uncertainties while preserving information on the history and structure of stellar populations. Radial profiles of global properties, as well as SFHs averaged over different regions are much more stable than for individual spaxels. Uncertainties related to the choice of base models are larger than those associated with data and method. Differences in mean age, mass and metallicity are ~ 0.15 to 0.25 dex, and 0.1 mag in A_V. Spectral residuals are ~ 1% on average, but with systematic features of up to 4%. The origin of these features is discussed. (Abridged)

We use deep integral field spectroscopy data from the CALIFA survey to study the warm interstellar medium (WIM) of 32 nearby early-type galaxies (ETGs). We propose a tentative subdivision of our sample ETGs into two groups, according to their Ha equivalent width (EW) and Lyman continuum (LyC) photon escape fraction (PLF). Type i ETGs show nearly constant EWs and a PLF~0, suggesting that photoionization by post-AGB stars is the main driver of their faint extranuclear nebular emission. Type ii ETGs are characterized by very low, outwardly increasing EWs, and a PLF as large as ~0.9 in their centers. Such properties point to a low, and inwardly decreasing WIM density and/or volume filling factor. We argue that, because of extensive LyC photon leakage, emission-line luminosities and EWs are reduced in type ii ETG nuclei by at least one order of magnitude. Consequently, the line weakness of these ETGs is by itself no compelling evidence for their containing merely "weak" (sub-Eddington accreting) active galactic nuclei (AGN). In fact, LyC photon escape, which has heretofore not been considered, may constitute a key element in understanding why many ETGs with prominent signatures of AGN activity in radio continuum and/or X-ray wavelengths show only faint emission lines and weak signatures of AGN activity in their optical spectra. The LyC photon escape, in conjunction with dilution of nuclear EWs by line-of-sight integration through a triaxial stellar host, can systematically impede detection of AGN in gas-poor galaxy spheroids through optical emission-line spectroscopy. We further find that type i and ii ETGs differ little (~0.4 dex) in their mean BPT line ratios, which in both cases are characteristic of LINERs. This potentially hints at a degeneracy of the projected, luminosity-weighted BPT ratios for the specific 3D properties of the WIM in ETGs. (abridged)

Fossil record methods based on spectral synthesis techniques have matured over the past decade, and their application to integrated galaxy spectra fostered substantial advances on the understanding of galaxies and their evolution. Yet, because of the lack of spatial resolution, these studies are limited to a global view, providing no information about the internal physics of galaxies. Motivated by the CALIFA survey, which is gathering Integral Field Spectroscopy over the full optical extent of 600 galaxies, we have developed an end-to-end pipeline which: (i) partitions the observed data cube into Voronoi zones in order to, when necessary and taking due account of correlated errors, increase the S/N, (ii) extracts spectra, including propagated errors and bad-pixel flags, (iii) feeds the spectra into the STARLIGHT spectral synthesis code, (iv) packs the results for all galaxy zones into a single file, (v) performs a series of post-processing operations, including zone-to-pixel image reconstruction and unpacking the spectral and stellar population properties into multi-dimensional time, metallicity, and spatial coordinates. This paper provides an illustrated description of this whole pipeline and its products. Using data for the nearby spiral NGC 2916 as a show case, we go through each of the steps involved, presenting ways of visualizing and analyzing this manifold. These include 2D maps of properties such as the v-field, stellar extinction, mean ages and metallicities, mass surface densities, star formation rates on different time scales and normalized in different ways, 1D averages in the temporal and spatial dimensions, projections of the stellar light and mass growth (x,y,t) cubes onto radius-age diagrams, etc. The results illustrate the richness of the combination of IFS data with spectral synthesis, providing a glimpse of what is to come from CALIFA and future surveys. (Abridged)

We present the results on the study of the global and local M-Z relation based on the first data available from the CALIFA survey (150 galaxies). This survey provides integral field spectroscopy of the complete optical extent of each galaxy (up to 2-3 effective radii), with enough resolution to separate individual HII regions and/or aggregations. Nearly $\sim$3000 individual HII regions have been detected. The spectra cover the wavelength range between [OII]3727 and [SII]6731, with a sufficient signal-to-noise to derive the oxygen abundance and star-formation rate associated with each region. In addition, we have computed the integrated and spatially resolved stellar masses (and surface densities), based on SDSS photometric data. We explore the relations between the stellar mass, oxygen abundance and star-formation rate using this dataset. We derive a tight relation between the integrated stellar mass and the gas-phase abundance, with a dispersion smaller than the one already reported in the literature ($\sigma_{\Delta{\rm log(O/H)}}=$0.07 dex). Indeed, this dispersion is only slightly larger than the typical error derived for our oxygen abundances. However, we do not find any secondary relation with the star-formation rate, other than the one induced due to the primary relation of this quantity with the stellar mass. We confirm the result using the $\sim$3000 individual HII regions, for the corresponding local relations. Our results agree with the scenario in which gas recycling in galaxies, both locally and globally, is much faster than other typical timescales, like that of gas accretion by inflow and/or metal loss due to outflows. In essence, late-type/disk dominated galaxies seem to be in a quasi-steady situation, with a behavior similar to the one expected from an instantaneous recycling/closed-box model.

This work investigates the effect of the aperture size on derived galaxy properties for which we have spatially-resolved optical spectra. We focus on some indicators of star formation activity and dust attenuation for spiral galaxies that have been widely used in previous work on galaxy evolution. We have used 104 spiral galaxies from the CALIFA survey for which 2D spectroscopy with complete spatial coverage is available. From the 3D cubes we have derived growth curves of the most conspicuous Balmer emission lines (Halpha, Hbeta) for circular apertures of different radii centered at the galaxy's nucleus after removing the underlying stellar continuum. We find that the Halpha flux (f(Halpha)) growth curve follows a well defined sequence with aperture radius showing low dispersion around the median value. From this analysis, we derive aperture corrections for galaxies in different magnitude and redshift intervals. Once stellar absorption is properly accounted for, the f(Halpha)/f(Hbeta) ratio growth curve shows a smooth decline, pointing towards the absence of differential dust attenuation as a function of radius. Aperture corrections as a function of the radius are provided in the interval [0.3,2.5]R_50. Finally, the Halpha equivalent width (EW(Halpha)) growth curve increases with the size of the aperture and shows a very large dispersion for small apertures. This large dispersion prevents the use of reliable aperture corrections for this quantity. In addition, this result suggests that separating star-forming and quiescent galaxies based on observed EW(Halpha) through small apertures is likely to result in low EW(Halpha) star-forming galaxies begin classified as quiescent.

We present the first public data release of the CALIFA survey. It consists of science-grade optical datacubes for the first 100 of eventually 600 nearby (0.005<z<0.03) galaxies, obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5m telescope at the Calar Alto observatory. The galaxies in DR1 already cover a wide range of properties in color-magnitude space, morphological type, stellar mass, and gas ionization conditions. This offers the potential to tackle a variety of open questions in galaxy evolution using spatially resolved spectroscopy. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the nominal wavelength range 3745-7500A with a spectral resolution of 6.0A (FWHM), and (ii) a medium-resolution V1200 setup covering the nominal wavelength range 3650-4840A with a spectral resolution of 2.3A (FWHM). We present the characteristics and data structure of the CALIFA datasets that should be taken into account for scientific exploitation of the data, in particular the effects of vignetting, bad pixels and spatially correlated noise. The data quality test for all 100 galaxies showed that we reach a median limiting continuum sensitivity of 1.0x10^-18erg/s/cm^2/A/arcsec^2 at 5635A and 2.2x10^-18erg/s/cm^2/A/arcsec^2 at 4500A for the V500 and V1200 setup respectively, which corresponds to limiting r and g band surface brightnesses of 23.6mag/arcsec^2 and 23.4mag/arcsec^2, or an unresolved emission-line flux detection limit of roughly 1x10^-17erg/s/cm^2/arcsec^2 and 0.6x10^-17erg/s/cm^2/arcsec^2, respectively. The median spatial resolution is 3.7", and the absolute spectrophotometric calibration is better than 15% (1sigma). We also describe the available interfaces and tools that allow easy access to this first public CALIFA data.

The Calar Alto Legacy Integral Field Area (CALIFA) is an ongoing 3D spectroscopic survey of 600 nearby galaxies of all kinds. This pioneer survey is providing valuable clues on how galaxies form and evolve. Processed through spectral synthesis techniques, CALIFA datacubes allow us to, for the first time, spatially resolve the star formation history of galaxies spread across the color-magnitude diagram. The richness of this approach is already evident from the results obtained for the first 107 galaxies. Here we show how the different galactic spatial sub-components ("bulge" and "disk") grow their stellar mass over time. We explore the results stacking galaxies in mass bins, finding that, except at the lowest masses, galaxies grow inside-out, and that the growth rate depends on a galaxy's mass. The growth rate of inner and outer regions differ maximally at intermediate masses. We also find a good correlation between the age radial gradient and the stellar mass density, suggesting that the local density is a main driver of galaxy evolution.

The growth of galaxies is one of the key problems in understanding the structure and evolution of the universe and its constituents. Galaxies can grow their stellar mass by accretion of halo or intergalactic gas clouds, or by merging with smaller or similar mass galaxies. The gas available translates into a rate of star formation, which controls the generation of metals in the universe. The spatially resolved history of their stellar mass assembly has not been obtained so far for any given galaxy beyond the Local Group. Here we demonstrate how massive galaxies grow their stellar mass inside-out. We report the results from the analysis of the first 105 galaxies of the largest to date three-dimensional spectroscopic survey of galaxies in the local universe (CALIFA). We apply the fossil record method of stellar population spectral synthesis to recover the spatially and time resolved star formation history of each galaxy. We show, for the first time, that the signal of downsizing is spatially preserved, with both inner and outer regions growing faster for more massive galaxies. Further, we show that the relative growth rate of the spheroidal component, nucleus and inner galaxy, that happened 5-7 Gyr ago, shows a maximum at a critical stellar mass ~10^10 Msun. We also find that galaxies less massive than ~10^10 Msun show a transition to outside-in growth, thus connecting with results from resolved studies of the growth of low mass galaxies.

[Abridged] Detecting cosmic ray hits (cosmics) in fiber-fed IFS data of single exposures is a challenging task, because of the complex signal recorded by IFS instruments. Existing detection algorithms are commonly found to be unreliable in the case of IFS data and the optimal parameter settings are usually unknown a-priori for a given dataset. The CALIFA survey generates hundreds of IFS datasets for which a reliable and robust detection algorithm for cosmics is required as an important part of the fully automatic CALIFA data reduction pipeline. We developed a novel algorithm, PyCosmic, which combines the edge-detection algorithm of L.A.Cosmic with a point-spread function convolution scheme. We generated mock data to compute the efficiency of different algorithms for a wide range of characteristic fibre-fed IFS datasets using the PMAS and VIMOS IFS instruments as representative cases. PyCosmic is the only algorithm that achieves an acceptable detection performance for CALIFA data. We find that PyCosmic is the most robust tool with a detection rate of >~90% and a false detection rate <5% for any of the tested IFS data. It has one less free parameter than the L.A.Cosmic algorithm. Only for strongly undersampled IFS data does L.A.Cosmic exceed the performance of PyCosmic by a few per cent. DCR never reaches the efficiency of the other two algorithms and should only be used if computational speed is a concern. Thus, PyCosmic appears to be the most versatile cosmics detection algorithm for IFS data. It is implemented in the new CALIFA data reduction pipeline as well as in recent versions of the multi-instrument IFS pipeline P3D.

(Abridged) We present a detailed analysis of a large sample of 31 low-redshift, mostly radio-quiet type 1 QSOs observed with integral field spectroscopy to study their extended emission-line regions (EELRs). We focus on the ionisation state of the gas, size and luminosity of extended narrow line regions (ENLRs), which corresponds to those parts of the EELR dominated by ionisation from the QSO, as well as the kinematics of the ionised gas. We detect EELRs around 19 of our 31 QSOs (61%) after deblending the unresolved QSO emission and the extended host galaxy light in the integral field data. We identify 13 EELRs to be entirely ionised by the QSO radiation, 3 EELRs are composed of HII regions and 3 EELRs display signatures of both ionisation mechanisms at different locations. The typical size of the ENLR is 10kpc at a median nuclear [OIII] luminosity of log(L([OIII])/[erg/s])=42.7+-0.15. We show that the ENLR sizes are least a factor of 2 larger than determined with HST, but are consistent with those of recently reported type 2 QSOs at matching [OIII] luminosities. The ENLR of type 1 and type 2 QSOs appear to follow the same size-luminosity relation. Furthermore, we show for the first time that the ENLR size is much better correlated with the QSO continuum luminosity than with the total/nuclear [OIII] luminosity. We show that ENLR luminosity and radio luminosity are correlated, and argue that radio jets even in radio-quiet QSOs are important for shaping the properties of the ENLR. Strikingly, the kinematics of the ionised gas is quiescent and likely gravitationally driven in the majority of cases and we find only 3 objects with radial gas velocities exceeding 400km/s in specific regions of the EELR that can be associate with radio jets. In general, these are significantly lower outflow velocities and detection rates compared to starburst galaxies or radio-loud QSOs.

In this work we analyze the spectroscopic properties of a large number of H ii regions, \sim2600, located in 38 galaxies. The sample of galaxies has been assembled from the face-on spirals in the PINGS survey and a sample described in M\'armol-Queralt\'o (2011, henceforth Paper I). All the galaxies were observed using Integral Field Spectroscopy with a similar setup, covering their optical extension up to \sim2.4 effective radii within a wavelength range from \sim3700 to \sim6900{\AA}. We develop a new automatic procedure to detect H ii regions, based on the contrast of the H{\alpha} intensity maps. Once detected, the procedure provides us with the integrated spectra of each individual segmented region. A well-tested automatic decoupling procedure has been applied to remove the underlying stellar population, deriving the main proper- ties of the strongest emission lines in the considered wavelength range (covering from [O ii] {\lambda}3727 to [S ii] {\lambda}6731). A final catalogue of the spectroscopic properties of these regions has been created for each galaxy. In the current study we focused on the understanding of the average properties of the H ii regions and their radial distributions. We find that the gas-phase oxygen abundance and the H{\alpha} equivalent width present negative and positive gradient, respectively. The distribution of slopes is statistically compatible with a random Gaussian distribution around the mean value, if the radial distances are measured in units of the respective effective radius. No difference in the slope is found for galaxies of different morphologies: barred/non-barred, grand-design/flocculent. Therefore, the effective radius is a universal scale length for gradients in the evolution of galaxies. Other properties have a larger variance across each object.

We demonstrate the existence of a -local- relation between galaxy surface mass density, gas metallicity, and star-formation rate density using spatially-resolved optical spectroscopy of HII regions in the local Universe. One of the projections of this distribution, -the local mass-metallicity relation- extends over three orders of magnitude in galaxy mass density and a factor of eight in gas metallicity. We explain the new relation as the combined effect of the differential radial distributions of mass and metallicity in the discs of galaxies, and a selective star-formation efficiency. We use this local relation to reproduce -with remarkable agreement- the total mass-metallicity relation seen in galaxies, and conclude that the latter is a scale-up integrated effect of a local relation, supporting the inside-out growth and downsizing scenarios of galaxy evolution.

We present here the Calar Alto Legacy Integral Field Area (CALIFA) survey, which has been designed to provide a first step in this direction.We summarize the survey goals and design, including sample selection and observational strategy.We also showcase the data taken during the first observing runs (June/July 2010) and outline the reduction pipeline, quality control schemes and general characteristics of the reduced data. This survey is obtaining spatially resolved spectroscopic information of a diameter selected sample of $\sim600$ galaxies in the Local Universe (0.005< z <0.03). CALIFA has been designed to allow the building of two-dimensional maps of the following quantities: (a) stellar populations: ages and metallicities; (b) ionized gas: distribution, excitation mechanism and chemical abundances; and (c) kinematic properties: both from stellar and ionized gas components. CALIFA uses the PPAK Integral Field Unit (IFU), with a hexagonal field-of-view of $\sim1.3\sq\arcmin'$, with a 100% covering factor by adopting a three-pointing dithering scheme. The optical wavelength range is covered from 3700 to 7000 {\AA}, using two overlapping setups (V500 and V1200), with different resolutions: R\sim850 and R\sim1650, respectively. CALIFA is a legacy survey, intended for the community. The reduced data will be released, once the quality has been guaranteed. The analyzed data fulfill the expectations of the original observing proposal, on the basis of a set of quality checks and exploratory analysis. We conclude from this first look at the data that CALIFA will be an important resource for archaeological studies of galaxies in the Local Universe.

[abridged] The metallicities of AGN are usually well above solar in their NLR, often reaching up to several times solar in their broad-line regions. Low-metallicity AGN are rare objects which have so far always been associated with low-mass galaxies hosting low-mass BHs (M_BH<10^6Msun). In this paper we present IFS data of the low-redshift QSO HE 2158-0107 for which we find strong evidence for sub-solar NLR metallicities associated with a massive BH (M_BH~3x10^8Msun). The QSO is surrounded by a large extended emission-line region reaching out to 30kpc from the QSO in a tail-like geometry. We present optical and near-IR images and investigate the properties of the host galaxy. The SED of the host is rather blue, indicative of a significant young age stellar population formed within the last 1Gyr. A 3sigma upper limit of L_bulge<4.5x10^10Lsun for the H band luminosity and a corresponding stellar mass upper limit of M_bulge<3.4x10^10Msun show that the host is offset from the local BH-bulge relations. This is independently supported by the kinematics of the gas. Although the stellar mass of the host galaxy is lower than expected, it cannot explain the exceptionally low metallicity of the gas. We suggest that the extended emission-line region and the galaxy growth are caused by the infall of nearly pristine gas from the environment of the QSO host. Minor mergers of dwarf galaxies or the theoretically predicted smooth accretion of cold gas are both potential drivers behind that process. Since the metallicity of the gas in the NLR is much lower than expected, we suspect that the external gas has already reached the galaxy centre and may even contribute to the current feeding of the BH. HE 2158-0107 appears to represent a particular phase of substantial BH and galaxy growth that can be observationally linked with the accretion of external material from its environment.

We present VIMOS integral field spectroscopy of the brightest radio-quiet QSO on the southern sky HE 1029-1401 at a redshift of z=0.086. Standard decomposition techniques for broad-band imaging are extended to integral field data in order to deblend the QSO and host emission. We perform a tentative analysis of the stellar continuum finding a young stellar population (<100Myr) or a featureless continuum embedded in an old stellar population (10Gyr) typical for a massive elliptical galaxy. The stellar velocity dispersion of sigma_*=320\pm90 km/s and the estimated black hole mass log(M_BH/M_sun)=8.7\pm0.3 are consistent with the local M_BH-sigma_* relation within the errors. For the first time we map the two-dimensional ionised gas distribution and the gas velocity field around HE 1029-1401. While the stellar host morphology is purely elliptical we find a highly structured distribution of ionised gas out to 16 kpc from the QSO. The gas is highly ionised solely by the QSO radiation and has a significantly lower metallicity than would be expected for the stellar mass of the host, indicating an external origin of the gas most likely due to minor mergers. We find a rotating gas disc around the QSO and a dispersion-dominated non-rotating gas component within the central 3 kpc. At larger distances the velocity field is heavily disturbed, which could be interpreted as another signature of past minor merger events. Alternatively, the arc-like structure seen in the ionised gas might also be indicative of a large-scale expanding bubble, centred on and possibly driven by the active nucleus.

We present a determination of the local (z=0) luminosity function of optically selected type 1 (broad-line) Active Galactic Nuclei. Our primary resource is the Hamburg/ESO Survey (HES), which provides a well-defined sample of more than 300 optically bright AGN with redshifts z<0.3 and blue magnitudes B<17.5. AGN luminosities were estimated in two ways, always taking care to minimise photometric biases due to host galaxy light contamination. Firstly, we measured broad-band B_J (blue) magnitudes of the objects over small apertures of the size of the seeing disk. Secondly, we extracted H alpha and H beta broad emission line luminosities from the spectra which should be entirely free of any starlight contribution. The resulting AGN luminosity function (AGNLF) is consistent with a single power law, also when considering the effects of number density evolution within the narrow redshift range. We compared our AGNLF with the H alpha luminosity function of lower luminosity Seyfert 1 galaxies by Hao et al. (2005) and found a smooth transition between both, with excellent agreement in the overlapping region. From the combination of HES and SDSS samples we constructed a single local AGNLF spanning more than 4 orders of magnitude in luminosity. It shows only mild curvature which can be well described as a double power law with slope indices of -2.0 for the faint end and -2.8 for the bright end. We predicted the local AGNLF in the soft X-ray domain and compared this to recent literature data. We also compared the local AGNLF with results obtained at higher redshifts and find strong evidence for luminosity-dependent evolution, in the sense that AGN with luminosities around M_B~-19 are as common in the local universe as they were at z=1.5, supporting the 'AGN downsizing' picture (abridged).

We use one of the deepest spectroscopic samples of broad line active galactic nuclei (AGN) currently available, extracted from the VIMOS VLT Deep Survey (VVDS), to compute MgII and CIV virial masses estimate of 120 super-massive black holes in the redshift range 1.0<z<1.9 and 2.6<z<4.3. We find that the mass-luminosity relation shows considerably enhanced dispersion towards low AGN luminosities (log L_bol ~ 45). At these luminosities, there is a substantial fraction of black holes accreting far below their Eddington limit (L_bol/L_Edd < 0.1), in marked contrast to what is generally found for AGN of higher luminosities. We speculate that these may be AGN on the decaying branch of their light-curves, well past their peak activity. This would agree with recent theoretical predictions of AGN evolution. In the electronic Appendix of this paper we publish an update of the VVDS type-1 AGN sample, including the first and most of the second epoch observations. This sample contains 298 objects of which 168 are new.

We searched for the presence of extended emission-line regions (EELRs) around low-redshift QSOs. We observed a sample of 20 mainly radio-quiet low-redshift quasars (z<0.3) by means of integral field spectroscopy. After decomposing the extended and nuclear emission components, we constructed [OIII] 5007 narrow-band images of the EELR to measure the total flux. From the same data we obtained high S/N (>50) nuclear spectra to measure properties such as [OIII]/Hbeta flux ratios, FeII equivalent widths and Hbeta line widths. A significant fraction of the quasars (8/20) show a luminous EELR, with detected linear sizes of several kpc. Whether or not a QSO has a luminous EELR is strongly related with nuclear properties, in the sense that an EELR was detected in objects with low FeII equivalent width and large Hbeta FWHM. The EELRs were detected preferentially in QSOs with larger black hole masses. There is no discernible relation, however, between EELR detection and QSO luminosity and Eddington ratio.

We present the analysis and results of a 20 ks XMM-Newton observation of RBS1423. X-ray spectral analysis is used to establish a significantly broadened relativistic iron K-alpha line from a highly ionised disk. A QSO at z=2.262 was considered to be the optical counterpart of this ROSAT Bright Survey X-ray source. Based on the improved XMM-Newton source position we identified a z=0.208 QSO as optical counterpart to RBS1423. The 0.2-12 keV X-ray luminosity of this radio-quiet QSO is 6x10^{44} erg/s. The XMM-EPIC spectra are well described by a power law with a significantly broadened iron K-alpha line. Disk line models for both Schwarzschild and Kerr black holes require hydrogen-like iron ions to fit the measured line profile. Significant ionisation of the reflection disk is confirmed by model fits with ionised disk models, resulting in an ionisation parameter xi~2000.