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.

Tidal effects in planetary systems are the main driver in the orbital migration of natural satellites. They result from physical processes occurring deep inside celestial bodies, whose effects are rarely observable from surface imaging. For giant planet systems, the tidal migration rate is determined by poorly understood dissipative processes in the planet, and standard theories suggest an orbital expansion rate inversely proportional to the power 11/2 in distance, implying little migration for outer moons such as Saturn's largest moon, Titan. Here, we use two independent measurements obtained with the Cassini spacecraft to measure Titan's orbital expansion rate. We find Titan migrates away from Saturn at 11.3 $\pm$ 2.0 cm/year, corresponding to a tidal quality factor of Saturn of Q $\simeq$ 100, and a migration timescale of roughly 10 Gyr. This rapid orbital expansion suggests Titan formed significantly closer to Saturn and has migrated outward to its current position. Our results for Titan and five other moons agree with the predictions of a resonance locking tidal theory, sustained by excitation of inertial waves inside the planet. The associated tidal expansion is only weakly sensitive to orbital distance, motivating a revision of the evolutionary history of Saturn's moon system. The resonance locking mechanism could operate in other systems such as stellar binaries and exoplanet systems, and it may allow for tidal dissipation to occur at larger orbital separations than previously believed.

We present observations of ZTF18abfcmjw (SN2019dge), a helium-rich supernova with a fast-evolving light curve indicating an extremely low ejecta mass ($\approx 0.3\,M_\odot$) and low kinetic energy ($\approx 1.2\times 10^{50}\,{\rm erg}$). Early-time (<4 d after explosion) photometry reveal evidence of shock cooling from an extended helium-rich envelope of $\sim0.1\,M_\odot$ located at $\sim 3\times 10^{12}\,{\rm cm}$ from the progenitor. Early-time He II line emission and subsequent spectra show signatures of interaction with helium-rich circumstellar material, which extends from $\gtrsim 5\times 10^{13}\,{\rm cm}$ to $\gtrsim 2\times 10^{16}\,{\rm cm}$. We interpret SN2019dge as a helium-rich supernova from an ultra-stripped progenitor, which originates from a close binary system consisting of a mass-losing helium star and a low-mass main sequence star or a compact object (i.e., a white dwarf, a neutron star, or a black hole). We infer that the local volumetric birth rate of 19dge-like ultra-stripped SNe is in the range of 1400--8200$\,{\rm Gpc^{-3}\, yr^{-1}}$ (i.e., 2--12% of core-collapse supernova rate). This can be compared to the observed coalescence rate of compact neutron star binaries that are not formed by dynamical capture.

The hyper-velocity star S5-HVS1, ejected 5 Myr ago from the Galactic Center at 1800 km/s, was most likely produced by tidal break-up of a tight binary by the supermassive black hole SgrA*. Taking a Monte Carlo approach, we show that the former companion of S5-HVS1 was likely a main-sequence star between 1.5 and 6 solar masses and was captured into a highly eccentric orbit with pericenter distance in the range 1-10 AU and semimajor axis about $10^3$ AU. We then explore the fate of the captured star. We find that the heat deposited by tidally excited stellar oscillation modes leads to runaway disruption if the pericenter distance is smaller than about 3 AU. Over the past 5 Myr, its angular momentum has been significantly modified by orbital relaxation, which may stochastically drive the pericenter inwards below 3 AU and cause tidal disruption. We find an overall survival probability in the range 5% to 50%, depending on the local relaxation time in the close environment of the captured star, and the initial pericenter at capture. The pericenter distance of the surviving star has migrated to 10-100 AU, making it potentially the most extreme member of the S-star cluster. From the ejection rate of S5-HVS1-like stars, we estimate that there may currently be a few stars in such highly eccentric orbits. They should be detectable (typically Ks < 18.5 mag) by the GRAVITY instrument and by future Extremely Large Telescopes and hence provide an extraordinary probe of the spin of SgrA*.

Asteroseismology is a powerful tool for probing the internal structures of stars by using their natural pulsation frequencies. It relies on identifying sequences of pulsation modes that can be compared with theoretical models, which has been done successfully for many classes of pulsators, including low-mass solar-type stars, red giants, high-mass stars and white dwarfs. However, a large group of pulsating stars of intermediate mass--the so-called delta Scuti stars--have rich pulsation spectra for which systematic mode identification has not hitherto been possible. This arises because only a seemingly random subset of possible modes are excited, and because rapid rotation tends to spoil the regular patterns. Here we report the detection of remarkably regular sequences of high-frequency pulsation modes in 60 intermediate-mass main-sequence stars, allowing definitive mode identification. Some of these stars have space motions that indicate they are members of known associations of young stars, and modelling of their pulsation spectra confirms that these stars are indeed young.

Compact white dwarf (WD) binaries are important sources for space-based gravitational-wave (GW) observatories, and an increasing number of them are being identified by surveys like ZTF. We study the effects of nonlinear dynamical tides in such binaries. We focus on the global three-mode parametric instability and show that it has a much lower threshold energy than the local wave-breaking condition studied previously. By integrating networks of coupled modes, we calculate the tidal dissipation rate as a function of orbital period. We construct phenomenological models that match these numerical results and use them to evaluate the spin and luminosity evolution of a WD binary. While in linear theory the WD's spin frequency can lock to the orbital frequency, we find that such a lock cannot be maintained when nonlinear effects are taken into account. Instead, as the orbit decays, the spin and orbit go in and out of synchronization. Each time they go out of synchronization, there is a brief but significant dip in the tidal heating rate. While most WDs in compact binaries should have luminosities that are similar to previous traveling-wave estimates, a few percent should be about ten times dimmer because they reside in heating rate dips. This offers a potential explanation for the low luminosity of the CO WD in J0651. Lastly, we consider the impact of tides on the GW signal and show that LISA and TianGO can constrain the WD's moment of inertia to better than 1% for deci-Hz systems.

Rotation and mass loss are crucially interlinked properties of massive stars, strongly affecting their evolution and ultimate fate. Massive stars rotating near their breakup limit shed mass centrifugally, creating Be stars with circumstellar disks and possibly driving outbursts. Using the MESA stellar evolution code, we examine the effects of efficient angular momentum transport on the main sequence and post-main sequence rotational evolution of massive stars. In rapid rotators, angular momentum transported from the contracting core to the expanding envelope can spin up the surface layers past the breakup rate, particularly for stars near (or beyond) the end of the main sequence and in low-metallicity environments. We also demonstrate that centrifugal instabilities could arise in rapidly rotating massive stars, potentially triggering the S Doradus outbursts observed in luminous blue variable stars. Prior mass accretion from a binary companion increases both the likelihood and the intensity of centrifugal mass loss. We discuss implications for massive stellar evolution, Be stars, and luminous blue variables.

CO~Cam (TIC 160268882) is the second ``single-sided pulsator'' to be discovered. These are stars where one hemisphere pulsates with a significantly higher amplitude than the other side of the star. CO~Cam is a binary star comprised of an Am $\delta$~Sct primary star with $T_{\rm eff} = 7070 \pm 150$\,K, and a spectroscopically undetected G main-sequence secondary star. The dominant pulsating side of the primary star is centred on the L$_1$ point. We have modelled the spectral energy distribution combined with radial velocities, and independently the {\em TESS} light curve combined with radial velocities. Both of these give excellent agreement and robust system parameters for both stars. The $\delta$~Sct star is an oblique pulsator with at least four low radial overtone (probably) f~modes with the pulsation axis coinciding with the tidal axis of the star, the line of apsides. Preliminary theoretical modelling indicates that the modes must produce much larger flux perturbations near the L$_1$ point, although this is difficult to understand because the pulsating star does not come near to filling its Roche lobe. More detailed models of distorted pulsating stars should be developed. These newly discovered single-sided pulsators offer new opportunities for astrophysical inference from stars that are oblique pulsators in close binary stars.

