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

We present PIPE3D, an analysis pipeline based on the FIT3D fitting tool, developed to explore the properties of the stellar populations and ionized gas of integral field spectroscopy (IFS) data. PIPE3D was created to provide 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). In 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 example 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. 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.

Pipe3D, a pipeline to analyse integral Field Spectroscopy data: I. New fitting philosophy 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. This tool was developed to analyze integral field spectroscopy data and it is the basis of PIPE3D, a pipeline used in the analysis of CALIFA, MaNGA, and SAMI data. We describe the philosophy and each step of the fitting procedure. We present an extensive set of simulations in order to estimate the precision and accuracy of the derived parameters for the stellar populations and the ionized gas. We report on the results of those simulations. Finally, we compare the results of the analysis using FIT3D with those provided by other widely used packages, and we find that the parameters derived by FIT3D are fully compatible with those derived using these other tools.

Integral field spectroscopy and multi-wavelength imaging of nearby spiral galaxies: NGC 5668 as a pilot case for MEGARA

MEGARA (Multi-Espectrógrafo en GTC de Alta Resolución para Astronomía) is an optical Integral-Field Unit (IFU) and Multi-Object Spectrograph (MOS) designed for the GTC 10.4 m telescope in La Palma. MEGARA will be a 3rd generation instrument for GTC. It is led by the University Complutense of Madrid with the collaboration of INAOE, IAA, UPM and comprises more than 50 researchers from a large number of institutions worldwide.

The Mice at play in the CALIFA survey. A case study of a gas-rich major merger between first passage and coalescence

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (≳10^11 M_⊙) gas-rich spirals NGC 4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties, and stellar population diagnostics, over the full optical extent of both galaxies with ~1.6 kpc spatial resolution. The Mice galaxies provide a perfect case study that highlights the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars most likely induced in both galaxies. The barred spiral NGC 4676B exhibits a strong twist in both its stellar and ionised gas disk. The edge-on disk galaxy NGC 4676A appears to be bulge free, with a strong bar causing its “boxy” light profile. On the other hand, the impact of the merger on the stellar populations has been minimal thus far. By combining the IFS data with archival multiwavelength observations we show that star formation induced by the recent close passage has not contributed significantly to the total star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC 4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks (vs ~ 350 km s^-1) extend to ~6.6 kpc above the disk plane. The measured ram pressure (P/k = 4.8 × 10^6 K cm^-3) and mass outflow rate (~8−20 M_⊙ yr^-1) are similar to superwinds from local ultra-luminous infrared galaxies, although NGC 4676A only has a moderate infrared luminosity of 3 × 10^10 L_⊙. Energy beyond what is provided by the mechanical energy of the starburst appears to be required to drive the outflow. Finally, we compare the observations to mock kinematic and stellar population maps extracted from a hydrodynamical merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

Mass-metallicity relation explored with CALIFA I. Is there a dependence on the star-formation rate?

We studied the global and local ℳ-Z relation based on the first data available from the CALIFA survey (150 galaxies). This survey provides integral field spectroscopy of the complete optical extent of each galaxy (up to 2−3 effective radii), with a resolution high enough to separate individual H II regions and/or aggregations. About 3000 individual H II regions have been detected. The spectra cover the wavelength range between [OII]3727 and [SII]6731, with a sufficient signal-to-noise ratio to derive the oxygen abundance and star-formation rate associated with each region. In addition, we computed the integrated and spatially resolved stellar masses (and surface densities) based on SDSS photometric data. We explore the relations between the stellar mass, oxygen abundance and star-formation rate using this dataset. We derive a tight relation between the integrated stellar mass and the gas-phase abundance, with a dispersion lower than the one already reported in the literature (σ_Δlog (O/H) = 0.07 dex). Indeed, this dispersion is only slightly higher than the typical error derived for our oxygen abundances. However, we found no secondary relation with the star-formation rate other than the one induced by the primary relation of this quantity with the stellar mass. The analysis for our sample of ~3000 individual H II   regions confirms (i) a local mass-metallicity relation and (ii) the lack of a secondary relation with the star-formation rate. The same analysis was performed with similar results for the specific star-formation rate. Our results agree with the scenario in which gas recycling in galaxies, both locally and globally, is much faster than other typical timescales, such like that of gas accretion by inflow and/or metal loss due to outflows. In essence, late-type/disk-dominated galaxies seem to be in a quasi-steady situation, with a behavior similar to the one expected from an instantaneous recycling/closed-box model.

The ionized gas in the CALIFA early-type galaxies I. Mapping two representative cases: NGC 6762 and NGC 5966

As part of the ongoing CALIFA survey, we have conducted a thorough bidimensional analysis of the ionized gas in two E/S0 galaxies, NGC 6762 and NGC 5966, aiming to shed light on the nature of their warm ionized ISM. Specifically, we present optical (3745–7300 Å) integral field spectroscopy obtained with the PMAS/PPAK integral field spectrophotometer. Its wide field-of-view (1′ × 1′) covers the entire optical extent of each galaxy down to faint continuum surface brightnesses. To recover the nebular lines, we modeled and subtracted the underlying stellar continuum from the observed spectra using the STARLIGHT spectral synthesis code. The pure emission-line spectra were used to investigate the gas properties and determine the possible sources of ionization. We show the advantages of IFU data in interpreting the complex nature of the ionized gas in NGC 6762 and NGC 5966. In NGC 6762, the ionized gas and stellar emission display similar morphologies, while the emission line morphology is elongated in NGC 5966, spanning ~6 kpc, and is oriented roughly orthogonal to the major axis of the stellar continuum ellipsoid. Whereas gas and stars are kinematically aligned in NGC 6762, the gas is kinematically decoupled from the stars in NGC 5966. A decoupled rotating disk or an “ionization cone” are two possible interpretations of the elongated ionized gas structure in NGC 5966. The latter would be the first “ionization cone” of such a dimension detected within a weak emission-line galaxy. Both galaxies have weak emission-lines relative to the continuum[EW(Hα) ≲ 3 Å] and have very low excitation, log([OIII]λ5007/Hβ) ≲ 0.5. Based on optical diagnostic ratios ([OIII]λ5007/Hβ, [NII]λ6584/Hα, [SII]λ6717, 6731/Hα, [OI]λ6300/Hα), both objects contain a LINER nucleus and an extended LINER-like gas emission. The emission line ratios do not vary significantly with radius or aperture, which indicates that the nebular properties are spatially homogeneous. The gas emission in NGC 6762 can be best explained by photoionization by pAGB stars without the need of invoking any other excitation mechanism. In the case of NGC 5966, the presence of a nuclear ionizing source seems to be required to shape the elongated gas emission feature in the “ionization cone” scenario, although ionization by pAGB stars cannot be ruled out. Further study of this object is needed to clarify the nature of its elongated gas structure.

