Showing 1–38 of 38 results for author: Bauer, B

Supernova Shocks Cannot Explain the Inflated State of Hypervelocity Runaways from White Dwarf Binaries

Authors: Aakash Bhat, Evan B. Bauer, Rüdiger Pakmor, Ken J. Shen, Ilaria Caiazzo, Abinaya Swaruba Rajamuthukumar, Kareem El-Badry, Wolfgang E. Kerzendorf

Abstract: 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… ▽ More 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. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: Submitted to A\&A. 15 pages, 17 figures

An Expanded Set of Los Alamos OPLIB Tables in MESA: Type-1 Rosseland-mean Opacities and Solar Models

Authors: Ebraheem Farag, Christopher J. Fontes, F. X. Timmes, Earl P. Bellinger, Joyce A. Guzik, Evan B. Bauer, Suzannah R. Wood, Katie Mussack, Peter Hakel, James Colgan, David P. Kilcrease, Manolo E. Sherrill, Tryston C. Raecke, Morgan T. Chidester

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 28 pages, 12 (13) figures. Accepted for Publication in ApJ

A Model for Eruptive Mass Loss in Massive Stars

Authors: Shelley J. Cheng, Jared A. Goldberg, Matteo Cantiello, Evan B. Bauer, Mathieu Renzo, Charlie Conroy

Abstract: 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 wit… ▽ More 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. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 16 pages, 9 figures, submitted to ApJ

OGLE-BLAP-009 -- A Case Study for the Properties and Evolution of Blue Large-Amplitude Pulsators

Authors: Corey W. Bradshaw, Matti Dorsch, Thomas Kupfer, Brad N. Barlow, Uli Heber, Evan B. Bauer, Lars Bildsten, Jan van Roestel

Abstract: 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 th… ▽ More 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}$. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS, 15 pages, 11 figures, 5 tables

Testing Modules for Experiments in Stellar Astrophysics (MESA)

Authors: William M. Wolf, Josiah Schwab, R. Farmer, Evan B. Bauer

Abstract: 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 compu… ▽ More 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. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 12 page, 7 figures, Accepted to ApJS

Modelling the AM CVn and Double Detonation Supernova Progenitor Binary System CD-30$^{\circ}$11223

Authors: Kunal Deshmukh, Evan B. Bauer, Thomas Kupfer, Matti Dorsch

Abstract: 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… ▽ More 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. △ Less

Submitted 25 October, 2023; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: 12 pages, 8 figures. Accepted for publication in MNRAS

Reinterpreting the Polluted White Dwarf SDSS J122859.93+104032.9 in Light of Thermohaline Mixing Models: More Polluting Material from a Larger Orbiting Solid Body

Authors: Arianna Dwomoh, Evan B. Bauer

Abstract: 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… ▽ More 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. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: 7 pages, 3 figures, 4 tables, accepted for publication in ApJ

Orbital decay in an accreting and eclipsing 13.7 minute orbital period binary with a luminous donor

Authors: Kevin B. Burdge, Kareem El-Badry, Saul Rappaport, Tin Long Sunny Wong, Evan B. Bauer, Lars Bildsten, Ilaria Caiazzo, Deepto Chakrabarty, Emma Chickles, Matthew J. Graham, Erin Kara, S. R. Kulkarni, Thomas R. Marsh, Melania Nynka, Thomas A. Prince, Robert A. Simcoe, Jan van Roestel, Zach Vanderbosch, Eric C. Bellm, Richard G. Dekany, Andrew J. Drake, George Helou, Frank J. Masci, Jennifer Milburn, Reed Riddle , et al. (2 additional authors not shown)

Abstract: 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 e… ▽ More 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σ$ 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. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 13 pages, 7 figures, 2 tables, submitted to ApJL

Carbon-Oxygen Phase Separation in MESA White Dwarf Models

Authors: Evan B. Bauer

Abstract: 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 cool… ▽ More 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. △ Less

Submitted 24 April, 2023; v1 submitted 17 March, 2023; originally announced March 2023.

Comments: 15 pages, 9 figures, accepted for publication in ApJ

Limiting the accretion disk light in two mass transferring hot subdwarf binaries

Authors: Kunal Deshmukh, Thomas Kupfer, Pasi Hakala, Evan B. Bauer, Andrei Berdyugin, Lars Bildsten, Thomas R. Marsh, Sandro Mereghetti, Vilppu Piirola

