Nebular phase spectra of core-collapse supernovae (SNe) provide critical and unique information on the progenitor massive star and its explosion. We present a set of 1-D steady-state non-local thermodynamic equilibrium radiative transfer calculations of type II SNe at 300d after explosion. Guided by results for a large set of stellar evolution simulations, we craft ejecta models for type II SNe from the explosion of a 12, 15, 20, and 25Msun star. The ejecta density structure and kinetic energy, the 56Ni mass, and the level of chemical mixing are parametrized. Our model spectra are sensitive to the adopted line Doppler width, a phenomenon we associate with the overlap of FeII and OI lines with Lyalpha and Lybeta. Our spectra show a strong sensitivity to 56Ni mixing since it determines where decay power is absorbed. Even at 300d after explosion, the H-rich layers reprocess the radiation from the inner metal rich layers. In a given progenitor model, variations in 56Ni mass and distribution impact the ejecta ionization, which can modulate the strength of all lines. Such ionization shifts can quench CaII line emission. In our set of models, the OI6300 doublet strength is the most robust signature of progenitor mass. However, we emphasize that convective shell merging in the progenitor massive star interior can pollute the O-rich shell with Ca, which will weaken the OI6300 doublet flux in the resulting nebular SN II spectrum. This process may occur in Nature, with a greater occurrence in higher mass progenitors, and may explain in part the preponderance of progenitor masses below 17Msun inferred from nebular spectra.

We present HST/STIS observations and analysis of two prominent nebular structures around the central source of Eta Carinae, the knots C and D. The former is brighter than the latter for emission lines from intermediate or high ionization potential ions. The brightness of lines from intermediate and high ionization potential ions significantly decreases at phases around periastron. We do not see conspicuous changes in the brightness of lines from low ionization potential (<13.6 eV) that the total extinction towards the Weigelt structures is that the total extinction towards the Weigelt structures is AsubV =2/0. that the total extinction towards the Weigelt structures is AV = 2.0. Weigelt C and D are characterized by an electron density of that the total extinction towards the Weigelt structures is AV = 2.0. Weigelt C and D are characterized by an electron density of 10exp6.9 cm-3 that does not significantly change throughout the orbital cycle. The electron temperature varies from 5500 K (around periastron) to 7200 K (around apastron). The relative changes in the brightness of He I lines are well reproduced by the variations in the electron temperature alone. We found that, at phases around periastron, the electron temperature seems to be higher for Weigelt C than that of D. The Weigelt structures are located close to the Homunculus equatorial plane, at a distance of about 1240 AU from the central source. From the analysis of proper motion and age, the Weigelt complex can be associated with the equatorial structure called the Butterfly Nebula surrounding the central binary system.

Eta Car is one of the most luminous and massive stars in our Galaxy and is the brightest mid-infrared (mid-IR) source in the sky, outside our solar system. Since the late 1990s the central source has dramatically brightened at ultraviolet and optical wavelengths. This might be explained by a decrease in circumstellar dust extinction. We aim to establish the mid-IR flux evolution and further our understanding of the star's ultraviolet and optical brightening. Mid-IR images from $8-20~\mu$m were obtained in 2018 with VISIR at the Very Large Telescope. Archival data from 2003 and 2005 are retrieved from the ESO Science Archive Facility and historical records are collected from publications. We present the highest angular resolution mid-IR images of $\eta$ Car to date at the corresponding wavelengths ($\geq 0.22''$). We reconstruct the mid-IR evolution of the spectral energy distribution of the spatially integrated Homunculus nebula from 1968 to 2018 and find no long-term changes. Eta Car's bolometric luminosity has been stable over the past five decades. We do not observe a long-term decrease in the mid-IR flux densities that could be associated with the brightening at ultraviolet and optical wavelengths, but circumstellar dust must be declining in our line-of-sight only. Short-term flux variations within about 25% of the mean levels could be present.

Hydrogen-rich (type II) supernovae (SNe) exhibit considerable photometric and spectroscopic diversity. Extending previous work that focused exclusively on photometry, we simultaneously model the multi-band light curves and optical spectra of Type II SNe using RSG progenitors that are characterized by their H-rich envelope masses or the mass and extent of an enshrouding cocoon at the star's surface. Reducing the H-rich envelope mass yields faster declining light curves, a shorter duration of the photospheric phase, broader line profiles at early times, but only a modest boost in early-time optical brightness. Increasing the mass of the circumstellar material (CSM) is more effective at boosting the early-time brightness and producing a fast-declining light curve while leaving the duration of the photospheric phase intact. It also makes the optical color bluer, delays the onset of recombination, and can severely reduce the speed of the fastest ejecta material. The early ejecta interaction with CSM is conducive to producing featureless spectra at $10-20$ d and a weak or absent H$\alpha$ absorption during the recombination phase. The slow decliners SNe 1999em, 2012aw, and 2004et can be explained with a $1.2 \times 10^{51}$ erg explosion in a compact ($\sim$600 R$_\odot$) RSG star from a 15 M$_\odot$ stellar evolution model. A small amount of CSM ($<0.2$ M$_\odot$) improves the match to the SN photometry at times $<$ 10 d. With more extended RSG progenitors, one predicts lower ejecta kinetic energies, but the SN color stays blue for too long and the spectral line widths are too narrow. The fast decliners SNe 2013ej and 2014G may require $0.5-1.0$ M$_\odot$ of CSM, although this depends on the CSM structure. A larger boost to the luminosity (as for fast decliners SNe 1979C or 1998S) requires interaction with a more spatially extended CSM, which might also be detached from the star.

In this study, we present one-dimensional, non-local-thermodynamic-equilibrium, radiative transfer simulations (using CMFGEN) in which we introduce micro-clumping at nebular times into two Type Ia supernova ejecta models. We use one sub-Chandrasekhar (sub-M$_{\rm Ch}$) ejecta with 1.02 M$_\odot$ and one M$_{\rm Ch}$ ejecta model with 1.40 M$_\odot$. We introduce clumping factors $f=$0.33,0.25, and 0.10 which are constant throughout the ejecta and compared to the unclumped $f=1.0$ case. We find that clumping is a natural mechanism to reduce the ionization of the ejecta, reducing emission from [Fe III], [Ar III], and [S III] by a factor of a few. For decreasing values of the clumping factor $f$, the [Ca II] $\lambda\lambda$7291,7324 doublet became a dominant cooling line for our M$_{\rm Ch}$ model but still weak in our sub-M$_{\rm Ch}$ model. Strong [Ca II] $\lambda\lambda$7291,7324 indicates non-thermal heating in that region and may constrain explosion modelling. Due to the low abundance of stable nickel, our sub-M$_{\rm Ch}$ model never showed the [Ni II] 1.939 micron diagnostic feature for all clumping values.

We present a new relativistic radiative-transfer code for $\gamma$-rays of energy less than 5 MeV in supernova (SN) ejecta. This code computes the opacities, the prompt emissivity (i.e. decay), and the scattering emissivity, and solves for the intensity in the co-moving frame. Because of the large expansion velocities of SN ejecta, we ignore redistribution effects associated with thermal motions. The energy deposition is calculated from the energy removed from the radiation field by scattering or photoelectric absorption. This new code yields comparable results to an independent Monte Carlo code. However, both yield non-trivial differences with the results from a pure absorption treatment of $\gamma$-ray transport. A synthetic observer's frame spectrum is also produced from the CMF intensity. At early times when the optical depth to $\gamma$-rays is large, the synthetic spectrum show asymmetric line profiles with redshifted absorption as seen in SN 2014J. This new code is integrated within CMFGEN and allows for an accurate and fast computation of the decay energy deposition in SN ejecta.

Much controversy surrounds the inferred progenitor masses of Type II-Plateau (II-P) supernovae (SNe). The debate is nourished by the discrepant results from radiation-hydrodynamics simulations, from pre-explosion imaging, and from studies of host stellar populations. Here, we present a controlled experiment using four solar metallicity models with zero-age main-sequence masses of 12, 15, 20, and 25Msun. Because of the effects of core burning and surface mass loss, these models reach core collapse as red-supergiant (RSG) stars with a similar H-rich envelope mass of 8 to 9Msun but with final masses in the range 11 to 16Msun. We explode the progenitors using a thermal bomb, adjusting the energy deposition to yield an asymptotic ejecta kinetic energy of 1.25 x 10^51 erg and an initial 56Ni mass of 0.04Msun. The resulting SNe produce similar photometric and spectroscopic properties from 10 to 200d. The spectral characteristics are degenerate. The scatter in early-time color results from the range in progenitor radii, while the differences in late-time spectra reflect the larger oxygen yields in more massive progenitors. Because the progenitors have a comparable H-rich envelope mass, the photospheric phase duration is comparable for all models; the difference in He-core mass is invisible. As different main-sequence masses can produce progenitors with a similar H-rich envelope mass, light curve modeling cannot provide a robust and unique solution for the ejecta mass of Type II-P SNe. The numerous uncertainties in massive star evolution and wind mass loss also prevent a robust association with a main-sequence star mass. Light curve modeling can at best propose compatibility.

