SHARDS: survey for high-Z absorption red & dead sources

SHARDS, an ESO/GTC Large Program, is an ultra-deep (26.5 mag) spectro-photometric survey with GTC/OSIRIS designed to select and study massive passively evolving galaxies at z=1.0-2.3 in the GOODS-N field using a set of 24 medium-band filters (FWHM ∼ 17 nm) covering the 500-950 nm spectral range. Our observing strategy has been planned to detect, for z>1 sources, the prominent Mg absorption feature (at rest-frame ∼ 280 nm), a distinctive, necessary, and sufficient feature of evolved stellar populations (older than 0.5 Gyr). These observations are being used to: (1) derive for the first time an unbiased sample of high-z quiescent galaxies, which extends to fainter magnitudes the samples selected with color techniques and spectroscopic surveys; (2) derive accurate ages and stellar masses based on robust measurements of spectral features such as the Mg_UV or D(4000) indices; (3) measure their redshift with an accuracy Δz/(1+z)<0.02; and (4) study emission-line galaxies (starbursts and AGN) up to very high redshifts. The well-sampled optical SEDs provided by SHARDS for all sources in the GOODS-N field are a valuable complement for current and future surveys carried out with other telescopes (e.g., Spitzer, HST, and Herschel).

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.

Spectrophotometric study in the near-ir of a sample of MIPS selected galaxies at Z~2

Our main objective is to determine what kind of galaxies dominate the cosmic SFR density at z~2. Our sample consists of 24 galaxies in Chandra Deep Field South, a unique field for the study of galaxy evolution (12 observed with GNIRS/GEMINI and 12 with ISAAC/VLT). We use H alpha together with the already merged X-ray, ultraviolet, optical, near and mid-infrared imaging data to obtain estimations of SFRs, metallicities, stellar and dynamical masses, AGN activity, and extinction properties. We have obtained 15 Hα detections, 4 rotation curves, and SFR relationship for 7 galaxies. The metallicities obtained for 8 galaxies of the sample are compatible with the metallicities of local galaxies.

Discovery of "Warm Dust" galaxies in clusters at z ~ 0.3: evidence for stripping of cool dust in the dense environment?

Using far-infrared imaging from the "Herschel Lensing Survey," we derive dust properties of spectroscopically confirmed cluster member galaxies within two massive systems at z ~ 0.3: the merging Bullet Cluster and the more relaxed MS2137.3-2353. Most star-forming cluster sources (~90%) have characteristic dust temperatures similar to local field galaxies of comparable infrared (IR) luminosity (T_dust ~ 30 K). Several sub-luminous infrared galaxy (LIRG; L_IR < 10^11 L_☉) Bullet Cluster members are much warmer (T_dust > 37 K) with far-infrared spectral energy distribution (SED) shapes resembling LIRG-type local templates. X-ray and mid-infrared data suggest that obscured active galactic nuclei do not contribute significantly to the infrared flux of these "warm dust" galaxies. Sources of comparable IR luminosity and dust temperature are not observed in the relaxed cluster MS2137, although the significance is too low to speculate on an origin involving recent cluster merging. "Warm dust" galaxies are, however, statistically rarer in field samples (>3σ), indicating that the responsible mechanism may relate to the dense environment. The spatial distribution of these sources is similar to the whole far-infrared bright population, i.e., preferentially located in the cluster periphery, although the galaxy hosts tend toward lower stellar masses (M_* < 10^10 M_☉). We propose dust stripping and heating processes which could be responsible for the unusually warm characteristic dust temperatures. A normal star-forming galaxy would need 30%-50% of its dust removed (preferentially stripped from the outer reaches, where dust is typically cooler) to recover an SED similar to a "warm dust" galaxy. These progenitors would not require a higher IR luminosity or dust mass than the currently observed normal star-forming population.

The relation between cool cluster cores and Herschel-detected star formation in brightest cluster galaxies

We present far-infrared (FIR) analysis of 68 brightest cluster galaxies (BCGs) at 0.08 < z < 1.0. Deriving total infrared luminosities directly from Spitzer and Herschel photometry spanning the peak of the dust component (24-500 μm), we calculate the obscured star formation rate (SFR). 22^+6.2 _–5.3% of the BCGs are detected in the far-infrared, with SFR = 1-150 M ☉ yr^–1. The infrared luminosity is highly correlated with cluster X-ray gas cooling times for cool-core clusters (gas cooling time <1 Gyr), strongly suggesting that the star formation in these BCGs is influenced by the cluster-scale cooling process. The occurrence of the molecular gas tracing Hα emission is also correlated with obscured star formation. For all but the most luminous BCGs (L_TIR > 2 × 10^11 L_☉), only a small (≤0.4 mag) reddening correction is required for SFR(Hα) to agree with SFR_FIR. The relatively low Hα extinction (dust obscuration), compared to values reported for the general star-forming population, lends further weight to an alternate (external) origin for the cold gas. Finally, we use a stacking analysis of non-cool-core clusters to show that the majority of the fuel for star formation in the FIR-bright BCGs is unlikely to originate from normal stellar mass loss.

