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.

An emission-line survey of the Milky Way /

Includes bibliographical references

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.

Non-LTE Monte Carlo radiative transfer. II. Nonisothermal solutions for viscous Keplerian disks

We discuss the basic hydrodynamics that determines the density structure of the disks around hot stars. Observational evidence supports the idea that these disks are Keplerian (rotationally supported) gaseous disks. A popular scenario in the literature, which naturally leads to the formation of Keplerian disks, is the viscous decretion model. According to this scenario, the disks are hydrostatically supported in the vertical direction, while the radial structure is governed by the viscous transport. This suggests that the temperature is one primary factor that governs the disk density structure. In a previous study we demonstrated, using three-dimensional non-LTE Monte Carlo simulations, that viscous Keplerian disks can be highly nonisothermal. In this paper we build on our previous work and solve the full problem of the steady state nonisothermal viscous diffusion and vertical hydrostatic equilibrium. We find that the self-consistent solution departs significantly from the analytic isothermal density, with potentially large effects on the emergent spectrum. This implies that nonisothermal disk models must be used for a detailed modeling of Be star disks.

TOWARD UNDERSTANDING THE B[e] PHENOMENON. IV. MODELING OF IRAS 00470+6429

FS CMa type stars are a recently described group of objects with the B[e] phenomenon which exhibits strong emission-line spectra and strong IR excesses. In this paper, we report the first attempt for a detailed modeling of IRAS 00470+6429, for which we have the best set of observations. Our modeling is based on two key assumptions: the star has a main-sequence luminosity for its spectral type (B2) and the circumstellar (CS) envelope is bimodal, composed of a slowly outflowing disklike wind and a fast polar wind. Both outflows are assumed to be purely radial. We adopt a novel approach to describe the dust formation site in the wind that employs timescale arguments for grain condensation and a self-consistent solution for the dust destruction surface. With the above assumptions we were able to satisfactorily reproduce many observational properties of IRAS 00470+6429, including the Hi line profiles and the overall shape of the spectral energy distribution. Our adopted recipe for dust formation proved successful in reproducing the correct amount of dust formed in the CS envelope. Possible shortcomings of our model, as well as suggestions for future improvements, are discussed.

TOWARD UNDERSTANDING THE B[e] PHENOMENON. III. PROPERTIES OF THE OPTICAL COUNTERPART OF IRAS 00470+6429

FS CMa type stars are a group of Galactic objects with the B[e] phenomenon. They exhibit strong emission-line spectra and infrared excesses, which are most likely due to recently formed circumstellar dust. The group content and identification criteria were described in the first two papers of the series. In this paper we report our spectroscopic and photometric observations of the optical counterpart of IRAS 00470+6429 obtained in 2003-2008. The optical spectrum is dominated by emission lines, most of which have P Cyg type profiles. We detected significant brightness variations, which may include a regular component, and variable spectral line profiles in both shape and position. The presence of a weak Li I 6708 angstrom line in the spectrum suggests that the object is most likely a binary system with a B2-B3 spectral-type primary companion of a luminosity log L/L(circle dot) = 3.9 +/- 0.3 and a late-type secondary companion. We estimate a distance toward the object to be 2.0 +/- 0.3 kpc from the Sun.

A cyclic bipolar wind in the interacting binary V 393 Scorpii

