Monte-Carlo methods for NLTE spectral synthesis of supernovae

We present JEKYLL, a new code for modelling of supernova (SN) spectra and lightcurves based on Monte-Carlo (MC) techniques for the radiative transfer. The code assumes spherical symmetry, homologous expansion and steady state for the matter, but is otherwise capable of solving the time-dependent radiative transfer problem in non-local-thermodynamic-equilibrium (NLTE). The method used was introduced in a series of papers by Lucy, but the full time-dependent NLTE capabilities of it have never been tested. Here, we have extended the method to include non-thermal excitation and ionization as well as charge-transfer and two-photon processes. Based on earlier work, the non-thermal rates are calculated by solving the Spencer-Fano equation. Using a method previously developed for the SUMO code, macroscopic mixing of the material is taken into account in a statistical sense. In addition, a statistical Markov-chain model is used to sample the emission frequency, and we introduce a method to control the sampling of the radiation field. Except for a description of JEKYLL, we provide comparisons with the ARTIS, SUMO and CMFGEN codes, which show good agreement in the calculated spectra as well as the state of the gas. In particular, the comparison with CMFGEN, which is similar in terms of physics but uses a different technique, shows that the Lucy method does indeed converge in the time-dependent NLTE case. Finally, as an example of the time-dependent NLTE capabilities of JEKYLL, we present a model of a Type IIb SN, taken from a set of models presented and discussed in detail in an accompanying paper. Based on this model we investigate the effects of NLTE, in particular those arising from non-thermal excitation and ionization, and find strong effects even on the bolometric lightcurve. This highlights the need for full NLTE calculations when simulating the spectra and lightcurves of SNe.Comment: Accepted for publication by Astronomy & Astrophysic

Astrophysics in 2006

The fastest pulsar and the slowest nova; the oldest galaxies and the youngest stars; the weirdest life forms and the commonest dwarfs; the highest energy particles and the lowest energy photons. These were some of the extremes of Astrophysics 2006. We attempt also to bring you updates on things of which there is currently only one (habitable planets, the Sun, and the universe) and others of which there are always many, like meteors and molecules, black holes and binaries.Comment: 244 pages, no figure

Spin transport properties of a quantum dot coupled to ferromagnetic leads with noncollinear magnetizations

A correct general formula for the spin current through an interacting quantum dot coupled to ferromagnetic leads with magnetization at an arbitrary angle $\theta$ is derived within the framework of the Keldysh formalism. Under asymmetric conditions, the spin current component J_{z} may change sign for $0<\theta<\pi$. It is shown that the spin current and spin tunneling magnetoresistance exhibit different angle dependence in the free and Coulomb blockade regimes. In the latter case, the competition of spin precession and the spin-valve effect could lead to an anomaly in the angle dependence of the spin current.Comment: 7 pages, 4 figures; some parts of the text has been revised in this version accepted by J. Phys.: Condens. Matte

Late-time spectral line formation in Type IIb supernovae, with application to SN 1993J, SN 2008ax, and SN 2011dh

We investigate line formation processes in Type IIb supernovae (SNe) from 100 to 500 days post-explosion using spectral synthesis calculations. The modeling identifies the nuclear burning layers and physical mechanisms that produce the major emission lines, and the diagnostic potential of these. We compare the model calculations with data on the three best observed Type IIb SNe to-date - SN 1993J, SN 2008ax, and SN 2011dh. Oxygen nucleosynthesis depends sensitively on the main-sequence mass of the star and modeling of the [O I] 6300, 6364 lines constrains the progenitors of these three SNe to the M_ZAMS=12-16 M_sun range (ejected oxygen masses 0.3-0.9 M_sun), with SN 2011dh towards the lower end and SN 1993J towards the upper end of the range. The high ejecta masses from M_ZAMS >= 17 M_sun progenitors give rise to brighter nebular phase emission lines than observed. Nucleosynthesis analysis thus supports a scenario of low/moderate mass progenitors for Type IIb SNe, and by implication an origin in binary systems. We demonstrate how oxygen and magnesium recombination lines may be combined to diagnose the magnesium mass in the SN ejecta. For SN 2011dh, a magnesium mass of of 0.02-0.14 M_sun is derived, which gives a Mg/O production ratio consistent with the solar value. Nitrogen left in the He envelope from CNO-burning gives strong [N II] 6548, 6583 emission lines that dominate over H-alpha emission in our models. The hydrogen envelopes of Type IIb SNe are too small and dilute to produce any noticeable H-alpha emission or absorption after ~150 days, and nebular phase emission seen around 6550 A is in many cases likely caused by [N II] 6548, 6583. Finally, the influence of radiative transport on the emergent line profiles is investigated...(abridged)Comment: Published versio

The nebular spectra of SN 2012aw and constraints on stellar nucleosynthesis from oxygen emission lines

We present nebular phase optical and near-infrared spectroscopy of the Type IIP supernova SN 2012aw combined with NLTE radiative transfer calculations applied to ejecta from stellar evolution/explosion models. Our spectral synthesis models generally show good agreement with the ejecta from a MZAMS = 15 Msun progenitor star. The emission lines of oxygen, sodium, and magnesium are all consistent with the nucleosynthesis in a progenitor in the 14 - 18 Msun range. We also demonstrate how the evolution of the oxygen cooling lines of [O I] 5577 A, [O I] 6300 A, and [O I] 6364 A can be used to constrain the mass of oxygen in the non-molecularly cooled ashes to < 1 Msun, independent of the mixing in the ejecta. This constraint implies that any progenitor model of initial mass greater than 20 Msun would be difficult to reconcile with the observed line strengths. A stellar progenitor of around MZAMS = 15 Msun can consistently explain the directly measured luminosity of the progenitor star, the observed nebular spectra, and the inferred pre-supernova mass-loss rate. We conclude that there is still no convincing example of a Type IIP explosion showing the nucleosynthesis expected from a MZAMS > 20 Msun progenitor.Comment: Accepted for publication in MNRA