Search

Recent observations have found a growing number of hypervelocity stars with speeds of $\approx 1500-2500\,$km\,s$^{-1}$ which could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been as white dwarfs filling their Roche lobe. While many simulations exist studying the dynamical phase leading to supernova detonation in these systems, no detailed calculations of the long-term structure of the runaways have yet been performed. We use an existing \textsc{Arepo} hydrodynamical simulation of a supernova in a white dwarf binary as a starting point for the evolution of these stars with the 1 dimensional stellar evolution code MESA. We show that the supernova shock is not enough to inflate the white dwarf over timescales longer than a few thousand years, significantly shorter than the $10^{5-6}$ year lifetimes inferred for observed hypervelocity runaways. Despite experiencing a shock from a supernova less than $\approx 0.02\,R_\odot$ away, our models do not experience significant interior heating, and all contract back to radii around $0.01\,R_\odot$ within about $10^4$\,years. Explaining the observed inflated states requires either an additional source of significant heating or some other physics that is not yet accounted for in the subsequent evolution., Comment: Submitted to A\&A. 15 pages, 17 figures

We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10$^{-6}$ to $Z$\,=\,0.2 which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity range $Z$\,=\,10$^{-4}$ to $Z$\,=\,0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions between $X$\,=\,0 and $X$\,=\,0.1 including $X$\,=\,$10^{-2}, 10^{-3}, 10^{-4}, 10^{-5}, 10^{-6}$ that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables in \MESA, and find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from $\approx$\,20--80\% across individual chemical mixtures, up to $\approx$\,8\% and $\approx$\,15\% at the bottom and top of the solar convection zone respectively, and $\approx$\,7\% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity., Comment: 28 pages, 12 (13) figures. Accepted for Publication in ApJ

Eruptive mass loss in massive stars is known to occur, but the mechanism(s) are not yet well-understood. One proposed physical explanation appeals to opacity-driven super-Eddington luminosities in stellar envelopes. Here, we present a 1D model for eruptive mass loss and implement this model in the MESA stellar evolution code. The model identifies regions in the star where the energy associated with a locally super-Eddington luminosity exceeds the binding energy of the overlaying envelope. The material above such regions is ejected from the star. Stars with masses $10-100~M_\odot$ at solar and SMC metallicities are evolved through core helium burning, with and without this new eruptive mass-loss scheme. We find that eruptive mass loss of up to $\sim10^{-2}~M_\odot \mathrm{yr}^{-1}$ can be driven by this mechanism, and occurs in a vertical band on the HR diagram between $3.5 \lesssim \log(T_\mathrm{eff}/\mathrm{K}) \lesssim 4.0$. This predicted eruptive mass loss prevents stars of initial masses $\gtrsim20~M_\odot$ from evolving to become red supergiants, with the stars instead ending their lives as blue supergiants, and therefore offers a possible explanation for the observed lack of red supergiants in that mass regime., Comment: 17 pages, 9 figures, accepted to ApJ

Blue large-amplitude pulsators (BLAPs) make up a rare class of hot pulsating stars with effective temperatures of $\approx$30,000 K and surface gravities of 4.0 - 5.0 dex (cgs). The evolutionary origin and current status of BLAPs is not well understood, largely based on a lack of spectroscopic observations and no available mass constraints. However, several theoretical models have been proposed that reproduce their observed properties, including studies that identify them as pulsating helium-core pre-white dwarfs (He-core pre-WDs). We present here follow-up high-speed photometry and phase-resolved spectroscopy of one of the original 14 BLAPs, OGLE-BLAP-009, discovered during the Optical Gravitational Lensing Experiment. We aim to explore its pulsation characteristics and determine stellar properties such as mass and radius in order to test the consistency of these results with He-core pre-WD models. Using the mean atmospheric parameters found using spectroscopy, we fit a spectral energy distribution to obtain a preliminary estimate of the radius, luminosity and mass by making use of the Gaia parallax. We then compare the consistency of these results to He-core pre-WD models generated using MESA, with predicted pulsation periods implemented using GYRE. We find that our mass constraints are in agreement with a low-mass He-core pre-WD of $\approx$0.30 M$_{\odot}$., Comment: Accepted for publication in MNRAS, 15 pages, 11 figures, 5 tables

