First science with SHARDS: emission line galaxies

SHARDS is an unbiased ultra-deep spectro-photometric survey with GTC@OSIRIS aimed at selecting and studying massive passively evolving galaxies at z=1.0-2.3 using a set of 24 medium-band filters (FWHM~17 nm) at 500-950 nm in GOODS-N. Our observing strategy is optimized to detect at z>1 the prominent Mg absorption feature at rest-frame ~ 280 nm, a distinctive, necessary, and sufficient feature of evolved stellar populations. Nonetheless, the data quality allow a plethora of studies on galaxy populations, including Emission Lines Galaxies (ELGs) about which we have started our first science verification project presented in this contribution.

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

Context. 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 (more commonly known as active galactic nucleus, AGN) of low luminosity, making LINER galaxies the largest AGN sub-population, which dominate 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. Aims. Using integral field spectroscopic data from the CALIFA survey, we compare the observed radial surface brightness profiles with what is expected from illumination by an AGN. Methods. Essential for this analysis is a proper extraction of emission lines, especially weak lines, such as Balmer H beta lines, which are superposed on an absorption trough. To accomplish this, we use the GANDALF code, which simultaneously fits the underlying stellar continuum and emission lines. Results. For 48 galaxies with LINER-like emission, we show 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. Conclusions. 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 the gas present and with 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 not a class defined by a property but rather by the absence of a property. It also explains why LINER emission is observed mostly in massive galaxies with old stars and little star formation.

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

Chromospheric activity and rotation of FGK stars in the solar vicinity. An estimation of the radial velocity jitter

Context. Chromospheric activity produces both photometric and spectroscopic variations that can be mistaken as planets. Large spots crossing the stellar disc can produce planet-like periodic variations in the light curve of a star. These spots clearly affect the spectral line profiles, and their perturbations alter the line centroids creating a radial velocity jitter that might “contaminate” the variations induced by a planet. Precise chromospheric activity measurements are needed to estimate the activity-induced noise that should be expected for a given star. Aims. We obtain precise chromospheric activity measurements and projected rotational velocities for nearby (d ≤ 25 pc) cool (spectral types F to K) stars, to estimate their expected activity-related jitter. As a complementary objective, we attempt to obtain relationships between fluxes in different activity indicator lines, that permit a transformation of traditional activity indicators, i.e., Ca II H & K lines, to others that hold noteworthy advantages. Methods. We used high resolution (~50 000) echelle optical spectra. Standard data reduction was performed using the IRAF ECHELLE package. To determine the chromospheric emission of the stars in the sample, we used the spectral subtraction technique. We measured the equivalent widths of the chromospheric emission lines in the subtracted spectrum and transformed them into fluxes by applying empirical equivalent width and flux relationships. Rotational velocities were determined using the cross-correlation technique. To infer activity-related radial velocity (RV) jitter, we used empirical relationships between this jitter and the R’_HK index. Results. We measured chromospheric activity, as given by different indicators throughout the optical spectra, and projected rotational velocities for 371 nearby cool stars. We have built empirical relationships among the most important chromospheric emission lines. Finally, we used the measured chromospheric activity to estimate the expected RV jitter for the active stars in the sample.

Optical and ultraviolet observations of a strong flare in the young, single K2 dwarf LQ Hya

We present high-resolution optical echelle spectra and IUE observations during a strong flare on 1993 December 22 in the very active, young, rapidly rotating, single K2 dwarf LQ Hya. The initial impulsive phase of the flare, which started sometime between 2:42 ut and 4:07 ut, was characterized by strong optical continuum enhancement and blueshifted emission lines with broad wings. The optical chromospheric lines reached their maximum intensity at ≈ 5:31 ut, by which time the blueshift vanished and the optical continuum enhancement had sharply decreased. Thereafter, the line emission slowly decreased and the lines redshift in a gradual phase that lasted at least two more hours. The Mg II lines behaved similarly. Quiescent C IV flux levels were not recovered until 21 h later, though a data gap and a possible second flare make the interpretation uncertain. In addition to the typically flare-enhanced emission lines (e.g., H α and H β), we observe He I D_3 going into emission, plus excess emission (after subtraction of the quiescent spectrum) in other He I and several strong neutral metal lines (e.g., Mg I b). Flare enhancement of the far-ultraviolet continuum generally agrees with an Si I recombination model. We estimate the total flare energy, and discuss the broad components, asymmetries and Doppler shifts seen in some of the emission lines.

Caught in the act: gas and stellar velocity dispersions in a fast Quenching compact star-forming galaxy at z ~ 1.7

