Cepheid masses derived from pulsations or binary dynamics are generally lower than those derived from stellar evolution models. Recent efforts have been dedicated to investigating the effects of abundances, mass loss, rotation, convection and overshooting prescriptions for modifying the evolution tracks to reduce or remove this Cepheid mass discrepancy. While these approaches are promising, either alone or in combination, more work is required to distinguish between possible solutions. Here we investigate nuclear reaction rate and opacity modifications on Cepheid evolution using the MESA code. We discuss the effects of opacity increases at envelope temperatures of 200,000-400,000 K proposed to explain the pulsation properties of hybrid main-sequence beta Cep/Slowly Pulsating B (SPB) variables which will evolve into Cepheids. We make use of the RSP nonlinear radial pulsation modeling capability in MESA to calculate periods and radial velocity amplitudes of Galactic Cepheids V1334 Cyg, Polaris, and delta Cep.

We explore the evolution of a select grid of solar metallicity stellar models from their pre-main sequence phase to near their final fates in a neutrino Hertzsprung-Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated \MESA\ solar model for the solar neutrino luminosity ($L_{\nu,\odot}$ = 0.02398 $\cdot$ $L_{\gamma,\odot}$ = 9.1795 $\times$ 10$^{31}$ erg s$^{-1}$) as a normalization, we identify $\simeq$ 0.3 MeV electron neutrino emission from helium burning during the helium flash (peak $L_{\nu} / L_{\nu,\odot} \simeq$ 10$^4$, flux $\Phi_{\nu, {\rm He \ flash}} \simeq$ 170 (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$ for a star located at a distance of $d$ parsec, timescale $\simeq$ 3 days) and the thermal pulse (peak $L_{\nu} / L_{\nu,\odot} \simeq$ 10$^9$, flux $\Phi_{\nu, {\rm TP}} \simeq$ 1.7$\times$10$^7$ (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$, timescale $\simeq$ 0.1 yr) phases of evolution in low mass stars as potential probes for stellar neutrino astronomy. We also delineate the contribution of neutrinos from nuclear reactions and thermal processes to the total neutrino loss along the stellar tracks in a neutrino Hertzsprung-Russell diagram. We find, broadly but with exceptions, that neutrinos from nuclear reactions dominate whenever hydrogen and helium burn, and that neutrinos from thermal processes dominate otherwise.

Los Alamos National Laboratory has calculated a new generation of radiative opacities (OPLIB data using the ATOMIC code) for elements with atomic number Z=1-30 with improved physics input, updated atomic data, and finer temperature grid to replace the Los Alamos LEDCOP opacities released in the year 2000. We calculate the evolution of standard solar models including these new opacities, and compare with models evolved using the Lawrence Livermore National Laboratory OPAL (Iglesias and Rogers 1996) opacities. We use the solar abundance mixture of Asplund et al. (2009). The new Los Alamos ATOMIC opacities have steeper opacity derivatives than those of OPAL for temperatures and densities of the solar interior radiative zone. We compare the calculated nonadiabatic solar oscillation frequencies and solar interior sound speed to observed frequencies and helioseismic inferences. The calculated sound-speed profiles are similar for models evolved using either the updated Iben evolution code (see \cite{Guzik2010}), or the MESA evolution code (Paxton et al., 2015). The LANL ATOMIC opacities partially mitigate the "solar abundance problem".