The interacting binary Eta Carinae remains one of the most enigmatic massive stars in our Galaxy despite over four centuries of observations. In this work, its light curve from the ultraviolet to the near-infrared is analysed using spatially resolved HST observations and intense monitoring at the La Plata Observatory, combined with previously published photometry. We have developed a method to separate the central stellar object in the ground-based images using HST photometry and applying it to the more numerous ground-based data, which supports the hypothesis that the central source is brightening faster than the almost-constant Homunculus. After detrending from long-term brightening, the light curve shows periodic orbital modulation ($\Delta V$ $\sim$ 0.6 mag) attributed to the wind-wind collision cavity as it sweeps around the primary star and it shows variable projected area to our line-of-sight. Two quasi-periodic components with time scales of 2-3 and 8-10 yr and low amplitude, $\Delta V$ $<$ 0.2 mag, are superimposed on the brightening light-curve, being the only stellar component of variability found, which indicates minimal stellar instability. Moreover, the light curve analysis shows no evidence of `shell ejections' at periastron. We propose that the long-term brightening of the stellar core is due to the dissipation of a dusty clump in front of the central star, which works like a natural coronagraph. Thus, the central stars appear to be more stable than previously thought since the dominant variability originates from a changing circumstellar medium. We predict that the brightening phase, due mainly to dust dissipation, will be completed around 2032 $\pm$ 4 yr, when the star will be brighter than in the 1600's by up to $\Delta V$ $\sim$ 1 mag.

Since the discovery of SN (supernova) 1987A, the number of Type II-peculiar SNe has grown, revealing a rich diversity in photometric and spectroscopic properties. In this study, using a single 15Msun low-metallicity progenitor that dies as a blue supergiant (BSG), we have generated explosions with a range of energies and 56Ni masses. We then performed the radiative transfer modeling with CMFGEN from 1d until 300d after explosion. Our models yield light curves that rise to optical maximum in ~100d, with a similar brightening rate, and with a peak absolute V-band magnitude spanning from -14 to -16.5mag. All models follow a similar color evolution, entering the recombination phase within a few days of explosion, and reddening further until the nebular phase. Their spectral evolution is analogous, mostly differing in line profile width. With this model set, we study the Type II-pec SNe 1987A, 2000cb, 2006V, 2006au, 2009E, and 2009mw. Their photometric and spectroscopic diversity suggest that there is no prototypical Type II-pec SN. These SNe brighten to maximum faster than our model set, except perhaps SN2009mw. The spectral evolution of SN1987A conflicts with other observations and with model predictions from 20d until maximum: Halpha narrows and weakens while BaII lines strengthen faster than predicted, which we interpret as signatures of clumping. SN2000cb rises to maximum in only 20d and shows weak BaII lines. Its spectral evolution is well matched by an energetic ejecta but the light curve may require asymmetry. The persistent blue color, narrow lines, and weak Halpha absorption, seen in SN2006V conflicts with expectations for a BSG explosion powered by 56Ni and may require an alternative power source. In addition to diversity arising from different BSG progenitors, we surmise that their ejecta are asymmetric, clumped, and, in some cases, not solely powered by 56Ni decay [abridged].

We present interferometric observations of six O-type stars that were made with the Precision Astronomical Visible Observations (PAVO) beam combiner at the Center for High Angular Resolution Astronomy (CHARA) Array. The observations include multiple brackets for three targets, $\lambda$~Ori~A, $\zeta$~Oph, and 10~Lac, but there are only preliminary, single observations of the other three stars, $\xi$~Per, $\alpha$~Cam, and $\zeta$~Ori~A. The stellar angular diameters range from 0.55 milliarcsec for $\zeta$~Ori~A down to 0.11 mas for 10~Lac, the smallest star yet resolved with the CHARA Array. The rotational oblateness of the rapidly rotating star $\zeta$ Oph is directly measured for the first time. We assembled ultraviolet to infrared flux measurements for these stars, and then derived angular diameters and reddening estimates using model atmospheres and an effective temperature set by published results from analysis of the line spectrum. The model-based angular diameters are in good agreement with observed angular diameters. We also present estimates for the effective temperatures of these stars derived by setting the interferometric angular size and fitting the spectrophotometry.

There is both observational and theoretical evidence that the ejecta of core-collapse supernovae (SNe) are structured. Rather than being smooth and homogeneous, the material is made of over-dense and under-dense regions of distinct composition. Here, we explore the effect of clumping on the SN radiation during the photospheric phase using 1-D non-local thermodynamic equilibrium radiative transfer and an ejecta model arising from a blue-supergiant explosion (yielding a Type II-peculiar SN). Neglecting chemical segregation, we adopt a velocity-dependent volume-filling factor approach that assumes that the clumps are small but does not change the column density along any sightline. We find that clumping boosts the recombination rate in the photospheric layers, leading to a faster recession of the photosphere, an increase in bolometric luminosity, and a reddening of the SN colors through enhanced blanketing. The SN bolometric light curve peaks earlier and transitions faster to the nebular phase. On the rise to maximum, the strongest luminosity contrast between our clumped and smooth models is obtained at the epoch when the photosphere has receded to ejecta layers where the clumping factor is only 0.5 -- this clumping factor may be larger in Nature. Clumping is seen to have a similar influence in a Type II-Plateau SN model. As we neglect both porosity and chemical segregation our models underestimate the true impact of clumping. These results warrant further study of the influence of clumping on the observables of other SN types during the photospheric phase.

$\eta$ Car is a massive, eccentric binary with a rich observational history. We obtained the first high-cadence, high-precision light curves with the BRITE-Constellation nanosatellites over 6 months in 2016 and 6 months in 2017. The light curve is contaminated by several sources including the Homunculus nebula and neighboring stars, including the eclipsing binary CPD$-$59$^\circ$2628. However, we found two coherent oscillations in the light curve. These may represent pulsations that are not yet understood but we postulate that they are related to tidally excited oscillations of $\eta$ Car's primary star, and would be similar to those detected in lower-mass eccentric binaries. In particular, one frequency was previously detected by van Genderen et al. and Sterken et al. through the time period of 1974 to 1995 through timing measurements of photometric maxima. Thus, this frequency seems to have been detected for nearly four decades, indicating that it has been stable in frequency over this time span. These pulsations could help provide the first direct constraints on the fundamental parameters of the primary star if confirmed and refined with future observations.

We report on the first detection of a global change in the X-ray emitting properties of a wind-wind collision, thanks to XMM-Newton observations of the massive SMC system HD5980. While its lightcurve had remained unchanged between 2000 and 2005, the X-ray flux has now increased by a factor of ~2.5, and slightly hardened. The new observations also extend the observational coverage over the entire orbit, pinpointing the lightcurve shape. It has not varied much despite the large overall brightening, and a tight correlation of fluxes with orbital separation is found, without any hysteresis effect. Moreover, the absence of eclipses and of absorption effects related to orientation suggests a large size for the X-ray emitting region. Simple analytical models of the wind-wind collision, considering the varying wind properties of the eruptive component in HD5980, are able to reproduce the recent hardening and the flux-separation relationship, at least qualitatively, but they predict a hardening at apastron and little change in mean flux, contrary to observations. The brightness change could then possibly be related to a recently theorized phenomenon linked to the varying strength of thin-shell instabilities in shocked wind regions.

While Chandrasekhar-mass ($M_{\rm Ch}$) models with a low $^{56}\rm{Ni}$ yield can match the peak luminosities of fast-declining, 91bg-like Type Ia supernovae (SNe Ia), they systematically fail to reproduce their faster light-curve evolution. Here we illustrate the impact of a low ejecta mass on the radiative display of low-luminosity SNe Ia, by comparing a sub-$M_{\rm Ch}$ model resulting from the pure central detonation of a C-O White Dwarf (WD) to a $M_{\rm Ch}$ delayed-detonation model with the same $^{56}\rm{Ni}$ yield of 0.12 M$_{\odot}$. Our sub-$M_{\rm Ch}$ model from a 0.90 M$_{\odot}$ WD progenitor has a $\sim$5 day shorter rise time in the integrated UV-optical-IR (uvoir) luminosity, as well as in the $B$-band, and a $\sim$20 per cent higher peak uvoir luminosity ($\sim$1 mag brighter peak $M_B$). This sub-$M_{\rm Ch}$ model also displays bluer maximum-light colours due to the larger specific heating rate, and larger post-maximum uvoir and $B$-band decline rates. The luminosity decline at nebular times is also more pronounced, reflecting the enhanced escape of gamma rays resulting from the lower density of the progenitor WD. The deficit of stable nickel in the innermost ejecta leads to a notable absence of forbidden lines of [Ni II] in the nebular spectra. In contrast, the $M_{\rm Ch}$ model displays a strong line due to [Ni II] 1.939 $\mu\rm{m}$, which could in principle serve to distinguish between different progenitor scenarios. Our sub-$M_{\rm Ch}$ model offers an unprecedented agreement with optical and near-infrared observations of the 91bg-like SN 1999by, making a strong case for a WD progenitor significantly below the Chandrasekhar-mass limit for this event and other low-luminosity SNe Ia.

We present one-dimensional non-local thermodynamic equilibrium time-dependent radiative transfer simulations (using CMFGEN) of two sub-Chandrasekhar (sub-$M_{\rm Ch}$), one $M_{\rm Ch}$ and one super-$M_{\rm Ch}$ Type Ia SN ejecta models. Three originate from $M_{\rm Ch}$ delayed detonation models, and the fourth is a sub-$M_{\rm Ch}$ detonation model. Ejecta masses are 1.02, 1.04, 1.40, and 1.70 M$_\odot$, and all models have 0.62 M$_\odot$ of $^{56}{\rm Ni}$. Sub-$M_{\rm Ch}$ model light curves evolve faster, reaching bolometric maximum 2--3 days earlier and having 3--4 days shorter bolometric half light widths. The models vary by $\sim$12 per cent at maximum bolometric luminosity and by 0.17 mag in $B_{\rm max}$. While $\Delta M_{15}(B)$ increases with ejecta mass it only varies by $\sim$5 per cent around 1 mag. Sub-$M_{\rm Ch}$ models are 0.25 mag bluer in $B-R$ at $B_{\rm max}$. Optical spectra share many similarities, but lower mass models exhibit less UV line blanketing during the photospheric phase. At nebular times, significant NIR spectroscopic differences are seen. In particular, emission lines of the Ca II NIR triplet; [S III] $\lambda\lambda$9068,9530; [Ca II] $\lambda\lambda$7291,7324; [Ar III] $\lambda\lambda$7135,7751; and [Ni II] 1.939 $\mu$m are stronger in higher mass models. The [Ni II] 1.939 $\mu$m line is absent in the sub-$M_{\rm Ch}$ detonation model, and provides a valuable potential tool to distinguish sub-$M_{\rm Ch}$ explosions from $M_{\rm Ch}$ explosions. In general, the nebular phase models are too highly ionized. We attribute this to the neglect of clumping and/or the distribution of intermediate mass and iron group elements. The two sub-$M_{\rm Ch}$ models, while exploded by different mechanisms, can be distinguished in the $J$ and $H$ bands at late times (e.g., $+200$ days).

The progenitors of low-luminosity Type II-Plateau supernovae (SNe II-P) are believed to be red supergiant (RSG) stars, but there is much disparity in the literature concerning their mass at core collapse and therefore on the main sequence. Here, we model the SN radiation arising from the low-energy explosion of RSG stars of 12, 25, and 27 M$_{\odot}$ on the main sequence and formed through single star evolution. Despite the narrow range in ejecta kinetic energy (2.5$-$4.2$\times$10$^{50}$ erg) in our model set, the SN observables from our three models are significantly distinct, reflecting the differences in progenitor structure (e.g., surface radius, H-rich envelope mass, He-core mass). Our higher mass RSG stars give rise to Type II SNe that tend to have bluer colors at early times, a shorter photospheric phase, and a faster declining $V$-band light curve (LC) more typical of Type II-linear SNe, in conflict with the LC plateau observed for low-luminosity SNe II. The complete fallback of the CO core in the low-energy explosions of our high mass RSG stars prevents the ejection of any ${}^{56}$Ni (nor any core O or Si), in contrast to low-luminosity SNe II-P, which eject at least 0.001 M$_{\odot}$ of ${}^{56}$Ni. In contrast to observations, type II SN models from higher mass RSGs tend to show an H$\alpha$ absorption that remains broad at late times (due to a larger velocity at the base of the H-rich envelope). In agreement with the analyses of pre-explosion photometry, we conclude that low-luminosity SNe II-P likely arise from low-mass rather than high-mass RSG stars.

Infrared observations of the dusty, massive Homunculus Nebula around the luminous blue variable $\eta$ Carinae are crucial to characterize the mass-loss history and help constrain the mechanisms leading to the Great Eruption. We present the 2.4 - 670 $\mu$m spectral energy distribution, constructed from legacy ISO observations and new spectroscopy obtained with the {\em{Herschel Space Observatory}}. Using radiative transfer modeling, we find that the two best-fit dust models yield compositions which are consistent with CNO-processed material, with iron, pyroxene and other metal-rich silicates, corundum, and magnesium-iron sulfide in common. Spherical corundum grains are supported by the good match to a narrow 20.2 $\mu$m feature. Our preferred model contains nitrides AlN and Si$_3$N$_4$ in low abundances. Dust masses range from 0.25 to 0.44 $M_\odot$ but $M_{\rm{tot}} \ge$ 45 $M_\odot$ in both cases due to an expected high Fe gas-to-dust ratio. The bulk of dust is within a 5$"$ $\times$ 7$"$ central region. An additional compact feature is detected at 390 $\mu$m. We obtain $L_{\rm{IR}}$ = 2.96 $\times$ 10$^6$ $L_\odot$, a 25\% decline from an average of mid-IR photometric levels observed in 1971-1977. This indicates a reduction in circumstellar extinction in conjunction with an increase in visual brightness, allowing 25-40\% of optical and UV radiation to escape from the central source. We also present an analysis of $^{12}$CO and $^{13}$CO $J = 5-4$ through $9-8$ lines, showing that the abundances are consistent with expectations for CNO-processed material. The [$^{12}$C~{\sc{ii}}] line is detected in absorption, which we suspect originates in foreground material at very low excitation temperatures.

The faster light-curve evolution of low-luminosity Type Ia supernovae (SNe Ia) suggests that they could result from the explosion of white dwarf (WD) progenitors below the Chandrasekhar mass ($M_{\rm Ch}$). Here we present 1D non-LTE time-dependent radiative transfer simulations of pure central detonations of carbon-oxygen WDs with a mass ($M_\rm{tot}$) between 0.88 M$_{\odot}$ and 1.15 M$_{\odot}$, and a $^{56}\rm{Ni}$ yield between 0.08 M$_{\odot}$ and 0.84 M$_{\odot}$. Their lower ejecta density compared to $M_{\rm Ch}$ models results in a more rapid increase of the luminosity at early times and an enhanced $\gamma$-ray escape fraction past maximum light. Consequently, their bolometric light curves display shorter rise times and larger post-maximum decline rates. Moreover, the higher $M(^{56}\rm{Ni})/M_\rm{tot}$ ratio at a given $^{56}\rm{Ni}$ mass enhances the temperature and ionization level in the spectrum-formation region for the less luminous models, giving rise to bluer colours at maximum light and a faster post-maximum evolution of the $B-V$ colour. For sub-$M_{\rm Ch}$ models fainter than $M_B\approx -18.5$ mag at peak, the greater bolometric decline and faster colour evolution lead to a larger $B$-band post-maximum decline rate, $\Delta M_{15}(B)$. In particular, all of our previously-published $M_{\rm Ch}$ models (standard and pulsational delayed detonations) are confined to $\Delta M_{15}(B) < 1.4$ mag, while the sub-$M_{\rm Ch}$ models with $M_\rm{tot}\lesssim 1$ M$_{\odot}$ extend beyond this limit to $\Delta M_{15}(B)\approx 1.65$ mag for a peak $M_B\approx -17$ mag, in better agreement with the observed width-luminosity relation (WLR). Regardless of the precise ignition mechanism, these simulations suggest that fast-declining SNe Ia at the faint end of the WLR could result from the explosion of WDs whose mass is significantly below the Chandrasekhar limit.

As part of a search for Wolf-Rayet (WR) stars in the Magellanic Clouds, we have discovered a new type of WR star in the Large Magellanic Cloud (LMC). These stars have both strong emission lines, as well as He ii and Balmer absorption lines and spectroscopically resemble a WN3 and O3V binary pair. However, they are visually too faint to be WN3+O3V binary systems. We have found nine of these WN3/O3s, making up ~6% of the population of LMC WRs. Using CMFGEN, we have successfully modeled their spectra as single stars and have compared the physical parameters with those of more typical LMC WNs. Their temperatures are around 100,000 K, a bit hotter than the majority of WN stars (by around 10,000 K) although a few hotter WNs are known. The abundances are what you would expect for CNO equilibrium. However, most anomalous are their mass-loss rates which are more like that of an O-type star than a WN star. While their evolutionary status is uncertain, their low mass-loss rates and wind velocities suggest that they are not products of homogeneous evolution. It is possible instead that these stars represent an intermediate stage between O stars and WNs. Since WN3/O3 stars are unknown in the Milky Way, we suspect that their formation depends upon metallicity, and we are investigating this further by a deep survey in M33, which posses a metallicity gradient.

Early-time observations of the Type II supernovae (SNe) 2013cu and 2013fs have revealed an interaction of ejecta with material near the star surface. Unlike the Type IIn SN2010jl, which interacts with a dense wind for ~1yr, the interaction ebbs after 2-3d, suggesting a dense and compact circumstellar envelope. Here, we use multi-group radiation-hydrodynamics and non-local-thermodynamic-equilibrium radiative transfer to explore the properties of red supergiant (RSG) star explosions embedded in a variety of dense envelopes. We consider the cases of an extended static atmosphere or a steady-state wind, adopting a range of mass loss rates. The shock-breakout signal, the SN radiation up to 10d, and the ejecta dynamics are strongly influenced by the properties of this nearby environment. This compromises the use of early-time observations to constrain Rstar. The presence of narrow lines for 2-3d in 2013fs and 2013cu require a cocoon of material of ~0.01Msun out to 5-10Rstar. Spectral lines evolve from electron-scattering to Doppler broadened, with a growing blueshift of their emission peaks. Recent studies propose a super-wind phase with a mass loss rate from 0.001 up to 1Msun/yr in the last months/years of the RSG life, although there is no observational constraint that this external material is a steady-state outflow. Alternatively, observations may be explained by the explosion of a RSG star inside its complex atmosphere. Indeed, spatially resolved observations reveal that RSG stars have extended atmospheres, with the presence of downflows and upflows out to several Rstar, even in a standard RSG like Betelgeuse. Mass loading in the region intermediate between star and wind can accommodate the 0.01Msun needed to explain the observations of 2013fs. A puzzling super-wind phase prior to core-collapse may be therefore superfluous [Abridged].

Using 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations, we study the ejecta properties required to match the early and late-time photometric and spectroscopic properties of supernovae (SNe) associated with long-duration gamma-ray bursts (LGRBs). To match the short rise time, narrow light curve peak, and extremely broad spectral lines of SN1998bw requires a model with <3Msun ejecta but a high explosion energy of a few 1e52erg and 0.5Msun of Ni56. However the relatively high luminosity, the presence of narrow spectral lines of intermediate mass elements, and the low ionization at the nebular stage are matched with a more standard C-rich Wolf-Rayet (WR) star explosion, with an ejecta of >10Msun, an explosion energy >1e51erg, and only 0.1Msun of Ni56. As the two models are mutually exclusive, the breaking of spherical symmetry is essential to match the early/late photometric/spectroscopic properties of SN1998bw. This conclusion confirms the notion that the ejecta of SN1998bw is aspherical on large scales. More generally, with asphericity, the energetics and Ni56 mass of LGRB/SNe are reduced and their ejecta mass is increased, favoring a massive fast-rotating Wolf-Rayet star progenitor. Contrary to persisting claims in favor of the proto-magnetar model for LGRB/SNe, such progenitor/ejecta properties are compatible with collapsar formation. Ejecta properties of LGRB/SNe inferred from 1D radiative-transfer modeling are fundamentally flawed.

As part of a multi-year survey for Wolf-Rayet stars in the Magellanic Clouds, we have discovered a new type of Wolf-Rayet star with both strong emission and absorption. While one might initially classify these stars as WN3+O3V binaries based on their spectra, such a pairing is unlikely given their faint visual magnitudes. Spectral modeling suggests effective temperatures and bolometric luminosities similar to those of other early-type LMC WNs but with mass-loss rates that are three to five times lower than expected. They additionally retain a significant amount of hydrogen, with nitrogen at its CNO-equilibrium value (10x enhanced). Their evolutionary status remains an open question. Here we discuss why these stars did not evolve through quasi-homogeneous evolution. Instead we suggest that based on a link with long-duration gamma ray bursts, they may form in lower metallicity environments. A new survey in M33, which has a large metallicity gradient, is underway.

Supernova (SN) 2008bk is a well observed low-luminosity Type II event visually associated with a low-mass red-supergiant progenitor. To model SN 2008bk, we evolve a 12Msun star from the main sequence until core collapse, when it has a total mass of 9.88Msun, a He-core mass of 3.22Msun, and a radius of 502Rsun. We then artificially trigger an explosion that produces 8.29Msun of ejecta with a total energy of 2.5x10^50erg and ~0.009Msun of 56Ni. We model the subsequent evolution of the ejecta with non-Local-Thermodynamic-Equilibrium time-dependent radiative transfer. Although somewhat too luminous and energetic, this model reproduces satisfactorily the multi-band light curves and multi-epoch spectra of SN 2008bk, confirming the suitability of a low-mass massive star progenitor. As in other low-luminosity SNe II, the structured Halpha profile at the end of the plateau phase is probably caused by BaII 6496.9A rather than asphericity. We discuss the sensitivity of our results to changes in progenitor radius and mass, as well as chemical mixing. A 15% increase in progenitor radius causes a 15% increase in luminosity and a 0.2mag V-band brightening of the plateau but leaves its length unaffected. An increase in ejecta mass by 10% lengthens the plateau by ~10d. Chemical mixing introduces slight changes to the bolometric light curve, limited to the end of the plateau, but has a large impact on colours and spectra at nebular times.

Context. The mass loss from massive stars is not understood well. Eta Car is a unique object for studying the massive stellar wind during the LBV phase. It is also an eccentric binary with a period of 5.54 yr. The nature of both stars is uncertain, although we know from X-ray studies that there is a wind-wind collision whose properties change with orbital phase. Methods. Observations of Eta Car were carried out with the ESO VLTI and the AMBER instrument between approximately five and seven months before the August 2014 periastron passage. Velocity-resolved aperture-synthesis images were reconstructed from the spectrally dispersed interferograms. Interferometric studies can provide information on the binary orbit, the primary wind, and the wind collision. Results. We present velocity-resolved aperture-synthesis images reconstructed in more than 100 different spectral channels distributed across the Br Gamma 2.166 micrometer emission line. The intensity distribution of the images strongly depends on wavelength. At wavelengths corresponding to radial velocities of approximately -140 to -376 km/s measured relative to line center, the intensity distribution has a fan-shaped structure. At the velocity of -277 km/s, the position angle of the symmetry axis of the fan is ~ 126 degree. The fan-shaped structure extends approximately 8.0 mas (~ 18.8 au) to the southeast and 5.8 mas (~ 13.6 au) to the northwest, measured along the symmetry axis at the 16% intensity contour. The shape of the intensity distributions suggests that the obtained images are the first direct images of the innermost wind-wind collision zone. Therefore, the observations provide velocity-dependent image structures that can be used to test three-dimensional hydrodynamical, radiative transfer models of the massive interacting winds of Eta Car.

We have monitored the massive binary star $\eta$ Carinae with the CTIO/SMARTS 1.5~m telescope and CHIRON spectrograph from the previous apastron passage of the system through the recent 2014.6 periastron passage. Our monitoring has resulted in a large, homogeneous data set with an unprecedented time-sampling, spectral resolving power, and signal-to-noise. This allowed us to investigate temporal variability previously unexplored in the system and discover a kinematic structure in the P Cygni absorption troughs of neutral helium wind lines. The features observed occurred prior to the periastron passage and are seen as we look through the trailing arm of the wind-wind collision shock cone. We show that the bulk of the variability is repeatable across the last five periastron passages, and that the absorption occurs in the inner 230 AU of the system. In addition, we found an additional, high-velocity absorption component super-imposed on the P Cygni absorption troughs that has been previously un-observed in these lines, but which bears resemblance to the observations of the He~I $\lambda$10830 \AA\ feature across previous cycles. Through a comparison of the current smoothed particle hydrodynamical simulations, we show that the observed variations are likely caused by instabilities in the wind-wind collision region in our line-of-sight, coupled with stochastic variability related to clumping in the winds.

The X-ray emission of $\eta$ Carinae shows multiple features at various spatial and temporal scales. The central constant emission (CCE) component is centred on the binary and arises from spatial scales much smaller than the bipolar Homunculus nebula, but likely larger than the central wind--wind collision region between the stars as it does not vary over the $\sim$2-3 month X-ray minimum when it can be observed. Using large-scale 3D smoothed particle hydrodynamics (SPH) simulations, we model both the colliding-wind region between the stars, and the region where the secondary wind collides with primary wind ejected from the previous periastron passage. The simulations extend out to one hundred semimajor axes and make two limiting assumptions (strong coupling and no coupling) about the influence of the primary radiation field on the secondary wind. We perform 3D radiative transfer calculations on the SPH output to synthesize the X-ray emission, with the aim of reproducing the CCE spectrum. For the preferred primary mass-loss rate $\dot{M}_A\approx8.5\times10^{-4}$ M$_\odot$ yr$^{-1}$, the model spectra well reproduce the observation as the strong- and no-coupling spectra bound the CCE observation for longitude of periastron $\omega\approx252^\circ$, and bound/converge on the observation for $\omega\approx90^\circ$. This suggests that $\eta$ Carinae has moderate coupling between the primary radiation and secondary wind, that both the region between the stars and the comoving collision on the backside of the secondary generate the CCE, and that the CCE cannot place constraints on the binary's line of sight. We also discuss comparisons with common X-ray fitting parameters.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875\AA\ at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772\AA\ line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$\mu$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.

Using radiation-hydrodynamics and radiative-transfer simulations, we explore the origin of the spectral diversity of interacting supernovae (SNe) of type IIn. We revisit SN1994W and investigate the dynamical configurations that can give rise to spectra with narrow lines at all times. We find that a standard ~10Msun 10^51erg SN ejecta ramming into a 0.4Msun dense CSM is inadequate for SN1994W, as it leads to the appearance of broad lines at late times. This structure, however, generates spectra that exhibit the key morphological changes seen in SN1998S. For SN1994W, we consider a completely different configuration, which involves the interaction at a large radius of a low mass inner shell with a high mass outer shell. Such a structure may arise in an 8-12Msun star from a nuclear flash (e.g., of Ne) followed within a few years by core collapse. Our simulations show that the large mass of the outer shell leads to the complete braking of the inner shell material, the formation of a slow dense shell, and the powering of a luminous SN IIn, even for a low inner shell energy. Early on, our model line profiles are typical of SNe IIn, exhibiting narrow cores and broad electron-scattering wings. As observed in SN1994W, they also remain narrow at late times. Our SN1994W model invokes two low energy ejections, both atypical of observed massive stars, and illustrates the diversity of configurations leading to SNe IIn. These results also highlight the importance of spectra to constrain the dynamical properties and understand the origin of SNe IIn.

{\eta} Carinae is an extremely massive binary system in which rapid spectrum variations occur near periastron. Most notably, near periastron the He II $\lambda 4686$ line increases rapidly in strength, drops to a minimum value, then increases briefly before fading away. To understand this behavior, we conducted an intense spectroscopic monitoring of the He II $\lambda 4686$ emission line across the 2014.6 periastron passage using ground- and space-based telescopes. Comparison with previous data confirmed the overall repeatability of EW(He II $\lambda 4686$), the line radial velocities, and the timing of the minimum, though the strongest peak was systematically larger in 2014 than in 2009 by 26%. The EW(He II $\lambda 4686$) variations, combined with other measurements, yield an orbital period $2022.7\pm0.3$ d. The observed variability of the EW(He II $\lambda 4686$) was reproduced by a model in which the line flux primarily arises at the apex of the wind-wind collision and scales inversely with the square of the stellar separation, if we account for the excess emission as the companion star plunges into the hot inner layers of the primary's atmosphere, and including absorption from the disturbed primary wind between the source and the observer. This model constrains the orbital inclination to $135^\circ$-$153^\circ$, and the longitude of periastron to $234^\circ$-$252^\circ$. It also suggests that periastron passage occurred on $T_0 = 2456874.4\pm1.3$ d. Our model also reproduced EW(He II $\lambda 4686$) variations from a polar view of the primary star as determined from the observed He II $\lambda 4686$ emission scattered off the Homunculus nebula.

We present a multi-wavelength (X-ray to optical) analysis, based on non-local thermodynamic equilibrium photospheric+wind models, of the B0 Ia-supergiant: $\epsilon$~Ori. The aim is to test the consistency of physical parameters, such as the mass-loss rate and CNO abundances, derived from different spectral bands. The derived mass-loss rate is $\dot{M}/\sqrt{f_\infty}\sim$1.6$\times$10$^{-6}$ M$_\odot$ yr$^{-1}$ where $f_\infty$ is the volume filling factor. However, the S IV $\lambda\lambda$1062,1073 profiles are too strong in the models; to fit the observed profiles it is necessary to use $f_\infty<$0.01. This value is a factor of 5 to 10 lower than inferred from other diagnostics, and implies $\dot{M} \lesssim1 \times 10^{-7}$ M$_\odot$ yr$^{-1}$. The discrepancy could be related to porosity-vorosity effects or a problem with the ionization of sulfur in the wind. To fit the UV profiles of N V and O VI it was necessary to include emission from an interclump medium with a density contrast ($\rho_{cl}/\rho_{ICM}$) of $\sim$100. X-ray emission in H-He like and Fe L lines was modeled using four plasma components located within the wind. We derive plasma temperatures from $1 \times 10^{6}$ to $7\times 10^{6}$ K, with lower temperatures starting in the outer regions (R$_0\sim$3-6 R$_*$), and a hot component starting closer to the star (R$_0\lesssim$2.9 R$_*$). From X-ray line profiles we infer $\dot{M} <\, 4.9\times10^{-7}$ M$_\odot$ yr$^{-1}$. The X-ray spectrum ($\geq$0.1 kev) yields an X-ray luminosity $L_{\rm X}\sim 2.0\times10^{-7} L_{\rm bol}$, consistent with the superion line profiles. X-ray abundances are in agreement with those derived from the UV and optical analysis: $\epsilon$ Ori is slightly enhanced in nitrogen and depleted in carbon and oxygen, evidence for CNO processed material.

As part of our ongoing Wolf-Rayet (WR) Magellanic Cloud survey, we have discovered 13 new WRs. However, the most exciting outcome of our survey is not the number of new WRs, but their unique characteristics. Eight of our discoveries appear to belong to an entirely new class of WRs. While one might naively classify these stars as WN3+O3V binaries, such a pairing is unlikely. Preliminary CMFGEN modeling suggests physical parameters similar to early-type WNs in the Large Magellanic Cloud except with mass-loss rates three to five times lower and slightly higher temperatures. The evolution status of these stars remains an open question.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".

The light curves of Type Ia supernovae (SNe Ia) are powered by the radioactive decay of $^{56}$Ni to $^{56}$Co at early times, and the decay of $^{56}$Co to $^{56}$Fe from ~60 days after explosion. We examine the evolution of the [Co III] 5892 A emission complex during the nebular phase for SNe Ia with multiple nebular spectra and show that the line flux follows the square of the mass of $^{56}$Co as a function of time. This result indicates both efficient local energy deposition from positrons produced in $^{56}$Co decay, and long-term stability of the ionization state of the nebula. We compile 77 nebular spectra of 25 SN Ia from the literature and present 17 new nebular spectra of 7 SNe Ia, including SN2014J. From these we measure the flux in the [Co III] 5892 A line and remove its well-behaved time dependence to infer the initial mass of $^{56}$Ni ($M_{Ni}$) produced in the explosion. We then examine $^{56}$Ni yields for different SN Ia ejected masses ($M_{ej}$ - calculated using the relation between light curve width and ejected mass) and find the $^{56}$Ni masses of SNe Ia fall into two regimes: for narrow light curves (low stretch s~0.7-0.9), $M_{Ni}$ is clustered near $M_{Ni}$ ~ 0.4$M_\odot$ and shows a shallow increase as $M_{ej}$ increases from ~1-1.4$M_\odot$; at high stretch, $M_{ej}$ clusters at the Chandrasekhar mass (1.4$M_\odot$) while $M_{Ni}$ spans a broad range from 0.6-1.2$M_\odot$. This could constitute evidence for two distinct SN Ia explosion mechanisms.

We present numerical simulations that include 1-D Eulerian multi-group radiation-hydrodynamics, 1-D non-LTE radiative transfer, and 2-D polarised radiative transfer for super-luminous interacting supernovae (SNe). Our reference model is a ~10Msun inner shell with 10^51erg ramming into a ~3Msun cold outer shell (the circumstellar-medium, or CSM) that extends from 10^15cm to 2x10^16cm and moves at 100km/s. We discuss the light curve evolution, which cannot be captured adequately with a grey approach. In these interactions, the shock-crossing time through the optically-thick CSM is much longer than the photon diffusion time. Radiation is thus continuously leaking from the shock through the CSM, in disagreement with the shell-shocked model that is often invoked. Our spectra redden with time, with a peak distribution in the near-UV during the first month gradually shifting to the optical range over the following year. Initially Balmer lines exhibit a narrow line core and the broad line wings that are characteristic of electron scattering in the SNe IIn atmospheres (CSM). At later times they also exhibit a broad blue shifted component which arises from the cold dense shell. Our model results are broadly consistent with the bolometric light curve and spectral evolution observed for SN2010jl. Invoking a prolate pole-to-equator density ratio in the CSM, we can also reproduce the ~2% continuum polarisation, and line depolarisation, observed in SN2010jl. By varying the inner shell kinetic energy and the mass and extent of the outer shell, a large range of peak luminosities and durations, broadly compatible with super-luminous SNe IIn like 2010jl or 2006gy, can be produced.

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.

We present 1D non-local thermodynamic equilibrium (non-LTE) time-dependent radiative-transfer simulations of a Chandrasekhar-mass delayed-detonation model which synthesizes 0.51 Msun of 56Ni, and confront our results to the Type Ia supernova (SN Ia) 2002bo over the first 100 days of its evolution. Assuming only homologous expansion, this same model reproduces the bolometric and multi-band light curves, the secondary near-infrared (NIR) maxima, and the optical and NIR spectra. The chemical stratification of our model qualitatively agrees with previous inferences by Stehle et al., but reveals significant quantitative differences for both iron-group and intermediate-mass elements. We show that +/-0.1 Msun (i.e., +/-20 per cent) variations in 56Ni mass have a modest impact on the bolometric and colour evolution of our model. One notable exception is the U-band, where a larger abundance of iron-group elements results in less opaque ejecta through ionization effects, our model with more 56Ni displaying a higher near-UV flux level. In the NIR range, such variations in 56Ni mass affect the timing of the secondary maxima but not their magnitude, in agreement with observational results. Moreover, the variation in the I, J, and K_s magnitudes is less than 0.1 mag within ~10 days from bolometric maximum, confirming the potential of NIR photometry of SNe Ia for cosmology. Overall, the delayed-detonation mechanism in single Chandrasekhar-mass white dwarf progenitors seems well suited for SN 2002bo and similar SNe Ia displaying a broad Si II 6355 A line. Whatever multidimensional processes are at play during the explosion leading to these events, they must conspire to produce an ejecta comparable to our spherically-symmetric model.

We present non-LTE time-dependent radiative transfer simulations for ejecta produced by the detonation of an helium shell at the surface of a low-mass carbon/oxygen white dwarf (WD). This mechanism is one possible origin for supernovae (SNe) with faint and fast-decaying light curves, such as .Ia SNe and Ca-rich transients. Our initial ejecta conditions at 1d are given by the 0.18B explosion model COp45HEp2 of Waldman et al.. The 0.2Msun ejecta initially contains 0.11Msun of He, 0.03Msun of Ca, and 0.03Msun of Ti. We obtain a ~5d rise to a bolometric maximum of 3.59x10^41 erg/s, primarily powered by 48V decay. Multi-band light curves show distinct morphologies, with a rise to maximum magnitude (-14.3 to -16.7mag) that varies between 3 to 9d from the U to the K bands. Near-IR light curves show no secondary maximum. Because of the presence of both HeI and SiII lines at early times we obtain a hybrid Type Ia/Ib classification. During the photospheric phase line blanketing is caused primarily by TiII. At nebular times, the spectra show strong CaII lines in the optical (but no [OI] 6300--6364A emission), and TiII in the near-IR. Overall, these results match qualitatively the very disparate properties of .Ia SNe and Ca-rich transients. Although the strong TiII blanketing and red colors that we predict are rarely observed, they are seen, for example, in OGLE-2013- SN-079. Furthermore, we obtain a faster light-curve evolution than, for example, PTF10iuv, indicating an ejecta mass >0.2Msun. An alternate scenario may be the merger of two WDs, one or both composed of He.

Massive stars show a variety of spectral variability: presence of discrete absorption components in UV P-Cygni profiles, optical line profile variability, X-ray variability, radial velocity modulations. Our goal is to study the spectral variability of single OB stars to better understand the relation between photospheric and wind variability. For that, we rely on high spectral resolution, high signal-to-noise ratio optical spectra collected with the spectrograph NARVAL on the Telescope Bernard Lyot at Pic du Midi. We investigate the variability of twelve spectral lines by means of the Temporal Variance Spectrum (TVS). The selected lines probe the radial structure of the atmosphere, from the photosphere to the outer wind. We also perform a spectroscopic analysis with atmosphere models to derive the stellar and wind properties, and to constrain the formation region of the selected lines. We show that variability is observed in the wind lines of all bright giants and supergiants, on a daily timescale. Lines formed in the photosphere are sometimes variable, sometimes not. The dwarf stars do not show any sign of variability. If variability is observed on a daily timescale, it can also (but not always) be observed on hourly timescales, albeit with lower amplitude. There is a very clear correlation between amplitude of the variability and fraction of the line formed in the wind. Strong anti-correlations between the different part of the temporal variance spectrum are observed. Our results indicate that variability is stronger in lines formed in the wind. A link between photospheric and wind variability is not obvious from our study, since wind variability is observed whatever the level of photospheric variability. Different photospheric lines also show different degrees of variability.

New spectroscopic observations of the LBV/WR multiple system HD5980 in the Small Magellanic Cloud are used to address the question of the masses and evolutionary status of the two very luminous stars in the 19.3d eclipsing binary system. Two distinct components of the N V 4944 A line are detected in emission and their radial velocity variations are used to derive masses of 61 and 66 Mo, under the assumption that binary interaction effects on this atomic transition are negligible. We propose that this binary system is the product of quasi-chemically homogeneous evolution with little or no mass transfer. Thus, both of these binary stars may be candidates for gamma-ray burst progenitors or even pair instability supernovae. Analysis of the photospheric absorption lines belonging to the third-light object in the system confirm that it consists of an O-type star in a 96.56d eccentric orbit (e=0.82) around an unseen companion. The 5:1 period ratio and high eccentricities of the two binaries suggest that they may constitute a hierarchical quadruple system.

Our new survey for Wolf-Rayet stars in the Magellanic Clouds is only 15% complete but has already found 9 new WRs in the LMC. This suggests that the total WR population in the LMC may be underestimated by 10-40%. Eight of the nine are WNs, demonstrating that the "observed" WC to WN ratio is too large, and is biased towards WC stars. The ninth is another rare WO star, the second we have found in the LMC in the past two years. Five (and possibly six) of the 8 WNs are of a new class of WRs, which pose a significant challenge to our understanding. Naively we would classify these stars as "WN3+O3V," but there are several reasons why such a pairing is unlikely, not the least of which is that the absolute visual magnitudes of these stars are faint, with M_V~ -2.3 to -3.1. We have performed a preliminary analysis with CMFGEN, and we find that (despite the faint visual magnitudes) the bolometric luminosities of these stars are normal for early-type WNs. Our fitting suggests that these stars are evolved, with significantly enriched N and He. Their effective temperatures are also normal for early-type WNs. What is unusual about these stars is that they have a surprisingly small mass-loss rate compared to other early-type WNs. How these stars got to be the way they are (single star evolution? binary evolution?) remains an open question. For now, we are designating this class as WN3/O3, in analogy to the late-type WN "slash" stars.

Over the years, directed surveys and incidental spectroscopy have identified 12 Wolf-Rayet (WR) stars in the SMC and 139 in the LMC, numbers which are often described as "essentially} complete." Yet, new WRs are discovered in the LMC almost yearly. We have therefore initiated a new survey of both Magellanic Clouds using the same interference-filter imaging technique previously applied to M31 and M33. We report on our first observing season, in which we have successfully surveyed ~15$% of our intended area of the SMC and LMC. Spectroscopy has confirmed 9 newly found WRs in the LMC (a 6% increase), including one of WO-type, only the third known in that galaxy and the second to be discovered recently. The other eight are WN3 stars that include an absorption component. In two, the absorption is likely from an O-type companion, but the other six are quite unusual. Five would be classified naively as "WN3+O3 V," but such a pairing is unlikely given the rarity of O3 stars, the short duration of this phase (which is incommensurate with the evolution of a companion to a WN star), and because these stars are considerably fainter than O3 V stars. The sixth star may also fall into this category. CMFGEN modeling suggests these stars are hot, bolometrically luminous, and N-rich like other WN3 stars, but lack the strong winds that characterize WNs. Finally, we discuss two rare Of?p stars and four Of supergiants we found, and propose that the B[e] star HD 38489 may have a WN companion.

We explore the physics of SN Ia light curves and spectra using the 1-D non-LTE time-dependent radiative-transfer code CMFGEN. Rather than adjusting ejecta properties to match observations, we select as input one "standard" 1-D Chandrasekhar-mass delayed-detonation hydrodynamical model, and then explore the sensitivity of radiation and gas properties on radiative-transfer modeling assumptions. The correct computation of SN Ia radiation is not exclusively a solution to an "opacity problem", characterized by the treatment of a large number of lines. It is also key to treat important atomic processes consistently. Besides handling line blanketing in non-LTE, we show that including forbidden line transitions of metals is increasingly important for the temperature and ionization of the gas beyond maximum light. Non-thermal ionization and excitation are also critical since they affect the color evolution and the Delta-M15 of our model. While impacting little the bolometric luminosity, a more complete treatment of decay routes leads to enhanced line blanketing, e.g., associated with 48Ti in the U and B bands. Overall, we find that SN Ia radiation properties are influenced in a complicated way by the atomic data we employ, so that obtaining converged results is a challenge. We nonetheless obtain a good match to the golden standard type Ia SN 2005cf in the optical and near-IR, from 5 to 60d after explosion, suggesting that assuming spherical symmetry is not detrimental to SN Ia radiative-transfer modeling at these times. Multi-D effects no doubt matter, but they are perhaps less important than accurately treating non-LTE processes [abridged].

Observations of SN 2011fe at early times reveal an evolution analogous to a fireball model of constant color. In contrast, our unmixed delayed detonations of Chandrasekhar-mass white dwarfs (DDC series) exhibit a faster brightening concomitant with a shift in color to the blue. In this paper, we study the origin of these discrepancies. We find that strong chemical mixing largely resolves the photometric mismatch at early times, but it leads to an enhanced line broadening that contrasts, for example, with the markedly narrow SiII6355A line of SN 2011fe. We also explore an alternative configuration with pulsational-delayed detonations (PDDEL model series). Because of the pulsation, PDDEL models retain more unburnt carbon, have little mass at high velocity, and have a much hotter outer ejecta after the explosion. The pulsation does not influence the inner ejecta, so PDDEL and DDC models exhibit similar radiative properties beyond maximum. However, at early times, PDDEL models show bluer optical colors and a higher luminosity, even for weak mixing. Their early-time radiation is derived primarily from the initial shock-deposited energy in the outer ejecta rather than radioactive decay heating. Furthermore, PDDEL models show short-lived CII lines, reminiscent of SN 2013dy. They typically exhibit lines that are weaker, narrower, and of near-constant width, reminiscent of SN 2011fe. In addition to multi-dimensional effects, varying configurations for such ``pulsations" offer a source of spectral diversity amongst SNe Ia. PDDEL and DDC models also provide one explanation for low- and high-velocity gradient SNe Ia.

The high metal content and fast expansion of supernova (SN) Ia ejecta lead to considerable line overlap in their optical spectra. Uncertainties in composition and ionization further complicate the process of line identification. In this paper, we focus on the 5900A emission feature seen in SN Ia spectra after bolometric maximum, a line which in the last two decades has been associated with [CoIII]5888A or NaID. Using non-LTE time-dependent radiative-transfer calculations based on Chandrasekhar-mass delayed-detonation models, we find that NaID line emission is extremely weak at all post-maximum epochs. Instead, we predict the presence of [CoIII]5888A after maximum in all our SN Ia models, which cover a range from 0.12 to 0.87Msun of 56Ni. We also find that the [CoIII]5888A forbidden line is present within days of bolometric maximum, and strengthens steadily for weeks thereafter. Both predictions are confirmed by observations. Rather than trivial taxonomy, these findings confirm that it is necessary to include forbidden-line transitions in radiative-transfer simulations of SNe Ia, both to obtain the correct ejecta cooling rate and to match observed optical spectra.

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 explore the properties of Type II-Plateau (II-P) supernovae (SNe) together with their red-supergiant (RSG) star progenitors. Using MESA STAR, we modulate the parameters (e.g., mixing length, overshoot, rotation, metallicity) that control the evolution of a 15Msun main-sequence star to produce a variety of physical pre-SN models and SN II-P ejecta. We extend previous modeling of SN II-P radiation to include photospheric and nebular phases, as well as multi-band light curves and spectra. Our treatment does not assume local thermodynamic equilibrium, is time dependent, treats explicitly the effects of line blanketing, and incorporates non-thermal processes. We find that the color properties of SNe II-P require large model atoms for FeI and FeII, much larger than adopted in Dessart & Hillier (2011). The color properties also imply RSG progenitors of limited extent (~500Rsun) --- larger progenitor stars produce a SN II-P radiation that remains too blue for too long. This finding calls for a reduction of RSG radii, perhaps through a strengthening of convective energy transport in RSG envelopes. Increased overshoot and rotation reduce the ratio of ejecta to helium-core mass, similarly to an increase in main-sequence mass, and thus complicate the inference ofprogenitor masses. In contrast to the great sensitivity on progenitor radius, SN II-P color evolution appears insensitive to variations in explosion energy. Finally, we document the numerous SN II-P signatures that vary with progenitor metallicity, revealing their potential for metallicity determinations in the nearby and distant Universe.

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.

The model atmosphere programs FASTWIND and CMFGEN are both elegantly designed to perform non-LTE analyses of the spectra of hot massive stars, and include sphericity and mass-loss. The two codes differ primarily in their approach towards line blanketing, with CMFGEN treating all of the lines in the co-moving frame and FASTWIND taking an approximate approach which speeds up execution times considerably. Although both have been extensively used to model the spectra of O-type stars, no studies have used the codes to independently model the same spectra of the same stars and compare the derived physical properties. We perform this task on ten O-type stars in the Magellanic Clouds. For the late-type O supergiants, both CMFGEN and FASTWIND have trouble fitting some of the He I lines, and we discuss causes and cures. We find that there is no difference in the average effective temperatures found by the two codes for the stars in our sample, although the dispersion is large, due primarily to the various difficulties each code has with He I. The surface gravities determined using FASTWIND are systematically lower by 0.12 dex compared to CMFGEN, a result we attribute to the better treatment of electron scattering by CMFGEN. This has implications for the interpretation of the origin of the so-called mass discrepancy, as the masses derived by FASTWIND are on average lower than inferred from stellar evolutionary models, while those found by CMFGEN are in better agreement.

The delayed-detonation explosion mechanism applied to a Chandrasekhar-mass white dwarf offers a very attractive model to explain the inferred characteristics of Type Ia supernovae (SNe Ia). The resulting ejecta are chemically stratified, have the same mass and roughly the same asymptotic kinetic energy, but exhibit a range in 56Ni mass. We investigate the contemporaneous photometric and spectroscopic properties of a sequence of delayed-detonation models, characterized by 56Ni masses between 0.18 and 0.81 Msun. Starting at 1d after explosion, we perform the full non-LTE, time-dependent radiative transfer with the code CMFGEN, with an accurate treatment of line blanketing, and compare our results to SNe Ia at bolometric maximum. Despite the 1D treatment, our approach delivers an excellent agreement to observations. We recover the range of SN Ia luminosities, colours, and spectral characteristics from the near-UV to 1 micron, for standard as well as low-luminosity 91bg-like SNe Ia. Our models predict an increase in rise time to peak with increasing 56Ni mass, from ~15 to ~21d, yield peak bolometric luminosities that match Arnett's rule to within 10%, and reproduce the much smaller scatter in near-IR magnitudes compared to the optical. We reproduce the morphology of individual spectral features, the stiff dependence of the R(Si) spectroscopic ratio on 56Ni mass, and the onset of blanketing from TiII/ScII in low-luminosity SNe Ia with a 56Ni mass <0.3 Msun. We find that ionization effects, which often dominate over abundance variations, can produce high-velocity features in CaII lines, even in 1D. Distinguishing between different SN Ia explosion mechanisms is a considerable challenge but the results presented here provide additional support to the viability of the delayed-detonation model.

We present non-LTE time-dependent radiative-transfer simulations of pair-instability supernovae (PISNe) stemming from red-supergiant (RSG), blue-supergiant (BSG) and Wolf-Rayet (WR) star rotation-free progenitors born in the mass range 160-230Msun, at 10^-4 Zsun. Although subject to uncertainties in convection and stellar mass-loss rates, our initial conditions come from physically-consistent models that treat evolution from the main-sequence, the onset of the pair-production instability, and the explosion phase. With our set of input models characterized by large 56Ni and ejecta masses, and large kinetic energies, we recover qualitatively the Type II-Plateau, II-peculiar, and Ib/c light-curve morphologies, although they have larger peak bolometric luminosities (~10^9 to 10^10 Lsun) and a longer duration (~200d). We discuss the spectral properties for each model during the photospheric and nebular phases, including Balmer lines in II-P and II-pec at early times, the dominance of lines from intermediate-mass-elements (IMEs) near the bolometric maximum, and the strengthening of metal line blanketing thereafter. Having similar He-core properties, all models exhibit similar post-peak spectra that are strongly blanketed by FeII and FeI lines, characterized by red colors, and that arise from photospheres/ejecta with a temperature of <4000K. Combined with the modest line widths after bolometric peak, these properties contrast with those of known super-luminous SNe suggesting that PISNe are yet to be discovered. Being reddish, PISNe will be difficult to observe at high redshift except when they stem from RSG explosions, in which case they could be used as metallicity probes and distance indicators.

We analyze spatially resolved spectroscopic observations of the Eta Carinae binary system obtained with HST/STIS. Eta Car is enshrouded by the dusty Homunculus nebula, which scatters light emitted by the central binary and provides a unique opportunity to study a massive binary system from different vantage points. We investigate the latitudinal and azimuthal dependence of H$\alpha$ line profiles caused by the presence of a wind-wind collision (WWC) cavity created by the companion star. Using two-dimensional radiative transfer models, we find that the wind cavity can qualitatively explain the observed line profiles around apastron. Regions of the Homunculus which scatter light that propagated through the WWC cavity show weaker or no H alpha absorption. Regions scattering light that propagated through a significant portion of the primary wind show stronger P Cygni absorption. Our models overestimate the H alpha absorption formed in the primary wind, which we attribute to photoionization by the companion, not presently included in the models. We can qualitatively explain the latitudinal changes that occur during periastron, shedding light on the nature of Eta Car's spectroscopic events. Our models support the idea that during the brief period of time around periastron when the primary wind flows unimpeded toward the observer, H alpha absorption occurs in directions toward the central object and Homunculus SE pole, but not toward equatorial regions close to the Weigelt blobs. We suggest that observed latitudinal and azimuthal variations are dominated by the companion star via the WWC cavity, rather than by rapid rotation of the primary star.

Aims. We report new results for the massive evolved and main sequence members of the young galactic cluster DBS2003 179. We determine the physical parameters and investigate the high-mass stellar content of the cluster, as well as of its close vicinity. Methods. Our analysis is based on ISAAC/VLT moderate-resolution (R\approx4000) infrared spectroscopy of the brightest cluster members. We derive stellar parameters for sixteen of the stellar members, using full non-LTE modeling of the obtained spectra. Results. The cluster contains three late WN or WR/LBV stars (Obj 4, Obj 15, and Obj 20:MDM32) and at least 5 OIf and 5 OV stars. According to the Hertzsprung-Russell diagram for DBS2003 179, the WR stars show masses above 85Msun, the OIf stars are between 40 and 80Msun, and the main sequence O stars are >20Msun. There are indications of binarity for Obj 4 and Obj 11, and Obj 3 shows a variable spectrum. The cluster is surrounded by a continuous protostar formation region most probably triggered by DBS2003 179. Conclusions. We confirm that DBS2003 179 is young massive cluster (2.5 10^4Msun) very close to the Galactic center at the distance of 7.9+-0.8 kpc.

Much uncertainty surrounds the origin of super-luminous supernovae (SNe). Motivated by the discovery of the Type Ic SN2007bi, we study its proposed association with a pair-instability SN (PISN). We compute stellar-evolution models for primordial ~200Msun stars, simulating the implosion/explosion due to the pair-production instability, and use them as inputs for detailed non-LTE time-dependent radiative-transfer simulations that include non-local energy deposition and non-thermal processes. We retrieve the basic morphology of PISN light curves from red-supergiant, blue-supergiant, and Wolf-Rayet (WR) star progenitors. Although we confirm that a progenitor 100Msun helium core (PISN model He100) fits well the SN2007bi light curve, the low ratios of its kinetic energy and 56Ni mass to the ejecta mass, similar to standard core-collapse SNe, conspire to produce cool photospheres, red spectra subject to strong line blanketing, and narrow line profiles, all conflicting with SN2007bi observations. He-core models of increasing 56Ni-to-ejecta mass ratio have bluer spectra, but still too red to match SN2007bi, even for model He125 -- the effect of 56Ni heating is offset by the associated increase in blanketing. In contrast, the delayed injection of energy by a magnetar represents a more attractive alternative to reproduce the blue, weakly-blanketed, and broad-lined spectra of super-luminous SNe. The extra heat source is free of blanketing and is not explicitly tied to the ejecta. Experimenting with a ~9Msun WR-star progenitor, initially exploded to yield a ~1.6B SN Ib/c ejecta but later influenced by tunable magnetar-like radiation, we produce a diversity of blue spectral morphologies reminiscent of SN2007bi, the peculiar Type Ib SN2005bf, and super-luminous SN2005ap-like events.

We investigate the formation of CIII 4647-51-50 and CIII 5696 in the atmosphere of O stars to see if they can be reliably used for abundance determinations. We use atmosphere models computed with the code CMFGEN. The key physical ingredients explaining the formation of the CIII lines are extracted from comparisons of models with different stellar parameters and through examining rates controlling the level populations. The strength of CIII 5696 critically depends on UV CIII lines at 386, 574 and 884 A. These lines control the CIII 5696 upper and lower level population. CIII 884 plays a key role in late O stars where it drains the lower level of CIII 5696. CIII 386 and CIII 574 are more important at early spectral types. The overlap of these UV lines with FeIV 386.262, FeIV 574.232 and SV 884.531 influences the radiative transfer at 386, 574 and 884 A, and consequently affects the strength of CIII 5696. CIII 4650 is mainly controlled by the CIII 538 line which acts as a drain on its lower level. FeIV 538.057 interacts with CIII 538 and has an impact on the CIII 4650 profile. Low temperature dielectronic recombinations have a negligible effect on the line profiles. Given our current understanding of the stellar and wind properties of O stars, and in view of the present results, the determination of accurate carbon abundances from CIII 4647-51-50 and CIII 5696 is an extremely challenging task. Uncertainties lower than a factor of two on C/H determinations based only on these two sets of lines should be regarder as highly doubtful. Our results provide a possible explanation of the variability of CIII 4650 in Of?p stars.

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 incorporate non-thermal excitation and ionization processes arising from non-thermal electrons that result from \gamma-ray energy deposition, into our radiative transfer code CMFGEN. The non-thermal electron distribution is obtained by solving the Spencer-Fano equation using the procedure of Kozma & Fransson (1992). We applied the non-thermal calculations to the blue supergiant explosion model whose early evolution was studied in Dessart & Hillier (2010). Non-thermal processes generally increase excitation and ionization and decrease the temperature of the ejecta. We confirm that non-thermal processes are crucial for modeling the nebular spectra. Both optical HI and HeI lines are significantly strengthened. While optical HeI lines are not easily discerned in observational spectra due to severe blending with other lines, HeI 2.058 \mu m provides an excellent opportunity to infer the influence of non-thermal processes. We also discuss the processes controlling the formation of the HeI lines during the nebular epoch. Most lines of other species are only slightly affected. We also show that the inclusion of FeI has substantial line-blanketing effects on the optical spectra. Our model spectra and synthetic light curves are compared to the observations of SN 1987A. The spectral evolution shows broad agreement with the observations, especially H\alpha. The uncertainties of the non-thermal solver are studied, and are expected to be small. With this new addition of non-thermal effects in CMFGEN, we now treat all known important processes controlling the radiative transfer of a supernova ejecta, whatever the type and the epoch.

Using non-LTE time-dependent radiative-transfer calculations, we study the impact of mixing and non-thermal processes associated with radioactive decay on SN IIb/Ib/Ic light curves (LCs) and spectra. Starting with short-period binary models of \leq5Msun He-rich stars (18-25Msun on the main-sequence), we produce 1.2B ejecta which we artificially mix to alter the chemical stratification. While the total 56Ni mass influences the LC peak, the spatial distribution of 56Ni, controlled by mixing processes, impacts both the multi-band LCs and spectra. With enhanced mixing, our synthetic LCs start their post-breakout re-brightening phase earlier, follow a more gradual rise to peak, appear redder, and fade faster after peak due to enhanced gamma-ray escape. Non-thermal electrons, crucial for the production of HeI lines, deposit a dominant fraction of their energy as heat. Because energy deposition is generally local well after the LC peak, the broad HeI lines characteristic of maximum-light SN IIb/Ib spectra require mixing that places 56Ni and helium nuclei to within a gamma-ray mean-free-path. This requirement indicates that SNe IIb and Ib most likely arise from the explosion of stripped-envelope massive stars (main-sequence masses \leq25Msun) that have evolved through mass-transfer in a binary system, rather than from more massive single WR stars. In contrast, the lack of HeI lines in SNe Ic may result from a variety of causes: A genuine helium deficiency; strongly-asymmetric mixing; weak mixing; or a more massive, perhaps single, progenitor characterized by a larger oxygen-rich core. Our models, subject to different mixing magnitudes, can produce a variety of SN types, including IIb, IIc, Ib, and Ic. As it is poorly constrained by explosion models, mixing challenges our ability to infer the progenitor and explosion properties of SNe IIb/Ib/Ic.

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.

We present 2D radiative transfer modeling of the Eta Carinae binary system accounting for the presence of a wind-wind collision (WWC) cavity carved in the optically-thick wind of the primary star. By comparing synthetic line profiles with HST/STIS spectra obtained near apastron, we show that the WWC cavity has a strong influence on multi-wavelength diagnostics. This influence is regulated by the modification of the optical depth in the continuum and spectral lines. We find that H-alpha, H-beta, and Fe II lines are the most affected by the WWC cavity, since they form over a large volume of the primary wind. These spectral lines depend on latitude and azimuth since, according to the orientation of the cavity, different velocity regions of a spectral line are affected. For 2D models with orientation corresponding to orbital inclination angle 110deg < i < 140deg and longitude of periastron 210deg < omega < 330deg, the blueshifted and zero-velocity regions of the line profiles are the most affected. These orbital orientations are required to simultaneously fit the UV and optical spectrum of Eta Car, for a half-opening angle of the cavity in the range 50-70deg. We find that the excess P-Cygni absorption seen in H-alpha, H-beta and optical Fe II lines in spherical models becomes much weaker or absent in the 2D models, in agreement with the observations. The observed UV spectrum of Eta Car, dominated by Fe II absorption lines, is superbly reproduced by our 2D cavity models. Small discrepancies still remain, as H-gamma and H-delta absorptions are overestimated by our models. We suggest that photoionization of the wind of the primary by the hot companion star is responsible for the weak absorption seen in these lines. Our CMFGEN models indicate that the primary star has a mass-loss rate of 8.5x10e-4 Msun/yr and wind terminal velocity of 420 km/s around the 2000 apastron.

In previous papers we discussed results from fully time-dependent radiative transfer models for core-collapse supernova (SN) ejecta, including the Type II-peculiar SN 1987A, the more "generic" SN II-Plateau, and more recently Type IIb/Ib/Ic SNe. Here we describe the modifications to our radiative modeling code, CMFGEN, which allowed those studies to be undertaken. The changes allow for time-dependent radiative transfer of SN ejecta in homologous expansion. In the modeling we treat the entire SN ejecta, from the innermost layer that does not fall back on the compact remnant out to the progenitor surface layers. From our non-LTE time-dependent line-blanketed synthetic spectra, we compute the bolometric and multi-band light curves: light curves and spectra are thus calculated simultaneously using the same physical processes and numerics. These upgrades, in conjunction with our previous modifications which allow the solution of the time dependent rate equations, will improve the modeling of SN spectra and light curves, and hence facilitate new insights into SN ejecta properties, the SN progenitors and the explosion mechanism(s). CMFGEN can now be applied to the modeling of all SN types

We present the first ultraviolet spectrum of the peculiar, magnetic Of?p star HD 108 obtained in its spectroscopic low state. The new data, obtained with the Space Telescope Imaging Spectrograph (STIS) on the Hubble Space Telescope, reveal significant changes compared to IUE spectra obtained in the high state: N V 1240, Si IV 1400, and C IV 1550 present weaker P-Cygni profiles (less absorption) in the new data, while N IV 1718 absorption is deeper, without the clear wind signature evident in the high state. Such changes contrast with those found in other magnetic massive stars, where more absorption is observed in the resonance doublets when the sightline is close to the plane of the magnetic equator. The new data show also that the photospheric Fe IV forest, at 1600--1700 angstroms, has strengthened compared to previous observations. The ultraviolet variability is large compared to that found in typical, non-magnetic O stars, but moderate when compared to the high-/low-state changes reported in the optical spectrum of HD 108 over several decades. We use non-LTE expanding-atmosphere models to analyze the new STIS observations. Overall, the results are in accord with a scenario in which the optical variability is mainly produced by magnetically constrained gas, close to the photosphere. The relatively modest changes found in the main ultraviolet wind lines suggest that the stellar wind is not substantially variable on a global scale. Nonetheless, multidimensional radiative-transfer models may be needed to understand some of the phenomena observed.

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths.

Our goal is to determine the stellar and wind properties of seven O stars in the cluster NGC2244 and three O stars in the OB association MonOB2. These properties give us insight into the mass loss rates of O stars, allow us to check the validity of rotational mixing in massive stars, and to better understand the effects of the ionizing flux and wind mechanical energy release on the surrounding interstellar medium and its influence on triggered star formation. We collect optical and UV spectra of the target stars which are analyzed by means of atmosphere models computed with the code CMFGEN. The spectra of binary stars are disentangled and the components are studied separately. All stars have an evolutionary age less than 5 million years, with the most massive stars being among the youngest. Nitrogen surface abundances show no clear relation with projected rotational velocities. Binaries and single stars show the same range of enrichment. This is attributed to the youth and/or wide separation of the binary systems in which the components have not (yet) experienced strong interaction. A clear trend of larger enrichment in higher luminosity objects is observed, consistent with what evolutionary models with rotation predict for a population of O stars at a given age. We confirm the weakness of winds in late O dwarfs. In general, mass loss rates derived from UV lines are lower than mass loss rates obtained from Ha. The UV mass loss rates are even lower than the single line driving limit in the latest type dwarfs. These issues are discussed in the context of the structure of massive stars winds. The evolutionary and spectroscopic masses are in agreement above 25 Msun but the uncertainties are large. Below this threshold, the few late-type O stars studied here indicate that the mass discrepancy still seems to hold.

Over the past 40 years, the massive LBV/WR system HD 5980 in the Small Magellanic Cloud has undergone a long-term S Doradus type variability cycle and two brief and violent eruptions in 1993 and 1994. In this paper we analyze a collection of UV and optical spectra obtained between 1979 and 2009 and perform CMFGEN model fits to spectra of 1994, 2000, 2002 and 2009. The results are as follows: a) The long term S Dor-type variability is associated with changes of the hydrostatic radius; b) The 1994 eruption involved changes in its bolometric luminosity and wind structure; c) the emission-line strength, the wind velocity and the continuum luminosity underwent correlated variations in the sense that a decreasing V$_\infty$ is associated with increasing emission line and continuum levels; and d) The spectrum of the third star in the system ({\it Star C}) is well-fit by a T$_{eff}$=32 K model atmosphere with SMC chemical abundances. For all epochs, the wind of the erupting star is optically thick at the sonic point and is thus driven mainly by the continuum opacity. We speculate that the wind switches between two stable regimes driven by the "hot" (during the eruption) and the "cool" (post-eruption) iron opacity bumps as defined by Lamers & Nugis (2002) and Gr\"afener and Hamann (2008), and thus the wind may undergo a bi-stability jump of a different nature from that which occurs in OB-stars.

The Type II-Plateau supernova (SN II-P) SN 2004dj was the first SN II-P for which spectropolarimetry data were obtained with fine temporal sampling before, during, and after the fall off of the photometric plateau -- the point that marks the transition from the photospheric to the nebular phase in SNe II-P. Unpolarized during the plateau, SN 2004dj showed a dramatic spike in polarization during the descent off of the plateau, and then exhibited a smooth polarization decline over the next two hundred days. This behavior was interpreted by Leonard et al. (2006) as evidence for a strongly non-spherical explosion mechanism that had imprinted asphericity only in the innermost ejecta. In this brief report, we compare nine similarly well-sampled epochs of spectropolarimetry of the Type II-P SN 2008bk to those of SN 2004dj. In contrast to SN 2004dj, SN 2008bk became polarized well before the end of the plateau and also retained a nearly constant level of polarization through the early nebular phase. Curiously, although the onset and persistence of polarization differ between the two objects, the detailed spectropolarimetric characteristics at the epochs of recorded maximum polarization for the two objects are extremely similar, feature by feature. We briefly interpret the data in light of non-Local-Thermodynamic Equilibrium, time-dependent radiative-transfer simulations specifically crafted for SN II-P ejecta.

We show that the significantly different effective temperatures (Teff) achieved by the luminous blue variable AG Carinae during the consecutive visual minima of 1985-1990 (Teff= 22,800 K) and 2000-2001 (Teff=17,000 K) place the star on different sides of the bistability limit, which occurs in line-driven stellar winds around Teff ~ 21,000 K. Decisive evidence is provided by huge changes in the optical depth of the Lyman continuum in the inner wind as Teff changes during the S Dor cycle. These changes cause different Fe ionization structures in the inner wind. The bistability mechanism is also related to the different wind parameters during visual minima: the wind terminal velocity was 2-3 times higher and the mass-loss rate roughly two times smaller in 1985-1990 than in 2000-2003. We obtain a projected rotational velocity of 220 +- 50 km/s during 1985-1990 which, combined with the high luminosity (L=1.5x10^6 Lsun), puts AG Car extremely close to the Eddington limit modified by rotation (Omega-Gamma limit): for an inclination angle of 90 deg, Gamma_Omega > 1.0 for M< 60 Msun. Based on evolutionary models and mass budget, we obtain an initial mass of ~ 100 Msun and a current mass of 60-70 Msun for AG Car. Therefore, AG Car is close to, if not at, the Omega-Gamma limit during visual minimum. Assuming M=70 Msun, we find that Gamma_Omega decreases from 0.93 to 0.72 as AG Car expands toward visual maximum, suggesting that the star is not above the Eddington limit during maximum phases.

We have developed a spectrum synthesis method for modeling the UV emission from the accretion disk from cataclysmic variables (CVs). The disk is separated into concentric rings, with an internal structure from the Wade & Hubeny disk-atmosphere models. For each ring, a wind atmosphere is calculated in the co-moving frame with a vertical velocity structure obtained from a solution of the Euler equation. Using simple assumptions, regarding rotation and the wind streamlines, these 1D models are combined into a single 2.5D model for which we compute synthetic spectra. We find that the resulting line and continuum behavior as a function of the orbital inclination is consistent with the observations, and verify that the accretion rate affects the wind temperature, leading to corresponding trends in the intensity of UV lines. In general, we also find that the primary mass has a strong effect on the P-Cygni absorption profiles, the synthetic emission line profiles are strongly sensitive to the wind temperature structure, and an increase in the mass loss rate enhances the resonance line intensities. Synthetic spectra were compared with UV data for two high orbital inclination nova-like CVs - RW Tri and V347 Pup. We needed to include disk regions with arbitrary enhanced mass loss to reproduce reasonably well widths and line profiles. This fact and a lack of flux in some high ionization lines may be the signature of the presence of density enhanced regions in the wind, or alternatively, may result from inadequacies in some of our simplifying assumptions.

We present synthetic single-line and continuum linear-polarisation signatures due to electron scattering in axially-symmetric Type II supernovae (SNe) which we calculate using a Monte Carlo and a long-characteristic radiative-transfer code. Aspherical ejecta are produced by prescribing a latitudinal scaling or stretching of SN ejecta inputs obtained from 1-D non-LTE time-dependent calculations. We study polarisation signatures as a function of inclination, shape factor, wavelength, line identity, post-explosion time. At early times, cancellation and optical-depth effects make the polarisation intrinsically low, causing complicated sign reversals with inclination or continuum wavelength, and across line profiles. While the line polarisation is positive (negative) for an oblate (prolate) morphology at the peak and in the red wing, the continuum polarisation may be of any sign. These complex polarisation variations are produced not just by the asymmetric distribution of scatterers but also of the flux. Our early-time signatures are in contradiction with predictions for a centrally illuminated aspherical nebula, although this becomes a better approximation at nebular times. For a fixed asymmetry, our synthetic continuum polarisation is generally low, may evolve non-monotonically during the plateau phase, but it systematically rises as the ejecta become optically thin. Changes in polarization over time do not necessarily imply a change in the asymmetry of the ejecta. The SN structure (e.g., density/ionization) critically influences the level of polarisation. Importantly, a low polarisation (<0.5%) at early times does not necessarily imply a low degree of asymmetry as usually assumed. Asphericity influences line-profile morphology and the luminosity, which may compromise the accuracy of SN characteristics inferred from these.

We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.

We have developed a radiative transfer code, CMFGEN, which allows us to model the spectra of massive stars and supernovae. Using CMFGEN we can derive fundamental parameters such as effective temperatures and surface gravities, derive abundances, and place constraints on stellar wind properties. The last of these is important since all massive stars are losing mass via a stellar wind that is driven from the star by radiation pressure, and this mass loss can substantially influence the spectral appearance and evolution of the star. Recently we have extended CMFGEN to allow us to undertake time-dependent radiative transfer calculations of supernovae. Such calculations will be used to place constraints on the supernova progenitor, to place constraints on the supernova explosion and nucleosynthesis, and to derive distances using a physical approach called the "Expanding Photosphere Method". We describe the assumptions underlying the code and the atomic processes involved. A crucial ingredient in the code is the atomic data. For the modeling we require accurate transition wavelengths, oscillator strengths, photoionization cross-sections, collision strengths, autoionization rates, and charge exchange rates for virtually all species up to, and including, cobalt. Presently, the available atomic data varies substantially in both quantity and quality.

We present the first non-LTE time-dependent radiative-transfer simulations of supernovae (SNe) II-Plateau (II-P) covering both the photospheric and nebular phases, from ~10 to >~1000d after the explosion, and based on 1.2B piston-driven ejecta produced from a 15Msun and a 25Msun non-rotating solar-metallicity star. The radial expansion of the gradually cooling photosphere gives rise to a near-constant luminosity up to >~100d after explosion. The photosphere remains in the outer 0.5Msun of the ejecta for up to ~50d after explosion. As the photosphere reaches the edge of the helium core, the SN luminosity drops by an amount mitigated by the progenitor radius and the 56Ni mass. Synthetic light-curves exhibit a bell-shape morphology, evolving faster for more compact progenitors, and with an earlier peak and narrower width in bluer filters. UV and U-band fluxes are very sensitive to line-blanketing, the metallicity, and the adopted model atoms. During the recombination epoch synthetic spectra are dominated by HI and metal lines, and are largely insensitive to the differing H/He/C/N/O composition of our two progenitor stars. In contrast, synthetic nebular-phase spectra reveal a broader/stronger OI doublet line in the higher-mass progenitor model, reflecting the larger masses of oxygen and nickel that are ejected. Our simulations overestimate the typical luminosity and the visual rise time of standard SNe II-P, likely a consequence of our progenitor stars being too big and/or too hydrogen rich. Comparison of our simulations with photospheric-phase observations of SN1999em of the same color are satisfactory. Our neglect of non-thermal processes leads to a fast disappearance of continuum radiation and Balmer-line emission at the end of the plateau phase. With the exception of HI lines, our nebular spectra show a striking similarity to contemporaneous observations of SN1999em.

We present a method for computing the net transmission of X-rays emitted by shock-heated plasma distributed throughout a partially optically thick stellar wind from a massive star. We find the transmission by an exact integration of the formal solution, assuming that the emitting plasma and absorbing plasma are mixed at a constant mass ratio above some minimum radius, below which there is assumed to be no emission. This model is more realistic than either the slab absorption associated with a corona at the base of the wind or the exospheric approximation that assumes that all observed X-rays are emitted without attenuation from above the radius of optical depth unity. Our model is implemented in XSPEC as a pre-calculated table that can be coupled to a user-defined table of the wavelength dependent wind opacity. We provide a default wind opacity model that is more representative of real wind opacities than the commonly used neutral interstellar medium (ISM) tabulation. Preliminary modeling of \textit{Chandra} grating data indicates that the X-ray hardness trend of OB stars with spectral subtype can largely be understood as a wind absorption effect.

We analyze interferometric measurements of the Luminous Blue Variable Eta Carinae with the goal of constraining the rotational velocity of the primary star and probing the influence of the companion. Using 2-D radiative transfer models of latitude-dependent stellar winds, we find that prolate wind models with a ratio of the rotational velocity (vrot) to the critical velocity (vcrit) of W=0.77-0.92, inclination angle of i=60-90 degrees, and position angle PA=108-142 degrees reproduce simultaneously K-band continuum visibilities from VLTI/VINCI and closure phase measurements from VLTI/AMBER. Interestingly, oblate models with W=0.73-0.90 and i=80-90 degrees produce similar fits to the interferometric data, but require PA=210-230 degrees. Therefore, both prolate and oblate models suggest that the rotation axis of the primary star is not aligned with the Homunculus polar axis. We also compute radiative transfer models of the primary star allowing for the presence of a cavity and dense wind-wind interaction region created by the companion star. We find that the wind-wind interaction has a significant effect on the K-band image mainly via free-free emission from the compressed walls and, for reasonable model parameters, can reproduce the VLTI/VINCI visibilities taken at phase 0.92-0.93. We conclude that the density structure of the primary wind can be sufficiently disturbed by the companion, thus mimicking the effects of fast rotation in the interferometric observables. Therefore, fast rotation may not be the only explanation for the interferometric observations. Intense temporal monitoring and 3-D modeling are needed to resolve these issues.

We report near-IR spectroscopic observations of the Eta Carinae massive binary system during 2008-2009 using VLT/CRIRES. We detect a strong, broad absorption wing in He I 10833 extending up to -1900 km/s across the 2009.0 spectroscopic event. Archival HST/STIS ultraviolet and optical data shows a similar high-velocity absorption (up to -2100 km/s) in the UV resonance lines of Si IV 1394, 1403 across the 2003.5 event. UV lines from low-ionization species, such as Si II 1527, 1533 and C II 1334, 1335, show absorption up to -1200 km/s, indicating that the absorption with v from -1200 to -2100 km/s originates in a region markedly faster and more ionized than the nominal wind of the primary star. Observations obtained at the OPD/LNA during the last 4 spectroscopic cycles (1989-2009) also display high-velocity absorption in He I 10833 during periastron. Based on the OPD/LNA dataset, we determine that material with v < -900 km/s is present in the phase range 0.976 < phi < 1.023 of the spectroscopic cycle, but absent in spectra taken at phi < 0.947 and phi > 1.049. Therefore, we constrain the duration of the high-velocity absorption to be 95 to 206 days (or 0.047 to 0.102 in phase). We suggest that the high-velocity absorption originates from shocked gas in the wind-wind collision zone, at distances of 15 to 45 AU in the line-of-sight to the primary star. Using 3-D hydrodynamical simulations of the wind-wind collision zone, we find that the dense high-velocity gas is in the line-of-sight to the primary star only if the binary system is oriented in the sky so that the companion is behind the primary star during periastron, corresponding to a longitude of periastron of omega ~ 240 to 270 degrees. We study a possible tilt of the orbital plane relative to the Homunculus equatorial plane and conclude that our data are broadly consistent with orbital inclinations in the range i=40 to 60 degrees.

We discuss a new one-dimensional non-LTE time-dependent radiative-transfer technique for the simulation of supernova (SN) spectra and light curves. Starting from a hydrodynamical input characterizing the homologously-expanding ejecta at a chosen post-explosion time, we model the evolution of the entire ejecta, including gas and radiation. Non-LTE, which holds in all regions at and above the photosphere, is accounted for. The effects of line blanketing on the radiation field are explicitly included, using complex model atoms and solving for all ion level populations appearing in the statistical-equilibrium equations. Here, we present results for SN1987A, evolving the model "lm18a7Ad" of Woosley from 0.27 to 20.8d. The fastest evolution occurs prior to day 1, with a spectral energy distribution peaking in the range 300-2000A, subject to line blanketing from highly ionized metal and CNO species. After day 1, our synthetic multi-band light curve and spectra reproduce the observations to within 10-20% in flux in the optical, with a greater mismatch for the faint UV flux. We do not encounter any of the former discrepancies associated with the HeI and HI lines in the optical, which can be fitted well with a standard Blue-supergiant-star surface composition and no contribution from radioactive decay. The effects of time dependence on the ionization structure, discussed in Dessart & Hillier, are recovered, and thus nicely integrated in this new scheme. Despite the 1D nature of our approach, its high physical consistency and accuracy will allow reliable inferences to be made on explosion properties and pre-SN star evolution.

We fit every emission line in the high-resolution Chandra grating spectrum of zeta Pup with an empirical line profile model that accounts for the effects of Doppler broadening and attenuation by the bulk wind. For each of sixteen lines or line complexes that can be reliably measured, we determine a best-fitting fiducial optical depth, tau_* = kappa*Mdot/4{pi}R_{\ast}v_{\infty}, and place confidence limits on this parameter. These sixteen lines include seven that have not previously been reported on in the literature. The extended wavelength range of these lines allows us to infer, for the first time, a clear increase in tau_* with line wavelength, as expected from the wavelength increase of bound-free absorption opacity. The small overall values of tau_*, reflected in the rather modest asymmetry in the line profiles, can moreover all be fit simultaneously by simply assuming a moderate mass-loss rate of 3.5 \pm 0.3 \times 10^{-6} Msun/yr, without any need to invoke porosity effects in the wind. The quoted uncertainty is statistical, but the largest source of uncertainty in the derived mass-loss rate is due to the uncertainty in the elemental abundances of zeta Pup, which affects the continuum opacity of the wind, and which we estimate to be a factor of two. Even so, the mass-loss rate we find is significantly below the most recent smooth-wind H-alpha mass-loss rate determinations for zeta Pup, but is in line with newer determinations that account for small-scale wind clumping. If zeta Pup is representative of other massive stars, these results will have important implications for stellar and galactic evolution.

We report optical observations of the Luminous Blue Variable (LBV) HR Carinae which show that the star has reached a visual minimum phase in 2009. More importantly, we detected absorptions due to Si IV 4088-4116 Angstroms. To match their observed line profiles from 2009 May, a high rotational velocity of vrot=150 +- 20 km/s is needed (assuming an inclination angle of 30 degrees), implying that HR Car rotates at ~0.88 +- 0.2 of its critical velocity for break-up (vcrit). Our results suggest that fast rotation is typical in all strong-variable, bona-fide galactic LBVs, which present S Dor-type variability. Strong-variable LBVs are located in a well-defined region of the HR diagram during visual minimum (the "LBV minimum instability strip"). We suggest this region corresponds to where vcrit is reached. To the left of this strip, a forbidden zone with vrot/vcrit>1 is present, explaining why no LBVs are detected in this zone. Since dormant/ex LBVs like P Cygni and HD 168625 have low vrot, we propose that LBVs can be separated in two groups: fast-rotating, strong-variable stars showing S-Dor cycles (such as AG Car and HR Car) and slow-rotating stars with much less variability (such as P Cygni and HD 168625). We speculate that SN progenitors which had S-Dor cycles before exploding (such as in SN 2001ig, SN 2003bg, and SN 2005gj) could have been fast rotators. We suggest that the potential difficulty of fast-rotating Galactic LBVs to lose angular momentum is an additional evidence that such stars could explode during the LBV phase.

We present a detailed spectroscopic analysis of the luminous blue variable AG Carinae during the last two visual minimum phases of its S-Dor cycle (1985-1990 and 2000-2003). The analysis reveals an overabundance of He, N, and Na, and a depletion of H, C, and O, on the surface of AG Car, indicating the presence of CNO-processed material. Furthermore, the ratio N/O is higher on the stellar surface than in the nebula. We found that the minimum phases of AG Car are not equal to each other, since we derived a noticeable difference between the maximum effective temperature achieved during 1985-1990 (22,800 K) and 2000-2001 (17,000 K). While the wind terminal velocity was 300 km/s in 1985-1990, it was as low as 105 km/s in 2001. The mass-loss rate, however, was lower from 1985-1990 (1.5 x 10^(-5) Msun/yr) than from 2000-2001 (3.7 x 10^(-5) Msun/yr). We found that the wind of AG Car is significantly clumped (f=0.10 - 0.25) and that clumps must be formed deep in the wind. We derived a bolometric luminosity of 1.5 x 10^6 Lsun during both minimum phases which, contrary to the common assumption, decreases to 1.0 x 10^6 Lsun as the star moves towards maximum flux in the V band. Assuming that the decrease in the bolometric luminosity of AG Car is due to the energy used to expand the outer layers of the star (Lamers 1995), we found that the expanding layers contain roughly 0.6 - 2 Msun. Such an amount of mass is an order of magnitude lower than the nebular mass around AG Car, but is comparable to the nebular mass found around lower-luminosity LBVs and to that of the Little Homunculus of Eta Car. If such a large amount of mass is indeed involved in the S Dor-type variability, we speculate that such instability could be a failed Giant Eruption, with several solar masses never becoming unbound from the star.(abridged)

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).

We derive the physical properties of three WNh stars in the SMC to constrain stellar evolution beyond the main sequence at low metallicity and to investigate the metallicity dependence of the clumping properties of massive stars. We compute atmosphere models to derive the stellar and wind properties of the three WNh targets. A FUV/UV/optical/near-infrared analysis gives access to temperatures, luminosities, mass loss rates, terminal velocities and stellar abundances. All stars still have a large hydrogen mass fraction in their atmosphere, and show clear signs of CNO processing in their surface abundances. One of the targets can be accounted for by normal stellar evolution. It is a star with initial mass around 40-50 Msun in, or close to, the core He burning phase. The other two objects must follow a peculiar evolution, governed by fast rotation. In particular, one object is likely evolving homogeneously due to its position blue-ward of the main sequence and its high H mass fraction. The clumping factor of one star is found to be 0.15+/-0.05. This is comparable to values found for Galactic Wolf-Rayet stars, indicating that within the uncertainties, the clumping factor does not seem to depend on metallicity.

We present a multi-epoch quantitative spectroscopic analysis of the Type IIn SN 1994W, an event interpreted by Chugai et al. as stemming from the interaction between the ejecta of a SN and a 0.4Msun circumstellar shell ejected 1.5yr before core collapse. During the brightening phase, our models suggest that the source of optical radiation is not unique, perhaps associated with an inner optically-thick Cold Dense Shell (CDS) and outer optically-thin shocked material. During the fading phase, our models support a single source of radiation, an hydrogen-rich optically-thick layer with a near-constant temperature of ~7000K that recedes from a radius of 4.3x10^15 at peak to 2.3x10^15cm 40 days later. We reproduce the hybrid narrow-core broad-wing line profile shapes of SN 1994W at all times, invoking an optically-thick photosphere exclusively (i.e., without any external optically-thick shell). In SN 1994W, slow expansion makes scattering with thermal electrons a key escape mechanism for photons trapped in optically-thick line cores, and allows the resulting broad incoherent electron-scattering wings to be seen around narrow line cores. In SNe with larger expansion velocities, the thermal broadening due to incoherent scattering is masked by the broad profile and the dominant frequency redshift occasioned by bulk motions. Given the absence of broad lines at all times and the very low 56Ni yields, we speculate whether SN 1994W could have resulted from an interaction between two ejected shells without core collapse. The high conversion efficiency of kinetic to thermal energy may not require a SN-like energy budget for SN1994W.

We present a measurement of metallicity in the Galactic center Quintuplet Cluster made using quantitative spectral analysis of two Luminous Blue Variables (LBVs). The analysis employs line-blanketed NLTE wind/atmosphere models fit to high-resolution near-infrared spectra containing lines of H, HeI, SiII, MgII, and FeII. We are able to break the H/He ratio vs. mass-loss rate degeneracy found in other LBVs and to obtain robust estimates of the He content of both objects. Our results indicate solar iron abundance and roughly twice solar abundance in the alpha-elements. These results are discussed within the framework of recent measurements of oxygen and carbon composition in the nearby Arches Cluster and iron abundances in red giants and supergiants within the central 30 pc of the Galaxy. The relatively large enrichment of alpha-elements with respect to iron is consistent with a history of more nucleosynthesis in high mass stars than the Galactic disk.

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 the serendipitous detection by GALEX of fast (<1 day) rising (>1 mag) UV emission from two Type II plateau (II-P) supernovae (SNe) at z=0.185 and 0.324 discovered by the Supernova Legacy Survey. Optical photometry and VLT spectroscopy 2 weeks after the GALEX detections link the onset of UV emission to the time of shock breakout. Using radiation hydrodynamics and non-LTE radiative transfer simulations, and starting from a standard red supergiant (RSG; Type II-P SN progenitor) star evolved self-consistently from the main sequence to iron core collapse, we model the shock breakout phase and the 55 hr that follow. The small scale height of our RSG atmosphere model suggests that the breakout signature is a thermal soft X-ray burst (lambda_peak ~ 90\AA) with a duration of <~ 2000 s. Longer durations are possible but require either an extended and tenuous non-standard envelope, or an unusually dense RSG wind with \dot{M} ~ 10^(-3) Msun yr^(-1). The GALEX observations miss the peak of the luminous (M_FUV ~ -20) UV burst but unambiguously capture the rise of the emission and a subsequent 2 day long plateau. The postbreakout, UV-bright plateau is a prediction of our model in which the shift of the peak of the spectral energy distribution (SED) from ~100 to ~1000\AA and the ejecta expansion both counteract the decrease in bolometric luminosity from ~10^11 to ~10^9 L_sun over that period. Based on the observed detection efficiency of our study we make predictions for the breakout detection rate of the GALEX Time Domain Survey.

A full description of the 5.5-yr low excitation events in Eta Carinae is presented. We show that they are not as simple and brief as previously thought, but a combination of two components. The first, the 'slow variation' component, is revealed by slow changes in the ionization level of circumstellar matter across the whole cycle and is caused by gradual changes in the wind-wind collision shock-cone orientation, angular opening and gaseous content. The second, the 'collapse' component, is restricted to around the minimum, and is due to a temporary global collapse of the wind-wind collision shock. High energy photons (E > 16 eV) from the companion star are strongly shielded, leaving the Weigelt objects at low ionization state for >6 months. High energy phenomena are sensitive only to the 'collapse', low energy only to the 'slow variation' and intermediate energies to both components. Simple eclipses and mechanisms effective only near periastron (e.g., shell ejection or accretion onto the secondary star) cannot account for the whole 5.5-yr cycle. We find anti-correlated changes in the intensity and the radial velocity of P Cygni absorption profiles in FeII 6455 and HeI 7065 lines, indicating that the former is associated to the primary and the latter to the secondary star. We present a set of light curves representative of the whole spectrum, useful for monitoring the next event (2009 January 11).

We analyze the chemical composition of the central star of the planetary nebula NGC 6543 based upon a detailed NLTE model of its stellar wind. The logarithmic abundances by number are H=12.00, He=11.00, C=9.03, N=8.36, O=9.02, Si=8.19, P=5.53, S=7.57 and Fe=7.24. Compared with the solar abundances, most of the elements have solar composition with respect to hydrogen except C which is overabundant by 0.28 dex and Fe which is depleted by $\sim 0.2$ dex. Contrary to most previous work, we find that the star is not H-poor and has a normal He composition. These abundances are compared with those found in the diffuse X-ray plasma and the nebular gas. Compared to the plasma emitting in diffuse X-rays, the stellar wind is much less depleted in iron. Since the iron depletions in the nebular gas and X-ray plasma are similar, we conclude that the plasma emitting diffuse X-rays is derived from the nebular gas rather than the stellar wind. Excellent agreement is obtained between the abundances in the stellar wind and the nebular recombination line abundances for He, C, and O relative to H. On the other hand, the derived stellar N abundance is smaller than the nebular N abundance derived from recombination lines and agrees with the abundance found from collisionally-excited lines. The mean temperature variation determined by five different methods indicates that the difference in the nebular abundances between the recombination lines and collisionally excited lines can be explained as due to the temperature variations in a chemically homogeneous medium.

We calculate the hydrogen and helium ionization in B[e] envelopes and explore their dependence on mass-loss and effective temperature. We also present simulated observations of the Halpha emission line and the C IV 1550 doublet, and study their behavior. This paper reports our first results in an ongoing study of B[e] supergiants, and provides a glimpse on the ionization of the most important elements in self-consistent numerical simulations. Our newly developed 2D stellar atmosphere code, ASTAROTH, was used for the numerical simulations. The code self-consistently solves for the continuum radiation, non-LTE level populations, and electron temperature in axi-symmetric stellar envelopes. Observed profiles were calculated by an auxiliary program developed separately from ASTAROTH. In all but one of our models, H remained fully ionized. Due to ionizations from excited states it is much more difficult to get a H neutral disk than indicated by previous analytical calculations. Near the poles, the ionization is high in all models, while helium recombined in the equatorial regions for all but our lowest mass-loss rate. Although the model parameters were not adjusted to provide fits to any particular star, the theoretical profiles show some features seen in the profiles of R126. These include the partially resolved double peaked profile of Halpha, and the weak emission associated with the UV C IV resonance line.

Extensive spectral observations of eta Carinae over the last cycle, and particularly around the 2003.5 low excitation event, have been obtained. The variability of both narrow and broad lines, when combined with data taken from two earlier cycles, reveal a common and well defined period. We have combined the cycle lengths derived from the many lines in the optical spectrum with those from broad-band X-rays, optical and near-infrared observations, and obtained a period length of 2022.7+-1.3 d. Spectroscopic data collected during the last 60 years yield an average period of 2020+-4 d, consistent with the present day period. The period cannot have changed by more than $\Delta$P/P=0.0007 since 1948. This confirms the previous claims of a true, stable periodicity, and gives strong support to the binary scenario. We have used the disappearance of the narrow component of HeI 6678 to define the epoch of the Cycle 11 minimum, T_0=JD 2,452,819.8. The next event is predicted to occur on 2009 January 11 (+-2 days). The dates for the start of the minimum in other spectral features and broad-bands is very close to this date, and have well determined time delays from the HeI epoch.

We analyze the Type II Plateau supernovae (SN II-P) 2005cs and 2006bp with the non-LTE model atmosphere code CMFGEN. We fit 13 spectra in the first month for SN 2005cs and 18 for SN 2006bp. {\sl Swift} ultraviolet photometry and ground-based optical photometry calibrate each spectrum. Our analysis shows both objects were discovered less than 3 days after they exploded, making these the earliest SN II-P spectra ever studied. They reveal broad and very weak lines from highly-ionized fast ejecta with an extremely steep density profile. We identify He{\sc ii} 4686\AA emission in the SN 2006bp ejecta. Days later, the spectra resemble the prototypical Type II-P SN 1999em, which had a supergiant-like photospheric composition. Despite the association of SN 2005cs with possible X-ray emission, the emergent UV and optical light comes from the photosphere, not from circumstellar emission. We surmise that the very steep density fall-off we infer at early times may be a fossil of the combined actions of the shock wave passage and radiation driving at shock breakout. Based on tailored CMFGEN models, the direct-fitting technique and the Expanding Photosphere Method both yield distances and explosion times that agree within a few percent. We derive a distance to NGC 5194, the host of SN 2005cs, of 8.9$\pm$0.5 Mpc and 17.5$\pm$0.8 Mpc for SN 2006bp in NGC 3953. The luminosity of SN 2006bp is 1.5 times that of SN 1999em, and 6 times that of SN 2005cs. Reliable distances to Type II-P supernovae that do not depend on a small range in luminosity provide an independent route to the Hubble Constant and improved constraints on other cosmological parameters.

We study a sample composed of 28 of the brightest stars in the Arches cluster. We analyze K-band spectra obtained with the integral field spectrograph SINFONI on the VLT. Atmosphere models computed with the code CMFGEN are used to derive the effective temperatures, luminosities, stellar abundances, mass loss rates and wind terminal velocities. We find that the stars in our sample are either H-rich WN7-9 stars (WN7-9h) or O supergiants, two being classified as OIf+. All stars are 2-4 Myr old. There is marginal evidence for a younger age among the most massive stars. The WN7-9h stars reach luminosities as large as 2 x 1e6 Lsun, consistent with initial masses of ~ 120 Msun. They are still quite H-rich, but show both N enhancement and C depletion. They are thus identified as core H-burning objects showing products of the CNO equilibrium at their surface. Their progenitors are most likely supergiants of spectral types earlier than O4-6 and initial masses > 60 Msun. Their winds follow a well defined modified wind momentum - luminosity relation (WLR): this is a strong indication that they are radiatively driven. Stellar abundances tend to favor a slightly super solar metallicity, at least for the lightest metals. We note however that the evolutionary models seem to under-predict the degree of N enrichment.

Spectroscopic modeling of Type II supernovae (SNe) generally assumes steady-state. Following the recent suggestion of Utrobin & Chugai, but using the 1D non-LTE line-blanketed model atmosphere code CMFGEN, we investigate the effects of including time-dependent terms that appear in the statistical and radiative equilibrium equations. We base our discussion on the ejecta properties and the spectroscopic signatures obtained from time-dependent simulations, investigating different ejecta configurations, and covering their evolution from one day to six weeks after shock breakout. Compared to equivalent steady-state models, our time-dependent models produce SN ejecta that are systematically over-ionized, affecting helium at one week after explosion, but ultimately affecting all ions after a few weeks. While the continuum remains essentially unchanged, time-dependence effects on observed spectral lines are large. At the recombination epoch, HI lines and NaID are considerably stronger and broader than in equivalent steady-state models, while CaII8500A is weakened. If time dependence is allowed for, the HeI lines at 5875A and 10830A appear about 3 times stronger at one week, and HeI10830A persists as a blue-shifted absorption feature even at 6 weeks after explosion. Time dependence operates through the energy gain from changes in ionization and excitation, and, perhaps more universally across SN types, from the competition between recombination and expansion, which in-turn, can be affected by optical-depth effects. Our time-dependent models compare well with observations of the low-luminosity low-velocity SN 1999br and the more standard SN 1999em, reproducing the Halpha line strength at the recombination epoch, and without the need for setting unphysical requirements on the magnitude of nickel mixing.

We study the stellar and wind properties of massive stars in the central cluster of the Galaxy. We use non-LTE atmosphere models including winds and line-blanketing to fit their H and K band spectra obtained with the 3D spectrograph SINFONI on the VLT. We derive the main stellar (Teff, L, abundances, ionizing flux) and wind (mass loss rate, terminal velocity) properties. They are found to be similar to other galactic massive stars. We show that a direct evolutionary link between Ofpe/WN9, WN8 and WN/C stars exists. Using individual SEDs for each massive star, we construct the total spectral energy distribution of the cluster and use it to compute photoionization models. We show that the nebular properties of the central HII region are well reproduced. We conclude that, contrary to previous claims, standard stellar evolution and atmosphere models are well suited to explain the properties of the central cluster. Our results indicate that massive stars in the central cluster do not have a peculiar evolution as could be expected from their proximity to the supermassive black hole SgrA*.

We present the results of a detailed spectroscopic analysis of 20 years of observations of AG Carinae using the radiative transfer code CMFGEN. Among the conclusions of this work, we highlight the importance of including time-dependent effects in the analysis of the full S Dor cycle. We obtained that the mass-loss rate is approximately constant during the cool phases, implying that the S Dor-type eruptions begin well earlier than the maximum seen in the visual lightcurve. We also determined that the S Dor cycles are ultimately a consequence of an increase/decrease of the hydrostatic radius in combination with the formation of a pseudo-photosphere.

We report the results of an X-ray observing campaign on the massive, evolved star Eta Carinae, concentrating on the 2003 X-ray minimum as seen by the XMM-Newton observatory. These are the first spatially-resolved X-ray monitoring observations of the stellar X-ray spectrum during the minimum. The hard X-ray emission, believed to be associated with the collision of Eta Carinae's wind with the wind from a massive companion star, varied strongly in flux on timescales of days, but not significantly on timescales of hours. The lowest X-ray flux in the 2-10 keV band seen by XMM-Newton was only 0.7% of the maximum seen by RXTE just before the X-ray minimum. The slope of the X-ray continuum above 5 keV did not vary in any observation, which suggests that the electron temperature of the hottest plasma associated with the stellar source did not vary significantly at any phase. Through the minimum, the absorption to the stellar source increased by a factor of 5-10 to NH ~3-4E23 cm-2. The thermal Fe XXV emission line showed significant excesses on both the red and blue sides of the line outside the minimum and exhibited an extreme red excess during the minimum. The Fe fluorescence line at 6.4 keV increased in equivalent width from 100 eV outside the minimum to 200 eV during the minimum. From these observed features, we discuss two possible causes of the X-ray minimum; the eclipse of the X-ray plasma and an intrinsic fading of the X-ray emissivity. The drop in the colliding wind X-ray emission also revealed the presence of an additional X-ray component which exhibited no variation on timescales of weeks to years. This component may be produced by the collision of high speed outflows at v \~1000-2000 km s-1 from Eta Carinae with ambient gas within a few thousand AU from the star.

We have analyzed high spatial, moderate spectral resolution observations of Eta Carinae obtained with the STIS from 1998.0 to 2004.3. The spectra show prominent P-Cygni lines in H I, Fe II and He I which are complicated by blends and contamination by nebular emission and absorption along the line-of-sight toward the observer. All lines show phase and species dependent variations in emission and absorption. For most of the cycle the He I emission is blueshifted relative to the H I and Fe II P-Cygni emission lines, which are approximately centered at system velocity. The blueshifted He I absorption varies in intensity and velocity throughout the 2024 day period. We construct radial velocity curves for the absorption component of the He I and H I lines. The He I absorption shows significant radial velocity variations throughout the cycle, with a rapid change of over 200 km/s near the 2003.5 event. The H I velocity curve is similar to that of the He I absorption, though offset in phase and reduced in amplitude. We interpret the complex line profile variations in He I, H I and Fe II to be a consequence of the dynamic interaction of the dense wind of Eta Car A with the less dense, faster wind plus the radiation field of a hot companion star, Eta Car B. During most of the orbit, Eta Car B and the He+ recombination zone are on the near side of Eta Car A, producing blueshifted He I emission. He I absorption is formed in the part of the He+ zone that intersects the line-of-sight toward Eta Car. We use the variations seen in He I and the other P-Cygni lines to constrain the geometry of the orbit and the character of Eta Car B.

We present far-UV to optical analyses of four hydrogen deficient central stars of planetary nebulae: BD+303639, NGC 40, NGC 5315 and NGC 6905. Using the radiative transfer code CMFGEN, we determined new physical parameters and chemical abundances for these stars. The results were analyzed in the context of the [WR] => PG 1159 evolution via the transformed radius(Rt)-temperature and HR diagrams. NGC 5315 showed itself as an odd object among the previously analyzed central stars. Its temperature (~76kK) is considerably lower than other early-type [WR] stars (~120-150kK). From our models for NGC 5315 and NGC 6905, it is unclear if early-type [WR] stars have smaller C/He mass ratios than other spectral classes, as claimed in the literature. A ratio of ~0.8 is found for NGC 6905. We analyzed FUSE spectra of these stars for the first time, and identified phosphorus in the spectra of BD+303639, NGC 40 and NGC 5315 through the transitions P V 1118,1128. The Fe, Si, P, S and Ne abundances were analyzed in the context of the nucleosynthesis occurring in previous evolutionary phases. We found evidence for Fe deficiency in BD +30 3639 and NGC 5315, and from fits to the Si IV lines we determined a solar Si abundance for BD+303639 and NGC 40. For phosphorus, an oversolar abundance in the NGC 5315 model was preferred, while in the other stars a solar phosphorus abundance cannot be discarded. Regarding sulfur, we estimated upper limits for its abundance, since no conspicuous lines can be seen in the observed spectra. We found that Ne is overabundant in BD +30 3639. In the other stars, Ne is weak or undetectable and upper limits for its abundance were estimated. Our results are in agreement with theoretical predictions and show the usefulness of [WR] stars as testbeds for nucleosynthesis calculations in the AGB and post-AGB phases.(abridged)

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 have made an extensive study of the UV spectrum of Eta Carinae, and find that we do not directly observe the star and its wind in the UV. Because of dust along our line of sight, the UV light that we observe arises from bound-bound scattering at large impact parameters. We obtain a reasonable fit to the UV spectrum by using only the flux that originates outside 0.033". This explains why we can still observe the primary star in the UV despite the large optical extinction -- it is due to the presence of an intrinsic coronagraph in the Eta Carinae system, and to the extension of the UV emitting region. It is not due to peculiar dust properties alone. We have computed the spectrum of the purported companion star, and show that it could only be directly detected in the UV spectrum preferentially in the Far Ultraviolet Spectroscopic Explorer (FUSE) spectral region (912-1175 Ang.). However, we find no direct evidence for a companion star, with the properties indicated by X-ray studies and studies of the Weigelt blobs, in UV spectra. This might be due to reprocessing of the companion's light by the dense stellar wind of the primary. Broad FeII and [FeII] emission lines, which form in the stellar wind, are detected in spectra taken in the SE lobe, 0.2" from the central star. The wind spectrum shows some similarities to the spectra of the B & D Weigelt blobs, but also shows some marked differences in that high excitation lines, and lines pumped by Ly-alpha, are not seen. The detection of the broad lines lends support to our interpretation of the UV spectrum, and to our model for Eta Carinae.

We report the detection of broad absorptions due to Si IV 4088-4116 A in the Luminous Blue Variable (LBV) AG Carinae during its last hot phase (2001-2003). Our NLTE spectral analysis, with the radiative transfer code CMFGEN, revealed the photospheric nature of these lines predicting, however, much narrower and deeper absorption profiles than observed. Using a recently-developed code to compute synthetic spectra in 2D geometry allowing for the effects of rotation, we could match the broad absorptions with a high projected rotational velocity of 190 +/- 30 km/s on 2001 April. Analysis of spectra obtained on 2002 March and 2003 January, when the star was cooling, yielded to a projected rotational velocity of 110 +/- 10 km/s and 85 +/- 10 km/s, respectively. The derived rotational velocities are proportional to R^-1, as expected from angular momentum conservation. We discuss the effects of such high rotation on the spectral analysis of AG Car, and on the wind terminal velocity. Our results show direct spectroscopic evidence, for the first time, that a LBV may rotate at a significant fraction of its break-up velocity. Thus, AG Car (and possibly other LBVs) is indeed close to the Gamma-Omega limit, as predicted by theoretical studies of LBVs.

We describe and analyze Hubble Space Telescope (HST) observations of transient emission near 4680 {\AA} in Eta Car, reported earlier by Steiner & Damineli (2004). If, as seems probable, this is He II $\lambda$4687, then it is a unique clue to Eta Car's 5.5-year cycle. According to our analysis, several aspects of this feature support a mass-ejection model of the observed spectroscopic events, and not an eclipse model. The He II emission appeared in early 2003, grew to a brief maximum during the 2003.5 spectroscopic event, and then abruptly disappeared. It did not appear in any other HST spectra before or after the event. The peak brightness was larger than previously reported, and is difficult to explain even if one allows for an uncertainty factor of order 3. The stellar wind must provide a temporary larger-than-normal energy supply, and we describe a special form of radiative amplification that may also be needed. These characteristics are consistent with a class of mass-ejection or wind-disturbance scenarios, which have implications for the physical structure and stability of Eta Car.

We present a new radiative transfer code for axi-symmetric stellar atmospheres and compare test results against 1D and 2D models with and without velocity fields. The code uses the short characteristic method with modifications to handle axi-symmetric and non-monotonic 3D wind velocities, and allows for distributed calculations. The formal solution along a characteristic is evaluated with a resolution that is proportional to the velocity gradient along the characteristic. This allows us to accurately map the variation of the opacities and emissivities as a function of frequency and spatial coordinates, but avoids unnecessary work in low velocity regions. We represent a characteristic with an impact-parameter vector p (a vector that is normal to the plane containing the characteristic and the origin) rather than the traditional unit vector in the direction of the ray. The code calculates the incoming intensities for the characteristics by a single latitudinal interpolation without any further interpolation in the radiation angles. Using this representation also provides a venue for distributed calculations since the radiative transfer can be done independently for each p.

We report the detection of a hot companion of $\eta$ Carinae using high resolution spectra (905 - 1180 \AA) obtained with the Far Ultraviolet Spectroscopic Explorer (\fuse) satellite. Observations were obtained at two epochs of the 2024-day orbit: 2003 June during ingress to the 2003.5 X-ray eclipse and 2004 April several months after egress. These data show that essentially all the far-UV flux from \etacar shortward of \lya disappeared at least two days before the start of the X-ray eclipse (2003 June 29), implying that the hot companion, \etaB, was also eclipsed by the dense wind or extended atmosphere of \etaA. Analysis of the far-UV spectrum shows that \etaB is a luminous hot star. The \nii \wll1084-1086 emission feature suggests that it may be nitrogen-rich. The observed far-UV flux levels and spectral features, combined with the timing of their disappearance, is consistent with \etacar\ being a massive binary system.

This work presents a detailed quantitative spectroscopic analysis of, and the determination of the distance to, the type II supernovae (SN) SN1999em with CMFGEN (Dessart & Hillier 2005a), based on spectrophotometric observations at eight dates up to 40 days after discovery. We use the same iron-group metal content for the ejecta, the same power-law density distribution (with exponent n~10), and a Hubble-velocity law at all times. We adopt a H/He/C/N/O abundance pattern compatible with CNO-cycle equilibrium values for a RSG/BSG progenitor, with C/O enhanced and N depleted at later times. Based on our synthetic fits to spectrophotometric observations of SN1999em, we obtain a distance of 11.5Mpc, similar to that of Baron et al. (2004) and the Cepheid distance to the galaxy host of 11.7Mpc (Leonard et al. 2003). Similarly, based on such models, the Expanding Photosphere Method (EPM) delivers a distance of 11.6Mpc, with negligible scatter between photometric bandpass sets; there is thus nothing wrong with the EPM as such. Previous determinations using the tabulated correction factors of Eastman et al. (1996) all led to 30-50% underestimates: we find that this is caused by 1) an underestimate of the correction factors compared to the only other study of the kind by Dessart & Hillier (2005b), 2) a neglect of the intrinsic >20% scatter of correction factors, and 3) the use of the EPM at late times when severe line blanketing makes the method inaccurate. The need of detailed model computations for reliable EPM distance estimates thus defeats the appeal and simplicity of the method. However, detailed fits to SN optical spectra, based on tailored models for individual SN observations, offers a promising approach to obtaining distances with 10-20% accuracy, either through the EPM or a la Baron et al. (2004).

Using Sloan Digital Sky Survey (SDSS) quasar spectra taken at multiple epochs, we find that the composite flux density differences in the rest frame wavelength range 1300-6000 AA can be fit by a standard thermal accretion disk model where the accretion rate has changed from one epoch to the next (without considering additional continuum emission components). The fit to the composite residual has two free parameters: a normalizing constant and the average characteristic temperature $\bar{T}^*$. In turn the characteristic temperature is dependent on the ratio of the mass accretion rate to the square of the black hole mass. We therefore conclude that most of the UV/optical variability may be due to processes involving the disk, and thus that a significant fraction of the UV/optical spectrum may come directly from the disk.

We apply the semi-analytical analysis of the steady nature of line-driven winds presented in two earlier papers to disk winds driven by the flux distribution of a standard Shakura & Sunyaev (1973) disk for typical cataclysmic variable (CV) parameters. We find that the wind critical point tends to be closer to the disk surface towards the inner disk regions. Our main conclusion, however, is that a line-driven wind, arising from a steady disk flux distribution of a standard Shakura-Sunyaev disk capable of locally supplying the corresponding mass flow, is steady. These results confirm the findings of an earlier paper that studied "simple" flux distributions that are more readily analyzable than those presented here. These results are consistent with the steady velocity nature of outflows observationally inferred for both CVs and quasi-stellar objects (QSOs). We find good agreement with the 2.5D CV disk wind models of Pereyra and collaborators. These results suggest that the likely scenario to account for the wind outflows commonly observed in CVs is the line-driven accretion disk wind scenario, as suggested early-on by Cordova & Mason (1982). For QSOs, these results show that the line-driven accretion disk wind continues to be a promising scenario to account for the outflows detected in broad absorption line (BAL) QSOs, as suggested early-on by Turnshek (1984), and analyzed in detail by Murray et al. (1995).

Due to their high intrinsic brightness, type II supernovae (SN) can be used as lighthouses to constrain distances in the Universe using variants of the Baade-Wesselink method. Based on a large set of CMFGEN models (Hillier & Miller 1998) covering the photospheric phase of type II SN, we study the various concepts entering one such technique, the Expanding Photosphere Method (EPM). We compute correction factors $\xi$ needed to approximate the synthetic Spectral Energy Distribution (SED) with that of a blackbody at temperature $T$. Our $\xi$, although similar, are systematically greater, by $\sim 0.1$, than the values obtained by Eastman et al. (1996). This translate into a systematic enhancement of 10-20% in EPM-distances. We find that line emission and absorption, not directly linked to color temperature variations, can considerably alter the synthetic magnitude, and cause above-unity $\xi$-values. Optically-thick lines, present in the optical at late-times, can also introduce a strong wavelength-dependence of the photospheric radius, invalidating the use of the Baade method. Both the impact of line-blanketing on the SED and the photospheric radius at low $T$ suggest that the EPM is best used at early times, when the outflow is fully ionized and line-opacity mostly confined to the UV range. We also investigate how reliably one can measure the photospheric velocity from P-Cygni line profiles. Contrary to the usually held belief, the velocity at maximum absorption in the P-Cygni trough of optically-thick lines can both {\it overestimate or underestimate} the photospheric velocity, with a magnitude that depends on the SN outflow density gradient and the optical thickness of the line.

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)

We present new calibrations of stellar parameters of O stars at solar metallicity taking non-LTE, wind, and line-blanketing effects into account. Gravities and absolute visual magnitudes are derived from results of recent spectroscopic analyses. Two types of effective temperature scales are derived: one from a compilation based on recent spectroscopic studies of a sample of massive stars - the ``observational scale'' - and the other from direct interpolations on a grid of non-LTE spherically extended line-blanketed models computed with the code CMFGEN (Hillier & Miller 1998) - the ``theoretical scale''. These Teff scales are then further used together with the grid of models to calibrate other parameters (bolometric correction, luminosity, radius, spectroscopic mass and ionising fluxes) as a function of spectral type and luminosity class. Compared to the earlier calibrations of Vacca et al. (1996) the main results are: 1) the effective temperature scales are cooler by 2000 to 8000 K, the reduction being the largest for early supergiants; 2) bolometric corrections are lower by 0.1 mag for a given temperature, and are lower by 0.3 - 0.6 mag for a given spectral type; 3) luminosities are reduced by 0.20 to 0.35 dex; 4) Lyman ionising fluxes are lower by 0.20 to 0.80 dex. We provide calibration tables for each luminosity class and for both types of Teff scales.

The core of the nebula surrounding Eta Carinae has been observed with the VLT Adaptive Optics system NACO and with the interferometer VLTI/MIDI to constrain spatially and spectrally the warm dusty environment and the central object. In particular, narrow-band images at 3.74 and 4.05 micron reveal the butterfly shaped dusty environment close to the central star with unprecedented spatial resolution. A void whose radius corresponds to the expected sublimation radius has been discovered around the central source. Fringes have been obtained in the Mid-IR which reveal a correlated flux of about 100Jy situated 0.3" south-east of the photocenter of the nebula at 8.7 micron, which corresponds with the location of the star as seen in other wavelengths. This correlated flux is partly attributed to the central object, and these observations provide an upper limit for the SED of the central source from 2.2 to 13.5 micron. Moreover, we have been able to spectrally disperse the signal from the nebula itself at PA=318 degree, i.e. in the direction of the bipolar nebula 310 degree) within the MIDI field of view of 3". A large amount of corundum (Al2O3) is discovered, peaking at 0.6-1.2" south-east from the star, whereas the dust content of the Weigelt blobs is dominated b silicates. We discuss the mechanisms of dust formation which are closely related to the geometry of this Butterfly nebulae.

We discuss work toward developing a 2.5D non-LTE radiative transfer code. Our code uses the short characteristic method with modifications to handle realistic 3D wind velocities. We also designed this code for parallel computing facilities by representing the characteristics with an impact parameter vector p. This makes the interpolation in the radiation angles unnecessary and allows for an independent calculation for characteristics represented by different p vectors. The effects of the velocity field are allowed for by increasing, as needed, the number of grid points along a short characteristic. This allows us to accurately map the variation of the opacities and emissivities as a function of frequency and spatial coordinates. In the future we plan to use this transfer code with a 2D non-LTE stellar atmosphere program to self-consistently solve for level populations, the radiation field and temperature structure for stars with winds and without spherical symmetry.

During $\eta$ Car's spectroscopic event in mid-2003, the stellar wind's bright H$\alpha$ and H$\beta$ emission lines temporarily had a distinctive shape unlike that reported on any previous occasion, and particularly unlike the 1997--98 event. Evidently the structure of the wind changed between 1997 and 2003. Combining this with other evidence, we suspect that the star may now be passing through a rapid stage in its recovery from the Great Eruption seen 160 years ago. In any case the data indicate that successive spectroscopic events differ, and the hydrogen line profiles are quantitative clues to the abnormal structure of the wind during a spectroscopic event.

We present a quantitative spectral analysis of five very massive stars in the Arches cluster, located near the Galactic center, to determine stellar parameters, stellar wind properties and, most importantly, metallicity content. The analysis uses a new technique, presented here for the first time, and uses line-blanketed NLTE wind/atmosphere models fit to high-resolution near-infrared spectra of late-type nitrogen-rich Wolf-Rayet stars and OfI+ stars in the cluster. It relies on the fact that massive stars reach a maximum nitrogen abundance that is related to initial metallicity when they are in the WNL phase. We determine the present-day nitrogen abundance of the WNL stars in the Arches cluster to be 1.6% (mass fraction) and constrain the stellar metallicity in the cluster to be solar. This result is invariant to assumptions about the mass-luminosity relationship, the mass-loss rates, and rotation speeds. In addition, from this analysis, we find the age of the Arches cluster to be 2-2.5Myr, assuming coeval formation.

We report the first detection of X-ray emission associated with the Homunculus Nebula which surrounds the supermassive star Eta Carinae. The emission is characterized by a temperature in excess of 100 MK, and is consistent with scattering of the time-delayed X-ray flux associated with the star. The nebular emission is bright in the northwestern lobe and near the central regions of the Homunculus, and fainter in the southeastern lobe. We also report the detection of an unusually broad Fe K fluorescent line, which may indicate fluorescent scattering off the wind of a companion star or some other high velocity outflow. The X-ray Homunculus is the nearest member of the small class of Galactic X-ray reflection nebulae, and the only one in which both the emitting and reflecting sources are distinguishable.

We present quantitative studies of 8 late O and early B-type supergiants in the Magellanic Clouds using far-ultraviolet FUSE, ultraviolet IUE/HST and optical VLT-UVES spectroscopy. Temperatures, mass-loss rates and CNO abundances are obtained using the non-LTE, spherical, line-blanketed model atmosphere code of Hillier & Miller (1998). We support recent results for lower temperatures of OB-type supergiants as a result of stellar winds and blanketing, which amounts to ~2000 K at B0 Ia. In general, H$\alpha$ derived mass-loss rates are consistent with UV and far-UV spectroscopy, although from consideration of the SIV $\lambda\lambda$1063-1073 doublet, clumped winds are preferred over homogenous models. AV 235 (B0 Iaw) is a notable exception, which has an unusually strong H$\alpha$ profile that is inconsistent with the other Balmer lines and UV wind diagnostics. We also derive CNO abundances for our sample, revealing substantial nitrogen enrichment, with carbon and oxygen depletion. Our results are supported by comparison with the Galactic supergiant HD 2905 (BC0.7 Ia) for which near-solar CNO abundances are obtained. This bolsters previous suggestions that ``normal'' OB-type supergiants exhibit atmospheric compositions indicative of partical CNO processing.

We present a quantitative study of the stellar and wind properties of massive stars in the compact star forming region SMC-N81. Non-LTE spherically extended atmosphere models including line-blanketing have been computed with the code CMFGEN (Hillier & Miller 1998) and the analysis of UV STIS spectra lead to the following results: 1) The SMC-N81 components are young (~ 0-4 Myrs) O stars with Teff compatible with medium to late subtypes and with luminosities lower than average Galactic O dwarfs, rendering them possible ZAMS candidates; 2) The winds are extremely weak: the mass loss rates (~1e-8/1e-9 Msol/yr) and the modified wind momenta are 1 to 2 orders of magnitude lower than observed so far for Galactic stars, and are lower than the predictions of the most recent hydrodynamical models. The reasons for such weak winds are investigated with special emphasis on the modified wind momenta: 1) There may be a break-down of the wind momentum - luminosity relation (WLR) for dwarf stars at low luminosity (log L/Lsol < 5.5). However, reasons for such a breakdown remain unknown; 2) The slope of the WLR may be steeper at low metallicity. This is predicted by the radiation driven wind theory, but the current hydrodynamical simulations do not show any change of the slope at SMC metallicity. Moreover, there are indications that some Galactic objects have wind momenta similar to those of the SMC stars; 3) Decoupling in the atmosphere may lead to multicomponent winds, but various tests indicate that this is not likely to be the case for the SMC-N81 stars. Hence, the origin of the weakness of the wind observed in these objects remains unknown. We suggest that this weakness may be linked with the youth of these stars and represents possibly the onset of stellar winds in recently formed massive stars.

We perform an analytic investigation of the stability of line-driven disk winds, independent of hydrodynamic simulations. Our motive is to determine whether or not line-driven disk winds can account for the wide/broad UV resonance absorption lines seen in cataclysmic variables (CVs) and quasi-stellar objects (QSOs). In both CVs and QSOs observations generally indicate that the absorption arising in the outflowing winds has a steady velocity structure on time scales exceeding years (for CVs) and decades (for QSOs). However, published results from hydrodynamic simulations of line-driven disk winds are mixed, with some researchers claiming that the models are inherently unsteady, while other models produce steady winds. The analytic investigation presented here shows that if the accretion disk is steady, then the line-driven disk wind emanating from it can also be steady. In particular, we show that a gravitational force initially increasing along the wind streamline, which is characteristic of disk winds, does not imply an unsteady wind. The steady nature of line-driven disk winds is consistent with the 1D streamline disk-wind models of Murray and collaborators and the 2.5D time-dependent models of Pereyra and collaborators. This paper emphasizes the underlying physics behind the steady nature of line-driven disk winds using mathematically simple models that mimic the disk environment

We present new high angular resolution observations at near-IR wavelengths of the core of the Luminous Blue Variable Eta Carinae, using NAOS-CONICA at the VLT and VINCI at the VLT Interferometer (VLTI). The latter observations provide spatial information on a scale of 5 milli-arcsec or ~11 AU at the distance of Eta Carinae. The present-day stellar wind of Eta Carinae is resolved on a scale of several stellar radii. Assuming spherical symmetry, we find a mass loss rate of 1.6x10^{-3} M_sun/yr and a wind clumping factor of 0.26. The VLTI data taken at a baseline of 24 meter show that the object is elongated with a de-projected axis ratio of approximately 1.5; the major axis is aligned with that of the large bi-polar nebula that was ejected in the 19th century. The most likely explanation for this observation is a counter-intuitive model in which stellar rotation near the critical velocity causes enhanced mass loss along the rotation axis. This results from the large temperature difference between pole and equator in rapidly rotating stars. Eta Carinae must rotate in excess of 90 per cent of its critical velocity to account for the observed shape. The large outburst may have been shaped in a similar way. Our observations provide strong support for the existence of a theoretically predicted rotational instability, known as the Omega limit.

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)

We use ultraviolet space-based (FUSE, HST) and optical/IR ground-based (2.3m MSSSO, NTT) spectroscopy to determine the physical parameters of six WC4-type Wolf-Rayet stars in the Large Magellanic Cloud. Stellar parameters are revised significantly relative to Grafener et al. (1998) based on improved observations and more sophisticated model atmosphere codes, which account for line blanketing and clumping. We find that stellar luminosities are revised upwards by up to 0.4 dex, with surface abundances spanning a lower range of 0.1<C/He<0.35 (20-45% carbon by mass) and O/He<0.06 (<10% oxygen by mass). Relative to Galactic WC5-8 stars at known distance, and analysed in a similar manner, LMC WC4 stars possess systematically higher stellar luminosities, 0.2dex lower wind densities, yet a similar range of surface chemistries. We illustrate how the classification CIII 5696 line is extremely sensitive to wind density, such that this is the principal difference between the subtype distribution of LMC and Galactic early-type WC stars. Temperature differences do play a role, but carbon abundance does not affect WC spectral types. We illustrate the effect of varying temperature and mass-loss rate on the WC spectral type for HD32257 (WC4, LMC) and HD156385 (WC7, Galaxy) which possess similar abundances and luminosities. Using the latest evolutionary models, pre-supernova stellar masses in the range 11-19 Mo are anticipated for LMC WC4 stars, with 7-14 Mo for Galactic WC stars with known distances. These values are consistent with pre-cursors of bright Type-Ic supernovae such as SN1998bw (alias GRB 980425) for which a minimum total mass of C and O of 14Mo has been independently derived.

We have undertaken quantitative analysis of four LMC and SMC O4--9.7 extreme supergiants using far-ultraviolet FUSE, ultraviolet IUE/HSTand optical VLT UVES spectroscopy. Extended, non-LTE model atmospheres that allow for the consistent treatment of line blanketing are used to analyse wind & photospheric spectral features simultaneously. Using Halpha to constrain dM/dt, HeI-II photospheric lines reveal T_eff's which are systematically (5-7.5kK) & substantially (15-20%) lower than previously derived from unblanketed, plane-parallel, non-LTE photospheric studies. We have confidence in these revisions, since derived T_eff's generally yield consistent fits across the entire 912-7000Ang observed spectral range. In particular, we are able to resolve previous UV-optical T_eff discrepancies. `Of' classification criteria are directly linked to (strong) nitrogen enrichment (via NIII 4097) & (weak) carbon depletion (via CIII 4647-51), providing evidence for mixing of unprocessed and CNO processed material at their stellar surfaces. Oxygen abundances are more difficult to constrain, except via OII lines in the O9.7 supergiant for which it is also found to be somewhat depleted. The effect of wind clumping is also investigated, for which PV 1118-28 potentially provides a useful diagnostic in O-star winds. Revised stellar properties affect existing calibrations of (i) Lyman continuum photons - a factor of two lower for the O4 supergiant; and (ii) kinetic energy released into the ISM by O supergiants. Our results also have importance for the calibration of the wind momentum-luminosity relationship for OB stars.

