Search

Context. Type Ia supernovae (SN Ia) have become an invaluable cosmological tool because their exceptional brightness makes them observable even at very large distances (up to redshifts around z ≈ 1). To investigate possible systematic differences between local and distant SN Ia requires detailed models whose synthetic spectra can be compared to observations and in which the solution of the radiative transfer is a key ingredient. One commonly employed method is the Feautrier scheme, which is generally very robust but can lead to wrong results under certain conditions that frequently occur when modeling supernova ejecta or even the radiatively driven expanding atmospheres of hot stars. Aims. We attempt to improve the procedure we have developed for simulating the radiative transfer of metal-rich, intermediate- and low-density, line-dominated atmospheres to allow the method to be applied successfully even under conditions of high expansion velocities. Methods. We use a sophisticated model atmosphere code that considers the non-LTE effects and large velocity gradients that strongly affect the physics of SN Ia atmospheres at all wavelengths to simulate the formation of SN Ia spectra by the thousands of strong spectral lines that intricately interact with the "pseudo-continuum" formed entirely by these Doppler-shifted lines themselves. We focus on investigating the behavior of the Feautrier scheme under these conditions. Results. Synthetic spectra of SN Ia, a complex product of computer models replicating numerous physical processes that determine the conditions of matter and radiation in the ejecta, are affected by large spatial jumps of the line-dominated opacities and source functions for which the application of even well established methods may harbor certain pitfalls. We analyse the conditions that can lead to a breakdown of conventional procedures and we derive a modified description that yields more accurate results in the given circumstances for the Feautrier radiative transfer solver. [ABSTRACT FROM AUTHOR]

Context. In type Ia supernova (SN Ia) envelopes a huge number of lines of different elements overlap within their thermal Doppler widths, and this problem is exacerbated by the circumstance that up to 20% of these lines can have a line optical depth higher than 1. The stagnation of the lambda iteration in such optically thick cases is one of the fundamental physical problems inherent in the iterative solution of the non-LTE problem, and the failure of a lambda iteration to converge is a point of crucial importance whose physical significance must be understood completely. Aims. We discuss a general problem related to radiative transfer under the physical conditions of supernova ejecta that involves a failure of the usual non-LTE iteration scheme to converge when multiple strong opacities belonging to different physical transitions come together, similar to the well-known situation where convergence is impaired even when only a single process attains high optical depths. The convergence problem is independent of the chosen frequency and depth grid spacing, independent of whether the radiative transfer is solved in the comoving or observer's frame, and independent of whether a common complete-linearization scheme or a conventional accelerated lambda iteration (ALI) is used. The problem appears when all millions of line transitions required for a realistic description of SN Ia envelopes are treated in the frame of a comprehensive non-LTE model. The only solution to this problem is a complete-linearization approach that considers all ions of all elements simultaneously, or an adequate generalization of the established ALI technique that accounts for the mutual interaction of the strong spectral lines of different elements and which thereby unfreezes the "stuck" state of the iteration. Methods. The physics of the atmospheres of SN Ia are strongly affected by the high-velocity expansion of the ejecta, which dominates the formation of the spectra at all wavelength ranges. Thus, hydrodynamic explosion models and realistic model atmospheres that take into account the strong deviation from local thermodynamic equilibrium (LTE) are necessary for the synthesis and analysis of the spectra. In this regard one of the biggest challenges we have found in modeling the radiative transfer in SN Ia is the fact that the radiative energy in the UV has to be transferred only via spectral lines into the optical regime to be able to leave the ejecta. However, convergence of the model toward a state where this is possible is impaired when using the standard procedures. We report on improvements in our approach of computing synthetic spectra for SN Ia with respect to (i) an improved and sophisticated treatment of many 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 arising from the radioactive decay of mostly 56Ni and 56Co. Results. We show that an ALI procedure we have developed for the mutual interaction of strong spectral lines appearing in the atmospheres of SNe Ia solves the long-standing problem of transferring the radiative energy from the UV into the optical regime. Our new method thus constitutes a foundation for more refined models, such as those including energy deposition. In this regard we furthermore show synthetic spectra obtained with various methods adopted for the released energy and compare them with observations. We discuss in detail applications of the diagnostic technique by example of a standard type Ia supernova, 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 emergent 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 knowledge provides an important insight into the process of extracting information from observed SN Ia spectra, since these spectra are a complex product of numerous unobservable SN Ia spectral features, which are thus analyzed in parallel to the observable SN Ia spectral features. [ABSTRACT FROM AUTHOR]

Context. Rotational mixing is known to significantly affect the evolution of massive stars; however, we still lack a consensus regarding the various possible modeling approaches and mixing recipes describing this process. The empirical investigation of surface abundances of carbon, nitrogen, and oxygen (CNO) in large samples of O- and B-type stars will be essential for providing meaningful observational constraints on the different available stellar evolution models. Aims. Setting up and testing adequate tools to perform CNO surface abundance determinations for large samples of O-type stars, by means of the fast performance, NLTE, unified model atmosphere code FASTWIND. Methods. We have developed a set of semi-automatic tools for measuring and analyzing the observed equivalent widths of strategic optical C, N, and O lines from different ions. Our analysis strategy is based on a χ2 minimization of weighted differences between observed and synthetic equivalent widths, the latter computed from tailored model grids. We have paid special attention to the (significant) errors introduced by typical uncertainties in stellar parameters. In this pilot study, we describe these tools, and test their performance and reliability using a set of high quality spectra of a sample of 18 presumably single Galactic O-type stars with low projected rotational velocities (v sin i≲100 km s−1), and previously determined stellar parameters. In particular, we have compared the outcome of our analysis with results from existing studies and theoretical stellar evolution models. Results. Most of our results for carbon and nitrogen agree, within the errors, with both theoretical expectations and literature values. While many cooler dwarfs display C and N abundances close to solar, some of the early- and mid-O dwarfs – and most supergiants – show significant enrichment in N and depletion in C. Our results for oxygen in late-O dwarfs are, however, unexpectedly low, possibly indicating deficiencies in the adopted oxygen model atom. For all other objects, no systematic problems in their oxygen content have been identified. Specific stars in our sample show peculiarities in their abundances, and we suggest hypotheses regarding their origin. Conclusions. Our method is (almost) ready to be applied to large samples of late and mid O-type stars – although the oxygen model atom needs to be improved and carefully tested first. For early O-type stars (O4 and hotter), a simultaneous UV analysis seems to be inevitable, due to the scarcity and weakness of optical C and O lines. This will necessarily imply a more complex modeling, additionally accounting for the effects of X-rays from wind-embedded shocks and wind inhomogeneities. [ABSTRACT FROM AUTHOR]

