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.

Decay of surface plasmons to hot carriers finds a wide variety of applications in energy conversion, photocatalysis and photodetection. However, a detailed theoretical description of plasmonic hot-carrier generation in real materials has remained incomplete. Here we report predictions for the prompt distributions of excited ‘hot’ electrons and holes generated by plasmon decay, before inelastic relaxation, using a quantized plasmon model with detailed electronic structure. We find that carrier energy distributions are sensitive to the electronic band structure of the metal: gold and copper produce holes hotter than electrons by 1–2 eV, while silver and aluminium distribute energies more equitably between electrons and holes. Momentum-direction distributions for hot carriers are anisotropic, dominated by the plasmon polarization for aluminium and by the crystal orientation for noble metals. We show that in thin metallic films intraband transitions can alter the carrier distributions, producing hotter electrons in gold, but interband transitions remain dominant.

Nanoscale localization of electromagnetic fields near metallic nanostructures underpins the fundamentals and applications of plasmonics. The unavoidable energy loss from plasmon decay, initially seen as a detriment, has now expanded the scope of plasmonic applications to exploit the generated hot carriers. However, quantitative understanding of the spatial localization of these hot carriers, akin to electromagnetic near-field maps, has been elusive. Here we spatially map hot-electron-driven reduction chemistry with 15 nm resolution as a function of time and electromagnetic field polarization for different plasmonic nanostructures. We combine experiments employing a six-electron photo-recycling process that modify the terminal group of a self-assembled monolayer on plasmonic silver nanoantennas, with theoretical predictions from first-principles calculations of non-equilibrium hot-carrier transport in these systems. The resulting localization of reactive regions, determined by hot-carrier transport from high-field regions, paves the way for improving efficiency in hot-carrier extraction science and nanoscale regio-selective surface chemistry.

The angular momentum (AM) evolution of stellar interiors, along with the resulting rotation rates of stellar remnants, remains poorly understood. Asteroseismic measurements of red giant stars reveal that their cores rotate much faster than their surfaces, but much slower than theoretically predicted, indicating an unidentified source of AM transport operates in their radiative cores. Motivated by this, we investigate the magnetic Tayler instability and argue that it saturates when turbulent dissipation of the perturbed magnetic field energy is equal to magnetic energy generation via winding. This leads to larger magnetic field amplitudes, more efficient AM transport, and smaller shears than predicted by the classic Tayler-Spruit dynamo. We provide prescriptions for the effective AM diffusivity and incorporate them into numerical stellar models, finding they largely reproduce (1) the nearly rigid rotation of the Sun and main sequence stars, (2) the core rotation rates of low-mass red giants during hydrogen shell and helium burning, and (3) the rotation rates of white dwarfs. We discuss implications for stellar rotational evolution, internal rotation profiles, rotational mixing, and the spins of compact objects.