Source author record

S. F. Sánchez

S. F. Sánchez appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

75works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

75 published item(s)

preprint2025arXiv

Beyond diagnostic-diagrams: A critical exploration on the classification of ionization processes

Optical emission-line diagnostic diagrams, such as the classical BPT, are widely used to identify ionisation mechanisms in galaxies but often suffer from degeneracies, especially when multiple sources coexist. We critically evaluate the effectiveness of these diagnostics in distinguishing star-forming galaxies, retired galaxies (RGs), and active galactic nuclei (AGNs), and propose refined methods to reduce misclassifications. Using a large sample of nearby galaxies from the NASA-Sloan Atlas cross-matched with SDSS spectroscopy, we define representative subsamples of late-type/star-forming galaxies, early-type/RGs, and multiwavelength-selected AGNs. Their distributions are analysed across traditional and modern diagnostics, including WHaN, WHaD, and the newly introduced WHaO diagram, which combine Hα equivalent width with [N II]/Hα, σ(Hα), and [O III]/[O II], respectively. Quantitative comparisons reveal that classical BPT diagrams overestimate star-forming galaxies by ~10% and misclassify up to 45% of AGNs and nearly all RGs. Diagnostics incorporating EW(Hα) achieve improved accuracy, reducing misclassifications to {\sim} 20\ % for AGNs and {\sim} 15\% for RGs. A revised classification scheme based on EW(Hα) thresholds and consistent WHaD/WHaO results yields the highest purity (8-25% misclassifications) and better reflects underlying physical conditions. Our analysis demonstrates that BPT-based methods fail to reliably separate ionisation mechanisms, particularly in galaxies hosting weak AGNs or evolved stellar populations. Updated EW(Hα)-based diagnostics offer a more accurate framework for studying galaxy ionisation and should replace traditional schemes in future population studies.

preprint2022arXiv

HII regions in CALIFA survey: II. The relation between their physical properties and galaxy evolution

We present in here the exploration of the physical properties of the sample of HII regions and aggregations of the last HII regions catalog of the CALIFA survey. This sample comprises the optical spectroscopic properties of more than ~26,000 ionized regions corresponding to 924 galaxies from the Integral Field Spectroscopy data, including the flux intensity and equivalent widths and the properties of their underlying stellar population. In the current study we derive a set of physical quantities for all these regions based on those properties, including (i) the fraction of young stars; (ii) the ionization strength (using six different estimations); (iii) the oxygen abundance (using 25 different calibrators); (iv) the nitrogen and nitrogen-to-oxygen abundance; (v) the dust extinction and (vi) the electron density. Using this dataset we explore how the loci in the classical diagnostic diagrams are connected with those quantities, the radial distributions of these parameters, and the inter-relations between themselves and with the properties of the underlying stellar populations. We conclude that many properties of the HII regions are tightly related to the galactic stellar evolution at the location where those regions are observed. Those properties are modulated only as a second-order effect by the properties of the ionizing stars and the ionized nebulae that do not depend on the astrophysical context in which they are formed. Our results highlight the importance of HII regions to explore the chemical evolution in galaxies, clarifying which of their properties can be used as proxies of that evolution.

preprint2022arXiv

Investigating the link between inner gravitational potential and star-formation quenching in CALIFA galaxies

It has been suggested that the gravitational potential can have a significant role in suppressing the star formation in the nearby galaxies. To establish observational constrains on this scenario, we investigate the connection between the dynamics, through the circular velocity curves (CVCs) as a proxy of the inner gravitational potential, and star formation quenching in 215 non-active galaxies across Hubble sequence from the Calar Alto Legacy Integral Field Area (CALIFA) survey. Our results show that galaxies with similar CVCs tend to have a certain star-formation quenching pattern. To explore these findings in more details, we construct kpc-resolved relations of the equivalent width of the H$α$ ($W_{{\rm H}α}$) versus the amplitude ($V_c$) and shape ($β= d\ln V_c/ d\ln R$) of the circular velocity at given radius. We find that the $W_{{\rm H}α}-V_c$ is a declining relationship, where the retired regions of the galaxies (the ones with $W_{{\rm H}α}$ values below 3 Å) tend to have higher $V_c$. Differently, $W_{{\rm H}α}-β$ is a bi-modal relationship, characterised by two peaks: concentration of the star forming regions at a positive $β$ (rising CVC) and another one of the retired regions with a negative $β$ (declining CVC). Our results show that both the amplitude of the CVC, driven by the mass of the galaxies, and the shape of the CVC, reflecting the internal structure of the galaxies, play an important role in galaxy's quenching history.

preprint2022arXiv

pyFIT3D and pyPipe3D -- The new version of the Integral Field Spectroscopy data analysis pipeline

We present a new version of the FIT3D and Pipe3D codes, two packages to derive properties of the stellar populations and the ionized emission lines from optical spectroscopy and integral field spectroscopy data respectively. The new codes have been fully transcribed to Python from the original Perl and C versions, modifying the algorithms when needed to make use of the unique capabilities of this language with the main goals of (1) respecting as much as possible the original philosophy of the algorithms, (2) maintaining a full compatibility with the original version in terms of the format of the required input and produced output files, and (3) improving the efficiency and accuracy of the algorithms, and solving known (and newly discovered) bugs. The complete package is freely distributed, with an available repository online. pyFIT3D and pyPipe3D are fully tested with data of the most recent IFS data surveys and compilations (e.g. CALIFA, MaNGA, SAMI and AMUSING++), and confronted with simulations. We describe here the code, its new implementation, its accuracy in recovering the parameters based on simulations, and a showcase of its implementation on a particular dataset.

preprint2022arXiv

pyHIIextractor: A tool to detect and extract physical properties of HII regions from Integral Field Spectroscopic data

We present a new code named pyHIIextractor, which detects and extracts the main features (positions and radii) of clumpy ionized regions, i.e. candidate HII regions, using Hα emission line images. Our code is optimized to be used on the dataproducts provided by the Pipe3D pipeline (or dataproducts with such a format), applied to high spatial resolution Integral Field Spectroscopy data (like that provided by the AMUSING++ compilation, using MUSE). The code provides the properties of both the underlying stellar population and the emission lines for each detected H ii candidate. Furthermore, the code delivers a novel estimation of the diffuse ionized gas (DIG) component, independent of its physical properties, which enables a decontamination of the properties of the HII regions from the DIG. Using simulated data, mimicking the expected observations of spiral galaxies, we characterise pyHIIextractor and its ability to extract the main properties of the H ii regions (and the DIG), including the line fluxes, ratios and equivalent widths. Finally, we compare our code with other such tools adopted in the literature, which have been developed or used for similar purposes: pyhIIexplorer, SourceExtractor, HIIphot, and astrodendro. We conclude that pyHIIextractor exceeds the performance of previous tools in aspects such as the number of recovered regions and the distribution of sizes and fluxes (an improvement that is especially noticeable for the faintest and smallest regions). pyHIIextractor is therefore an optimals tool to detect candidate HII regions, offering an accurate estimation of their properties and a good decontamination of the DIG component.

preprint2022arXiv

Revisiting the Mass-Excitation (MEx) diagram using the MaNGA dataset

The diagram comparing the flux ratio of the [OIII] and H$β$ emission lines with the total stellar mass of galaxies (also known as the mass-excitation diagram, MEx) has been widely used to classify the ionization mechanism in high redshift galaxies between star formation and active galactic nuclear ones. This diagram was mainly derived using single-fiber spectroscopy from the SDSS-DR7 survey. In this study, we revise this diagram using the central and integrated spectral measurement from the entire Integral Field Spectroscopic MaNGA sample. Our results suggest that along with the physical parameters of this diagram, the equivalent width of the H$α$ emission line is also required to constrain the ionization mechanism of a high-redshifted galaxy. Furthermore, the location of a galaxy in the excitation-mass diagram varies depending on the use of central or integrated properties.

preprint2022arXiv

SDSS-IV MaNGA: The radial distribution of physical properties within galaxies in the nearby universe

