Search

Stars in close binaries are tidally distorted, and this has a strong effect on their pulsation modes. We compute the mode frequencies and geometries of tidally distorted stars using perturbation theory, accounting for the effects of the Coriolis force and the coupling between different azimuthal orders $m$ of a multiplet induced by the tidal distortion. For tidally coupled dipole pressure modes, the tidal coupling dominates over the Coriolis force and the resulting pulsations are ``triaxial", with each of the three modes in a multiplet ``tidally tilted" to be aligned with the one of the three principal axes of the star. The observed amplitudes and phases of the dipole modes aligned orthogonal to the spin axis are modulated throughout the orbit, producing doublets in the power spectrum that are spaced by exactly twice the orbital frequency. Quadrupole modes have similar but slightly more complex behavior. This amplitude modulation allows for mode identification which can potentially enable detailed asteroseismic analyses of tidally tilted pulsators. Pressure modes should exhibit this behavior in stellar binaries close enough to be tidally synchronized, while gravity modes should remain aligned with the star's spin axis. We discuss applications to various types of pulsating stars, and the relationship between tidal tilting of pulsations and the ``single-sided" pulsations sometimes observed in very tidally distorted stars., Comment: Submitted to ApJ, comments welcome!

The origin of magnetic fields in white dwarfs (WDs) remains mysterious. Magnetic WDs are traditionally associated with field strengths $\gtrsim1\,\mathrm{MG}$, set by the sensitivity of typical spectroscopic magnetic field measurements. Informed by recent developments in red giant magnetoasteroseismology, we revisit the use of WD pulsations as a seismic magnetometer. WD pulsations primarily probe near-surface magnetic fields, whose effect on oscillation mode frequencies is to asymmetrize rotational multiplets and, if strong enough, suppress gravity-mode propagation altogether. The sensitivity of seismology to magnetic fields increases strongly with mode period and decreases quickly with the depth of the partial ionization-driven surface convective zone. We place upper limits for magnetic fields in $24$ pulsating WDs: $20$ hydrogen-atmosphere (DAV) and three helium-atmosphere (DBV) carbon-oxygen WDs, and one extremely low-mass (helium-core) pulsator. These bounds are typically $\sim1$-$10\,\mathrm{kG}$, although they can reach down to $\sim10$-$100\,\mathrm{G}$ for DAVs and helium-core WDs in which lower-frequency modes are excited. Seismic magnetometry may enable new insights into the formation and evolution of WD magnetism., Comment: 21 pages, 9 figures, submitted to ApJ

Failed supernovae (SNe), which are likely the main channel for forming stellar-mass black holes, are predicted to accompany mass ejections much weaker than typical core-collapse SNe. We conduct a grid of one-dimensional radiation hydrodynamical simulations to explore the emission of failed SNe from red supergiant progenitors, leveraging recent understanding of the weak explosion and the dense circumstellar matter (CSM) surrounding these stars. We find from these simulations and semi-analytical modeling that diffusion in the CSM prolongs the early emission powered by shock breakout/cooling. The early emission has peak luminosities of $\sim 10^7$-$10^8~L_\odot$ in optical and UV, and durations of days to weeks. The presence of dense CSM aids detection of the early bright peak from these events via near-future wide-field surveys such as Rubin Observatory, ULTRASAT and UVEX., Comment: 17 pages, 8 figures, submitted to ApJ. Model light curves available here: https://github.com/DTsuna/failedSNeLCs

The Transiting Exoplanet Survey Satellite (TESS) has enabled the discovery of numerous tidally tilted pulsators (TTPs), which are pulsating stars in close binaries where the presence of a tidal bulge has the effect of tilting the primary star's pulsation axes into the orbital plane. Recently, the modeling framework developed to analyze TTPs has been applied to the emerging class of tri-axial pulsators, which exhibit nonradial pulsations about three perpendicular axes. In this work, we report on the identification of the second-ever discovered tri-axial pulsator, with sixteen robustly-detected pulsation multiplets, of which fourteen are dipole doublets separated by 2$\nu_{\rm orb}$. We jointly fit the spectral energy distribution (SED) and TESS light curve of the star, and find that the primary is slightly evolved off the zero-age main sequence, while the less massive secondary still lies on the zero-age main sequence. Of the fourteen doublets, we associate eight with $Y_{10x}$ modes and six with novel $Y_{10y}$ modes. We exclude the existence of $Y_{11x}$ modes in this star and show that the observed pulsation modes must be $Y_{10y}$. We also present a toy model for the tri-axial pulsation framework in the context of this star. The techniques presented here can be utilized to rapidly analyze and confirm future tri-axial pulsator candidates., Comment: 19 pages, 11 figures, 4 tables. Accepted for publication in ApJ

Hot subdwarf B (sdB) stars are stripped helium-burning stars that are often found in close binaries, where they experience strong tidal interactions. The dissipation of tidally excited gravity waves alters their rotational evolution throughout the sdB lifetime. While many sdB binaries have well-measured rotational and orbital frequencies, there have been few theoretical efforts to accurately calculate the tidal torque produced by gravity waves. In this work, we directly calculate the tidal excitation of internal gravity waves in realistic sdB stellar models and integrate the coupled spin-orbit evolution of sdB binaries. We find that for canonical sdB ($M_\mathrm{sdB}=0.47\,M_\odot$) binaries, the transitional orbital period below which they could reach tidal synchronization in the sdB lifetime is $\sim \! 0.2\;\mathrm{days}$, with weak dependence on the companion masses. For low-mass sdBs ($M_\mathrm{sdB}=0.37\,M_\odot$) formed from more massive progenitor stars, the transitional orbital period becomes $\sim \! 0.15\;\mathrm{days}$. These values are very similar to the tidal synchronization boundary ($\sim \! 0.2\;\mathrm{days}$) evident from observations. We discuss the dependence of tidal torques on stellar radii, and we make predictions for the rapidly rotating white dwarfs formed from synchronized sdB binaries., Comment: Accepted to ApJ

Binary mass transfer can occur at high rates due to rapid expansion of the donor's envelope. In the case where mass transfer is unstable, the binary can rapidly shrink its orbit and lead to a merger. In this work we consider the appearance of the system preceding merger, specifically for the case of a low-mass ($\approx 2.5$-$3~M_\odot$) helium star with a neutron star (NS) companion. Modeling the mass transfer history as well as the wind launched by super-Eddington accretion onto the NS, we find that such systems can power slowly rising transients with timescales as long as years, and luminosities of $\sim 10^{40}$-$10^{41}$ erg s$^{-1}$ from optical to UV. The final explosion following the merger (or core-collapse of the helium star in some cases) leads to an interaction-powered transient with properties resembling Type Ibn supernovae (SNe), possibly with a bright early peak powered by shock cooling emission for merger-powered explosions. We apply our model to the Type Ibn SN 2023fyq, that displayed a long-term precursor activity from years before the terminal explosion., Comment: Accepted to the Open Journal of Astrophysics. A typographic error in equation 34 has been corrected, and correspondingly Figure 5 has been slightly updated (other results are unaffected). Light curve source code available here: https://github.com/DTsuna/merger_precursor.git

