Self-consistent modelling of line-driven hot-star winds with Monte Carlo radiation hydrodynamics

Radiative pressure exerted by line interactions is a prominent driver of outflows in astrophysical systems, being at work in the outflows emerging from hot stars or from the accretion discs of cataclysmic variables, massive young stars and active galactic nuclei. In this work, a new radiation hydrodynamical approach to model line-driven hot-star winds is presented. By coupling a Monte Carlo radiative transfer scheme with a finite volume fluid dynamical method, line-driven mass outflows may be modelled self-consistently, benefiting from the advantages of Monte Carlo techniques in treating multiline effects, such as multiple scatterings, and in dealing with arbitrary multidimensional configurations. In this work, we introduce our approach in detail by highlighting the key numerical techniques and verifying their operation in a number of simplified applications, specifically in a series of self-consistent, one-dimensional, Sobolev-type, hot-star wind calculations. The utility and accuracy of our approach are demonstrated by comparing the obtained results with the predictions of various formulations of the so-called CAK theory and by confronting the calculations with modern sophisticated techniques of predicting the wind structure. Using these calculations, we also point out some useful diagnostic capabilities our approach provides. Finally, we discuss some of the current limitations of our method, some possible extensions and potential future applications.

5.9-keV Mn K-shell X-ray luminosity from the decay of 55Fe in Type Ia supernova models

We show that the X-ray line flux of the Mn Kα line at 5.9 keV from the decay of 55Fe is a promising diagnostic to distinguish between Type Ia supernova (SN Ia) explosion models. Using radiation transport calculations, we compute the line flux for two three-dimensional explosion models: a near-Chandrasekhar mass delayed detonation and a violent merger of two (1.1 and 0.9 M⊙) white dwarfs. Both models are based on solar metallicity zero-age main-sequence progenitors. Due to explosive nuclear burning at higher density, the delayed-detonation model synthesizes ˜3.5 times more radioactive 55Fe than the merger model. As a result, we find that the peak Mn Kα line flux of the delayed-detonation model exceeds that of the merger model by a factor of ˜4.5. Since in both models the 5.9-keV X-ray flux peaks five to six years after the explosion, a single measurement of the X-ray line emission at this time can place a constraint on the explosion physics that is complementary to those derived from earlier phase optical spectra or light curves. We perform detector simulations of current and future X-ray telescopes to investigate the possibilities of detecting the X-ray line at 5.9 keV. Of the currently existing telescopes, XMM-Newton/pn is the best instrument for close (≲1-2 Mpc), non-background limited SNe Ia because of its large effective area. Due to its low instrumental background, Chandra/ACIS is currently the best choice for SNe Ia at distances above ˜2 Mpc. For the delayed-detonation scenario, a line detection is feasible with Chandra up to ˜3 Mpc for an exposure time of 106 s. We find that it should be possible with currently existing X-ray instruments (with exposure times ≲5 × 105 s) to detect both of our models at sufficiently high S/N to distinguish between them for hypothetical events within the Local Group. The prospects for detection will be better with future missions. For example, the proposed Athena/X-IFU instrument could detect our delayed-detonation model out to a distance of ˜5 Mpc. This would make it possible to study future events occurring during its operational life at distances comparable to those of the recent supernovae SN 2011fe (˜6.4 Mpc) and SN 2014J (˜3.5 Mpc).

The ejected mass distribution of Type Ia supernovae: a significant rate of non-Chandrasekhar-mass progenitors

The ejected mass distribution of Type Ia supernovae (SNe Ia) directly probes progenitor evolutionary history and explosion mechanisms, with implications for their use as cosmological probes. Although the Chandrasekhar mass is a natural mass scale for the explosion of white dwarfs as SNe Ia, models allowing SNe Ia to explode at other masses have attracted much recent attention. Using an empirical relation between the ejected mass and the light-curve width, we derive ejected masses Mej and 56Ni masses MNi for a sample of 337 SNe Ia with redshifts z <0.7 used in recent cosmological analyses. We use hierarchical Bayesian inference to reconstruct the joint Mej-MNi distribution, accounting for measurement errors. The inferred marginal distribution of Mej has a long tail towards sub-Chandrasekhar masses, but cuts off sharply above 1.4 M⊙. Our results imply that 25-50 per cent of normal SNe Ia are inconsistent with Chandrasekhar-mass explosions, with almost all of these being sub-Chandrasekhar mass; super-Chandrasekhar-mass explosions make up no more than 1 per cent of all spectroscopically normal SNe Ia. We interpret the SN Ia width-luminosity relation as an underlying relation between Mej and MNi, and show that the inferred relation is not naturally explained by the predictions of any single known explosion mechanism.