Abstract: 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 accre… ▽ More 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. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 10 pages, 8 figures. Accepted for publication in MNRAS

Modules for Experiments in Stellar Astrophysics (MESA): Time-Dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure

Authors: Adam S. Jermyn, Evan B. Bauer, Josiah Schwab, R. Farmer, Warrick H. Ball, Earl P. Bellinger, Aaron Dotter, Meridith Joyce, Pablo Marchant, Joey S. G. Mombarg, William M. Wolf, Tin Long Sunny Wong, Giulia C. Cinquegrana, Eoin Farrell, R. Smolec, Anne Thoul, Matteo Cantiello, Falk Herwig, Odette Toloza, Lars Bildsten, Richard H. D. Townsend, F. X. Timmes

Abstract: 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 MES… ▽ More 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. △ Less

Submitted 30 December, 2022; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 50 pages, 29 figures; Accepted to ApJS

Moosinesq Convection in the Cores of Moosive Stars

Authors: Evan H. Anders, Evan B. Bauer, Adam S. Jermyn, Samuel J. Van Kooten, Benjamin P. Brown, Eric W. Hester, Mindy Wilkinson, Jared A. Goldberg, Tania Varesano, Daniel Lecoanet

Abstract: 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… ▽ More 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. △ Less

Submitted 30 March, 2022; originally announced April 2022.

Comments: 7 pages, 3 figures, 1 moose, no logic in this place

Accurate Diffusion Coefficients for Dense White Dwarf Plasma Mixtures

Authors: M. E. Caplan, E. B. Bauer, I. F. Freeman

Abstract: Diffusion coefficients are essential microphysics input for modeling white dwarf evolution, as they impact phase separation at crystallization and sedimentary heat sources. Present schemes for computing diffusion coefficients are accurate at weak coupling ($Γ\ll 1$), but they have errors as large as a factor of two in the strongly coupled liquid regime ($1 \lesssim Γ\lesssim 200$). With modern mol… ▽ More Diffusion coefficients are essential microphysics input for modeling white dwarf evolution, as they impact phase separation at crystallization and sedimentary heat sources. Present schemes for computing diffusion coefficients are accurate at weak coupling ($Γ\ll 1$), but they have errors as large as a factor of two in the strongly coupled liquid regime ($1 \lesssim Γ\lesssim 200$). With modern molecular dynamics codes it is possible to accurately determine diffusion coefficients in select systems with percent-level precision. In this work, we develop a theoretically motivated law for diffusion coefficients which works across the wide range of parameters typical for white dwarf interiors. We perform molecular dynamics simulations of pure systems and two mixtures that respectively model a typical-mass C/O white dwarf and a higher-mass O/Ne white dwarf, and resolve diffusion coefficients for several trace neutron-rich nuclides. We fit the model to the pure systems and propose a physically motivated generalization for mixtures. We show that this model is accurate to roughly 15% when compared to molecular dynamics for many individual elements under conditions typical of white dwarfs, and is straightforward to implement in stellar evolution codes. △ Less

Submitted 24 March, 2022; originally announced March 2022.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS Letters

The Forward Physics Facility at the High-Luminosity LHC

Authors: Jonathan L. Feng, Felix Kling, Mary Hall Reno, Juan Rojo, Dennis Soldin, Luis A. Anchordoqui, Jamie Boyd, Ahmed Ismail, Lucian Harland-Lang, Kevin J. Kelly, Vishvas Pandey, Sebastian Trojanowski, Yu-Dai Tsai, Jean-Marco Alameddine, Takeshi Araki, Akitaka Ariga, Tomoko Ariga, Kento Asai, Alessandro Bacchetta, Kincso Balazs, Alan J. Barr, Michele Battistin, Jianming Bian, Caterina Bertone, Weidong Bai , et al. (211 additional authors not shown)

