Search for: All records

Stars and planets move supersonically in a gaseous medium during planetary engulfment, stellar interactions, and within protoplanetary discs. For a nearly uniform medium, the relevant parameters are the Mach number and the size of the body, R, relative to its accretion radius, RA. Over many decades, numerical and analytical work has characterized the flow, the drag on the body, and the possible suite of instabilities. Only a limited amount of work has treated the stellar boundary as it is in many of these astrophysical settings, a hard sphere at R. Thus, we present new 3D athena++  hydrodynamic calculations for a large range of parameters. For RA ≪ R, the results are as expected for pure hydrodynamics with minimal impact from gravity, which we verify by comparing to experimental wind tunnel data in air. When RA ≈ R, a hydrostatically supported separation bubble forms behind the gravitating body, exerting significant pressure on the sphere and driving a recompression shock, which intersects with the bow shock. For RA ≫ R, the bubble transitions into an isentropic, spherically symmetric halo, as seen in earlier works. These two distinct regimes of flow morphology may be treated separately in terms of their shock stand-off distance and drag coefficients. Most importantly for astrophysical applications, we propose a new formula for the dynamical friction, which depends on the ratio of the shock stand-off distance to RA. That exploration also reveals the minimum size of the simulation domain needed to accurately capture the deflection of incoming streamlines due to gravity.

 

Although stellar radii from asteroseismic scaling relations agree at the per cent level with independent estimates for main sequence and most first-ascent red giant branch (RGB) stars, the scaling relations over-predict radii at the tens of per cent level for the most luminous stars ($R \gtrsim 30 \, \mathrm{R}_{\odot }$). These evolved stars have significantly superadiabatic envelopes, and the extent of these regions increase with increasing radius. However, adiabaticity is assumed in the theoretical derivation of the scaling relations as well as in corrections to the large frequency separation. Here, we show that a part of the scaling relation radius inflation may arise from this assumption of adiabaticity. With a new reduction of Kepler asteroseismic data, we find that scaling relation radii and Gaia radii agree to within at least 2 per cent for stars with $R \lesssim 30\, \mathrm{R}_{\odot }$, when treated under the adiabatic assumption. The accuracy of scaling relation radii for stars with $50\, \mathrm{R}_{\odot }\lesssim R \lesssim 100\, \mathrm{R}_{\odot }$, however, is not better than $10~{{\ \rm per \, cent}}-15~{{\ \rm per \, cent}}$ using adiabatic large frequency separation corrections. We find that up to one third of this disagreement for stars with $R \approx 100\, \mathrm{R}_{\odot }$ could be caused by the adiabatic assumption, and that this adiabatic error increases with radius to reach 10 per cent at the tip of the RGB. We demonstrate that, unlike the solar case, the superadiabatic gradient remains large very deep in luminous stars. A large fraction of the acoustic cavity is also in the optically thin atmosphere. The observed discrepancies may therefore reflect the simplified treatment of convection and atmospheres.

 

The detonation of an overlying helium layer on a 0.8–1.1M_⊙carbon–oxygen (CO) white dwarf (WD) can detonate the CO WD and create a thermonuclear supernova (SN). Many authors have recently shown that when the mass of the He layer is low (≲0.03M_⊙), the ashes from its detonation minimally impact the spectra and light curve from the CO detonation, allowing the explosion to appear remarkably similar to Type Ia SNe. These new insights motivate our investigation of dynamical He shell burning and our search for a binary scenario that stably accumulates thermally unstable He shells in the 0.01–0.08M_⊙range, thick enough to detonate, but also often thin enough for minimal impact on the observables. We first show that our improved nonadiabatic evolution of convective He shell burning in this range of shell mass leads to conditions ripe for a He detonation. We also find that a stable mass transfer scenario with a high-entropy He WD donor of mass 0.15–0.25M_⊙yields the He shell masses needed to achieve the double detonations. This scenario also predicts that the surviving He donor leaves with a spatial velocity consistent with the unusual runaway object, D6-2. We find that hot He WD donors originate in common-envelope events when a 1.3–2.0M_⊙star fills its Roche lobe at the base of the red giant branch at orbital periods of 1–10 days with the CO WD.

 

