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.
more »
« less
Synthesizing Spectra from 3D Radiation Hydrodynamic Models of Massive Stars Using Monte Carlo Radiation Transport
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.
more »
« less
- Award ID(s):
- 1663688
- NSF-PAR ID:
- 10414668
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 945
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 58
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Context. Asymptotic giant branch (AGB) stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics, which amplifies the uncertainties on stellar parameters and distances. Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics (RHD) simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code O PTIM 3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius using the average 3D intensity profile and, finally, compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 μ m, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths (3.299 ± 0.005 mas and 3.053 ± 0.006 mas at 1.52 and 1.70 μ m, respectively) and recovered the radius ( R = 307 ± 41 and R = 284 ± 38 R ⊙ ) using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.more » « less
-
more » « less
Abstract We explore the three-dimensional properties of convective, luminous (
L ≈ 104.5–105L ⊙), 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. -
more » « less
Abstract Massive stars die in catastrophic explosions that seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion’s character and the remnant mass depend on models of the star’s evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometric observations reveal a ubiquitous red noise signal on massive main sequence stars; a hypothesized source of this noise is gravity waves driven by core convection. We present three-dimensional simulations of massive star convection extending from the star’s centre to near its surface, with realistic stellar luminosities. Using these simulations, we predict the photometric variability due to convectively driven gravity waves at the surfaces of massive stars, and find that gravity waves produce photometric variability of a lower amplitude and lower characteristic frequency than the observed red noise. We infer that the photometric signal of gravity waves excited by core convection is below the noise limit of current observations, and thus the red noise must be generated by an alternative process.
-
more » « less
ABSTRACT We analyse two binary systems containing giant stars, V723 Mon (‘the Unicorn’) and 2M04123153+6738486 (‘the Giraffe’). Both giants orbit more massive but less luminous companions, previously proposed to be mass-gap black holes. Spectral disentangling reveals luminous companions with star-like spectra in both systems. Joint modelling of the spectra, light curves, and spectral energy distributions robustly constrains the masses, temperatures, and radii of both components: the primaries are luminous, cool giants ($T_{\rm eff,\, giant} = 3800$ and $4000\, \rm K$, $R_{\rm giant}= 22.5$ and $25\, {\rm R}_{\odot }$) with exceptionally low masses ($M_{\rm giant} \approx 0.4\, {\rm M}_{\odot }$) that likely fill their Roche lobes. The secondaries are only slightly warmer subgiants ($T_{\rm eff,\, 2} = 5800$ and $5150\, \rm K$, $R_2= 8.3$ and $9\, {\rm R}_{\odot }$) and thus are consistent with observed UV limits that would rule out main-sequence stars with similar masses ($M_2 \approx 2.8$ and ${\approx}1.8\, {\rm M}_{\odot }$). In the Unicorn, rapid rotation blurs the spectral lines of the subgiant, making it challenging to detect even at wavelengths where it dominates the total light. Both giants have surface abundances indicative of CNO processing and subsequent envelope stripping. The properties of both systems can be reproduced by binary evolution models in which a $1{-}2\, {\rm M}_{\odot }$ primary is stripped by a companion as it ascends the giant branch. The fact that the companions are also evolved implies either that the initial mass ratio was very near unity, or that the companions are temporarily inflated due to rapid accretion. The Unicorn and Giraffe offer a window into into a rarely observed phase of binary evolution preceding the formation of wide-orbit helium white dwarfs, and eventually, compact binaries containing two helium white dwarfs.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/1538-4357/acb701