Context. The diffuse ionized gas (DIG) constitutes the largest fraction of the total ionized interstellar matter in star-forming galaxies, but it is still unclear whether the ionization is driven predominantly by the ionizing radiation of hot massive stars, as in H II regions, or whether additional sources of ionization have to be considered. Key to understanding the ionization mechanisms in the DIG is the line emission by the ionized gas. Aims. We systematically explore a plausible subset of the parameter space involving effective temperatures and metallicities of the ionizing sources, the effects of the hardening of their radiation by surrounding "leaky" H II 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 computed 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 interstellar gas we used spherically symmetric photoionization models and state-of-the-art three-dimensional (3D) non-LTE radiative transfer simulations, considering hydrogen, helium, and the most abundant metals. We first applied these methods to classical H II regions around hot stars, using the model SEDs at different metallicities and effective temperatures as ionizing fluxes, and compute the SEDs of the escaping radiation for different escape fractions of hydrogen-ionizing photons. In a second step, we studied the effects of the escaping radiation on the more dilute and extended DIG. Using 3D models simulating a section of a galactic spiral arm, we computed the ionization structure in the DIG for different scenarios for the inhomogeneity of the gas, assuming ionization by a stellar population SED based on plausible parameters. Results. We provide quantitative predictions of how the line ratios from H II regions and the DIG vary as a function of metallicity Z, stellar effective temperature Teff, and escape fraction fesc from the H II region. The range of predicted line ratios reinforces the hypothesis that the DIG is ionized by (filtered) radiation from hot stars. At one-tenth solar metallicity, radiation hardening is mostly due to hydrogen and helium, whereas at solar metallicity absorption by metals plays a significant role. The effects of hardening are seen primarily in the increase in the emission line ratios of the most important cooling lines of the gas, [N II]∕Hβ and [O II]∕Hβ at lower Teff, and [O III]∕Hβ at higher Teff. For low Teff nearly the entire He I-ionizing radiation is absorbed in the H II regions, thereby preventing the formation of high ionization stages such as O III in the DIG. The ionization structure of the DIG depends strongly on both the clumping factor fcl = 〈nH2〉/〈nH2〉 ${f_{\textrm{cl}}}\,{=}\,\langle{n_{\textrm{H}}}^2\rangle/\langle{n_{\textrm{H}}}\rangle^2$ f cl   =   〈 n H 2 〉 / 〈 n H 〉 2 and the large-scale distribution of the gas. In our simulations about 10% of the ionizing radiation produced by hot massive stars in a spiral arm is sufficient to ionize the DIG up to a height of approximately 1 kpc above the galactic plane for a clumping factor close to the observed value of fcl ~ 5. Even small changes in simulation parameters such as 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. 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. [ABSTRACT FROM AUTHOR]

Context. The determination of chemical abundances constitutes a fundamental requirement for obtaining a complete picture of a star. Particularly in massive stars, CNO abundances are of prime interest, due to the nuclear CNO-cycle, and various mixing processes which bring these elements to the surface. The precise determination of carbon abundances, together with N and O, is thus a key ingredient for understanding the different phases of stellar evolution. Aims. We aim to enable a reliable carbon spectroscopy for our unified non-LTE atmosphere code FASTWIND. Methods. We have developed a new carbon model atom including C II/III/IV/V, and we discuss specific problems related to carbon spectroscopy in O-type stars. We describe different tests we have performed to examine the reliability of our implementation, and investigate which mechanisms influence the carbon ionization balance. By comparing with high-resolution spectra from six O-type stars, we verified to what extent observational constraints can be reproduced by our new carbon line synthesis. Results. Carbon lines are even more sensitive to a variation of Teff, log g, and Ṁ, than hydrogen and helium lines. We are able to reproduce most of the observed lines from our stellar sample, and to estimate those specific carbon abundances which bring the lines from different ions into agreement (three stages in parallel for cool objects, two for intermediate O-types). For hot dwarfs and supergiants earlier than O7, X-rays from wind-embedded shocks can have an impact on the synthesized line strengths, particularly for C IV, potentially affecting the abundance determination. Dielectronic recombination has a significant impact on the ionization balance in the wind. Conclusions. We demonstrate our capability to derive realistic carbon abundances by means of FASTWIND, using our recently developed model atom. We find that complex effects can have a strong influence on the carbon ionization balance in hot stars. For a further understanding, the UV range needs to be explored as well. By means of detailed and available nitrogen and oxygen model atoms, we will be able to perform a complete CNO abundance analysis for larger samples of massive stars, and to provide constraints on corresponding evolutionary models and aspects. [ABSTRACT FROM AUTHOR]