Regular, automated testing is a foundational principle of modern software development. Numerous widely-used continuous integration systems exist, but they are often not suitable for the unique needs of scientific simulation software. Here we describe the testing infrastructure developed for and used by the Modules for Experiments in Stellar Astrophysics (MESA) project. This system allows the computationally-demanding MESA test suite to be regularly run on a heterogeneous set of computers and aggregates and displays the testing results in a form that allows for the rapid identification and diagnosis of regressions. Regularly collecting comprehensive testing data also enables longitudinal studies of the performance of the software and the properties of the models it generates., Comment: 12 page, 7 figures, Accepted to ApJS

We present a detailed modelling study of CD-30$^{\circ}$11223 (CD-30), a hot subdwarf (sdB)-white dwarf (WD) binary identified as a double detonation supernova progenitor, using the open-source stellar evolution software MESA. We focus on implementing binary evolution models carefully tuned to match the observed characteristics of the system including $\log g$ and $T_{\rm eff}$. For the first time, we account for the structure of the hydrogen envelope throughout the modelling, and find that the inclusion of element diffusion is important for matching the observed radius and temperature. We investigate the two sdB mass solutions (0.47 and 0.54 $M_{\odot}$) previously proposed for this system, strongly favouring the 0.47 $M_{\odot}$ solution. The WD cooling age is compared against the sdB age using our models, which suggest an sdB likely older than the WD, contrary to the standard assumption for compact sdB-WD binaries. Subsequently, we propose a possible alternate formation channel for CD-30. We also perform binary evolution modelling of the system to study various aspects such as mass transfer, orbital period evolution and luminosity evolution. Our models confirm CD-30 as a double detonation supernova progenitor, expected to explode $\approx55$ Myr from now. The WD accretes a $\approx0.17$ $M_{\odot}$ thick helium shell that causes a detonation, leaving a 0.30 $M_{\odot}$ sdB ejected at $\approx$750 km/s. The final 15 Myr of the system are characterised by helium accretion which dominates the system luminosity, possibly resembling an AM CVn-type system., Comment: 12 pages, 8 figures. Accepted for publication in MNRAS

The polluted white dwarf (WD) system SDSS J122859.93+104032.9 (SDSS J1228) shows variable emission features interpreted as originating from a solid core fragment held together against tidal forces by its own internal strength, orbiting within its surrounding debris disk. Estimating the size of this orbiting solid body requires modeling the accretion rate of the polluting material that is observed mixing into the WD surface. That material is supplied via sublimation from the surface of the orbiting solid body. The sublimation rate can be estimated as a simple function of the surface area of the solid body and the incident flux from the nearby hot WD. On the other hand, estimating the accretion rate requires detailed modeling of the surface structure and mixing in the accreting WD. In this work, we present MESA WD models for SDSS J1228 that account for thermohaline instability and mixing in addition to heavy element sedimentation to accurately constrain the sublimation and accretion rate necessary to supply the observed pollution. We derive a total accretion rate of $\dot M_{\rm acc}=1.8\times 10^{11}\,\rm g\,s^{-1}$, several orders of magnitude higher than the $\dot M_{\rm acc}=5.6\times 10^{8}\,\rm g\,s^{-1}$ estimate obtained in earlier efforts. The larger mass accretion rate implies that the minimum estimated radius of the orbiting solid body is r$_{\rm{min}}$ = 72 km, which, although significantly larger than prior estimates, still lies within upper bounds (a few hundred km) for which the internal strength could no longer withstand tidal forces from the gravity of the WD., Comment: 7 pages, 3 figures, 4 tables, accepted for publication in ApJ

We report the discovery of ZTF J0127+5258, a compact mass-transferring binary with an orbital period of 13.7 minutes. The system contains a white dwarf accretor, which likely originated as a post-common envelope carbon-oxygen (CO) white dwarf, and a warm donor ($T_{\rm eff,\,donor}= 16,400\pm1000\,\rm K$). The donor probably formed during a common envelope phase between the CO white dwarf and an evolving giant which left behind a helium star or helium white dwarf in a close orbit with the CO white dwarf. We measure gravitational wave-driven orbital inspiral with $\sim 35\sigma$ significance, which yields a joint constraint on the component masses and mass transfer rate. While the accretion disk in the system is dominated by ionized helium emission, the donor exhibits a mixture of hydrogen and helium absorption lines. Phase-resolved spectroscopy yields a donor radial-velocity semi-amplitude of $771\pm27\,\rm km\, s^{-1}$, and high-speed photometry reveals that the system is eclipsing. We detect a {\it Chandra} X-ray counterpart with $L_{X}\sim 3\times 10^{31}\,\rm erg\,s^{-1}$. Depending on the mass-transfer rate, the system will likely evolve into either a stably mass-transferring helium CV, merge to become an R Crb star, or explode as a Type Ia supernova in the next million years. We predict that the Laser Space Interferometer Antenna (LISA) will detect the source with a signal-to-noise ratio of $24\pm6$ after 4 years of observations. The system is the first \emph{LISA}-loud mass-transferring binary with an intrinsically luminous donor, a class of sources that provide the opportunity to leverage the synergy between optical and infrared time domain surveys, X-ray facilities, and gravitational-wave observatories to probe general relativity, accretion physics, and binary evolution., Comment: 13 pages, 7 figures, 2 tables, submitted to ApJL