It has long been suspected that tidal forces in close binary stars could modify the orientation of the pulsation axis of the constituent stars. Such stars have been searched for, but until now never detected. Here we report the discovery of tidally trapped pulsations in the ellipsoidal variable HD 74423 in TESS space photometry data. The system contains a Delta Scuti pulsator in a 1.6-d orbit, whose pulsation mode amplitude is strongly modulated at the orbital frequency, which can be explained if the pulsations have a much larger amplitude in one hemisphere of the star. We interpret this as an obliquely pulsating distorted dipole oscillation with a pulsation axis aligned with the tidal axis. This is the first time that oblique pulsation along a tidal axis has been recognized. It is unclear whether the pulsations are trapped in the hemisphere directed towards the companion or in the side facing away from it, but future spectral measurements can provide the solution. In the meantime, the single-sided pulsator HD 74423 stands out as the prototype of a new class of obliquely pulsating stars in which the interactions of stellar pulsations and tidal distortion can be studied.

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.

Robust evidence of an ice giant planet shedding its atmosphere around the white dwarf WD J0914+1914 represents a milestone in exoplanetary science, allowing us to finally supplement our knowledge of white dwarf metal pollution, debris discs and minor planets with the presence of a major planet. Here, we discuss the possible dynamical origins of this planet, WD J0914+1914 b. The very young cooling age of the host white dwarf (13 Myr) combined with the currently estimated planet-star separation of about 0.07 au imposes particularly intriguing and restrictive coupled constraints on its current orbit and its tidal dissipation characteristics. The planet must have been scattered from a distance of at least a few au to its current location, requiring the current or former presence of at least one more major planet in the system in the absence of a hidden binary companion. We show that WD J0914+1914 b could not have subsequently shrunk its orbit through chaotic f-mode tidal excitation (characteristic of such highly eccentric orbits) unless the planet was or is highly inflated and possibly had partially thermally self-disrupted from mode-based energy release. We also demonstrate that if the planet is currently assumed to reside on a near-circular orbit at 0.07 au, then non-chaotic equilibrium tides impose unrealistic values for the planet's tidal quality factor. We conclude that WD~J0914+1914~b either (i) actually resides interior to 0.07 au, (ii) resembles a disrupted "Super-Puff" whose remains reside on a circular orbit, or (iii) resembles a larger or denser ice giant on a currently eccentric orbit. Distinguishing these three possibilities strongly motivates follow-up observations.