Using the largest sample of galaxies observed with an optical integral field unit (IFU, the SDSS-IV MaNGA survey, $\sim$10000 targets), we derive the radial distribution of the physical properties obtained from the stellar continuum and the ionized-gas emission lines. Given the large sample, we are able to explore the impact of the total stellar mass and morphology by averaging those radial distributions for different bins of both global properties. We use a piece-wise analysis to characterize the slopes of the gradients from those properties at different galactocentric distances. In general we find that most of the properties -- derived from both the stellar continuum and the ionized gas emission lines -- exhibit a negative gradient with a secondary impact by global properties such as the total stellar mass or morphology. Our results confirm the intimate interplay between the properties of the stellar component and those of the ionized gas at local (kpc) scales in order to set the observed gradients. Furthermore, the resemblance of the gradients for similar global properties (in particular for the stellar parameters) indicates statistical similar histories of star formation and chemical enrichment with an initial radial gas distribution following the potential of the galaxy.

preprint2022arXiv

Which galaxy property is the best gauge of the oxygen abundance?

We present an extensive exploration of the impact of 29 physical parameters in the oxygen abundance for a sample of 299 star-forming galaxies extracted from the extended CALIFA sample. We corroborate that the stellar mass is the physical parameter that better traces the observed oxygen abundance (i.e., the mass-metallicity relation, MZR), while other physical parameters could play a potential role in shaping this abundance, but with a lower significant impact. We find that the functional form that best describes the MZR is a third-order polynomial function. From the residuals between this best functional form and the MZR, we find that once considered the impact of the mass in the oxygen abundance, the other physical parameters do not play a significant secondary role in shaping the oxygen abundance in these galaxies (including the gas fraction or the star formation rate). Our analysis suggests that the origin of the MZR is related to the chemical enrichment evolution of the interstellar medium due, most likely, to the build-up of stellar mass in these star-forming galaxies.

preprint2021arXiv

SDSS-IV MaNGA: global and local stellar population properties of elliptical galaxies

We study the spatially resolved properties of 343 elliptical galaxies with the MaNGA/SDSS-IV survey. We used the DESI Legacy Imaging Surveys for morphological classification. Based on integrated spectroscopic properties and colors, we classified seven classes of ellipticals. We inferred the stellar age and metallicity ($Z$) gradients out to a 1.5 effective radius of classical "red and dead", recently quenched, and blue star-forming ellipticals (CLEs, RQEs, and BSFs), corresponding to 73%, 10%, and 4% of the sample, respectively. We reconstructed their global and radial histories of star formation and mass growth. We find the mass- and luminosity-weighted age gradients of CLEs are nearly flat or mildly negative. The respective $Z$ gradients are negative, being flatter as the mass is smaller. The more massive CLEs formed stars earlier and quenched faster than the less massive ones. The CLEs show a weak inside-out growth and a clear inside-out quenching. They finished their quenching globally 3.8 $\pm$ 1.2 Gyr ago, with quenching time-scales of 3.4 $\pm$ 0.8 Gyr. At masses < 10$^{11}$ $M_{\bigodot}$, the age and $Z$ gradients of the RQEs and BSFs are flatter than those of the CLEs but with larger scatters. They show very weak inside-out growth and quenching, which is slow and even not completed at $z\sim0$ for the BSFs. Instead, the massive RQEs show an outside-in quenching and positive gradients in the luminosity-weighted age and stellar metallicities. The RQEs of all masses quenched 1.2 $\pm$ 0.9 Gyr ago on average. Our results for the CLEs are consistent with a two-phase scenario where their inner parts formed by an early and coeval dissipative collapse with a consequent burst of star formation and further quenching, whereas the outer parts continued their assembly likely by dry mergers. We also discuss some evolutionary scenarios for the RQE and BSF galaxies.

preprint2021arXiv

The differences between mass- and light-derived structural parameters over time for MaNGA Elliptical galaxies

We apply stellar population synthesis analysis to obtain spatially-resolved archaeological inferences for a large sample of "red and dead" Elliptical galaxies (Classical Ellipticals; CLEs) from the MaNGA/SDSS-IV DR15 survey. From their 2D stellar light and mass maps, we explore the differences between the radial mass and light distributions in the rest-frame bands $g,$ $r,$ and $i$ as functions of look-back time, $t_{\rm lb}$, or redshift, $z$. We characterize these differences through the ratios between the following mass- and light-derived global properties: sizes, concentrations, and effective surface densities. We find that the mass-to-light ratios of these properties change with $t_{\rm lb}$, more the more massive the galaxies are. The CLE galaxy archaeological progenitors are, on average, less compact, concentrated, and dense in light than in mass as $z$ decreases. However, at later times, when also the evolution of the progenitors becomes passive at all radii, there is an upturn in these trends and the differences between mass and light in compactness/concentration decrease towards $z\sim 0$. The trends in the ratios of mass to light sizes agree qualitatively with results from direct observations in galaxy surveys at different redshifts. We discuss the caveats and interpretations of our results, and speculate that the strong structural evolution found in some previous studies for early-type galaxies could be explained partially by photometric changes rather than by intrinsic structural changes.

preprint2020arXiv

A Universal fundamental plane and the $M_{dyn}-M_{\star}$ relation for galaxies with CALIFA and MaNGA

We use the stellar kinematics for $2458$ galaxies from the MaNGA survey to explore dynamical scaling relations between the stellar mass $M_{\star}$ and the total velocity parameter at the effective radius, $R_e$, defined as $S_{K}^{2}=KV_{R_e}^{2}+σ_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $σ_{\star_e}$. We confirm that spheroidal and spiral galaxies follow the same $M_{\star}-S_{0.5}$ scaling relation with lower scatter than the $M_{\star}-V_{R_e}$ and $M_{\star}-σ_{\star_e}$ ones. We also explore a more general Universal Fundamental Plane described by the equation $log(Υ_{e}) = log (S_{0.5}^{2}) - log (I_{e}) - log (R_{e}) + C$, which in addition to kinematics, $S_{0.5}$, and effective radius, $R_e$, it includes surface brightness, $I_e$, and dynamical mass-to-light ratio, $Υ_e$. We use sophisticated Schwarzschild dynamical models for a sub-sample of 300 galaxies from the CALIFA survey to calibrate the so called Universal Fundamental Plane. That calibration allows us to propose both: (i) a parametrization to estimate the difficult-to-measure dynamical mass-to-light ratio at the effective radius; and (ii) a new dynamical mass proxy consistent with dynamical models within $0.09\ dex$. We reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. We use the estimated dynamical mass-to-light ratio from our analysis, $Υ_{e}^{fit}$, to explore the Universal Fundamental Plane with the MaNGA data set. We find that all classes of galaxies, from spheroids to disks, follow this Universal Fundamental Plane with a scatter significantly smaller $(0.05\ dex)$ than the one reported for the $M_{\star}-S_{0.5}$ relation $(0.1\ dex)$, the Fundamental Plane $(\sim 0.09\ dex)$ and comparable with Tully-Fisher studies $(\sim 0.05\ dex)$, but for a wider range of galaxy types.

preprint2020arXiv

HII regions in the CALIFA survey: I. Catalog presentation

We present a new catalog of HII regions based on the integral field spectroscopy (IFS) data of the extended CALIFA and PISCO samples. The selection of HII regions was based on two assumptions: a clumpy structure with high contrast of H$α$ emission and an underlying stellar population comprising young stars. The catalog provides the spectroscopic information of 26,408 individual regions corresponding to 924 galaxies, including the flux intensities and equivalent widths of 51 emission lines covering the wavelength range between 3745-7200A. To our knowledge, this is the largest catalog of spectroscopic properties of HII regions. We explore a new approach to decontaminate the emission lines from diffuse ionized gas contribution. This diffuse gas correction was estimated to correct every emission line within the considered spectral range. With the catalog of HII regions corrected, new demarcation lines are proposed for the classical diagnostic diagrams. Finally, we study the properties of the underlying stellar populations of the HII regions. It was found that there is a direct relationship between the ionization conditions on the nebulae and the properties of stellar populations besides the physicals condition on the ionized regions.

preprint2020arXiv

Larger $λ_R$ in the disc of isolated active spiral galaxies than in their non-active twins