We present Keck I MOSFIRE spectroscopy in the Y and H bands of GDN-8231, a massive, compact, star-forming galaxy at a redshift of z ~ 1.7. Its spectrum reveals both Hα and [Nii] emission lines and strong Balmer absorption lines. The Hα and Spitzer MIPS 24 μm fluxes are both weak, thus indicating a low star-formation rate of SFR≲5-10 M_⨀ yr−1. This, added to a relatively young age of ~700 Myr measured from the absorption lines, provides the first direct evidence for a distant galaxy being caught in the act of rapidly shutting down its star formation. Such quenching allows GDN-8231 to become a compact, quiescent galaxy, similar to three other galaxies in our sample, by z ~ 1.5. Moreover, the color profile of GDN-8231 shows a bluer center, consistent with the predictions of recent simulations for an early phase of inside-out quenching. Its line-of-sight velocity dispersion for the gas, σ_LOG^gas = 127 ± 32 km s^−1, is nearly 40% smaller than that of its stars, σ_LOG^* = 215 ± 35 km s^−1. High-resolution hydro-simulations of galaxies explain such apparently colder gas kinematics of up to a factor of ~1.5 with rotating disks being viewed at different inclinations and/or centrally concentrated star-forming regions. A clear prediction is that their compact, quiescent descendants preserve some remnant rotation from their star-forming progenitors.

Stars and brown dwarfs in the σ Orionis cluster (Research Note) III. OSIRIS/GTC low-resolution spectroscopy of variable sources

Context. Although many studies have been performed so far, there are still dozens of low-mass stars and brown dwarfs in the young σ Orionis open cluster without detailed spectroscopic characterisation. Aims. We look for unknown strong accretors and disc hosts that were undetected in previous surveys. Methods. We collected low-resolution spectroscopy (R ~ 700) of ten low-mass stars and brown dwarfs in σ Orionis with OSIRIS at the Gran Telescopio Canarias under very poor weather conditions. These objects display variability in the optical, infrared, Hα, and/or X-rays on time scales of hours to years. We complemented our spectra with optical and near-/mid-infrared photometry. Results. For seven targets, we detected lithium in absorption, identified Hα, the calcium doublet, and forbidden lines in emission, and/or determined spectral types for the first time. We characterise in detail a faint, T Tauri-like brown dwarf with an 18 h-period variability in the optical and a large Hα equivalent width of –125  ±  15 Å, as well as two M1-type, X-ray-flaring, low-mass stars, one with a warm disc and forbidden emission lines, the other with a previously unknown cold disc with a large inner hole. Conclusions. New unrevealed strong accretors and disc hosts, even below the substellar limit, await discovery among the list of known σ Orionis stars and brown dwarfs that are variable in the optical and have no detailed spectroscopic characterisation yet.

An X-ray study of the SNR G344.7-0.1 and the central object CXOU J170357.8-414302

Aims. We report results of an X-ray study of the supernova remnant (SNR) G344.7-0.1 and the point-like X-ray source located at the geometrical center of the SNR radio structure. Methods. The morphology and spectral properties of the remnant and the central X-ray point-like source were studied using data from the XMM-Newton and Chandra satellites. Archival radio data and infrared Spitzer observations at 8 and 24 mu m were used to compare and study its multi-band properties at different wavelengths. Results. The XMM-Newton and Chandra observations reveal that the overall X-ray emission of G344.7-0.1 is extended and correlates very well with regions of bright radio and infrared emission. The X-ray spectrum is dominated by prominent atomic emission lines. These characteristics suggest that the X-ray emission originated in a thin thermal plasma, whose radiation is represented well by a plane-parallel shock plasma model (PSHOCK). Our study favors the scenario in which G344.7-0.1 is a 6 x 10^3 year old SNR expanding in a medium with a high density gradient and is most likely encountering a molecular cloud on the western side. In addition, we report the discovery of a soft point-like X-ray source located at the geometrical center of the radio SNR structure. The object presents some characteristics of the so-called compact central objects (CCO). However, its neutral hydrogen absorption column (N_H) is inconsistent with that of the SNR. Coincident with the position of the source, we found infrared and optical objects with typical early-K star characteristics. The X-ray source may be a foreground star or the CCO associated with the SNR. If this latter possibility were confirmed, the point-like source would be the farthest CCO detected so far and the eighth member of the new population of isolated and weakly magnetized neutron stars.

Local luminous infrared galaxies: Spatially resolved mid-infrared observations with Spitzer/IRS