Abstract: High energy collisions at the High-Luminosity Large Hadron Collider (LHC) produce a large number of particles along the beam collision axis, outside of the acceptance of existing LHC experiments. The proposed Forward Physics Facility (FPF), to be located several hundred meters from the ATLAS interaction point and shielded by concrete and rock, will host a suite of experiments to probe Standard Mod… ▽ More High energy collisions at the High-Luminosity Large Hadron Collider (LHC) produce a large number of particles along the beam collision axis, outside of the acceptance of existing LHC experiments. The proposed Forward Physics Facility (FPF), to be located several hundred meters from the ATLAS interaction point and shielded by concrete and rock, will host a suite of experiments to probe Standard Model (SM) processes and search for physics beyond the Standard Model (BSM). In this report, we review the status of the civil engineering plans and the experiments to explore the diverse physics signals that can be uniquely probed in the forward region. FPF experiments will be sensitive to a broad range of BSM physics through searches for new particle scattering or decay signatures and deviations from SM expectations in high statistics analyses with TeV neutrinos in this low-background environment. High statistics neutrino detection will also provide valuable data for fundamental topics in perturbative and non-perturbative QCD and in weak interactions. Experiments at the FPF will enable synergies between forward particle production at the LHC and astroparticle physics to be exploited. We report here on these physics topics, on infrastructure, detector, and simulation studies, and on future directions to realize the FPF's physics potential. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: 429 pages, contribution to Snowmass 2021

Report number: UCI-TR-2022-01, CERN-PBC-Notes-2022-001, FERMILAB-PUB-22-094-ND-SCD-T, INT-PUB-22-006, BONN-TH-2022-04

Masses of White Dwarf Binary Companions to Type Ia Supernovae Measured from Runaway Velocities

Authors: Evan B. Bauer, Vedant Chandra, Ken J. Shen, J. J. Hermes

Abstract: 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… ▽ More 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. △ Less

Submitted 14 December, 2021; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: 8 pages, 4 figures, accepted for publication in ApJL

Discovery of a double detonation thermonuclear supernova progenitor

Authors: Thomas Kupfer, Evan B. Bauer, Jan van Roestel, Eric C. Bellm, Lars Bildsten, Jim Fuller, Thomas A. Prince, Ulrich Heber, Stephan Geier, Matthew J. Green, Shrinivas R. Kulkarni, Steven Bloemen, Russ R. Laher, Ben Rusholme, David Schneider

Abstract: 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 we… ▽ More 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. △ Less

Submitted 6 January, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: accepted for publication in ApJL, 14 pages, 5 figures, 1 table

The SNIa Runaway LP 398-9: Detection of Circumstellar Material and Surface Rotation

Authors: Vedant Chandra, Hsiang-Chih Hwang, Nadia L. Zakamska, Simon Blouin, Andrew Swan, Thomas R. Marsh, Ken J. Shen, Boris T. Gänsicke, J. J. Hermes, Odelia Putterman, Evan B. Bauer, Evan Petrosky, Vikram S. Dhillon, Stuart P. Littlefair, Richard P. Ashley

Abstract: 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 r… ▽ More 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. △ Less

Submitted 13 April, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: Accepted to MNRAS

A 99-minute Double-lined White Dwarf Binary from SDSS-V

Authors: Vedant Chandra, Hsiang-Chih Hwang, Nadia L. Zakamska, Boris T. Gaensicke, J. J. Hermes, Axel Schwope, Carles Badenes, Gagik Tovmassian, Evan B. Bauer, Dan Maoz, Matthias R. Schreiber, Odette F. Toloza, Keith P. Inight, Hans-Walter Rix, Warren R. Brown

Abstract: 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 near… ▽ More 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. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Comments: Accepted to ApJ

The final fates of close hot subdwarf - white dwarf binaries: mergers involving He/C/O white dwarfs and the formation of unusual giant stars with C/O-dominated envelopes

Authors: Josiah Schwab, Evan B. Bauer

Abstract: 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… ▽ More 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. △ Less

Submitted 5 August, 2021; originally announced August 2021.

Comments: 7 pages, 4 figures; Accepted to ApJ

Phases of Mass Transfer from Hot Subdwarfs to White Dwarf Companions and Their Photometric Properties

Authors: Evan B. Bauer, Thomas Kupfer

Abstract: 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.… ▽ More 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. △ Less

Submitted 9 September, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: 16 pages, 8 figures, accepted for publication in ApJ

Horizontal spreading of planetary debris accreted by white dwarfs

Authors: Tim Cunningham, Pier-Emmanuel Tremblay, Evan B. Bauer, Odette Toloza, Elena Cukanovaite, Detlev Koester, Jay Farihi, Bernd Freytag, Boris T. Gänsicke, Hans-Günter Ludwig, Dimitri Veras