A newly developed 3-D Monte Carlo model is used, in conjunction with a multi-line non-LTE radiative transfer model, to determine the mass-loss rate of the Wolf-Rayet (W-R) star in the massive binary \object{V444 Cyg} (WN5+O6). This independent estimate of mass-loss rate is attained by fitting the observed \HeI (5876) \AA and \HeII (5412) \AA line profiles, and the continuum light curves of three Stokes parameters ((I, Q, U)) in the (V) band simultaneously. The high accuracy of our determination arises from the use of many observational constraints, and the sensitivity of the continuum polarization to the mass-loss rate. Our best fit model suggests that the mass-loss rate of the system is (\dot{M}_{\WR}=0.6(\pm 0.2) \times 10^{-5} M_{\sun} \mathrm{yr}^{-1} ), and is independent of the assumed distance to \object{V444 Cyg}. The fits did not allow a unique value for the radius of the W-R star to be derived. The range of the volume filling factor for the W-R star atmosphere is estimated to be in the range of 0.050 (for $R_{\WR}=5.0 R_{\sun}$) to 0.075 (for $R_{\WR}=2.5 R_{\sun}$). We also found that the blue-side of \HeI (5876 ) \AA and \HeII (5412) \AA lines at phase 0.8 is relatively unaffected by the emission from the wind-wind interaction zone and the absorption by the O-star atmosphere; hence, the profiles at this phase are suitable for spectral line fittings using a spherical radiative transfer model.

