BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Lars Bildsten

Lars Bildsten contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.HEastro-ph.GA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they undergo core-collapse supernovae (CCSNe). The resulting dense circumstellar medium (CSM) may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomena were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2-4 km/s convective motions seed the intrinsic convective instability associated with the high entropy material deep in the envelope, enabling gas from deep within the envelope to escape, increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by 10x along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Discovery of a double detonation thermonuclear supernova progenitor

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical Simulations of Convective 3-Dimensional Red Supergiant Envelopes

We explore the three-dimensional properties of convective, luminous ($L\approx10^{4.5}-10^{5}L_\odot$), Hydrogen-rich envelopes of Red Supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry using $\texttt{Athena++}$. These computations comprise $\approx30\%$ of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer $\approx 3M_\odot$ of the simulated $M\approx15M_\odot$ stars. This work reveals a radius, $R_\mathrm{corr}$, around which the nature of the convection changes. For $r>R_\mathrm{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 outwards 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 II-P supernovae. Inside of $R_\mathrm{corr}$, we find a nearly flat entropy profile as expected in the efficient regime of mixing-length-theory (MLT). Radiation pressure provides $\approx1/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 ($\approx300$ days), with amplitudes consistent with those observed photometrically.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shock Breakout in 3-Dimensional Red Supergiant Envelopes

Using Athena++, we perform 3D Radiation-Hydrodynamic calculations of the radiative breakout of the shock wave in the outer envelope of a red supergiant (RSG) which 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 $\approx$1-2 hours. 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 $\approx$3-6 hours 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_\mathrm{bol}\sim10^{44}\mathrm{erg\ s^{-1}}$), 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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Still Brighter than Pre-Explosion, SN 2012Z Did Not Disappear: Comparing Hubble Space Telescope Observations a Decade Apart