Abstract: 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 whic… ▽ More 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. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: Accepted for publication in MNRAS. 22 pages, 17 figures

On The Impact Of 22Ne On The Pulsation Periods Of Carbon-Oxygen White Dwarfs With Helium Dominated Atmospheres

Authors: Morgan T. Chidester, F. X. Timmes, Josiah Schwab, Richard H. D. Townsend, Ebraheem Farag, Anne Thoul, C. E. Fields, Evan B. Bauer, Michael H. Montgomery

Abstract: 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 wi… ▽ More 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 $σ_{rms}$ $\lesssim$ 0.3 s, for a fit precision of $σ_{rms}/ P$ $\lesssim$ 0.3 %. We find average relative period shifts of $Δ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 $Δ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. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 23 pages, 19 figures total. Accepted to the ApJ. Presented results at 2021 AAS January meeting

Multi-Gigayear White Dwarf Cooling Delays from Clustering-Enhanced Gravitational Sedimentation

Authors: Evan B. Bauer, Josiah Schwab, Lars Bildsten, Sihao Cheng

Abstract: 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… ▽ More 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. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: 23 pages, 14 figures, accepted for publication in ApJ

Cool, Luminous, and Highly Variable Stars in the Magellanic Clouds from ASAS-SN: Implications for Thorne-Żytkow Objects and Super-Asymptotic Giant Branch Stars

Authors: Anna J. G. O'Grady, Maria R. Drout, B. J. Shappee, Evan B. Bauer, Jim Fuller, C. S. Kochanek, T. Jayasinghe, B. M. Gaensler, K. Z. Stanek, Thomas W. -S. Holoien, J. L. Prieto, Todd A. Thompson

Abstract: 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-Żytkow Object… ▽ More 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-Żytkow Object (TŻO) 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 $Δ$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ŻOs 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. △ Less

Submitted 14 August, 2020; originally announced August 2020.

Comments: 29 pages, 14 figures, 4 tables, plus Appendix. Accepted for publication in the Astrophysical Journal

A new class of Roche lobe-filling hot subdwarf binaries

Authors: Thomas Kupfer, Evan B. Bauer, Kevin B. Burdge, Jan van Roestel, Eric C. Bellm, Jim Fuller, JJ Hermes, Thomas R. Marsh, Lars Bildsten, Shrinivas R. Kulkarni, E. S. Phinney, Thomas A. Prince, Paula Szkody, Yuhan Yao, Andreas Irrgang, Ulrich Heber, David Schneider, Vik S. Dhillon, Gabriel Murawski, Andrew J. Drake, Dmitry A. Duev, Michael Feeney, Matthew J. Graham, Russ R. Laher, S. P. Littlefair , et al. (8 additional authors not shown)

Abstract: 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 lig… ▽ More 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. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Comments: accepted for publication in ApJL, 11 pages, 5 figures and 1 tables. arXiv admin note: text overlap with arXiv:2002.01485

Initial results from the New Horizons exploration of 2014 MU69, a small Kuiper Belt Object

Authors: S. A. Stern, H. A. Weaver, J. R. Spencer, C. B. Olkin, G. R. Gladstone, W. M. Grundy, J. M. Moore, D. P. Cruikshank, H. A. Elliott, W. B. McKinnon, J. Wm. Parker, A. J. Verbiscer, L. A. Young, D. A. Aguilar, J. M. Albers, T. Andert, J. P. Andrews, F. Bagenal, M. E. Banks, B. A. Bauer, J. A. Bauman, K. E. Bechtold, C. B. Beddingfield, N. Behrooz, K. B. Beisser , et al. (180 additional authors not shown)

Abstract: The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a fl… ▽ More The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes. △ Less

Submitted 2 April, 2020; originally announced April 2020.

Comments: 43 pages, 8 figure

Journal ref: Science 364, eaaw9771 (2019)

Intensity Mapping as a Probe of Axion Dark Matter

Authors: Jurek B. Bauer, David J. E. Marsh, Renée Hložek, Hamsa Padmanabhan, Alex Laguë

