Model atmospheres for type Ia supernovae: Basic steps towards realistic synthetic spectra

Type Ia supernovae are an important tool for studying the expansion history of the universe. Advancing our yet incomplete understanding of the explosion scenario requires detailed and realistic numerical models in order to interpret and analyze the growing amount of observational data. Here we present first results of our new NLTE model calculations for the expanding atmospheres of type Ia supernovae that employ a detailed and consistent treatment of all important NLTE effects as well as line blocking and blanketing. The comparison of the synthetic spectra resulting from these models with observed data shows that the employed methods represent an important step towards a more realistic description of the atmospheres of supernovae Ia.Comment: 4 pages, 1 figure, to appear in: Proceedings of the 11th Workshop on Nuclear Astrophysics, Ringberg Castle, Germany, 200

Numerical Models for the Diffuse Ionized Gas in Galaxies. II. Three-dimensional radiative transfer in inhomogeneous interstellar structures as a tool for analyzing the diffuse ionized gas

Aims: We systematically explore a plausible subset of the parameter space involving effective temperatures and metallicities of the ionizing stellar sources, the effects of the hardening of their radiation by surrounding leaky HII regions with different escape fractions, as well as different scenarios for the clumpiness of the DIG, and compute the resulting line strength ratios for a number of diagnostic optical emission lines. Methods: For the ionizing fluxes we compute a grid of stellar spectral energy distributions (SEDs) from detailed, fully non-LTE model atmospheres that include the effects of stellar winds and line blocking and blanketing. To calculate the ionization and temperature structure in the HII regions and the diffuse ionized gas we use spherically symmetric photoionization models as well as state-of-the-art three-dimensional (3D) non-LTE radiative transfer simulations, considering hydrogen, helium, and the most abundant metals. Results: We provide quantitative predictions of how the line ratios from HII regions and the DIG vary as a function of metallicity, stellar effective temperature, and escape fraction from the HII region. The range of predicted line ratios reinforces the hypothesis that the DIG is ionized by (filtered) radiation from hot stars; however, comparison of observed and predicted line ratios indicates that the DIG is typically ionized with a softer SED than predicted by the chosen stellar population synthesis model. Even small changes in simulation parameters like the clumping factor can lead to considerable variation in the ionized volume. Both for a more homogeneous gas and a very inhomogeneous gas containing both dense clumps and channels with low gas density, the ionized region in the dilute gas above the galactic plane can cease to be radiation-bounded, allowing the ionizing radiation to leak into the intergalactic medium.Comment: 21 pages, 9 figures, accepted by A&

Radiation-driven winds of hot luminous stars. XVI. Expanding atmospheres of massive and very massive stars and the evolution of dense stellar clusters

Context: Starbursts, and particularly their high-mass stars, play an essential role in the evolution of galaxies. The winds of massive stars not only significantly influence their surroundings, but the mass loss also profoundly affects the evolution of the stars themselves. In addition to the evolution of each star, the evolution of the dense cores of massive starburst clusters is affected by N-body interactions, and the formation of very massive stars via mergers may be decisive for the evolution of the cluster. Aims: To introduce an advanced diagnostic method of O-type stellar atmospheres with winds, including an assessment of the accuracy of the determinations of abundances, stellar and wind parameters. Methods: We combine consistent models of expanding atmospheres with detailed stellar evolutionary calculations of massive and very massive single stars with regard to the evolution of dense stellar clusters. Accurate predictions of the mass loss rates of very massive stars requires a highly consistent treatment of the statistical equilibrium and the hydrodynamic and radiative processes in the expanding atmospheres. Results: We present computed mass loss rates, terminal wind velocities, and spectral energy distributions of massive and very massive stars of different metallicities, calculated from atmospheric models with an improved level of consistency. Conclusions: Stellar evolutionary calculations using our computed mass loss rates show that low-metallicity very massive stars lose only a very small amount of their mass, making it unlikely that very massive population III stars cause a significant helium enrichment of the interstellar medium. Solar-metallicity stars have higher mass-loss rates, but these are not so high to exclude very massive stars formed by mergers in dense clusters from ending their life massive enough to form intermediate-mass black holes.Comment: Accepted by A&

Ionizing Photon Emission Rates from O- and Early B-type Stars and Clusters

