Context. Accreted material onto CTTSs is expected to form a hot quasi-periodic plasma structure that radiates in X-rays. Simulations of this phenomenon only partly match with observations. They all rely on a static model for the chromosphere model and on the assumption that radiation and matter are decoupled. Aims. We explore the effects on the structure and on the dynamics of the accretion flow of both a shock-heated chromosphere and of the coupling between radiation and hydrodynamics. Methods. We simulate accretion columns falling onto a stellar chromosphere using the 1D ALE code AstroLabE. This code solves the hydrodynamics equations along with the two first momenta equations for radiation transfer, with the help of a dedicated opacity table for the coupling between matter and radiation. We derive the total electron and ions densities from collisional-radiative NLTE ionisation equilibrium. Results. The chromospheric acoustic heating has an impact on the duration of the cycle and on the structure of the heated slab. In addition, the coupling between radiation and hydrodynamics leads to a heating of the accretion flow and the chromosphere, inducing a possible unburial of the whole column. These two last conclusions are in agreement with the computed monochromatic intensity. Both effects (acoustic heating and radiation coupling) have an influence on the amplitude and temporal variations of the net X-ray luminosity, which varies between 30 and 94% of the incoming mechanical energy flux, depending on the model considered.

We have performed a census of the UV-bright population in 78 globular clusters using wide-field UV telescopes. This population includes a variety of phases of post-horizontal branch (HB) evolution, including hot post-asymptotic giant branch (AGB) stars, and post-early AGB stars. There are indications that old stellar systems like globular clusters produce fewer post-(early) AGB stars than currently predicted by evolutionary models, but observations are still scarce. We obtained FORS2 spectroscopy of eleven of these UV-selected objects (covering a range of -2.3<[Fe/H]<-1.0), which we (re-)analysed together with previously observed data. We used model atmospheres of different metallicities, including super-solar ones. Where possible, we verified our atmospheric parameters using UV spectrophotometry and searched for metal lines in the optical spectra. We calculated evolutionary sequences for four metallicity regimes and used them together with information about the HB morphology of the globular clusters to estimate the expected numbers of post-AGB stars. Seven of the eleven new luminous UV-bright stars are post-AGB or post-early AGB stars, two are evolving away from the HB, one is a foreground white dwarf, and one is a white dwarf merger. So spectroscopy is clearly required to identify the evolutionary status of hot UV-bright stars. For hotter stars, metal-rich model spectra are required to reproduce their optical and UV spectra, which may affect the flux contribution of hot post-AGB stars to the UV spectra of evolved populations. Adding published information on other hot UV-bright stars in globular clusters, we find that the number of observed hot post-AGB stars generally agrees with the predicted values, although the numbers are still low.

Using a sample of 31 main-sequence OB stars located between galactocentric distances 8.4 - 15.6 kpc, we aim to probe the present-day radial abundance gradients of the Galactic disk. The analysis is based on high-resolution spectra obtained with the MIKE spectrograph on the Magellan Clay 6.5-m telescope on Las Campanas. We used a non-NLTE analysis in a self-consistent semi-automatic routine based on TLUSTY and SYNSPEC to determine atmospheric parameters and chemical abundances. Stellar parameters (effective temperature, surface gravity, projected rotational velocity, microturbulence, and macroturbulence) and silicon and oxygen abundances are presented for 28 stars located beyond 9 kpc from the Galactic centre plus three stars in the solar neighborhood. The stars of our sample are mostly on the main-sequence, with effective temperatures between 20800 - 31300 K, and surface gravities between 3.23 - 4.45 dex. The radial oxygen and silicon abundance gradients are negative and have slopes of -0.07 dex/kpc and -0.09 dex/kpc, respectively, in the region $8.4 \leq R_G \leq 15.6$\,kpc. The obtained gradients are compatible with the present-day oxygen and silicon abundances measured in the solar neighborhood and are consistent with radial metallicity gradients predicted by chemodynamical models of Galaxy Evolution for a subsample of young stars located close to the Galactic plane.

This paper presents a detailed operational manual for TLUSTY. It provides a guide for understanding the essential features and the basic modes of operation of the program. To help the user, it is divided into two parts. The first part describes the most important input parameters and available numerical options. The second part covers additional details and a comprehensive description of all physical and numerical options, and a description of all input parameters, many of which needed only in special cases.

This is the first of three papers that present a detailed guide for working with the codes {\sc tlusty} and {\sc synspec} to generate model stellar atmospheres or accretion disks, and to produce detailed synthetic spectra. In this paper, we present a very brief manual intended for casual users who intend to use these codes for simple, well defined tasks. This paper does not present any background theory, or a description of the adopted numerical approaches, but instead uses simple examples to explain how to employ these codes. In particular, it shows how to produce a simple model atmosphere from the scratch, or how to improve an existing model by considering more extended model atoms. This paper also presents a brief guide to the spectrum synthesis program {\sc synspec}.

This is the second part of a three-volume guide to TLUSTY and SYNSPEC. It presents a detailed reference manual for TLUSTY, which contains a detailed description of basic physical assumptions and equations used to model an atmosphere, together with an overview of the numerical methods to solve these equations.

