Source author record

D. John Hillier

D. John Hillier appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.HE astro-ph.CO astro-ph.GA astro-ph

Catalog footprint

What is connected

53works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Eta Carinae: an evolving view of the central binary, its interacting winds and its foreground ejecta

FUV spectra of Eta Car, recorded across two decades with HST/STIS, document multiple changes in resonant lines caused by dissipating extinction in our line of sight. The FUV flux has increased nearly ten-fold which has led to increased ionization of the multiple shells within the Homunculus and photo-destruction of molecular hydrogen. Comparison of observed resonant line profiles with CMFGEN model profiles allows separation of wind-wind collision and shell absorptions from the primary wind, P Cygni profiles.The dissipating occulter preferentially obscured the central binary and interacting winds relative to the very extended primary wind. We are now able to monitor changes in the colliding winds with orbital phase. High velocity transient absorptions occurred across the most recent periastron passage, indicating acceleration of the primary wind by the secondary wind which leads to a downstream, high velocity bowshock that is newly generated every orbital period. There is no evidence of changes in the properties of the binary winds.

preprint2022arXiv

Modeling the signatures of interaction in Type II supernovae: UV emission, high-velocity features, broad-boxy profiles

Because mass loss is a fundamental phenomenon in massive stars, interaction with circumstellar material (CSM) should be universal in core-collapse supernovae (SNe). Leaving aside the extreme CSM density, extent, or mass typically encountered in Type IIn SNe, we investigate the diverse long-term radiative signatures of interaction between a Type II SN ejecta and CSM corresponding to mass loss rates up to 10$^{-3}$ $M_{\odot}$ yr$^{-1}$. Because these CSM are relatively tenuous and optically-thin to electron-scattering beyond a few stellar radii, radiation hydrodynamics is not essential and one may treat the interaction directly as an additional power source in the non-local thermodynamic equilibrium radiative transfer problem. The CSM accumulated since shock breakout forms a dense shell in the outer ejecta and leads to high-velocity absorption features in spectral lines, even for negligible shock power. Besides Balmer lines, such features may appear in NaID, HeI lines etc. A stronger interaction strengthens the continuum flux (preferentially in the UV), quenches the absorption of P-Cygni profiles, boosts the MgII $λλ$ $2795,2802$ doublet, and fosters the production of a broad boxy H$α$ emission component. The rise in ionization in the outer ejecta may quench some lines (e.g., the CaII near-infrared triplet). The interaction power emerges preferentially in the UV, in particular at later times, shifting the optical color to the blue, but increasing modestly the optical luminosity. Strong thermalization and clumping seem to be required to make an interaction superluminous in the optical. The UV range contains essential signatures that provide critical constraints to infer the mass loss history and inner workings of core-collapse SN progenitors at death.

preprint2022arXiv

New hydrodynamic solutions for line-driven winds of hot massive stars using Lambert $W$-function

Hot massive stars present strong stellar winds that are driven by absorption, scattering and re\-emission of photons by the ions of the atmosphere (\textit{line-driven winds}). A better comprehension of this phenomenon, and a more accurate calculation of hydrodynamics and radiative acceleration is required to reduce the number of free parameters in spectral fitting, to determine accurate wind parameters such as mass-loss rates and velocity profiles. We use the non-LTE model-atmosphere code CMFGEN to numerically solve the radiative transfer equation in the stellar atmosphere and to calculate the radiative acceleration $g_\text{rad}(r)$. Under the assumption that the radiative acceleration depends only on the radial coordinate, we solve analytically the equation of motion by means of the Lambert $W$-function. An iterative procedure between the solution of the radiative transfer and the equation of motion is executed in order to obtain a final self-consistent velocity field that is no longer based on any $β$-law. We apply the Lambert-procedure to three O supergiant stars ($ζ$-Puppis, HD~165763 and $α$-Cam) and discuss the Lambert-solutions for the velocity profiles. It is found that, even without recalculation of the mass-loss rate, the Lambert-procedure allows the calculation of consistent velocity profiles that reduce the number of free parameters when a spectral fitting using CMFGEN is performed. Synthetic spectra calculated from our Lambert-solutions show significant differences compared to the initial $β$-law CMFGEN models. The results indicate the importance of consistent velocity profile calculation in the CMFGEN code and its usage in a fitting procedure and interpretation of observed spectra.

preprint2022arXiv

NICER X-ray Observations of Eta Carinae During its Most Recent Periastron Passage