Locations of peculiar supernovae as a diagnostic of their origins

We put constraints on the properties of the progenitors of peculiar calcium-rich transients using the distribution of locations within their host galaxies. We confirm that this class of transients do not follow the galaxy stellar mass profile and are more likely to be found in remote locations of their apparent hosts. We test the hypothesis that these transients are from low-metallicity progenitors by comparing their spatial distributions with the predictions of self-consistent cosmological simulations that include star formation and chemical enrichment. We find that while metal-poor stars and our transient sample show a consistent preference for large offsets, metallicity alone cannot explain the extreme cases. Invoking a lower age limit on the progenitor helps to improve the match, indicating these events may result from a very old metal-poor population. We also investigate the radial distribution of globular cluster systems, and show that they too are consistent with the class of calcium-rich transients. Because photometric upper limits exist for globular clusters for some members of the class, a production mechanism related to the dense environment of globular clusters is not favoured for the calcium-rich events. However, the methods developed in this paper may be used in the future to constrain the effects of low metallicity on radially distant core-collapse events or help establish a correlation with globular clusters for other classes of peculiar explosions.

SN 2011fu: a type IIb supernova with a luminous double-peaked light curve

We present optical and near-infrared observations of the type IIb supernova (SN) 2011fu from a few days to similar to 300 d after explosion. The SN presents a double-peaked light curve (LC) similar to that of SN 1993J, although more luminous and with a longer cooling phase after the primary peak. The spectral evolution is also similar to SN 1993J's, with hydrogen dominating the spectra to similar to 40 d, then helium gaining strength, and nebular emission lines appearing from similar to 60 d post-explosion. The velocities derived from the P-Cygni absorptions are overall similar to those of other type IIb SNe. We have found a strong similarity between the oxygen and magnesium line profiles at late times, which suggests that these lines are forming at the same location within the ejecta. The hydrodynamical modelling of the pseudo-bolometric LC and the observed photospheric velocities suggest that SN 2011fu was the explosion of an extended star (R similar to 450 R-circle dot), in which 1.3 x 10(51) erg of kinetic energy were released and 0.15 M-circle dot of Ni-56 were synthesized. In addition, a better reproduction of the observed early pseudo-bolometric LC is achieved if a more massive H-rich envelope than for other type IIb SNe is considered (0.3 M-circle dot). The hydrodynamical modelling of the LC and the comparison of our late-time spectra with nebular spectral models for type IIb SNe, point to a progenitor for SN 2011fu with a Zero Age Main Sequence (ZAMS) mass of 13-18 M-circle dot.

SN 2009ip at late times - an interacting transient at+2 years

We present photometric and spectroscopic observations of the interacting transient SN 2009ip taken during the 2013 and 2014 observing seasons. We characterize the photometric evolution as a steady and smooth decline in all bands, with a decline rate that is slower than expected for a solely Co-56-powered supernova at late phases. No further outbursts or eruptions were seen over a two year period from 2012 December until 2014 December. SN 2009ip remains brighter than its historic minimum from pre-discovery images. Spectroscopically, SN 2009ip continues to be dominated by strong, narrow (less than or similar to 2000 km s(-1)) emission lines of H, He, Ca, and Fe. While we make tenuous detections of [Fe II] lambda 7155 and [O I] lambda lambda 6300, 6364 lines at the end of 2013 June and the start of 2013 October, respectively, we see no strong broad nebular emission lines that could point to a core-collapse origin. In general, the lines appear relatively symmetric, with the exception of our final spectrum in 2014 May, when we observe the appearance of a redshifted shoulder of emission at +550 km s(-1). The lines are not blueshifted, and we see no significant near-or mid-infrared excess. From the spectroscopic and photometric evolution of SN 2009ip until 820 d after the start of the 2012a event, we still see no conclusive evidence for core-collapse, although whether any such signs could be masked by ongoing interaction is unclear.

PTF11iqb: cool supergiant mass-loss that bridges the gap between Type IIn and normal supernovae