We enhance the treatment of crystallization for models of white dwarfs (WDs) in the stellar evolution software MESA by implementing carbon-oxygen (C/O) phase separation. The phase separation process during crystallization leads to transport of oxygen toward the center of WDs, resulting in a more compact structure that liberates gravitational energy as additional heating that modestly slows WD cooling timescales. We quantify this cooling delay in MESA C/O WD models over the mass range 0.5-1.0 $M_\odot$, finding delays of 0.5-0.8 Gyr for typical C/O interior profiles. MESA WD cooling timescales including this effect are generally comparable to other WD evolution models that make similar assumptions about input physics. When considering phase separation alongside $^{22}$Ne sedimentation, however, we find that both MESA and BaSTI WD cooling models predict a more modest sedimentation delay than the latest LPCODE models, and this may therefore require a re-evaluation of previously proposed solutions to some WD cooling anomalies that were based on LPCODE models of $^{22}$Ne sedimentation. Our implementation of C/O phase separation in the open-source stellar evolution software MESA provides an important tool for building realistic grids of WD cooling models, as well as a framework for expanding on our implementation to explore additional physical processes related to phase transitions and associated fluid motions in WD interiors., Comment: 15 pages, 9 figures, accepted for publication in ApJ

We report the results from follow-up observations of two Roche-lobe filling hot subdwarf binaries with white dwarf companions predicted to have accretion disks. ZTF J213056.71+442046.5 (ZTF J2130) with a 39-minute period and ZTF J205515.98+465106.5 (ZTF J2055) with a 56-minute period were both discovered as subdwarf binaries with light curves that could only be explained well by including an accretion disk in their models. We performed a detailed high-resolution spectral analysis using Keck/ESI to search for possible accretion features for both objects. We also employed polarimetric analysis using the Nordic Optical Telescope (NOT) for ZTF J2130. We did not find any signatures of an accretion disk in either object, and placed upper limits on the flux contribution and variation in degree of polarisation due to the disk. Owing to the short 39-minute period and availability of photometric data over six years for ZTF J2130, we conducted an extensive $O - C$ timing analysis in an attempt to look for orbital decay due to gravitational wave radiation. No such decay was detected conclusively, and a few more years of data paired with precise and consistent timing measurements were deemed necessary to constrain $\dot P$ observationally., Comment: 10 pages, 8 figures. Accepted for publication in MNRAS

We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron degenerate ignition events. We strengthen MESA's implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator split nuclear burning mode. We close by discussing major updates to MESA's software infrastructure that enhance source code development and community engagement., Comment: 50 pages, 29 figures; Accepted to ApJS

Stars with masses $\gtrsim 4 \times 10^{27}M_{\rm{moose}} \approx 1.1 M_\odot$ have core convection zones during their time on the main sequence. In these moosive stars, convection introduces many uncertainties in stellar modeling. In this Letter, we build upon the Boussinesq approximation to present the first-ever simulations of Moosinesq convection, which captures the complex geometric structure of the convection zones of these stars. These flows are bounded in a manner informed by the majestic terrestrial Alces alces (moose) and could have important consequences for the evolution of these stars. We find that Moosinesq convection results in very interesting flow morphologies and rapid heat transfer, and posit this as a mechanism of biomechanical thermoregulation., Comment: 7 pages, 3 figures, 1 moose, no logic in this place

The recently proposed "dynamically driven double-degenerate double-detonation" (D6) scenario posits that Type Ia supernovae (SNe) may occur during dynamically unstable mass transfer between two white dwarfs (WDs) in a binary. This scenario predicts that the donor WD may then survive the explosion and be released as a hypervelocity runaway, opening up the exciting possibility of identifying remnant stars from D6 SNe and using them to study the physics of detonations that produce Type Ia SNe. Three candidate D6 runaway objects have been identified in Gaia data. The observable runaway velocity of these remnant objects represents their orbital speed at the time of SN detonation. The orbital dynamics and Roche lobe geometry required in the D6 scenario place specific constraints on the radius and mass of the donor WD that becomes the hypervelocity runaway. In this letter, we calculate the radii required for D6 donor WDs as a function of the runaway velocity. Using mass-radius relations for WDs, we then constrain the masses of the donor stars as well. With measured velocities for each of the three D6 candidate objects based on Gaia EDR3, this work provides a new probe of the masses and mass ratios in WD binary systems that produce SN detonations and hypervelocity runaways., Comment: 8 pages, 4 figures, accepted for publication in ApJL