V393 Scorpii is a double periodic variable characterized by a relatively stable non-orbital photometric cycle of 253 d. Mennickent et al. argue for the presence of a massive optically thick disc around the more massive B-type component and describe the evolutionary stage of the system. In this paper, we analyse the behaviour of the main spectroscopic optical lines during the long non-orbital photometric cycle. We study the radial velocity of the donor determining its orbital elements and find a small but significant orbital eccentricity (e = 0.04). The donor spectral features are modelled and removed from the spectrum at every observing epoch using the light-curve model given by Mennickent et al. We find that the line emission is larger during eclipses and mostly comes from a bipolar wind. We also find that the long cycle is explained in terms of a modulation of the wind strength; the wind has a larger line and continuum emissivity at the high state. We report the discovery of highly variable chromospheric emission in the donor, as revealed by the Doppler maps of the emission lines Mg II 4481 and C I 6588. We discuss notable and some novel spectroscopic features like discrete absorption components, especially visible at blue depressed O I 7773 absorption wings during the second half-cycle, Balmer double emission with V/R curves showing 'Z-type' and 'S-type' excursions around secondary and main eclipses, respectively, and H beta emission wings extending up to +/- 2000 km s(-1). We also discuss possible causes for these phenomena and for their modulations with the long cycle.

Probing the properties of Be star discs with spectroastrometry and NLTE radiative transfer modelling: beta CMi

While the presence of discs around classical Be stars is well established, their origin is still uncertain. To understand what processes result in the creation of these discs and how angular momentum is transported within them, their physical properties must be constrained. This requires comparing high spatial and spectral resolution data with detailed radiative transfer modelling. We present a high spectral resolution, R similar to 80 000, sub-milliarcsecond precision, spectroastrometric study of the circumstellar disc around the Be star beta CMi. The data are confronted with 3D, non-local thermodynamic equilibrium radiative transfer calculations to directly constrain the properties of the disc. Furthermore, we compare the data to disc models featuring two velocity laws: Keplerian, the prediction of the viscous disc model, and angular momentum conserving rotation. It is shown that the observations of beta CMi can only be reproduced using Keplerian rotation. The agreement between the model and the observed spectral energy distribution, polarization and spectroastrometric signature of beta CMi confirms that the discs around Be stars are well modelled as viscous decretion discs.

THE FIRST DETERMINATION OF THE VISCOSITY PARAMETER IN THE CIRCUMSTELLAR DISK OF A Be STAR

Be stars possess gaseous circumstellar decretion disks, which are well described using standard alpha-disk theory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous decretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter alpha = 1.0 +/- 0.2, corresponding to a mass injection rate (M) over dot = (3.5 +/- 1.3) x 10(-8) M-circle dot yr(-1). Such a large value of a suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass-loss rate inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, but some other mechanism.

Stochastic gravito-inertial modes discovered by CoRoT in the hot Be star HD51452

Context. Be stars are rapidly rotating stars with a circumstellar decretion disk. They usually undergo pressure and/or gravity pulsation modes excited by the kappa-mechanism, i.e. an effect of the opacity of iron-peak elements in the envelope of the star. In the Milky Way, p-modes are observed in stars that are hotter than or equal to the B3 spectral type, while g-modes are observed at the B2 spectral type and cooler. Aims. We observed a B0IVe star, HD51452, with the high-precision, high-cadence photometric CoRoT satellite and high-resolution, ground-based HARPS and SOPHIE spectrographs to study its pulsations in great detail. We also used the lower resolution spectra available in the BeSS database. Methods. We analyzed the CoRoT and spectroscopic data with several methods: CLEAN-NG, FREQFIND, and a sliding window method. We also analyzed spectral quantities, such as the violet over red (V/R) emission variations, to obtain information about the variation in the circumstellar environment. We calculated a stellar structure model with the ESTER code to test the various interpretation of the results. Results. We detect 189 frequencies of variations in the CoRoT light curve in the range between 0 and 4.5 c d(-1). The main frequencies are also recovered in the spectroscopic data. In particular we find that HD51452 undergoes gravito-inertial modes that are not in the domain of those excited by the kappa-mechanism. We propose that these are stochastic modes excited in the convective zones and that at least some of them are a multiplet of r-modes (i.e. subinertial modes mainly driven by the Coriolis acceleration). Stochastically excited gravito-inertial modes had never been observed in any star, and theory predicted that their very low amplitudes would be undetectable even with CoRoT. We suggest that the amplitudes are enhanced in HD51452 because of the very rapid stellar rotation. In addition, we find that the amplitude variations of these modes are related to the occurrence of minor outbursts. Conclusions. Thanks to CoRoT data, we have detected a new kind of pulsations in HD51452, which are stochastically excited gravito-inertial modes, probably due to its very rapid rotation. These modes are probably also present in other rapidly rotating hot Be stars.