PTF11iqb was initially classified as a TypeIIn event caught very early after explosion. It showed narrow Wolf-Rayet (WR) spectral features on day 2, but the narrow emission weakened quickly and the spectrum morphed to resemble those of Types II-L and II-P. At late times, Halpha emission exhibited a complex, multipeaked profile reminiscent of SN1998S. In terms of spectroscopic evolution, we find that PTF11iqb was a near twin of SN~1998S, although with weaker interaction with circumstellar material (CSM) at early times, and stronger CSM interaction at late times. We interpret the spectral changes as caused by early interaction with asymmetric CSM that is quickly (by day 20) enveloped by the expanding SN ejecta photosphere, but then revealed again after the end of the plateau when the photosphere recedes. The light curve can be matched with a simple model for weak CSM interaction added to the light curve of a normal SN~II-P. This plateau requires that the progenitor had an extended H envelope like a red supergiant, consistent with the slow progenitor wind speed indicated by narrow emission. The cool supergiant progenitor is significant because PTF11iqb showed WR features in its early spectrum --- meaning that the presence of such WR features in an early SN spectrum does not necessarily indicate a WR-like progenitor. [abridged] Overall, PTF11iqb bridges SNe~IIn with weaker pre-SN mass loss seen in SNe II-L and II-P, implying a continuum between these types.

Electron-Impact Uncertainty Analysis and its Impact on Certain Temperature Diagnostics

In this study we calculate the electron-impact uncertainties in atomic data for direct ionization and recombination and investigate the role of these uncertainties on spectral diagnostics. We outline a systematic approach to assigning meaningful uncertainties that vary with electron temperature. Once these uncertainty parameters have been evaluated, we can then calculate the uncertainties on key diagnostics through a Monte Carlo routine, using the Astrophysical Emission Code (APEC) [Smith et al. 2001]. We incorporate these uncertainties into well known temperature diagnostics, such as the Lyman alpha versus resonance line ratio and the G ratio. We compare these calculations to a study performed by [Testa et al. 2004], where significant discrepancies in the two diagnostic ratios were observed. We conclude that while the atomic physics uncertainties play a noticeable role in the discrepancies observed by Testa, they do not explain all of them. This indicates that there is another physical process occurring in the system that is not being taken into account. This work is supported in part by the National Science Foundation REU and Department of Defense ASSURE programs under NSF Grant no. 1262851 and by the Smithsonian Institution.

Valence and L-shell photoionization of Cl-like argon using R-matrix techniques

Photoionization cross-sections are obtained using the relativistic DiracAtomic R-matrix Codes (DARC) for all valence and L-shell energy ranges between 27 and 270 eV. A total of 557 levels arising from the dominant configurations 3s²3p⁴, 3s³p5, 3p⁶, 3s²3p³[3d, 4s, 4p], 3p⁵3d, 3s²3p²3d², 3s³p4³d, 3s³p3³d2 and 2s²2p⁵3s²3p⁵ have been included in the targetwavefunction representation of the Ar III ion, including up to 4p in the orbital basis. We also performed a smaller Breit-Pauli (BP) calculation containing the lowest 124 levels. Direct comparisons are made with previous theoretical and experimental work for both valence shell and L-shell photoionization. Excellent agreement was found for transitions involving the ²P^o initial state to all allowed final states for both calculations across a range of photon energies. A number of resonant states have been identified to help analyse and explain the nature of the spectra at photon energies between 250 and 270 eV.

Ultraviolet emission lines of Si II in cool star and solar spectra

Recent atomic physics calculations for Si II are employed within the CLOUDY modelling code to analyse Hubble Space Telescope (HST) STIS ultraviolet spectra of three cool stars, β Geminorum, α Centauri A and B, as well as previously published HST/GHRS observations of α Tau, plus solar quiet Sun data from the High Resolution Telescope and Spectrograph. Discrepancies found previously between theory and observation for line intensity ratios involving the 3s²3p ²PJ-3s3p^{2 4}P^J' intercombination multiplet of Si II at ~ 2335 Å are significantly reduced, as are those for ratios containing the 3s²3p ²PJ-3s3p^{2 2}DJ ~ transitions at ~1816 Å. This is primarily due to the effect of the new Si II transition probabilities. However, these atomic data are not only very different from previous calculations, but also show large disagreements with measurements, specifically those of Calamai et al. for the intercombination lines. New measurements of transition probabilities for Si II are hence urgently required to confirm (or otherwise) the accuracy of the recently calculated values. If the new calculations are confirmed, then a long-standing discrepancy between theory and observation will have finally been resolved. However, if the older measurements are found to be correct, then the agreement between theory and observation is simply a coincidence and the existing discrepancies remain.