Infrared Excess sources: Compton thick QSOs, low-luminosity Seyferts or starbursts?

We explore the nature of Infrared Excess sources (IRX), which are proposed as candidates for luminous [L_X(2–10 keV) > 10^43 erg s^−1] Compton thick (NH > 2 × 1024 cm−2) QSOs at z≈ 2. Lower redshift, z≈ 1, analogues of the distant IRX population are identified by first redshifting to z= 2 the spectral energy distributions (SEDs) of all sources with secure spectroscopic redshifts in the AEGIS (6488) and the GOODS-North (1784) surveys and then selecting those that qualify as IRX sources at that redshift. A total of 19 galaxies are selected. The mean redshift of the sample is z≈ 1. We do not find strong evidence for Compton thick QSOs in the sample. For nine sources with X-ray counterparts, the X-ray spectra are consistent with Compton thin active galactic nucleus (AGN). Only three of them show tentative evidence for Compton thick obscuration. The SEDs of the X-ray undetected population are consistent with starburst activity. There is no evidence for a hot dust component at the mid-infrared associated with AGN heated dust. If the X-ray undetected sources host AGN, an upper limit of L_X(2–10 keV) = 10^43 erg s^−1 is estimated for their intrinsic luminosity. We propose that a large fraction of the z≈ 2 IRX population is not Compton thick quasi-stellar objects (QSOs) but low-luminosity [L_X(2–10 keV) < 10^43 erg s^−1], possibly Compton thin, AGN or dusty starbursts. It is shown that the decomposition of the AGN and starburst contribution to the mid-IR is essential for interpreting the nature of this population, as star formation may dominate this wavelength regime.

Improving the identification of high-z Herschel sources with position priors and optical/NIR and FIR/mm photometric redshifts

We present preliminary results about the detection of high redshift (U)LIRGs in the Bullet cluster field by the PACS and SPIRE instruments within the Herschel Lensing Survey (HLS) Program. We describe in detail a photometric procedure designed to recover robust fluxes and deblend faint Herschel sources near the confusion noise. The method is based on the use of the positions of Spitzer/MIPS 24 μm sources as priors. Our catalogs are able to reliably (5σ) recover galaxies with fluxes above 6 and 10 mJy in the PACS 100 and 160 μm channels, respectively, and 12 to 18 mJy in the SPIRE bands. We also obtain spectral energy distributions covering the optical through the far-infrared/millimeter spectral ranges of all the Herschel detected sources, and analyze them to obtain independent estimations of the photometric redshift based on either stellar population or dust emission models. We exemplify the potential of the combined use of Spitzer position priors plus independent optical and IR photometric redshifts to robustly assign optical/NIR counterparts to the sources detected by Herschel and other (sub-)mm instruments.

Local luminous infrared galaxies. I. Spatially resolved observations with the Spitzer infrared spectrograph

