ABSTRACT
The spectra of Wolf–Rayet (WR) stars exhibit strong, broad emission lines that originate in the wind. These winds are radiatively driven and are susceptible to hydrodynamic instabilities ...that result in the formation of clumps. When modelling spectra of WR stars the volume-filling factor (VFF) approach is usually employed to treat clumpy winds. However, it is based on the assumption that the entire wind mass resides in optically thin clumps, which is not necessarily justifiable in dense winds. To test the validity of the VFF approach, we use a previously described method of treating clumping, the ‘shell’ approach, to study line and continuum formation in the dense wind of the WN4 star, HD 50896. Our models indicate that fully intact spherical shells are in tension with observed spectra; a persistent ‘dip’ in emission lines occurs at line centre. Removing this dip requires our models to use ‘broken’ shells – shells that are highly decoherent laterally. This insinuates that the wind of HD 50896, and by extension the winds of other WR stars, are comprised of small laterally confined and radially compressed clumps – clumps smaller than the Sobolev length. We discuss some of the conditions necessary for the VFF approach to be valid.
Over the years, directed surveys and incidental spectroscopy have identified 12 Wolf-Rayet (WR) stars in the Small Magellanic Cloud (SMC) and 139 in the Large Magellanic Cloud (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 nine 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 Be star HD 38489 may have a WN companion.
ABSTRACT
An analysis of the Large Magellanic Cloud (LMC) WC4 star BAT99-9 (HD 32125, FD 4, Brey 7, WS 3) shows that the star still contains photospheric nitrogen. Three N emission features (N v ...λλ1238, 1242, N iv λ1719, and N iv λλ3479–3485) are unambiguously identified in the spectrum. cmfgen models of the star yield an N/C ratio of 0.004 ± 0.002 (by number) and a C/He ratio of $0.15_{-0.05}^{+0.10}$. Due to the similarity of BAT99-9 to other WC4 stars, and the good fit achieved by cmfgen to both the classic WC4 spectrum and the N lines, we argue that the N lines are intrinsic to BAT99-9. An examination of a limited set of rotating models for single-star evolution at LMC and Galactic metallicities shows that a model with a Galactic metallicity (z = 0.014) and a progenitor mass of around 50 M⊙ can have an N/C ratio similar to, or larger than, what we observe for a significant fraction of its lifetime. However, the LMC models (z = 0.006) are inconsistent with the observations. Both the single and binary BPASS models predict that many WC stars can have an N/C ratio similar to, or larger than, what we observe for a significant fraction of their lifetime. While the binary models cover a wider range of luminosities and provide a somewhat better match to BAT99-9, it is not currently possible to rule out BAT99-9 being formed through single-star evolution, given the uncertainties in mass-loss rates, and the treatment of convection and mixing processes.
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 one-dimensional steady-state ...non-local thermodynamic equilibrium radiative transfer calculations of type II SNe at 300 d after explosion. Guided by the results obtained from a large set of stellar evolution simulations, we craft ejecta models for type II SNe from the explosion of a 12, 15, 20, and 25 M⊙ 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 Fe II and O I lines with Ly α and Ly β. Our spectra show a strong sensitivity to 56Ni mixing since it determines where decay power is absorbed. Even at 300 d 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 Ca II line emission. In our set of models, the O I λλ 6300, 6364 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 would weaken the O I doublet flux in the resulting nebular SN II spectrum. This process may occur in nature, with a greater occurrence in higher mass progenitors, and this may explain in part the preponderance of progenitor masses below 17 M⊙ that are inferred from nebular spectra.
Abstract
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), though 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 possesses a metallicity gradient.
We present non-Local Thermodynamic Equilibrium (LTE) time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ∼1051 erg piston-driven ejecta, with ...and without 56Ni, produced from single and binary Wolf-Rayet (WR) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-d long post-breakout plateau with a luminosity of 1-5 × 107 L⊙, visually brighter by ≳10 mag than the progenitor WR star. In our 56Ni-rich models, with ∼3 M⊙ ejecta masses, this plateau precedes a 20 to 30 d long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. A larger ejecta mass or a deeper 56Ni location increases the heat diffusion time and acts to both delay and broaden the light-curve peak. Discriminating between the two effects requires spectroscopic modelling. In low ejecta-mass models with moderate mixing, γ-ray leakage starts as early as ∼50 d after explosion and causes the nebular luminosity to steeply decline by ∼0.02 mag d−1. 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 ≲20 M⊙ 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.001 M⊙, yield the characteristic SN IIb spectral morphology at early times. However at later times, ∼15 d after the explosion, only Hα may remain as a weak absorption feature. Our binary models, characterized by helium surface mass fractions of ≳0.85, systematically show He i 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 present 1D 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
e of 1.7–5.2 M⊙, kinetic energies E
kin of 0.6–5.0 × 1051 erg, and 56Ni masses of 0.05–0.30 M⊙. We find a strong correlation between the 56Ni 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 10 d 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 He i 5875 Å at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The O i 7772 Å 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 ∼1.05 μm feature in all models. The reversal is due to He i alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionization, 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.
We present 1D non-local thermodynamic equilibrium time-dependent radiative transfer simulations of a Chandrasekhar-mass delayed-detonation model which synthesizes 0.51 M⊙ of 56Ni, and confront our ...results to the Type Ia supernova (SN Ia) 2002bo over the first 100 d of its evolution. Assuming only homologous expansion, this same model reproduces the bolometric and multiband 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 M⊙ (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-ultraviolet 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 d 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 Å 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 1D 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.2 × 1051 erg, but characterized by different ejecta masses (2–5 M⊙), 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.32 M⊙ of helium in a CO-rich core of 5.11 M⊙. 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 per cent while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong C i lines around 9000Å prior to maximum is an indicator that the pre-SN star was underabundant in helium. As noted by others, the 1.08μm feature is a complex blend of C i, Mg ii, and He i lines, which makes the identification of He uncertain in SNe Ibc unless other He i lines can be identified. Our models show little scatter in (V − R) colour 10 d 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’.
Infrared observations of the dusty, massive Homunculus Nebula around the luminous blue variable
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
m spectral energy distribution, constructed from legacy ISO observations and new spectroscopy obtained with the
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
m feature. Our preferred model contains nitrides AlN and Si
N
in low abundances. Dust masses range from 0.25 to 0.44
but
≥ 45
in both cases due to an expected high Fe gas-to-dust ratio. The bulk of dust is within a 5″ × 7″ central region. An additional compact feature is detected at 390
m. We obtain
= 2.96 × 10
, 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
CO and
CO
= 5 - 4 through 9 - 8 lines, showing that the abundances are consistent with expectations for CNO-processed material. The
C II line is detected in absorption, which we suspect originates in foreground material at very low excitation temperatures.