New activity in the large circumstellar disk of the Be-shell star 48 Librae

Aims. Spectroscopic, polarimetric, and high spectral resolution interferometric data covering the period 1995-2011 are analyzed to document the transition into a new phase of circumstellar disk activity in the classical Be-shell star 48 Lib. The objective is to use this broad data set to additionally test disk oscillations as the basic underlying dynamical process. Methods. The long-term disk evolution is described using the V/R ratio of the violet and red emission components of H alpha and Br gamma, radial velocities and profiles of He I and optical metal shell lines, as well as multi-band BVRI polarimetry. Single-epoch broad-band and high-resolution interferometric visibilities and phases are discussed with respect to a classical disk model and the given baseline orientations. Results. Spectroscopic signatures of disk asymmetries in 48 Lib vanished in the late nineties but recovered some time between 2004 and 2007, as shown by a new large-amplitude and long-duration V/R cycle. Variations in the radial velocity and line profile of conventional shell lines correlate with the V/R behavior. They are shared by narrow absorption cores superimposed on otherwise seemingly photospheric He I lines, which may form in high-density gas at the inner disk close to the photosphere. Large radial velocity variations continued also during the V/R-quiet years, suggesting that V/R may not always be a good indicator of global density waves in the disk. The comparison of the polarization after the recovery of the V/R activity shows a slight increase, while the polarization angle has been constant for more than 20 years, placing tight limits on any 3-D precession or warping of the disk. The broad H-band interferometry gives a disk diameter of (1.72 +/- 0.2) mas (equivalent to 15 stellar radii), position angle of the disk (50 +/- 9)degrees and a relatively low disk flattening of 1.66 +/- 0.3. Within the errors the same disk position angle is derived from polarimetric observations and from photocenter shifts across Br gamma. The high-resolution interferometric visibility and phase profiles show a double or even multiple-component structure. A preliminary estimate based on the size of the Br gamma emitting region indicates a large diameter for the disk (tens of stellar radii). Overall, no serious contradiction between the observations and the disk-oscillation model could be construed.

Stochastic gravito-inertial modes discovered by CoRoT in the hot Be star HD51452

Context. Be stars are rapidly rotating stars with a circumstellar decretion disk. They usually undergo pressure and/or gravity pulsation modes excited by the κ-mechanism, i.e. an effect of the opacity of iron-peak elements in the envelope of the star. In the Milky Way, p-modes are observed in stars that are hotter than or equal to the B3 spectral type, while g-modes are observed at the B2 spectral type and cooler. Aims. We observed a B0IVe star, HD51452, with the high-precision, high-cadence photometric CoRoT satellite and high-resolution, ground-based HARPS and SOPHIE spectrographs to study its pulsations in great detail. We also used the lower resolution spectra available in the BeSS database. Methods. We analyzed the CoRoT and spectroscopic data with several methods: Clean-NG, FreqFind, and a sliding window method. We also analyzed spectral quantities, such as the violet over red (V/R) emission variations, to obtain information about the variation in the circumstellar environment. We calculated a stellar structure model with the ESTER code to test the various interpretation of the results. Results. We detect 189 frequencies of variations in the CoRoT light curve in the range between 0 and 4.5 c d−1. The main frequencies are also recovered in the spectroscopic data. In particular we find that HD51452 undergoes gravito-inertial modes that are not in the domain of those excited by the κ-mechanism. We propose that these are stochastic modes excited in the convective zones and that at least some of them are a multiplet of r-modes (i.e. subinertial modes mainly driven by the Coriolis acceleration). Stochastically excited gravito-inertial modes had never been observed in any star, and theory predicted that their very low amplitudes would be undetectable even with CoRoT. We suggest that the amplitudes are enhanced in HD51452 because of the very rapid stellar rotation. In addition, we find that the amplitude variations of these modes are related to the occurrence of minor outbursts. Conclusions. Thanks to CoRoT data, we have detected a new kind of pulsations in HD51452, which are stochastically excited gravito-inertial modes, probably due to its very rapid rotation. These modes are probably also present in other rapidly rotating hot Be stars.