We present the discovery of a new double detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with an P=76.34179(2) min orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits we find a low mass sdB star, $M_{\rm sdB}=0.383\pm0.028$ M$_\odot$ with a massive white dwarf companion, $M_{\rm WD}=0.725\pm0.026$ M$_\odot$. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of $\approx$25 Myrs whereas our MESA model predicts an sdB age of $\approx$170 Myrs. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion $\approx$25 Myrs ago. Using the MESA stellar evolutionary code we find that the sdB star will start mass transfer in $\approx$6 Myrs and in $\approx$60 Myrs the white dwarf will reach a total mass of $0.92$ M$_\odot$ with a thick helium layer of $0.17$ M$_\odot$. This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 2238+7430 is only the second confirmed candidate for a double detonation thermonuclear supernova. Using both systems we estimate that at least $\approx$1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions., Comment: accepted for publication in ApJL, 14 pages, 5 figures, 1 table

A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back $\approx 10^5$ years to a known supernova remnant. Here we report the discovery of carbon-rich circumstellar material around LP 398-9, revealed by a strong infrared excess and analyzed with follow-up spectroscopy. The circumstellar material is most plausibly composed of inflated layers from the star itself, mechanically and radioactively heated by the past companion's supernova. We also detect a 15.4 hr periodic signal in the UV and optical light curves of LP 398-9, which we interpret as surface rotation. The rotation rate is consistent with theoretical predictions from this supernova mechanism, and the brightness variations could originate from surface inhomogeneity deposited by the supernova itself. Our observations strengthen the case for this double-degenerate SNIa progenitor channel, and motivate the search for more runaway SNIa donors., Comment: Accepted to MNRAS

We report the discovery of SDSS J133725.26+395237.7 (hereafter SDSS J1337+3952), a double-lined white dwarf (WD+WD) binary identified in early data from the fifth generation Sloan Digital Sky Survey (SDSS-V). The double-lined nature of the system enables us to fully determine its orbital and stellar parameters with follow-up Gemini spectroscopy and Swift UVOT ultraviolet fluxes. The system is nearby ($d = 113$ pc), and consists of a $0.51\, M_\odot$ primary and a $0.32\, M_\odot$ secondary. SDSS J1337+3952 is a powerful source of gravitational waves in the millihertz regime, and will be detectable by future space-based interferometers. Due to this gravitational wave emission, the binary orbit will shrink down to the point of interaction in $\approx 220$ Myr. The inferred stellar masses indicate that SDSS J1337+3952 will likely not explode as a Type Ia supernova (SN Ia). Instead, the system will probably merge and evolve into a rapidly rotating helium star, and could produce an under-luminous thermonuclear supernova along the way. The continuing search for similar systems in SDSS-V will grow the statistical sample of double-degenerate binaries across parameter space, constraining models of binary evolution and SNe Ia., Comment: Accepted to ApJ

Recently, a class of Roche-lobe-filling binary systems consisting of hot subdwarf stars and white dwarfs with sub-hour periods has been discovered. At present, the hot subdwarf is in a shell He burning phase and is transferring some of its remaining thin H envelope to its white dwarf companion. As the evolution of the hot subdwarf continues, it is expected to detach, leaving behind a low mass C/O core white dwarf secondary with a thick He layer. Then, on a timescale of $\sim 10$ Myr, gravitational wave radiation will again bring the systems into contact. If the mass transfer is unstable and results in a merger and a catastrophic thermonuclear explosion is not triggered, it creates a remnant with a C/O-dominated envelope, but one still rich enough in He to support an R Corona Borealis-like shell burning phase. We present evolutionary calculations of this phase and discuss its potential impact on the cooling of the remnant white dwarf., Comment: 7 pages, 4 figures; Accepted to ApJ