We present new computations of the ionizing spectral energy distributions (SEDs) and Lyman continuum (Lyc) and HeI continuum photon emission rates, for hot O-type and early B-type stars. We consider solar metallicity stars, with effective temperatures ranging from 25,000 to 55,000 K and surface gravities (cm s^-2) logg ranging from 3 to 4, covering the full range of spectral types and luminosity classes for hot stars. We use our updated (WM-basic) code to construct radiation-driven wind atmosphere models for hot stars. Our models include the coupled effects of hydrodynamics and non-LTE radiative transfer in spherically outflowing winds, including the detailed effects of metal line blocking and line blanketing on the radiative transfer and energy balance. We incorporate our hot-star models into our population synthesis code (STARS), and we compute the time-dependent SEDs and resulting Lyc and HeI emission rates for evolving star clusters. We present results for continuous and impulsive star formation for a range of assumed stellar initial mass functions.Comment: 23 pages, 7 figures. To appear in the Astrophysical Journal. For grid of star models see ftp://wise3.tau.ac.il/pub/star

Non-LTE models for synthetic spectra of type Ia supernovae. III. An accelerated lambda iteration procedure for the mutual interaction of strong spectral lines in SN Ia models with and without energy deposition

Context. Spectroscopic analyses to interpret the spectra of the brightest supernovae from the UV to the near-IR provide a powerful tool with great astrophysical potential for the determination of the physical state of the ejecta, their chemical composition, and the SNe distances even at significant redshifts. Methods. We report on improvements of computing synthetic spectra for SNIa with respect to i) an improved and sophisticated treatment of thousands of strong lines that interact intricately with the "pseudo-continuum" formed entirely by Doppler- shifted spectral lines, ii) an improved and expanded atomic database, and iii) the inclusion of energy deposition within the ejecta. Results. We show that an accelerated lambda iteration procedure we have developed for the mutual interaction of strong spectral lines appearing in the atmospheres of SNeIa solves the longstanding problem of transferring the radiative energy from the UV into the optical regime. In detail we discuss applications of the diagnostic technique by example of a standard SNIa, where the comparison of calculated and observed spectra revealed that in the early phases the consideration of the energy deposition within the spectrum-forming regions of the ejecta does not qualitatively alter the shape of the spectra. Conclusions. The results of our investigation lead to an improved understanding of how the shape of the spectrum changes radically as function of depth in the ejecta, and show how different emergent spectra are formed as a result of the particular physical properties of SNe Ia ejecta and the resulting peculiarities in the radiative transfer. This provides an important insight into the process of extracting information from observed SNIa spectra, since these spectra are a complex product of numerous unobservable SNIa spectral features which are thus analyzed in parallel to the observable spectral features.Comment: 27 pages, 19 figures. Submitted to A&A, revised versio

Radiation-driven winds of hot luminous stars XVII. Parameters of selected central stars of PN from consistent optical and UV spectral analysis and the universality of the mass-luminosity relation

Context: The commonly accepted mass-luminosity relation of central stars of planetary nebulae (CSPNs) might not be universally valid. While earlier optical analyses could not derive masses and luminosities independently (instead taking them from theoretical evolutionary models) hydrodynamically consistent modelling of the stellar winds allows using fits to the UV spectra to consistently determine also stellar radii, masses, and luminosities without assuming a mass-luminosity relation. Recent application to a sample of CSPNs raised questions regarding the validity of the theoretical mass-luminosity relation of CSPNs. Aims: The results of the earlier UV analysis are reassessed by means of a simultaneous comparison of observed optical and UV spectra with corresponding synthetic spectra. Methods: Using published stellar parameters (a) from a consistent UV analysis and (b) from fits to optical H and He lines, we calculate simultaneous optical and UV spectra with our model atmosphere code, which has been improved by implementing Stark broadening for H and He lines. Results: Spectra computed with the parameter sets from the UV analysis yield good agreement to the observations, but spectra computed with the stellar parameters from the published optical analysis and using corresponding consistent wind parameters show large discrepancies to both the observed optical and UV spectra. The published optical analyses give good fits to the observed spectrum only because the wind parameters assumed in these analyses are inconsistent with their stellar parameters. By enforcing consistency between stellar and wind parameters, stellar parameters are obtained which disagree with the core-mass-luminosity relation for the objects analyzed. This disagreement is also evident from a completely different approach: an investigation of the dynamical wind parameters.Comment: 22 pages, 18 fugre