The mass loss mechanism of red supergiant stars is not well understood, even though it has crucial consequences for their stellar evolution and the appearance of supernovae that occur upon core-collapse. We argue that outgoing shock waves launched near the photosphere can support a dense chromosphere between the star's surface and the dust formation radius at several stellar radii. We derive analytic expressions for the time-averaged density profile of the chromosphere, and we use these to estimate mass loss rates due to winds launched by radiation pressure at the dust formation radius. These mass loss rates are similar to recent observations, possibly explaining the upward kink in mass loss rates of luminous red supergiants. Our models predict that low-mass red supergiants lose less mass than commonly assumed, while high-mass red supergiants lose more. The chromospheric mass of our models is $\sim$0.01 solar masses, most of which lies within a few stellar radii. This can help explain the early light curves and spectra of type-II P supernovae without requiring extreme pre-supernova mass loss. We discuss implications for stellar evolution, type II-P supernovae, SN 2023ixf, and Betelgeuse., Comment: Accepted to the Open Journal of Astrophysics

Mergers between helium white dwarfs and main-sequence stars are likely common, producing red giant-like remnants making up roughly a few percent of all low-mass ($\lesssim2M_\odot$) red giants. Through detailed modeling, we show that these merger remnants possess distinctive photometric, asteroseismic, and surface abundance signatures through which they may be identified. During hydrogen shell burning, merger remnants reach higher luminosities and possess pulsations which depart from the usual degenerate sequence on the asteroseismic $\Delta\nu$--$\Delta\Pi$ diagram for red giant branch stars. For sufficiently massive helium white dwarfs, merger remnants undergo especially violent helium flashes which can dredge up a large amount of core material (up to $\sim0.1M_\odot$) into the envelope. Such post-dredge-up remnants are more luminous than normal red clump stars, are surface carbon-, helium-, and possibly lithium-rich, and possess a wider range of asteroseismic g-mode period spacings and mixed-mode couplings. Recent asteroseismically determined low-mass ($\lesssim0.8M_\odot$) red clump stars may be core helium-burning remnants of mergers involving lower-mass helium white dwarfs., Comment: 24 pages, 16 figures, accepted to the Open Journal of Astrophysics!

We present SN 2023zaw $-$ a sub-luminous ($\mathrm{M_r} = -16.7$ mag) and rapidly-evolving supernova ($\mathrm{t_{1/2,r}} = 4.9$ days), with the lowest nickel mass ($\approx0.002$ $\mathrm{M_\odot}$) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He I and Ca NIR emission lines with velocities of $\sim10\ 000 - 12\ 000$ $\mathrm{km\ s^{-1}}$. The late-time spectra show prominent narrow He I emission lines at $\sim$1000$\ \mathrm{km\ s^{-1}}$, indicative of interaction with He-rich circumstellar material. SN 2023zaw is located in the spiral arm of a star-forming galaxy. We perform radiation-hydrodynamical and analytical modeling of the lightcurve by fitting with a combination of shock-cooling emission and nickel decay. The progenitor has a best-fit envelope mass of $\approx0.2$ $\mathrm{M_\odot}$ and an envelope radius of $\approx50$ $\mathrm{R_\odot}$. The extremely low nickel mass and low ejecta mass ($\approx0.5$ $\mathrm{M_\odot}$) suggest an ultra-stripped SN, which originates from a mass-losing low mass He-star (ZAMS mass $<$ 10 $\mathrm{M_\odot}$) in a close binary system. This is a channel to form double neutron star systems, whose merger is detectable with LIGO. SN 2023zaw underscores the existence of a previously undiscovered population of extremely low nickel mass ($< 0.005$ $\mathrm{M_\odot}$) stripped-envelope supernovae, which can be explored with deep and high-cadence transient surveys., Comment: Accepted for publication in ApJL

Tidal interactions between moons and planets can have major effects on the orbits, spins, and thermal evolution of the moons. In the Saturn system, tidal dissipation in the planet transfers angular momentum from Saturn to the moons, causing them to migrate outwards. The rate of migration is determined by the mechanism of dissipation within the planet, which is closely tied to the planet's uncertain structure. We review current knowledge of giant planet internal structure and evolution, which has improved thanks to data from the \textit{Juno} and \textit{Cassini} missions. We discuss general principles of tidal dissipation, describing both equilibrium and dynamical tides, and how dissipation can occur in a solid core or a fluid envelope. Finally, we discuss the possibility of resonance locking, whereby a moon can lock into resonance with a planetary oscillation mode, producing enhanced tidal migration relative to classical theories, and possibly explaining recent measurements of moon migration rates., Comment: Accepted for Space Science Reviews. Chapter in the book based on the ISSI workshop "New Vision of the Saturnian System in the Context of a Highly Dissipative Saturn" (9-13 May 2022)

Double white dwarf (WD) binaries are increasingly being discovered at short orbital periods where strong tidal effects and significant tidal heating signatures may occur. We assume the tidal potential of the companion excites outgoing gravity waves within the WD primary, the dissipation of which leads to an increase in the WD's surface temperature. We compute the excitation and dissipation of the waves in cooling WD models in evolving MESA binary simulations. Tidal heating is self-consistently computed and added to the models at every time step. As a binary inspirals to orbital periods less than $\sim$20 minutes, the WD's behavior changes from cooling to heating, with temperature enhancements that can exceed 10,000 K compared with non-tidally heated models. We compare a grid of tidally heated WD models to observed short-period systems with hot WD primaries. While tidal heating affects their $T_{\rm eff}$, it is likely not the dominant luminosity. Instead these WDs are probably intrinsically young and hot, implying the binaries formed at short orbital periods. The binaries are consistent with undergoing common envelope evolution with a somewhat low efficiency $\alpha_{\rm CE}$. We delineate the parameter space where the traveling wave assumption is most valid, noting that it breaks down for WDs that cool sufficiently, where standing waves may instead be formed., Comment: 14 pages, 13 figures. Accepted to ApJ

A fraction of core-collapse supernovae (SNe) with signs of interaction with a dense circumstellar matter are preceded by bright precursor emission. While the precursors are likely caused by a mass ejection before core-collapse, their mechanism to power energetic bursts, sometimes reaching $10^{48}$--$10^{49}\ {\rm erg}$ that are larger than the binding energies of red supergiant envelopes, is still under debate. Remarkably, such a huge energy-deposition should result in an almost complete envelope ejection and hence a strong sign of interaction, but the observed SNe with precursors show in fact typical properties among the interacting SNe. More generally, the observed luminosity of $10^{40-42}\,\rm erg\,s^{-1}$ is shown to be challenging for a single SN progenitor. To resolve these tensions, we propose a scenario where the progenitor is in a binary system with a compact object (CO), and an outburst from the star leads to a super-Eddington accretion onto the CO. We show that for sufficiently short separations, outbursts with moderate initial kinetic energies of $10^{46}$--$10^{47}$ erg can be energized by the accreting CO so that their radiative output can be consistent with the observed precursors. We discuss the implications of our model in relation to CO binaries detectable with Gaia and gravitational wave detectors., Comment: Accepted by ApJ; 18 pages, 5 Figures (+ appendix). Light curve source code available here: https://github.com/DTsuna/binary_precursor

Context. Slowly pulsating B (SPB) stars display multi-periodic variability in the gravito-inertial mode regime with indications of non-linear resonances between modes. Several have undergone asteroseismic modeling in the past few years to infer their internal properties in a linear setting. Rotation is typically included in the modeling by means of the traditional approximation of rotation (TAR). Aims. We aim to extend the set of tools available for asteroseismology, by describing time-independent (stationary) resonant non-linear coupling among three gravito-inertial modes within the TAR. Such coupling offers the opportunity to use mode amplitude ratios in the asteroseismic modeling process, instead of only relying on frequencies of linear eigenmodes. Methods. Following observational detections, we derive expressions for the resonant stationary non-linear coupling between three gravito-inertial modes in rotating stars. We assess selection rules and stability domains for stationary solutions and predict non-linear frequencies and amplitude ratio observables that can be compared with their observed counterparts. Results. The non-linear frequency shifts of stationary couplings are negligible compared to typical frequency errors derived from observations. The theoretically predicted amplitude ratios of combination frequencies match with some of their observational counterparts in the SPB targets. Other observed ratios could be linked to other saturation mechanisms, to interactions between different modes, or to different opacity gradients in the driving zone. Conclusions. Our non-linear mode coupling formalism can explain some of the stationary amplitude ratios of observed resonant mode couplings in single SPB stars monitored during 4 years by Kepler., Comment: After language revision; Accepted by A&A (26 pages, 6 figures, 6 tables, 7 appendices). Abstract abridged for arXiv submission