Binary systems of a hot subdwarf B (sdB) star + a white dwarf (WD) with orbital periods less than 2-3 hours can come into contact due to gravitational waves and transfer mass from the sdB star to the WD before the sdB star ceases nuclear burning and contracts to become a WD. Motivated by the growing class of observed systems in this category, we study the phases of mass transfer in these systems. We find that because the residual outer hydrogen envelope accounts for a large fraction of an sdB star's radius, sdB stars can spend a significant amount of time ($\sim$10s of Myr) transferring this small amount of material at low rates ($\sim 10^{-10}$-$10^{-9}\ M_\odot\,\rm yr^{-1}$) before transitioning to a phase where the bulk of their He transfers at much faster rates ($\gtrsim 10^{-8}\ M_\odot\,\rm yr^{-1}$). These systems therefore spend a surprising amount of time with Roche-filling sdB donors at orbital periods longer than the range associated with He star models without an envelope. We predict that the envelope transfer phase should be detectable by searching for ellipsoidal modulation of Roche-filling objects with $P_{\rm orb}=30$-$100$ min and $T_{\rm eff}=20{,}000$-$30{,}000$ K, and that many ($\geq$10) such systems may be found in the Galactic plane after accounting for reddening. We also argue that many of these systems may go through a phase of He transfer that matches the signatures of AM CVn systems, and that some AM CVn systems associated with young stellar populations likely descend from this channel., Comment: 16 pages, 8 figures, accepted for publication in ApJ

White dwarfs with metal-polluted atmospheres have been studied widely in the context of the accretion of rocky debris from evolved planetary systems. One open question is the geometry of accretion and how material arrives and mixes in the white dwarf surface layers. Using the 3D radiation-hydrodynamics code CO$^5$BOLD, we present the first transport coefficients in degenerate star atmospheres which describe the advection-diffusion of a passive scalar across the surface-plane. We couple newly derived horizontal diffusion coefficients with previously published vertical diffusion coefficients to provide theoretical constraints on surface spreading of metals in white dwarfs. Our grid of 3D simulations probes the vast majority of the parameter space of convective white dwarfs, with pure-hydrogen atmospheres in the effective temperature range 6000-18000 K and pure-helium atmospheres in the range 12000-34000 K. Our results suggest that warm hydrogen-rich atmospheres (DA; $\gtrsim$13000 K) and helium-rich atmospheres (DB, DBA; $\gtrsim$30000 K) are unable to efficiently spread the accreted metals across their surface, regardless of the time dependence of accretion. This result may be at odds with the current non-detection of surface abundance variations at white dwarfs with debris discs. For cooler hydrogen- and helium-rich atmospheres, we predict a largely homogeneous distribution of metals across the surface within a vertical diffusion timescale. This is typically less than 0.1 per cent of disc lifetime estimates, a quantity which is revisited in this paper using the overshoot results. These results have relevance for studies of the bulk composition of evolved planetary systems and models of accretion disc physics., Comment: Accepted for publication in MNRAS. 22 pages, 17 figures

We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models without $^{22}$Ne are fit to the observed periods and other properties. The root-mean-square residuals to the $\simeq$ 150$-$400 s low-order g-mode periods are typically in the range of $\sigma_{rms}$ $\lesssim$ 0.3 s, for a fit precision of $\sigma_{rms}/ P$ $\lesssim$ 0.3 %. We find average relative period shifts of $\Delta P/P$ $\simeq$ $\pm$ 0.5 % for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of $\Delta P/P$ depending on the specific g-mode, abundance of $^{22}$Ne, effective temperature, and mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when $^{22}$Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of $^{22}$Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including $^{22}$Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model fitting processes., Comment: 23 pages, 19 figures total. Accepted to the ApJ. Presented results at 2021 AAS January meeting

Cooling white dwarfs (WDs) can yield accurate ages when theoretical cooling models fully account for the physics of the dense plasma of WD interiors. We use MESA to investigate cooling models for a set of massive and ultra-massive WDs (0.9-1.3 $M_\odot$) for which previous models fail to match kinematic age indicators based on Gaia DR2. We find that the WDs in this population can be explained as C/O cores experiencing unexpectedly rapid $^{22}$Ne sedimentation in the strongly liquid interior just prior to crystallization. We propose that this rapid sedimentation is due to the formation of solid clusters of $^{22}$Ne in the liquid C/O background plasma. We show that these heavier solid clusters sink faster than individual $^{22}$Ne ions and enhance the sedimentation heating rate enough to dramatically slow WD cooling. MESA models including our prescription for cluster formation and sedimentation experience cooling delays of $\approx$4 Gyr on the WD Q branch, alleviating tension between cooling ages and kinematic ages. This same model then predicts cooling delays coinciding with crystallization of 6 Gyr or more in lower mass WDs (0.6-0.8 $M_\odot$). Such delays are compatible with, and perhaps required by, observations of WD populations in the local 100 pc WD sample and the open cluster NGC 6791. These results motivate new investigations of the physics of strongly coupled C/O/Ne plasma mixtures in the strongly liquid state near crystallization and tests through comparisons with observed WD cooling., Comment: 23 pages, 14 figures, accepted for publication in ApJ

