Stellar populations in bulges of spiral galaxies

We present a few results of a work in progress tackling the radial spectroscopic properties of bulges of spiral galaxies.

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.

The star formation history of CALIFA galaxies: Radial structures

We have studied 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 applied the fossil record method based on spectral synthesis techniques to recover the star formation history (SFH), resolved in space and time, in spheroidal and disk dominated galaxies with masses from 10^9 to 10^12 M_⊙. We derived the half-mass radius, and we found that galaxies are on average 15% more compact in mass than in light. The ratio of half-mass radius to half-light radius (HLR) shows a dual dependence with galaxy stellar mass; it decreases with increasing mass for disk galaxies, but is almost constant in spheroidal galaxies. In terms of integrated versus spatially resolved properties, 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 are present in most of the galaxies, supporting an inside-out formation. The larger inner (≤1 HLR) age gradients occur in the most massive (10^11 M_⊙) 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 (~7 × 10^2 M_⊙/pc^2 at 1 HLR) that is independent of the galaxy mass. Thus, all the massive spheroidal galaxies have similar local μ_∗ at the same distance (in HLR units) from the nucleus. 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 μ_∗; this suggests 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.

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.

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.

DUst around NEarby Stars. The survey observational results

Context. Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system counterparts are the asteroid and Edgeworth-Kuiper belts. Aims. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. Methods. We used Herschel/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 μm were obtained, complemented in some cases with observations at 70 μm, and at 250, 350 and 500 μm using SPIRE. The observing strategy was to integrate as deep as possible at 100 μm to detect the stellar photosphere. Results. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of ~12.1% ± 5% before Herschel to ~20.2% ± 2%. A significant fraction (~52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70–160 μm range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

Chemically tagging the Hyades Supercluster. A homogeneous sample of F6-K4 kinematically selected northern stars

Stellar kinematic groups are kinematical coherent groups of stars that might have a common origin. These groups are dispersed throughout the Galaxy over time by the tidal effects of both Galactic rotation and disc heating, although their chemical content remains unchanged. The aim of chemical tagging is to establish that the abundances of every element in the analysis are homogeneus among the members. We study the case of the Hyades Supercluster to compile a reliable list of members (FGK stars) based on our chemical tagging analysis. For a total of 61 stars from the Hyades Supercluster, stellar atmospheric parameters (T_eff, log g, ξ, and [Fe/H]) are determined using our code called StePar, which is based on the sensitivity to the stellar atmospheric parameters of the iron EWs measured in the spectra. We derive the chemical abundances of 20 elements and find that their [X/Fe] ratios are consistent with Galactic abundance trends reported in previous studies. The chemical tagging method is applied with a carefully developed differential abundance analysis of each candidate member of the Hyades Supercluster, using a well-known member of the Hyades cluster as a reference (vB 153). We find that only 28 stars (26 dwarfs and 2 giants) are members, i.e. that 46% of our candidates are members based on the differential abundance analysis. This result confirms that the Hyades Supercluster cannot originate solely from the Hyades cluster.

The impact of a future solar minimum on climate change projections in the Northern Hemisphere

Solar variability represents a source of uncertainty in the future forcings used in climate model simulations. Current knowledge indicates that a descent of solar activity into an extended minimum state is a possible scenario. With aid of experiments from a state-of-the-art Earth system model, we investigate the impact of a future solar minimum on Northern Hemisphere climate change projections. This scenario is constructed from recent 11 year solar-cycle minima of the solar spectral irradiance, and is therefore more conservative than the 'grand' minima employed in some previous modeling studies. Despite the small reduction in total solar irradiance (0.36 W m^-2), relatively large responses emerge in the winter Northern Hemisphere, with a reduction in regional-scale projected warming by up to 40%. To identify the origin of the enhanced regional signals, we assess the role of the different mechanisms by performing additional experiments forced only by irradiance changes at different wavelengths of the solar spectrum. We find that a reduction in visible irradiance drives changes in the stationary wave pattern of the North Pacific and sea-ice cover. A decrease in UV irradiance leads to smaller surface signals, although its regional effects are not negligible. These results point to a distinct but additive role of UV and visible irradiance in the Earth's climate, and stress the need to account for solar forcing as a source of uncertainty in regional scale projections.

Pathways to quiescence: SHARDS view on the star formation histories of massive quiescent galaxies at 1.0 < z < 1.5

We present star formation histories (SFHs) for a sample of 104 massive (stellar mass M > 10^10 M_⊙) quiescent galaxies (MQGs) at z = 1.0–1.5 from the analysis of spectrophotometric data from the Survey for High-z Absorption Red and Dead Sources (SHARDS) and HST/WFC3 G102 and G141 surveys of the GOODS-North field, jointly with broad-band observations from ultraviolet (UV) to far-infrared (far-IR). The sample is constructed on the basis of rest-frame UVJ colours and specific star formation rates (sSFRs = SFR/Mass). The spectral energy distributions (SEDs) of each galaxy are compared to models assuming a delayed exponentially declining SFH. A Monte Carlo algorithm characterizes the degeneracies, which we are able to break taking advantage of the SHARDS data resolution, by measuring indices such as MgUV and D4000. The population of MQGs shows a duality in their properties. The sample is dominated (85 per cent) by galaxies with young mass-weighted ages, t_M t_M < 2 Gyr, short star formation time-scales, 〈τ〉 ∼ 60–200 Myr, and masses log(M/M_⊙) ∼ 10.5. There is an older population (15 per cent) with t_M t_M = 2–4 Gyr, longer star formation time-scales, 〈τ〉∼ 400 Myr, and larger masses, log(M/M_⊙) ∼ 10.7. The SFHs of our MQGs are consistent with the slope and the location of the main sequence of star-forming galaxies at z > 1.0, when our galaxies were 0.5–1.0 Gyr old. According to these SFHs, all the MQGs experienced a luminous infrared galaxy phase that lasts for ∼500 Myr, and half of them an ultraluminous infrared galaxy phase for ∼100 Myr. We find that the MQG population is almost assembled at z ∼ 1, and continues evolving passively with few additions to the population.

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.

EChO. Exoplanet characterisation observatory

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO—the Exoplanet Characterisation Observatory—is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. The use of passive cooling, few moving parts and well established technology gives a low-risk and potentially long-lived mission. EChO will build on observations by Hubble, Spitzer and ground-based telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. However, EChO’s configuration and specifications are designed to study a number of systems in a consistent manner that will eliminate the ambiguities affecting prior observations. EChO will simultaneously observe a broad enough spectral region—from the visible to the mid-infrared—to constrain from one single spectrum the temperature structure of the atmosphere, the abundances of the major carbon and oxygen bearing species, the expected photochemically-produced species and magnetospheric signatures. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules and retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures T_ eq up to 2,000 K, to those of a few Earth masses, with T _eq \u223c 300 K. The list will include planets with no Solar System analog, such as the recently discovered planets GJ1214b, whose density lies between that of terrestrial and gaseous planets, or the rocky-iron planet 55 Cnc e, with day-side temperature close to 3,000 K. As the number of detected exoplanets is growing rapidly each year, and the mass and radius of those detected steadily decreases, the target list will be constantly adjusted to include the most interesting systems. We have baselined a dispersive spectrograph design covering continuously the 0.4–16 μm spectral range in 6 channels (1 in the visible, 5 in the InfraRed), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to \u223c45 K. EChO will be placed in a grand halo orbit around L2. This orbit, in combination with an optimised thermal shield design, provides a highly stable thermal environment and a high degree of visibility of the sky to observe repeatedly several tens of targets over the year. Both the baseline and alternative designs have been evaluated and no critical items with Technology Readiness Level (TRL) less than 4–5 have been identified. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

Magnetic activity and differential rotation in the young Sun-like stars KIC 7985370 and KIC 7765135

Aims. We present a detailed study of the two Sun-like stars KIC 7985370 and KIC 7765135, to determine their activity level, spot distribution, and differential rotation. Both stars were previously discovered by us to be young stars and were observed by the NASA Kepler mission. Methods. The fundamental stellar parameters (vsini, spectral type, T_eff, log g, and [Fe/H]) were derived from optical spectroscopy by comparison with both standard-star and synthetic spectra. The spectra of the targets allowed us to study the chromospheric activity based on the emission in the core of hydrogen Hα and Ca ii infrared triplet (IRT) lines, which was revealed by the subtraction of inactive templates. The high-precision Kepler photometric data spanning over 229 days were then fitted with a robust spot model. Model selection and parameter estimation were performed in a Bayesian manner, using a Markov chain Monte Carlo method. Results. We find that both stars are Sun-like (of G1.5 V spectral type) and have an age of about 100–200 Myr, based on their lithium content and kinematics. Their youth is confirmed by their high level of chromospheric activity, which is comparable to that displayed by the early G-type stars in the Pleiades cluster. The Balmer decrement and flux ratio of their Ca ii-IRT lines suggest that the formation of the core of these lines occurs mainly in optically thick regions that are analogous to solar plages. The spot model applied to the Kepler photometry requires at least seven persistent spots in the case of KIC 7985370 and nine spots in the case of KIC 7765135 to provide a satisfactory fit to the data. The assumption of the longevity of the star spots, whose area is allowed to evolve with time, is at the heart of our spot-modelling approach. On both stars, the surface differential rotation is Sun-like, with the high-latitude spots rotating slower than the low-latitude ones. We found, for both stars, a rather high value of the equator-to-pole differential rotation (dΩ ≈ 0.18 rad d^-1), which disagrees with the predictions of some mean-field models of differential rotation for rapidly rotating stars. Our results agree instead with previous works on solar-type stars and other models that predict a higher latitudinal shear, increasing with equatorial angular velocity, that can vary during the magnetic cycle.

No surviving evolved companions of the progenitor of SN 1006

Type Ia supernovae are thought to occur when a white dwarf made of carbon and oxygen accretes sufficient mass to trigger a thermonuclear explosion(1). The accretion could be slow, from an unevolved (main-sequence) or evolved (subgiant or giant) star(2,3) (the single-degenerate channel), or rapid, as the primary star breaks up a smaller orbiting white dwarf(3,4) (the double-degenerate channel). A companion star will survive the explosion only in the single-degenerate channel(5). Both channels might contribute to the production of type Ia supernovae(6,7), but the relative proportions of their contributions remain a fundamental puzzle in astronomy. Previous searches for remnant companions have revealed one possible case for SN 1572 (refs 8, 9), although that has been questioned(10). More recently, observations have restricted surviving companions to be small, main-sequence stars(11-13), ruling out giant companions but still allowing the single-degenerate channel. Here we report the results of a search for surviving companions of the progenitor of SN 1006 (ref. 14). None of the stars within 4 arc minutes of the apparent site of the explosion is associated with the supernova remnant, and we can firmly exclude all giant and subgiant stars from being companions of the progenitor. In combination with previous results, our findings indicate that fewer than 20 per cent of type Ia supernovae occur through the single-degenerate channel.