Radiation-driven winds of hot stars. XIII. A description of NLTE line blocking and blanketing towards realistic models of expanding atmospheres

We present significant improvements of our approach in constructing detailed atmospheric models and synthetic spectra for hot luminous stars: 1. A sophisticated and consistent description of line blocking and blanketing. Our solution concept renders the line blocking influence on the ionizing fluxes (mainly the EUV and UV are affected) in identical quality as the synthetic high resolution spectra representing the observable region. Line blanketing is properly accounted for in the energy balance. 2. A considerably improved and enhanced atomic data archive providing the basis for a detailed multilevel NLTE treatment of the metal ions (from C to Zn) and an adequate representation of line blocking and the radiative line acceleration. 3. A revised inclusion of EUV and X-ray radiation produced by cooling zones originating from shock heated matter. This new tool not only provides an easy to use method for O-star diagnostics, whereby physical constraints on the properties of stellar winds, stellar parameters, and abundances can be obtained via a comparison of observed and synthetic spectra, but also allows the astrophysically important information about the ionizing fluxes of hot stars to be determined automatically. Results illustrating this are discussed by means of a basic model grid calculated for O-stars with solar metallicity. To further demonstrate the astrophysical potential of our new method we provide a first detailed spectral diagnostic determination of the stellar parameters, the wind parameters, and the abundances by an exemplary application to the O9.5Ia supergiant alpha Cam.Comment: 42 Pages, 39 Figures. LaTeX documentclass aa, packages graphicx and times. Submitted to Astronomy and Astrophysic

Models of X-ray Photoionization in LMC X-4: Slices of a Stellar Wind

We show that the orbital variation in the UV P Cygni lines of the X-ray binary LMC X-4 results when X-rays photoionize nearly the entire region outside of the X-ray shadow of the normal star. We fit models to HST GHRS observations of N V and C IV P Cygni line profiles. Analytic methods assuming a spherically symmetric wind show that the wind velocity law is well-fit by v~(1-1/r)^beta, where beta is likely 1.4-1.6 and definitely <2.5. Escape probability models can fit the observed P Cygni profiles, and provide measurements of the stellar wind parameters. The fits determine Lx/Mdot=2.6+/-0.1 x10^43 erg/s/Msun yr, where Lx is the X-ray luminosity and Mdot is the mass-loss rate of the star. Allowing an inhomogeneous wind improves the fits. IUE spectra show greater P Cygni absorption during the second half of the orbit than during the first. We discuss possible causes of this effect.Comment: 56 pages, 12 figures, to be published in the Astrophysical Journa

On the Light Curve and Spectrum of SN 2003dh Separated from the Optical Afterglow of GRB 030329

The net optical light curves and spectra of the supernova (SN) 2003dh are obtained from the published spectra of GRB 030329, covering about 6 days before SN maximum to about 60 days after. The bulk of the U-band flux is subtracted from the observed spectra using early-time afterglow templates, because strong line blanketing greatly depresses the UV and U-band SN flux in a metal-rich, fast-moving SN atmosphere. The blue-end spectra of the gamma-ray burst (GRB)connected hypernova SN 1998bw is used to determine the amount of subtraction. The subtraction of a host galaxy template affects the late-time results. The derived SN 2003dh light curves are narrower than those of SN 1998bw, rising as fast before maximum, reaching a possibly fainter maximum, and then declining ~ 1.2-1.4 times faster. We then build UVOIR bolometric SN light curve. Allowing for uncertainties, it can be reproduced with a spherical ejecta model of Mej ~ 7+/-3 Msun, KE ~ (3.5+/-1.5)E52 ergs, with KE/Mej ~ 5 following previous spectrum modelling, and M(Ni56) ~ (0.4 +0.15/-0.1) Msun. This suggests a progenitor main-sequence mass of about 25-40 Msun, lower than SN 1998bw but significantly higher than normal Type Ic SNe and the GRB-unrelated hypernova SN 2002ap.Comment: 18 pages, 7 figures, published by Ap