Researcher profile

J. O. Sundqvist

J. O. Sundqvist contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
4topics
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

9 published item(s)

preprint2022arXiv

First 3D Radiation-Hydrodynamic Simulations of Wolf-Rayet Winds

Classical Wolf Rayet (WR) stars are direct supernova progenitors undergoing vigorous mass-loss. Understanding the dense and fast outflows of such WR stars is thus crucial for understanding advanced stages of stellar evolution, the dynamical feedback of massive stars on their environments, and characterizing the distribution of black hole masses. In this paper, we develop first time-dependent, multi-dimensional, radiation-hydrodynamical models of the extended optically thick atmospheres and wind outflows of hydrogen-free classical WR stars. A flux limiting radiation hydrodynamics approach is used on a finite volume mesh to model WR outflows. The opacities are described using a combination of tabulated Rosseland mean opacities and the enhanced line opacities expected within a supersonic flow. For high-luminosity models, a radiation-driven, dense, supersonic wind is launched from deep sub-surface regions associated with peaks in the Rosseland mean opacity. For a model with lower luminosity, on the other hand, the Rosseland mean opacity is not sufficient to sustain a net-radial outflow in the sub-surface regions. Rather, what develops in this case is a "standard" line-driven wind launched from the optically thin regions above an extended and highly turbulent atmosphere. We thus find here a natural transition from optically thick outflows of classical WR stars to optically thin winds of hot, compact sub-dwarfs; in our simulations this transition occurs approximately at a luminosity that is about 40% of the Eddington luminosity. Because of the changing character of the wind-launching mechanism, this transition is also accompanied by a large drop (on the low-luminosity end) in average mass-loss rate.

preprint2022arXiv

New predictions for radiation-driven, steady-state mass-loss and wind-momentum from hot, massive stars III. Updated mass-loss rates for stellar evolution

Massive stars lose a large fraction of their mass to radiation-driven winds throughout their entire life. These outflows impact both the life and death of these stars and their surroundings. Theoretical mass-loss rates of hot, massive stars are derived to be used in applications such as stellar evolution. The behaviour of these rates in the OB-star regime is analysed, and their effects on massive-star evolution predictions is studied. Dynamically-consistent models are computed by solving the spherically symmetric, steady-state equation-of-motion for a large grid of hot, massive stars. The radiative acceleration is derived from non-local thermodynamic equilibrium radiative transfer in the co-moving frame. The resulting mass-loss rates are used to derive a simple scaling recipe with stellar parameters and to evaluate some first impacts upon evolution tracks. We provide a new prescription for steady-state, radiation-driven mass-loss from hot, massive stars depending on their fundamental parameters. The rates for O-stars are lower by about a factor ~3 than the rates typically used in previous stellar-evolution calculations, where differences generally decrease with increasing luminosity and temperature. For cooler B-giants/supergiants we find larger discrepancies. This arises because we do not find any systematic increase in mass-loss rates below the so-called bi-stability region; indeed, our results do not show any sign of a significant bi-stability jump within the parameter range covered by the grid. Due to the lower mass-loss rates we find that envelopes are not easily stripped by means of standard steady-state winds, making it difficult to create classical Wolf-Rayet stars via this channel. However, a remaining key uncertainty regarding these predictions concerns unsteady mass loss for very high-luminosity stars close to the Eddington limit as well as the impact of non-line-driven winds.

preprint2022arXiv

Theoretical wind clumping predictions from 2D LDI models of O-star winds at different metallicities

Hot, massive (OB) stars experience strong line-driven stellar winds and mass loss. As the majority of efficient driving lines are metallic, the amount of wind driving and mass loss is dependent on the stellar metallicity Z. In addition, line-driven winds are intrinsically inhomogeneous and clumpy. However, to date, neither theoretical nor empirical studies of line-driven winds have investigated how such wind clumping may also depend on Z. We theoretically investigated the degree of wind clumping due to the line-deshadowing instability (LDI) as a function of Z. We performed two-dimensional hydrodynamic simulations of the LDI with an assumed one-dimensional radiation line force for a grid of O-star wind models with fixed luminosity, but with different metal contents by varying the accumulative line strength Qbar describing the total ensemble of driving lines. We find that, for this fixed luminosity, the amount of wind clumping decreases with metallicity. The decrease is clearly seen in the statistical properties of our simulations, but is nonetheless rather weak; a simple power-law fit for the dependence of the clumping factor f_cl = <rho^2>/<rho>^2 on metallicity yields f_cl ~ Z^(0.15 +/- 0.01). This implies that empirically derived power-law dependencies of mass-loss rate Mdot on metallicity -- which were previously inferred from spectral diagnostics effectively depending on Mdot*sqrt(f_cl) but without having any constraints on f_cl(Z) -- should be only modestly altered by clumping. We expect that this prediction can be directly tested using new data from the Hubble Space Telescope Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) project.