We report high-precision X-ray monitoring observations in the 0.4-10 keV band of the luminous, long-period colliding-wind binary Eta Carinae up to and through its most recent X-ray minimum/periastron passage in February 2020. Eta Carinae reached its observed maximum X-ray flux on 7 January 2020, at a flux level of $3.30 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$, followed by a rapid plunge to its observed minimum flux, $0.03 \times 10^{-10}$ ergs s$^{-1}$ cm$^{-2}$ near 17 February 2020. The NICER observations show an X-ray recovery from minimum of only $\sim$16 days, the shortest X-ray minimum observed so far. We provide new constraints of the "deep" and "shallow" minimum intervals. Variations in the characteristic X-ray temperature of the hottest observed X-ray emission indicate that the apex of the wind-wind "bow shock" enters the companion's wind acceleration zone about 81 days before the start of the X-ray minimum. There is a step-like increase in column density just before the X-ray minimum, probably associated with the presence of dense clumps near the shock apex. During recovery and after, the column density shows a smooth decline, which agrees with previous $N_{H}$ measurements made by SWIFT at the same orbital phase, indicating that changes in mass-loss rate are only a few percent over the two cycles. Finally, we use the variations in the X-ray flux of the outer ejecta seen by NICER to derive a kinetic X-ray luminosity of the ejecta of $\sim 10^{41}$ ergs s$^{-1}$ near the time of the "Great Eruption'.

preprint2022arXiv

Stable nickel production in Type Ia supernovae: A smoking gun for the progenitor mass?

At present, there are strong indications that white dwarf (WD) stars with masses well below the Chandrasekhar limit (MCh ~ 1.4 Msun) contribute a significant fraction of SN Ia progenitors. The relative fraction of stable iron-group elements synthesized in the explosion has been suggested as a possible discriminant between MCh and sub-MCh events. In particular, it is thought that the higher-density ejecta of MCh WDs, which favours the synthesis of stable isotopes of nickel, results in prominent [Ni II] lines in late-time spectra. We study the explosive nucleosynthesis of stable nickel in SNe Ia resulting from MCh and sub-MCh progenitors. We explore the potential for lines of [Ni II] at 7378 Å and 1.94 microns in late-time spectra to serve as a diagnostic of the exploding WD mass, using nonlocal thermodynamic equilibrium radiative-transfer simulations with the CMFGEN code. We find that the radiative proton-capture reaction 57Co(p,gamma)58Ni is the dominant production mode for 58Ni in both MCh and sub-MCh models, while the alpha-capture reaction on 54Fe has a negligible impact on the final 58Ni yield. More importantly, we demonstrate that the lack of [Ni II] lines in late-time spectra of sub-MCh events is not always due to an under-abundance of stable Ni; rather, it results from the higher ionization of Ni in the inner ejecta. Conversely, the strong [Ni II] lines predicted in our 1D MCh models are completely suppressed when 56Ni is sufficiently mixed with the innermost layers, which are rich in stable iron-group elements. [Ni II] lines in late-time SN Ia spectra have a complex dependency on the abundance of stable Ni, which limits their use in distinguishing among MCh and sub-MCh progenitors. However, we argue that a low-luminosity SN Ia displaying strong [Ni II] lines would most likely result from a Chandrasekhar-mass progenitor. [Abridged]

preprint2022arXiv

The orbital kinematics of eta Carinae over three periastra with a possible detection of the elusive secondary's motion

The binary eta Carinae is the closest example of a very massive star, which may have formed through a merger during its Great Eruption in the mid-nineteenth century. We aimed to confirm and improve the kinematics using a spectroscopic data set taken with the CTIO 1.5 m telescope over the time period of 2008-2020, covering three periastron passages of the highly eccentric orbit. We measure line variability of H-alpha and H-beta, where the radial velocity and orbital kinematics of the primary star were measured from the H-beta emission line using a bisector method. At phases away from periastron, we observed the He II 4686 emission moving opposite the primary star, consistent with a possible Wolf-Rayet companion, although with a seemingly narrow emission line. This could represent the first detection of emission from the companion.

preprint2022arXiv

WO-Type Wolf-Rayet Stars: the Last Hurrah of Massive Star Evolution