We rediscuss the temperature of O dwarfs based on new non-LTE line blanketed atmosphere models including stellar winds computed with the CMFGEN code of Hillier & Miller (1998). Compared to the latest calibration of Vacca et al. (1996), the inclusion of line blanketing leads to lower effective temperatures, typically by 4000 to 1500 K for O3 to O9.5 dwarf stars. The dependence of the Teff-scale on stellar and model parameters - such as mass loss, microturbulence, and metallicity - is explored, and model predictions are compared to optical observations of O stars. Even for an SMC metallicity we find a non-negligible effect of line blanketing on the Teff-scale. The temperature reduction implies downward revisions of luminosities by 0.1 dex and Lyman continuum fluxes Q0 by approximately 40% for dwarfs of a given spectral type.

We have developed a 3-D Monte Carlo radiative transfer model which computes line and continuum polarization variability for a binary system with an optically thick non-axisymmetric envelope. This allows us to investigate the complex (phase-locked) line and continuum polarization variability features displayed by many massive binaries: W-R+O, O+O, etc. An 8-way tree data structure constructed via a ``cell-splitting'' method allows for high precision with efficient use of computer resources. The model is not restricted to binary systems; it can easily be adapted to a system with an arbitrary density distribution and large density gradients. As an application to a real system, the phase dependent Stokes parameters (I, Q, U) and the phase dependent He I (5876) profiles of the massive binary system V444 Cyg (WN5+O6 III-V) are computed.

