The distribution of stars in the Hertzsprung Russell diagram (HRD) for a stellar conglomeration represents a snapshot of its evolving stellar population. Some of the supergiant stars may transit the HRD from blue to red and then again to blue during their late evolutionary stages, as exemplified by the progenitor of SN 1987A. Others may transit a given part of the HRD more than twice in a "blue loop" and end up as red supergiants before they explode. Since stars in blue loops spend a considerable part of their lives there, these stages may change the relative number of modeled supergiants in the HRD. Their lifetimes in turn depend upon the initial mass of the star, how convection in its interior is modeled, and how much mass loss takes place during its evolution. The observed ratio of the number of blue to red supergiants and yellow to red supergiants sensitively test the stellar evolution theory. We compare modeled number ratios of these supergiants with observed data from the Large Magellanic Cloud as it has a metallicity very similar to that of the environment of SN 2013ej. We successfully model these by taking into account moderate (exponential) convective overshooting. We explore its effect on the final radius and mass of the star prior to core collapse. The radius differs dramatically with overshoot. These factors controlling pre-supernova structure may affect the post-explosion optical/IR light curves and spectral development.

It has been well established from a variety of observations that red supergiants (RSGs) loose a lot of mass in stellar wind. Dust formed in this emitted gas over a few decades before core-collapse can lead to substantial extinction and obscure the intrinsic luminosity of the progenitor RSG. This may lead to a difficulty in determining the range of progenitor masses that lead to the different classes of supernovae. Even the nearby, well studied supernovae with pre-explosion observations, such as SN 2013ej may suffer from this uncertainty in the progenitor mass. We explore here two different masses proposed for its progenitor. We compute their pre-supernova characteristics using Modules for Experiments in Stellar Astrophysics (MESA). We show that a non-rotating star with the initial mass of 26 M$_{\odot}$ would require a considerable amount of circum-stellar medium (A$_V \sim$ 3) to obscure its high luminosity given the observed pre-explosion magnitudes detected by the Hubble Space Telescope (HST). Such a high value of visual extinction appears to be inconsistent with that derived for SN 2013ej as well as SN 2003gd in the same host galaxy M74. In contrast, the evolutionary models of lower mass (13 M$_{\odot}$) star are easily accommodated within the observed HST magnitudes. Some of the 26 M$_{\odot}$ simulations show luminosity variation in the last few years which could be discriminated by high cadence and multi-band monitoring of supernova candidates in nearby galaxies. We demonstrate that our calculations are well-resolved with adequate zoning and evolutionary time-steps.

We present the evolution of massive star progenitors of supernovae of type IIP. We take the example of the nearby and well-studied SN 2013ej. We explore how convective overshoot affects the stellar structure, surface abundances, and effective temperature of massive stars, using the Modules for Experiments in Stellar Astrophysics (MESA). In particular models with moderate overshoot ($f$ = 0.02 to 0.031) show the presence of blue loops in the Hertzsprung-Russell diagram with a red to blue [$log_{10}(T_{eff}/\rm K)$ from $< 3.6$ to $> 4.0$] excursion and transition back to red, during core helium burning phase. Models with overshoot outside this range of $f$ values kept the star in the red supergiant state throughout the post helium ignition phases. The surface CNO abundance shows enrichment post-main-sequence and again around the time when helium is exhausted in core. These evolutionary changes in surface CNO abundance are indistinguishable in the currently available observations due to large observational uncertainties. However, these observations may distinguish between the ratios of surface nitrogen to oxygen at different evolutionary stages of the star. We also compare the effects of convective overshoot on various parameters related to likelihood of explosion of a star as opposed to collapse to a black hole. These parameters are the compactness parameter, M$_4$, and $\mu_4$. Combination $\mu_4 \times $M$_4$, and $\mu_4$ have similar variation with $f$ and both peak out at $f$ = 0.032. We find that all of our 13 M$_{\odot}$ models are likely to explode.

The European Space Agency has invested heavily in two cornerstones missions; Herschel and Planck. The legacy data from these missions provides us with an unprecedented opportunity to study cosmic dust in galaxies so that we can answer fundamental questions about, for example: the origin of the chemical elements, physical processes in the interstellar medium (ISM), its effect on stellar radiation, its relation to star formation and how this relates to the cosmic far infrared background. In this paper we describe the DustPedia project, which is enabling us to develop tools and computer models that will help us relate observed cosmic dust emission to its physical properties (chemical composition, size distribution, temperature), to its origins (evolved stars, super novae, growth in the ISM) and the processes that destroy it (high energy collisions and shock heated gas). To carry out this research we will combine the Herschel/Planck data with that from other sources of data, providing observations at numerous wavelengths (< 41) across the spectral energy distribution, thus creating the DustPedia database. To maximise our spatial resolution and sensitivity to cosmic dust we limit our analysis to 4231 local galaxies (v < 3000 km/s) selected via their near infrared luminosity (stellar mass). To help us interpret the data we have developed a new physical model for dust (THEMIS), a new Bayesian method of fitting and interpreting spectral energy distributions (HerBIE) and a state-of-the-art Monte Carlo photon tracing radiative transfer model (SKIRT). In this the first of the DustPedia papers we describe the project objectives, data sets used and provide an insight into the new scientific methods we plan to implement.

We present newly processed archival Herschel images of molecular cloud MCLD 123.5+24.9 in the Polaris Flare. This cloud contains five starless cores. Using the spectral synthesis code Cloudy, we explore uncertainties in the derivation of column densities, hence, masses of molecular cores from Herschel data. We first consider several detailed grain models that predict far-IR grain opacities. Opacities predicted by the models differ by more than a factor of two, leading to uncertainties in derived column densities by the same factor. Then we consider uncertainties associated with the modified blackbody fitting process used by observers to estimate column densities. For high column density clouds (N(H) $\gg$ 10$^{22}$ cm$^{-2}$), this fitting technique can underestimate column densities by about a factor of three. Finally, we consider the virial stability of the five starless cores in MCLD 123.5+24.9. All of these cores appear to have strongly sub-virial masses, assuming, as we argue, that $^{13}$CO line data provide reliable estimates of velocity dispersions. Evidently, they are not self-gravitating, so it is no surprise that they are starless.

Dielectronic recombination (DR) is the dominant recombination process for most heavy elements in photoionized clouds. Accurate DR rates for a species can be predicted when the positions of autoionizing states are known. Unfortunately such data are not available for most third and higher-row elements. This introduces an uncertainty that is especially acute for photoionized clouds, where the low temperatures mean that DR occurs energetically through very low-lying autoionizing states. This paper discusses S$^{2+} \rightarrow$ S$^+$ DR, the process that is largely responsible for establishing the [S~III]/[S~II] ratio in nebulae. We derive an empirical rate coefficient using a novel method for second-row ions, which do have accurate data. Photoionization models are used to reproduce the [O~III] / [O~II] / [O~I] / [Ne~III] intensity ratios in central regions of the Orion Nebula. O and Ne have accurate atomic data and can be used to derive an empirical S$^{2+} \rightarrow$ S$^+$ DR rate coefficient at $\sim 10^{4}$~K. We present new calculations of the DR rate coefficient for S$^{2+} \rightarrow$ S$^+$ and quantify how uncertainties in the autoionizing level positions affect it. The empirical and theoretical results are combined and we derive a simple fit to the resulting rate coefficient at all temperatures for incorporation into spectral synthesis codes. This method can be used to derive empirical DR rates for other ions, provided that good observations of several stages of ionization of O and Ne are available.