We present a comparison of the spin parameter $λ_R$, measured in a region dominated by the galaxy disc, between 20 pairs of nearby (0.005$<$z$<$0.03) seemingly isolated twin galaxies differing in nuclear activity. We find that 80--82% of the active galaxies show higher values of $λ_R$ than their corresponding non-active twin(s), indicating larger rotational support in the AGN discs. This result is driven by the 11 pairs of unbarred galaxies, for which 100% of the AGN show larger $λ_R$ than their twins. These results can be explained by a more efficient angular momentum transfer from the inflowing gas to the disc baryonic matter in the case of the active galaxies. This gas inflow could have been induced by disc or bar instabilities, although we cannot rule out minor mergers if these are prevalent in our active galaxies. This result represents the first evidence of galaxy-scale differences between the dynamics of active and non-active isolated spiral galaxies of intermediate stellar masses (10$^{10}<M_*<10^{11}$ M$_{\odot}$) in the Local Universe.

preprint2020arXiv

MUSE reveals extended circumnuclear outflows in the Seyfert 1 NGC 7469

NGC 7469 is a well known Luminous IR Galaxy, with a circumnuclear star formation ring ($\sim 830$ pc radius) surrounding a Seyfert 1 AGN. Nuclear unresolved winds were previously detected in X-rays and UV, as well as an extended biconical outflow in IR coronal lines. We search for extended outflows by measuring the kinematics of the $\mathrm{Hβ}$ and [O III] $λ5007$ optical emission lines, in data of the VLT/MUSE integral field spectrograph. We find evidence of two outflow kinematic regimes: one slower regime extending across most of the star formation ring -- possibly driven by the massive star formation -- and a faster regime (with a maximum velocity of $-715 \ \mathrm{km \ s^{-1}}$), only observed in [O III], in the western region between the AGN and the massive star forming regions of the ring, likely AGN-driven. This work shows a case where combined AGN/star-formation feedback can be effectively spatially-resolved, opening up a promising path toward a deeper understanding of feedback processes in the central kiloparsec of AGN.

preprint2019arXiv

SDSS-IV MaNGA: Bar pattern speed estimates with the Tremaine-Weinberg method and their error sources

Estimating the bar pattern speed (\Om{}) is one of the main challenges faced in understanding the role of stellar bars in galaxy dynamical evolution. This work aims to characterise different uncertainty sources affecting the Tremaine Weinberg (TW)-method to study the correlation between bar and galaxies physical parameters. We use a sample of 15 MaNGA SDSS-IV galaxies and 3 CALIFA galaxies from \cite{Aguerri2015}. We studied the errors related with (i) galaxy centre determination, (ii) disc position angle (PA) emphasising the difficulties triggered by outer non-axisymmetric structures besides the bar, (iii) the slits length and (iv) the spatial resolution. In average, the PA uncertainties range $\sim 15 \%$, the slit length $\sim 9 \%$ and the centring error $\sim 5 \%$. Reducing the spatial resolution increases the sensitivity to the PA error. Through Monte Carlo simulations, we estimate the probability distribution of the \R{} bar speed parameter. The present sample is composed of 7 slow, 4 fast and 7 ultrafast bars, with no trend with morphological types. Although uncertainties and low sample numbers may mask potential correlations between physical properties, we present a discussion of them: We observe an anti-correlation of \Om{} with the bar length and the stellar mass, suggesting that massive galaxies tend to host longer and slower bars. We also observe a correlation of the molecular gas fraction with \R{}, and a weak anti-correlation with \Om{}, suggesting that bars rotate slower in gaseous discs. Confirmation of such trends awaits future studies.

preprint2016arXiv

Aperture effects on the oxygen abundance determinations from CALIFA data

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.

preprint2016arXiv

CALIFA, the Calar Alto Legacy Integral Field Area survey: IV. Third Public data release

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.

preprint2016arXiv

Evidence of ongoing radial migration in NGC 6754: Azimuthal variations of the gas properties

Understanding the nature of spiral structure in disk galaxies is one of the main, and still unsolved questions in galactic astronomy. However, theoretical works are proposing new testable predictions whose detection is becoming feasible with recent development in instrumentation. In particular, streaming motions along spiral arms are expected to induce azimuthal variations in the chemical composition of a galaxy at a given galactic radius. In this letter we analyse the gas content in NGC 6754 with VLT/MUSE data to characterise its 2D chemical composition and H$α$ line-of-sight velocity distribution. We find that the trailing (leading) edge of the NGC 6754 spiral arms show signatures of tangentially-slower, radially-outward (tangentially-faster, radially-inward) streaming motions of metal-rich (poor) gas over a large range of radii. These results show direct evidence of gas radial migration for the first time. We compare our results with the gas behaviour in a $N$-body disk simulation showing spiral morphological features rotating with a similar speed as the gas at every radius, in good agreement with the observed trend. This indicates that the spiral arm features in NGC 6754 may be transient and rotate similarly as the gas does at a large range of radii.

preprint2016arXiv

First survey of Wolf-Rayet star populations over the full extension of nearby galaxies observed with CALIFA

The search of extragalactic regions with conspicuous presence of Wolf-Rayet (WR) stars outside the Local Group is challenging task due to the difficulties in detecting their faint spectral features. In this exploratory work, we develop a methodology to perform an automated search of WR signatures through a pixel-by-pixel analysis of integral field spectroscopy (IFS) data belonging to the Calar Alto Legacy Integral Field Area survey, CALIFA. This technique allowed us to build the first catalogue of Wolf-Rayet rich regions with spatially-resolved information, allowing to study the properties of these complexes in a 2D context. The detection technique is based on the identification of the blue WR bump (around He II 4686 Å, mainly associated to nitrogen-rich WR stars, WN) and the red WR bump (around C IV 5808 Å and associated to carbon-rich WR stars, WC) using a pixel-by-pixel analysis. We identified 44 WR-rich regions with blue bumps distributed in 25 galaxies of a total of 558. The red WR bump was identified only in 5 of those regions. We found that the majority of the galaxies hosting WR populations in our sample are involved in some kind of interaction process. Half of the host galaxies share some properties with gamma-ray burst (GRB) hosts where WR stars, as potential candidates to being the progenitors of GRBs, are found. We also compared the WR properties derived from the CALIFA data with stellar population synthesis models, and confirm that simple star models are generally not able to reproduce the observations. We conclude that other effects, such as the binary star channel (which could extend the WR phase up to 10 Myr), fast rotation or other physical processes that causes the loss of observed Lyman continuum photons, are very likely affecting the derived WR properties, and hence should be considered when modelling the evolution of massive stars.

preprint2016arXiv

IMF shape constraints from stellar populations and dynamics from CALIFA

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, $Υ_\mathrm{dyn}$ and $Υ_{\ast}$, over a homogenous aperture of 0.5~$R_{e}$. We use the constraint $Υ_\mathrm{dyn} \ge Υ_{\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$σ$) 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.

preprint2016arXiv

MUSE Reveals a Recent Merger in the Post-starburst Host Galaxy of the TDE ASASSN-14li