Type Iax supernovae represent the largest class of peculiar white-dwarf supernovae. The type Iax SN~2012Z in NGC 1309 is the only white dwarf supernova with a detected progenitor system in pre-explosion observations. Deep \textit{Hubble Space Telescope} images taken before SN~2012Z show a luminous, blue source that we have interpreted as a helium-star companion (donor) to the exploding white dwarf. We present here late-time \textit{HST} observations taken $\sim$1400 days after the explosion to test this model. We find the SN light curve can empirically be fit by an exponential decay model in magnitude units. The fitted asymptotic brightness is within $10\%$ of our latest measurements and approximately twice the brightness of the pre-explosion source. The decline of the light curve is too slow to be powered by $^{56}$Co or $^{57}$Co decay: if radioactive decay is the dominate power source, it must be from longer half-life species like $^{55}$Fe. Interaction with circumstellar material may contribute to the light curve, as may shock heating of the companion star. Companion-star models underpredict the observed flux in the optical, producing most of their flux in the UV at these epochs. A radioactively-heated bound remnant, left after only a partial disruption of the white dwarf, is also capable of producing the observed excess late-time flux. Our analysis suggests that the total ejecta + remnant mass is consistent with the Chandrasekhar mass for a range of type Iax supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Stochastic Low Frequency Variability in 3-Dimensional Radiation Hydrodynamical Models of Massive Star Envelopes

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$_\odot$ star envelopes (the outer $\approx$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 ($\approx$3-7$\,$days$^{-1}$). These simulations, which have no free parameters, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Massive Star's Dying Breaths: Pulsating Red Supergiants and Their Resulting Type IIP Supernovae

Massive stars undergo fundamental-mode and first-overtone radial pulsations with periods of 100-1000 days as Red Supergiants (RSGs). At large amplitudes, these pulsations substantially modify the outer envelope's density structure encountered by the outgoing shock wave from the eventual core collapse of these $M>9M_\odot$ stars. Using Modules for Experiments in Stellar Astrophysics (MESA), we model the effects of fundamental-mode and first-overtone pulsations in the RSG envelopes, and the resulting Type IIP supernovae (SNe) using MESA+STELLA. We find that, in the case of fundamental mode pulsations, SN plateau observables such as the luminosity at day 50, $L_{50}$, time-integrated shock energy $ET$, and plateau duration $t_{\rm p}$ are consistent with radial scalings derived considering explosions of non-pulsating stars. Namely, most of the effect of the pulsation is consistent with the behavior expected for a star of a different size at the time of explosion. However, in the case of overtone pulsations, the Lagrangian displacement is not monotonic. Therefore, in such cases, excessively bright or faint SN emission at different times reflects the underdense or overdense structure of the emitting region near the SN photosphere.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A new class of Roche lobe-filling hot subdwarf binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Comparing Moment-Based and Monte Carlo Methods of Radiation Transport Modeling for Type II-Plateau Supernova Light Curves

Time-dependent electromagnetic signatures from core-collapse supernovae are the result of detailed transport of the shock-deposited and radioactively-powered radiation through the stellar ejecta. Due to the complexity of the underlying radiative processes, considerable approximations are made to simplify key aspects of the radiation transport problem. We present a systematic comparison of the moment-based radiation hydrodynamical code STELLA and the Monte Carlo radiation transport code Sedona in the 1D modeling of Type II-Plateau supernovae. Based on explosion models generated from the Modules for Experiments in Stellar Astrophysics (MESA) instrument, we find remarkable agreements in the modeled light curves and the ejecta structure thermal evolution, affirming the fidelity of both radiation transport modeling approaches. The radiative moments computed directly by the Monte Carlo scheme in Sedona also verify the accuracy of the moment-based scheme. We find that the coarse resolutions of the opacity tables and the numerical approximations in STELLA have insignificant impact on the resulting bolometric light curves, making it an efficient tool for the specific task of optical light curve modeling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The first ultracompact Roche lobe-filling hot subdwarf binary

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Value of Progenitor Radius Measurements for Explosion Modeling of Type II-Plateau Supernovae

Using Modules for Experiments in Stellar Astrophysics (MESA)+STELLA, we show that very different physical models can adequately reproduce a specific observed Type II-Plateau Supernova (SN). We consider SN2004A, SN2004et, SN2009ib, SN2017eaw, and SN2017gmr, Nickel-rich ($M_\mathrm{Ni}>0.03M_\odot$) events with bolometric lightcurves and a well-sampled decline from the plateau. These events also have constraints on the progenitor radius, via a progenitor image, or, in the case of SN2017gmr, a radius from fitting shock-cooling models. In general, many explosions spanning the parameter space of progenitors can yield excellent lightcurve and Fe line velocity agreement, demonstrating the success of scaling laws in motivating models which match plateau properties for a given radius and highlighting the degeneracy between plateau luminosity and velocity in models and observed events, which can span over 50% in ejecta mass, radius, and explosion energy. This can help explain disagreements in explosion properties reported for the same event using different model calculations. Our calculations yield explosion properties when combined with pre-explosion progenitor radius measurements or a robust understanding of the outermost $<0.1\,M_\odot$ of material that quantifies the progenitor radius from SN observations a few days after explosion.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Variability of massive stars in M31 from the Palomar Transient Factory

Using data from the (intermediate) Palomar Transient Factory (iPTF), we characterize the time variability of ~500 massive stars in M31. Our sample is those stars which are spectrally typed by Massey and collaborators, including Luminous Blue Variables, Wolf-Rayets, and warm and cool supergiants. We use the high-cadence, long-baseline (~5 years) data from the iPTF survey, coupled with data-processing tools that model complex features in the light curves. We find widespread photometric (R-band) variability in the upper Hertzsprung Russell diagram (or CMD) with an increasing prevalence of variability with later spectral types. Red stars (V-I>1.5) exhibit larger amplitude fluctuations than their bluer counterparts. We extract a characteristic variability timescale, tch, via wavelet transformations that are sensitive to both continuous and localized fluctuations. Cool supergiants are characterized by longer timescales (>100 days) than the hotter stars. The latter have typical timescales of tens of days but cover a wider range, from our resolution limit of a few days to longer than 100 days timescales. Using a 60-night block of data straddling two nights with a cadence of around 2 minutes, we extracted tch in the range 0.1--10 days with amplitudes of a few percent for 13 stars. Though there is broad agreement between the observed variability characteristics in the different parts of the upper CMD with theoretical predictions, detailed comparison requires models with a more comprehensive treatment of the various physical processes operating in these stars such as pulsation, subsurface convection, and the effect of binary companions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Remnants of Subdwarf Helium Donor Stars Ejected from Close Binaries with Thermonuclear Supernovae

Some binary systems composed of a white dwarf (WD) and a hot subdwarf (sdB) helium star will make contact within the helium burning lifetime of the sdB star. The accreted helium on the WD inevitably undergoes a thermonuclear instability, causing a detonation that is expected to transition into the WD core and lead to a thermonuclear supernova while the donor orbits nearby with high velocity. Motivated by the recent discovery of fast-moving objects that occupy unusual locations on the HR diagram, we explore the impact of the thermonuclear supernovae on the donors in this specific double detonation scenario. We use MESA to model the binary up to the moment of detonation, then 3D Athena++ to model the hydrodynamic interaction of the supernova ejecta with the donor star, calculating the amount of mass that is stripped and the entropy deposited in the deep stellar interior by the strong shock that traverses it. We show that these donor remnants are ejected with velocities primarily set by their orbital speeds: $700-900\ \rm km\ s^{-1}$. We model the long-term thermal evolution of remnants by introducing the shock entropy into MESA models. In response to this entropy change, donor remnants expand and brighten for timescales ranging from $10^6-10^8$ years, giving ample time for these runaway stars to be observed in their inflated state before they leave the galaxy. Even after surface layers are stripped, some donors retain enough mass to resume core helium burning and further delay fading for more than $10^8$ years.

0 citations0 reviewsNo code linkedOpen access