The orbital architectures of short-period exoplanet systems are shaped by tidal dissipation in their host stars. For low-mass M-dwarfs whose dynamical tidal response comprises a dense spectrum of inertial modes at low frequencies, resolving the frequency dependence of tidal dissipation is crucial to capturing the effect of tides on planetary orbits throughout the evolutionary stages of the host star. We use non-perturbative spectral methods to calculate the normal mode oscillations of a fully-convective M-dwarf modeled using realistic stellar profiles from MESA. We compute the dissipative tidal response composed of contributions from each mode as well as non-adiabatic coupling between the modes, which we find to be an essential component of the dissipative calculations. Using our results for dissipation, we then compute of the evolution of circular, coplanar planetary orbits under the influence of tides in the host star. We find that orbital migration driven by resonance locking affects the orbits of Earth-mass planets at orbital periods $P_{\rm orb} \lesssim 1.5$ day and of Jupiter-mass planets at $P_{\rm orb} \lesssim 2.5$ day. Due to resonantly-driven orbital decay and outward migration, we predict a dearth of small planets closer than $P_{\rm orb} \sim 1$ day and similarly sparse numbers of more massive planets out to $P_{\rm orb} \sim 3$ day., Comment: Accepted in ApJ

Post-common-envelope binaries (PCEBs) containing a white dwarf (WD) and a main-sequence (MS) star can constrain the physics of common envelope evolution and calibrate binary evolution models. Most PCEBs studied to date have short orbital periods ($P_{\rm orb} \lesssim 1\,$d), implying relatively inefficient harnessing of binaries' orbital energy for envelope expulsion. Here, we present follow-up observations of five binaries from {\it Gaia} DR3 containing solar-type MS stars and probable ultramassive WDs ($M\gtrsim 1.2\,M_{\odot}$) with significantly wider orbits than previously known PCEBs, $P_{\rm orb} = 18-49\,$d. The WD masses are much higher than expected for systems formed via stable mass transfer at these periods, and their near-circular orbits suggest partial tidal circularization when the WD progenitors were giants. These properties strongly suggest that the binaries are PCEBs. Forming PCEBs at such wide separations requires highly efficient envelope ejection, and we find that the observed periods can only be explained if a significant fraction of the energy released when the envelope recombines goes into ejecting it. Our 1D stellar models including recombination energy confirm prior predictions that a wide range of PCEB orbital periods, extending up to months or years, can potentially result from Roche lobe overflow of a luminous AGB star. This evolutionary scenario may also explain the formation of several wide WD+MS binaries discovered via self-lensing, as well as a significant fraction of post-AGB binaries and barium stars., Comment: 21 pages, 10 figures, accepted to MNRAS

Recent asteroseismic measurements have revealed a small population of stars in close binaries, containing primaries with extremely slow rotation rates. Such stars defy the standard expectation of tidal synchronization in such systems, but they can potentially be explained if they are trapped in a spin-orbit equilibrium known as Cassini state 2 (CS2). This state is maintained by orbital precession due to an outer tertiary star, and it typically results in a very sub-synchronous rotation rate and high degree of spin-orbit misalignment. We examine how CS2 is affected by magnetic braking and different types of tidal dissipation. Magnetic braking results in a slower equilibrium rotation rate, while tidal dissipation via gravity waves can result in a slightly higher rotation rate than predicted by equilibrium tidal theory, and dissipation via inertial waves can result in much slower rotation rates. For seven binary systems with slowly rotating primaries, we predict the location of the outer tertiary predicted by the CS2 theory. In five of these systems, a tertiary companion has already been detected, although it closer than expected in three of these, potentially indicating tidal dissipation via inertial waves. We also identify a few new candidate systems among a population of eclipsing binaries with rotation measurements via spot modulation., Comment: 14 pages, 13 figures

Stars in short-period binaries typically have spins that are aligned and synchronized with the orbit of their companion. In triple systems, however, the combination of spin and orbital precession can cause the star's rotation to evolve to a highly misaligned and sub-synchronous equilibrium known as a Cassini state. We identify a population of recently discovered stars that exhibit these characteristics and which are already known to have tertiary companions. These third bodies have a suitable orbital period to allow the inner binary to evolve into the sub-synchronous Cassini state, which we confirm with orbital evolution models. We also compute the expected stellar obliquity and spin period, showing that the observed rotation rates are often slower than expected from equilibrium tidal models. However, we show that tidal dissipation via inertial waves can alter the expected spin-orbit misalignment angle and rotation rate, potentially creating the very slow rotation rates in some systems. Finally, we show how additional discoveries of such systems can be used to constrain the tidal physics and orbital evolution histories of stellar systems., Comment: Revised after submission to MNRAS

When main-sequence stars expand into red giants, they are expected to engulf close-in planets. Until now, the absence of planets with short orbital periods around post-expansion, core-helium-burning red giants has been interpreted as evidence that short-period planets around Sun-like stars do not survive the giant expansion phase of their host stars. Here we present the discovery that the giant planet 8 Ursae Minoris b orbits a core-helium-burning red giant. At a distance of only 0.5 au from its host star, the planet would have been engulfed by its host star, which is predicted by standard single-star evolution to have previously expanded to a radius of 0.7 au. Given the brief lifetime of helium-burning giants, the nearly circular orbit of the planet is challenging to reconcile with scenarios in which the planet survives by having a distant orbit initially. Instead, the planet may have avoided engulfment through a stellar merger that either altered the evolution of the host star or produced 8 Ursae Minoris b as a second-generation planet. This system shows that core-helium-burning red giants can harbour close planets and provides evidence for the role of non-canonical stellar evolution in the extended survival of late-stage exoplanetary systems., Comment: Published in Nature on 28 June 2023. In press

Eruptive mass loss of massive stars prior to supernova (SN) explosion is key to understanding their evolution and end fate. An observational signature of pre-SN mass loss is the detection of an early, short-lived peak prior to the radioactive-powered peak in the lightcurve of the SN. This is usually attributed to the SN shock passing through an extended envelope or circumstellar medium (CSM). Such an early peak is common for double-peaked Type IIb SNe with an extended Hydrogen envelope but is uncommon for normal Type Ibc SNe with very compact progenitors. In this paper, we systematically study a sample of 14 double-peaked Type Ibc SNe out of 475 Type Ibc SNe detected by the Zwicky Transient Facility. The rate of these events is ~ 3-9 % of Type Ibc SNe. A strong correlation is seen between the peak brightness of the first and the second peak. We perform a holistic analysis of this sample's photometric and spectroscopic properties. We find that six SNe have ejecta mass less than 1.5 Msun. Based on the nebular spectra and lightcurve properties, we estimate that the progenitor masses for these are less than ~ 12 Msun. The rest have an ejecta mass > 2.4 Msun and a higher progenitor mass. This sample suggests that the SNe with low progenitor masses undergo late-time binary mass transfer. Meanwhile, the SNe with higher progenitor masses are consistent with wave-driven mass loss or pulsation-pair instability-driven mass loss simulations., Comment: Accepted for publication in ApJ