We present results from the Spitzer Infrared Spectrograph spectral mapping observations of 15 local luminous infrared galaxies (LIRGs). In this paper, we investigate the spatial variations of the mid-IR emission which includes fine structure lines, molecular hydrogen lines, polycyclic aromatic features (PAHs), continuum emission, and the 9.7 μm silicate feature. We also compare the nuclear and integrated spectra. We find that the star formation takes place in extended regions (several kpc) as probed by the PAH emission, as well as the [Ne II]12.81 μm and [Ne III]15.56 μm emissions. The behavior of the integrated PAH emission and 9.7 μm silicate feature is similar to that of local starburst galaxies. We also find that the minima of the [Ne III]15.56 μm/[Ne II]12.81 μm ratio tends to be located at the nuclei and its value is lower than that of H II regions in our LIRGs and nearby galaxies. It is likely that increased densities in the nuclei of LIRGs are responsible for the smaller nuclear [Ne III]15.56 μm/[Ne II]12.81 μm ratios. This includes the possibility that some of the most massive stars in the nuclei are still embedded in ultracompact H II regions. In a large fraction of our sample, the 11.3 μm PAH emission appears more extended than the dust 5.5 μm continuum emission. We find a dependency of the 11.3 μm PAH/7.7 μm PAH and [Ne II]12.81 μm/11.3 μm PAH ratios with the age of the stellar populations. Smaller and larger ratios, respectively, indicate recent star formation. The estimated warm (300 K

MID-IR luminosities and UV/optical star formation rates at z < 1.4

Ultraviolet (UV) nonionizing continuum and mid-infrared (IR) emission constitute the basis of two widely used star formation (SF) indicators at intermediate and high redshifts. We study 2430 galaxies with z < 1.4 in the Extended Groth Strip with deep MIPS 24 μm observations from FIDEL, spectroscopy from DEEP2, and UV, optical, and near-IR photometry from the AEGIS. The data are coupled with dust-reddened stellar population models and Bayesian spectral energy distribution (SED) fitting to estimate dust-corrected star formation rates (SFRs). In order to probe the dust heating from stellar populations of various ages, the derived SFRs were averaged over various timescales—from 100 Myr for "current" SFR (corresponding to young stars) to 1-3 Gyr for long-timescale SFRs (corresponding to the light-weighted age of the dominant stellar populations). These SED-based UV/optical SFRs are compared to total IR luminosities extrapolated from 24 μm observations, corresponding to 10-18 μm rest frame. The total IR luminosities are in the range of normal star-forming galaxies and luminous IR galaxies (10^10-10^12 L_☉). We show that the IR luminosity can be estimated from the UV and optical photometry to within a factor of 2, implying that most z < 1.4 galaxies are not optically thick. We find that for the blue, actively star-forming galaxies the correlation between the IR luminosity and the UV/optical SFR shows a decrease in scatter when going from shorter to longer SFR-averaging timescales. We interpret this as the greater role of intermediate age stellar populations in heating the dust than what is typically assumed. Equivalently, we observe that the IR luminosity is better correlated with dust-corrected optical luminosity than with dust-corrected UV light. We find that this holds over the entire redshift range. Many so-called green valley galaxies are simply dust-obscured actively star-forming galaxies. However, there exist 24 μm detected galaxies, some with L_IR>10^11 L_☉, yet with little current SF. For them a reasonable amount of dust absorption of stellar light (but presumably higher than in nearby early-type galaxies) is sufficient to produce the observed levels of IR, which includes a large contribution from intermediate and old stellar populations. In our sample, which contains very few ultraluminous IR galaxies, optical and X-ray active galactic nuclei do not contribute on average more than ~50% to the mid-IR luminosity, and we see no evidence for a large population of "IR excess" galaxies.

The host galaxies and black holes of typical z~0.5-1.4 AGNs

We study the stellar and star formation properties of the host galaxies of 58 X-ray-selected AGNs in the GOODS portion of the Chandra Deep Field South (CDF-S) region at z ~ 0.5-1.4. The AGNs are selected such that their rest-frame UV to near-infrared spectral energy distributions (SEDs) are dominated by stellar emission; i.e., they show a prominent 1.6 μm bump, thus minimizing the AGN emission "contamination." This AGN population comprises approximately 50% of the X-ray-selected AGNs at these redshifts. We find that AGNs reside in the most massive galaxies at the redshifts probed here. Their characteristic stellar masses (M_* ~ 7.8 × 10^10 and M_* ~ 1.2 × 10^11 M_☉ at median redshifts of 0.67 and 1.07, respectively) appear to be representative of the X-ray-selected AGN population at these redshifts and are intermediate between those of local type 2 AGNs and high-redshift (z ~ 2) AGNs. The inferred black hole masses (M_BH ~ 2 × 10^8 M_☉) of typical AGNs are similar to those of optically identified quasars at similar redshifts. Since the AGNs in our sample are much less luminous (L_2–10 keV < 10^44 erg s^−1) than quasars, typical AGNs have low Eddington ratios (η ~ 0.01-0.001). This suggests that, at least at intermediate redshifts, the cosmic AGN "downsizing" is due to both a decrease in the characteristic stellar mass of typical host galaxies and less efficient accretion. Finally, there is no strong evidence in AGN host galaxies for either highly suppressed star formation (expected if AGNs played a role in quenching star formation) or elevated star formation when compared to mass-selected (i.e., IRAC-selected) galaxies of similar stellar masses and redshifts.

The star formation and extinction coevolution of UV-selected galaxies over 0.05 < z < 1.2

We use a new stacking technique to obtain mean mid-IR and far-IR to far-UV flux ratios over the rest-frame near-UV, near-IR color-magnitude diagram. We employ COMBO-17 redshifts and COMBO-17 optical, GALEX far- and near-UV, and Spitzer IRAC and MIPS mid-IR photometry. This technique permits us to probe the infrared excess (IRX), the ratio of far-IR to far-UV luminosity, and the specific star formation rate (SSFR) and their coevolution over 2 orders of magnitude of stellar mass and over redshift 0.1 < z < 1.2. We find that the SSFR and the characteristic mass (Script M_0) above which the SSFR drops increase with redshift (downsizing). At any given epoch, the IRX is an increasing function of mass up to Script M_0. Above this mass the IRX falls, suggesting gas exhaustion. In a given mass bin below Script M_0, the IRX increases with time in a fashion consistent with enrichment. We interpret these trends using a simple model with a Schmidt-Kennicutt law and extinction that tracks gas density and enrichment. We find that the average IRX and SSFR follow a galaxy age parameter ξ, which is determined mainly by the galaxy mass and time since formation. We conclude that blue-sequence galaxies have properties which show simple, systematic trends with mass and time such as the steady buildup of heavy elements in the interstellar media of evolving galaxies and the exhaustion of gas in galaxies that are evolving off the blue sequence. The IRX represents a tool for selecting galaxies at various stages of evolution.

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.

Exploitation du potentiel sismique des étoiles naines blanches

Le but de cette thèse est d’explorer le potentiel sismique des étoiles naines blanches pulsantes, et en particulier celles à atmosphères riches en hydrogène, les étoiles ZZ Ceti. La technique d’astérosismologie exploite l’information contenue dans les modes normaux de vibration qui peuvent être excités lors de phases particulières de l’évolution d’une étoile. Ces modes modulent le flux émergent de l’étoile pulsante et se manifestent principalement en termes de variations lumineuses multi-périodiques. L’astérosismologie consiste donc à examiner la luminosité d’étoiles pulsantes en fonction du temps, afin d’en extraire les périodes, les amplitudes apparentes, ainsi que les phases relatives des modes de pulsation détectés, en utilisant des méthodes standards de traitement de signal, telles que des techniques de Fourier. L’étape suivante consiste à comparer les périodes de pulsation observées avec des périodes générées par un modèle stellaire en cherchant l’accord optimal avec un modèle physique reconstituant le plus fidèlement possible l’étoile pulsante. Afin d’assurer une recherche optimale dans l’espace des paramètres, il est nécessaire d’avoir de bons modèles physiques, un algorithme d’optimisation de comparaison de périodes efficace, et une puissance de calcul considérable. Les périodes des modes de pulsation de modèles stellaires de naines blanches peuvent être généralement calculées de manière précise et fiable sur la base de la théorie linéaire des pulsations stellaires dans sa version adiabatique. Afin de définir dans son ensemble un modèle statique de naine blanche propre à l’analyse astérosismologique, il est nécessaire de spécifier la gravité de surface, la température effective, ainsi que différents paramètres décrivant la disposition en couche de l’enveloppe. En utilisant parallèlement les informations obtenues de manière indépendante (température effective et gravité de surface) par la méthode spectroscopique, il devient possible de vérifier la validité de la solution obtenue et de restreindre de manière remarquable l’espace des paramètres. L’exercice astérosismologique, s’il est réussi, mène donc à la détermination précise des paramètres de la structure globale de l’étoile pulsante et fournit de l’information unique sur sa structure interne et l’état de sa phase évolutive. On présente dans cette thèse l’analyse complète réussie, de l’extraction des fréquences à la solution sismique, de quatre étoiles naines blanches pulsantes. Il a été possible de déterminer les paramètres structuraux de ces étoiles et de les comparer remarquablement à toutes les contraintes indépendantes disponibles dans la littérature, mais aussi d’inférer sur la dynamique interne et de reconstruire le profil de rotation interne. Dans un premier temps, on analyse le duo d’étoiles ZZ Ceti, GD 165 et Ross 548, afin de comprendre les différences entre leurs propriétés de pulsation, malgré le fait qu’elles soient des étoiles similaires en tout point, spectroscopiquement parlant. L’analyse sismique révèle des structures internes différentes, et dévoile la sensibilité de certains modes de pulsation à la composition interne du noyau de l’étoile. Afin de palier à cette sensibilité, nouvellement découverte, et de rivaliser avec les données de qualité exceptionnelle que nous fournissent les missions spatiales Kepler et Kepler2, on développe une nouvelle paramétrisation des profils chimiques dans le coeur, et on valide la robustesse de notre technique et de nos modèles par de nombreux tests. Avec en main la nouvelle paramétrisation du noyau, on décroche enfin le ”Saint Graal” de l’astérosismologie, en étant capable de reproduire pour la première fois les périodes observées à la précision des observations, dans le cas de l’étude sismique des étoiles KIC 08626021 et de GD 1212.

Investigating the potential magnetic origin of wind variability in OB stars

In this thesis, the origin of large-scale structures in hot star winds, believed to be responsible for the presence of discrete absorption components (DACs) in the absorption troughs of ultraviolet resonance lines, is constrained using both observations and numerical simulations. These structures are understood as arising from bright regions on the stellar surface, although their physical cause remains unknown. First, we use high quality circular spectropolarimetric observations of 13 well-studied OB stars to evaluate the potential role of dipolar magnetic fields in producing DACs. We perform longitudinal field measurements and place limits on the field strength using Bayesian inference, assuming that it is dipolar. No magnetic field was detected within this sample. The derived constraints statistically refute any significant dynamical influence from a magnetic dipole on the wind for all of these stars, ruling out such fields as a cause for DACs. Second, we perform numerical simulations using bright spots constrained by broadband optical photometric observations. We calculate hydrodynamical wind models using three sets of spot sizes and strengths. Co-rotating interaction regions are yielded in each model, and radiative transfer shows that the properties of the variations in the UV resonance lines synthesized from these models are consistent with those found in observed UV spectra, establishing the first consistent link between UV spectroscopic line profile variability and photometric variations and thus supporting the bright spot paradigm (BSP). Finally, we develop and apply a phenomenological model to quantify the measurable effects co-rotating bright spots would have on broadband optical photometry and on the profiles of photopheric lines in optical spectra. This model can be used to evaluate the existence of these spots, and, in the event of their detection, characterize them. Furthermore, a tentative spot evolution model is presented. A preliminary analysis of its output, compared to the observed photometric variations of xi Persei, suggests the possible existence of “active longitudes” on the surface of this star. Future work will expand the range of observational diagnostics that can be interpreted within the BSP, and link phenomenology (bright spots) to physical processes (magnetic spots or non-radial pulsations).

The Rotational Evolution and Magnetospheric Emission of the Magnetic Early B-type Stars

How do the magnetic fields of massive stars evolve over time? Are their gyrochronological ages consistent with ages inferred from evolutionary tracks? Why do most stars predicted to host Centrifugal Magnetospheres (CMs) display no H$\alpha$ emission? Does plasma escape from CMs via centrifugal breakout events, or by a steady-state leakage mechanism? This thesis investigates these questions via a population study with a sample of 51 magnetic early B-type stars. The longitudinal magnetic field \bz~was measured from Least Squares Deconvolution profiles extracted from high-resolution spectropolarimetric data. New rotational periods $P_{\rm rot}$ were determined for 15 stars from \bz, leaving only 3 stars for which $P_{\rm rot}$ is unknown. Projected rotational velocities \vsini~were measured from multiple spectral lines. Effective temperatures and surface gravities were measured via ionization balances and line profile fitting of H Balmer lines. Fundamental physical parameters, \bz, \vsini, and $P_{\rm rot}$ were then used to determine radii, masses, ages, dipole oblique rotator model, stellar wind, magnetospheric, and spindown parameters using a Monte Carlo approach that self-consistently calculates all parameters while accounting for all available constraints on stellar properties. Dipole magnetic field strengths $B_{\rm d}$ follow a log-normal distribution similar to that of Ap stars, and decline over time in a fashion consistent with the expected conservation of fossil magnetic flux. $P_{\rm rot}$ increases with fractional main sequence age, mass, and $B_{\rm d}$, as expected from magnetospheric braking. However, comparison of evolutionary track ages to maximum spindown ages $t_{\rm S,max}$ shows that initial rotation fractions may be far below critical for stars with $M_*>10 M_\odot$. Computing $t_{\rm S,max}$ with different mass-loss prescriptions indicates that the mass-loss rates of B-type stars are likely much lower than expected from extrapolation from O-type stars. Stars with H$\alpha$ in emission and absorption occupy distinct regions in the updated rotation-magnetic confinement diagram: H$\alpha$-bright stars are found to be younger, more rapidly rotating, and more strongly magnetized than the general population. Emission strength is sensitive both to the volume of the CM and to the mass-loss rate, favouring leakage over centrifugal breakout.