Source author record

N. D. Kee

N. D. Kee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.SR astro-ph.GA astro-ph.HE

Catalog footprint

What is connected

8works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

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.

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' 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.

preprint2020arXiv

A Detailed View of the Circumstellar Environment and Disk of the Forming O-star AFGL 4176

We present a detailed analysis of the disk and circumstellar environment of the forming O-type star AFGL 4176 mm1, placing results from the Atacama Large Millimeter/submillimeter Array (ALMA) into context with multiwavelength data. With ALMA, we detect seventeen 1.2 mm continuum sources within 5$''$ (21,000 au) of AFGL 4176 mm1. We find that mm1 has a spectral index of 3.4$\pm$0.2 across the ALMA band, with $>$87% of its 1.2 mm continuum emission from dust. The source mm2, projected 4200 au from mm1, may be a companion or a blueshifted knot in a jet. We also explore the morphological differences between the molecular lines detected with ALMA, finding 203 lines from 25 molecules, which we categorize into several morphological types. Our results show that AFGL 4176 mm1 provides an example of a forming O-star with a large and chemically complex disk, which is mainly traced by nitrogen-bearing molecules. Lines that show strong emission on the blueshifted side of the disk are predominantly oxygen-bearing, which we suggest are tracing a disk accretion shock. The molecules C$^{34}$S, H$_2$CS and CH$_{3}$CCN trace a slow wide-angle wind or dense structures in the outflow cavity walls. With the Australia Telescope Compact Array (ATCA), we detect a compact continuum source ($<$2000 $\times$ 760 au) at 1.2 cm, associated with mm1, of which $>$96% is from ionized gas. The ATCA NH$_3$(1,1) and (2,2) emission traces a large-scale (r$\sim$0.5 pc) rotating toroid with the disk source mm1 in the blueshifted part of this structure offset to the northwest.

preprint2020arXiv

Spiral arms and instability within the AFGL 4176 mm1 disc

We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

preprint2016arXiv

Line-driven ablation of circumstellar disks: I. Optically thin decretion disks of classical Oe/Be stars

The extreme luminosities of hot, massive stars drive strong stellar winds through UV line-scattering. For OB stars with an orbiting circumstellar disk, we explore the effect of such line-scattering in ablating disk material, initially focusing on the marginally optically thin decretion disks of classical Oe and Be stars. For this we apply a multi-dimensional radiation-hydrodynamics code, assuming optically thin ray tracing for the stellar continuum and a multi-ray Sobolev treatment of the line transfer. This accounts for desaturation of line-absorption by Keplerian shear in the disk, and associated driving by non-radial photons. Results show dense, intermediate-speed surface ablation, consistent with the strong, blue-shifted absorption seen in UV wind lines of Be shell stars. The asymptotic ablation rate is typically an order-unity factor times the stellar wind mass loss rate, leading to disk destruction times of order months to years for Be disks, consistent with observations. The much stronger radiative forces of O stars reduce this time to order days, making sustaining a disk difficult, and so providing a natural explanation for the rarity of Galactic Oe stars. Additionally, the weakened line-driving at lower metallicity implies both a reduction in the winds that help spin-down stars from near-critical rotation, and in the ablation of decretion disks, thus providing a natural explanation for the higher fraction of Classical Be stars, and the presence of Oe stars, in the Magellanic Clouds. We conclude with a discussion of future extensions to study line-driven ablation of denser, optically thick, accretion disks around pre-main-sequence massive stars.

preprint2015arXiv

Radiative Ablation of Disks Around Massive Stars

Hot, massive stars (spectral types O and B) have extreme luminosities ($10^4 -10^6 L_\odot$) that drive strong stellar winds through UV line-scattering. Some massive stars also have disks, formed by either decretion from the star (as in the rapidly rotating "Classical Be stars"), or accretion during the star's formation. This dissertation examines the role of stellar radiation in driving (ablating) material away from these circumstellar disks. A key result is that the observed month to year decay of Classical Be disks can be explained by line-driven ablation without, as previously done, appealing to anomalously strong viscous diffusion. Moreover, the higher luminosity of O stars leads to ablation of optically thin disks on dynamical timescales of order a day, providing a natural explanation for the lack of observed Oe stars. In addition to the destruction of Be disks, this dissertation also introduces a model for their formation by coupling observationally inferred non-radial pulsation modes and rapid stellar rotation to launch material into orbiting Keplerian disks. In contrast to such Be decretion disks, star-forming accretion disks are much denser and so are generally optically thick to continuum processes. To circumvent the computational challenges associated with radiation hydrodynamics through optically thick media, we develop an approximate method for treating continuum absorption in the limit of geometrically thin disks. The comparison of ablation with and without continuum absorption shows that accounting for disk optical thickness leads to less than a 50$\%$ reduction in ablation rate, implying that ablation rate depends mainly on stellar properties like luminosity.

preprint2014arXiv

Suppression of X-rays from radiative shocks by their thin-shell instability

We examine X-rays from radiatively cooled shocks, focusing on how their thin-shell instability reduces X-ray emission. For 2D simulations of collision between equal expanding winds, we carry out a parameter study of such instability as a function of the ratio of radiative vs. adiabatic-expansion cooling lengths. In the adiabatic regime, the extended cooling layer suppresses instability, leading to planar shock compression with X-ray luminosity that follows closely the expected (Lx ~ M^2) quadratic scaling with mass-loss rate M . In the strongly radiative limit, the X-ray emission now follows an expected linear scaling with mass loss (Lx ~ M), but the instability deforms the shock compression into extended shear layers with oblique shocks along fingers of cooled, dense material. The spatial dispersion of shock thermalization limits strong X-ray emission to the tips and troughs of the fingers, and so reduces the X-ray emission (here by about a factor 1/50) below what is expected from analytic radiative-shock models without unstable structure. Between these two limits, X-ray emission can switch between a high-state associated with extended shock compression, and a low-state characterized by extensive shear. Further study is needed to clarify the origin of this "shear mixing reduction factor" in X-ray emission, and its dependence on parameters like the shock Mach number.