Are WO-type Wolf Rayet (WR) stars in the final stage of massive star evolution before core-collapse? Although WC- and WO-type WRs have very similar spectra, WOs show a much stronger O VI $λλ$3811,34 emission-line feature. This has usually been interpreted to mean that WOs are more oxygen rich than WCs, and thus further evolved. However, previous studies have failed to model this line, leaving the relative abundances uncertain, and the relationship between the two types unresolved. To answer this fundamental question, we modeled six WCs and two WOs in the LMC using UV, optical, and NIR spectra with the radiative transfer code CMFGEN in order to determine their physical properties. We find that WOs are not richer in oxygen; rather, the O VI feature is insensitive to the abundance. However, the WOs have a significantly higher carbon and lower helium content than the WCs, and hence are further evolved. A comparison of our results with single-star Geneva and binary BPASS evolutionary models show that while many properties match, there is more carbon and less oxygen in the WOs than either set of evolutionary model predicts. This discrepancy may be due to the large uncertainty in the $^{12}$C+$^4$He$\rightarrow^{16}$O nuclear reaction rate; we show that if the Kunz et al. rate is decreased by a factor of 25-50%, then there would be a good match with the observations. It would also help explain the LIGO/VIRGO detection of black holes whose masses are in the theoretical upper mass gap.

preprint2021arXiv

The Physical Parameters of Four WC-type Wolf-Rayet Stars in the Large Magellanic Cloud: Evidence of Evolution

We present a spectral analysis of four LMC WC-type Wolf-Rayet (WR) stars (BAT99-8, BAT99-9, BAT99-11, and BAT99-52) to shed light on two evolutionary questions surrounding massive stars. The first is: are WO-type WR stars more oxygen enriched than the WC-type stars, indicating further chemical evolution, or are the strong high-excitation oxygen lines in the WO-type stars an indication of higher temperatures. This study will act as a baseline for answering the question of where WO-type stars fall in WR evolution. Each star's spectrum, extending from 1100~Å to 25000~Å, was modeled using \cmfgen\ to determine the star's physical properties such as luminosity, mass-loss rate, and chemical abundances. The oxygen abundance is a key evolutionary diagnostic, and with higher resolution data and an improved stellar atmosphere code, we found the oxygen abundance to be up to a factor of 5 lower than previous studies. The second evolutionary question revolves around the formation of WR stars: Do they evolve by themselves or is a close companion star necessary for their formation? Using our derived physical parameters, we compared our results to the Geneva single-star evolutionary models and the BPASS binary evolutionary models. We found that both the Geneva solar metallicity models and BPASS LMC metallicity models are in agreement with the four WC-type stars, while the Geneva LMC metallicity models are not. Therefore, these four WC4 stars could have been formed either via binary or single-star evolution.

preprint2019arXiv

Understanding nebular spectra of Type Ia supernovae

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] $λλ$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] $λλ$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.

preprint2016arXiv

Inferring supernova IIb/Ib/Ic ejecta properties from light curves and spectra: Correlations from radiative-transfer models

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Å at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772Å 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$μ$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.

preprint2016arXiv

Modelling the Central Constant Emission X-ray component of eta Carinae

The X-ray emission of $η$ 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 $ω\approx252^\circ$, and bound/converge on the observation for $ω\approx90^\circ$. This suggests that $η$ 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.

preprint2016arXiv

Models of interacting supernovae and their spectral diversity

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.

preprint2016arXiv

To $v_\infty$ and beyond! The He I absorption variability across the 2014.6 periastron passage of $η$ Carinae

We have monitored the massive binary star $η$ 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 $λ$10830 Å 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.

preprint2015arXiv

A one-dimensional Chandrasekhar-mass delayed-detonation model for the broad-lined Type Ia supernova 2002bo

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.

preprint2015arXiv

Measuring nickel masses in Type Ia supernovae using cobalt emission in nebular phase spectra

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.

preprint2015arXiv

Numerical simulations of super-luminous supernovae of type IIn

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.

preprint2015arXiv

Radiative-transfer models for supernovae IIb/Ib/Ic from binary-star progenitors

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

preprint2015arXiv

X-ray, UV and optical analysis of supergiants: $ε$ Ori