preprint2021arXiv

Analytic, dust-independent mass-loss rates for red supergiant winds initiated by turbulent pressure

Context. Red supergiants are observed to undergo vigorous mass-loss. However, to date, no theoretical model has succeeded in explaining the origins of these objects&#39; winds. This strongly limits our understanding of red supergiant evolution and Type II-P and II-L supernova progenitor properties. Aims. We examine the role that vigorous atmospheric turbulence may play in initiating and determining the mass-loss rates of red supergiant stars. Methods. We analytically and numerically solve the equations of conservation of mass and momentum, which we later couple to an atmospheric temperature structure, to obtain theoretically motivated mass-loss rates. We then compare these to state-of-the-art empirical mass-loss rate scaling formulae as well as observationally inferred mass-loss rates of red supergiants. Results. We find that the pressure due to the characteristic turbulent velocities inferred for red supergiants is sufficient to explain the mass-loss rates of these objects in the absence of the normally employed opacity from circumstellar dust. Motivated by this initial success, we provide a first theoretical and fully analytic mass-loss rate prescription for red supergiants. We conclude by highlighting some intriguing possible implications of these rates for future studies of stellar evolution, especially in light of the lack of a direct dependence on metallicity.

preprint2021arXiv

Binary-object spectral-synthesis in 3D (BOSS-3D) -- Modelling H-alpha emission in the enigmatic multiple system LB-1

Context: To quantitatively decode the information stored within an observed spectrum, detailed modelling of the physical state and accurate radiative transfer solution schemes are required. In the analysis of stellar spectra, the numerical model often needs to account for binary companions and 3D structures in the stellar envelopes. The enigmatic binary (or multiple) system LB-1 constitutes a perfect example of such a complex multi-D problem. Aims: To improve our understanding of the LB-1 system, we directly modelled the phase-dependent H-alpha line profiles of this system. To this end, we developed a multi-purpose binary-object spectral-synthesis code in 3D (BOSS-3D). Methods: BOSS-3D calculates synthetic line profiles for a given state of the circumstellar material. The standard pz-geometry commonly used for single stars is extended by defining individual coordinate systems for each involved object and by accounting for the appropriate coordinate transformations. The code is then applied to the LB-1 system, considering two main hypotheses, a binary containing a stripped star and Be star, or a B star and a black hole with a disc. Results: Comparing these two scenarios, neither model can reproduce the detailed phase-dependent shape of the H-alpha line profiles. A satisfactory match with the observations, however, is obtained by invoking a disc around the primary object in addition to the Be-star disc or the black-hole accretion disc. Conclusions: The developed code can be used to model synthetic line profiles for a wide variety of binary systems, ranging from transit spectra of planetary atmospheres, to post-asymptotic giant branch binaries including circumstellar and circumbinary discs and massive-star binaries with stellar winds and disc systems. For the LB-1 system, our modelling provides strong evidence that each object in the system contains a disc-like structure.

preprint2021arXiv

Radiation-Hydrodynamics with MPI-AMRVAC: Flux-Limited Diffusion