We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities $\gtrsim 1300\,\rm km\,s^{-1}$. These are most likely the surviving donors from double-degenerate binaries in which the other WD exploded. The other two objects have lower minimum velocities, $\gtrsim 600\,\rm km\,s^{-1}$, and may have formed through a different mechanism, such as pure deflagration of a WD in a Type Iax supernova. The four fastest stars are hotter and smaller than the previously known "D$^6$ stars," with effective temperatures ranging from $\sim$20,000 to $\sim$130,000 K and radii of $\sim 0.02-0.10\,R_{\odot}$. Three of these have carbon-dominated atmospheres, and one has a helium-dominated atmosphere. Two stars have RVs of $-1694$ and $-2285\rm \,km\,s^{-1}$ -- the fastest systemic stellar RVs ever measured. Their inferred birth velocities, $\sim 2200-2500\,\rm km\,s^{-1}$, imply that both WDs in the progenitor binary had masses $>1.0\,M_{\odot}$. The high observed velocities suggest that a dominant fraction of the observed hypervelocity WD population comes from double-degenerate binaries whose total mass significantly exceeds the Chandrasekhar limit. However, the two nearest and faintest D$^6$ stars have the lowest velocities and masses, suggesting that observational selection effects favor rarer, higher-mass stars. A significant population of fainter low-mass runaways may still await discovery. We infer a birth rate of D$^6$ stars that is consistent with the SN Ia rate. The birth rate is poorly constrained, however, because the luminosities and lifetimes of $\rm D^6$ stars are uncertain., Comment: 26 pages, 17 figures. Accepted to OJA

The detections of four apparently young radio pulsars in the Milky Way globular clusters are difficult to reconcile with standard neutron star formation scenarios associated with massive star evolution. Here we discuss formation of these young pulsars through white dwarf mergers in dynamically-old clusters that have undergone core collapse. Based on observed properties of magnetic white dwarfs, we argue neutron stars formed via white dwarf merger are born with spin periods of roughly $10-100\,$ms and magnetic fields of roughly $10^{11}-10^{13}\,$G. As these neutron stars spin down via magnetic dipole radiation, they naturally reproduce the four observed young pulsars in the Milky Way clusters. Rates inferred from $N$-body cluster simulations as well as the binarity, host cluster properties, and cluster offsets observed for these young pulsars hint further at a white dwarf merger origin. These young pulsars may be descendants of neutron stars capable of powering fast radio bursts analogous to the bursts observed recently in a globular cluster in M81., Comment: 5 pages, 2 figures, 1 table. Accepted for publication in MNRAS

Gravitational wave observations indicate the existence of merging black holes (BHs) with high spin ($a\gtrsim0.3$), whose formation pathways are still an open question. A possible way to form those binaries is through the tidal spin-up of a Wolf-Rayet (WR) star by its BH companion. In this work, we investigate this scenario by directly calculating the tidal excitation of oscillation modes in WR star models, determining the tidal spin-up rate, and integrating the coupled spin-orbit evolution for WR-BH binaries. We find that for short-period orbits and massive WR stars, the tidal interaction is mostly contributed by standing gravity modes, in contrast to Zahn's model of travelling waves which is frequently assumed in the literature. The standing modes are less efficiently damped than traveling waves, meaning that prior estimates of tidal spin-up may be overestimated. We show that tidal synchronization is rarely reached in WR-BH binaries, and the resulting BH spins have $a \lesssim 0.4$ for all but the shortest period ($P_{\rm orb} \! \lesssim 0.5 \, {\rm d}$) binaries. Tidal spin-up in lower-mass systems is more efficient, providing an anti-correlation between the mass and spin of the BHs, which could be tested in future gravitational wave data. Nonlinear damping processes are poorly understood but may allow for more efficient tidal spin-up. We also discuss a new class of gravito-thermal modes that appear in our calculations., Comment: Accepted to ApJ, corrected by the erratum published at https://iopscience.iop.org/article/10.3847/1538-4357/ad372c/meta

Normal mode seismology is a promising means of measuring rotation in gas giant interiors, and ring seismology presents a singular opportunity to do so at Saturn. We calculate Saturn's normal modes of oscillation and zonal gravity field, using nonperturbative methods for normal modes in the rigidly rotating approximation, and perturbative methods for the shifts that Saturn's deep winds induce in the mode frequencies and zonal gravity harmonics. The latter are calculated by solving the thermo-gravitational wind equation in an oblate geometry. Comparing many such models to gravity data and the frequencies of ring patterns excited by Saturn normal modes, we use statistical methods to estimate that Saturn's cloud-level winds extend inward along cylinders before decaying at a depth 0.125-0.138 times Saturn's equatorial radius, or 7,530-8,320 km, consistent with analyses of Cassini gravity and magnetic field data. The seismology is especially useful for pinning down Saturn's poorly constrained deep rotation period, which we estimate at 634.7 min (median) with a 5/95% quantile range 633.8-635.5 min. Outstanding residuals in mode frequencies at low angular degree suggest a more complicated deep interior than has been considered to date. Smaller but still significant residuals at high angular degree also show that our picture for the thermal, composition, and/or rotation profile in Saturn's envelope is not yet complete., Comment: Accepted to PSJ

Strong magnetic fields in the cores of stars are expected to significantly modify the behavior of gravity waves: this is likely the origin of suppressed dipole modes observed in many red giants. However, a detailed understanding of how such fields alter the spectrum and spatial structure of magnetogravity waves has been elusive. For a dipole field, we analytically characterize the horizontal eigenfunctions of magnetogravity modes, assuming that the wavevector is primarily radial. For axisymmetric modes ($m=0$), the magnetogravity wave eigenfunctions become Hough functions, and they have a radial turning point for sufficiently strong magnetic fields. For non-axisymmetric modes ($m\neq0$), the interaction between the discrete $g$ mode spectrum and a continuum of Alfv\'en waves produces nearly discontinuous features in the fluid displacements at critical latitudes associated with a singularity in the fluid equations. effects. We find that magnetogravity modes cannot propagate in regions with sufficiently strong magnetic fields, instead becoming evanescent. When encountering strong magnetic fields, ingoing gravity waves are likely refracted into outgoing slow magnetic waves. These outgoing waves approach infinite radial wavenumbers, which are likely to be damped efficiently. However, it may be possible for a small fraction of the wave power to escape the stellar core as pure Alfv\'en waves or magnetogravity waves confined to a very narrow equatorial band. The artificially sharp features in the WKB-separated solutions suggest the need for global mode solutions which include small terms neglected in our analysis., Comment: 22 pages, 10 figures, 1 table, comments welcome! (submitted to MNRAS)

Strong magnetic fields are observed in a substantial fraction of upper main sequence stars and white dwarfs. Many such stars are observed to exhibit photometric modulations as the magnetic poles rotate in and out of view, which could be a consequence of magnetic perturbations to the star's thermal structure. The magnetic pressure is typically larger than the gas pressure at the star's photosphere, but much smaller than the gas pressure in the star's interior, so the expected surface flux perturbations are not clear. We compute magnetically perturbed stellar structures of young $3 \, M_\odot$ stars that are in both hydrostatic and thermal equilibrium, and which contain both poloidal and toroidal components of a dipolar magnetic field as expected for stable fossil fields. This provides semi-analytical models of such fields in baroclinic stably stratified regions. The star's internal pressure, temperature, and flux perturbations can have a range of magnitudes, though we argue the most likely configurations exhibit flux perturbations much smaller than the ratio of surface magnetic pressure to surface gas pressure, but much larger than the ratio of surface magnetic pressure to central gas pressure. The magnetic pole is hotter than the equator in our models, but a cooler magnetic pole is possible depending on the magnetic field configuration. The expected flux variations for observed field strengths are $\delta L/L \! \lesssim \! 10^{-6}$, much smaller than those observed in magnetic stars, suggesting that observed perturbations stem from changes to the emergent spectrum rather than changes to the bolometric flux., Comment: Submitted to MNRAS, comments welcome!