We present MUSE integral field spectroscopic observations of the host galaxy (PGC 043234) of one of the closest ($z=0.0206$, $D\simeq 90$ Mpc) and best-studied tidal disruption events (TDE), ASASSN-14li. The MUSE integral field data reveal asymmetric and filamentary structures that extend up to $\gtrsim 10$ kpc from the post-starburst host galaxy of ASASSN-14li. The structures are traced only through the strong nebular [O III] $λ$5007, [N II] $λ$6584, and H$α$ emission lines. The total off nuclear [O III] $λ$5007 luminosity is luminosity is $4.7\times 10^{39}$ erg s$^{-1}$ and the ionized H mass is $\rm \sim 10^4(500/n_e)\,M_{\odot}$. Based on the BPT diagram, the nebular emission can be driven by either AGN photoionization or shock excitation, with AGN photoionization favored given the narrow intrinsic line widths. The emission line ratios and spatial distribution strongly resemble ionization nebulae around fading AGNs such as IC 2497 (Hanny's Voorwerp) and ionization "cones" around Seyfert 2 nuclei. The morphology of the emission line filaments strongly suggest that PGC 043234 is a recent merger, which likely triggered a strong starburst and AGN activity leading to the post-starburst spectral signatures and the extended nebular emission line features we see today. We briefly discuss the implications of these observations in the context of the strongly enhanced TDE rates observed in post-starburst galaxies and their connection to enhanced theoretical TDE rates produced by supermassive black-hole binaries.

preprint2016arXiv

Photoionization models of the CALIFA HII regions. I. Hybrid models

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.

preprint2016arXiv

Pipe3D, a pipeline to analyse integral field spectroscopy data: II. Analysis sequence and CALIFA dataproducts

We present Pipe3D, an analysis pipeline based on the FIT3D fitting tool, devel- oped to explore the properties of the stellar populations and ionized gas of Integral Field Spectroscopy data. Pipe3D was created to provide with coherent, simple to distribute, and comparable dataproducts, independently of the origin of the data, focused on the data of the most recent IFU surveys (e.g., CALIFA, MaNGA, and SAMI), and the last generation IFS instruments (e.g., MUSE). Along this article we describe the different steps involved in the analysis of the data, illustrating them by showing the dataproducts derived for NGC 2916, observed by CALIFA and P-MaNGA. As a practical use of the pipeline we present the complete set of dataproducts derived for the 200 datacubes that comprises the V500 setup of the CALIFA Data Release 2 (DR2), making them freely available through the network (ftp://ftp.caha.es/CALIFA/dataproducts/DR2/Pipe3D). Finally, we explore the hypothesis that the properties of the stellar populations and ionized gas of galaxies at the effective radius are representative of the overall average ones, finding that this is indeed the case.

preprint2016arXiv

Shape of the oxygen abundance profiles in CALIFA face-on spiral galaxies

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 $α_{O/H} = -\,0.075\,\rm{dex}/r_e$ with a scatter of $σ= 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.

preprint2016arXiv

Simultaneous spectroscopic and photometric analysis of galaxies with STARLIGHT: CALIFA $+$ GALEX

We present an extended version of the spectral synthesis code STARLIGHT designed to incorporate both $λ$-by-$λ$ spectra and photometric fluxes in the estimation of stellar population properties of galaxies. The code is tested with simulations and data for 260 galaxies culled from the CALIFA survey, spatially matching the 3700--7000 Å optical datacubes to GALEX near and far UV images. The sample spans E--Sd galaxies with masses from $10^9$ to $10^{12} M_\odot$ and stellar populations all the way from star-forming to old, passive systems. Comparing results derived from purely optical fits with those which also consider the NUV and FUV data we find that: (1) The new code is capable of matching the input UV data within the errors while keeping the quality of the optical fit essentially unchanged. (2) Despite being unreliable predictors of the UV fluxes, purely optical fits yield stellar population properties which agree well with those obtained in optical+UV fits for nearly 90% of our sample. (3) The addition of UV constraints has little impact on properties such as stellar mass and dust optical depth. Mean stellar ages and metallicities also remain nearly the same for most galaxies, the exception being low-mass, late-type galaxies, which become older and less enriched due to rearrangements of their youngest populations. (4) The revised ages are better correlated with observables such as the 4000 Å break index, and the $NUV - r$ and $u - r$ colours, an empirical indication that the addition of UV constraints helps mitigating the effects of age-metallicity-extinction degeneracies.

preprint2016arXiv

Spatially-Resolved Star Formation Main Sequence of Galaxies in the CALIFA Survey

The "main sequence of galaxies" $-$ defined in terms of the total star formation rate $ψ$ 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 ($σ= 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.

preprint2016arXiv

Star formation along the Hubble sequence: Radial structure of the star formation of CALIFA galaxies

The aim of this paper is to characterize the radial structure of the star formation rate (SFR) in galaxies in the nearby Universe as represented by the CALIFA survey. The sample under study contains 416 galaxies observed with IFS, covering a wide range of Hubble types and stellar masses. Spectral synthesis techniques are applied to obtain radial profiles of the intensity of the star formation rate in the recent past, and the local sSFR. To emphasize the behavior of these properties for galaxies that are on and off the main sequence of star formation (MSSF) we stack the individual radial profiles in bins of galaxy morphology and stellar masses. Our main results are: a) The intensity of SFR shows declining profiles that exhibit very little differences between spirals. The dispersion between the profiles is significantly smaller in late type spirals. This confirms that the MSSF is a sequence of galaxies with nearly constant intensity of SFR b) sSFR values scale with Hubble type and increase radially outwards, with a steeper slope in the inner 1 HLR. This behavior suggests that galaxies are quenched inside-out, and that this process is faster in the central, bulge-dominated part than in the disks. c) As a whole, and at all radii, E and S0 are off the MSSF. d) Applying the volume-corrections for the CALIFA sample, we obtain a density of star formation in the local Universe of 0.0105 Msun/yr/Mpc^{-3}. Most of the star formation is occurring in the disks of spirals. e) The volume averaged birthrate parameter, b'=0.39, suggests that the present day Universe is forming stars at 1/3 of its past average rate. E, S0, and the bulge of early type spirals contribute little to the recent SFR of the Universe, which is dominated by the disks of later spirals. f) There is a tight relation between the intensity of the SFR and stellar mass, defining a local MSSF relation with a logarithmic slope of 0.8.

preprint2016arXiv

Stellar kinematics across the Hubble sequence in the CALIFA survey: General properties and aperture corrections

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.

preprint2016arXiv

The Califa and Hipass velocity function for all morphological galaxy types

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.

preprint2016arXiv

