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 nonlinear 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, as well as 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.

WASP-12b is a transiting hot Jupiter on a 1.09 day orbit around a late-F star. Since the planet’s discovery in 2008, the time interval between transits has been decreasing by 29 ± 2 ms yr⁻¹. This is a possible sign of orbital decay, although the previously available data left open the possibility that the planet’s orbit is slightly eccentric and is undergoing apsidal precession. Here, we present new transit and occultation observations that provide more decisive evidence for orbital decay, which is favored over apsidal precession by a of 22.3 or Bayes factor of 70,000. We also present new radial-velocity data that rule out the Rømer effect as the cause of the period change. This makes WASP-12 the first planetary system for which we can be confident that the orbit is decaying. The decay timescale for the orbit is . Interpreting the decay as the result of tidal dissipation, the modified stellar tidal quality factor is .

We present new 3.6 and 4.5 μm Spitzer phase curves for the highly irradiated hot Jupiter WASP-33b and the unusually dense Saturn-mass planet HD 149026b. As part of this analysis, we develop a new variant of pixel-level decorrelation that is effective at removing intrapixel sensitivity variations for long observations (>10 hr) where the position of the star can vary by a significant fraction of a pixel. Using this algorithm, we measure eclipse depths, phase amplitudes, and phase offsets for both planets at 3.6 and 4.5 μm. We use a simple toy model to show that WASP-33b’s phase offset, albedo, and heat recirculation efficiency are largely similar to those of other hot Jupiters despite its very high irradiation. On the other hand, our fits for HD 149026b prefer a very high albedo. We also compare our results to predictions from general circulation models, and we find that while neither planet matches the models well, the discrepancies for HD 149026b are especially large. We speculate that this may be related to its high bulk metallicity, which could lead to enhanced atmospheric opacities and the formation of reflective cloud layers in localized regions of the atmosphere. We then place these two planets in a broader context by exploring relationships between the temperatures, albedos, heat transport efficiencies, and phase offsets of all planets with published thermal phase curves. We find a striking relationship between phase offset and irradiation temperature: the former drops with increasing temperature until around 3400 K and rises thereafter. Although some aspects of this trend are mirrored in the circulation models, there are notable differences that provide important clues for future modeling efforts.