About ten percent of Sun-like (1–2M_⊙) stars will engulf a 1–10M_Jplanet as they expand during the red giant branch (RGB) or asymptotic giant branch (AGB) phase of their evolution. Once engulfed, these planets experience a strong drag force in the star’s convective envelope and spiral inward, depositing energy and angular momentum. For these mass ratios, the inspiral takes ∼10–10²yr (∼10²–10³orbits); the planet undergoes tidal disruption at a radius of ∼1R_⊙. We use the Modules for Experiments in Stellar Astrophysics (MESA) software instrument to track the stellar response to the energy deposition while simultaneously evolving the planetary orbit. For RGB stars, as well as AGB stars withM_p≲ 5M_Jplanets, the star responds quasi-statically but still brightens measurably on a timescale of years. In addition, asteroseismic indicators, such as the frequency spacing or rotational splitting, differ before and after engulfment. For AGB stars, engulfment of anM_p≳ 5M_Jplanet drives supersonic expansion of the envelope, causing a bright, red, dusty eruption similar to a “luminous red nova.” Based on the peak luminosity, color, duration, and expected rate of these events, we suggest that engulfment events on the AGB could be a significant fraction of low-luminosity red novae in the Galaxy. We do not find conditions where the envelope is ejected prior to the planet’s tidal disruption, complicating the interpretation of short-period giant planets orbiting white dwarfs as survivors of common envelope evolution.

 

Abstract Observations indicate that turbulent motions are present on most massive star surfaces. Starting from the observed phenomena of spectral lines with widths that are much larger than their thermal broadening (e.g., micro- and macroturbulence), and considering the detection of stochastic low-frequency variability (SLFV) in the Transiting Exoplanet Survey Satellite photometry, these stars clearly have large-scale turbulent motions on their surfaces. The cause of this turbulence is debated, with near-surface convection zones, core internal gravity waves, and wind variability being proposed. Our 3D gray radiation hydrodynamic (RHD) models previously characterized the convective dynamics of the surfaces, driven by near-surface convection zones, and provided reasonable matches to the observed SLFV of the most luminous massive stars. We now explore the complex emitting surfaces of these 3D RHD models, which strongly violate the 1D assumption of a plane-parallel atmosphere. By post-processing the gray RHD models with the Monte Carlo radiation transport code Sedona , we synthesize stellar spectra and extract information from the broadening of individual photospheric lines. The use of Sedona enables the calculation of the viewing angle and temporal dependence of spectral absorption line profiles. By combining uncorrelated temporal snapshots together, we compare the turbulent broadening from the 3D RHD models to the thermal broadening of the extended emitting region, showing that our synthesized spectral lines closely resemble the observed macroturbulent broadening from similarly luminous stars. More generally, the new techniques that we have developed will allow for systematic studies of the origins of turbulent velocity broadening from any future 3D simulations. 

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

 

UsingAthena++, we perform 3D radiation-hydrodynamic calculations of the radiative breakout of the shock wave in the outer envelope of a red supergiant (RSG) that has suffered core collapse and will become a Type IIP supernova. The intrinsically 3D structure of the fully convective RSG envelope yields key differences in the brightness and duration of the shock breakout (SBO) from that predicted in a 1D stellar model. First, the lower-density “halo” of material outside of the traditional photosphere in 3D models leads to a shock breakout at lower densities than 1D models. This would prolong the duration of the shock breakout flash at any given location on the surface to ≈1–2 hr. However, we find that the even larger impact is the intrinsically 3D effect associated with large-scale fluctuations in density that cause the shock to break out at different radii at different times. This substantially prolongs the SBO duration to ≈3–6 hr and implies a diversity of radiative temperatures, as different patches across the stellar surface are at different stages of their radiative breakout and cooling at any given time. These predicted durations are in better agreement with existing observations of SBO. The longer durations lower the predicted luminosities by a factor of 3–10 (L_bol∼ 10⁴⁴erg s⁻¹), and we derive the new scalings of brightness and duration with explosion energies and stellar properties. These intrinsically 3D properties eliminate the possibility of using observed rise times to measure the stellar radius via light-travel time effects.

 