A quantitative study of the Luminous Blue Variable NGC 2363-V1 in the Magellanic galaxy NGC 2366 (D = 3.44 Mpc) is presented, based on ultraviolet and optical HST/STIS spectroscopy. Contemporary WFPC2 and WHT imaging reveals a modest V-band brightness increase of ~ 0.2 mag per year between 1996 January and 1997 November, reaching V=17.4 mag, corresponding to Mv=-10.4 mag. Subsequently, V1 underwent a similar decrease in V-band brightness, together with a UV brightening of 0.35 mag from 1997 November to 1999 November. The optical spectrum of V1 is dominated by H emission lines, with Fe II, He I and Na I also detected. In the ultraviolet, a forest of Fe absorption features and numerous absorption lines typical of mid-B supergiants are observed. From a spectral analysis with the non-LTE, line-blanketed code of Hillier & Miller (1998), we derive stellar parameters of T*=11kK, R*=420Ro, log(L/Lo)=6.35 during 1997 November, and T*=13kK, R*=315Ro, log(L/Lo)=6.4 for 1999 July. The wind properties of V1 are also exceptional, with Mdot ~ 4.4 x 10e-4 Mo/yr and $v_{\infty} ~ $300 km/s, allowing for a clumped wind (filling factor = 0.3), and assuming H/He ~ 4 by number. The presence of Fe lines in the UV and optical spectrum of V1 permits an estimate of the heavy elemental abundance of NGC 2363 from our spectral synthesis. Although some deficiencies remain, allowance for charge exchange reactions in our calculations supports a SMC-like metallicity, that has previously been determined for NGC 2363 from nebular oxygen diagnostics. Considering a variety of possible progenitor stars, V1 has definitely undergone a giant eruption, with a substantial increase in stellar luminosity, radius, and almost certainly mass-loss rate.

