Context. A currently unsolved question in supernova (SN) research is the origin of stripped-envelope supernovae (SESNe). Such SNe lack spectral signatures of hydrogen (Type Ib), or hydrogen and helium (Type Ic), indicating that the outer stellar layers have been stripped during their evolution. The mechanism for this is not well understood, and to disentangle the different scenarios’ determination of nucleosynthesis yields from observed spectra can be attempted. However, the interpretation of observations depends on the adopted spectral models. A previously missing ingredient in these is the inclusion of molecular effects, which can be significant.

Aims. We aim to investigate how the molecular chemistry in SESNe affect physical conditions and optical spectra, and produce ro-vibrational emission in the mid-infrared (MIR). We also aim to assess the diagnostic potential of observations of such MIR emission with JWST.

Methods. We coupled a chemical kinetic network including carbon, oxygen, silicon, and sulfur-bearing molecules into the nonlocal thermal equilibrium (NLTE) spectral synthesis code SUMO. We let four species – CO, SiO, SiS, and SO – participate in NLTE cooling of the gas to achieve self-consistency between the molecule formation and the temperature. We applied the new framework to model the spectrum of a Type Ic SN in the 100–600 days time range.

Results. Molecules are predicted to form in SESN ejecta in significant quantities (typical mass 10⁻³ M_⊙) throughout the 100–600 days interval. The impact on the temperature and optical emission depends on the density of the oxygen zones and varies with epoch. For example, the [O I] 6300, 6364 feature can be quenched by molecules from 200 to 450 days depending on density. The MIR predictions show strong emission in the fundamental bands of CO, SiO, and SiS, and in the CO and SiO overtones.

Conclusions. Type Ibc SN ejecta have a rich chemistry and considering the effect of molecules is important for modeling the temperature and atomic emission in the nebular phase. Observations of SESNe with JWST hold promise to provide the first detections of SiS and SO, and to give information on zone masses and densities of the ejecta. Combined optical, near-infrared, and MIR observations can break degeneracies and achieve a more complete picture of the nucleosynthesis, chemistry, and origin of Type Ibc SNe.

Iron is one of the most important elements in-stellar astrophysics. However, spectroscopic measurements of its abundance are prone to systematic modelling errors. We present three dimensional non-local thermodynamic equilibrium (3D non-LTE) calculations across 32 STAGGER-grid models with effective temperatures from 5000 K to 6500 K, surface gravities of 4.0 dex and 4.5 dex, and metallicities from −3 dex to 0 dex, and we study the effects on 171 Fe I and 12 Fe II optical lines. In warm metal-poor stars, the 3D non-LTE abundances are up to 0.5 dex larger than 1D LTE abundances inferred from Fe I lines of an intermediate excitation potential. In contrast, the 3D non-LTE abundances can be 0.2 dex smaller in cool metal-poor stars when using Fe I lines of a low excitation potential. The corresponding abundance differences between 3D non-LTE and 1D non-LTE are generally less severe but can still reach ±0.2 dex. For Fe II lines, the 3D abundances range from up to 0.15 dex larger to 0.10 dex smaller than 1D abundances, with negligible departures from 3D LTE except for the warmest stars at the lowest metallicities. The results were used to correct 1D LTE abundances of the Sun and Procyon (HD 61421), and of the metal-poor stars HD 84937 and HD 140283, using an interpolation routine based on neural networks. The 3D non-LTE models achieve an improved ionisation balance in all four stars. In the two metal-poor stars, they removed excitation imbalances amounting to 250 K to 300 K errors in effective temperature. For Procyon, the 3D non-LTE models suggest [Fe/H] = 0.11 ± 0.03, which is significantly larger than literature values based on simpler models. We make the 3D non-LTE interpolation routine for FG-type dwarfs publicly available, in addition to 1D non-LTE departure coefficients for standard MARCS models of FGKM-type dwarfs and giants. These tools, together with an extended 3D LTE grid for Fe II from 2019, can help improve the accuracy of stellar parameter and iron abundance determinations for late-type stars.

Context. The inclusion of molecular physics is an important piece that tends to be missing from the puzzle when modeling the spectra of supernovae (SNe). Molecules have both a direct impact on the spectra, particularly in the infrared, and an indirect one as a result of their influence on certain physical conditions, such as temperature.

Aims. In this paper, we aim to investigate molecular formation and non-local thermodynamic equilibrium (NLTE) cooling, with a particular focus on CO, the most commonly detected molecule in supernovae. We also aim to determine the dependency of supernova chemistry on physical parameters and the relative sensitivity to rate uncertainties.

Methods. We implemented a chemical kinetic description of the destruction and formation of molecules into the SN spectral synthesis code SUMO. In addition, selected molecules were coupled into the full NLTE level population framework and, thus, we incorporated molecular NLTE cooling into the temperature equation. We produced a test model of the CO formation in SN 1987A between 150 and 600 days and investigated the sensitivity of the resulting molecular masses to the input parameters.

Results. We find that there is a close inter-dependency between the thermal evolution and the amount of CO formed, mainly through an important temperature-sensitive CO destruction process with O+. After a few hundred days, CO completely dominates the cooling of the oxygen-carbon zone of the supernova which, therefore, contributes little optical emission. The uncertainty of the calculated CO mass scales approximately linearly with the typical uncertainty factor for individual rates. We demonstrate how molecular masses can potentially be used to constrain various physical parameters of the supernova.

Context. The stellar winds of asymptotic giant branch (AGB) stars are commonly attributed to radiation pressure on dust grains, formed in the wake of shock waves that arise in the stellar atmospheres. The mass loss due to these outflows is substantial, and modelling the dynamical properties of the winds is essential both for studies of individual stars and for understanding the evolution of stellar populations with low to intermediate mass.

Aims. The purpose of this work is to present an extensive grid of dynamical atmosphere and wind models for M-type AGB stars, covering a wide range of relevant stellar parameters.

Methods. We used the DARWIN code, which includes frequency-dependent radiation-hydrodynamics and a time-dependent description of dust condensation and evaporation, to simulate the dynamical atmosphere. The wind-driving mechanism is photon scattering on submicron-sized Mg2SiO4 grains. The grid consists of similar to 4000 models, with luminosities from L-* = 890 L-circle dot to L-* = 40 000 L-circle dot and effective temperatures from 2200 to 3400 K. For the first time different current stellar masses are explored with M-type DARWIN models, ranging from 0.75 M-circle dot to 3 M-circle dot. The modelling results are radial atmospheric structures, dynamical properties such as mass-loss rates and wind velocities, and dust properties (e.g. grain sizes, dust-to-gas ratios, and degree of condensed Si).

Results. We find that the mass-loss rates of the models correlate strongly with luminosity. They also correlate with the ratio L-*/M-* : increasing L-*/M-* by an order of magnitude increases the mass-loss rates by about three orders of magnitude, which may naturally create a superwind regime in evolution models. There is, however, no discernible trend of mass-loss rate with effective temperature, in contrast to what is found for C-type AGB stars. We also find that the mass-loss rates level off at luminosities higher than similar to 14 000 L-circle dot, and consequently at pulsation periods longer than similar to 800 days. The final grain radii range from 0.25 to 0.6 mu m. The amount of condensed Si is typically between 10 and 40%, with gas-to-dust mass ratios between 500 and 4000.