In early-type stars a fossil magnetic field may be generated during the star formation process or be the result of a stellar merger event. Surface magnetic fields are thought to be erased by (sub)surface convection layers, which typically leave behind weak disordered fields. However, if the fossil field is strong enough it can prevent the onset of (sub)surface convection and so be preserved onto the main sequence. We calculate the critical field strength at which this occurs, and find that it corresponds well with the lower limit amplitude of observed fields in strongly magnetised Ap/Bp stars ($\approx$ 300 G). The critical field strength is predicted to increase slightly during the main sequence evolution, which could also explain the observed decline in the fraction of magnetic stars. This supports the conclusion that the bimodal distribution of observed magnetic fields in early-type stars reflects two different field origin stories: strongly magnetic fields are fossils fields inherited from star formation or a merger event, and weak fields are the product of on-going dynamo action.

A common situation in galactic and intergalactic gas involves cold dense gas in motion relative to hot diffuse gas. Kelvin-Helmholtz instability creates a turbulent mixing layer and populates the intermediate-temperature phase, which often cools rapidly. The energy lost to cooling is balanced by the advection of hot high enthalpy gas into the mixing layer, resulting in growth and acceleration of the cold phase. This process may play a major role in determining the interstellar medium and circumgalactic medium phase structure, and accelerating cold gas in galactic winds and cosmic filaments. Cooling in these mixing layers occurs in a thin corrugated sheet, which we argue has an area with fractal dimension $D=5/2$ and a thickness that adjusts to match the hot phase mixing time to the cooling time. These cooling sheet properties form the basis of a new model for how the cooling rate and hot gas inflow velocity depend on the size $L$, cooling time $t_{\rm cool}$, relative velocity $v_{\rm rel}$, and density contrast $\rho_{\rm cold}/\rho_{\rm hot}$ of the system. Entrainment is expected to be enhanced in environments with short $t_{\rm cool}$, large $v_{\rm rel}$, and large $\rho_{\rm cold}/\rho_{\rm hot}$. Using a large suite of three dimensional hydrodynamic simulations, we demonstrate that this fractal cooling layer model accurately captures the energetics and evolution of turbulent interfaces and can therefore be used as a foundation for understanding multiphase mixing with strong radiative cooling.

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.

Over time, tides synchronize the rotation periods of stars in a binary system to the orbital period. However, if the star exhibits differential rotation then only a portion of it can rotate at the orbital period, so the rotation period at the surface may not match the orbital period. The difference between the rotation and orbital periods can therefore be used to infer the extent of the differential rotation. We use a simple parameterization of differential rotation in stars with convective envelopes in circular orbits to predict the difference between the surface rotation period and the orbital period. Comparing this parameterization to observed eclipsing binary systems, we find that in the surface convection zones of stars in short-period binaries there is very little radial differential rotation, with $|r\partial_r \ln \Omega| < 0.02$. This holds even for longer orbital periods, though it is harder to say which systems are synchronized at long periods, and larger differential rotation is degenerate with asynchronous rotation.

We argue that the interplay between cosmic rays, the initial mass function, and star formation plays a crucial role in regulating the star-forming "main sequence". To explore these phenomena we develop a toy model for galaxy evolution in which star formation is regulated by a combination of a temperature-dependent initial mass function and heating due to starlight, cosmic rays and, at very high redshift, the cosmic microwave background. This produces an attractor, near-equilibrium solution which is consistent with observations of the star-forming main sequence over a broad redshift range. Additional solutions to the same equations may correspond to other observed phases of galaxy evolution including quiescent galaxies. This model makes several falsifiable predictions, including higher metallicities and dust masses than anticipated at high redshift and isotopic abundances in the Milky Way. It also predicts that stellar mass-to-light ratios are lower than produced using a Milky Way-derived IMF, so that inferences of stellar masses and star formation rates for high redshift galaxies are overestimated. In some cases, this may also transform inferred dark matter profiles from core-like to cusp-like.

Sulfur is one of the most abundant elements in the Universe, with important roles in astro-, geo-, and biochemistry. Its main reservoirs in planet-forming disks have previously eluded detection: gaseous molecules only account for $<1\,$\% of total elemental sulfur, with the rest likely in either ices or refractory minerals. Mechanisms such as giant planets can filter out dust from gas accreting onto disk-hosting stars. For stars above 1.4 solar masses, this leaves a chemical signature on the stellar photosphere that can be used to determine the fraction of each element that is locked in dust. Here, we present an application of this method to sulfur, zinc, and sodium. We analyse the accretion-contaminated photospheres of a sample of young stars and find $(89\pm8)\,$\% of elemental sulfur is in refractory form in their disks. The main carrier is much more refractory than water ice, consistent with sulfide minerals such as FeS.

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.

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.