Common envelope (CE) evolution, which is crucial in creating short period binaries and associated astrophysical events, can be constrained by reverse modeling of such binaries' formation histories. Through analysis of a sample of well-constrained white dwarf (WD) binaries with low-mass primaries (7 eclipsing double WDs, 2 non-eclipsing double WDs, 1 WD-brown dwarf), we estimate the CE energy efficiency $\alpha_{\rm{CE}}$ needed to unbind the hydrogen envelope. We use grids of He- and CO-core WD models to determine the masses and cooling ages that match each primary WD's radius and temperature. Assuming gravitational wave-driven orbital decay, we then calculate the associated ranges in post-CE orbital period. By mapping WD models to a grid of red giant progenitor stars, we determine the total envelope binding energies and possible orbital periods at the point CE evolution is initiated, thereby constraining $\alpha_{\rm CE}$. Assuming He-core WDs with progenitors of 0.9 - 2.0 $M_\odot$, we find $\alpha_{\rm CE} \! \sim \! 0.2-0.4$ is consistent with each system we model. Significantly higher values of $\alpha_{\rm{CE}}$ are required for higher mass progenitors and for CO-core WDs, so these scenarios are deemed unlikely. Our values are mostly consistent with previous studies of post-CE WD binaries, and they suggest a nearly constant and low envelope ejection efficiency for CE events that produce He-core WDs., Comment: 20 pages, 15 figures; Added references, corrected typos in discussion, updated figure resolution; Accepted by MNRAS

Planet engulfment can be inferred from enhancement of refractory elements in the photosphere of the engulfing star following accretion of rocky planetary material. Such refractory enrichments are subject to stellar interior mixing processes, namely thermohaline mixing induced by an inverse mean-molecular-weight gradient between the convective envelope and radiative core. Using MESA stellar models, we quantified the strength and duration of engulfment signatures following planet engulfment. We found that thermohaline mixing dominates during the first $\sim$5$-$45 Myr post-engulfment, weakening signatures by a factor of $\sim$2 before giving way to depletion via gravitational settling on longer timescales. Solar metallicity stars in the 0.5-1.2 $M_{\odot}$ mass range have observable signature timescales of $\sim$1 Myr$-$8 Gyr, depending on the engulfing star mass and amount of material engulfed. Early type stars exhibit larger initial refractory enhancements but more rapid depletion. Solar-like stars ($M$ = 0.9$-$1.1 $M_{\odot}$) maintain observable signatures ($>$0.05 dex) over timescales of $\sim$20 Myr$-$1.7 Gyr for nominal 10 $M_{\oplus}$ engulfment events, with longer-lived signatures occurring for low-metallicity and/or hotter stars (1 $M_{\odot}$, $\sim$2$-$3 Gyr). Engulfment events occurring well after the zero-age main sequence produce larger signals due to suppression of thermohaline mixing by gravitational settling of helium (1 $M_{\odot}$, $\sim$1.5 Gyr). These results indicate that it may be difficult to observe engulfment signatures in solar-like stars that are several Gyr old., Comment: 13 pages, 8 figures; accepted for publication in MNRAS

Many core collapse supernovae (SNe) with hydrogen-poor and low-mass ejecta, such as ultra-stripped SNe and type Ibn SNe, are observed to interact with dense circumstellar material (CSM). These events likely arise from the core-collapse of helium stars which have been heavily stripped by a binary companion and ejected significant mass during the last weeks to years of their lives. In helium star models run to days before core-collapse, we identify a range of helium core masses $\approx 2.5 -3 M_{\odot}$ whose envelopes expand substantially due to helium shell burning while the core undergoes neon and oxygen burning. When modeled in binary systems, the rapid expansion of these helium stars induces extremely high rates of late-stage mass transfer ($\dot{M} \gtrsim 10^{-2} M_\odot/{\rm yr}$) beginning weeks to decades before core-collapse. We consider two scenarios for producing CSM in these systems: either mass transfer remains stable and mass loss is driven from the system in the vicinity of the accreting companion, or mass transfer becomes unstable and causes a common envelope event (CEE) through which the helium envelope is unbound. The ensuing CSM properties are consistent with the CSM masses ($\sim 10^{-2}-1 M_\odot$) and radii ($\sim 10^{13}-10^{16} {\rm cm}$) inferred for ultra-stripped SNe and several type Ibn SNe. Furthermore, systems that undergo a CEE could produce short-period NS binaries that merge in less than 100 Myr.

In over a thousand known cataclysmic variables (CVs), where a white dwarf is accreting from a hydrogen-rich star, only a dozen have orbital periods below 75 minutes. One way to achieve these short periods requires the donor star to have undergone substantial nuclear evolution prior to interacting with the white dwarf, and it is expected that these objects will transition to helium accretion. These transitional CVs have been proposed as progenitors of helium CVs. However, no known transitional CV is expected to reach an orbital period short enough to account for most of the helium CV population, leaving the role of this evolutionary pathway unclear. Here we report observations of ZTF J1813+4251, a 51-minute orbital period, fully eclipsing binary system consisting of a star with a temperature comparable to that of the Sun but a density 100 times greater due to its helium-rich composition, accreting onto a white dwarf. Phase-resolved spectra, multi-band light curves and the broadband spectral energy distribution allow us to obtain precise and robust constraints on the masses, radii and temperatures of both components. Evolutionary modeling shows that ZTF J1813+4251 is destined to become a helium CV binary, reaching an orbital period under 20 minutes, rendering ZTF J1813+4251 a previously missing link between helium CV binaries and hydrogen-rich CVs., Comment: 48 Pages, 12 figures, 2 tables, Published online by Nature on Oct 5, 2022

The Tayler instability is an important but poorly studied magnetohydrodynamic instability that likely operates in stellar interiors. The nonlinear saturation of the Tayler instability is poorly understood and has crucial consequences for dynamo action and angular momentum transport in radiative regions of stars. We perform three-dimensional MHD simulations of the Tayler instability in a cylindrical geometry, including strong buoyancy and Coriolis forces as appropriate for its operation in realistic rotating stars. The linear growth of the instability is characterized by a predominantly $m=1$ oscillation with growth rates roughly following analytical expectations. The non-linear saturation of the instability appears to be caused by secondary shear instabilities and is also accompanied by a morphological change of the flow. We argue, however, that non-linear saturation likely occurs via other mechanisms in real stars where the separation of scales is larger than those reached by our simulations. We also observe dynamo action via the amplification of the axisymmetric poloidal magnetic field, suggesting that Tayler instability could be important for magnetic field generation and angular momentum transport in the radiative regions of evolving stars., Comment: 11 pages, 10 figures, published by MNRAS