Radiation controls the dynamics and energetics of many astrophysical environments. To capture the coupling between the radiation and matter, however, is often a physically complex and computationally expensive endeavour. We develop a numerical tool to perform radiation-hydrodynamics simulations in various configurations at an affordable cost. We build upon the finite volume code MPI-AMRVAC to solve the equations of hydrodynamics on multi-dimensional adaptive meshes and introduce a new module to handle the coupling with radiation. A non-equilibrium, flux-limiting diffusion approximation is used to close the radiation momentum and energy equations. The time-dependent radiation energy equation is then solved within a flexible framework, accounting fully for radiation forces and work terms and further allowing the user to adopt a variety of descriptions for the radiation-matter interaction terms (the &#39;opacities&#39;). We validate the radiation module on a set of standard testcases for which different terms of the radiative energy equation predominate. As a preliminary application to a scientific case, we calculate spherically symmetric models of the radiation-driven and optically thick supersonic outflows from massive Wolf-Rayet stars. This also demonstrates our code&#39;s flexibility, as the illustrated simulation combines opacities typically used in static stellar structure models with a parametrised form for the enhanced line-opacity expected in supersonic flows. This new module provides a convenient and versatile tool to perform multi-dimensional and high resolution radiative-hydrodynamics simulations in optically thick environments with the MPI-AMRVAC code. The code is ready to be used for a variety of astrophysical applications, where a first target for us will be multi-dimensional simulations of stellar outflows from Wolf-Rayet stars.

preprint2021arXiv

The inner circumstellar dust of the red supergiant Antares as seen with VLT/SPHERE/ZIMPOL

The processes by which red supergiants lose mass are not fully understood thus-far and their mass-loss rates lack theoretical constraints. The ambient surroundings of the nearby M0.5 Iab star Antares offers an ideal environment to obtain detailed empirical information on the outflow properties at its onset, and hence indirectly, on the mode(s) of mass loss. We present and analyse optical VLT/SPHERE/ZIMPOL polarimetric imaging with angular resolution down to 23 milli-arcsec, sufficient to spatially resolve both the stellar disk and its direct surroundings. We detect a conspicuous feature in polarised intensity that we identify as a clump containing dust, which we characterise through 3D radiative transfer modelling. The clump is positioned behind the plane of the sky, therefore has been released from the backside of the star, and its inner edge is only 0.3 stellar radii above the surface. The current dust mass in the clump is $1.3^{+0.2}_{-1.0} \times 10^{-8}$ M$_{\odot}$, though its proximity to the star implies that dust nucleation is probably still ongoing. The ejection of clumps of gas and dust makes a non-negligible contribution to the total mass lost from the star which could possibly be linked to localised surface activity such as convective motions or non-radial pulsations.

preprint2021arXiv

Upper Mass-Loss Limits and Clumping in the Intermediate and Outer Wind Regions of OB stars

We probe the radial clumping stratification of OB stars in the intermediate and outer wind regions (r>~2 R*) to derive upper limits for mass-loss rates, and compare to current mass-loss implementation. Together with archival multi-wavelength data, our new far-infrared continuum observations for a sample of 25 OB stars (including 13 B Supergiants) uniquely constrain the clumping properties of the intermediate wind region. We derive the minimum radial stratification of the clumping factor through the stellar wind, fclmin(r), and the corresponding maximum mass-loss rate, Mdotmax, normalising clumping factors to the outermost wind region (clfar=1). The clumping degree for r>~2 R* decreases or stays constant with increasing radius for almost the whole sample. There is a dependence on luminosity class and spectral type at the intermediate region relative to the outer ones: O Supergiants (OSGs) present a factor 2 larger clumping factors than B Supergiants (BSGs). The maximum clumping of roughly 1/3 of the OB Supergiants occurs close to the wind base (r<~2 R*) and then decreases monotonically. This contrasts with the more frequent case where the lowermost clumping increases towards a maximum, and needs to be addressed by theoretical models. Additionally, the estimated Mdotmax for BSGs is at least one order of magnitude lower than theoretical values, whereas for OSGs our results and predictions agree within errors. Assuming values of clfar=4-9 from hydrodynamical models would imply a reduction of mass-loss rates included in stellar evolution models by a factor 2-3 for OSGs and by factors 6-200 for BSGs below the first bi-stability jump. This implies large reductions of mass-loss rates applied in evolution-models for BSGs, independently of the actual clumping properties of these winds, and a thorough re-investigation of BSG mass-loss rates and their effects on stellar evolution.

preprint2019arXiv

New predictions for radiation-driven, steady-state mass-loss and wind-momentum from hot, massive stars. I. Method and first results

[Abridged] Context: Radiation-driven mass loss plays a key role in the life-cycles of massive stars. However, basic predictions of such mass loss still suffer from significant quantitative uncertainties. Aims: We develop new radiation-driven, steady-state wind models for massive stars with hot surfaces, suitable for quantitative predictions of global parameters like mass-loss and wind-momentum rates. Methods: The simulations presented here are based on a self-consistent, iterative grid-solution to the spherically symmetric, steady-state equation of motion, using full NLTE radiative transfer solutions in the co-moving frame to derive the radiative acceleration. We do not rely on any distribution functions or parametrization for computation of the line force responsible for the wind driving. Results: In this first paper, we present models representing two prototypical O-stars in the Galaxy, one with a higher stellar mass M/Msun=59 and luminosity log10(L/Lsun)= 5.87 and one with a lower M/Msun= 27 and log10(L/Lsun)= 5.1. For these simulations, basic predictions for global mass-loss rates and velocity laws are given, and the influence from additional parameters like wind clumping and microturbulent speeds discussed. A key result is that our mass-loss rates are significantly lower than those predicted by the mass-loss recipes normally included in models of massive-star evolution. Conclusions: Our results support previous suggestions that Galactic O-star mass-loss rates may be overestimated in present-day stellar evolution models, and that new rates thus might be needed. Indeed, future papers in this series will incorporate our new models into such simulations of stellar evolution, extending the very first simulations presented here toward larger grids covering a range of metallicities, B supergiants across the bistability jump, and possibly also Wolf-Rayet stars.