Many Type II supernovae (SNe) show hot early (~30 days) emission, and a diversity in their light curves extending from the Type IIP to the Type IIL, which can be explained by interaction with dense and confined circumstellar material (CSM). We perform hydrodynamical simulations of red supergiants to model the ejection of CSM caused by wave heating during late-stage nuclear burning. Even a small amount of deposited energy ($10^{46}-10^{47}$ erg), which is roughly that expected due to waves excited by convection in the core, is sufficient to change the shapes of SN light curves and bring them into better agreement with observations. As a test case, we consider the specific example of SN 2017eaw, which shows that a nuclear burning episode is able to explain the light curve if it occurs ~150-450 days prior to core-collapse. Due to the long timescale it takes for the low energy shock to traverse the star, this would manifest as a pre-SN outburst ~50-350 days prior to the full-fledged SN. Applying work like this to other SNe will provide a direct connection between the SN and pre-SN outburst properties, which can be tested by future wide field surveys. In addition, we show that our models can qualitatively explain the short lived `flash-ionization' lines seen in the early spectra of many Type II SNe.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

We report the discovery of a detached double white dwarf binary with an orbital period of $\approx20.6\,\rm minutes$, PTF J053332.05+020911.6. The visible object in this binary, PTF J0533+0209B, is a $\approx 0.17\,M_\odot$ mass white dwarf with a helium-dominated atmosphere containing traces of hydrogen (DBA). This object exhibits ellipsoidal variations due to tidal deformation, and is the visible component in a single-lined spectroscopic binary with a velocity semi-amplitude of $K_B=618.7\pm6.9 \, \rm km \, s^{-1}$. We have detected significant orbital decay due to the emission of gravitational radiation, and we expect that the Laser Interferometer Space Antenna (LISA) will detect this system with a signal to noise of $8.4^{+4.2}_{-3.0}$ after four years of operation. Because this system already has a well determined orbital period, radial velocity semi-amplitude, temperature, atmospheric composition, surface gravity, and orbital decay rate, a LISA signal will help fully constrain the properties of this system by providing a direct measurement of its inclination. Thus, this binary demonstrates the synergy between electromagnetic and gravitational radiation for constraining the physical properties of an astrophysical object.

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

Our understanding of the core-collapse supernova explosion mechanism is incomplete. While the favoured scenario is delayed revival of the stalled shock by neutrino heating, it is difficult to reliably compute explosion outcomes and energies, which depend sensitively on the complex radiation hydrodynamics of the post-shock region. The dynamics of the (non-)explosion depend sensitively on how energy is transported from inside and near the proto-neutron star (PNS) to material just behind the supernova shock. Although most of the PNS energy is lost in the form of neutrinos, hydrodynamic and hydromagnetic waves can also carry energy from the PNS to the shock. We show that gravity waves excited by core PNS convection can couple with outgoing acoustic waves that present an appreciable source of energy and pressure in the post-shock region. Using one-dimensional simulations, we estimate the gravity wave energy flux excited by PNS convection and the fraction of this energy transmitted upward to the post-shock region as acoustic waves. We find wave energy fluxes near $10^{51}\,\mathrm{erg}\,\mathrm{s}^{-1}$ are likely to persist for $\sim1\,\mathrm{s}$ post-bounce. The wave pressure on the shock may exceed $10\%$ of the thermal pressure, potentially contributing to shock revival and, subsequently, a successful and energetic explosion. We also discuss how future simulations can better capture the effects of waves, and more accurately quantify wave heating rates.

We present the characterization of KIC 4142768, an eclipsing binary with two evolved A-type stars in an eccentric orbit with a period of 14 days. We measure the fundamental parameters of the two components ($M_1=2.05M_{\odot}, R_1=2.96R_{\odot}$ and $M_2=2.05M_{\odot}, R_2=2.51R_{\odot}$) by combining {\it Kepler} photometry and spectra from {\it Keck} HIRES. The measured surface rotation rates are only one-fifth of the pseudo-synchronous rate of the eccentric orbit. Fourier spectrum of the light curve reveals hybrid pulsations of $\delta$ Scuti and $\gamma$ Doradus type, with pulsation frequencies at about $15-18$ day$^{-1}$ for p modes and about $0.2-1.2$ day$^{-1}$ for low-frequency g modes. Some of the g modes are exact orbital harmonics and are likely tidally excited. Their pulsation amplitudes and phases both agree with predictions from the linear tidal theory for $l=2, m=2$ prograde modes. We examine the period spacing patterns in the free oscillating g modes and identify them mostly as prograde sectoral dipole modes. The unstable frequency range and frequency spacing of p modes and the inferred asymptotic g-mode period spacings both agree with the stellar model for the primary star evolved to a late stage of the main sequence. The inferred rotation rate of the convective core boundary is very slow, similar to the small surface rotation rate inferred from the spectroscopy. The measured surface and near-core rotation rates provide constraints for testing the mechanism of angular momentum transfer and tidal synchronization in evolved eccentric binary star systems.

A gas giant planet which survives the giant branch stages of evolution at a distance of many au and then is subsequently perturbed sufficiently close to a white dwarf will experience orbital shrinkage and circularization due to star-planet tides. The circularization timescale, when combined with a known white dwarf cooling age, can place coupled constraints on the scattering epoch as well as the active tidal mechanisms. Here, we explore this coupling across the entire plausible parameter phase space by computing orbit shrinkage and potential self-disruption due to chaotic f-mode excitation and heating in planets on orbits with eccentricities near unity, followed by weakly dissipative equilibrium tides. We find that chaotic f-mode evolution activates only for orbital pericentres which are within twice the white dwarf Roche radius, and easily restructures or destroys ice giants but not gas giants. This type of internal thermal destruction provides an additional potential source of white dwarf metal pollution. Subsequent tidal evolution for the surviving planets is dominated by non-chaotic equilibrium and dynamical tides which may be well-constrained by observations of giant planets around white dwarfs at early cooling ages.

General relativity predicts that short orbital period binaries emit significant gravitational radiation, and the upcoming Laser Interferometer Space Antenna (LISA) is expected to detect tens of thousands of such systems; however, few have been identified, and only one is eclipsing--the double white dwarf binary SDSS J065133.338+284423.37, which has an orbital period of 12.75 minutes. Here, we report the discovery of an eclipsing double white dwarf binary system with an orbital period of only 6.91 minutes, ZTF J153932.16+502738.8. This system has an orbital period close to half that of SDSS J065133.338+284423.37 and an orbit so compact that the entire binary could fit within the diameter of the planet Saturn. The system exhibits a deep eclipse, and a double-lined spectroscopic nature. We observe rapid orbital decay, consistent with that expected from general relativity. ZTF J153932.16+502738.8 is a significant source of gravitational radiation close to the peak of LISA's sensitivity, and should be detected within the first week of LISA observations.

The age of gravitational wave (GW) astronomy has begun, and black hole (BH) mergers detected by LIGO are providing novel constraints on massive star evolution. A major uncertainty in stellar theory is the angular momentum (AM) transport within the star that determines its core rotation rate and the resulting BH's spin. Internal rotation rates of low-mass stars measured from asteroseismology prove that AM transport is efficient, suggesting that massive stellar cores may rotate slower than prior expectations. We investigate AM transport via the magnetic Tayler instability, which can largely explain the rotation rates of low-mass stars and white dwarfs. Implementing an updated AM transport prescription into models of high-mass stars, we compute the spins of their BH remnants. We predict that BHs born from single stars rotate very slowly, with $a \sim 10^{-2}$, regardless of initial rotation rate, possibly explaining the low $\chi_{\rm eff}$ of most BH binaries detected by LIGO thus far. A limited set of binary models suggests slow rotation for many binary scenarios as well, although homogeneous evolution and tidal spin-up of post-common envelope helium stars can create moderate or high BH spins. We make predictions for the values of $\chi_{\rm eff}$ in future LIGO events, and we discuss implications for engine-powered transients.

The internal rotational dynamics of massive stars are poorly understood. If angular momentum (AM) transport between the core and the envelope is inefficient, the large core AM upon core-collapse will produce rapidly rotating neutron stars (NSs). However, observations of low-mass stars suggest an efficient AM transport mechanism is at work, which could drastically reduce NS spin rates. Here we study the effects of the baroclinic instability and the magnetic Tayler instability in differentially rotating radiative zones. Although the baroclinic instability may occur, the Tayler instability is likely to be more effective for AM transport. We implement Tayler torques as prescribed by Fuller et al. 2019 into models of massive stars, finding they remove the vast majority of the core's AM as it contracts between the main sequence and helium-burning phases of evolution. If core AM is conserved during core-collapse, we predict natal NS rotation periods of $P_{\rm NS} \approx 50-200 \, {\rm ms}$, suggesting these torques help explain the relatively slow rotation rates of most young NSs, and the rarity of rapidly rotating engine-driven supernovae. Stochastic spin-up via waves just before core-collapse, asymmetric explosions, and various binary evolution scenarios may increase the initial rotation rates of many NSs.

$\varepsilon$ Lupi A is a binary system consisting of two main sequence early B-type stars Aa and Ab in a short period, moderately eccentric orbit. The close binary pair is the only doubly-magnetic massive binary currently known. Using photometric data from the BRITE-Constellation we identify a modest heartbeat variation. Combining the photometry with radial velocities of both components we determine a full orbital solution including empirical masses and radii. These results are compared with stellar evolution models as well as interferometry and the differences discussed. We also find additional photometric variability at several frequencies, finding it unlikely these frequencies can be caused by tidally excited oscillations. We do, however, determine that these signals are consistent with gravity mode pulsations typical for slowly pulsating B stars. Finally we discuss how the evolution of this system will be affected by magnetism, determining that tidal interactions will still be dominant.

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.

Asteroseismology is the only observational tool in astronomy that can probe the interiors of stars, and is a benchmark method for deriving fundamental properties of stars and exoplanets. Over the coming decade, space-based and ground-based observations will provide a several order of magnitude increase of solar-like oscillators, as well as a dramatic increase in the number and quality of classical pulsator observations, providing unprecedented possibilities to study stellar physics and galactic stellar populations. In this white paper, we describe key science questions and necessary facilities to continue the asteroseismology revolution into the 2020's.

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

Type Iax supernovae may arise from failed explosions of white dwarfs that leave behind a bound remnant (i.e., a "postgenitor" star) that could be identified in wide field surveys. To understand their observational signatures, we simulate these white dwarf (WD) postgenitors from shortly after explosion until they move back down the WD cooling track, and we consider several possible WD masses and explosion energies. To predict the peculiar surface abundances of the WD postgenitors, our models take into account gravitational settling and radiative levitation. We find that radiative levitation is significant at temperatures above a mass-dependent critical temperature, typically in the range Teff ~ 50-100 * 10^3 K, significantly increasing surface abundances of iron-group elements. Due to enhanced iron group opacity compared to normal WDs, the postgenitor peak luminosity and cooling timescale depend sensitively on mass, with more massive WDs becoming brighter but cooling much faster. We discuss our results in light of recently discovered hypervelocity white dwarfs with peculiar surface compositions, finding that our low-mass postgenitor models match many of their observational characteristics. Finally, we explore the effects of thermohaline diffusion, tentatively finding that it strongly suppresses abundance enhancements created by radiative levitation, but more realistic modeling is required to reach a firm conclusion.

Type I rapidly fading supernovae (RFSNe) appear to originate from hydrogen-free stars with large radii that produce predominantly shock-cooling light curves, in contrast with more typical nickel-rich SNe Ibc. However, it remains to be determined what types of stars would produce bright shock-cooling light curves without significant contribution from radioactive nickel. Bare helium stars in the mass range ~2-4 solar masses are known to hydrostatically develop radii as large as 100 solar radii or more due to strong He and C shell burning outside of a core with a sharp density gradient. We produce several such stellar models and demonstrate that, when exploded, these helium giants can naturally produce RFSN light curves. Since many prototypical SNe Ibc should come from large-radius stars in this mass range as well, we predict that these RFSNe may be distinct from SNe Ibc solely due to the absence of substantial nickel.

Context. While rotation has a major impact on stellar structure and evolution, its effects are not well understood. Thanks to high- quality and long timebase photometric observations obtained with recent space missions, we are now able to study stellar rotation more precisely. Aims. We aim to constrain radial differential rotation profiles in gamma Doradus (gamma Dor) stars, and to develop new theoretical seismic diagnosis for such stars with rapid and potentially non-uniform rotation. Methods. We derive a new asymptotic description which accounts for the impact of weak differential near-core rotation on gravity- mode period spacings. The theoretical predictions are illustrated from pulsation computations with the code GYRE and compared with observations of gamma Dor stars. When possible, we also derive the surface rotation rates in these stars by detecting and analysing signatures of rotational modulation, and compute the core-to-surface rotation ratios. Results. Stellar rotation has to be strongly differential before its effects on period spacing patterns can be detected, unless multiple period spacing patterns can be compared. Six stars in our sample exhibit a single unexplained period spacing pattern of retrograde modes. We hypothesise that these are Yanai modes. Finally, we find signatures of rotational spot modulation in the photometric data of eight targets. Conclusions. If only one period spacing pattern is detected and analysed for a star, it is difficult to detect differential rotation. A rigidly rotating model will often provide the best solution. Differential rotation can only be detected when multiple period spacing patterns have been found for a single star or its surface rotation rate is known as well. This is the case for eight stars in our sample, revealing surface-to-core rotation ratios between 0.95 and 1.05.

Pre-supernova (SN) outbursts from massive stars may be driven by hydrodynamical wave energy emerging from the core of the progenitor star during late nuclear burning phases. Here, we examine the effects of wave heating in stars containing little or no hydrogen, i.e., progenitors of type IIb/Ib SNe. Because there is no massive hydrogen envelope, wave energy is thermalized near the stellar surface where the overlying atmospheric mass is small but the optical depth is large. Wave energy can thus unbind this material, driving an optically thick, super-Eddington wind. Using 1D hydodynamic MESA simulations of $\sim \! 5 \, M_\odot$ He stars, we find that wave heating can drive pre-SN outbursts composed of a dense wind whose mass loss rate can exceed $\sim \! 0.1 \, M_\odot/{\rm yr}$. The wind terminal velocities are a few $100 \, {\rm km}/{\rm s}$, and outburst luminosities can reach $\sim \! 10^6 \, L_\odot$. Wave-driven outbursts may be linked with observed or inferred pre-SN outbursts of type Ibn/transitional/transformational SNe, and pre-SN wave-driven mass loss is a good candidate to produce these types of SNe. However, we also show that non-linear wave breaking in the core of the star may prevent such outbursts in stars with thick convective helium-burning shells. Hence, only a limited subset of SN progenitors are likely to experience wave-driven pre-SN outbursts.

Tidal dissipation inside giant planets is important for the orbital evolution of their natural satellites. It is conventionally treated by parameterized equilibrium tidal theory, in which the tidal torque declines rapidly with distance, and orbital expansion was faster in the past. However, Lainey et al. (2017) find that some Saturnian satellites are currently migrating outward faster than predicted by equilibrium tidal theory. Resonance locking between satellites and internal oscillations of Saturn, proposed by Fuller et al. (2016), naturally matches the observed migration rates. Here, we show that the resonance locking theory predicts dynamical tidal perturbations to Saturn's gravitational field in addition to those produced by equilibrium tidal bulges. We show that these perturbations can likely be detected during Cassini's proximal orbits if migration of satellites results from resonant gravity modes, but will likely be undetectable if migration results from inertial wave attractors or dissipation of the equilibrium tide. Additionally, we show that the detection of gravity modes would place constraints on the size of the hypothetical stably stratified region in Saturn.

Tidal dissipation in binary star and planetary systems is poorly understood. Fortunately, eccentric binaries known as heartbeat stars often exhibit tidally excited oscillations, providing observable diagnostics of tidal circularization mechanisms and timescales. We apply tidal theories to observations of the heartbeat star KIC 8164262, which contains an F-type primary in a very eccentric orbit that exhibits a prominent tidally excited oscillation. We demonstrate that the prominent oscillation is unlikely to result from a chance resonance between tidal forcing and a stellar oscillation mode. However, the oscillation has a frequency and amplitude consistent with the prediction of resonance locking, a mechanism in which coupled stellar and orbital evolution maintain a stable resonance between tidal forcing and a stellar oscillation mode. The resonantly excited mode produces efficient tidal dissipation (corresponding to an effective tidal quality factor $Q \sim 5 \times 10^4$), such that tidal orbital decay/circularization proceeds on a stellar evolution time scale.

We present the analysis of KIC 8164262, a heartbeat star with a high-amplitude (~1 mmag), tidally resonant pulsation (a mode in resonance with the orbit) at 229 times the orbital frequency and a plethora of tidally induced g-mode pulsations (modes excited by the orbit). The analysis combines Kepler light curves with follow-up spectroscopic data from the Keck telescope, KPNO (Kitt Peak National Observatory) 4-m Mayal telescope and the 2.7-m telescope at the McDonald observatory. We apply the binary modelling software, PHOEBE, to the Kepler light curve and radial velocity data to determine a detailed binary star model that includes the prominent pulsation and Doppler boosting, alongside the usual attributes of a binary star model (including tidal distortion and reflection). The results show that the system contains a slightly evolved F star with an M secondary companion in a highly eccentric orbit (e = 0.886). We use the results of the binary star model in a companion paper (Fuller et al., 2017) where we show that the prominent pulsation can be explained by a tidally excited oscillation mode held near resonance by a resonance locking mechanism.

Heartbeat stars are eccentric binary stars in short period orbits whose light curves are shaped by tidal distortion, reflection, and Doppler beaming. Some heartbeat stars exhibit tidally excited oscillations and present new opportunities for understanding the physics of tidal dissipation within stars. We present detailed methods to compute the forced amplitudes, frequencies, and phases of tidally excited oscillations in eccentric binary systems. Our methods i) factor out the equilibrium tide for easier comparison with observations, ii) account for rotation using the traditional approximation, iii) incorporate non-adiabatic effects to reliably compute surface luminosity perturbations, iv) allow for spin-orbit misalignment, and v) correctly sum over contributions from many oscillation modes. We also discuss why tidally excited oscillations are more visible in hot stars with surface temperatures $T \! \gtrsim \! 6500 \, {\rm K}$, and we derive some basic probability theory that can be used to compare models with data in a statistical manner. Application of this theory to heartbeat systems can be used to determine whether observed tidally excited oscillations can be explained by chance resonances with stellar oscillation modes, or whether a resonance locking process is operating.

Early observations of supernovae (SNe) indicate that enhanced mass loss and pre-SN outbursts may occur in progenitors of many types of SNe. We investigate the role of energy transport via waves driven by vigorous convection during late-stage nuclear burning of otherwise typical $15 \, M_\odot$ red supergiant SNe progenitors. Using MESA stellar evolution models including 1D hydrodynamics, we find that waves carry $\sim \! 10^7 \, L_\odot$ of power from the core to the envelope during core neon/oxygen burning in the final years before core collapse. The waves damp via shocks and radiative diffusion at the base of the hydrogen envelope, which heats up fast enough to launch a pressure wave into the overlying envelope that steepens into a weak shock near the stellar surface, causing a mild stellar outburst and ejecting a small ($\lesssim 1 \, M_\odot$) amount of mass at low speed ($\lesssim \! 50 \, {\rm km}/{\rm s}$) roughly one year before the SN. The wave heating inflates the stellar envelope but does not completely unbind it, producing a non-hydrostatic pre-SN envelope density structure different from prior expectations. In our models, wave heating is unlikely to lead to luminous type IIn SNe, but it may contribute to flash-ionized SNe and some of the diversity seen in II-P/II-L SNe.

Eccentric binaries known as heartbeat stars experience strong dynamical tides as the stars pass through periastron, providing a laboratory to study tidal interactions. We measure the rotation periods of 24 heartbeat systems, using the Kepler light curves to identify rotation peaks in the Fourier transform. Where possible, we compare the rotation period to the pseudosynchronization period derived by Hut 1981. Few of our heartbeat stars are pseudosynchronized with the orbital period. For four systems, we were able to identify two sets of rotation peaks, which we interpret as the rotation from both stars in the binary. The majority of the systems have a rotation period that is approximately 3/2 times the pseudosynchronization period predicted by Hut 1981, suggesting that other physical mechanisms influence the stars' rotation, or that stars typically reach tidal spin equilibrium at a rotation period slightly longer than predicted.

{\iota} Ori is a well studied massive binary consisting of an O9 III + B1 III/IV star. Due to its high eccentricity (e = 0.764) and short orbital period (P orb = 29.13376 d), it has been considered to be a good candidate to show evidence of tidal effects; however, none have previously been identified. Using photometry from the BRITE-Constellation space photometry mission we have confirmed the existence of tidal distortions through the presence of a heartbeat signal at periastron. We combine spectroscopic and light curve analyses to measure the masses and radii of the components, revealing {\iota} Ori to be the most massive heartbeat system known to date. In addition, using a thorough frequency analysis, we also report the unprecedented discovery of multiple tidally induced oscillations in an O star. The amplitudes of the pulsations allow us to empirically estimate the tidal circularization rate, yielding an effective tidal quality factor Q $\approx 4 \times 10^{4}$ .

Extrasolar planets on eccentric short-period orbits provide a laboratory in which to study radiative and tidal interactions between a planet and its host star under extreme forcing conditions. Studying such systems probes how the planet's atmosphere redistributes the time-varying heat flux from its host and how the host star responds to transient tidal distortion. Here, we report the insights into the planet-star interactions in HAT-P-2's eccentric planetary system gained from the analysis of 350 hr of 4.5 micron observations with the Spitzer Space Telescope. The observations show no sign of orbit-to-orbit variability nor of orbital evolution of the eccentric planetary companion, HAT-P-2 b. The extensive coverage allows us to better differentiate instrumental systematics from the transient heating of HAT-P-2 b's 4.5 micron photosphere and yields the detection of stellar pulsations with an amplitude of approximately 40 ppm. These pulsation modes correspond to exact harmonics of the planet's orbital frequency, indicative of a tidal origin. Transient tidal effects can excite pulsation modes in the envelope of a star, but, to date, such pulsations had only been detected in highly eccentric stellar binaries. Current stellar models are unable to reproduce HAT-P-2's pulsations, suggesting that our understanding of the interactions at play in this system is incomplete.

Asteroseismology probes the interiors of stars by studying oscillation modes at a star's surface. Although pulsation spectra are well understood for solar-like oscillators, a substantial fraction of red giant stars observed by Kepler exhibit abnormally low-amplitude dipole oscillation modes. Fuller et al. (2015) suggests this effect is produced by strong core magnetic fields that scatter dipole internal gravity waves (IGWs) into higher multipole IGWs or magnetic waves. In this paper, we study the interaction of IGWs with a magnetic field to test this mechanism. We consider two background stellar structures: one with a uniform magnetic field, and another with a magnetic field that varies both horizontally and vertically. We derive analytic solutions to the wave propagation problem and validate them with numerical simulations. In both cases, we find perfect conversion from IGWs into magnetic waves when the IGWs propagate into a region exceeding a critical magnetic field strength. Downward propagating IGWs cannot reflect into upward propagating IGWs because their vertical wavenumber never approaches zero. Instead, they are converted into upward propagating slow (Alfvenic) waves, and we show they will likely dissipate as they propagate back into weakly magnetized regions. Therefore, strong internal magnetic fields can produce dipole mode suppression in red giants, and gravity modes will likely be totally absent from the pulsation spectra of sufficiently magnetized stars.

A white dwarf (WD) captured into a high-eccentricity orbit around a massive black hole (MBH) may undergo many pericenter passages before tidal disruption. During these passages, the tidal potential of the MBH excites internal oscillations or waves in the WD, and the dissipation of these oscillations can significantly influence the physical properties of the WD prior to its disruption. We calculate the amplitude of the tidally excited gravity (buoyancy) waves in the WD as a function of the pericenter distance and eccentricity for realistic WD models, under the assumption that these outgoing gravity waves are efficiently dissipated in outer layers of the WD by non-linear effects or radiative damping. We obtain fitting formulae for the tidal energy and angular momentum transfer rates as well as the tidal heating rate. We find that these dynamical tides are much weaker than gravitational radiation in driving the orbital decay of the WD-MBH binary, and they are also inefficient in changing the WD spin during the orbital evolution. Incorporating our computed tidal dissipation rate into a MESA-based WD evolution code, we find that tidal heating can lead to appreciable brightening of the WD and may induce runaway fusion in the hydrogen envelope well before the WD undergoes tidal disruption.

KIC 3230227 is a short period ($P\approx 7.0$ days) eclipsing binary with a very eccentric orbit ($e=0.6$). From combined analysis of radial velocities and {\it Kepler} light curves, this system is found to be composed of two A-type stars, with masses of $M_1=1.84\pm 0.18M_{\odot}$, $M_2=1.73\pm 0.17M_{\odot}$ and radii of $R_1=2.01\pm 0.09R_{\odot}$, $R_2=1.68\pm 0.08 R_{\odot}$ for the primary and secondary, respectively. In addition to an eclipse, the binary light curve shows a brightening and dimming near periastron, making this a somewhat rare eclipsing heartbeat star system. After removing the binary light curve model, more than ten pulsational frequencies are present in the Fourier spectrum of the residuals, and most of them are integer multiples of the orbital frequency. These pulsations are tidally driven, and both the amplitudes and phases are in agreement with predictions from linear tidal theory for $l=2, m=-2$ prograde modes.

The inner moons of Jupiter and Saturn migrate outwards due to tidal energy dissipation within the planets, the details of which remain poorly understood. We demonstrate that resonance locking between moons and internal oscillation modes of the planet can produce rapid tidal migration. Resonance locking arises due to the internal structural evolution of the planet and typically produces an outward migration rate comparable to the age of the solar system. Resonance locking predicts a similar migration timescale but a different effective tidal quality factor $Q$ governing the migration of each moon. It also predicts nearly constant migration timescales a function of semi-major axis, such that effective $Q$ values were larger in the past. Recent measurements of Jupiter and Saturn's moon systems find effective $Q$ values that are smaller than expected (and are different between moons), and which correspond to migration timescales of $\sim$10 Gyr. If confirmed, the measurements are broadly consistent with resonance locking as the dominant source of tidal dissipation in Jupiter and Saturn. Resonance locking also provides solutions to several problems posed by current measurements: it naturally explains the exceptionally small $Q$ governing Rhea's migration, it allows the large heating rate of Enceladus to be achieved in an equilibrium eccentricity configuration, and it resolves evolutionary problems arising from present-day migration/heating rates.

Heartbeat stars (HB stars) are a class of eccentric binary stars with close periastron passages. The characteristic photometric HB signal evident in their light curves is produced by a combination of tidal distortion, heating, and Doppler boosting near orbital periastron. Many HB stars continue to oscillate after periastron and along the entire orbit, indicative of the tidal excitation of oscillation modes within one or both stars. These systems are among the most eccentric binaries known, and they constitute astrophysical laboratories for the study of tidal effects. We have undertaken a radial velocity (RV) monitoring campaign of Kepler HB stars in order to measure their orbits. We present our first results here, including a sample of 21 Kepler HB systems, where for 19 of them we obtained the Keplerian orbit and for 3 other systems we did not detect a statistically significant RV variability. Results presented here are based on 218 spectra obtained with the Keck/HIRES spectrograph during the 2015 Kepler observing season, and they have allowed us to obtain the largest sample of HB stars with orbits measured using a single instrument, which roughly doubles the number of HB stars with an RV measured orbit. The 19 systems measured here have orbital periods from 7 to 90 d and eccentricities from 0.2 to 0.9. We show that HB stars draw the upper envelope of the eccentricity - period distribution. Therefore, HB stars likely represent a population of stars currently undergoing high eccentricity migration via tidal orbital circularization, and they will allow for new tests of high eccentricity migration theories.

Heartbeat stars are eccentric (e > 0.2) ellipsoidal variables whose light curves resemble a cardiogram. We present the observations and corresponding model of KIC 3749404, a highly eccentric (e = 0.66), short period (P = 20.3 d) heartbeat star with tidally induced pulsations. A binary star model was created using PHOEBE, which we modified to include tidally induced pulsations and Doppler boosting. The morphology of the photometric periastron variation (heartbeat) depends strongly on the eccentricity, inclination and argument of periastron. We show that the inclusion of tidally induced pulsations in the model significantly changes the parameter values, specifically the inclination and those parameters dependent on it. Furthermore, we determine the rate of apsidal advance by modelling the periastron variation at the beginning and end of the 4-yr Kepler data set and dividing by the elapsed time. We compare the model with the theoretical expectations for classical and general relativistic apsidal motion and find the observed rate to be two orders of magnitude greater than the theoretical rate. We find that the observed rate cannot be explained by tidally induced pulsations alone and consequently hypothesise the presence of a third body in the system.

The asteroseismology of red giant stars has continued to yield surprises since the onset of high-precision photometry from space-based observations. An exciting new theoretical result shows that the previously observed suppression of dipole oscillation modes in red giants can be used to detect strong magnetic fields in the stellar cores. A fundamental facet of the theory is that nearly all the mode energy leaking into the core is trapped by the magnetic greenhouse effect. This results in clear predictions for how the mode visibility changes as a star evolves up the red giant branch, and how that depends on stellar mass, spherical degree, and mode lifetime. Here, we investigate the validity of these predictions with a focus on the visibility of different spherical degrees. We find that mode suppression weakens for higher degree modes with an average reduction in the quadrupole mode visibility of up to 49% for the least evolved stars in our sample, and no detectable suppression of octupole modes, in agreement with the theoretical predictions. We furthermore find evidence for the influence of increasing mode lifetimes on the measured visibilities along the red giant branch, in agreement with previous independent observations. These results support the theory that strong internal magnetic fields are responsible for the observed suppression of non-radial modes in red giants. We also find preliminary evidence that stars with suppressed dipole modes on average have slightly lower metallicity than normal stars.

Recent asteroseismic analyses have revealed the presence of strong (B $\gtrsim 10^5$ G) magnetic fields in the cores of many red giant stars. Here, we examine the implications of these results for the evolution of stellar magnetic fields, and we make predictions for future observations. Those stars with suppressed dipole modes indicative of strong core fields should exhibit moderate but detectable quadrupole mode suppression. The long magnetic diffusion times within stellar cores ensure that dynamo-generated fields are confined to mass coordinates within the main sequence convective core, and the observed sharp increase in dipole mode suppression rates above $1.5 \, M_\odot$ may be explained by the larger convective core masses and faster rotation of these more massive stars. In clump stars, core fields of $\sim10^5 \, {\rm G}$ can suppress dipole modes, whose visibility should be equal to or less than the visibility of suppressed modes in ascending red giants. High dipole mode suppression rates in low-mass ($M \lesssim 2 \, M_\odot$) clump stars would indicate that magnetic fields generated during the main sequence can withstand subsequent convective phases and survive into the compact remnant phase. Finally, we discuss implications for observed magnetic fields in white dwarfs and neutron stars, as well as the effects of magnetic fields in various types of pulsating stars.

Magnetic fields play a role in almost all stages of stellar evolution. Most low-mass stars, including the Sun, show surface fields that are generated by dynamo processes in their convective envelopes. Intermediate-mass stars do not have deep convective envelopes, although 10% exhibit strong surface fields that are presumed to be residuals from the stellar formation process. These stars do have convective cores that might produce internal magnetic fields, and these might even survive into later stages of stellar evolution, but information has been limited by our inability to measure the fields below the stellar surface. Here we use asteroseismology to study the occurrence of strong magnetic fields in the cores of low- and intermediate-mass stars. We have measured the strength of dipolar oscillation modes, which can be suppressed by a strong magnetic field in the core, in over 3600 red giant stars observed by Kepler. About 20% of our sample show mode suppression but this fraction is a strong function of mass. Strong core fields only occur in red giants above 1.1 solar masses (1.1Msun), and the occurrence rate is at least 60% for intermediate-mass stars (1.6--2.0Msun), indicating that powerful dynamos were very common in the convective cores of these stars.

Internal stellar magnetic fields are inaccessible to direct observations and little is known about their amplitude, geometry and evolution. We demonstrate that strong magnetic fields in the cores of red giant stars can be identified with asteroseismology. The fields can manifest themselves via depressed dipole stellar oscillation modes, which arises from a magnetic greenhouse effect that scatters and traps oscillation mode energy within the core of the star. The Kepler satellite has observed a few dozen red giants with depressed dipole modes which we interpret as stars with strongly magnetized cores. We find field strengths larger than $\sim\! 10^5 \,{\rm G}$ may produce the observed depression, and in one case we infer a minimum core field strength of $\approx \! \! 10^7 \,{\rm G}$.

The core rotation rates of massive stars have a substantial impact on the nature of core-collapse supernovae and their compact remnants. We demonstrate that internal gravity waves (IGW), excited via envelope convection during a red supergiant phase or during vigorous late time burning phases, can have a significant impact on the rotation rate of the pre-SN core. In typical ($10 \, M_\odot \lesssim M \lesssim 20 \, M_\odot$) supernova progenitors, IGW may substantially spin down the core, leading to iron core rotation periods $P_{\rm min,Fe} \gtrsim 30 \, {\rm s}$. Angular momentum (AM) conservation during the supernova would entail minimum NS rotation periods of $P_{\rm min,NS} \gtrsim 3 \, {\rm ms}$. In most cases, the combined effects of magnetic torques and IGW AM transport likely lead to substantially longer rotation periods. However, the stochastic influx of AM delivered by IGW during shell burning phases inevitably spin up a slowly rotating stellar core, leading to a maximum possible core rotation period. We estimate maximum iron core rotation periods of $P_{\rm max,Fe} \lesssim 5 \times 10^3 \, {\rm s}$ in typical core-collapse supernova progenitors, and a corresponding spin period of $P_{\rm max, NS} \lesssim 500 \, {\rm ms}$ for newborn neutron stars. This is comparable to the typical birth spin periods of most radio pulsars. Stochastic spin-up via IGW during shell O/Si burning may thus determine the initial rotation rate of most neutron stars. For a given progenitor, this theory predicts a Maxwellian distribution in pre-collapse core rotation frequency that is uncorrelated with the spin of the overlying envelope.

We study low frequency waves that propagate in a region of layered semi-convection. Layered semi-convection is predicted to be present in stellar and planetary interiors and can significantly modify the rate of thermal and compositional mixing. We derive a series of analytical dispersion relations for plane-parallel layered semi-convection in the Boussinesq approximation using a matrix transfer formalism. We find that like a continuously stratified medium, a semi-convective staircase -- in which small convective regions are separated by sharp density jumps -- supports internal gravity waves (g-modes). When the wavelength is much longer than the distance between semi-convective steps, these behave nearly like g-modes in a continuously stratified medium. However, the g-mode period spacing in a semi-convective region is systematically {\em smaller} than in a continuously stratified medium, and it decreases with decreasing mode frequency. When the g-mode wavelength becomes comparable to the distance between semi-convective steps, the g-mode frequencies deviate significantly from those of a continuously stratified medium (the frequencies are higher). G-modes with vertical wavelengths smaller than the distance between semi-convective steps are evanescent and do not propagate in the staircase. Thus, there is a lower cutoff frequency for a given horizontal wavenumber. We generalize our results to gravito-inertial waves relevant for rapidly rotating stars and planets. Finally, we assess the prospects for detecting layered semi-convection using astero/planetary seismology.

Fast radio bursts (FRBs) are an emerging class of short and bright radio transients whose sources remain enigmatic. Within the galactic center, the non-detection of pulsars within the inner $\sim \!10\,{\rm pc}$ has created a missing pulsar problem that has intensified with time. With all reserve, we advance the notion that the two problems could be linked by a common solution: the collapse of neutron stars (NS) due to capture and sedimentation of dark matter (DM) within their cores. Bramante \& Linden (2014), Phys.\ Rev.\ Lett.~19, 191301 showed that certain DM properties allow for rapid NS collapse within the high DM density environments near galactic centers while permitting NS survival elsewhere. Each DM-induced collapse could generate an FRB as the NS magnetosphere is suddenly expelled. This scenario could explain several features of FRBs: their short time scales, large energies, locally produced scattering tails, and high event rates. We predict that FRBs are localized to galactic centers, and that our own galactic center harbors a large population of NS-mass ($M\sim1.4 M_\odot$) black holes. The DM-induced collapse scenario is intrinsically unlikely because it can only occur in a small region of allowable DM parameter space. However, if observed to occur, it would place tight constraints on DM properties.

Gravitational waves (GW) generated during a core-collapse supernova open a window into the heart of the explosion. At core bounce, progenitors with rapid core rotation rates exhibit a characteristic GW signal which can be used to constrain the properties of the core of the progenitor star. We investigate the dynamics of rapidly rotating core collapse, focusing on hydrodynamic waves generated by the core bounce and the GW spectrum they produce. The centrifugal distortion of the rapidly rotating proto-neutron star (PNS) leads to the generation of axisymmetric quadrupolar oscillations within the PNS and surrounding envelope. Using linear perturbation theory, we estimate the frequencies, amplitudes, damping times, and GW spectra of the oscillations. Our analysis provides a qualitative explanation for several features of the GW spectrum and shows reasonable agreement with nonlinear hydrodynamic simulations, although a few discrepancies due to non-linear/rotational effects are evident. The dominant early postbounce GW signal is produced by the fundamental quadrupolar oscillation mode of the PNS, at a frequency $0.70 \, {\rm kHz} \lesssim f \lesssim 0.80\,{\rm kHz}$, whose energy is largely trapped within the PNS and leaks out on a $\sim\!10$ ms timescale. Quasi-radial oscillations are not trapped within the PNS and quickly propagate outwards until they steepen into shocks. Both the PNS structure and Coriolis/centrifugal forces have a strong impact on the GW spectrum, and a detection of the GW signal can therefore be used to constrain progenitor properties.

We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations.

Recent asteroseismic advances have allowed for direct measurements of the internal rotation rates of many sub-giant and red giant stars. Unlike the nearly rigidly rotating Sun, these evolved stars contain radiative cores that spin faster than their overlying convective envelopes, but slower than they would in the absence of internal angular momentum transport. We investigate the role of internal gravity waves in angular momentum transport in evolving low mass stars. In agreement with previous results, we find that convectively excited gravity waves can prevent the development of strong differential rotation in the radiative cores of Sun-like stars. As stars evolve into sub-giants, however, low frequency gravity waves become strongly attenuated and cannot propagate below the hydrogen burning shell, allowing the spin of the core to decouple from the convective envelope. This decoupling occurs at the base of the sub-giant branch when stars have surface temperatures of roughly 5500 K. However, gravity waves can still spin down the upper radiative region, implying that the observed differential rotation is likely confined to the deep core near the hydrogen burning shell. The torque on the upper radiative region may also prevent the core from accreting high-angular momentum material and slow the rate of core spin-up. The observed spin-down of cores on the red giant branch cannot be totally attributed to gravity waves, but the waves may enhance shear within the radiative region and thus increase the efficacy of viscous/magnetic torques.

Seismology allows for direct observational constraints on the interior structures of stars and planets. Recent observations of Saturn's ring system have revealed the presence of density waves within the rings excited by oscillation modes within Saturn, allowing for precise measurements of a limited set of the planet's mode frequencies. We construct interior structure models of Saturn, compute the corresponding mode frequencies, and compare them with the observed mode frequencies. The fundamental mode frequencies of our models match the observed frequencies (of the largest amplitude waves) to an accuracy of $\sim 1 \%$, confirming that these waves are indeed excited by Saturn's f-modes. The presence of the lower amplitude waves (finely split in frequency from the f-modes) can only be reproduced in models containing gravity modes that propagate in a stably stratified region of the planet. The stable stratification must exist deep within the planet near the large density gradients between the core and envelope. Our models cannot easily reproduce the observed fine splitting of the $m=-3$ modes, suggesting that additional effects (e.g., significant latitudinal differential rotation) may be important.

KIC 8560861 (HD 183648) is a marginally eccentric (e=0.05) eclipsing binary with an orbital period of P_orb=31.973d, exhibiting mmag amplitude pulsations on time scales of a few days. We present the results of the complex analysis of high and medium-resolution spectroscopic data and Kepler Q0 -- Q16 long cadence photometry. The iterative combination of spectral disentangling, atmospheric analysis, radial velocity and eclipse timing variation studies, separation of pulsational features of the light curve, and binary light curve analysis led to the accurate determination of the fundamental stellar parameters. We found that the binary is composed of two main sequence stars with an age of 0.9\+-0.2 Gyr, having masses, radii and temperatures of M_1=1.93+-0.12 M_sun, R_1=3.30+-0.07 R_sun, T_eff1=7650+-100 K for the primary, and M_2=1.06+-0.08 M_sun, R_2=1.11+-0.03 R_sun, T_eff2=6450+-100 K for the secondary. After subtracting the binary model, we found three independent frequencies, two of which are separated by twice the orbital frequency. We also found an enigmatic half orbital period sinusoidal variation that we attribute to an anomalous ellipsoidal effect. Both of these observations indicate that tidal effects are strongly influencing the luminosity variations of HD 183648. The analysis of the eclipse timing variations revealed both a parabolic trend, and apsidal motion with a period of (P_apse)_obs=10,400+-3,000 y, which is three times faster than what is theoretically expected. These findings might indicate the presence of a distant, unseen companion.

Tidal interactions play an important role in the evolution and ultimate fate of compact white dwarf (WD) binaries. Not only do tides affect the pre-merger state (such as temperature and rotation rate) of the WDs, but they may also determine which systems merge and which undergo stable mass transfer. In this paper, we attempt to quantify the effects of rotation on tidal angular momentum transport in binary stars, with specific calculations applied to WD stellar models. We incorporate the effect of rotation using the traditional approximation, in which the dynamically excited gravity waves within the WDs are transformed into gravito-inertial Hough waves. The Coriolis force has only a minor effect on prograde gravity waves, and previous results predicting the tidal spin-up and heating of inspiraling WDs are not significantly modified. However, rotation strongly alters retrograde gravity waves and inertial waves, with important consequences for the tidal spin-down of accreting WDs. We identify new dynamical tidal forcing terms that arise from a proper separation of the equilibrium and dynamical tide components; these new forcing terms are very important for systems near synchronous rotation. Additionally, we discuss the impact of Stokes drift currents on the wave angular momentum flux. Finally, we speculate on how tidal interactions will affect super-synchronously rotating WDs in accreting systems.

Recent observations have revealed evidence for the global oscillations of Jupiter and Saturn, which can potentially provide a new window into the interior structure of giant planets. Motivated by these observations, we study the non-radial oscillation modes of giant planets containing a solid core. Our calculations include the elastic response of the core and consider a wide range of possible values of the core shear modulus. While the elasticity of the core only slightly changes the frequencies of acoustic modes, which reside mostly in the fluid envelope, it adds two new classes of shear modes that are largely confined to the core. We also calculate the effects of the Coriolis force on the planetary oscillation modes. In addition to changing the mode frequencies, the Coriolis force can cause the shear modes to mix with the f-modes. Such mixing occurs when the frequencies of the shear mode and the f-mode are close to each other, and results in "mixed modes" with similar properties that are slightly split in frequency. We discuss our results in light of the recent work by Hedman & Nicholson (2013), which revealed the presence of density waves in Saturn's C-ring that appear to be excited by the gravitational perturbations associated with the f-mode oscillations within Saturn. We find that the fine splitting in wave frequencies observed in the rings can in principle be explained by the rotation-induced mixing between core shear modes and f-modes, possibly indicating the presence of a solid core within Saturn. However, in our current calculations, which assume rigid-body rotation and include only first-order rotational effects, significant fine-tuning in the planetary model parameters is needed in order to achieve these mode mixings and to explain the observed fine frequency splitting. We briefly discuss other effects that may modify the f-modes and facilitate mode mixing.

We present Kepler photometry and ground based spectroscopy of KIC 4544587, a short-period eccentric eclipsing binary system with self-excited pressure and gravity modes, tidally excited modes, tidally influenced p modes, and rapid apsidal motion of 182 y per cycle. The primary and secondary components of KIC 4544587 reside within the delta Scuti and gamma Dor instability region of the Hurtzsprung-Russell diagram, respectively. By applying the binary modelling software PHOEBE to prewhitened Kepler photometric data and radial velocity data obtained using the William Herschel Telescope and 4-m Mayall telescope at KPNO, the fundamental parameters of this important system have been determined, including the stellar masses, 1.98+/-0.07 MSun and 1.60+/-0.06 MSun, and radii, 1.76+/-0.03 RSun and 1.42+/-0.02 RSun, for the primary and secondary components, respectively. Frequency analysis of the residual data revealed 31 modes, 14 in the gravity mode region and 17 in the pressure mode region. Of the 14 gravity modes 8 are orbital harmonics: a signature of tidal resonance. While the measured amplitude of these modes may be partially attributed to residual signal from binary model subtraction, we demonstrate through consideration of the folded light curve that these frequencies do in fact correspond to tidally excited pulsations. Furthermore, we present an echelle diagram of the pressure mode frequency region (modulo the orbital frequency) and demonstrate that the tides are also influencing the p modes. A first look at asteroseismology hints that the secondary component is responsible for the p modes, which is contrary to our expectation that the hotter star should pulsate in higher radial overtone, higher frequency p modes.

We report that Kepler Object of Interest 256 (KOI-256) is a mutually eclipsing post-common envelope binary (ePCEB), consisting of a cool white dwarf (M = 0.592 +/- 0.089 MSun, R = 0.01345 +/- 0.00091 RSun, Teff = 7100 +/- 700 K) and an active M3 dwarf (M = 0.51 +/- 0.16 MSun, R = 0.540 +/- 0.014 RSun, Teff = 3450 +/- 50 K) with an orbital period of 1.37865 +/- 0.00001 days. KOI-256 is listed as hosting a transiting planet-candidate by Borucki et al. and Batalha et al.; here we report that the planet-candidate transit signal is in fact the occultation of a white dwarf as it passes behind the M dwarf. We combine publicly-available long- and short-cadence Kepler light curves with ground-based measurements to robustly determine the system parameters. The occultation events are readily apparent in the Kepler light curve, as is spin-orbit synchronization of the M dwarf, and we detect the transit of the white dwarf in front of the M dwarf halfway between the occultation events. The size of the white dwarf with respect to the Einstein ring during transit (REin = 0.00473 +/- 0.00055 RSun) causes the transit depth to be shallower than expected from pure geometry due to gravitational lensing. KOI-256 is an old, long-period ePCEB and serves as a benchmark object for studying the evolution of binary star systems as well as white dwarfs themselves, thanks largely to the availability of near-continuous, ultra-precise Kepler photometry.

We investigate the nature of tidal effects in compact triple-star systems. The hierarchical structure of a triple system produces tidal forcing at high frequencies unobtainable in binary systems, allowing for the tidal excitation of high frequency p-modes in the stellar components. The tidal forcing exists even for circular, aligned, and synchronized systems. We calculate the magnitude and frequencies of three-body tidal forcing on the central primary star for circular and coplanar orbits, and we estimate the amplitude of the tidally excited oscillation modes. We also calculate the secular orbital changes induced by the tidally excited modes, and show that they can cause significant orbital decay. During certain phases of stellar evolution, the tidal dissipation may be greatly enhanced by resonance locking. We then compare our theory to observations of HD 181068, which is a hierarchical triply eclipsing star system in the Kepler field of view. The observed oscillation frequencies in HD 181068 can be naturally explained by three-body tidal effects. We then compare the observed oscillation amplitudes and phases in HD 181068 to our predictions, finding mostly good agreement. Finally, we discuss the past and future evolution of compact triple systems like HD 181068.

Tidal dissipation in compact white dwarf (WD) binary systems significantly influences the physical conditions (such as surface temperature and rotation rate) of the WDs prior to mass transfer or merger. In these systems, the dominant tidal effects involve the excitation of gravity waves and their dissipation in the outer envelope of the star. We calculate the amplitude of tidally excited gravity waves in low-mass (0.3M_\odot) helium-core (He) WDs as a function of the tidal forcing frequency \omega. Like carbon-oxygen (CO) WDs studied in our previous paper, we find that the dimensionless tidal torque F(\omega) (inversely proportional to the effective tidal quality factor) has an erratic dependence on \omega. On average, F(\omega) scales approximately as \omega^6, and is several orders of magnitude smaller for He WDs than for CO WDs. We find that tidal torques can begin to synchronize the WD rotation when the orbital period is less than about a hour, although a nearly constant asynchronization is maintained even at small periods. We examine where the tidally excited gravity waves experience non-linear breaking or resonant absorption at a critical layer, allowing us to estimate the location and magnitude of tidal heating in the WD envelope. We then incorporate tidal heating in the MESA stellar evolution code, calculating the physical conditions of the WD as a function of orbital period for different WD models. We find that tidal heating makes a significant contribution to the WD luminosity for short-period (~10 min) systems such as SDSS J0651+2844. We also find that for WDs containing a hydrogen envelope, tidal heating can trigger runaway hydrogen shell burning, leading to a nova-like event before the onset of mass transfer.

Compact binary white dwarfs (WDs) undergoing orbital decay due to gravitational radiation can experience significant tidal heating prior to merger. In these WDs, the dominant tidal effect involves the excitation of outgoing gravity waves in the inner stellar envelope and the dissipation of these waves in the outer envelope. As the binary orbit decays, the WDs are synchronized from outside in (with the envelope synchronized first, followed by the core). We examine the deposition of tidal heat in the envelope of a Carbon-Oxygen WD and study how such tidal heating affects the structure and evolution of the WD. We show that significant tidal heating can occur in the star's degenerate hydrogen layer. This layer heats up faster than it cools, triggering runaway nuclear fusion. Such "tidal novae" may occur in all WD binaries containing a CO WD, at orbital periods between 5 min and 20 min, and precede the final merger by 10^5-10^6 years.

In compact white dwarf (WD) binary systems (with periods ranging from minutes to hours), dynamical tides involving the excitation and dissipation of gravity waves play a dominant role in determining the physical conditions of the WDs prior to mass transfer or binary merger. We calculate the amplitude of the tidally excited gravity waves as a function of the tidal forcing frequency \omega=2(\Omega-\Omega_s) (where \Omega is the orbital frequency and \Omega_s is the spin frequency) for several realistic carbon-oxygen WD models, assuming that the waves are efficiently dissipated in the outer layer of the star by nonlinear effects or radiative damping. The mechanism of wave excitation in WDs is complex due to the sharp features associated with composition changes inside the WD, and in our WD models gravity waves are launched just below the helium-carbon boundary. We find that the tidal torque on the WD and the related tidal energy transfer rate, \dot E_{\rm tide}, depend on \omega in an erratic way. On average, \dot E_{\rm tide} scales approximately as \Omega^5\omega^5 for a large range of tidal frequencies. We also study the effects of dynamical tides on the long-term evolution of WD binaries. Above a critical orbital frequency \Omega_c, corresponding to an orbital period of order one hour (depending on WD models), dynamical tides efficiently drive \Omega_s toward \Omega, although a small, almost constant degree of asynchronization (\Omega-\Omega_s\sim {\rm constant}) is maintained even at the smallest binary periods. While the orbital decay is always dominated by gravitational radiation, the tidal energy transfer can induce significant phase error in the low-frequency gravitational waveforms, detectable by the planned LISA project. Tidal dissipation may also lead to significant heating of the WD envelope and brightening of the system long before binary merger.

We study the tidal excitation of gravity modes (g-modes) in compact white dwarf binary systems with periods ranging from minutes to hours. As the orbit of the system decays via gravitational radiation, the orbital frequency increases and sweeps through a series of resonances with the g-modes of the white dwarf. At each resonance, the tidal force excites the g-mode to a relatively large amplitude, transferring the orbital energy to the stellar oscillation. We calculate the eigenfrequencies of g-modes and their coupling coefficients with the tidal field for realistic non-rotating white dwarf models. Using these mode properties, we numerically compute the excited mode amplitude in the linear approximation as the orbit passes though the resonance, including the backreaction of the mode on the orbit. We also derive analytical estimates for the mode amplitude and the duration of the resonance, which accurately reproduce our numerical results for most binary parameters. We find that the g-modes can be excited to a dimensionless (mass-weighted) amplitude up to 0.1, with the mode energy approaching $10^{-3}$ of the gravitational binding energy of the star. This suggests that thousands of years prior to the binary merger, the white dwarf may be heated up significantly by tidal interactions. However, more study is needed since the physical amplitudes of the excited oscillation modes become highly nonlinear in the outer layer of the star, which can reduce the mode amplitude attained by tidal excitation.

Recent observation of the tidally-excited stellar oscillations in the main-sequence binary KOI-54 by the KEPLER satellite provides a unique opportunity for studying dynamical tides in eccentric binary systems. We develop a general theory of tidal excitation of oscillation modes of rotating binary stars, and apply our theory to tidally excited gravity modes (g-modes) in KOI-54. The strongest observed oscillations, which occur at 90 and 91 times the orbital frequency, are likely due to prograde m=2 modes (relative to the stellar spin axis) locked in resonance with the orbit. The remaining flux oscillations with frequencies that are integer multiples of the orbital frequency are likely due to nearly resonant m=0 g-modes; such axisymmetric modes generate larger flux variations compared to the m=2 modes, assuming that the spin inclination angle of the star is comparable to the orbital inclination angle. We examine the process of resonance mode locking under the combined effects of dynamical tides on the stellar spin and orbit and the intrinsic stellar spindown. We show that KOI-54 can naturally evolve into a state in which at least one m=2 mode is locked in resonance with the orbital frequency. Our analysis provides an explanation for the fact that only oscillations with frequencies less than 90-100 times the orbital frequency are observed. We have also found evidence from the published KEPLER result that three-mode nonlinear coupling occurs in the KOI-54 system. We suggest that such nonlinear mode coupling may explain the observed oscillations that are not harmonics of the orbital frequency.