The multiphase starburst-driven galactic wind in NGC 5394

We present a detailed study of the neutral and ionized gas phases in the galactic wind for the nearby starburst galaxy NGC 5394 based on new integral field spectroscopy obtained with the INTEGRAL fibre system at the William Herschel Telescope. The neutral gas phase in the wind is detected via the interstellar Na I D doublet absorption. After a careful removal of the stellar contribution to these lines, a significant amount of neutral gas (∼10^7 M_⊙) is detected in a central region of ∼1.75 kpc size. This neutral gas is blueshifted by ∼165 km s^−1 with respect to the underlying galaxy. The mass outflow of neutral gas is comparable to the star formation rate of the host galaxy. Simultaneously, several emission lines (Hα, [N II], [S II]) are also analysed looking for the ionized warm phase counterpart of the wind. A careful kinematic decomposition of the line profiles reveals the presence of a secondary, broader, kinematic component. This component is found roughly in the same region where the Na I D absorption is detected. It presents higher [N II]/Hα and [S II]/Hα line ratios than the narrow component at the same locations, indicative of contamination by shock ionization. This secondary component also presents blueshifted velocities, although smaller than those measured for the neutral gas, averaging to ∼−30 km s^−1. The mass and mass outflow rate of the wind is dominated by the neutral gas, of which a small fraction might be able to escape the gravitational potential of the host galaxy. The observations in this system can be readily understood within a bipolar gas flow scenario.

A spectroscopy study of nearby late-type stars, possible members of stellar kinematic groups

Context. Nearby late-type stars are excellent targets for seeking young objects in stellar associations and moving groups. The origin of these structures is still misunderstood, and lists of moving group members often change with time and also from author to author. Most members of these groups have been identified by means of kinematic criteria, leading to an important contamination of previous lists by old field stars. Aims. We attempt to identify unambiguous moving group members among a sample of nearby-late type stars by studying their kinematics, lithium abundance, chromospheric activity, and other age-related properties. Methods. High-resolution echelle spectra (R ~ 57 000) of a sample of nearby late-type stars are used to derive accurate radial velocities that are combined with the precise Hipparcos parallaxes and proper motions to compute galactic-spatial velocity components. Stars are classified as possible members of the classical moving groups according to their kinematics. The spectra are also used to study several age-related properties for young late-type stars, i.e., the equivalent width of the lithium Li i 6707.8 Å line or the R'_HK index. Additional information like X-ray fluxes from the ROSAT All-Sky Survey or the presence of debris discs is also taken into account. The different age estimators are compared and the moving group membership of the kinematically selected candidates are discussed. Results. From a total list of 405 nearby stars, 102 have been classified as moving group candidates according to their kinematics. i.e., only ~25.2% of the sample. The number reduces when age estimates are considered, and only 26 moving group candidates (25.5% of the 102 candidates) have ages in agreement with the star having the same age as an MG member.

Large-scale clustering measurements with photometric redshifts: comparing the dark matter haloes of X-ray AGN, star-forming and passive galaxies at z ≈ 1

We combine multi-wavelength data in the AEGIS-XD and C-COSMOS surveys to measure the typical dark matter halo mass of X-ray selected active galactic nuclei (AGN) [L_X(2–10 keV) > 10^42 erg s^− 1] in comparison with far-infrared selected star-forming galaxies detected in the Herschel/PEP survey (PACS Evolutionary Probe; L_IR > 10^11 L_⊙) and quiescent systems at z ≈ 1. We develop a novel method to measure the clustering of extragalactic populations that uses photometric redshift probability distribution functions in addition to any spectroscopy. This is advantageous in that all sources in the sample are used in the clustering analysis, not just the subset with secure spectroscopy. The method works best for large samples. The loss of accuracy because of the lack of spectroscopy is balanced by increasing the number of sources used to measure the clustering. We find that X-ray AGN, far-infrared selected star-forming galaxies and passive systems in the redshift interval 0.6 < z < 1.4 are found in haloes of similar mass, log M_DMH/(M_⊙ h^−1) ≈ 13.0. We argue that this is because the galaxies in all three samples (AGN, star-forming, passive) have similar stellar mass distributions, approximated by the J-band luminosity. Therefore, all galaxies that can potentially host X-ray AGN, because they have stellar masses in the appropriate range, live in dark matter haloes of log M_DMH/(M_⊙ h^−1) ≈ 13.0 independent of their star formation rates. This suggests that the stellar mass of X-ray AGN hosts is driving the observed clustering properties of this population. We also speculate that trends between AGN properties (e.g. luminosity, level of obscuration) and large-scale environment may be related to differences in the stellar mass of the host galaxies.