Abstract We explore the three-dimensional properties of convective, luminous ( L ≈ 10 4.5 –10 5 L ⊙ ), hydrogen-rich envelopes of red supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry using Athena ++. These computations comprise ≈30% of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer ≈3 M ⊙ of the simulated M ≈ 15 M ⊙ stars. This work reveals a radius, R corr , around which the nature of the convection changes. For r > R corr , though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well studied in less luminous stars, but in RSGs, the near- (or above-)Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outward rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This “halo” of material will impact predictions for both shock breakout and early lightcurves of Type IIP supernovae. Inside of R corr , we find a nearly flat entropy profile as expected in the efficient regime of mixing-length theory (MLT). Radiation pressure provides ≈1/3 of the support against gravity in this region. Our comparisons to MLT suggest a mixing length of α = 3–4, consistent with the sizes of convective plumes seen in the simulations. The temporal variability of these 3D models is mostly on the timescale of the convective plume lifetimes (≈300 days), with amplitudes consistent with those observed photometrically. 

Abstract Increasing main-sequence stellar luminosity with stellar mass leads to the eventual dominance of radiation pressure in stellar-envelope hydrostatic balance. As the luminosity approaches the Eddington limit, additional instabilities (beyond conventional convection) can occur. These instabilities readily manifest in the outer envelopes of OB stars, where the opacity increase associated with iron yields density and gas-pressure inversions in 1D models. Additionally, recent photometric surveys (e.g., TESS) have detected excess broadband low-frequency variability in power spectra of OB star lightcurves, called stochastic low-frequency variability (SLFV). This motivates our novel 3D Athena++ radiation hydrodynamical (RHD) simulations of two 35 M ⊙ star envelopes (the outer ≈15% of the stellar radial extent), one on the zero-age main sequence and the other in the middle of the main sequence. Both models exhibit turbulent motion far above and below the conventional iron-opacity peak convection zone (FeCZ), obliterating any “quiet” part of the near-surface region and leading to velocities at the photosphere of 10–100 km s −1 , directly agreeing with spectroscopic data. Surface turbulence also produces SLFV in model lightcurves with amplitudes and power-law slopes that are strikingly similar to those of observed stars. The characteristic frequencies associated with SLFV in our models are comparable to the thermal time in the FeCZ (≈3–7 day −1 ). These ab initio simulations are directly validated by observations and, though more models are needed, we remain optimistic that 3D RHD models of main-sequence O-star envelopes exhibit SLFV originating from the FeCZ. 

We report the discovery of ZTF J0127+5258, a compact mass-transferring binary with an orbital period of 13.7 minutes. The system contains a white dwarf accretor, which likely originated as a post–common envelope carbon–oxygen (CO) white dwarf, and a warm donor (T_eff,donor= 16,400 ± 1000 K). The donor probably formed during a common envelope phase between the CO white dwarf and an evolving giant that left behind a helium star or white dwarf in a close orbit with the CO white dwarf. We measure gravitational wave–driven orbital inspiral with ∼51σsignificance, which yields a joint constraint on the component masses and mass transfer rate. While the accretion disk in the system is dominated by ionized helium emission, the donor exhibits a mixture of hydrogen and helium absorption lines. Phase-resolved spectroscopy yields a donor radial velocity semiamplitude of 771 ± 27 km s⁻¹, and high-speed photometry reveals that the system is eclipsing. We detect a Chandra X-ray counterpart withL_X∼ 3 × 10³¹erg s⁻¹. Depending on the mass transfer rate, the system will likely either evolve into a stably mass-transferring helium cataclysmic variable, merge to become an R CrB star, or explode as a Type Ia supernova in the next million years. We predict that the Laser Space Interferometer Antenna (LISA) will detect the source with a signal-to-noise ratio of 24 ± 6 after 4 yr of observations. The system is the first LISA-loud mass-transferring binary with an intrinsically luminous donor, a class of sources that provide the opportunity to leverage the synergy between optical and infrared time domain surveys, X-ray facilities, and gravitational-wave observatories to probe general relativity, accretion physics, and binary evolution.