We report discovery of a bright, nearby ($G = 13.8;\,\,d = 480\,\rm pc$) Sun-like star orbiting a dark object. We identified the system as a black hole candidate via its astrometric orbital solution from the Gaia mission. Radial velocities validated and refined the Gaia solution, and spectroscopy ruled out significant light contributions from another star. Joint modeling of radial velocities and astrometry constrains the companion mass to $M_2 = 9.62\pm 0.18\,M_{\odot}$. The spectroscopic orbit alone sets a minimum companion mass of $M_2>5\,M_{\odot}$; if the companion were a $5\,M_{\odot}$ star, it would be $500$ times more luminous than the entire system. These constraints are insensitive to the mass of the luminous star, which appears as a slowly-rotating G dwarf ($T_{\rm eff}=5850\,\rm K$, $\log g = 4.5$, $M=0.93\,M_{\odot}$), with near-solar metallicity ($\rm [Fe/H] = -0.2$) and an unremarkable abundance pattern. We find no plausible astrophysical scenario that can explain the orbit and does not involve a black hole. The orbital period, $P_{\rm orb}=185.6$ days, is longer than that of any known stellar-mass black hole binary. The system's modest eccentricity ($e=0.45$), high metallicity, and thin-disk Galactic orbit suggest that it was born in the Milky Way disk with at most a weak natal kick. How the system formed is uncertain. Common envelope evolution can only produce the system's wide orbit under extreme and likely unphysical assumptions. Formation models involving triples or dynamical assembly in an open cluster may be more promising. This is the nearest known black hole by a factor of 3, and its discovery suggests the existence of a sizable population of dormant black holes in binaries. Future Gaia releases will likely facilitate the discovery of dozens more., Comment: 30 pages, 15 figures, 1-2 black holes. Accepted to MNRAS