This paper presents the sharpest near-IR images of the massive cluster R136 to date, based on the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO/VLT and operated in its IRDIS imaging mode. Stacking-up a few hundreds of short exposures in J and Ks spectral bands over a FoV of 10.9"x12.3" centered on the R136a1 stellar component, enabled us to carry a refined photometric analysis of the core of R136. We detected 1110 and 1059 sources in J and Ks images respectively with 818 common sources. Thanks to better angular resolution and dynamic range, we found that more than 62.6% (16.5%) of the stars, detected both in J and Ks data, have neighbours closer than 0.2" (0.1"). Among resolved and detected sources R136a1 and R136c have visual companions and R136a3 is resolved as two stars separated by 59mas. The new set of detected sources were used to re-assess the age and extinction of R136 based on 54 spectroscopically stars that have been recently studied with HST slit-spectroscopy. Over 90% of these 54 sources identified visual companions (closer than 0.2"). We found the most probable age and extinction for these sources within the photometric and spectroscopic error-bars. Additionally, using PARSEC evolutionary isochrones and tracks, we estimated the stellar mass range for each detected source (common in J and K data) and plotted the generalized histogram of mass (MF with error-bars). Using SPHERE data, we have gone one step further and partially resolved and studied the IMF covering mass range of (3-300) Msun at the age of 1 and 1.5 Myr. The density in the core of R136 is estimated and extrapolated in 3D and larger radii (up to 6pc). We show that the stars in the core are still unresolved due to crowding, and the results we obtained are upper limits. Higher angular resolution is mandatory to overcome these difficulties.

We present mid-infrared observations for a sample of 16 O-type stars. The data were acquired with the NASA Spitzer Space Telescope, using the IRS instrument at moderate resolution (R $\sim$ 600), covering the range $\sim 10-37$ microns. Our sample includes early, mid and late O supergiants and dwarfs. We explore for the first time their mid-IR spectral morphology in a quantitative way. We use NLTE expanding atmosphere models to help with line identifications, analyze profile contributions and line-formation regions. The O supergiants present a rich emission line spectra. The most intense features are from hydrogen - $6\alpha$, 7$\alpha$, and $8\alpha$ - which have non-negligible contributions of HeI or HeII lines, depending on the spectral type. The spectrum of early O supergiants is a composite of HI and HeII lines, HeI lines being absent. On the other hand, late O supergiants present features composed mainly by HI and HeI lines. All emission lines are formed throughout the stellar wind. We found that O dwarfs exhibit a featureless mid-IR spectrum. Two stars of our sample exhibit very similar mid-IR features, despite having a very different optical spectral classification. The analysis of O-type stars based on mid-IR spectra alone to infer spectral classes or to estimate physical parameters may thus be prone to substantial errors. Our results may therefore inform spectroscopic observations of massive stars located in heavily obscured regions and help establish an initial framework for observations of massive stars using the Mid-Infrared Instrument on the James Webb Space Telescope.

The globular cluster (GC) system of our Galaxy contains four planetary nebulae (PNe): K 648 (or Ps 1) in M15, IRAS 18333-2357 in M22, JaFu 1 in Pal 6, and JaFu 2 in NGC 6441. Because single-star evolution at the low stellar mass of present-epoch GCs was considered incapable of producing visible PNe, their origin presented a puzzle. We imaged the PN JaFu 1 with the Hubble Space Telescope (HST) to obtain photometry of its central star (CS) and high-resolution morphological information. We imaged IRAS 18333-2357 with better depth and resolution, and we analyzed its archival HST spectra to constrain its CS temperature and luminosity. All PNe in Galactic GCs now have high-quality HST data, allowing us to improve CS mass estimates. We find reasonably consistent masses between 0.53 and 0.58 Msun for all four objects, though estimates vary when adopting different stellar evolutionary calculations. The CS mass of IRAS 18333-2357, though, depends strongly on its temperature, which remains elusive due to reddening uncertainties. For all four objects, we consider their CS and nebular masses, their morphologies, and other incongruities to assess the likelihood that these objects formed from binary stars. Although generally limited by uncertainties (~0.02 Msun) in post-AGB tracks and core mass vs. luminosity relations, the high-mass CS in K 648 indicates a binary origin. The CS of JaFu 1 exhibits compact bright [O III] and Halpha emission, like EGB 6, suggesting a binary companion or disk. Evidence is weaker for a binary origin of JaFu 2.

We present an overview of the scientific potential of MATISSE, the Multi Aperture mid-Infrared SpectroScopic Experiment for the Very Large Telescope Interferometer. For this purpose we outline selected case studies from various areas, such as star and planet formation, active galactic nuclei, evolved stars, extrasolar planets, and solar system minor bodies and discuss strategies for the planning and analysis of future MATISSE observations. Moreover, the importance of MATISSE observations in the context of complementary high-angular resolution observations at near-infrared and submillimeter/millimeter wavelengths is highlighted.

MATISSE is the second-generation mid-infrared spectrograph and imager for the Very Large Telescope Interferometer (VLTI) at Paranal. This new interferometric instrument will allow significant advances by opening new avenues in various fundamental research fields: studying the planet-forming region of disks around young stellar objects, understanding the surface structures and mass loss phenomena affecting evolved stars, and probing the environments of black holes in active galactic nuclei. As a first breakthrough, MATISSE will enlarge the spectral domain of current optical interferometers by offering the L and M bands in addition to the N band. This will open a wide wavelength domain, ranging from 2.8 to 13 um, exploring angular scales as small as 3 mas (L band) / 10 mas (N band). As a second breakthrough, MATISSE will allow mid-infrared imaging - closure-phase aperture-synthesis imaging - with up to four Unit Telescopes (UT) or Auxiliary Telescopes (AT) of the VLTI. Moreover, MATISSE will offer a spectral resolution range from R ~ 30 to R ~ 5000. Here, we present one of the main science objectives, the study of protoplanetary disks, that has driven the instrument design and motivated several VLTI upgrades (GRA4MAT and NAOMI). We introduce the physical concept of MATISSE including a description of the signal on the detectors and an evaluation of the expected performances. We also discuss the current status of the MATISSE instrument, which is entering its testing phase, and the foreseen schedule for the next two years that will lead to the first light at Paranal.

MATISSE represents a great opportunity to image the environment around massive and evolved stars. This will allow one to put constraints on the circumstellar structure, on the mass ejection of dust and its reorganization , and on the dust-nature and formation processes. MATISSE measurements will often be pivotal for the understanding of large multiwavelength datasets on the same targets collected through many high-angular resolution facilities at ESO like sub-millimeter interferometry (ALMA), near-infrared adaptive optics (NACO, SPHERE), interferometry (PIONIER, GRAVITY), spectroscopy (CRIRES), and mid-infrared imaging (VISIR). Among main sequence and evolved stars, several cases of interest have been identified that we describe in this paper.

We present new near-infrared photometric measurements of the core of the young massive cluster NGC 3603 obtained with extreme adaptive optics. The data were obtained with the SPHERE instrument mounted on ESO Very Large Telescope, and cover three fields in the core of this cluster. We applied a correction for the effect of extinction to our data obtained in the J and K broadband filters and estimated the mass of detected sources inside the field of view of SPHERE/IRDIS, which is 13.5"x13.5". We derived the mass function (MF) slope for each spectral band and field. The MF slope in the core is unusual compared to previous results based on Hubble space telescope (HST) and very large telescope (VLT) observations. The average slope in the core is estimated as -1.06^{+0.26}_{-0.26} for the main sequence stars with 3.5 Msun < M < 120 Msun.Thanks to the SPHERE extreme adaptive optics, 814 low-mass stars were detected to estimate the MF slope for the pre-main sequence stars with 0.6 Msun< M < 3.5 Msun , Gamma = -0.54^{+0.11}_{-0.11} in the K-band images in two fields in the core of the cluster. For the first time, we derive the mass function of the very core of the NGC 3603 young cluster for masses in the range 0.6 - 120 Msun. Previous studies were either limited by crowding, lack of dynamic range, or a combination of both.

We present a spectroscopic analysis of HST/COS observations of three massive stars in the low metallicity dwarf galaxies IC 1613 and WLM. These stars, were previously observed with VLT/X-shooter by Tramper et al. (2011, 2014) who claimed that their mass-loss rates are higher than expected from theoretical predictions for the underlying metallicity. A comparison of the FUV spectra with those of stars of similar spectral types/luminosity classes in the Galaxy, and the Magellanic Clouds provides a direct, model-independent check of the mass-loss - metallicity relation. Then, a quantitative spectroscopic analysis is carried out using the NLTE stellar atmosphere code CMFGEN. We derive the photospheric and wind characteristics, benefiting from a much better sensitivity of the FUV lines to wind properties than Ha. Iron and CNO abundances are measured, providing an independent check of the stellar metallicity. The spectroscopic analysis indicates that Z/Zsun = 1/5, similar to a SMC-type environment, and higher than usually quoted for IC 1613 and WLM. The mass-loss rates are smaller than the empirical ones by Tramper et al. (2014), and those predicted by the widely used theoretical recipe by Vink et al. (2001). On the other hand, we show that the empirical, FUV-based, mass-loss rates are in good agreement with those derived from mass fluxes computed by Lucy (2012). We do not concur with Tramper et al. (2011, 2014) that there is a breakdown in the mass-loss - metallicity relation.

The aim of this study is to analyse and determine elemental abundances for a large sample of distant B stars in the outer Galactic disk in order to constrain the chemical distribution of the Galactic disk and models of chemical evolution of the Galaxy. Here, we present preliminary results on a few stars along with the adopted methodology based on securing simultaneous O and Si ionization equilibria with consistent NLTE model atmospheres.

We explore a method for metallicity determinations based on quantitative spectroscopy of type II-Plateau (II-P) supernovae (SNe). For consistency, we first evolve a set of 15Msun main sequence stars at 0.1, 0.4, 1, and 2 x the solar metallicity. At the onset of core collapse, we trigger a piston-driven explosion and model the resulting ejecta and radiation. Our theoretical models of such red-supergiant-star explosions at different metallicity show that synthetic spectra of SNe II-P possess optical signatures during the recombination phase that are sensitive to metallicity variations. This sensitivity can be quantified and the metallicity inferred from the strength of metal-line absorptions. Furthermore, these signatures are not limited to O, but also include Na, Ca, Sc, Ti, or Fe. When compared to a sample of SNe II-P from the Carnegie SN Project and previous SN followup programs, we find that most events lie at a metallicity between 0.4 and 2 x solar, with a marked scarcity of SN II-P events at SMC metallicity. This most likely reflects the paucity of low metallicity star forming regions in the local Universe. SNe II-P have high plateau luminosities that make them observable spectroscopically at large distances. Because they exhibit signatures of diverse metal species, in the future they may offer a means to constrain the evolution of the composition (e.g., the O/Fe ratio) in the Universe out to a redshift of one and beyond.

We present the results of a Hubble Space Telescope program to search for pulsating hot subdwarfs in the core of NGC 2808. These observations were motivated by the recent discovery of such stars in the outskirts of omega Cen. Both NGC 2808 and omega Cen are massive globular clusters exhibiting complex stellar populations and large numbers of extreme horizontal branch stars. Our far-UV photometric monitoring of over 100 hot evolved stars has revealed six pulsating subdwarfs with periods ranging from 85 to 149 s and UV amplitudes of 2.0 to 6.8%. In the UV color-magnitude diagram of NGC 2808, all six of these stars lie immediately below the canonical horizontal branch, a region populated by the subluminous "blue-hook" stars. For three of these six pulsators, we also have low-resolution far-UV spectroscopy that is sufficient to broadly constrain their atmospheric abundances and effective temperatures. Curiously, and in contrast to the omega Cen pulsators, the NGC 2808 pulsators do not exhibit the spectroscopic or photometric uniformity one might expect from a well-defined instability strip, although they all fall within a narrow band (0.2 mag) of far-UV luminosity.

(abridged) AIMS. We investigate the dynamics and stability of post-shock plasma streaming along nonuniform stellar magnetic fields at the impact region of accretion columns. We study how the magnetic field configuration and strength determine the structure, geometry, and location of the shock-heated plasma. METHODS. We model the impact of an accretion stream onto the chromosphere of a CTTS by 2D axisymmetric magnetohydrodynamic simulations. Our model takes into account the gravity, the radiative cooling, and the magnetic-field-oriented thermal conduction. RESULTS. The structure, stability, and location of the shocked plasma strongly depend on the configuration and strength of the magnetic field. For weak magnetic fields, a large component of B may develop perpendicular to the stream at the base of the accretion column, limiting the sinking of the shocked plasma into the chromosphere. An envelope of dense and cold chromospheric material may also develop around the shocked column. For strong magnetic fields, the field configuration determines the position of the shock and its stand-off height. If the field is strongly tapered close to the chromosphere, an oblique shock may form well above the stellar surface. In general, a nonuniform magnetic field makes the distribution of emission measure vs. temperature of the shocked plasma lower than in the case of uniform magnetic field. CONCLUSIONS. The initial strength and configuration of the magnetic field in the impact region of the stream are expected to influence the chromospheric absorption and, therefore, the observability of the shock-heated plasma in the X-ray band. The field strength and configuration influence also the energy balance of the shocked plasma, its emission measure at T > 1 MK being lower than expected for a uniform field. The above effects contribute in underestimating the mass accretion rates derived in the X-ray band.

Context. Theoretical arguments and numerical simulations of radiative shocks produced by the impact of the accreting gas onto young stars predict quasi-periodic oscillations in the emitted radiation. However, observational data do not show evidence of such periodicity. Aims. We investigate whether physically plausible perturbations in the accretion column or in the chromosphere could disrupt the shock structure influencing the observability of the oscillatory behavior. Methods. We performed local 2D magneto-hydrodynamical simulations of an accretion shock impacting a chromosphere, taking optically thin radiation losses and thermal conduction into account. We investigated the effects of several perturbation types, such as clumps in the accretion stream or chromospheric fluctuations, and also explored a wide range of plasma-\beta values. Results. In the case of a weak magnetic field, the post-shock region shows chaotic motion and mixing, smoothing out the perturbations and retaining a global periodic signature. On the other hand, a strong magnetic field confines the plasma in flux tubes, which leads to the formation of fibrils that oscillate independently. Realistic values for the amplitude, length, and time scales of the perturbation are capable of bringing the fibril oscillations out of phase, suppressing the periodicity of the emission. Conclusions. The strength of a locally uniform magnetic field in YSO accretion shocks determines the structure of the post-shock region, namely, whether it will be somewhat homogeneous or if it will split up to form a collection of fibrils. In the second case, the size and shape of the fibrils is found to depend strongly on the plasma-\beta value but not on the perturbation type. Therefore, the actual value of the protostellar magnetic field is expected to play a critical role in the time dependence of the observable emission.

We study the evolution, rotation, and surface abundances of O-type dwarfs in the Small Magellanic Cloud. We analyzed the UV and optical spectra of twenty-three objects and derived photospheric and wind properties. The observed binary fraction of the sample is ~ 26%, which is compatible with more systematic studies, if one considers that the actual binary fraction is potentially larger owing to low-luminosity companions and that the sample excluded obvious spectroscopic binaries. The location of the fastest rotators in the H-R diagram indicates that these could be several Myr old. The offset in the position of these fast rotators compared with the other stars confirms the predictions of evolutionary models that fast-rotating stars tend to evolve more vertically in the H-R diagram. Only one star of luminosity-class Vz, expected to best characterize extreme youth, is located on the ZAMS, the other two stars are more evolved. The distribution of nitrogen abundance of O and B stars suggests that the mechanisms responsible for the chemical enrichment of slowly rotating massive stars depends only weakly on the star's mass. We confirm that the group of slowly rotating N-rich stars is not reproduced by the evolutionary tracks. Our results call for stronger mixing in the models to explain the range of observed N abundances. All stars have an N/C ratio as a function of stellar luminosity that matches the predictions of the stellar evolution models well. More massive stars have a higher N/C ratio than the less massive stars. Faster rotators show on average a higher N/C ratio than slower rotators. The N/O versus N/C ratios agree qualitatively well with those of stellar evolution models. The only discrepant behavior is observed for the youngest two stars of the sample, which both show very strong signs of mixing, which is unexpected for their evolutionary status.

We present IRIS, a new generic three-dimensional (3D) spectral radiative transfer code that generates synthetic spectra, or images. It can be used as a diagnostic tool for comparison with astrophysical observations or laboratory astrophysics experiments. We have developed a 3D short-characteristic solver that works with a 3D nonuniform Cartesian grid. We have implemented a piecewise cubic, locally monotonic, interpolation technique that dramatically reduces the numerical diffusion effect. The code takes into account the velocity gradient effect resulting in gradual Doppler shifts of photon frequencies and subsequent alterations of spectral line profiles. It can also handle periodic boundary conditions. This first version of the code assumes Local Thermodynamic Equilibrium (LTE) and no scattering. The opacities and source functions are specified by the user. In the near future, the capabilities of IRIS will be extended to allow for non-LTE and scattering modeling. IRIS has been validated through a number of tests. We provide the results for the most relevant ones, in particular a searchlight beam test, a comparison with a 1D plane-parallel model, and a test of the velocity gradient effect. IRIS is a generic code to address a wide variety of astrophysical issues applied to different objects or structures, such as accretion shocks, jets in young stellar objects, stellar atmospheres, exoplanet atmospheres, accretion disks, rotating stellar winds, cosmological structures. It can also be applied to model laboratory astrophysics experiments, such as radiative shocks produced with high power lasers.

We study the formation of photospheric emission lines in O stars and show that the rectangular profiles, sometimes double peaked, that are observed for some stars are a direct consequence of rotation, and it is unnecessary to invoke an enhanced density structure in the equatorial regions. Emission lines, such as N IV 4058 and the N III 4634-4640-4642 multiplet, exhibit non-standard "limb darkening" laws. The lines can be in absorption for rays striking the center of the star and in emission for rays near the limb. Weak features in the flux spectrum do not necessarily indicate an intrinsically weak feature -- instead the feature can be weak because of cancellation between absorption in "core" rays and emission from rays near the limb. Rotation also modifies line profiles of wind diagnostics such as He II 4686 and Halpha and should not be neglected when inferring the actual stratification, level and nature of wind structures.

We aim to constrain the properties and evolutionary status of early and mid-spectral type supergiants (from O4 to O7.5). These posses the highest mass-loss rates among the O stars, and exhibit conspicuous wind profiles. Using the non-LTE wind code CMFGEN, we simultaneously analyzed the FUV-UV and optical spectral range to determine the photospheric properties and wind parameters. We derived effective temperatures, luminosities, surface gravities, surface abundances, mass-loss rates, wind terminal velocities, and clumping filling factors. The supergiants define a very clear evolutionary sequence, in terms of ages and masses, from younger and more massive stars to older stars with lower initial masses. O4 supergiants cluster around the 3 Myr isochrone and are more massive than 60 Msun, while the O5 to O7.5 stars have masses in the range 50 - 40 Msun and are 4 +/- 0.3 Myr old. The surface chemical composition is typical of evolved O supergiants (nitrogen-rich, carbon- and oxygen-poor). While the observed ranges of carbon and nitrogen mass-fractions are compatible with those expected from evolutionary models for the measured stellar masses, the N/C ratios as a function of age are inconsistent with the theoretical predictions for the four earliest (O4 spectral type) stars of the sample. We question the efficiency of rotational mixing as a function of age for these stars and suggest that another mechanism may be needed to explain the observed abundance patterns. Mass-loss rates derived with clumped-models range within a factor of three of the theoretical mass-loss rates. The corresponding volume-filling factors associated with small-scale clumping are 0.05 +/- 0.02. Clumping is found to start close to the photosphere for all but three stars, two of which are fast rotators.

Blue hook stars are a class of subluminous extreme horizontal branch stars that were discovered in UV images of the massive globular clusters omega Cen and NGC 2808. These stars occupy a region of the HR diagram that is unexplained by canonical stellar evolution theory. Using new theoretical evolutionary and atmospheric models, we have shown that the blue hook stars are very likely the progeny of stars that undergo extensive internal mixing during a late helium-core flash on the white dwarf cooling curve. This "flash mixing" produces hotter-than-normal EHB stars with atmospheres significantly enhanced in helium and carbon. The larger bolometric correction, combined with the decrease in hydrogen opacity, makes these stars appear subluminous in the optical and UV. Flash mixing is more likely to occur in stars born with a high helium abundance, due to their lower mass at the main sequence turnoff. For this reason, the phenomenon is more common in those massive globular clusters that show evidence for secondary populations enhanced in helium. However, a high helium abundance does not, by itself, explain the presence of blue hook stars in massive globular clusters. Here, we present new observational evidence for flash mixing, using recent HST observations. These include UV color-magnitude diagrams of six massive globular clusters and far-UV spectroscopy of hot subdwarfs in one of these clusters (NGC 2808).

[Abridged] We present new HST FUV spectroscopy of 24 hot evolved stars in NGC2808, a massive globular cluster with a large population of "blue-hook" (BHk) stars. The BHk stars are found in UV color-magnitude diagrams of the most massive globular clusters, where they fall at luminosities immediately below the hot end of the horizontal branch (HB), in a region of the HR diagram unexplained by canonical stellar evolution theory. Using new evolutionary and atmospheric models, we have shown that these subluminous HB stars are very likely the progeny of stars that undergo extensive internal mixing during a late He-core flash on the white dwarf cooling curve. This flash mixing leads to hotter temperatures and an enormous enhancement of the surface He and C; these hotter temperatures, together with the decrease in H opacity shortward of the Lyman limit, make the BHk stars brighter in the EUV while appearing subluminous in the UV and optical. Our FUV spectroscopy demonstrates that, relative to normal HB stars at the same color, the BHk stars of NGC2808 are hotter and greatly enhanced in He and C, thus providing unambiguous evidence of flash mixing. Although the C abundance in the BHk stars is orders of magnitude larger than that in the normal HB stars, the C abundance in both the BHk and normal HB stars appears to be affected by gravitational settling. The variations seen in Si and the Fe-peak elements also indicate that atmospheric diffusion is at play in our sample, with all of our hot subdwarfs at 25,000 to 50,000 K exhibiting large enhancements of the Fe-peak elements. The hottest subdwarfs in our BHk sample may be pulsators, given that they fall in the temperature range of newly-discovered pulsating subdwarfs in omega Cen. In addition to the normal hot HB and BHk stars, we obtain spectra of 5 blue HB stars, a post-HB star, and 3 unclassified stars with unusually blue UV colors.

We present non-LTE oxygen abundances for a sample of B stars in the Orion association. The abundance calculations included non-LTE line formation and used fully blanketed non-LTE model atmospheres. The stellar parameters were the same as adopted in the previous study by Cunha & Lambert (1994). We find that the young Orion stars in this sample of 10 stars are described by a single oxygen abundance with an average value of A(O)=8.78 and a small dispersion of +/- 0.05 dex, which is of the order of the uncertainties in the analysis. This average oxygen abundance compares well with the average oxygen abundance obtained previously in Cunha & Lambert (1994): A(O) = 8.72 +/- 0.13 although this earlier study, based upon non-blanketed model atmospheres in LTE, displayed larger scatter. Small scatter of chemical abundances in Orion B stars had also been found in our previous studies for neon and argon; all based on the same effective temperature scale. The derived oxygen abundance distribution for the Orion association compares well with other results for the oxygen abundance in the solar neighborhood.

UV observations of some massive globular clusters have revealed a significant population of stars hotter and fainter than the hot end of the horizontal branch (HB), the so-called blue hook stars. This feature might be explained either by the late hot flasher scenario where stars experience the helium flash while on the white dwarf cooling curve or by the progeny of the helium-enriched sub-population postulated to exist in some clusters. Previous spectroscopic analyses of blue hook stars in omega Cen and NGC 2808 support the late hot flasher scenario, but the stars contain much less helium than expected and the predicted C and N enrichment cannot be verified. We compare the observed effective temperatures, surface gravities, helium abundances, and carbon line strengths (where detectable) of our targets stars with the predictions of the two aforementioned scenarios. Moderately high resolution spectra of hot HB stars in the globular cluster omega Cen were analysed for radial velocity variations, atmospheric parameters, and abundances using LTE and non-LTE model atmospheres. We find no evidence of close binaries among our target stars. All stars below 30,000K are helium-poor and very similar to HB stars observed in that temperature range in other globular clusters. In the temperature range 30,000K to 50,000K, we find that 28% of our stars are helium-poor (log(He/H) < -1.6), while 72% have roughly solar or super-solar helium abundance (log(He/H) >= -1.5). We also find that carbon enrichment is strongly correlated with helium enrichment, with a maximum carbon enrichment of 3% by mass. A strong carbon enrichment in tandem with helium enrichment is predicted by the late hot flasher scenario, but not by the helium-enrichment scenario. We conclude that the helium-rich HB stars in omega Cen cannot be explained solely by the helium-enrichment scenario invoked to explain the blue main sequence.

We present new HST ultraviolet color-magnitude diagrams of 5 massive Galactic globular clusters: NGC 2419, NGC 6273, NGC 6715, NGC 6388, and NGC 6441. These observations were obtained to investigate the "blue hook" phenomenon previously observed in UV images of the globular clusters omega Cen and NGC 2808. Blue hook stars are a class of hot (approximately 35,000 K) subluminous horizontal branch stars that occupy a region of the HR diagram that is unexplained by canonical stellar evolution theory. By coupling new stellar evolution models to appropriate non-LTE synthetic spectra, we investigate various theoretical explanations for these stars. Specifically, we compare our photometry to canonical models at standard cluster abundances, canonical models with enhanced helium (consistent with cluster self-enrichment at early times), and flash-mixed models formed via a late helium-core flash on the white dwarf cooling curve. We find that flash-mixed models are required to explain the faint luminosity of the blue hook stars, although neither the canonical models nor the flash-mixed models can explain the range of color observed in such stars, especially those in the most metal-rich clusters. Aside from the variation in the color range, no clear trends emerge in the morphology of the blue hook population with respect to metallicity.

We present a comprehensive multiwavelength analysis of the hot DB white dwarf PG 0112+104. Our analysis relies on newly-acquired FUSE observations, on medium-resolution FOS and GHRS data, on archival high-resolution GHRS observations, on optical spectrophotometry both in the blue and around Halpha, as well as on time-resolved photometry. From the optical data, we derive a self-consistent effective temperature of 31,300+-500 K, a surface gravity of log g = 7.8 +- 0.1 (M=0.52 Msun), and a hydrogen abundance of log N(H)/N(He) < -4.0. The FUSE spectra reveal the presence of CII and CIII lines that complement the previous detection of CII transitions with the GHRS. The improved carbon abundance in this hot object is log N(C)/N(He) = -6.15 +- 0.23. No photospheric features associated with other heavy elements are detected. We reconsider the role of PG 0112+104 in the definition of the blue edge of the V777 Her instability strip in light of our high-speed photometry, and contrast our results with those of previous observations carried out at the McDonald Observatory.

We report on the investigation of strong radiative shocks generated with the high energy, sub-nanosecond iodine laser at PALS. These shock waves are characterized by a developed radiative precursor and their dynamics is analyzed over long time scales (~50 ns), approaching a quasi-stationary limit. We present the first preliminary results on the rear side XUV spectroscopy. These studies are relevant to the understanding of the spectroscopic signatures of accretion shocks in Classical T Tauri Stars.

Advances in laser and Z-pinch technology, coupled with the development of plasma diagnostics and the availability of high-performance computers, have recently stimulated the growth of high-energy density laboratory astrophysics. In particular a number of experiments have been designed to study radiative shocks and jets with the aim of shedding new light on physical processes linked to the ejection and accretion of mass by newly born stars. Although general scaling laws are a powerful tools to link laboratory experiments with astrophysical plasmas, the phenomena modelled are often too complicated for simple scaling to remain relevant. Nevertheless, the experiments can still give important insights into the physics of astrophysical systems and can be used to provide the basic experimental validation of numerical simulations in regimes of interest to astrophysics. We will illustrate the possible links between laboratory experiments, numerical simulations and astrophysics in the context of stellar jets. First we will discuss the propagation of stellar jets in a cross-moving interstellar medium and the scaling to Z-pinch produced jets. Our second example focuses on slab-jets produced at the PALS (Prague Asterix Laser System) laser installation and their practical applications to astrophysics. Finally, we illustrate the limitations of scaling for radiative shocks, which are found at the head of the most rapid stellar jets.

Aims. Soft X-ray excesses have been detected in several Be/X-ray binaries and interpreted as the signature of hard X-ray reprocessing in the inner accretion disk. The system XMMU J054134.7-682550, located in the LMC, featured a giant Type II outburst in August 2007. The geometry of this system can be understood by studying the response of the soft excess emission to the hard X-ray pulses. Methods. We have analyzed series of simultaneous observations obtained with XMM-Newton/EPIC-MOS and RXTE/PCA in order to derive spectral and temporal characteristics of the system, before, during and after the giant outburst. Spectral fits were performed and a timing analysis has been carried out. Spectral variability, spin period evolution and energy dependent pulse shapes are analysed. Results. The outburst (L_X = 3* 10^38 erg/s \approx L_EDD) spectrum could be modeled successfully using a cutoff powerlaw, a cold disk emission, a hot blackbody, and a cyclotron absorption line. The magnetic field and magnetospheric radius could be constrained. The thickness of the inner accretion disk is broadened to a width of 75 km. The hot blackbody component features sinusoidal modulations indicating that the bulk of the hard X-ray emission is emitted preferentially along the magnetic equator. The spin period of the pulsar decreased very significantly during the outburst. This is consistent with a variety of neutron star equations of state and indicates a very high accretion rate.

White dwarf spectra have been widely used as a calibration source for X-ray and EUV instruments. The in-flight effective area calibration of the reflection grating spectrometers (RGS) of XMM-Newton depend upon the availability of reliable calibration sources. We investigate how well these white dwarf spectra can be used as standard candles at the lowest X-ray energies in order to gauge the absolute effective area scale of X-ray instruments. We calculate a grid of model atmospheres for Sirius B and HZ 43A, and adjust the parameters using several constraints until the ratio of the spectra of both stars agrees with the ratio as observed by the low energy transmission grating spectrometer (LETGS) of Chandra. This ratio is independent of any errors in the effective area of the LETGS. We find that we can constrain the absolute X-ray spectrum of both stars with better than 5 % accuracy. The best-fit model for both stars is close to a pure hydrogen atmosphere, and we put tight limits to the amount of helium or the thickness of a hydrogen layer in both stars. Our upper limit to the helium abundance in Sirius B is 4 times below the previous detection based on EUVE data. We also find that our results are sensitive to the adopted cut-off in the Lyman pseudo-continuum opacity in Sirius B. We get best agreement with a long wavelength cut-off. White dwarf model atmospheres can be used to derive the effective area of X-ray spectrometers in the lowest energy band. An accuracy of 3-4 % in the absolute effective area can be achieved.

We have investigated the stellar and wind properties of a sample of five late-type O dwarfs in order to address the weak wind problem. A grid of TLUSTY models was used to obtain the stellar parameters, and the wind parameters were determined by using the CMFGEN code. We found that the spectra have mainly a photospheric origin. A weak wind signature is seen in CIV 1549, from where mass-loss rates consistent with previous CMFGEN results regarding O8-9V stars were obtained. A discrepancy of roughly 2 orders of magnitude is found between these mass-loss rates and the values predicted by theory (Mdot(Vink)), confirming a breakdown or a steepening of the modified wind momentum-luminosity relation at log L/Lsun < 5.2. We have estimated the carbon abundance for the stars of our sample and concluded that its uncertainty cannot cause the weak wind problem. Upper limits on Mdot were established for all objects using lines of different ions, namely, PV 1118,28, CIII 1176, NV 1239,43, Si IV 1394,03, and NIV 1718. All the values obtained are also in disagreement with theoretical predictions, bringing support to the reality of weak winds. Together with CIV 1549, the use of NV 1239,43 results in the lowest mass-loss rates: the upper limits indicate that Mdot must be less than about -1.0 dex Mdot(Vink). Regarding the other transitions, the upper limits still point to low rates: Mdot must be less than about $(-0.5 \pm 0.2)$ dex Mdot(Vink). We have studied the behavior of the Halpha line with different mass-loss rates. We have also explored ways to fit the observed spectra with Mdot(Vink). By using large amounts of X-rays, we verified that few wind emissions take place, as in weak winds. However, unrealistic X-rays luminosities had to be used (log Lx/Lbol > -3.5) (abridged).

The Gaia satellite will be launched at the end of 2011. It will observe at least 1 billion stars, and among them several million emission line stars and hot stars. Gaia will provide parallaxes for each star and spectra for stars till V magnitude equal to 17. After a general description of Gaia, we present the codes and methods, which are currently developed by our team. They will provide automatically the astrophysical parameters and spectral classification for the hot and emission line stars in the Milky Way and other close Local Group galaxies such as the Magellanic Clouds.

We have studied superionization and X-ray line formation in the spectra of Zeta Pup using our new stellar atmosphere code (XCMFGEN) that can be used to simultaneously analyze optical, UV, and X-ray observations. Here, we present results on the formation of the O VI ll1032, 1038 doublet. Our simulations, supported by simple theoretical calculations, show that clumped wind models that assume void in the interclump space cannot reproduce the observed O VI profiles. However, enough O VI can be produced if the voids are filled by a low density gas. The recombination of O VI is very efficient in the dense material but in the tenuous interclump region an observable amount of O VI can be maintained. We also find that different UV resonance lines are sensitive to different density regimes in Zeta Pup : C IV is almost exclusively formed within the densest regions, while the majority of O VI resides between clumps. N V is an intermediate case, with contributions from both the tenuous gas and clumps.

We report here the detection of a weak magnetic field of 50 - 100 G on the O9.7 supergiant zeta Ori A, using spectropolarimetric observations obtained with NARVAL at the 2m Telescope Bernard Lyot atop Pic du Midi (France). zeta Ori A is the third O star known to host a magnetic field (along with theta^1 Ori C and HD 191612), and the first detection on a 'normal' rapidly-rotating O star. The magnetic field of zeta Ori A is the weakest magnetic field ever detected on a massive star. The measured field is lower than the thermal equipartition limit (about 100 G). By fitting NLTE model atmospheres to our spectra, we determined that zeta Ori A is a 40 Msun star with a radius of 25 Rsun and an age of about 5 - 6 Myr, showing no surface nitrogen enhancement and losing mass at a rate of about 2x10^(-6) Msol/yr. The magnetic topology of zeta Ori A is apparently more complex than a dipole and involves two main magnetic polarities located on both sides of the same hemisphere; our data also suggest that zeta Ori A rotates in about 7.0 d and is about 40 degrees away from pole-on to an Earth-based observer. Despite its weakness, the detected magnetic field significantly affects the wind structure; the corresponding Alfven radius is however very close to the surface, thus generating a different rotational modulation in wind lines than that reported on the two other known magnetic O stars. The rapid rotation of zeta Ori A with respect to theta^1 Ori C appears as a surprise, both stars having similar unsigned magnetic fluxes (once rescaled to the same radius); it may suggest that the sub-equipartition field detected on zeta Ori A is not a fossil remnant (as opposed to that of theta^1 Ori C and HD 191612), but the result of an exotic dynamo action produced through MHD instabilities.

Argon abundances have been derived for a sample of B main-sequence stars in the Orion association. The abundance calculations are based on NLTE metal line-blanketed model atmospheres calculated with the NLTE code TLUSTY and an updated and complete argon model atom. We derive an average argon abundance for this young population of A(Ar) = 6.66 +- 0.06. While our result is in excellent agreement with a recent analysis of the Orion nebula, it is significantly higher than the currently recommended solar value which is based on abundance measurements in the solar corona. Moreover, the derived argon abundances in the Orion B stars agree very well with a measurement from a solar impulsive flare during which unmodified solar photospheric material was brought to flare conditions. We therefore argue that the argon abundances obtained independently for both the Orion B stars and the Orion nebula are representative of the disk abundance value in the solar neighborhood. The lower coronal abundance may reflect a depletion related to the FIP effect. We propose a new reference value for the abundance of argon in the solar neighborhood, A(Ar) = 6.63 +- 0.10, corresponding to Ar/O = 0.009.

We present FUSE observations of the hot white dwarfs in the post-common envelope binaries Feige 24, EUVE J0720-317, BPM 6502, and EUVE J2013+400. The spectra show numerous photospheric absorption lines which trace the white dwarf orbital motion. We report the detection of C III, O VI, P V, and Si IV in the spectra of Feige 24, EUVE J0720-317 and EUVE J2013+400, and the detection of C III, N II, Si III, Si IV, and Fe III in the spectra of BPM 6502. Abundance measurements support the possibility that white dwarfs in post-common envelope binaries accrete material from the secondary star wind. The FUSE observations of BPM 6502 and EUVE J2013+400 cover a complete binary orbit. We used the FUSE spectra to measure the radial velocities traced by the white dwarf in the four binaries, where the zero-point velocity were fixed using the ISM velocities in the line of sight of the stellar systems. For BPM 6502 we determined a white dwarf velocity semi-amplitude of K_WD = 18.6+/-0.5km/s, and with the velocity semi-amplitude of the red dwarf companion (K_RD = 75.2+/-3.1 km/s), we estimate the mass ratio to be q = 0.25+/-0.01. Adopting a spectroscopic mass determination for the white dwarf, we infer a low secondary mass of M_RD = 0.14+/-0.01 M_solar. For EUVE J2013+400 we determine a white dwarf velocity semi-amplitude of K_WD = 36.7+/-0.7 km/s. The FUSE observations of EUVE J0720-317 cover approximately 30% of the binary period and combined with the HST GHRS measurements (Vennes et al. 1999, ApJ 523, 386), we update the binary properties. FUSE observations of Feige 24 cover approximately 60% of the orbit and we combine this data set with HST STIS (Vennes et al. 2000, ApJ, 544, 423) data to update the binary properties.

UV observations of some massive globular clusters uncovered a significant population of very hot stars below the hot end of the horizontal branch (HB), the so-called blue hook stars. This feature might be explained either as results of the late hot flasher scenario where stars experience the helium flash while on the white dwarf cooling curve or by the progeny of the helium-enriched sub-population recently postulated to exist in some clusters. Moderately high resolution spectra of stars at the hot end of the blue HB in omega Cen were analysed for atmospheric parameters and abundances using LTE and Non-LTE model atmospheres. In the temperature range 30,000K to 50,000K we find that 35% of our stars are helium-poor (log(n_He/n_H) < -2), 51% have solar helium abundance within a factor of 3 (-1.5 <= log(n_He/n_H) <= -0.5) and 14% are helium-rich (log(n_He/n_H)> -0.4). We also find carbon enrichment in step with helium enrichment, with a maximum carbon enrichment of 3% by mass. At least 14% of the hottest HB stars in omega Cen show helium abundances well above the highest predictions from the helium enrichment scenario (Y = 0.42 corresponding to log(n_He/n_H) ~ -0.74). In addition, the most helium-rich stars show strong carbon enrichment as predicted by the late hot flasher scenario. We conclude that the helium-rich HB stars in omega Cen cannot be explained solely by the helium-enrichment scenario invoked to explain the blue main sequence. (Abridged)

We report on initial results of our Spitzer Cycle 2 program to observe the young massive star R Mon and its associated HH 39 Herbig-Haro object in the mid-infrared. Our program used all instruments on-board Spitzer to obtain deep images with IRAC of the HH 39 complex and of R Mon and its surroundings, a deep image of HH 39 at 24 and 70 $\mu$m with MIPS, and mid-infrared spectra with the SH, LH, and LL modules of IRS. The aim of this program is to study the physical links in a young massive star between accretion disk, outflows and jets, and sh ocks in the associated HH object. Our preliminary analysis reveals that several knots of HH 39 are clearly detected in most IRAC bands. In IRAC4 (8 $\mu$m), diffuse emission, probably from PAHs, appears as foreground emission covering the HH 39 emission. The HH 39 knots are detected at 24 microns, despite the fact that dust continuum emission covers the knots and shows the same structure as observed with IRAC4. The IRS spectra of HH 39 show weak evidence of [Ne II] 12.8 $\mu$m and 0--0 S(1) H$_2$ 17.0 $\mu$m lines. A more detailed analysis is, however, required due to the faintness of the Herbig-Haro knots. Finally, we obtained the SH and MIPS SED spectra of R Mon. A PAH emission feature at 11.3 $\mu$m is detected on top of the strong continuum; although no strong emission or absorption lines are observed, we will seek to detect faint lines. The combined IRAC, IRS, and MIPS data of the R Mon/HH 39 system will help us to understand circumstellar disk processing, and the connection between jets, outflows, and HH objects.

Radiative shocks (also called supercritical shocks) are high Mach number shock waves that photoionize the medium ahead of the shock front and give rise to a radiative precursor. They are generated in the laboratory using high-energy or high-power lasers and are frequently present in a wide range of astronomical objects. Their modelisation in one dimension has been the subject of numerous studies, but generalization to three dimensions is not straightforward. We calculate analyticaly the absorption of radiation in a grey uniform cylinder and show how it decreases with chi R, the product of the opacity chi and of the cylinder radius R. Simple formulas, whose validity range increases when chi R diminishes, are derived for the radiation field on the axis of symmetry. Numerical calculations in three dimensions of the radiative energy density, flux and pressure created by a stationary shock wave show how the radiation decreases whith R. Finally, the bidimensional structures of both the precursor and the radiation field are calculated with time-dependent radiation hydrodynamics numerical simulations and the influence of two-dimensional effects on the electron density, the temperature, the shock velocity and the shock geometry are exhibited. These simulations show how the radiative precursor shortens, cools and slows down when R is decreased.

We have constructed a comprehensive grid of 1540 metal line-blanketed, NLTE, plane-parallel, hydrostatic model atmospheres for the basic parameters appropriate to early B-type stars. The BSTAR2006 grid considers 16 values of effective temperatures, 15,000 K <= Teff <= 30,000 K with 1,000 K steps, 13 surface gravities, 1.75 <= log g <= 4.75 with 0.25 dex steps, 6 chemical compositions, and a microturbulent velocity of 2 km/s. The lower limit of log g for a given effective temperature is set by an approximate location of the Eddington limit. The selected chemical compositions range from twice to one tenth of the solar metallicity and metal-free. Additional model atmospheres for B supergiants (log g <= 3.0) have been calculated with a higher microturbulent velocity (10 km/s) and a surface composition that is enriched in helium and nitrogen, and depleted in carbon. This new grid complements our earlier OSTAR2002 grid of O-type stars (Lanz & Hubeny, 2003, ApJS, 146, 417). The paper contains a description of the BSTAR2006 grid and some illustrative examples and comparisons. NLTE ionization fractions, bolometric corrections, radiative accelerations, and effective gravities are obtained over the parameter range covered by the grid. By extrapolating radiative accelerations, we have determined an improved estimate of the Eddington limit in absence of rotation between 55,000 and 15,000 K. The complete BSTAR2006 grid is available at the TLUSTY website (http://nova.astro.umd.edu).

We present an analysis of high-resolution FLAMES spectra of approximately 50 early B-type stars in three young clusters at different metallicities, NGC6611 in the Galaxy, N11 in the Large Magellanic Cloud (LMC) and NGC346 in the Small Magellanic Cloud (SMC). Using the TLUSTY non-LTE model atmospheres code, atmospheric parameters and photospheric abundances (C, N, O, Mg and Si) of each star have been determined. These results represent a significant improvement on the number of Magellanic Cloud B-type stars with detailed and homogeneous estimates of their atmospheric parameters and chemical compositions. The relationships between effective temperature and spectral type are discussed for all three metallicity regimes, with the effective temperature for a given spectral type increasing as one moves to a lower metallicity regime. Additionally the difficulties in estimating the microturbulent velocity and the anomalous values obtained, particularly in the lowest metallicity regime, are discussed. Our chemical composition estimates are compared with previous studies, both stellar and interstellar with, in general, encouraging agreement being found. Abundances in the Magellanic Clouds relative to the Galaxy are discussed and we also present our best estimates of the base-line chemical composition of the LMC and SMC as derived from B-type stars. Additionally we discuss the use of nitrogen as a probe of the evolutionary history of stars, investigating the roles of rotational mixing, mass-loss, blue loops and binarity on the observed nitrogen abundances and making comparisons with stellar evolutionary models where possible.

We present a study of the iron abundance pattern in hot hydrogen-rich (DA) white dwarfs. The study is based on new and archival far ultraviolet spectroscopy of a sample of white dwarfs in the temperature range 30,000 K < T_eff < 64,000 K. The spectra obtained with the Far Ultraviolet Spectroscopic Explorer along with spectra obtained with the Hubble Space Telescope Imaging Spectrograph and the International Ultraviolet Explorer sample FeIII to FeVI absorption lines enabling a detailed iron abundance analysis over a wider range of effective temperatures than previously afforded. The measurements reveal abundance variations in excess of two orders of magnitude between the highest and the lowest temperatures probed, but also show considerable variations (over one order of magnitude) between objects with similar temperatures and surface gravities. Such variations in cooler objects may be imputed to accretion from unseen companions or so-called circumstellar debris although the effect of residual mass-loss and selective radiation pressure in the hottest objects in the sample remain dominant.

We report on non-LTE Ne abundances for a sample of B-type stellar members of the Orion Association. The abundances were derived by means of non-LTE fully metal-blanketed model atmospheres and extensive model atoms with updated atomic data. We find that these young stars have a very homogeneous abundance of A(Ne) = 8.27 +/- 0.05. This abundance is higher by ~0.4 dex than currently adopted solar value, A(Ne)=7.84, which is derived from lines produced in the corona and active regions. The general agreement between the abundances of C, N, and O derived for B stars with the solar abundances of these elements derived from 3-D hydrodynamical models atmospheres strongly suggests that the abundance patterns of the light elements in the Sun and B stars are broadly similar. If this hypothesis is true, then the Ne abundance derived here is the same within the uncertainties as the value required to reconcile solar models with helioseismological observations.

Recent calculations and analyses of O star spectra have revealed discrepancies between theory and observations, and between different theoretical calculations, for the strength of optical HeI singlet transitions.We investigate the source of these discrepancies. Using a non-LTE radiative transfer code we have undertaken detailed test calculations for a range of O star properties. Our principal test model has parameters similar to those of the O9V star, 10 Lac. We show that the discrepancies arise from uncertainties in the radiation field in the HeI resonance transition near 584Angs. The radiation field at 584Angs. is influenced by model assumptions, such as the treatment of line-blanketing and the adopted turbulent velocity, and by the FeIV atomic data. It isshown that two FeIV transitions near 584Angs can have a substantial influence on the strength of the HeI singlet transitions. Because of the difficulty of modeling the HeI singlet lines, particularly in stars with solar metalicity, the HeI triplet lines should be preferred in spectral analyses. These lines are much less sensitive to model assumptions.

We present a spectroscopic analysis of HST/STIS and FOS low- and intermediate-resolution spectroscopy of 55 stars (turn-off stars, horizontal branch stars and blue stragglers) in four globular clusters (47 Tucanae, M3, NGC6752, and NGC6397). Stars were analyzed with non-Local Thermodynamic Equilibrium model atmospheres, and values for their effective temperatures and gravities and some rotation rates were obtained. Using photometric fluxes, we also obtained radii, luminosities and spectroscopic masses.

We have analyzed the far-UV spectrum of two Galactic O4 stars, the O4If+ supergiant HD190429A and the O4V((f)) dwarf HD96715, using archival FUSE and IUE data. We have conducted a quantitative analysis based on the two NLTE model atmosphere and wind codes, TLUSTY and CMFGEN. We have derived the stellar and wind parameters and the surface composition of the two stars. The surface of HD190429A has a composition typical of an evolved O supergiant (N-rich, C and O-poor), while HD96715 exhibits surface N enhancement similar to the enrichment found in SMC O dwarfs and attributed to rotationally-induced mixing. We find that homogeneous wind models could not match the observed profile of O V1371 and require very low phosphorus abundance to fit the P V1118-1128 resonance lines. However, we are able to match the O V and P V lines using clumped wind models. We find that N IV1718 is also sensitive to wind clumping. For both stars, we have calculated clumped wind models that match well all these lines from different species and that remain consistent with Halpha data. These fits therefore provide a coherent and thus much stronger evidence of wind clumping in O stars than earlier claims. We find that the wind of these two stars is highly clumped, as expressed by very small volume filling factors, namely f=0.04 for HD190429A and f=0.02 for HD96715. In agreement with our analysis of SMC stars, clumping starts deep in the wind, just above the sonic point. The most crucial consequence of our analysis is that the mass loss rates of O stars need to be revised downward significantly, by a factor of 3 and more. Accounting for wind clumping is essential when determining the wind properties of O stars. Our study therefore calls for a fundamental revision in our understanding of mass loss and of O-type star winds. (abridged)

An analysis of high-resolution VLT/UVES spectra of two B-type main sequence stars, NGC 346-11 and AV 304, in the Small Magellanic Cloud (SMC), has been undertaken, using the non-LTE TLUSTY model atmospheres to derive the stellar parameters and chemical compositions of each star. The chemical compositions of the two stars are in reasonable agreement. Moreover, our stellar analysis agrees well with earlier analyses of H II regions. The results derived here should be representative of the current base-line chemical composition of the SMC interstellar medium as derived from B-type stars.

High-resolution UCLES/AAT spectra are presented for nine B-type supergiants in the SMC, chosen on the basis that they may show varying amounts of nuclear-synthetically processed material mixed to their surface. These spectra have been analysed using a new grid of approximately 12 000 non-LTE line blanketed {\sc tlusty} model atmospheres to estimate atmospheric parameters and chemical composition. The abundance estimates for O, Mg and Si are in excellent agreement with those deduced from other studies, whilst the low estimate for C may reflect the use of the \ion{C}{ii} doublet at 4267\AA. The N estimates are approximately an order of magnitude greater than those found in unevolved B-type stars or \ion{H}{ii} regions but are consistent with the other estimates in AB-type supergiants. These results have been combined with results from a unified model atmosphere analysis of UVES/VLT spectra of B-type supergiants (Trundle et al. \cite{Tru04}) to discuss the evolutionary status of these objects. For two stars that are in common with those discussed by Trundle et al., we have undertaken a careful comparison in order to try to understand the relative importance of the different uncertainties present in such analyses, including observational errors and the use of static or unified models. We find that even for these relatively luminous supergiants, {\sc tlusty} models yield atmospheric parameters and chemical compositions similar to those deduced from the unified code {\sc fastwind}.

We present a comprehensive analysis of high-resolution, far-UV HST/STIS, FUSE, and optical spectra of 17 O stars in the SMC. Our analysis is based on NLTE metal line-blanketed model atmospheres calculated with our NLTE code TLUSTY. We systematically explore the sensitivity of various UV and optical lines to different stellar parameters. We have obtained consistent fits of the UV and the optical spectrum to derive the effective temperature, surface gravity, surface composition, and microturbulent velocity of each star. Stellar radii, masses, luminosities and ages then follow. Similarly to more limited recent studies, we derive cooler temperatures than the standard Teff calibration of O stars. We propose a new calibration between the spectral type and effective temperature based on our results from UV metal lines as well as optical hydrogen and helium lines. For stars of the same spectral subtype, we find a general good agreement between Teff determinations obtained with TLUSTY, CMFGEN, and FASTWIND models. We derive ionizing luminosities that are smaller by a factor of 3 compared to luminosities inferred from previous standard calibrations. The chemical composition analysis reveals that the surface of about 3/4 of the program stars is moderately to strongly enriched in nitrogen, while showing the original helium, carbon, and oxygen abundances. Our results support the new stellar evolution models that predict that the surface of fast rotating stars becomes N-rich during the main sequence phase because of rotationally-induced mixing. Most stars exhibit the ``mass discrepancy'' problem. This discrepancy too is a result of fast rotation which lowers the measured effective gravity. Our study thus emphasizes the importance of rotation in our understanding of the properties of massive stars. (abridged)

We present a non-LTE model atmosphere analysis of Chandra HRC-S/LETG and XMM-Newton RGS spectroscopy of the prototypical supersoft source CAL 83 in the Large Magellanic Cloud. Taken with a 16-month interval, the Chandra and XMM-Newton spectra are very similar. They reveal a very rich absorption line spectrum from the hot white dwarf photosphere, but no spectral signatures of a wind. We also report a third X-ray off-state during a later Chandra observation, demonstrating the recurrent nature of CAL 83. Moreover, we found evidence of short-timescale variability in the soft X-ray spectrum. We completed the analysis of the LETG and RGS spectra of CAL 83 with new NLTE line-blanketed model atmospheres that explicitly include 74 ions of the 11 most abundant species. We successfully matched the Chandra and XMM-Newton spectra assuming a model composition with LMC metallicity. We derived the basic stellar parameters of the hot white dwarf, but the current state of atomic data in the soft X-ray domain precludes a detailed chemical analysis. We have obtained the first direct spectroscopic evidence that the white dwarf is massive (Mwd > 1 Msun). The short timescale of the X-ray off-states is consistent with a high white dwarf mass. Our analysis thus provides direct support for supersoft sources as likely progenitors of SN Ia.

We present an analysis of optical HST/STIS and HST/FOS spectroscopy of 6 blue stragglers found in the globular clusters M3, NGC6752 and NGC6397. These stars are a subsample of a set of ~50 blue stragglers and stars above the main sequence turn-off in four globular clusters which will be presented in an forthcoming paper. All but the 6 stars presented here can be well fitted with non-LTE model atmospheres. The 6 misfits, on the other hand, possess Balmer jumps which are too large for the effective temperatures implied by their Paschen continua. We find that our data for these stars are consistent with models only if we account for extra absorption of stellar Balmer photons by an ionized circumstellar disk. Column densities of HI and CaII are derived as are the the disks' thicknesses. This is the first time that a circumstellar disk is detected around blue stragglers. The presence of magnetically-locked disks attached to the stars has been suggested as a mechanism to lose the large angular momentum imparted by the collision event at the birth of these stars. The disks implied by our study might not be massive enough to constitute such an angular momentum sink, but they could be the leftovers of once larger disks.

We reinvestigate the UV spectrum of NGC1535 by means of recently developed fully line-blanketed non-LTE models. These new models account for the wind in spherical geometry while handling the atomic data in a very similar way to the TLUSTY code. This approach ensures at the same time realistic predictions of the photospheric absorption lines and of the emission lines formed in the wind. Our analysis confirms the results of previous studies. We derive T_* = 70kK, log M_dot = -7.8M_sun/year and v_infty = 2000km/s.

We present FUSE spectra of 3 He-rich sdB stars. Two of these stars, PG1544+488 and JL87, reveal extremely strong C III lines at 977 and 1176A, while the carbon lines are quite weak in the third star, LB1766. We have analyzed the FUSE data using TLUSTY NLTE line-blanketed model atmospheres, and find that PG1544+488 has a surface composition of 96% He, 2% C, and 1% N. JL87 shows a similar surface enrichment of carbon and nitrogen, but some significant fraction of hydrogen still remains in its atmosphere. LB1766 has a surface composition devoid of hydrogen and strongly depleted of carbon, indicating that its surface material has undergone CN-cycle processing. We interpret these observations with new evolutionary calculations which suggest that He-rich sdB stars with C-rich compositions are the progeny of stars which underwent a delayed He-core flash on the white-dwarf cooling curve. During such a flash the interior convection zone will penetrate into the H envelope, thereby mixing the envelope with the He- and C-rich core. Such `flash-mixed' stars will arrive on the extreme horizontal branch (EHB) with He- and C-rich surface compositions and will be hotter than the hottest canonical (i.e., unmixed) EHB stars. Two types of flash mixing are possible: `deep' and `shallow', depending on whether the H envelope is mixed deeply into the site of the He flash or only with the outer layers of the core. Based on both their stellar parameters and surface compositions, we suggest that PG1544+488 and JL87 are examples of `deep' and `shallow' flash mixing, respectively. Flash mixing may therefore represent a new evolutionary channel for producing the hottest EHB stars. However, flash mixing cannot explain the abundance pattern in LB1766, which remains a challenge to current evolutionary models.

We present FUSE spectra of three He-rich sdB stars. Two of these stars, PG1544+488 and JL87, reveal extremely strong C III lines, suggesting that they have mixed triple-alpha carbon from the deep interior out to their surfaces. Using TLUSTY NLTE line-blanketed model atmospheres, we find that PG1544+488 has a surface composition of 96% He, 2% C, and 1% N. JL87 shows a similar surface enrichment of C and N but still retains a significant amount of hydrogen. In contrast, the third star, LB1766, is devoid of hydrogen and strongly depleted of carbon, indicating that its surface material has undergone CN-cycle processing. We interpret these observations with new evolutionary calculations which suggest that He-rich sdB stars with C-rich compositions arise from a delayed helium-core flash on the white-dwarf cooling curve. During such a flash the interior convection zone will penetrate into the stellar envelope, thereby mixing the envelope with the He- and C-rich core. Such "flash-mixed" stars will arrive on the extreme horizontal branch (EHB) with He- and C-rich surface compositions and will be hotter than the hottest canonical EHB stars. Two types of flash mixing are possible: "deep" and "shallow", depending on whether the hydrogen envelope is mixed deeply into the site of the helium flash or only with the outer layers of the core. Based on both their stellar parameters and surface compositions, we suggest that PG1544+488 and JL87 are examples of "deep" and "shallow" flash mixing, respectively.

We present the results of a detailed UV and optical spectral analysis of the properties of 6 dwarf O-type stars in the SMC H II region NGC 346. Stellar parameters, chemical abundances, and wind parameters have been determined using NLTE line blanketed models calculated with the photospheric code, Tlusty, and with the wind code, CMFGEN. The results, in particular iron abundances, obtained with the two NLTE codes compare very favorably, demonstrating that basic photospheric parameters of O dwarfs can be reliably determined using NLTE static model atmospheres. The two NLTE codes require a microturbulent velocity to match the observed spectra. Our results hint at a decrease of the microturbulent velocity from early O stars to late O stars. Similarly to several recent studies of galactic, LMC and SMC stars, we derive effective temperatures lower than predicted from the widely-used relation between spectral type and Teff, resulting in lower stellar luminosities and lower ionizing fluxes. From evolutionary tracks in the HR diagram, we find an age of 3 10^6 years for NGC 346. A majority of the stars in our sample reveal CNO-cycle processed material at their surface during the MS stage, indicating thus fast stellar rotation and/or very efficient mixing processes. We obtain an overall metallicity, Z = 0.2 Zsun, in good agreement with other recent analyses of SMC stars. The derived mass loss rate of the three most luminous stars agrees with recent theoretical predictions. However, the three other stars of our sample reveal very weak wind signatures. We obtain mass loss rates that are significantly lower than 10^{-8} Msun/yr, which is below the predictions of radiative line-driven wind theory by an order of magnitude or more. (abridged version)

With the goal of producing a reliable set of model atoms and singly-ionized ions for use in building NLTE model atmospheres, we have combined measured energy levels, critically-compiled line transition probabilities, and resonance-averaged calculations of photoionization cross-sections. A majority of the elements from Li to Ca are considered, covering most of the important species in late-type atmospheres. These include elements which contribute free electrons and/or continuous opacity in the ultraviolet (e.g., Mg, and Si), as well as trace elements whose abundance determinations rely on ultraviolet lines (e.g., B from B I lines). The new data complement and, for the species in common, supersede a previous collection of model atoms originally designed for use in studies of early-type stars.

We have constructed a comprehensive grid of 690 metal line-blanketed, NLTE, plane-parallel, hydrostatic model atmospheres for the basic parameters appropriate to O-type stars. The OSTAR2002 grid considers 12 values of effective temperatures, 27500K < Teff < 55000K with 2500K steps, 8 surface gravities, 3.0 < log g< 4.75 with 0.25 dex steps, and 10 chemical compositions, from metal-rich relative to the Sun to metal-free. The lower limit of log g for a given effective temperature is set by an approximate location of the Eddington limit. The selected chemical compositions have been chosen to cover a number of typical environments of massive stars: the galactic center, the Magellanic Clouds, Blue Compact Dwarf galaxies like I Zw 18, and galaxies at high redshifts. The paper contains a description of the OSTAR2002 grid and some illustrative examples and comparisons. The complete OSTAR2002 grid is available at our website at http://tlusty.gsfc.nasa.gov.

We present UV images and spectra of the starburst galaxy I Zw 18, taken with the Space Telescope Imaging Spectrograph. The high spatial resolution of these data allows us to isolate clusters containing Wolf-Rayet stars of the subtype WC. Our far-UV spectra clearly show CIV 1548,1551 and HeII 1640 emission of WC stars in two clusters: one within the bright (NW) half of I Zw 18, and one on the outskirts of this region. The latter spectrum is unusual, because the CIV is seen only in emission, indicating a spectrum dominated by WC stars. These data also demonstrate that the HI column in I Zw 18 is strongly peaked in the fainter (SE) half of I Zw 18, with a column depth far larger than that reported in previous analyses.

Stars with very large mass loss on the red-giant branch can undergo the helium flash while descending the white-dwarf cooling curve. Under these conditions the flash convection zone will mix the hydrogen envelope with the hot helium- burning core. Such ``flash-mixed'' stars will arrive on the extreme horizontal branch (EHB) with helium- and carbon-rich envelopes and will lie at higher temperatures than the hottest canonical (i.e., unmixed) EHB stars. Flash mixing provides a new evolutionary channel for populating the hot end of the EHB and may explain the origin of the high gravity, helium-rich sdO and sdB stars.

Hot subluminous stars lying up to 0.7 mag below the extreme horizontal branch (EHB) are found in the UV color-magnitude diagrams of omega Cen and NGC 2808. Such stars are unexplained by canonical HB theory. In order to explore the origin of these subluminous stars, we evolved a set of low-mass stars from the main sequence through the helium-core flash to the HB for a wide range in the mass loss along the red-giant branch (RGB). Stars with the largest mass loss evolve off the RGB to high effective temperatures before igniting helium in their cores. Our results indicate that the subluminous EHB stars, as well as the gap within the EHB of NGC 2808, can be explained if these stars undergo a late helium-core flash on the white-dwarf cooling curve. Under these conditions the flash convection will penetrate into the stellar envelope, thereby mixing most, if not all, of the envelope hydrogen into the hot helium- burning interior. This phenomenon is analogous to the "born-again" scenario for producing hydrogen-deficient stars during a very late helium-shell flash. "Flash mixing" greatly enhances the envelope helium and carbon abundances and, as a result, leads to an abrupt increase in the HB effective temperature. We argue that the EHB gap in NGC 2808 is caused by this theoretically predicted dichotomy in the HB morphology. Using new helium- and carbon-rich stellar atmospheres, we show that the flash-mixed stars have the same reduced UV flux as the subluminous EHB stars. Moreover, we demonstrate that models without flash mixing lie, at most, ~0.1 mag below the EHB and hence fail to explain the observations. Flash mixing may also provide a new evolutionary channel for producing the high gravity, He-rich sdO and sdB stars.

We present a UV CMD spanning the hot horizontal branch (HB), blue straggler, and white dwarf populations of the globular cluster NGC 2808. These data were obtained with the Space Telescope Imaging Spectrograph (STIS). Although previous optical CMDs of NGC 2808 show a high temperature gap within the hot HB population, no such gap is evident in our UV CMD. Instead, we find a population of hot subluminous HB stars, an anomaly only previously reported for the globular cluster omega Cen. Our modeling indicates that the location of these subluminous stars in the UV CMD, as well as the high temperature gap along the HB in optical CMDs, can be explained if these stars underwent a late He-core flash while descending the white dwarf cooling curve. We show that the convection zone produced by such a late He flash will penetrate into the H envelope, thereby mixing H into the hot He-burning interior, where it is rapidly consumed. This phenomenon is analogous to the "born again" scenario for producing H-deficient stars following a late He-shell flash. The flash mixing of the envelope greatly enhances the envelope He and C abundances, and leads, in turn, to a discontinuous increase in the HB effective temperatures at the transition between canonical and flash-mixed stars. We argue that the hot HB gap is associated with this theoretically predicted dichotomy in the HB properties. Moreover, the changes in the emergent spectral energy distribution caused by these abundance changes are primarily responsible for explaining the hot subluminous HB stars. Although further evidence is needed to confirm that a late He-core flash can account for the subluminous HB stars and the hot HB gap, we demonstrate that an understanding of these stars requires the use of appropriate theoretical models for their evolution, atmospheres, and spectra.

We present new coronagraphic images of Beta Pictoris obtained with the Space Telescope Imaging Spectrograph (STIS) in September 1997. The high-resolution images (0.1") clearly detect the circumstellar disk as close as 0.75" to the star, corresponding to a projected radius of 15 AU. The images define the warp in the disk with greater precision and at closer radii to Beta Pic than do previous observations. They show that the warp can be modelled by the projection of two components: the main disk, and a fainter component, which is inclined to the main component by 4-5 degrees, and which extends only as far as ~4" from the star. We interpret the main component as arising primarily in the outer disk and the tilted component as defining the inner region of the disk. The observed properties of the warped inner disk are inconsistent with a driving force from stellar radiation. However, warping induced by the gravitational potential of one or more planets is consistent with the data. Using models of planet-warped disks constructed by Larwood and Papaloizou (1997), we derive possible masses of the perturbing object.

The recent spectral analysis of LBDS 53W091 by Spinrad and his collaborators has suggested that this red galaxy at z=1.55 is at least 3.5 Gyr old. This imposes an important constraint on cosmology, suggesting that this galaxy formed at z > 6.5, assuming recent estimates of cosmological parameters. We have performed chi^2 tests to the continuum of this galaxy using its UV spectrum and photometric data (RJHK). We have used the updated Yi models that are based on the Yale tracks. We find it extremely difficult to reproduce such large age estimates, under the assumption of the most probable input parameters. Using the same configuration as in Spinrad et al. (solar abundance models), our analysis suggests an age of approximately 1.4 -- 1.8 Gyr. The discrepancy between Spinrad et al.'s age estimate (based on the 1997 Jimenez models) and ours originates from the large difference in the model integrated spectrum: the Jimenez models are much bluer than the Yi models and the Bruzual \& Charlot (BC) models. Preliminary tests favor the Yi and BC models. The updated age estimate of LBDS 53W091 would suggest that this galaxy formed approximately at z=2-3.

We have analyzed the far UV-spectrum of the globular cluster M79. We show that the nearly Gaussian mass distribution of zero-age horizontal-branch stars, as derived by Dixon et al.(1996), is able to reproduce the far-UV Hopkins Utraviolet Telescope (HUT) spectrum, if there is a luminous UV-bright star of about an effective temperature of 9,500 K within the HUT entrance slit, or, more likely, if the horizontal branch morphology becomes considerably redder in the core of M79, as observed in some other centrally condensed globular clusters. Agreement between the synthetic and observed far-UV spectra for M79 would also be improved if the surface abundances of the heavy elements in the hot horizontal-branch stars were enhanced by radiative diffusion. Contrary to Dixon et al. (1996) we do not need extremely low gravities to reproduce the width of Ly alpha.

We present high-quality HST/GHRS spectra in the Hydrogen L alpha spectral region of Vega and Sirius-A. Thanks to the signal-to-noise ratio achieved in these observations and to the similarity of the two spectra, we found clear evidence of emission features in the low flux region, 1190-1222 A. These emission lines can be attributed unambiguously to Fe II and Cr II transitions. In this spectral range, silicon lines are observed in absorption. We built a series of non-LTE model atmospheres with different, prescribed temperature stratification in the upper atmosphere and treating Fe II with various degrees of sophistication in non-LTE. Emission lines are produced by the combined effect of the Schuster mechanism and radiative interlocking, and can be explained without the presence of a chromosphere. Silicon absorption lines and the L alpha profile set constraints on the presence of a chromosphere, excluding a strong temperature rise in layers deeper than tau_R ~ 10^{-4}.

An echelle spectrogram (R = 30,000) of the 2300-3100 A region in the ultraviolet spectrum of the F8V star 9 Comae is presented. The observation is used to calibrate features in the mid-ultraviolet spectra of similar stars according to age and metal content. In particular, the spectral break at 2640 A is interpreted using the spectral synthesis code SYNSPEC. We use this feature to estimate the time since the last major star formation episode in the z=1.55 early-type galaxy LBDS 53W091, whose rest frame mid-ultraviolet spectrum, observed with the Keck Telescope, is dominated by the flux from similar stars that are at or near the main-sequence turnoff in that system (Spinrad et al. 1997). Our result, 1 Gyr if the flux-dominating stellar population has a metallicity twice solar, or 2 Gyr for a more plausible solar metallicity, is significantly lower than the previous estimate and thereby relaxes constraints on cosmological parameters that were implied by the earlier work.