Stars with unusual properties can provide a wealth of information about rare stages of stellar evolution and exotic physics. However, determining the true nature of peculiar stars is often difficult. In this work, we conduct a systematic search for cool and luminous stars in the Magellanic Clouds with extreme variability, motivated by the properties of the unusual SMC star and Thorne-\.Zytkow Object (T\.ZO) candidate HV2112. Using light curves from ASAS-SN we identify 38 stars with surface temperatures T $<$ 4800K, luminosities $\log$(L/L$_\odot$) $>$ 4.3, variability periods $>$ 400 days, and variability amplitudes $\Delta$V $>$ 2.5 mag. Eleven of these stars possess the distinctive double-peaked light curve morphology of HV2112. We use the pulsation properties and derived occurrence rates for these 12 objects to constrain their nature. From comparisons to stellar populations and models, we find that one star may be a red supergiant with large amplitude pulsations. For the other 11 stars we derive current masses of $\sim$5-10 M$_{\odot}$, below the theoretical minimum mass of $\sim$15 M$_{\odot}$ for T\.ZOs to be stable, casting doubt on this interpretation. Instead, we find that the temperatures, luminosities, mass-loss rates, and periods of these stars are consistent with predictions for super-Asymptotic Giant Branch (s-AGB) stars that have begun carbon burning but have not reached the superwind phase. We infer lifetimes in this phase of $\sim($1$-$7) $\times$ 10$^{4}$ years, also consistent with an s-AGB interpretation. If confirmed, these objects would represent the first identified population of s-AGB stars, illuminating the transition between low- and high-mass stellar evolution., Comment: 29 pages, 14 figures, 4 tables, plus Appendix. Accepted for publication in the Astrophysical Journal

We present the discovery of the second binary with a Roche lobe-filling hot subdwarf transferring mass to a white dwarf (WD) companion. This 56 minute binary was discovered using data from the Zwicky Transient Facility. Spectroscopic observations reveal an He-sdOB star with an effective temperature of $T_{\rm eff}=33,700\pm1000$ K and a surface gravity of $log(g)=5.54\pm0.11$. The GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the He-sdOB star and shows an eclipse of the He-sdOB by an accretion disk as well as a weak eclipse of the WD. We infer a He-sdOB mass of $M_{\rm sdOB}=0.41\pm0.04$ M$_\odot$ and a WD mass of $M_{\rm WD}=0.68\pm0.05$ M$_\odot$. The weak eclipses imply a WD black-body temperature of $63,000\pm10,000$ K and a radius $R_{\rm WD}=0.0148\pm0.0020$ M$_\odot$ as expected for a WD of such high temperature. The He-sdOB star is likely undergoing hydrogen shell burning and will continue transferring mass for $\approx1$ Myrs at a rate of $10^{-9} M_\odot {\rm yr}^{-1}$ which is consistent with the high WD temperature. The hot subdwarf will then turn into a WD and the system will merge in $\approx30$ Myrs. We suggest that Galactic reddening could bias discoveries towards preferentially finding Roche lobe-filling systems during the short-lived shell burning phase. Studies using reddening corrected samples should reveal a large population of helium core-burning hot subdwarfs with $T_{\rm eff}\approx25,000$ K in binaries of 60-90 minutes with WDs. Though not yet in contact, these binaries would eventually come into contact through gravitational wave emission and explode as a sub-luminous thermonuclear supernova or evolve into a massive single WD., Comment: accepted for publication in ApJL, 11 pages, 5 figures and 1 tables. arXiv admin note: text overlap with arXiv:2002.01485