We constrain the orbital period ($P_{\rm orb}$) distribution of low-mass detached main-sequence eclipsing binaries (EBs) with light curves from the Zwicky Transient Facility (ZTF), which provides a well-understood selection function and sensitivity to faint stars. At short periods ($P_{\rm orb}\lesssim 2$ days), binaries are predicted to evolve significantly due to magnetic braking (MB), which shrinks orbits and ultimately brings detached binaries into contact. The period distribution is thus a sensitive probe of MB. We find that the intrinsic period distribution of low-mass ($0.1\lesssim M_1/M_{\odot} < 0.9$) binaries is basically flat (${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^0$), from $P_{\rm orb}=10$ days down to the contact limit. This is strongly inconsistent with predictions of classical MB models based on the Skumanich relation, which are widely used in binary evolution calculations and predict ${\rm d}N/{\rm d}P_{\rm orb} \propto P_{\rm orb}^{7/3}$ at short periods. The observed distributions are best reproduced by models in which the magnetic field saturates at short periods, with a MB torque that scales roughly as $\dot{J}\propto P_{\rm orb}^{-1}$, as opposed to $\dot{J} \propto P_{\rm orb}^{-3}$ in the standard Skumanich law. We also find no significant difference between the period distributions of binaries containing fully and partially convective stars. Our results confirm that a saturated MB law, which was previously found to describe the spin-down of rapidly rotating isolated M dwarfs, also operates in tidally locked binaries. We advocate using saturated MB models in binary evolution calculations. Our work supports previous suggestions that MB in cataclysmic variables (CVs) is much weaker than assumed in the standard evolutionary model, unless mass transfer leads to significant additional angular momentum loss in CVs., Comment: Table 1 included in source files. Accepted to MNRAS

Planetary engulfment events can occur while host stars are on the main sequence. The addition of rocky planetary material during engulfment will lead to refractory abundance enhancements in the host star photosphere, but the level of enrichment and its duration will depend on mixing processes that occur within the stellar interior, such as convection, diffusion, and thermohaline mixing. We examine engulfment signatures by modeling the evolution of photospheric lithium abundances. Because lithium can be burned before or after the engulfment event, it produces unique signatures that vary with time and host star type. Using MESA stellar models, we quantify the strength and duration of these signatures following the engulfment of a 1, 10, or 100 $M_{\oplus}$ planetary companion with bulk Earth composition, for solar-metallicity host stars with masses ranging from 0.5$-$1.4 $M_{\odot}$. We find that lithium is quickly depleted via burning in low-mass host stars ($\lesssim 0.7 \, M_\odot$) on a time scale of a few hundred Myrs, but significant lithium enrichment signatures can last for Gyrs in G-type stars ($\sim \! 0.9 \, M_{\odot}$). For more massive stars (1.3$-$1.4 $M_{\odot}$), engulfment can enhance internal mixing and diffusion processes, potentially decreasing the surface lithium abundance. Our predicted signatures from exoplanet engulfment are consistent with observed lithium-rich solar-type stars and abundance enhancements in chemically inhomogeneous binary stars., Comment: 13 pages, 9 figures, in review for MNRAS

The excitation of density and bending waves in Saturn's C ring by planetary oscillation modes presents a unique opportunity to learn about gas giant interiors and rotation. However, theoretical complications related to Saturn's rapid and differential rotation pose a barrier to the full utilization of ring wave detections. We calculate oscillation modes using a complete, non-perturbative treatment of differential rotation modelled after Saturn's zonal winds in self-consistently computed, polytropic equilibria. We find that previous, approximate treatments of the effects of differential rotation in Saturn overestimate shifts in the frequencies of fundamental modes (f-modes) thought to be responsible for the majority of the waves detected in the C ring, due to an omitted modification of the equilibrium shape and structure of the planet by differential rotation. The bias introduced by these frequency overestimates is small, but significant relative to the uncertainties afforded by Cassini data. We additionally consider the non-perturbative effects of Saturn-like differential rotation on the rotational mixing of f-modes and internal gravity modes (g-modes), which is relevant to detections of multiple density waves with very closely split pattern speeds. We find that higher order rotational effects can produce orders-of-magnitude enhancements in the surface gravitational perturbations of g-modes dominated by large spherical harmonic degrees $\ell$, regardless of frequency separation from the sectoral f-mode. Despite this enhancement, we find that the observed fine-splitting of density waves is unlikely to involve g-modes dominated by $\ell\gtrsim 10$. This restriction may aid in the inference of possible internal structures for Saturn., Comment: 15 pages, 14 figures. Published in MNRAS. Update includes minor corrections made while in press

Red Clump stars have been found to be enhanced in lithium relative to stars at the tip of the Red Giant Branch (TRGB), which is unexpected in current stellar models. At the TRGB, stars undergo the helium flash, during which helium burning briefly generates roughly $10^9 \, L_\odot$ of power and drives vigorous convection within the star's core. The helium-burning shell excites large fluxes of internal gravity waves. Here we investigate whether or not these waves can deposit enough heat to destabilize the hydrogen-burning shell, generate a convection zone there, and thereby drive the Cameron-Fowler process to enhance surface $^{7}\mathrm{Li}$. We study this with detailed stellar evolution models, and find that while the waves deposit $\sim 10^6 \, L_\odot$ near the hydrogen-burning shell, and while this generally does produce a convection zone, the resulting convection does not reach high enough to merge with the envelope, and so cannot explain enhancements to surface $^{7}\mathrm{Li}$., Comment: 18 pages, 7 figures

In a substantial number of core-collapse supernovae, early-time interaction indicates a dense circumstellar medium (CSM) that may be produced by outbursts from the progenitor star. Wave-driven mass loss is a possible mechanism to produce these signatures, with previous work suggesting this mechanism is most effective for low mass ($ \sim \! 11 \, M_\odot$) SN progenitors. Using one-dimensional hydrodynamic simulations with MESA, we study the effects of this wave heating in SN progenitors of masses $M_{\rm ZAMS} = 10-13\, M_{\odot}$. This range encompasses stars that experience semi-degenerate, central neon burning and more degenerate, off-center neon ignition. We find that central Ne ignition at $M_{\rm ZAMS} = 11\, M_{\odot}$ produces a burst of intense wave heating that transmits $\sim10^{47}$ erg of energy at $10$ years before core collapse, whereas other masses experience smaller levels of wave heating. Wave heating does not hydrodynamically drive mass loss in any of our models and is unlikely to produce a very massive CSM on its own. However, wave heating can cause large radial expansion (by more than an order of magnitude), photospheric cooling, and luminosity brightening by up to $\sim 10^6\, L_{\odot}$ in hydrogen-poor stripped star models. Some type Ib/c progenitors could drastically change their appearance in the final years of their lives, with brightness in visual bands increasing by nearly 3 mags. Moreover, interaction with a close binary companion could drive intense mass loss, with implications for type Ibn and other interaction-powered SNe.

Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution., Comment: 63 pages, 15 figures, 2 tables, published in Nature on May 4, 2022

In some semi-detached binary systems, the donor star may transfer mass to the companion at a very high rate. We propose that, at sufficiently high mass-transfer rates such that the accretion disk around the companion becomes geometrically thick (or advection-dominated) near the disk outer radius, a large fraction of the transferred mass will be lost through the outer Lagrangian (L2) point, as a result of the excessive energy generated by viscous heating that cannot be efficiently radiated away. A physical model is constructed where the L2 mass loss fraction is given by the requirement that the remaining material in the disk has Bernoulli number equal to the L2 potential energy. Our model predicts significant L2 mass loss at mass-transfer rates exceeding a few times 10^{-4} Msun/yr. An equatorial circum-binary outflow (CBO) is formed in these systems. Implications for the orbital evolution and the observational appearance are discussed. In particular, (1) rapid angular momentum loss from the system tends to shrink the orbit and hence may increase the formation rate of mergers and gravitational-wave sources; (2) photons from the hot disk wind are reprocessed by the CBO into longer wavelength emission in the infrared bands, consistent with Spitzer observations of some ultra-luminous X-ray sources., Comment: MNRAS in press, 14 pages, 7 figures, plus appendix. Codes available here: https://github.com/wenbinlu/L2massloss

Massive stars exhibit a variety of instabilities, many of which are poorly understood. We explore instabilities induced by centrifugal forces and angular momentum transport in massive rotating stars. First, we derive and numerically solve linearized oscillation equations for adiabatic radial modes in polytropic stellar models. In the presence of differential rotation, we show that centrifugal and Coriolis forces combined with viscous angular momentum transport can excite stellar pulsation modes, under both low- or high-viscosity conditions. In the low-viscosity limit, which is common in real stars, we demonstrate how to compute mode growth/damping rates via a work integral. Finally, we build realistic rotating $30\,M_\odot$ star models and show that overstable (growing) radial modes are predicted to exist for most of the star's life, in the absence of non-adiabatic effects. Peak growth rates are predicted to occur while the star is crossing the Hertzsprung-Russell gap, though non-adiabatic damping may dominate over viscous driving, depending on the effective viscosity produced by convective and/or magnetic torques. Viscous instability could be a new mechanism to drive massive star pulsations and is possibly related to instabilities of luminous blue variable stars., Comment: 14 pages, 12 figures; revised in response to comments from a referee

Tidally tilted pulsators (TTPs) are an intriguing new class of oscillating stars in binary systems; in such stars, the pulsation axis coincides with the line of apsides, or semi-major axis, of the binary. All three TTPs discovered so far have been $\delta$~Scuti stars. In this Letter, we report the first conclusive discovery of tidally tilted pulsations in a subdwarf B (sdB) star. HD 265435 is an sdB--white dwarf binary with a 1.65-hr period that has been identified and characterized as the nearest potential Type Ia supernova progenitor. Using TESS 20-s cadence data from Sectors 44 and 45, we show that the pulsation axis of the sdB star has been tidally tilted into the orbital plane and aligned with the tidal axis of the binary. We identify 31 independent pulsation frequencies, 27 of which have between 1 and 7 sidebands separated by the orbital frequency ($\nu_{\rm orb}$), or multiples thereof. Using the observed amplitude and phase variability due to tidal tilting, we assign $\ell$ and $m$ values to most of the observed oscillation modes and use these mode identifications to generate preliminary asteroseismic constraints. Our work significantly expands our understanding of TTPs, as we now know that (i) they can be found in stars other than $\delta$~Scuti pulsators, especially highly-evolved stars that have lost their H-rich envelopes, and (ii) tidally tilted pulsations can be used to probe the interiors of stars in very tight binaries., Comment: 17 pages, 5 figures + 2 appendix figures, 1 table; accepted for publication in ApJL

It has been recently suggested that white dwarfs generate magnetic fields in a process analogous to the Earth. The crystallization of the core creates a compositional inversion that drives convection, and combined with rotation, this can sustain a magnetic dynamo. We reanalyse the dynamo mechanism, arising from the slow crystallization of the core, and find convective turnover times $t_{\rm conv}$ of weeks to months - longer by orders of magnitude than previously thought. With white dwarf spin periods $P\ll t_{\rm conv}$, crystallization-driven dynamos are almost always in the fast rotating regime, where the magnetic field $B$ is at least in equipartition with the convective motion and is possibly further enhanced by a factor of $B\propto (t_{\rm conv}/P)^{1/2}$, depending on the assumed dynamo scaling law. We track the growth of the crystallized core using MESA and compute the magnetic field $B(T_{\rm eff})$ as a function of the white dwarf's effective temperature $T_{\rm eff}$. We compare this prediction with observations and show that crystallization-driven dynamos can explain some - but not all - of the $\sim$MG magnetic fields measured for single white dwarfs, as well as the stronger fields measured for white dwarfs in cataclysmic variables, which were spun up by mass accretion to short $P$. Our $B(T_{\rm eff})$ curves might also explain the clustering of white dwarfs with Balmer emission lines around $T_{\rm eff}\approx 7500\textrm{ K}$., Comment: Accepted version (MNRAS): improved calculation of the convection zone and the magnetic diffusion time, and a larger observational sample

The angular momentum (AM) content of massive stellar cores helps to determine the natal spin rates of neutron stars and black holes. Asteroseismic measurements of low-mass stars have proven that stellar cores rotate slower than predicted by most prior work, so revised models are necessary. In this work, we apply an updated AM transport model based on the Tayler instability to massive helium stars in close binaries, in which tidal spin-up can greatly increase the star's AM. Consistent with prior work, these stars can produce highly spinning black holes upon core-collapse if the orbital period is less than $P_{\rm orb} \lesssim \! 1 \, {\rm day}$. For neutron stars, we predict a strong correlation between the pre-explosion mass and the neutron star rotation rate, with millisecond periods ($P_{\rm NS} \lesssim 5 \, {\rm ms}$) only achievable for massive ($M \gtrsim 10 \, M_\odot$) helium stars in tight ($P_{\rm orb} \lesssim 1 \, {\rm day}$) binaries. Finally, we discuss our models in relation to type Ib/c supernovae, superluminous supernove, gamma-ray bursts, and LIGO measurements of black hole spins. Our models are roughly consistent with the rates and energetics of these phenomena, with the exception of broad-lined Ic supernovae, whose high rates and ejecta energies are difficult to explain., Comment: Revised after response to MNRAS referee. Comments welcome!

The tidally tilted pulsators are a new type of oscillating star in close binary systems that have their pulsation axis in the orbital plane because of the tidal distortion caused by their companion. We describe this group of stars on the basis of the first three representatives discovered and illustrate the basic methods used for their analysis. Their value for astrophysical study is rooted in the combination of the strengths of binary star and asteroseismic analyses; pulsational mode identifications can be achieved because the oscillations are visible over nearly 360 degrees of aspect throughout the orbital cycle. An illustrative case of a particularly interesting system is presented., Comment: 6 pages, 6 figures, Proceedings of the 40th meeting of the Polish Astronomical Society, submitted

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

Recent discoveries of double neutron star (DNS) mergers and ultra-stripped supernovae (SNe) raise the questions of their origin and connection. We present the first 1D~model of a DNS progenitor system which is calculated self-consistently until an ultra-stripped iron core collapse. We apply the \texttt{MESA} code starting from a post-common envelope binary consisting of a $1.35\;M_\odot$ NS and a $3.20\;M_\odot$ zero-age main-sequence helium star and continue the modelling via Case~BB Roche-lobe overflow until the infall velocity of the collapsing iron core exceeds $1000\;{\rm km\,s}^{-1}$. The exploding star has a total mass of $\sim 1.90\;M_\odot$, consisting of a $\sim 0.29\;M_\odot$ He-rich envelope embedding a CO core of $\sim 1.61\;M_\odot$ and an iron-rich core of $\sim 1.50\;M_\odot$. The resulting second-born NS has an estimated mass of $\sim 1.44\;M_\odot$ and we discuss the fate of the post-SN system, as well as the mild recycling of the first-born NS. Depending on the initial conditions, this family of systems is anticipated to reproduce the DNS mergers detected by the LIGO-network., Comment: 8 pages, 5 figures, ApJL in press

The discovery of rapidly rising and fading supernovae powered by circumstellar interaction has suggested the pre-supernova mass eruption phase as a critical phenomenon in massive star evolution. It is important to understand the mass and radial extent of the circumstellar medium (CSM) from theoretically predicted mass ejection mechanisms. In this work, we study the wave heating process in massive hydrogen-poor stars, running a suite of stellar models in order to predict the wave energy and pre-explosion time scale of surface energy deposition. We survey stellar models with main sequence progenitor masses from 20--70 $M_{\odot}$ and metallicity from 0.002 to 0.02. Most of these models predict that less than $\sim \! 10^{47} \, {\rm erg}$ is deposited in the envelope, with the majority of the energy deposited in the last week of stellar evolution. This translates to CSM masses less than $\sim \! 10^{-2} \, M_\odot$ that extend to less than $\sim \! 10^{14} \, {\rm cm}$, too small to greatly impact the light curves or spectra of the subsequent supernovae, except perhaps during the shock breakout phase. However, a few models predict somewhat higher wave energy fluxes, for which we perform hydrodynamical simulations of the mass ejection process. Radiative transfer simulations of the subsequent supernovae predict a bright but brief shock-cooling phase that could be detected in some type Ib/c supernovae if they are discovered within a couple days of explosion., Comment: 19 pages, 36 figures. Accepted for publication in Astrophysical Journal, submitted at 21 September 2021, accepted at 1 October 2021, published at 9 December 2021

Stellar mergers are important processes in stellar evolution, dynamics, and transient science. However, it is difficult to identify merger remnant stars because they cannot easily be distinguished from single stars based on their surface properties. We demonstrate that merger remnants can potentially be identified through asteroseismology of red giant stars using measurements of the gravity mode period spacing together with the asteroseismic mass. For mergers that occur after the formation of a degenerate core, remnant stars have over-massive envelopes relative to their cores, which is manifested asteroseismically by a g~mode period spacing smaller than expected for the star's mass. Remnants of mergers which occur when the primary is still on the main sequence or whose total mass is less than $\approx\! 2 \, M_\odot$ are much harder to distinguish from single stars. Using the red giant asteroseismic catalogs of Vrard et al. 2016 and Yu et al. 2018, we identify $24$ promising candidates for merger remnant stars. In some cases, merger remnants could also be detectable using only their temperature, luminosity, and asteroseismic mass, a technique that could be applied to a larger population of red giants without a reliable period spacing measurement., Comment: 13 pages, 6 figures, 1 table, submitted to MNRAS

Normal mode oscillations in Saturn excite density and bending waves in the C Ring, providing a valuable window into the planet's interior. Saturn's fundamental modes (f modes) excite the majority of the observed waves, while gravito-inertial modes (rotationally modified g modes) associated with stable stratification in the deep interior provide a compelling explanation for additional density waves with low azimuthal wavenumbers m. However, multiplets of density waves with nearly degenerate frequencies, including an m=3 triplet, still lack a definitive explanation. We investigate the effects of rapid and differential rotation on Saturn's oscillations, calculating normal modes for independently constrained interior models. We use a non-perturbative treatment of rotation that captures the full effects of the Coriolis and centrifugal forces, and consequently the mixing of sectoral f modes with g modes characterized by very different spherical harmonic degrees. Realistic profiles for differential rotation associated with Saturn's zonal winds can enhance these mode interactions, producing detectable oscillations with frequencies separated by less than 1%. Our calculations demonstrate that a three-mode interaction involving an f mode and two g modes can feasibly explain the finely split m=3 triplet, although the fine-tuning required to produce such an interaction generally worsens agreement with seismological constraints provided by m=2 density waves. Our calculations additionally demonstrate that sectoral f mode frequencies are measurably sensitive to differential rotation in Saturn's convective envelope. Finally, we find that including realistic equatorial antisymmetry in Saturn's differential rotation profile couples modes with even and odd equatorial parity, producing oscillations that could in principle excite both density and bending waves simultaneously., Comment: 23 pages, 16 figures, published in PSJ. Replacement fixes typos identified while in press

A large fraction of known exoplanets have short orbital periods where tidal excitation of gravity waves within the host star causes the planets' orbits to decay. We study the effects of tidal resonance locking, in which the planet locks into resonance with a tidally excited stellar gravity mode. Because a star's gravity mode frequencies typically increase as the star evolves, the planet's orbital frequency increases in lockstep, potentially causing much faster orbital decay than predicted by other tidal theories. Due to nonlinear mode damping, resonance locking in Sun-like stars likely only operates for low-mass planets ($M \lesssim 0.1 \, M_{\rm Jup}$), but in stars with convective cores it can likely operate for all planetary masses. The orbital decay timescale with resonance locking is typically comparable to the star's main-sequence lifetime, corresponding to a wide range in effective stellar quality factor ($10^3 \lesssim Q' \lesssim 10^9$), depending on the planet's mass and orbital period. We make predictions for several individual systems and examine the orbital evolution resulting from both resonance locking and nonlinear wave dissipation. Our models demonstrate how short-period massive planets can be quickly destroyed by nonlinear mode damping, while short-period low-mass planets can survive, even though they undergo substantial inward tidal migration via resonance locking., Comment: Accepted to ApJ

AM CVn systems are a rare type of accreting binary that consists of a white dwarf and a helium-rich, degenerate donor star. Using the Zwicky Transient Facility (ZTF), we searched for new AM CVn systems by focusing on blue, outbursting stars. We first selected outbursting stars using the ZTF alerts. We cross-matched the candidates with $Gaia$ and Pan-STARRS catalogs. The initial selection of candidates based on the $Gaia$ $BP$-$RP$ contains 1751 unknown objects. We used the Pan-STARRS $g$-$r$ and $r$-$i$ color in combination with the $Gaia$ color to identify 59 high-priority candidates. We obtained identification spectra of 35 sources, of which 18 are high priority candidates, and discovered 9 new AM CVn systems and one magnetic CV which shows only He-II lines. Using the outburst recurrence time, we estimate the orbital periods which are in the range of 29 to 50 minutes. We conclude that targeted followup of blue, outbursting sources is an efficient method to find new AM CVn systems, and we plan to followup all candidates we identified to systematically study the population of outbursting AM CVn systems., Comment: Submitted, comments and feedback welcome