Abstract: Intensity mapping (IM) of spectral lines has the potential to revolutionize cosmology by increasing the total number of observed modes by several orders of magnitude compared to the cosmic microwave background (CMB) anisotropies. In this paper, we consider IM of neutral hydrogen (HI) in the redshift range $0 \lesssim z \lesssim 3$ employing a halo model approach where HI is assumed to follow the d… ▽ More Intensity mapping (IM) of spectral lines has the potential to revolutionize cosmology by increasing the total number of observed modes by several orders of magnitude compared to the cosmic microwave background (CMB) anisotropies. In this paper, we consider IM of neutral hydrogen (HI) in the redshift range $0 \lesssim z \lesssim 3$ employing a halo model approach where HI is assumed to follow the distribution of dark matter (DM) halos. If a portion of the DM is composed of ultralight axions then the abundance of halos is changed compared to cold dark matter below the axion Jeans mass. With fixed total HI density, $Ω_{\rm HI}$, assumed to reside entirely in halos, this effect introduces a scale-independent increase in the HI power spectrum on scales above the axion Jeans scale, which our model predicts consistent with N-body simulations. Lighter axions introduce a scale-dependent feature even on linear scales due to its suppression of the matter power spectrum near the Jeans scale. We use the Fisher matrix formalism to forecast the ability of future HI surveys to constrain the axion fraction of DM and marginalize over astrophysical and model uncertainties. We find that a HIRAX-like survey is a very reliable IM survey configuration, being affected minimally by uncertainties due to non-linear scales, while the SKA1MID configuration is the most constraining as it is sensitive to non-linear scales. Including non-linear scales and combining a SKA1MID-like IM survey with the Simons Observatory CMB, the benchmark "fuzzy DM" model with $m_a = 10^{-22}\text{ eV}$ can be constrained at the 10% level. For lighter ULAs this limit improves below 1%, and allows the possibility to test the connection between axion models and the grand unification scale across a wide range of masses. △ Less

Submitted 19 December, 2020; v1 submitted 21 March, 2020; originally announced March 2020.

Comments: 16 pages, 14 figures, 2 tables, code to reproduce our results is publicly available at https://github.com/JurekBauer/axion21cmIM.git

Journal ref: MNRAS 500.3, 3162 (2021)

The first ultracompact Roche lobe-filling hot subdwarf binary

Authors: Thomas Kupfer, Evan B. Bauer, Thomas R. Marsh, Jan van Roestel, Eric C. Bellm, Kevin B. Burdge, Michael W. Coughlin, Jim Fuller, JJ Hermes, Lars Bildsten, Shrinivas R. Kulkarni, Thomas A. Prince, Paula Szkody, Vik S. Dhillon, Gabriel Murawski, Rick Burruss, Richard Dekany, Alex Delacroix, Andrew J. Drake, Dmitry A. Duev, Michael Feeney, Matthew J. Graham, David L. Kaplan, Russ R. Laher, S. P. Littlefair , et al. (7 additional authors not shown)

Abstract: 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 kno… ▽ More 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. △ Less

Submitted 4 February, 2020; originally announced February 2020.

Comments: accepted for publication in ApJ, 21 pages, 12 figures and 4 tables

Remnants of Subdwarf Helium Donor Stars Ejected from Close Binaries with Thermonuclear Supernovae

Authors: Evan B. Bauer, Christopher J. White, Lars Bildsten

Abstract: 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. Motiv… ▽ More 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. △ Less

Submitted 6 October, 2019; v1 submitted 21 June, 2019; originally announced June 2019.

Comments: 12 pages, 8 figures, 3 tables, accepted for publication in ApJ. Movies for figures 1 and 2 available in source

A new class of pulsating hot subdwarfs

Authors: Thomas Kupfer, Evan B. Bauer, Kevin B. Burdge, Eric C. Bellm, Lars Bildsten, Jim Fuller, JJ Hermes, Shrinivas R. Kulkarni, Thomas A. Prince, Jan van Roestel, Richard Dekany, Dmitry A. Duev, Michael Feeney, Matteo Giomi, Matthew J. Graham, Stephen Kaye, Russ R. Laher, Frank J. Masci, Michael Porter, Reed Riddle, David L. Shupe, Roger M. Smith, Maayane T. Soumagnac, Paula Szkody, Charlotte Ward

Abstract: 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 spect… ▽ More 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. △ Less

Submitted 3 June, 2019; originally announced June 2019.

Comments: accepted for publication in ApJL, 9 pages, 3 figures and 2 tables

Polluted White Dwarfs: Mixing Regions and Diffusion Timescales

Authors: Evan B. Bauer, Lars Bildsten

Abstract: 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… ▽ More 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$. △ Less