The dependence of oxygen and nitrogen abundances on stellar mass from the CALIFA survey

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érez-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 $ρ_s$ = 0.80) than between O/H and N/O slopes ($ρ_s$ = 0.39) or for O/H and N/O in the individual \hii\ regions ($ρ_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.

preprint2015arXiv

Bar pattern speeds in CALIFA galaxies: I. Fast bars across the Hubble sequence

The bar pattern speed ($Ω_{\rm b}$) is defined as the rotational frequency of the bar, and it determines the bar dynamics. Several methods have been proposed for measuring $Ω_{\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 $Ω_{\rm b}$ on galaxy properties, such as the Hubble type. We measured $Ω_{\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 $Ω_{\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.

preprint2015arXiv

Central star formation and metallicity in CALIFA interacting galaxies

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.

preprint2015arXiv

IMF - metallicity: a tight local relation revealed by the CALIFA survey

Variations in the stellar initial mass function (IMF) have been invoked to explain the spectroscopic and dynamical properties of early-type galaxies. However, no observations have yet been able to disentangle the physical driver. We analyse here a sample of 24 early-type galaxies drawn from the CALIFA survey, deriving in a homogeneous way their stellar population and kinematic properties. We find that the local IMF is tightly related to the local metallicity, becoming more bottom-heavy towards metal-rich populations. Our result, combined with the galaxy mass-metallicity relation, naturally explains previous claims of a galaxy mass-IMF relation, derived from non-IFU spectra. If we assume that - within the star formation environment of early-type galaxies - metallicity is the main driver of IMF variations, a significant revision of the interpretation of galaxy evolution observables is necessary.

preprint2015arXiv

No direct coupling between bending of galaxy disc stellar age and light profiles

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.

preprint2015arXiv

On the morphology of dust lanes in galactic bars

The aim of our study is to use dynamical simulations to explore the influence of two important dynamical bar parameters, bar strength and bar pattern speed, on the shape of the bar dust lanes. To quantify the shape of the dust lanes we have developed a new systematic method to measure the dust lane curvature. Previous numerical simulations have compared the curvature of bar dust lanes with the bar strength, predicting a relation between both parameters which has been supported by observational studies but with a large spread. We take into account the bar pattern speed to explore, simultaneously, the effect of both parameters on the dust lane shape. To that end, we separate our galactic bars in fast bars $\left(1 < \mathcal{R} < 1.4 \right)$ and slow bars $\left(\mathcal{R} > 1.4 \right)$, obtaining, as previous simulations, an inverse relation between the dust lane curvature and the bar strength for fast bars. For the first time, we extend the study to slow bars, finding a constant curvature as a function of the bar strength. As a result, we conclude that weak bars with straight dust lanes are candidates for slow bars. Finally, we have analysed a pilot sample of ten S$^4$G galaxies, obtaining dust lane curvatures lying within the range covered by the simulations.

preprint2015arXiv

Outer-disk reddening and gas-phase metallicities: The CALIFA connection

We study, for the first time in a statistically significant and well-defined sample, the relation between the outer-disk ionized-gas metallicity gradients and the presence of breaks in the surface brightness profiles of disk galaxies. SDSS g'- and r'-band surface brightness, (g'- r') color, and ionized-gas oxygen abundance profiles for 324 galaxies within the CALIFA survey are used for this purpose. We perform a detailed light-profile classification finding that 84% of our disks show down- or up-bending profiles (Type II and Type III, respectively) while the remaining 16% are well fitted by one single exponential (Type I). The analysis of the color gradients at both sides of this break shows a U-shaped profile for most Type II galaxies with an average minimum (g'- r') color of ~0.5 mag and a ionized-gas metallicity flattening associated to it only in the case of low-mass galaxies. More massive systems show a rather uniform negative metallicity gradient. The correlation between metallicity flattening and stellar mass results in p-values as low as 0.01. Independently of the mechanism having shaped the outer light profiles of these galaxies, stellar migration or a previous episode of star formation in a shrinking star-forming disk, it is clear that the imprint in their ionized-gas metallicity was different for low- and high-mass Type II galaxies. In the case of Type III disks, a positive correlation between the change in color and abundance gradient is found (the null hypothesis is ruled out with a p-value of 0.02), with the outer disks of Type III galaxies with masses $\leq$10$^{10}$ M$_{\odot}$ showing a weak color reddening or even a bluing. This is interpreted as primarily due to a mass down-sizing effect on the population of Type III galaxies having recently experienced an enhanced inside-out growth.

preprint2015arXiv

P-MaNGA: Emission Lines Properties - Gas Ionisation and Chemical Abundances from Prototype Observations

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a SDSS-IV survey that will obtain spatially resolved spectroscopy from 3600 Å to 10300 Å for a representative sample of over 10000 nearby galaxies. In this paper we present the analysis of nebular emission line properties in 14 galaxies obtained with P-MaNGA, a prototype of the MaNGA instrument. Using spatially resolved diagnostic diagrams we find extended star formation in galaxies that are centrally dominated by Seyfert/LINER-like emission, which illustrates that galaxy characterisations based on single fibre spectra are necessarily incomplete. We observe extended LINER-like emission (up to $\rm 1 R_{e}$) in three galaxies. We make use of the $\rm EW(H α)$ to argue that the observed emission is consistent with ionisation from hot evolved stars. We derive stellar population indices and demonstrate a clear correlation between $\rm D_n(4000)$ and $\rm EW(H δ_A)$ and the position in the ionisation diagnostic diagram: resolved galactic regions which are ionised by a Seyfert/LINER-like radiation field are also devoid of recent star formation and host older and/or more metal rich stellar populations. We also detect extraplanar LINER-like emission in two highly inclined galaxies, and identify it with diffuse ionised gas. We investigate spatially resolved metallicities and find a positive correlation between metallicity and star formation rate (SFR) surface density. We study the relation between N/O vs O/H on resolved scales. We find that, at given N/O, regions within individual galaxies are spread towards lower metallicities, deviating from the sequence defined by galactic central regions as traced by Sloan $3''$ fibre spectra. We suggest that the observed dispersion can be a tracer for gas flows in galaxies: infalls of pristine gas and/or the effect of a galactic fountain.

preprint2015arXiv

Pipe3D, a pipeline to analyze Integral Field Spectroscopy data: I. New fitting phylosophy of FIT3D

We present an improved version of FIT3D, a fitting tool for the analysis of the spectroscopic properties of the stellar populations and the ionized gas derived from moderate resolution spectra of galaxies. FIT3D is a tool developed to analyze Integral Field Spectroscopy data and it is the basis of Pipe3D, a pipeline already used in the analysis of datasets like CALIFA, MaNGA, and SAMI. We describe the philosophy behind the fitting procedure, and in detail each of the different steps in the analysis. We present an extensive set of simulations in order to estimate the precision and accuracy of the derived parameters for the stellar populations. In summary, we find that using different stellar population templates we reproduce the mean properties of the stellar population (age, metallicity, and dust attenuation) within ~0.1 dex. A similar approach is adopted for the ionized gas, where a set of simulated emission- line systems was created. Finally, we compare the results of the analysis using FIT3D with those provided by other widely used packages for the analysis of the stellar population (Starlight, Steckmap, and analysis based on stellar indices) using real high S/N data. In general we find that the parameters for the stellar populations derived by FIT3D are fully compatible with those derived using these other tools.

preprint2015arXiv

Spectroscopic aperture biases in inside-out evolving early-type galaxies from CALIFA

Integral field spectroscopy studies based on CALIFA data have recently revealed the presence of ongoing low-level star formation (SF) in the periphery of ~10% of local early-type galaxies (ETGs), witnessing a still ongoing inside-out galaxy growth process. A distinctive property of the nebular component in these ETGs, classified i+, is a two-radial-zone structure, with the inner zone displaying LINER emission with a Hαequivalent width EW~1Å, and the outer one (3Å<EW<~20Å) showing HII-region characteristics. Using CALIFA IFS data, we empirically demonstrate that the confinement of nebular emission to the galaxy periphery leads to a strong aperture (or, redshift) bias in spectroscopic single-fiber studies of type i+ ETGs: At low redshift (<~0.45), SDSS spectroscopy is restricted to the inner (SF-devoid LINER) zone, thereby leading to their erroneous classification as "retired" galaxies (systems lacking SF and whose faint emission is powered by pAGB stars). Only at higher z's the SDSS aperture can encompass the outer SF zone, permitting their unbiased classification as "composite SF/LINER". We also demonstrate that the principal effect of a decreasing aperture on the classification of i+ ETGs via standard BPT emission-line ratios consists in a monotonic up-right shift precisely along the upper-right wing of the "seagull" distribution. Motivated by these insights, we also investigate theoretically these biases in aperture-limited studies of inside-out growing galaxies as a function of z. To this end, we devise a simple model, which involves an outwardly propagating SF process, that reproduces the radial extent and two-zone EW distribution of i+ ETGs. By simulating on this model the spectroscopic SDSS aperture, we find that SDSS studies at z<~1 are progressively restricted to the inner LINER-zone, and miss an increasingly large portion of the Hα-emitting periphery.

preprint2015arXiv

Spiral-like star-forming patterns in CALIFA early-type galaxies

Based on a combined analysis of SDSS imaging and CALIFA integral field spectroscopy data, we report on the detection of faint (24 < μ$_r$ mag/arcsec$^2$ < 26) star-forming spiral-arm-like features in the periphery of three nearby early-type galaxies (ETGs). These features are of considerable interest because they document the still ongoing inside-out growth of some local ETGs and may add valuable observational insight into the origin and evolution of spiral structure in triaxial stellar systems. A characteristic property of the nebular component in the studied ETGs, classified i+, is a two-radial-zone structure, with the inner zone that displays faint (EW(Hα)$\simeq$1Å) low-ionization nuclear emission-line region (LINER) properties, and the outer one (3Å<EW(Hα)<~20Å) HII-region characteristics. This spatial segregation of nebular emission in two physically distinct concentric zones calls for an examination of aperture effects in studies of type i+ ETGs with single-fiber spectroscopic data.

preprint2015arXiv

Star Formation in the Local Universe from the CALIFA sample. I. Calibrating the SFR using IFS data

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α, 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α. The extinction-corrected Hα luminosity agrees with the updated hybrid SFR estimators based on either UV or Hα 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$α$ 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α 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.

preprint2015arXiv

Testing the Modern Merger Hypothesis via the Assembly of Massive Blue Elliptical Galaxies in the Local Universe

The modern merger hypothesis offers a method of forming a new elliptical galaxy through merging two equal-mass, gas-rich disk galaxies fuelling a nuclear starburst followed by efficient quenching and dynamical stabilization. A key prediction of this scenario is a central concentration of young stars during the brief phase of morphological transformation from highly-disturbed remnant to new elliptical galaxy. To test this aspect of the merger hypothesis, we use integral field spectroscopy to track the stellar Balmer absorption and 4000Å break strength indices as a function of galactic radius for 12 massive (${\rm M_{*}}\ge10^{10}{\rm M_{\odot}}$), nearby (${\rm z}\le0.03$), visually-selected plausible new ellipticals with blue-cloud optical colours and varying degrees of morphological peculiarities. We find that these index values and their radial dependence correlate with specific morphological features such that the most disturbed galaxies have the smallest 4000Å break strengths and the largest Balmer absorption values. Overall, two-thirds of our sample are inconsistent with the predictions of the modern merger hypothesis. Of these eight, half exhibit signatures consistent with recent minor merger interactions. The other half have star formation histories similar to local, quiescent early-type galaxies. Of the remaining four galaxies, three have the strong morphological disturbances and star-forming optical colours consistent with being remnants of recent, gas-rich major mergers, but exhibit a weak, central burst consistent with forming $\sim5\%$ of their stars. The final galaxy possesses spectroscopic signatures of a strong, centrally-concentrated starburst and quiescent core optical colours indicative of recent quenching (i.e., a post-starburst signature) as prescribed by the modern merger hypothesis.

preprint2015arXiv

The CALIFA survey across the Hubble sequence: Spatially resolved stellar population properties in galaxies

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.

preprint2015arXiv

The incidence of bar-like kinematic flows in CALIFA galaxies

We carry out a direct search for bar-like non-circular flows in intermediate-inclination, gas-rich disk galaxies with a range of morphological types and photometric bar classifications from the first data release (DR1) of the CALIFA survey. We use the DiskFit algorithm to apply rotation only and bisymmetric flow models to H$α$ velocity fields for 49/100 CALIFA DR1 systems that meet our selection criteria. We find satisfactory fits for a final sample of 37 systems. DiskFit is sensitive to the radial or tangential components of a bar-like flow with amplitudes greater than $15\,$km$\,$s$^{-1}$ across at least two independent radial bins in the fit, or ~2.25 kpc at the characteristic final sample distance of ~75 Mpc. The velocity fields of 25/37 $(67.6^{+6.6}_{-8.5}\%)$ galaxies are best characterized by pure rotation, although only 17/25 $(68.0^{+7.7}_{-10.4}\%)$ of them have sufficient H$α$ emission near the galaxy centre to afford a search for non-circular flows. We detect non-circular flows in the remaining 12/37 $(32.4^{+8.5}_{-6.6}\%)$ galaxies. We conclude that the non-circular flows detected in 11/12 $(91.7^{+2.8}_{-14.9}\%)$ systems stem from bars. Galaxies with intermediate (AB) bars are largely undetected, and our detection thresholds therefore represent upper limits to the amplitude of the non-circular flows therein. We find 2/23 $(8.7^{+9.6}_{-2.9}\%)$ galaxies that show non-circular motions consistent with a bar-like flow, yet no photometric bar is evident. This suggests that in ~10% of galaxies either the existence of a bar may be missed completely in photometry or other processes may drive bar-like flows and thus secular galaxy evolution.

preprint2015arXiv

The warm ionized gas in CALIFA early-type galaxies: 2D emission-line patterns and kinematics for 32 galaxies

The morphological, spectroscopic and kinematical properties of the warm interstellar medium (wim) in early-type galaxies (ETGs) hold key observational constraints to nuclear activity and the buildup history of these massive, quiescent systems. High-quality integral field spectroscopy (IFS) data with a wide spectral and spatial coverage, such as those from the CALIFA survey, offer an unprecedented opportunity for advancing our understanding of the wim in ETGs. This article centers on a 2D investigation of the wim component in 32 nearby (<~150Mpc) ETGs from CALIFA, complementing a previous 1D analysis of the same sample (Papaderos et al. 2013; P13). We include here Hα intensity and equivalent width (EW) maps and radial profiles, diagnostic emission-line ratios, besides ionized-gas and stellar kinematics. This study is supplemented by τ-ratio maps as an efficient means to quantify the role of photoionization by pAGB stars, as compared to other mechanisms (e.g., AGN, low-level star formation). Additionally, we extend the tentative classification proposed in P13 by the type i+, which is assigned to a subset of type i ETGs exhibiting ongoing low-level star-formation (SF) in their periphery. This finding along with faint traces of localized SF in the extranuclear component of several of our sample ETGs points to a non-negligible contribution by OB stars to the total ionizing budget. We also demonstrate that, at the typical emission-line detection threshold of ~2Åin previous studies, most of the extranuclear wim emission in an ETG may evade detection, which could in turn prompt its classification as an entirely gas-devoid system. This study adds further observational evidence for a considerable heterogeneity among ETGs with regard to the physical properties and 2D kinematics of the wim component, and underscores the importance of IFS studies over their entire optical extent.

preprint2015arXiv

Tracing kinematic (mis)alignments in CALIFA merging galaxies: Stellar and ionized gas kinematic orientations at every merger stage

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 ($δ$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 $δ$PA$_{\mathrm{kin}}$ are large for both samples. However, for a large fraction of interacting galaxies the ionized gas $δ$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.

preprint2014arXiv

CALIFA across the Hubble types: Spatially resolved properties of the stellar populations

We analyze the spatially resolved star formation history of 300 nearby galaxies from the CALIFA integral field spectroscopic survey to investigate the radial structure and gradients of the present day stellar populations properties as a function of Hubble type and galaxy stellar mass. A fossil record method based on spectral synthesis techniques is used to recover spatially and temporally resolved maps of stellar population properties of spheroidal and spiral galaxies with masses $10^9$ to $7 \times 10^{11}$ M$_\odot$. The results show that galaxy-wide spatially averaged stellar population properties (stellar mass, mass surface density, age, metallicity, and extinction) match those obtained from the integrated spectrum, and that these spatially averaged properties match those at $R = 1$ HLR (half light radius), proving that the effective radii are really effective. Further, the individual radial profiles of the stellar mass surface density ($μ_\star$), luminosity weighted ages ($< {\rm log}\,age>_L$), and mass weighted metallicity ($< \log Z_\star >_M$) are stacked in bins of galaxy morphology (E, S0, Sa, Sb, Sbc, Sc and Sd). All these properties show negative gradients as a sign of the inside-out growth of massive galaxies. However, the gradients depend on the Hubble type in different ways. For the same galaxy mass, E and S0 galaxies show the largest inner gradients in $μ_\star$; while MW-like galaxies (Sb with $M_\star \sim 10^{11} M_\odot$) show the largest inner age and metallicity gradients. The age and metallicity gradients suggest that major mergers have a relevant role in growing the center (within 3 HLR) of massive early type galaxies; and radial mixing may play a role flattening the radial metallicity gradient in MW-like disks.

preprint2014arXiv

CALIFA, the Calar Alto Legacy Integral Field Area survey. III. Second public data release

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 Å with a spectral resolution of 6.0 Å (FWHM), and (ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 Å with a spectral resolution of 2.3 Å (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.

preprint2014arXiv

CALIFA: a diameter-selected sample for an integral field spectroscopy galaxy survey

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.

preprint2014arXiv

Census of HII regions in NGC 6754 derived with MUSE: Constraints on the metal mixing scale

We present a study of the HII regions in the galaxy NGC 6754 from a two pointing mosaic comprising 197,637 individual spectra, using Integral Field Spectrocopy (IFS) recently acquired with the MUSE instrument during its Science Verification program. The data cover the entire galaxy out to ~2 effective radii (re ), sampling its morphological structures with unprecedented spatial resolution for a wide-field IFU. A complete census of the H ii regions limited by the atmospheric seeing conditions was derived, comprising 396 individual ionized sources. This is one of the largest and most complete catalogue of H ii regions with spectroscopic information in a single galaxy. We use this catalogue to derive the radial abundance gradient in this SBb galaxy, finding a negative gradient with a slope consistent with the characteristic value for disk galaxies recently reported. The large number of H ii regions allow us to estimate the typical mixing scale-length (rmix ~0.4 re ), which sets strong constraints on the proposed mechanisms for metal mixing in disk galaxies, like radial movements associated with bars and spiral arms, when comparing with simulations. We found evidence for an azimuthal variation of the oxygen abundance, that may be related with the radial migration. These results illustrate the unique capabilities of MUSE for the study of the enrichment mechanisms in Local Universe galaxies.

preprint2014arXiv

Insights on the stellar mass-metallicity relation from the CALIFA survey

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 ($μ$*) 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 $μ$* and Z* (the $μ$ZR), betraying the influence of local factors in determining Z*. This shows that both local ($μ$*-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, $μ$* regulates Z*, producing a strong $μ$ZR whose amplitude is modulated by M*. In spheroids it is M* who dominates the physics of star formation and chemical enrichment, with $μ$* 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.

preprint2014arXiv

Integral field spectroscopy of nearby QSOs: I. ENLR size-luminosity relation, ongoing star formation & resolved gas-phase metallicities

[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.

preprint2014arXiv

Ionizing stellar population in the disc of NGC 3310 - II. The Wolf-Rayet population

We use integral field spectroscopy to study in detail the Wolf-Rayet (WR) population in NGC 3310, spatially resolving 18 star-forming knots with typical sizes of 200-300 pc in the disc of the galaxy hosting a substantial population of WRs. The detected emission in the so-called blue bump is attributed mainly to late-type nitrogen WRs (WNL), ranging from a few dozens to several hundreds of stars per region. Our estimated WNL/(WNL+O) ratio is comparable to reported empirical relations once the extinction-corrected emission is further corrected by the presence of dust grains inside the nebula that absorb a non-negligible fraction of UV photons. Comparisons of observables with stellar population models show disagreement by factors larger than 2-3. However, if the effects of interacting binaries and/or photon leakage are taken into account, observations and predictions tend to converge. We estimate the binary fraction of the \hii regions hosting WRs to be significant in order to recover the observed X-ray flux, hence proving that the binary channel can be critical when predicting observables. We also explore the connection of the environment with the current hypothesis that WRs can be progenitors to long-duration gamma-ray bursts (GRBs). Galaxy interactions, which can trigger strong episodes of star formation in the central regions, may be a plausible environment where WRs may act as progenitors of GRBs. Finally, even though the chemical abundance is generally homogeneous, we also find weak evidence for rapid N pollution by WR stellar winds at scales of ~ 200 pc.

preprint2014arXiv

Ionizing stellar population in the disk of NGC 3310 $-$ I. The impact of a minor merger on galaxy evolution

Numerical simulations of minor mergers predict little enhancement in the global star formation activity. However, it is still unclear the impact they have on the chemical state of the whole galaxy and on the mass build-up in the galaxy bulge and disc. We present a two-dimensional analysis of NCG 3310, currently undergoing an intense starburst likely caused by a recent minor interaction, using data from the PPAK Integral Field Spectroscopy (IFS) Nearby Galaxies Survey (PINGS). With data from a large sample of about a hundred HII regions identified throughout the disc and spiral arms we derive, using strong-line metallicity indicators and direct derivations, a rather flat gaseous abundance gradient. Thus, metal mixing processes occurred, as in observed galaxy interactions. Spectra from PINGS data and additional multiwavelength imaging were used to perform a spectral energy distribution fitting to the stellar emission and a photoionization modelling of the nebulae. The ionizing stellar population is characterized by single populations with a narrow age range (2.5-5 Myr) and a broad range of masses ($10^4-6\times10^6 M_\odot$). The effect of dust grains in the nebulae is important, indicating that 25-70% of the ultraviolet photons can be absorbed by dust. The ionizing stellar population within the HII regions represents typically a few percent of the total stellar mass. This ratio, a proxy to the specific star formation rate, presents a flat or negative radial gradient. Therefore, minor interactions may indeed play an important role in the mass build-up of the bulge.

preprint2014arXiv

Kinematic alignment of non-interacting CALIFA galaxies: Quantifying the impact of bars on stellar and ionised gas velocity field orientations

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.

preprint2014arXiv

Spectral synthesis of stellar populations in the 3D era: The CALIFA experience

Methods to recover the fossil record of galaxy evolution encoded in their optical spectra have been instrumental in processing the avalanche of data from mega-surveys along the last decade, effectively transforming observed spectra onto a long and rich list of physical properties: from stellar masses and mean ages to full star formation histories. This promoted progress in our understanding of galaxies as a whole. Yet, the lack of spatial resolution introduces undesirable aperture effects, and hampers advances on the internal physics of galaxies. This is now changing with 3D surveys. The mapping of stellar populations in data-cubes allows us to figure what comes from where, unscrambling information previously available only in integrated form. This contribution uses our starlight-based analysis of 300 CALIFA galaxies to illustrate the power of spectral synthesis applied to data-cubes. The selected results highlighted here include: (a) The evolution of the mass-metallicity and mass-density-metallicity relations, as traced by the mean stellar metallicity. (b) A comparison of star formation rates obtained from Hα to those derived from full spectral fits. (c) The relation between star formation rate and dust optical depth within galaxies, which turns out to mimic the Schmidt-Kennicutt law. (d) PCA tomography experiments.

preprint2014arXiv

The CALIFA survey: Status Report

We present here a brief summary of the status of the on-going CALIFA survey with an em- phasis on the results that have been recently published. In particular, we make a summary of the most relevant results found regarding the properties of Hii regions discovered using this survey, and the evidence uncovered for an inside-out growth of galaxies. In particular, we present an updated version of our results on the exploration of the Mass-Metallicity relation using spatial resolved information, with almost the double number of objects analysed, and an enlarged number of galaxies in the lower-mass range. We confirm our previous results that shows a lack of secondary relation with the SFR.

preprint2013arXiv

Aperture corrections for disk galaxy properties derived from the CALIFA survey. Balmer emission lines in spiral galaxies

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.

preprint2013arXiv

CALIFA, the Calar Alto Legacy Integral Field Area survey. II. First public data release

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.

preprint2013arXiv

Integral field spectroscopy of M1-67. A Wolf-Rayet nebula with LBVN appearance

This work aims to disentangle the morphological, kinematic, and chemical components of the nebula M1-67 to shed light on its process of formation around the central Wolf-Rayet (WR) star WR124. We have carried out integral field spectroscopy observations over two regions of M1-67, covering most of the nebula in the optical range. Maps of electron density, line ratios, and radial velocity were created to perform a detailed analysis of the two-dimensional structure. We studied the physical and chemical properties by means of integrated spectra selected over the whole nebula. Photoionization models were performed to confirm the empirical chemical results theoretically. In addition, we analysed infrared spectroscopic data and the MIPS 24micron image of M1-67 from Spitzer. We find that the ionized gas of M1-67 is condensed in knots aligned in a preferred axis along the NE-SW direction, like a bipolar structure. Both electron density and radial velocity decrease in this direction when moving away from the central star. From the derived electron temperature, Te~8200 K, we have estimated chemical abundances, obtaining that nitrogen appears strongly enriched and oxygen depleted. From the last two results, we infer that this bipolarity is the consequence of an ejection of an evolved stage of WR124 with material processed in the CNO cycle. The infrared study has revealed that the bipolar axis is composed of ionized gas with a low ionization degree that is well mixed with warm dust and of a spherical bubble surrounding the ejection at 24micron. Taking the evolution of a 60 Mo star and the temporal scale of the bipolar ejection into account, we propose that the observed gas was ejected during an eruption in the luminous blue variable. The star has entered the WR phase recently without apparent signs of interaction between WR-winds and interstellar material.

preprint2013arXiv

The evolution of galaxies resolved in space and time: an inside-out growth view from the CALIFA survey

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.

preprint2013arXiv

The nature of LINER galaxies: Ubiquitous hot old stars and rare accreting black holes

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.

preprint2013arXiv

The O3N2 and N2 abundance indicators revisited: improved calibrations based on CALIFA and Te-based literature data

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érez-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.

preprint2013arXiv

The star formation history of CALIFA galaxies: Radial structures

We study the radial structure of the stellar mass surface density ($μ$) 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 $μ$ 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 $μ$ gradients that steepen with stellar mass. In spheroidal galaxies $μ$ saturates at a critical value that is independent of the galaxy mass. Thus, all the massive spheroidal galaxies have similar local $μ$ at the same radius (in HLR units). The SFH of the regions beyond 1 HLR are well correlated with their local $μ$, and follow the same relation as the galaxy-averaged age and $μ$; 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 $μ$. Thus, galaxy mass is a more fundamental property in spheroidal systems while the local stellar mass surface density is more important in disks.

preprint2012arXiv

A new scaling relation for HII regions in spiral galaxies: unveiling the true nature of the mass-metallicity relation

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.

preprint2012arXiv

Ionization structure and chemical abundances of the Wolf-Rayet nebula NGC6888 with integral field spectroscopy

This work aims to search for the observational footprints of the interactions between the interstellar medium (ISM) and stellar winds in the Wolf-Rayet (WR) nebula NGC6888 in order to understand its ionization structure, chemical composition, and kinematics. We have collected a set of integral field spectroscopy observations across NGC6888, obtained with PPAK in the optical range performing both 2D and 1D analyses. Attending to the 2D analysis in the northeast part of NGC6888, we have generated maps of the extinction structure and electron density. We produced statistical frequency distributions of the radial velocity and diagnostic diagrams. We have found that the spectra of a localized region to the southwest of this pointing can be represented well by shock models. Furthermore, we performed a thorough study of integrated spectra in nine regions over the whole nebula. We derived electron densities ranging from <100 to 360 cm^(-3). The electron temperature varies from ~7700 K to ~10200 K. A strong variation of up to a factor 10 between different regions in the nitrogen abundance has been found: N/H appears lower than the solar abundance in those positions observed at the edges and very enhanced in the observed inner parts. Oxygen appears slightly underabundant with respect to solar value, whereas the helium abundance is found to be above it. Finally, we provide a scenario for the evolution of NGC6888 to explain the features observed. This scheme consists of a structure of multiple shells.

preprint2012arXiv

PyCosmic: a robust method to detect cosmics in CALIFA and other fiber-fed integral-field spectroscopy datasets

[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.

preprint2012arXiv

The oxygen abundance in the IFU era

Spatially-resolved information of gas-phase emission provided by integral field units (IFUs) are allowing us to perform a new generation of emission-line surveys, based on large samples of HII regions and full two-dimensional coverage. Here we present two highlights of our current studies employing this technique: 1) A statistical approach to the abundance gradients of spiral galaxies, which indicates an -universal- radial gradient for oxygen abundance; and 2) The discovery of a new scaling relation of HII regions in spiral galaxies, the "local" mass-metallicity relation of star-forming galaxies.

preprint2012arXiv

The properties of the extended warm ionised gas around low-redshift QSOs and the lack of extended high-velocity outflows

(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.

preprint2011arXiv

PPAK Wide-field Integral Field Spectroscopy of NGC 628: II. Emission line abundance analysis

In this second paper of the series, we present the 2-dimensional (2D) emission line abundance analysis of NGC 628, the largest object within the PPAK Integral Field Spectroscopy (IFS) Nearby Galaxies Survey: PINGS. We introduce the methodology applied to the 2D IFS data in order to extract and deal with large spectral samples, from which a 2D abundance analysis can be later performed. We obtain the most complete and reliable abundance gradient of the galaxy up-to-date, by using the largest number of spectroscopic points sampled in the galaxy, and by comparing the statistical significance of different strong-line metallicity indicators. We find features not previously reported for this galaxy that imply a multi-modality of the abundance gradient consistent with a nearly flat-distribution in the innermost regions of the galaxy, a steep negative gradient along the disc and a shallow gradient or nearly-constant metallicity beyond the optical edge of the galaxy. The N/O ratio seems to follow the same radial behaviour. We demonstrate that the observed dispersion in metallicity shows no systematic dependence with the spatial position, signal-to-noise or ionization conditions, implying that the scatter in abundance for a given radius is reflecting a true spatial physical variation of the oxygen content. Furthermore, by exploiting the 2D IFS data, we were able to construct the 2D metallicity structure of the galaxy, detecting regions of metal enhancement, and showing that they vary depending on the choice of the metallicity estimator. The analysis of axisymmetric variations in the disc of NGC 628 suggest that the physical conditions and the star formation history of different-symmetric regions of the galaxy have evolved in a different manner.

preprint2011arXiv

The low-metallicity QSO HE 2158-0107: A massive galaxy growing by the accretion of nearly pristine gas from its environment?

[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.

preprint2010arXiv

CALIFA, the Calar Alto Legacy Integral Field Area survey: Early Report

We present the Calar Alto Legacy Integral Field Area survey (CALIFA). CALIFA's main aim is to obtain spatially resolved spectroscopic information for ~600 galaxies of all Hubble types in the Local Universe (0.005< z <0.03). The survey has been designed to allow three key measurements to be made: (a) Two-dimensional maps of stellar populations (star formation histories, chemical elements); (b) The distribution of the excitation mechanism and element abundances of the ionized gas; and (c) Kinematic properties (velocity fields, velocity dispersion), both from emission and from absorption lines. To cover the full optical extension of the target galaxies (i.e. out to a 3sigma depth of ~23 mag/arcsec2), CALIFA uses the exceptionally large field of view of the PPAK/PMAS IFU at the 3.5m telescope of the Calar Alto observatory. We use two grating setups, one covering the wavelength range between 3700 and 5000 AA at a spectral resolution R~1650, and the other covering 4300 to 7000 AA at R~850. The survey was allocated 210 dark nights, distributed in 6 semesters and starting in July 2010 and is carried out by the CALIFA collaboration, comprising ~70 astronomers from 8 different countries. As a legacy survey, the fully reduced data will be made publically available, once their quality has been verified. We showcase here early results obtained from the data taken so far (21 galaxies).

preprint2010arXiv

Integral Field Spectroscopy of Local LCBGs: NGC 7673, a case study. Physical properties of star-forming regions

Physical properties of the star-forming regions in the local Luminous Compact Blue Galaxy NGC 7673 are studied in detail using 3D spectroscopic data taken with the PPAK IFU at the 3.5-m telescope in CAHA. We derive integrated and spatially resolved properties such as extinction, star formation rate and metallicity for this galaxy. Our data show an extinction map with maximum values located at the position of the main clumps of star formation showing small spatial variations (E(B-V)_{t}=0.12-0.21 mag). We derive a Hα-based SFR for this galaxy of 6.2 \pm 0.8 M_{\odot}/yr in agreement with the SFR derived from infrared and radio continuum fluxes. The star formation is located mainly in clumps A, B, C and F. Different properties measured in clump B makes this region peculiar. We find the highest Hαluminosity with a SFR surface density of 0.5 M_{\odot}yr^{-1}kpc^{-2} in this clump. In our previous work, the kinematic analysis for this galaxy shows an asymmetrical ionized gas velocity field with a kinematic decoupled component located at the position of clump B. This region shows the absence of strong absorption features and the presence of a Wolf-Rayet stellar population indicating this is a young burst of massive stars. Furthermore, we estimate a gas metallicity of 12+log(O/H)=8.20\pm0.15 for the integrated galaxy using the R23 index. The values derived for the different clumps with this method show small metallicity variations in this galaxy, with values in the range 8.12 (for clump A) - 8.23 (for clump B) for 12+log(O/H). The analysis of the emission line ratios discards the presence of any AGN activity or shocks as the ionization source in this galaxy. Between the possible mechanisms to explain the starburst activity in this galaxy, our 3D spectroscopic data support the scenario of an on-going interaction with the possibility for clump B to be the dwarf satellite galaxy.

preprint2010arXiv

Mapping the ionised gas around the luminous QSO HE 1029-1401: Evidence for minor merger events?

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.

preprint2009arXiv

3D Spectroscopy of Local Luminous Compact Blue Galaxies: Kinematics of NGC 7673

The kinematic properties of the ionized gas of local Luminous Compact Blue Galaxy (LCBG) NGC 7673 are presented using three dimensional data taken with the PPAK integral field unit at the 3.5-m telescope in the Centro Astronómico Hispano Alemán. Our data reveal an asymmetric rotating velocity field with a peak to peak difference of 60 km s$^{-1}$. The kinematic centre is found to be at the position of a central velocity width maximum ($σ=54\pm1$ km s$^{-1}$), which is consistent with the position of the luminosity-weighted centroid of the entire galaxy. The position angle of the minor rotation axis is 168$^{\circ}$ as measured from the orientation of the velocity field contours. At least two decoupled kinematic components are found. The first one is compact and coincides with the position of the second most active star formation region (clump B). The second one is extended and does not have a clear optical counterpart. No evidence of active galactic nuclei activity or supernovae galactic winds powering any of these two components has been found. Our data, however, show evidence in support of a previously proposed minor merger scenario in which a dwarf galaxy, tentatively identified with clump B, is falling into NGC 7673. and triggers the starburst. Finally, it is shown that the dynamical mass of this galaxy may be severely underestimated when using the derived rotation curve or the integrated velocity width, under the assumption of virialization.