We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001%. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen-and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development.

We examine the properties of white dwarfs (WDs) accreting hydrogen-rich matter in and near the stable burning regime of accretion rates as modeled by time-dependent calculations done with Modules for Experiments in Stellar Astrophysics (MESA). We report the stability boundary for WDs of masses between 0.51 M and 1.34 M as found via time-dependent calculations. We also examine recurrent novae that are accreting at rates close to, but below, the stable burning limit and report their recurrence times. Our dense grid in accretion rates finds the expected minimum possible recurrence times as a function of the WD mass. This enables inferences to be made about the minimum WD mass possible to reach a specific recurrence time. We compare our computational models of post-outburst novae to the stably burning WDs and explicitly calculate the duration and effective temperature (T eff ) of the post-novae WD in the supersoft phase. We agree with the measured turnoff time -T eff relation in M31 by Henze and collaborators, infer WD masses in the 1.0-1.3 M range, and predict ejection masses consistent with those observed. We close by commenting on the importance of the hot helium layer generated by stable or unstable hydrogen burning for the short-and long-term evolution of accreting white dwarfs.

We present observations of four rapidly rising (t rise ≈ 10 d) transients with peak luminosities between those of supernovae (SNe) and superluminous SNe (M peak ≈ −20) -one discovered and followed by the Palomar Transient Factory (PTF) and three by the Supernova Legacy Survey (SNLS). The light curves resemble those of SN 2011kl, recently shown to be associated with an ultra-long-duration gamma ray burst (GRB), though no GRB was seen to accompany our SNe. The rapid rise to a luminous peak places these events in a unique part of SN phase space, challenging standard SN emission mechanisms. Spectra of the PTF event formally classify it as a Type II SN due to broad Hα emission, but an unusual absorption feature, which can be interpreted as either high velocity Hα (though deeper than in previously known cases) or Si II (as seen in Type Ia SNe), is also observed. We find that existing models of white dwarf detonations, CSM interaction, shock breakout in a wind (or steeper CSM) and magnetar spindown can not readily explain the observations. We consider the possibility that a "Type 1.5 SN" scenario could be the origin of our events. More detailed models for these kinds of transients and more constraining observations of future such events should help better determine their nature.

The critical role in surface reactions and heterogeneous catalysis of metal atoms with low coordination numbers, such as found at atomic steps and surface defects, is firmly established. But despite the growing availability of tools that enable detailed in situ characterization, so far it has not been possible to document this role directly. Surface properties can be mapped with high spatial resolution, and catalytic conversion can be tracked with a clear chemical signature; however, the combination of the two, which would enable high-spatial-resolution detection of reactions on catalytic surfaces, has rarely been achieved. Single-molecule fluorescence spectroscopy has been used to image and characterize single turnover sites at catalytic surfaces, but is restricted to reactions that generate highly fluorescing product molecules. Herein the chemical conversion of N-heterocyclic carbene molecules attached to catalytic particles is mapped using synchrotron-radiation-based infrared nanospectroscopy with a spatial resolution of 25 nanometres, which enabled particle regions that differ in reactivity to be distinguished. These observations demonstrate that, compared to the flat regions on top of the particles, the peripheries of the particles-which contain metal atoms with low coordination numbers-are more active in catalysing oxidation and reduction of chemically active groups in surface-anchored N-heterocyclic carbene molecules.