We report the discovery of the first short period binary in which a hot subdwarf star (sdOB) fills its Roche lobe and started mass transfer to its companion. The object was discovered as part of a dedicated high-cadence survey of the Galactic Plane named the Zwicky Transient Facility and exhibits a period of $P_{\rm orb}=39.3401(1)$ min, making it the most compact hot subdwarf binary currently known. Spectroscopic observations are consistent with an intermediate He-sdOB star with an effective temperature of $T_{\rm eff}=42,400\pm300$ K and a surface gravity of $\log(g)=5.77\pm0.05$. A high-signal-to noise GTC+HiPERCAM light curve is dominated by the ellipsoidal deformation of the sdOB star and an eclipse of the sdOB by an accretion disk. We infer a low-mass hot subdwarf donor with a mass $M_{\rm sdOB}=0.337\pm0.015$ M$_\odot$ and a white dwarf accretor with a mass $M_{\rm WD}=0.545\pm0.020$ M$_\odot$. Theoretical binary modeling indicates the hot subdwarf formed during a common envelope phase when a $2.5-2.8$ M$_\odot$ star lost its envelope when crossing the Hertzsprung Gap. To match its current $P_{\rm orb}$, $T_{\rm eff}$, $\log(g)$, and masses, we estimate a post-common envelope period of $P_{\rm orb}\approx150$ min, and find the sdOB star is currently undergoing hydrogen shell burning. We estimate that the hot subdwarf will become a white dwarf with a thick helium layer of $\approx0.1$ M$_\odot$ and will merge with its carbon/oxygen white dwarf companion after $\approx17$ Myr and presumably explode as a thermonuclear supernova or form an R CrB star., Comment: accepted for publication in ApJ, 21 pages, 12 figures and 4 tables

Some binary systems composed of a white dwarf (WD) and a hot subdwarf (sdB) helium star will make contact within the helium burning lifetime of the sdB star. The accreted helium on the WD inevitably undergoes a thermonuclear instability, causing a detonation that is expected to transition into the WD core and lead to a thermonuclear supernova while the donor orbits nearby with high velocity. Motivated by the recent discovery of fast-moving objects that occupy unusual locations on the HR diagram, we explore the impact of the thermonuclear supernovae on the donors in this specific double detonation scenario. We use MESA to model the binary up to the moment of detonation, then 3D Athena++ to model the hydrodynamic interaction of the supernova ejecta with the donor star, calculating the amount of mass that is stripped and the entropy deposited in the deep stellar interior by the strong shock that traverses it. We show that these donor remnants are ejected with velocities primarily set by their orbital speeds: $700-900\ \rm km\ s^{-1}$. We model the long-term thermal evolution of remnants by introducing the shock entropy into MESA models. In response to this entropy change, donor remnants expand and brighten for timescales ranging from $10^6-10^8$ years, giving ample time for these runaway stars to be observed in their inflated state before they leave the galaxy. Even after surface layers are stripped, some donors retain enough mass to resume core helium burning and further delay fading for more than $10^8$ years., Comment: 12 pages, 8 figures, 3 tables, accepted for publication in ApJ. Movies for figures 1 and 2 available in source

Using high-cadence observations from the Zwicky Transient Facility at low Galactic latitudes, we have discovered a new class of pulsating, hot, compact stars. We have found four candidates, exhibiting blue colors ($g-r\leq-0.1$ mag), pulsation amplitudes of $>5\%$, and pulsation periods of $200 - 475$ sec. Fourier transforms of the lightcurves show only one dominant frequency. Phase-resolved spectroscopy for three objects reveals significant radial velocity, T$_{\rm eff}$ and log(g) variations over the pulsation cycle, consistent with large amplitude radial oscillations. The mean T$_{\rm eff}$ and log(g) for these stars are consistent with hot subdwarf B (sdB) effective temperatures and surface gravities. We calculate evolutionary tracks using MESA and adiabatic pulsations using GYRE for low-mass helium-core pre-white dwarfs and low mass helium-burning stars. Comparison of low-order radial oscillation mode periods with the observed pulsation periods show better agreement with the pre-white dwarf models. Therefore, we suggest that these new pulsators and Blue Large-Amplitude Pulsators (BLAPs) could be members of the same class of pulsators, composed of young $\approx0.25-0.35$ M$_\odot$ helium-core pre-white dwarfs., Comment: accepted for publication in ApJL, 9 pages, 3 figures and 2 tables

Many isolated white dwarfs (WDs) show spectral evidence of atmospheric metal pollution. Since heavy element sedimentation timescales are short, this most likely indicates ongoing accretion. Accreted metals encounter a variety of mixing processes at the WD surface: convection, gravitational sedimentation, overshoot, and thermohaline instability. We present MESA WD models that explore each of these processes and their implications for inferred accretion rates. We provide diffusion timescales for many individual metals, and we quantify the regimes in which thermohaline mixing dominates over gravitational sedimentation in setting the effective settling rate of the heavy elements. We build upon and confirm earlier work finding that accretion rates as high as $10^{13} \, \rm g \, s^{-1}$ are needed to explain observed pollution in DA WDs for $T_{\rm eff} > 15,000 \, \rm K$, and we provide tabulated results from our models that enable accretion rate inferences from observations of polluted DA WDs. If these rates are representative of young WDs, we estimate that the total mass of planetesimal material accreted over a WD lifetime may be as high as $10^{28} \, \rm g$, though this estimate is susceptible to potential selection biases and uncertainties about the nature of disk processes that supply accretion to the WD surface. We also find that polluted DB WDs experience much less thermohaline mixing than DA WDs, and we do not expect thermohaline instability to be active for polluted DB WDs with $T_{\rm eff} < 18,000 \, \rm K$., Comment: 15 pages, 11 figures, 3 tables, Accepted for publication in ApJ

KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $\mu$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $\mu$Hz. A mean frequency difference of 30 $\mu$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction., Comment: 7 pages, 5 figures; Accepted for publication in ApJ Letters

Many isolated, old white dwarfs (WDs) show surprising evidence of metals in their photospheres. Given that the timescale for gravitational sedimentation is astronomically short, this is taken as evidence for ongoing accretion, likely of tidally disrupted planetesimals. The rate of such accretion, $\dot M_{\rm acc}$, is important to constrain, and most modeling of this process relies on assuming an equilibrium between diffusive sedimentation and metal accretion supplied to the WD's surface convective envelope. Building on earlier work of Deal and collaborators, we show that high $\dot M_{\rm acc}$ models with only diffusive sedimentation are unstable to thermohaline mixing and that models which account for the enhanced mixing from the active thermohaline instability require larger accretion rates, sometimes reaching $\dot M_{\rm acc} \approx 10^{13} \, \rm g \, s^{-1}$ to explain observed Calcium abundances. We present results from a grid of MESA models that include both diffusion and thermohaline mixing. These results demonstrate that both mechanisms are essential for understanding metal pollution across the range of polluted WDs with hydrogen atmospheres. Another consequence of active thermohaline mixing is that the observed metal abundance ratios are identical to accreted material., Comment: 7 pages, 4 figures; Accepted for publication in ApJ Letters

Most massive stars end their lives as Red Supergiants (RSGs), a short-lived evolution phase when they are known to pulsate with varying amplitudes. The RSG period-luminosity (PL) relation has been measured in the Milky Way, the Magellanic Clouds and M33 for about 120 stars in total. Using over 1500 epochs of R-band monitoring from the Palomar Transient Factory (PTF) survey over a five-year period, we study the variability of 255 spectroscopically cataloged RSGs in M31. We find that all RGSs brighter than M_K~ -10 mag (log(L/L_sun)>4.8) are variable at dm_R>0.05 mag. Our period analysis finds 63 with significant pulsation periods. Using the periods found and the known values of M_K for these stars, we derive the RSG PL relation in M31 and show that it is consistent with those derived earlier in other galaxies of different metallicities. We also detect, for the first time, a sequence of likely first-overtone pulsations. Comparison to stellar evolution models from MESA confirms the first overtone hypothesis and indicates that the variable stars in this sample have 12 M_sun

We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact., Comment: 64 pages, 61 figures; Accepted to AAS Journals

White dwarfs (WDs) that accrete helium at rates $\sim 10^{-8} \, M_\odot \, \rm yr^{-1}$, such as those in close binaries with sdB stars, can accumulate large ($\gtrsim 0.1 \, M_\odot$) helium envelopes which are likely to detonate. We perform binary stellar evolution calculations of sdB+WD binary systems with MESA, incorporating the important reaction chain $^{14}{\rm N}(e^-, \nu) {^{14}{\rm C}} ( \alpha, \gamma) {^{18}{\rm O}}$ (NCO), including a recent measurement for the ${^{14}{\rm C}} ( \alpha, \gamma) {^{18}{\rm O}}$ rate. In large accreted helium shells, the NCO reaction chain leads to ignitions at the dense base of the freshly accreted envelope, in contrast to $3\alpha$ ignitions which occur away from the base of the shell. In addition, at these accretion rates, the shells accumulate on a timescale comparable to their thermal time, leading to an enhanced sensitivity of the outcome on the accretion rate history. Hence, time dependent accretion rates from binary stellar evolution are necessary to determine the helium layer mass at ignition. We model the observed sdB+WD system ${\rm CD}\,-30^\circ 11223$ and find that the inclusion of these effects predicts ignition of a $0.153 \, M_\odot$ helium shell, nearly a factor of two larger than previous predictions. A shell with this mass will ignite dynamically, a necessary condition for a helium shell detonation., Comment: 11 pages, 9 figures; Accepted for publication in ApJ

We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development., Comment: Published in the ApJ Supplement Series; updated to account for published erratum