Luminous Infrared (IR) Galaxies (LIRGs, L_IR=10^11-10 L_⨀) are an important cosmological class of galaxies as they are the main contributors to the co-moving star formation rate density of the universe at z=1. In this paper we present a guaranteed time observation (GTO) Spitzer InfraRed Spectrograph (IRS) program aimed to obtain spectral mapping of a sample of 14 local d<76Mpc LIRGs. The data cubes map, at least, the central 20arcsec X 20arcsec to 30 arcsec X 30 arcsec regions of the galaxies, and use all four IRS modules covering the full 5-38 μ m spectral range. The final goal of this project is to characterize fully the mid-IR properties of local LIRGs as a first step to understanding their more distant counterparts. In this paper we present the first results of this GTO program. The IRS spectral mapping data allow us to build spectral maps of the bright mid-IR emission lines (e.g., [Ne II] 12.81 μ m, [Ne III]15.56 μ m, [S III] 18.71 μ m, H_2 at 17 μ m), continuum, the 6.2 and 11.3 μ m polycyclic aromatic hydrocarbon (PAH) features, and the 9.7 μ m silicate feature, as well as to extract 1D spectra for regions of interest in each galaxy. The IRS data are used to obtain spatially resolved measurements of the extinction using the 9.7 μ m silicate feature, and to trace star forming regions using the neon lines and the PAH features. We also investigate a number of active galactic nuclei (AGN) indicators, including the presence of high excitation emission lines and a strong dust continuum emission at around 6 9.7 μ m . We finally use the integrated Spitzer/IRS spectra as templates of local LIRGs. We discuss several possible uses for these templates, including the calibration of the star formation rate of IR-bright galaxies at high redshift. We also predict the intensities of the brightest mid-IR emission lines for LIRGs as a function of redshift, and compare them with the expected sensitivities of future space IR missions.

Are long gamma-ray bursts biased tracers of star formation? Clues from the host galaxies of the Swift/BAT6 complete sample of bright LGRBs II. Star formation rates and metallicities at z < 1

Aims. Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and might therefore be a potentially powerful tool for tracing cosmic star formation. However, especially at low redshifts (z< 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the effect of the environment on GRB formation. Methods. We studied host galaxy spectra of the Swift/BAT6 complete sample of 14 z< 1 bright LGRBs. We used the detected nebular emission lines to measure the dust extinction, star formation rate (SFR), and nebular metallicity (Z) of the hosts and supplemented the data set with previously measured stellar masses M_*. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M_* relations) are compared to samples of field star-forming galaxies. Results. We find that LGRB hosts at z< 1 have on average lower SFRs than if they were direct star formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to 12 +log (O/H)~8.4−8.5, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are consistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of an increased incidence of starburst galaxies in the Swift/BAT6 z< 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless, the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevents us from investigating more quantitatively whether the starburst incidence is such as expected after taking into account the high-metallicity aversion of LGRB host galaxies.

Aperture effects on the oxygen abundance determinations from Califa data

This paper aims to provide 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([O III] λ5007/Hβ)/([N II] λ6583/Hα) (O3N2) and log[N II] lambda 6583/Hα (N2) flux ratios since they are closely connected to different empirical calibrations of the oxygen abundances in star-forming galaxies. We compute the median growth curves of Hα, Hα/Hβ, O3N2, and N-2 up to 2.5R(50) and 1.5 disk R-eff. These distances cover most of the optical spatial extent of the CALIFA galaxies. The growth curves simulate the effect of observing galaxies through apertures of varying radii. We split these growth curves by morphological types and stellar masses to check if there is any dependence on these properties. The median growth curve of the Hα flux shows a monotonous increase with radius with no strong dependence on galaxy inclination, morphological type, and stellar mass. The median growth curve of the Hα/HβH ratio monotonically decreases from the center toward larger radii, showing for small apertures a maximum value of ≈10% larger than the integrated one. It does not show any dependence on inclination, morphological type, and stellar mass. The median growth curve of N-2 shows a similar behavior, decreasing from the center toward larger radii. No strong dependence is seen on the inclination, morphological type, and stellar mass. Finally, the median growth curve of O3N2 increases monotonically with radius, and it does not show dependence on the inclination. 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.5 R_eff) given the high dispersion shown around the median growth curves. Thus, the application of these median aperture corrections to derive abundances for individual galaxies is not recommended when their fluxes come from radii much smaller than either R_50 or R_eff.

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.

Emission-line galaxy surveys with GTC: The Hα Universe

Emission line galaxies are the most easily detected and studied objects in the high redshift Universe. They are being used to trace the evolution of critical observables of the Universe such as Star Formation Rate densities, starburst properties and abundances. Most of the research is being done using [OII]3727 and UV lines, but Hα is still one the best tracers for Star Formation Rate and physical properties of current star-forming galaxies. As a complementary contribution to studies of galaxy evolution, our team has focused into a long-term project to study the population of Hα-selected star-forming galaxies of the Universe at different redshifts. In 1995 we first determined the local Hα luminosity function, and from it the Star Formation Rate density (SFRd) of the local Universe. We then, using narrow-band imaging in the optical, extended this measurement to z ≃ 0.24 and z ≃ 0.4. Working in the near-infrared, GTC will be a very powerful tool to study the evolution of the Hα emission-line galaxy populations at different redshifts. We will both quantify the SFRd evolution and characterize the star-forming galaxy populations by directly comparing the same observables at all redshifts up to z ≃ 2.5.