Submitted 17 January, 2019; v1 submitted 22 December, 2018; originally announced December 2018.

Comments: 15 pages, 11 figures, 3 tables, Accepted for publication in ApJ

The Impact of White Dwarf Luminosity Profiles on Oscillation Frequencies

Authors: F. X. Timmes, Richard H. D. Townsend, Evan B. Bauer, Anne Thoul, C. E. Fields, William M. Wolf

Abstract: 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 proportion… ▽ More 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 $μ$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 $μ$Hz. A mean frequency difference of 30 $μ$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. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 7 pages, 5 figures; Accepted for publication in ApJ Letters

Increases to Inferred Rates of Planetesimal Accretion Due to Thermohaline Mixing in Metal Accreting White Dwarfs

Authors: Evan B. Bauer, Lars Bildsten

Abstract: 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… ▽ More 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. △ Less

Submitted 14 May, 2018; originally announced May 2018.

Comments: 7 pages, 4 figures; Accepted for publication in ApJ Letters

Variability of Red Supergiants in M31 from the Palomar Transient Factory

Authors: Monika D. Soraisam, Lars Bildsten, Maria R. Drout, Evan B. Bauer, Marat Gilfanov, Thomas Kupfer, Russ R. Laher, Frank Masci, Thomas A. Prince, Shrinivas R. Kulkarni, Thomas Matheson, Abhijit Saha

Abstract: 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,… ▽ More 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<M<24 M_sun. As these RSGs are the immediate progenitors to Type II-P core-collapse supernovae (SNe), we also explore the implication of their variability in the initial-mass estimates for SN progenitors based on archival images of the progenitors. We find that this effect is small compared to the present measurement errors. △ Less

Submitted 27 March, 2018; originally announced March 2018.

Comments: 17 pages, 10 figures

Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions

Authors: Bill Paxton, Josiah Schwab, Evan B. Bauer, Lars Bildsten, Sergei Blinnikov, Paul Duffell, R. Farmer, Jared A. Goldberg, Pablo Marchant, Elena Sorokina, Anne Thoul, Richard H. D. Townsend, F. X. Timmes

Abstract: 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… ▽ More 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. △ Less

Submitted 3 January, 2018; v1 submitted 23 October, 2017; originally announced October 2017.

Comments: 64 pages, 61 figures; Accepted to AAS Journals

Electron Captures on $^{14}{\rm N}$ as a Trigger for Helium Shell Detonations

Authors: Evan B. Bauer, Josiah Schwab, Lars Bildsten

Abstract: 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… ▽ More 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^-, ν) {^{14}{\rm C}} ( α, γ) {^{18}{\rm O}}$ (NCO), including a recent measurement for the ${^{14}{\rm C}} ( α, γ) {^{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α$ 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. △ Less

Submitted 17 July, 2017; originally announced July 2017.

Comments: 11 pages, 9 figures; Accepted for publication in ApJ

The Pluto system: Initial results from its exploration by New Horizons

Authors: S. A. Stern, F. Bagenal, K. Ennico, G. R. Gladstone, W. M. Grundy, W. B. McKinnon, J. M. Moore, C. B. Olkin, J. R. Spencer, H. A. Weaver, L. A. Young, T. Andert, J. Andrews, M. Banks, B. Bauer, J. Bauman, O. S. Barnouin, P. Bedini, K. Beisser, R. A. Beyer, S. Bhaskaran, R. P. Binzel, E. Birath, M. Bird, D. J. Bogan , et al. (126 additional authors not shown)

Abstract: The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly ext… ▽ More The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected. △ Less

Submitted 26 October, 2015; originally announced October 2015.

Comments: 8 pages - Initial Science paper from NASA's New Horizons Pluto Encounter

Journal ref: Science, vol 350, 292, 2015

Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions

Authors: Bill Paxton, Pablo Marchant, Josiah Schwab, Evan B. Bauer, Lars Bildsten, Matteo Cantiello, Luc Dessart, R. Farmer, H. Hu, N. Langer, R. H. D. Townsend, Dean M. Townsley, F. X. Timmes

Abstract: 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… ▽ More 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. △ Less

Submitted 8 January, 2017; v1 submitted 9 June, 2015; originally announced June 2015.

Comments: Published in the ApJ Supplement Series; updated to account for published erratum

Journal ref: ApJS, 220, 15 (2015)