Astrophysical parameters of LS 2883 and implications for the PSR B1259-63 gamma-ray binary

Only a few binary systems with compact objects display TeV emission. The physical properties of the companion stars represent basic input for understanding the physical mechanisms behind the particle acceleration, emission, and absorption processes in these so-called gamma-ray binaries. Here we present high-resolution and high signal-to-noise optical spectra of LS 2883, the Be star forming a gamma-ray binary with the young non-accreting pulsar PSR B1259-63, showing it to rotate faster and be significantly earlier and more luminous than previously thought. Analysis of the interstellar lines suggests that the system is located at the same distance as (and thus is likely a member of) Cen OB1. Taking the distance to the association, d = 2.3 kpc, and a color excess of E(B – V) = 0.85 for LS 2883 results in MV ≈ –4.4. Because of fast rotation, LS 2883 is oblate (R eq sime 9.7 R ☉ and R pole sime 8.1 R ☉) and presents a temperature gradient (T eq≈ 27,500 K, log g eq = 3.7; T pole≈ 34,000 K, log g pole = 4.1). If the star did not rotate, it would have parameters corresponding to a late O-type star. We estimate its luminosity at log(L */L ☉) sime 4.79 and its mass at M * ≈ 30 M ☉. The mass function then implies an inclination of the binary system i orb ≈ 23°, slightly smaller than previous estimates. We discuss the implications of these new astrophysical parameters of LS 2883 for the production of high-energy and very high-energy gamma rays in the PSR B1259-63/LS 2883 gamma-ray binary system. In particular, the stellar properties are very important for prediction of the line-like bulk Comptonization component from the unshocked ultrarelativistic pulsar wind.

NGC 7419 as a template for red supergiant clusters

Context. The open cluster NGC 7419 is known to contain five red supergiants and a very high number of Be stars. However, there are conflicting reports about its age and distance that prevent a useful comparison with other clusters. Aims. We intend to obtain more accurate parameters for NGC 7419, using techniques different from those of previous authors, so that it may be used as a calibrator for more obscured clusters. Methods. We obtained Strömgren photometry of the open cluster NGC 7419, as well as classification spectroscopy of ~20 stars in the area. We then applied standard analysis and classification techniques. Results. We find a distance of 4 ± 0.4 kpc and an age of 14 ± 2 Myr for NGC 7419. The main-sequence turn-off is found at spectral type B1, in excellent agreement. We identify 179 B-type members, implying that there are more than 1200 M⊙ in B stars at present. Extrapolating this to lower masses indicates an initial cluster mass of between 7000 and 10 000 M⊙, depending on the shape of the initial mass function. We find a very high fraction (≈40%) of Be stars around the turn-off, but very few Be stars at lower masses. We also report for the first time a strong variability in the emission characteristics of Be stars. We verified that the parameters of the red supergiant members can be used to obtain accurate cluster parameters. Conclusions. NGC 7419 is sufficiently massive to serve as a testbed for theoretical predictions and as a template to compare more obscured clusters. The distribution of stars above the main-sequence turn-off is difficult to accommodate with current evolutionary tracks. Though the presence of five red supergiants is marginally consistent with theoretical expectations, the high number of Be stars and very low number of luminous evolved B stars hint at some unknown physical factor that is not considered in current synthesis models.