Effective collision strengths for forbidden transitions among the 3s23p3 fine-structure levels of Cl III

Effective collision strengths for the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^o, ²D^o and ²P^o levels in the 3s²3p³ configuration of Cl III are presented. The calculation employs the multichannel R-matrix method to compute the electron-impact excitation collision strengths in a close-coupling expansion, which incorporates the lowest 23 LS target eigenstates of Cl III. These states are formed from the 3s²3p³, 3s3p⁴, 3s²3p²3d and 3s²3p²4s configurations. The Maxwellian-averaged effective collision strengths are presented graphically for all 10 fine-structure transitions over a wide range of electron temperatures appropriate for astrophysical applications [log T(K) = 3.3 - log T(K) = 5.9]. Comparisons are made with the earlier seven-state close-coupling calculation of Butler & Zeippen, and in general excellent agreement is found in the low-temperature region where a comparison is possible [log T(K) = 3.3 - log T(K) = 4.7]. However, discrepancies of up to 30 per cent are found to occur for the forbidden transitions which involve the ⁴S^o ground state level, particularly for the lowest temperatures considered. At the higher temperatures, the present data are the only reliable results currently available.

Effective collision strengths for fine-structure forbidden transitions among the 2s22p3 levels of S x

Effective collision strengths for electron-impact excitation of the N-like ion S x are calculated in the close-coupling approximation using the multichannel R-matrix method. Specific attention is given to the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^O, ²D^o and ²P^o levels in the 2s²2p³ ground configuration. The total (e^- + ion) wavefunction is expanded in terms of the 11 lowest LS eigenstates of S x, and each eigenstate is represented by extensive configuration-interaction wavefunctions. The collision strengths obtained are thermally averaged over a Maxwellian distribution of velocities, for all 10 fine-structure transitions, over the range of electron temperatures log T(K) = 4.6-6.7 (the range appropriate for astrophysical applications). The present effective collision strengths are the only results currently available for these fine-structure transition rates.

Effective collision strengths for fine-structure forbidden transitions among the 2s22p3 levels of Ne IV

Effective collision strengths for electron-impact excitation of the N-like ion NeIV are calculated in the close-coupling approximation using the multichannel R-matrix method. Specific attention is given to the 10 astrophysically important fine-structure forbidden transitions among the ⁴S^o, ²D^o and ²P^o levels in the 2s²2p³ ground-state configuration. The expansion of the total wavefunction incorporates the lowest 11 LS eigenstates of NeIV, consisting of eight n = 2 terms with configurations 2s²2p³, 2s2p⁴ and 2p⁵, together with three n = 3 states of configuration 2s²2p²3s. We present in graphical form the effective collision strengths obtained by thermally averaging the collision strengths over a Maxwellian distribution of velocities, for all 10 fine-structure transitions, over the range of electron temperatures log T(K) = 3.6 to log T(K) = 6.1 (the range appropriate for astrophysical applications). Comparisons are made with the earlier, less sophisticated close-coupling calculation of Giles, and excellent agreement is found in the limited temperature region where a comparison is possible [log T(K) = 3.7 to log 7(K) = 4.3]. At higher temperatures the present data are the only reliable results currently available.

Effective collision strengths for fine-structure forbidden transitions among the 3s 23p 3 levels of Ar IV

The multichannel R-matrix method is used to compute electron impact excitation collision strengths in Ar IV for all fine-structure transitions among the ⁴S°, ²D° and ²P° levels in the 3s ²3p ³ ground configuration. Included in the expansion of the total wavefunction are the lowest 13 LS target eigenstates of Ar iv formed from the 3s ²3p ³, 3s3p ⁴ and 3s ²3p ²3d configurations. The effective collision strengths, obtained by averaging the electron collision strengths over a Maxwellian distribution of electron velocities, are presented for all 10 fine-structure transitions over a wide range of electron temperatures of astrophysical interest (T _e = 2000-100 000 K). Comparisons are made with an earlier 7-state close-coupling calculation by Zeippen, Butler & Le Bourlot, and significant differences are found to occur for many of the forbidden transitions considered, in particular those involving the ⁴S° ground state, where discrepancies of up to a factor of 3 are found in the low-temperature region. © 1997 RAS.