We present a multi-wavelength (X-ray to optical) analysis, based on non-local thermodynamic equilibrium photospheric+wind models, of the B0 Ia-supergiant: $ε$~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 $λλ$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 ($ρ_{cl}/ρ_{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: $ε$ Ori is slightly enhanced in nitrogen and depleted in carbon and oxygen, evidence for CNO processed material.

preprint2014arXiv

[CoIII] versus NaID in type Ia supernova spectra

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.

preprint2014arXiv

A Modern Search for Wolf-Rayet Stars in the Magellanic Clouds: First Results

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.

preprint2014arXiv

A New Class of Wolf-Rayet Stars: WN3/O3s

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.

preprint2014arXiv

Constraints on the explosion mechanism and progenitors of type Ia supernovae

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.

preprint2014arXiv

Critical ingredients of supernova Ia radiative-transfer modeling

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

preprint2014arXiv

One-dimensional non-LTE time-dependent radiative transfer of an He-detonation model and the connection to faint and fast-decaying supernovae

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.

preprint2014arXiv

The HD5980 multiple system: Masses and evolutionary status

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.

preprint2013arXiv

A Bake-Off Between CMFGEN and FASTWIND: Modeling the Physical Properties of SMC and LMC O-type Stars

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.

preprint2013arXiv

Massive stars at low metallicity: Evolution and surface abundances of O dwarfs in the SMC

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.

preprint2013arXiv

Type II-Plateau supernova radiation: dependencies on progenitor and explosion properties

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.

preprint2012arXiv

Non-thermal excitation and ionization in supernovae

We incorporate non-thermal excitation and ionization processes arising from non-thermal electrons that result from γ-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 μ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α. 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.

preprint2012arXiv

On the influence of the companion star in Eta Carinae: 2D radiative transfer modeling of the ultraviolet and optical spectra

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.

preprint2012arXiv

On the nature of supernovae Ib and Ic

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.

preprint2012arXiv

One-dimensional delayed-detonation models of Type Ia supernovae: Confrontation to observations at bolometric maximum

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.

preprint2012arXiv

Radiative Properties of Pair-instability Supernova Explosions

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.

preprint2012arXiv

Super-luminous supernovae: 56Ni power versus magnetar radiation

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.

preprint2012arXiv

The influence of rotation on optical emission profiles of O 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.

preprint2012arXiv

Time Dependent Radiative Transfer Calculations for Supernovae

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

preprint2011arXiv

A Method for the Study of Accretion Disk Emission in Cataclysmic Variables I: The Model

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.

preprint2011arXiv

A Spectropolarimetric Comparison of the Type II-Plateau Supernovae SN 2008bk and SN 2004dj

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.

preprint2011arXiv

Core-collapse explosions of Wolf-Rayet stars and the connection to type IIb/Ib/Ic supernovae

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.

preprint2011arXiv

Multi-Wavelength Implications of the Companion Star in Eta Carinae

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.

preprint2011arXiv

On the nature of the prototype LBV AG Carinae II. Witnessing a massive star evolving close to the Eddington and bistability limits

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.

preprint2011arXiv

Synthetic line and continuum linear-polarisation signatures of axisymmetric type II supernova ejecta

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.

preprint2011arXiv

The Atomic Physics Underlying the Spectroscopic Analysis of Massive Stars and Supernovae

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.

preprint2010arXiv

A mass-loss rate determination for zeta Puppis from the quantitative analysis of X-ray emission line profiles

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.

preprint2010arXiv

Is Eta Carinae a fast rotator, and how much does the companion influence the inner wind structure?

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.

preprint2010arXiv

Non-LTE time-dependent spectroscopic modelling of type II-plateau supernovae from the photospheric to the nebular phase: case study for 15 and 25Msun progenitor stars

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.

preprint2010arXiv

Supernova radiative-transfer modeling: A new approach using non-LTE and full time dependence

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.

preprint2009arXiv

Bona-fide, strong-variable galactic Luminous Blue Variable stars are fast rotators: detection of a high rotational velocity in HR Carinae

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.

preprint2009arXiv

On the nature of the prototype LBV AG Carinae I. Fundamental parameters during visual minimum phases and changes in the bolometric luminosity during the S-Dor cycle

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)

preprint2008arXiv

Metallicity in the Galactic Center: The Quintuplet cluster

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.

preprint2007arXiv

X-ray Spectral Variation of Eta Carinae through the 2003 X-ray Minimum

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.

preprint2004arXiv

Metallicity in the Galactic Center: The Arches cluster

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.

preprint2003arXiv

On the Steady Nature of Line-Driven Disk Winds

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

D. John Hillier

What is connected

Connect this record

See the researcher in context

Building this map preview

53 published item(s)

Eta Carinae: an evolving view of the central binary, its interacting winds and its foreground ejecta

Modeling the signatures of interaction in Type II supernovae: UV emission, high-velocity features, broad-boxy profiles

New hydrodynamic solutions for line-driven winds of hot massive stars using Lambert $W$-function

NICER X-ray Observations of Eta Carinae During its Most Recent Periastron Passage

Stable nickel production in Type Ia supernovae: A smoking gun for the progenitor mass?

The orbital kinematics of eta Carinae over three periastra with a possible detection of the elusive secondary's motion

WO-Type Wolf-Rayet Stars: the Last Hurrah of Massive Star Evolution

The Physical Parameters of Four WC-type Wolf-Rayet Stars in the Large Magellanic Cloud: Evidence of Evolution

Understanding nebular spectra of Type Ia supernovae

Inferring supernova IIb/Ib/Ic ejecta properties from light curves and spectra: Correlations from radiative-transfer models

Modelling the Central Constant Emission X-ray component of eta Carinae

Models of interacting supernovae and their spectral diversity

To $v_\infty$ and beyond! The He I absorption variability across the 2014.6 periastron passage of $η$ Carinae

A one-dimensional Chandrasekhar-mass delayed-detonation model for the broad-lined Type Ia supernova 2002bo

Measuring nickel masses in Type Ia supernovae using cobalt emission in nebular phase spectra

Numerical simulations of super-luminous supernovae of type IIn

Radiative-transfer models for supernovae IIb/Ib/Ic from binary-star progenitors

X-ray, UV and optical analysis of supergiants: $ε$ Ori

[CoIII] versus NaID in type Ia supernova spectra

A Modern Search for Wolf-Rayet Stars in the Magellanic Clouds: First Results

A New Class of Wolf-Rayet Stars: WN3/O3s

Constraints on the explosion mechanism and progenitors of type Ia supernovae

Critical ingredients of supernova Ia radiative-transfer modeling

One-dimensional non-LTE time-dependent radiative transfer of an He-detonation model and the connection to faint and fast-decaying supernovae

The HD5980 multiple system: Masses and evolutionary status

A Bake-Off Between CMFGEN and FASTWIND: Modeling the Physical Properties of SMC and LMC O-type Stars

Massive stars at low metallicity: Evolution and surface abundances of O dwarfs in the SMC

Type II-Plateau supernova radiation: dependencies on progenitor and explosion properties

Non-thermal excitation and ionization in supernovae

On the influence of the companion star in Eta Carinae: 2D radiative transfer modeling of the ultraviolet and optical spectra

On the nature of supernovae Ib and Ic

One-dimensional delayed-detonation models of Type Ia supernovae: Confrontation to observations at bolometric maximum

Radiative Properties of Pair-instability Supernova Explosions

Super-luminous supernovae: 56Ni power versus magnetar radiation

The influence of rotation on optical emission profiles of O stars

Time Dependent Radiative Transfer Calculations for Supernovae

A Method for the Study of Accretion Disk Emission in Cataclysmic Variables I: The Model

A Spectropolarimetric Comparison of the Type II-Plateau Supernovae SN 2008bk and SN 2004dj

Core-collapse explosions of Wolf-Rayet stars and the connection to type IIb/Ib/Ic supernovae

Multi-Wavelength Implications of the Companion Star in Eta Carinae

On the nature of the prototype LBV AG Carinae II. Witnessing a massive star evolving close to the Eddington and bistability limits

Synthetic line and continuum linear-polarisation signatures of axisymmetric type II supernova ejecta

The Atomic Physics Underlying the Spectroscopic Analysis of Massive Stars and Supernovae

A mass-loss rate determination for zeta Puppis from the quantitative analysis of X-ray emission line profiles

Is Eta Carinae a fast rotator, and how much does the companion influence the inner wind structure?

Non-LTE time-dependent spectroscopic modelling of type II-plateau supernovae from the photospheric to the nebular phase: case study for 15 and 25Msun progenitor stars

Supernova radiative-transfer modeling: A new approach using non-LTE and full time dependence

Bona-fide, strong-variable galactic Luminous Blue Variable stars are fast rotators: detection of a high rotational velocity in HR Carinae

On the nature of the prototype LBV AG Carinae I. Fundamental parameters during visual minimum phases and changes in the bolometric luminosity during the S-Dor cycle

Metallicity in the Galactic Center: The Quintuplet cluster

X-ray Spectral Variation of Eta Carinae through the 2003 X-ray Minimum

Metallicity in the Galactic Center: The Arches cluster

On the Steady Nature of Line-Driven Disk Winds