We present Far-Ultraviolet Spectroscopic Explorer (FUSE) spectroscopy of Sand 2, a LMC WO-type Wolf-Rayet star, revealing the OVI resonance P Cygni doublet at 1032-38A. These data are combined with HST/FOS ultraviolet and Mt Stromlo 2.3m optical spectroscopy, and analysed using a spherical, non-LTE, line-blanketed code. Our study reveals exceptional stellar parameters: T*=150,000K, v_inf=4100 km/s, log (L/Lo)=5.3, and Mdot=10^-5 Mo/yr if we adopt a volume filling factor of 10%. Elemental abundances of C/He=0.7+-0.2 and O/He=0.15(-0.05+0.10) by number qualitatively support previous recombination line studies. We confirm that Sand 2 is more chemically enriched in carbon than LMC WC stars, and is expected to undergo a supernova explosion within the next 50,000 yr.

In this paper we derive stellar parameters for the Wolf-Rayet star in the gamma Velorum binary system (WR11), from a detailed non-LTE model of its optical and infrared spectra. Compared to the study of Schaerer et al., the parameters of the WC8 star are revised to a hotter effective temperature (Teff~57kK),a higher luminosity (log(L/Lsun = 5.00), and a lower mass loss rate (log(Mdot / Msun/yr) = -5.0 - using a 10% clumping filling factor). These changes lead to a significant decrease in wind efficiency number, from 144 to 7, so that the driving mechanism of the wind of this WR star may be simply radiation pressure on lines. The derived spectroscopic luminosity is found to be 40% lower than that derived by De Marco & Schmutz through the mass-luminosity relationship for WR stars (log(L/Lsun = 5.2). The paper furthermore presents a comparison of the independently-developed modelling programs, CMFGEN and ISA-Wind . Overall, there seems to be very reasonable agreement between the derived parameters for WR11, except for the carbon content, which is 2 times higher for CMFGEN (C/He=0.15 vs. 0.06, by number).The comparison also confirms a disparity in the predicted flux at lambda<$400A, found by Crowther et al., which will have effects on several nebular line strengths. The paper also presents the first independent check of the photon loss mechanism proposed by Schmutz. We conclude that, not only is it important to include very many lines to realistically model line blanketing, but in particular those ones that critically interact with strong resonance lines (e.g. HeII lambda303.78). The inclusion of these latter lines may significantly alter the wind ionization structure.

Neon abundances are derived in four Galactic WC stars -- gamma Vel (WR11, WC8+O7.5III), HD156385 (WR90, WC7), HD192103 (WR135, WC8), and WR146 (WC5+O8) - using mid-infrared fine structure lines obtained with ISO/SWS. Stellar parameters for each star are derived using a non-LTE model atmospheric code (Hillier & Miller 1998) together with ultraviolet (IUE), optical (INT, AAT) and infrared (UKIRT, ISO) spectroscopy. In the case of gamma Vel, we adopt results from De Marco et al. (2000), who followed an identical approach. ISO/SWS datasets reveal the [NeIII] 15.5um line in each of our targets, while [NeII] 12.8um, [SIV] 10.5um and [SIII] 18.7um are observed solely in gamma Vel. Using a method updated from Barlow et al. (1988) to account for clumped winds, we derive Ne/He=3-4x10^-3 by number, plus S/He=6x10^-5 for gamma Vel. Neon is highly enriched, such that Ne/S in gamma Vel is eight times higher than cosmic values. However, observed Ne/He ratios are a factor of two times lower than predictions of current evolutionary models of massive stars. An imprecise mass-loss and distance were responsible for the much greater discrepancy in neon content identified by Barlow et al. Our sample of WC5--8 stars span a narrow range in T* (=55--71kK), with no trend towards higher temperature at earlier spectral type, supporting earlier results for a larger sample by Koesterke & Hamann (1995). Stellar luminosities range from 100,000 to 500,000 Lo, while 10^-5.1 < Mdot/(Mo/yr) < 10^-4.5, adopting clumped winds, in which volume filling factors are 10%. In all cases, wind performance numbers are less than 10, significantly lower than recent estimates. Carbon abundances span 0.08 < C/He < 0.25 by number, while oxygen abundances remain poorly constrained.

We use WR124 (WN8h) and its associated nebula M1-67, to test theoretical non-LTE models for Wolf-Rayet (WR) stars. Lyman continuum ionizing flux distributions derived from a stellar analysis of WR124, are compared with nebular properties via photo-ionization modelling. Our study demonstrates the significant role that line blanketing plays in affecting the Lyman ionizing energy distribution of WR stars, of particular relevance to the study of HII regions containing young stellar populations. We confirm previous results that non-line blanketed WR energy distributions fail to explain the observed nebular properties of M1-67, such that the predicted ionizing spectrum is too hard. A line blanketed analysis of WR124 is carried out using the method of Hillier & Miller (1998), with stellar properties in accord with previous results, except that the inclusion of clumping in the stellar wind reduces its wind performance factor to only approx2. The ionizing spectrum of the line blanketed model is much softer than for a comparable temperature unblanketed case, such that negligible flux is emitted with energy above the HeI 504 edge. Photo-ionization modelling, incorporating the observed radial density distribution for M1-67 reveals excellent agreement with the observed nebular electron temperature, ionization balance and line strengths. An alternative stellar model of WR124 is calculated, following the technique of de Koter et al. (1997), augmented to include line blanketing following Schmutz et al. (1991). Good consistency is reached regarding the stellar properties of WR124, but agreement with the nebular properties of M1-67 is somewhat poorer than for the Hillier & Miller code.