Convection in the cores of massive stars becomes anisotropic when they rotate. This anisotropy leads to a misalignment of the thermal gradient and the thermal flux, which in turn results in baroclinicity and circulation currents in the upper radiative zone. We show that this induces a much stronger meridional flow in the radiative zone than previously thought. This drives significantly enhanced mixing, though this mixing does not necessarily reach the surface. The extra mixing takes on a similar form to convective overshooting, and is relatively insensitive to the rotation rate above a threshold, and may help explain the large overshoot distances inferred from observations. This has significant consequences for the evolution of these stars by enhancing core-envelope mixing.

We present an approach to turbulence closure based on mixing length theory with three-dimensional fluctuations against a two-dimensional background. This model is intended to be rapidly computable for implementation in stellar evolution software and to capture a wide range of relevant phenomena with just a single free parameter, namely the mixing length. We incorporate magnetic, rotational, baroclinic and buoyancy effects exactly within the formalism of linear growth theories with nonlinear decay. We treat differential rotation effects perturbatively in the corotating frame using a novel controlled approximation which matches the time evolution of the reference frame to arbitrary order. We then implement this model in an efficient open source code and discuss the resulting turbulent stresses and transport coefficients. We demonstrate that this model exhibits convective, baroclinic and shear instabilities as well as the magnetorotational instability (MRI). It also exhibits non-linear saturation behaviour, and we use this to extract the asymptotic scaling of various transport coefficients in physically interesting limits.

Post-asymptotic giant branch (post-AGB) stars with discs are all binaries. Many of these binaries have orbital periods between 100 and 1000 days so cannot have avoided mass transfer between the AGB star and its companion, likely through a common-envelope type interaction. We report on preliminary results of our project to model circumbinary discs around post-AGB stars using our binary population synthesis code binary_c. We combine a simple analytic thin-disc model with binary stellar evolution to estimate the impact of the disc on the binary, and vice versa, fast enough that we can model stellar populations and hence explore the rather uncertain parameter space involved with disc formation. We find that, provided the discs form with sufficient mass and angular momentum, and have an inner edge that is relatively close to the binary, they can both prolong the life of their parent post-AGB star and pump the eccentricity of orbits of their inner binaries.

We argue that an increased temperature in star-forming clouds alters the stellar initial mass function to be more bottom-light than in the Milky Way. At redshifts $z \gtrsim 6$, heating from the cosmic microwave background radiation produces this effect in all galaxies, and it is also present at lower redshifts in galaxies with very high star formation rates (SFRs). A failure to account for it means that at present, photometric template fitting likely overestimates stellar masses and star formation rates for the highest-redshift and highest-SFR galaxies. In addition this may resolve several outstanding problems in the chemical evolution of galactic halos.

Protoplanetary disks, debris disks, and disrupted or evaporating planets can all feed accretion onto stars. The photospheric abundances of such stars may then reveal the composition of the accreted material. This is especially likely in B to mid-F type stars, which have radiative envelopes and hence less bulk--photosphere mixing. We present a theoretical framework (\texttt{CAM}) considering diffusion, rotation, and other stellar mixing mechanisms, to describe how the accreted material interacts with the bulk of the star. This allows the abundance pattern of the circumstellar material to be calculated from measured stellar abundances and parameters ($v_{\rm rot}$, $T_{\rm eff}$). We discuss the \lboo\ phenomenon and the application of \texttt{CAM} on stars hosting protoplanetary disks (HD~100546, HD~163296), debris disks (HD~141569, HD~21997), and evaporating planets (HD~195689/KELT-9).

Telescopes are much more expensive than astronomers, so it is essential to minimize required sample sizes by using the most data-efficient statistical methods possible. However, the most commonly used model-independent techniques for finding the relationship between two variables in astronomy are flawed. In the worst case they can lead without warning to subtly yet catastrophically wrong results, and even in the best case they require more data than necessary. Unfortunately, there is no single best technique for nonparametric regression. Instead, we provide a guide for how astronomers can choose the best method for their specific problem and provide a python library with both wrappers for the most useful existing algorithms and implementations of two new algorithms developed here.

We study the interaction between stellar irradiation and tidal heating in gaseous planets with short orbital periods. The intentionally simplified atmospheric model we employ makes the problem analytically tractable and permits the derivation of useful scaling relations. We show that many tidal models provide thermal feedback, producing interior radiative zones and leading to enhanced g-mode dissipation with a wide spectrum of resonances. These resonances are dynamically tuned by the thermal feedback, and so represent a novel form of thermomechanical feedback, coupling vibrational modes to the very slow thermal evolution of the planet. We then show that stellar irradiation allows the heat produced by these modes to be trapped at depth with high efficiency, leading to entropy increase in the central convective region, as well as expansion of the planet's radius sufficient to match observed swelling. We find that thermally driven winds play an essential role in this process by making the thermal structure of the atmosphere spherically symmetric within a few scale heights of the photosphere. We characterise the relationship between the swelling factor, the orbital period and the host star and determine the timescale for swelling. We show that these g-modes suffice to produce bloating on the order of the radius of the planet over $\mathrm{Gyr}$ timescales when combined with significant insolation and we provide analytic relations for the relative magnitudes of tidal heating and insolation.