Search for: All records

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.

 

We explore the three-dimensional properties of convective, luminous (L≈ 10^4.5–10⁵L_⊙), hydrogen-rich envelopes of red supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry usingAthena++. These computations comprise ≈30% of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer ≈3M_⊙of the simulatedM≈ 15M_⊙stars. This work reveals a radius,R_corr, around which the nature of the convection changes. Forr>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 ofR_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.

 

We present the discovery of a new double-detonation progenitor system consisting of a hot subdwarf B (sdB) binary with a white dwarf companion with aP_orb= 76.34179(2) minutes orbital period. Spectroscopic observations are consistent with an sdB star during helium core burning residing on the extreme horizontal branch. Chimera light curves are dominated by ellipsoidal deformation of the sdB star and a weak eclipse of the companion white dwarf. Combining spectroscopic and light curve fits, we find a low-mass sdB star,M_sdB= 0.383 ± 0.028M_⊙with a massive white dwarf companion,M_WD= 0.725 ± 0.026M_⊙. From the eclipses we find a blackbody temperature for the white dwarf of 26,800 K resulting in a cooling age of ≈25 Myr whereas ourMESAmodel predicts an sdB age of ≈170 Myr. We conclude that the sdB formed first through stable mass transfer followed by a common envelope which led to the formation of the white dwarf companion ≈25 Myr ago. Using theMESAstellar evolutionary code we find that the sdB star will start mass transfer in ≈6 Myr and in ≈60 Myr the white dwarf will reach a total mass of 0.92M_⊙with a thick helium layer of 0.17M_⊙. This will lead to a detonation that will likely destroy the white dwarf in a peculiar thermonuclear supernova. PTF1 J2238+7430 is only the second confirmed candidate for a double-detonation thermonuclear supernova. Using both systems we estimate that at least ≈1% of white dwarf thermonuclear supernovae originate from sdB+WD binaries with thick helium layers, consistent with the small number of observed peculiar thermonuclear explosions.

 

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. 

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. 

Abstract Binary systems of a hot subdwarf B (sdB) star + a white dwarf (WD) with orbital periods less than 2–3 hr can come into contact due to gravitational waves and transfer mass from the sdB star to the WD before the sdB star ceases nuclear burning and contracts to become a WD. Motivated by the growing class of observed systems in this category, we study the phases of mass transfer in these systems. We find that because the residual outer hydrogen envelope accounts for a large fraction of an sdB star’s radius, sdB stars can spend a significant amount of time (∼tens of megayears) transferring this small amount of material at low rates (∼10 −10 –10 −9 M ⊙ yr −1 ) before transitioning to a phase where the bulk of their He transfers at much faster rates ( ≳10 −8 M ⊙ yr −1 ). These systems therefore spend a surprising amount of time with Roche-filling sdB donors at orbital periods longer than the range associated with He star models without an envelope. We predict that the envelope transfer phase should be detectable by searching for ellipsoidal modulation of Roche-filling objects with P orb = 30–100 minutes and T eff = 20,000–30,000 K, and that many (≥10) such systems may be found in the Galactic plane after accounting for reddening. We also argue that many of these systems may go through a phase of He transfer that matches the signatures of AM CVn systems, and that some AM CVn systems associated with young stellar populations likely descend from this channel. 

Abstract Helium star–carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code MESA , we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the time of explosion. We convert these model properties to near-UV/optical magnitudes assuming a blackbody spectrum and support this approach using a matched stellar atmosphere model. These models will be valuable to compare with pre-explosion imaging for future supernovae, though we emphasize the observational difficulty of detecting extremely blue companions. The pre-explosion source detected in association with SN 2012Z has been interpreted as a helium star binary containing an initially ultra-massive WD in a multiday orbit. However, extending our binary models to initial CO WD masses of up to 1.2 M ⊙ , we find that these systems undergo off-center carbon ignitions and thus are not expected to produce thermonuclear supernovae. This tension suggests that, if SN 2012Z is associated with a helium star–WD binary, then the pre-explosion optical light from the system must be significantly modified by the binary environment and/or the WD does not have a carbon-rich interior composition. 

Abstract In this work we derive the minimum allowed orbital periods of H-rich bodies ranging in mass from Saturn’s mass to 1 M ⊙ , emphasizing gas giants and brown dwarfs (BDs) over the range 0.0003–0.074 M ⊙ . Analytic fitting formulae for as a function of the mass of the body and as a function of the mean density are presented. We assume that the density of the host star is sufficiently high so as not to limit the minimum period. In many instances this implies that the host star is a white dwarf. This work is aimed, in part, toward distinguishing BDs from planets that are found transiting the host white dwarf without recourse to near-infrared or radial velocity measurements. In particular, orbital periods of ≲100 minutes are very likely to be BDs. The overall minimum period over this entire mass range is ≃37 minutes. 

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