Dissecting galaxies with quantitative spectroscopy of the brightest stars in the Universe

Measuring distances to galaxies, determining their chemical composition, investigating the nature of their stellar populations and the absorbing properties of their interstellar medium are fundamental activities in modern extragalactic astronomy helping to understand the evolution of galaxies and the expanding universe. The optically brightest stars in the universe, blue supergiants of spectral A and B, are unique tools for these purposes. With absolute visual magnitudes up to M_V ≃ ‐9.5 they are ideal to obtain accurate quantitative information about galaxies through the powerful modern methods of quantitative stellar spectroscopy. The spectral analysis of individual blue supergiant targets provides invaluable information about chemical abundances and abundance gradients, which is more comprehensive than the one obtained from H_II regions, as it includes additional atomic species, and which is also more accurate, since it avoids the systematic uncertainties inherent in the strong line studies usually applied to the H_II regions of spiral galaxies beyond the Local Group. Simultaneously, the spectral analysis yields stellar parameters and interstellar extinction for each individual supergiant target, which provides an alternative very accurate way to determine extragalactic distances through a newly developed method, called the Flux‐weighted Gravity–Luminosity Relationship (FGLR). With the present generation of 10 m‐class telescopes these spectroscopic studies can reach out to distances of 10 Mpc. The new generation of 30 m‐class telescopes will allow to extend this work out to 30 Mpc, a substantial volume of the local universe (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Diffuse Interstellar Band at 862 nm as a Reddening Tracer for GAIA

We investigate the plausibility of using diffuse interstellar band at862 nm for tracing interstellar extinction with the ESA's astrometric space mission GAIA. For this purpose we perform numerical tests to simulate the conditions of real observations, covering a wide range of stellar parameters and different amounts of interstellar extinction. Our simulations indicate that with the present Radial Velocity Spectrometer setup the uncertainty in color excess of σ_E(B-V)≤ 0.05 can be achieved only for the interstellar reddening tracers brighter than V ∼ 13. None of the plausible tracers can provide accurate color excesses (σ _E(B-V) ≤ 0.05) at the distances beyond 2 kpc. We therefore conclude that with the currently planned instrumentation onboard GAIA this method can not be used as a stand-alone approach for probing interstellar extinction on the Galactic distance scales within the framework of the GAIA mission. This revised version was published online in July 2006 with corrections to the Cover Date.     

Globular Clusters in the Large Magellanic Cloud: An Impact from GAIA Photometry

The goal of this work is to assess the expected scientific output from the photometric studies of globular clusters in the Large Magellanic Cloud with ESA's astrometric space mission GAIA. For this purpose we simulate GAIA photometry of individual stars in synthetic cluster populations, covering a large range of cluster ages and metallicities. We find that accurate effective temperatures (Δ T _eff<10%) can be obtained from GAIA photometry down to V ∼ 18 for stars in populations within the studied metallicity range ([M/H] = -0.4 ... -1.7). GAIA will also provide photometric metallicities (Δ [M/H] ≲ 0.3 dex) for the cluster giants brighter than V ∼ 17.5. The knowledge of the effective temperature sand metallicities will allow to obtain accurate ages of stellar populations younger than about 1 Gyr using the usual procedure of main sequence turn-off point fitting. Ages of older stellar populations (≳ 1 Gyr) may be constrained from the isochrone fits to the giant branches in the observed CMDs. We conclude that GAIA will provide excellent opportunities for studying star formation histories far beyond the Milky Way, providing means for better understanding of stellar and galactic evolution in different astrophysical environments. This revised version was published online in July 2006 with corrections to the Cover Date.     

Calculated atmospheric color refraction and observed stellar positions

Astrometric observations at different zenith distances have been performed in Dresden in an area centered atNGC 6791 where there are some stars with reliable color information (widely dispersed spectral types in the MK systemand color indices B_T – V_T) as well as with accurate positions from Tycho‐2 catalog. The results are used to estimate how significant improvements in stellar positions may be when accurate corrections for color refraction are taken into account. We have treated two cases for refraction calculations: (1) a photometric case for color indices and (2) a spectral case for spectral types and luminosity classes. To calculate refraction we use Stone's modified computer code (Malyuto & Meinel 2000). To treat the photometric case we have calculated the synthetic color indices for the spectral energy distributions of Sviderskiene (1988). The positional improvements due to including color refraction corrections are significant and slightly larger in the spectral case. An improvement of about 15% is reached at a zenith distance of 65°. Our basic conclusion is that color refraction should be taken into account for obtaining accurate stellar positions from ground based observations at larger zenith distances. Reliable refraction corrections may be calculated from spectral and/or photometric data.     

The ratio of total to selective absorption in the Carina OB2 association

Several previous studies in and around the Great Carina nebula (NGC 3372) have strongly indicated an abnormal interstellar extinction law with a high ratio of total to selective absorptionR. In the present study, newUBV photometric data and accurate MK spectral-types of stars in the region of the Carina OB2 association are used to show that (1) Car OB2 is a genuine stellar association located at a distance of 3.2 kpc, and (2) the interstellar extinction law seems to be normal throughout most of the region. A representative value ofR = A _v/E(B – V) = 3.0 is derived for Car OB2 from the variable extinction method.Visiting Astronomer, Cerro Tololo Inter-American Observatory, operated by AURA, Inc. under contract with the National Science Foundation.     

3-D Structure of the Galactic Interstellar Matter: A Contribution from GAIA

We investigate the possibilities for tracing interstellar extinction with the ESA's astrometric space mission GAIA. The analysis is based on detailed simulations of the GAIA photometry, which are used to derive the distribution of interstellar matter in a modelled Galaxy. We find that `small' diffuse clouds (diameter D = 4 pc, E <sub>B-V</sub> = 0.06) will be easily traced with GAIA up to the distances of ∼ 800 pc. `Large' diffuse interstellar clouds (D = 10 pc, E _B-V = 0.13) will be located up to the distances of ∼ 2.5 kpc. This holds for the reddening tracers of spectral types O – K2 brighter than V = 17. Inmost cases, due to their low spatial density, the early type stars (O– A2) cannot provide reliable information about the distribution of interstellar matter. None of the reddening tracers measured by GAIA will provide reliable identification of the individual interstellar clouds beyond the distances of ∼ 3 kpc. Therefore, we conclude that the information available from photometric observations will be not sufficient for the detailed reconstruction of the 3-D distribution of Galactic interstellar matter. It is therefore extremely important to define the new strategies which would allow to combine all the available information, including the earlier space- and/or ground-based investigations, together with the information which will be provided by GAIA itself (parallaxes, E _B-V etc.). This revised version was published online in July 2006 with corrections to the Cover Date.     

Asymptotic Giant Branch Stars as Tracers of Star Formation Histories: the GAIA Context

We discuss a possible use of the asymptotic giant branch (AGB) stars for tracing star formation histories on the Galactic and extragalactic distance scales with the ESA's astrometric space mission GAIA. Extensive numerical simulations demonstrate that metallicities (Δ [M/H] ≲ 0.3) can be obtained for the AGB stars with GAIA up to the distances of ∼ 200 kpc, if no interstellar extinction is present. Reliable population ages can be also obtained from the AGB stars if their T _eff are constrained precisely. We show that precise effective temperatures can be obtained by fitting observed spectral energy distributions of the AGB stars with theoretical fluxes calculated from the synthetic spectra. A combination of the derived effective temperatures with the bolometric luminosities allows to derive precise population ages for a wide range of ages and metallicities over the large distance scales. This demonstrates that AGB stars can be employed very effectively for tracing star formation histories with GAIA, allowing to refine the global evolutionary scenarios of stellar populations in the Milky Way and the galaxies beyond. This revised version was published online in July 2006 with corrections to the Cover Date.     

Near-IR spectral evolution of dusty starburst galaxies

We propose a multicomponent analysis of starburst galaxies, based on a model that takes into account the young and evolved stellar components and the gas emission, with their respective extinction, in the frame of a coherent dust distribution pattern. Near-IR signatures are preferentially investigated, in order to penetrate as deep as possible into the dusty starburst cores. We computed the 1.4-2.5 μm spectra of synthetic stellar populations evolving through strong, short timescale bursts of star formation (continuum and lines, R ? 500). The evolution model is specifically sensitive to cool stellar populations (AGB and red supergiant stars). It takes advantage of the stellar library of Lançon & Rocca-Volmerange (1992) [A&ASS, 96, 593], observed with the same instrument (FTS/CFHT) as the analysed galaxy sample, so that the instrumental effects are minimised. The main near-IR observable constraints are the molecular signatures of CO and H₂O and the slope of the continuum, observed over a range exceptionally broad for spectroscopic data. The H - K colour determined from the spectra measures the intrinsic stellar energy distribution but also differential extinction, which is further constrained by optical emission line ratios. Other observational constraints are the near-IR emission lines (Brγ, He I 2.06 μm, [Fe II] 1.64 μm, H₂ 2.12 μm) and the far-IR luminosity. The coherence of the results relies on the interpretation in terms of stellar populations from which all observable properties are derived, so that the link between the various wavelength ranges is secured. The luminosity L_K is used for the absolute calibration.We apply this approach to the typical spectrum of the core of NGC 1614. Consistent solutions for the starburst characteristics (star-formation rate, IMF, burst age, morphology) are found and the role of each observational constraint in deriving satisfactory models is extensively discussed. The acceptable contamination of the K band light by the underlying population amounts ≥ 15% even through a 5 arcsec aperture. The model leads to a limit on the direct absorption of Lyman continuum photons by dust situated inside the ionised areas, which in turn, with standard gas-to-dust ratios, translates into small characteristic sizes for the individual coexisting H II regions of the massive starburst area (clusters containing ∼ 10² ionising stars). We show that room is left for IMFs extending to 120 M_⊙, rather than truncated at ∼ 60 M_⊙ as most conservative studies conclude. High internal velocity dispersions (≥ 20 km s⁻¹) are then needed for the H II regions. An original feature of this work is to base the analysis of near-infrared spectral galaxy observations on a large wavelength range, using models constructed with spectral stellar data observed with the same instrument. However a broader use of this spectral evolution model on other spectral or photometric data samples is possible if the spectral resolution of the model is adapted to observations or if colours are derived from the energy distributions.Catherine J. Cesarsky     

Two new LBV candidates in the M 33 galaxy

We present two new luminous blue variable (LBV) candidate stars discovered in the M33 galaxy. We identified these stars as massive star candidates at the final stages of evolution, presumably with a notable interstellar extinction. The candidates were selected from the Massey et al. catalog based on the following criteria: emission in H _α , V<18^./m 5 and 0_.^m 35 < (B - V) < 1_.^m 2. The spectra of both stars reveal a broad and strong H _α emission with extended wings (770 and 1000 kms⁻¹). Based on the spectra we estimated the main parameters of the stars. Object N45901 has a bolometric luminosity log(L/L_⊙) = 6.0–6.2 with the value of interstellar extinction A _V = 2.3 ± 0.1. The temperature of the star’s photosphere is estimated as T⋆ ∼ 13000–15000 K, its probable mass on the Zero Age Main Sequence is M∼ 60–80 M_⊙. The infrared excess in N 45901 corresponds to the emission of warm dust with the temperature Twarm ∼ 1000 K, and amounts to 0.1%of the bolometric luminosity. A comparison of stellar magnitude estimates from different catalogs points to the probable variability of the object N45901. Bolometric luminosity of the second object, N125093, is log(L/L_⊙) = 6.3 − 6.6, the value of interstellar extinction is A _V = 2.75 ± 0.15. We estimate its photosphere’s temperature as T⋆∼ 13000–16000K, the initial mass as M ∼ 90–120 M_⊙. The infrared excess in N125093 amounts to 5–6% of the bolometric luminosity. Its spectral energy distribution reveals two thermal components with the temperatures T_warm ∼ 1000K and T_cold ∼ 480 K. The [Ca II] λλ7291, 7323 lines, observed in LBV-like stars Var A and N93351 in M33 are also present in the spectrum of N 125093. These lines indicate relatively recent gas eruptions and dust activity linked with them. High bolometric luminosity of these stars and broad H α emissions allow classifying the studied objects as LBV candidates.     

CDM, Feedback and the Hubble Sequence

We have performed TreeSPH simulations of galaxy formation in a standard ΛCDM cosmology, including effects of star formation, energetic stellar feedback processes and a meta-galactic UV field, and obtain a mix of disk, lenticular and elliptical galaxies. The disk galaxies are deficient in angular momentum by only about a factor of two compared to observed disk galaxies. The stellar disks have approximately exponential surface density profiles, and those of the bulges range from exponential to r ^1/4, as observed. The bulge-to-disk ratios of the disk galaxies are consistent with observations and likewise are their integrated B-V colours, which have been calculated using stellar population synthesis techniques. Furthermore, we can match the observed I-band Tully-Fisher (TF) relation, provided that the mass-to-light ratio of disk galaxies is (M/L _I) ≃ 0.6–0.7. The ellipticals and lenticulars have approximately r ^1/4 stellar surface density profiles, are dominated by non-disklike kinematics and flattened due to non-isotropic stellar velocity distributions, again consistent with observations. This revised version was published online in September 2006 with corrections to the Cover Date.     

Surface brightness fluctuations as primary and secondary distance indicators

The surface brightness fluctuations (SBF) method measures the variance in a galaxy’s light distribution arising from fluctuations in the numbers and luminosities of individual stars per resolution element. Once calibrated for stellar population effects, SBF measurements with HST provide distances to early-type galaxies with unrivaled precision. Optical SBF data from HST for the Virgo and Fornax clusters give the relative distances of these nearby fiducial clusters with 2% precision and constrain their internal structures. Observations in hand will allow us to tie the Coma cluster, the standard of comparison for distant cluster studies, into the same precise relative distance scale. The SBF method can be calibrated in an absolute sense either empirically from Cepheids or theoretically from stellar population models. The agreement between the model and empirical zero points has improved dramatically, providing an independent confirmation of the Cepheid distance scale. SBF is still brighter in the near-IR, and an ongoing program to calibrate the method for the F110W and F160W passbands of the Wide Field Camera 3 IR channel will enable accurate distance derivation whenever a large early-type galaxy or bulge is observed in these passbands at distances reaching well out into the Hubble flow.     

WR 38/38a and the ratio of total-to-selective extinction in Carina

A reanalysis of the (seemingly very distant) open cluster Shorlin 1, the group of stars associated with WR 38 and WR 38a, is made on the basis of existing UBV and JHK _s observations for cluster members. The 2MASS observations, in particular, imply a mean cluster reddening of E _B−V =1.45±0.07 and a distance of 2.94±0.12 kpc. The reddening agrees with the UBV results provided that the local reddening slope is described by E _U−B /E _B−V =0.64±0.01, but the distance estimates in the 2MASS and UBV systems agree only if the ratio of total-to-selective extinction for the associated dust is R=A _V /E _B−V =4.0±0.1. Both results are similar to what has been obtained for adjacent clusters in the Eta Carinae region by similar analyses, which suggests that “anomalous” dust extinction is widespread through the region, particularly for groups reddened by relatively nearby dust. Dust associated with the Eta Carinae complex itself appears to exhibit more “normal” qualities. The results have direct implications for the interpretation of distances to optical spiral arm indicators for the Galaxy at ℓ=287°–291°, in particular the Carina arm here is probably little more than ∼2 kpc distant, rather than 2.5–3 kpc distant as implied in previous studies. Newly-derived intrinsic parameters for the two cluster Wolf-Rayet stars WR 38 (WC4) and WR 38a (WN5) are in good agreement with what is found for other WR stars in Galactic open clusters, which was not the case previously.     

Photometric variability and evolutionary status of the supergiant with an infrared excess HD 179821=V1427 aquilae

We present new results of our UBV photometry for HD 179821=V1427 Aql, an F supergiant with an infrared excess, from 2000 to 2008. The semiregular low-amplitude (ΔV = 0_. ^m 05−0_. ^m 20) photometric variability of the star with a cycle period from 130 to 200 days is caused by pulsations, along with the instability of a variable stellar wind. V1427 Aql also exhibits a long-term trend in the brightness and colors that is probably attributable to a change in the stellar temperature as a result of mass loss episodes, which cause variations in the continuum formation level. We present the results of our JHKLM photometry for V1427 Aql in 1992–2008. We trace the trend in the near-infrared brightness, which agrees with the long-term variability in the V band. Based on broadband photometry, we have determined the color excess for V1427 Aql: E(B−V) = 0.7. Based on low-resolution spectroscopy, we have estimated the stellar temperature and revealed variability of the Hα line caused by a change in the contribution from the emission component. The hypotheses of whether the star belongs to post-AGB objects or to massive yellow hypergiants are discussed.     

The Tully-Fisher relations of the Eridanus group of galaxies

The Tully-Fisher (TF) or the luminosity-linewidth relations of the galaxies in the Eridanus group are constructed using the HI rotation curves and the luminosities in the optical and in the near-infrared bands. The slopes of the TF relations (absolute magnitudevs log2V _flat) are −8.6 ± 1.1, −10.0 ±1.5, −10.7 ±2.1, and −9.7 ±1.3 in the R, J, H, and K bands respectively for galaxies having flat HI rotation curves. These values of the slopes are consistent with those obtained from studies of other groups and clusters. The scatter in the TF relations is in the range 0.5-1.1 mag in different bands. This scatter is considerably larger compared to those observed in other groups and clusters. It is suggested that the larger scatter in the TF relations for the Eridanus group is related to the loose structure of the group. If the TF relations are constructed using the baryonic mass (stellar +HI + Helium mass) instead of the stellar luminosity, nearly identical slopes are obtained in the R and in the near-infrared bands. The baryonic TF (baryonic massvs log2V _flat) slope is in the range 3.5–4.1.     

Pulsational instability of yellow hypergiants

Instability of population I (X = 0.7, Z = 0.02) massive stars against radial oscillations during the post-main-sequence gravitational contraction of the helium core is investigated. Initial stellar masses are in the range 65M _⊙ ≤ M _ZAMS ≤ 90M _⊙. In hydrodynamic computations of self-exciting stellar oscillations we assumed that energy transfer in the envelope of the pulsating star is due to radiative heat conduction and convection. The convective heat transfer was treated in the framework of the theory of time-dependent turbulent convection. During evolutionary expansion of outer layers after hydrogen exhaustion in the stellar core the star is shown to be unstable against radial oscillations while its effective temperature is T _eff > 6700 K for M _ZAMS = 65M _⊙ and T _eff > 7200 K for M _ZAMS = 90M _⊙. Pulsational instability is due to the κ-mechanism in helium ionization zones and at lower effective temperature oscillations decay because of significantly increasing convection. The upper limit of the period of radial pulsations on this stage of evolution does not exceed ≈200 day. Radial oscillations of the hypergiant resume during evolutionary contraction of outer layers when the effective temperature is T _eff > 7300 K for M _ZAMS = 65M _⊙ and T _eff > 7600 K for M _ZAMS = 90M _⊙. Initially radial oscillations are due to instability of the first overtone and transition to fundamental mode pulsations takes place at higher effective temperatures (T _eff > 7700 K for M _ZAMS = 65M _⊙ and T _eff > 8200 K for M _ZAMS = 90M _⊙). The upper limit of the period of radial oscillations of evolving blueward yellow hypergiants does not exceed ≈130 day. Thus, yellow hypergiants are stable against radial stellar pulsations during the major part of their evolutionary stage.     

The Restricted Three-Body Problem: Comments on the Rotational Effect

The synchronization between the orbital motion and axial rotation of the two component stars of a binary system is reviewed. Some previous published papers are mentioned and the general conclusion is outlined: If we shall use a rotating coordinate system synchronous with one of the two stellar axial rotations, it is not possible to obtain a Jacobi integral and the Roche geometry cannot be further analyzed. In addition, a theoretical approach is summarized in order to use the axial rotations of the two component stars, even if the constants of the stellar structure (k₂)₁, (k₂)₂must be taken into consideration. So it is found that if the stellar angular velocities are higher than the corresponding Keplerian angular velocity (ω_i≫ ω_k, i= ), the problem of the rotational effect could be of practical consideration. Finally, a theoretical relationship between the two constants (k₂)₁and (k₂)₂of the stellar structure is established.     

Candidate red clump giants in the PPMX catalogue

This paper presents the RCGP catalogue of more than 0.5 million candidate red clump stars with the limiting magnitude K _s = 9.5 ^m . These stars are selected from the PPMX catalogue as the most probable red clump members by analyzing the color-reduced proper motion diagrams built from the proper motions given in PPMX and J, K _s -photometry given in the 2MASS catalogue. Reddening of the selected stars is used to find extinction in the K _s -band and to consider it in the further analysis. The two-dimensional galactic rotation model generalized by Ogorodnikov is used to investigate the tangential velocity field of the selected red clump members, most of which are thin disk stars located within 1.5 kpc from the sun. The values of kinematic parameters and solar components are determined as a function of stellar heights above the galactic equatorial plane and their heliocentric distances.     

On the instability of weakly radially anisotropic star clusters

We discuss contradictions existing in the literature in the problem on the stability of collisionless spherical stellar systems, which are the simplest anisotropic generalization of the well-known polytropic models. On the one hand, calculations of the growth rates within the framework of a linear stability theory and N-body simulations suggest that these systems should become stable when the parameter s characterizing the degree of anisotropy of the stellar velocity distribution becomes lower than some critical value s _crit > 0. On the other hand, according to Palmer and Papaloizou, the growth rate should be nonzero up to the isotropic limit s = 0. Using our method of determining the eigenmodes of stellar systems, we show that even though the mode growth rates in weakly radially anisotropic systems of this type are nonzero, they are exponentially small, i.e., decrease as γ ∝ exp(−a/s) when s → 0. For slightly radially anisotropic systems with a finite lifetime, this actually implies stability.     

Using the reduced proper motions of Tycho-2 stars with <Emphasis Type="Italic">B-V</Emphasis> from 0.75 to 1.25: Monte Carlo simulations

The stellar composition of the Tycho-2 Catalogue in the range B-V = 0_· ^m 75–1_· ^m 25 has been reproduced through Monte Carlo simulations. For young and old stars of the red giant clump (RGC), the red giant branch, subgiants, red dwarfs, and thick-disk giants, we have specified the distributions in coordinates, velocities, B-V, and M _V as a function of B-V and calculated their reduced proper motions, photometric distances from the (B-V)-M _V calibration, and photoastrometric distances from the reduced proper motion-M _V calibration. Our simulations have shown the following: (1) a sample of thin-disk giants within 500 pc with an admixture of less than 10% of other stars can be produced; (2) a sample of dwarfs within 100 pc almost without any admixture of other stars can be produced; (3) the Local Spiral Arm affects the RGC composition of any magnitude-limited catalog in favor of giants younger than 2 Gyr; (4) the samples produced using reduced proper motions can be used for kinematic studies, provided that the biases of the quantities being determined are simulated and taken into account; (5) the photometric distances correlate with the photoastrometric ones because of the correlation between the proper motion and magnitude; (6) the photometric distances are closer to the true ones for the red giant branch and red dwarfs as the categories of stars with a clear (B-V)-M _V relation, while the photoastrometric distances are closer to the true ones for the RGC, subgiants, and thick-disk giants; (7) the calculated distances differ systematically from the true ones, but they can be used to analyze the three-dimensional distribution of stars. Our simulations confirm the validity of our previous selection of RGC stars from Tycho-2.     

Three-component dust models for interstellar extinction

Interstellar extinction curves obtained from the ‘extinction without standard’ method were used to constrain the dust characteristics in the mean ISM (R _V = 3.1), along the lines of sight through a high latitude diffuse molecular cloud towards HD 210121 (R _V = 2.1) and in a dense interstellar environment towards the cluster NGC 1977 (R _V = 6.42). We have used three-component dust models comprising silicate, graphite and very small carbonaceous grains (polycyclic aromatic hydrocarbons) following the grain size distributions introduced by Li & Draine in 2001. It is shown that oxygen, carbon and silicon abundances derived from our models are closer with the available elemental abundances for the dust grains in the ISM if F & G type stars atmospheric abundances are taken for the ISM than the solar. The importance of very small grains in modelling the variation of interstellar extinction curves has been investigated. Grain size distributions and elemental abundances locked up in dust are studied and compared at different interstellar environments using these three extinction curves. We present the albedo and the scattering asymmetry parameter evaluated from optical to extreme-UV wavelengths for the proposed dust models.