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.

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.

Context: Knowledge about hot, massive stars is usually inferred from quantitative spectroscopy. To analyse non-spherical phenomena, the existing 1D codes must be extended to higher dimensions, and corresponding tools need to be developed. Aims: We present a 3D radiative transfer code that is capable of calculating continuum and line scattering problems in the winds of hot stars. By considering spherically symmetric test models, we discuss potential error sources, and indicate advantages and disadvantages by comparing different solution methods. Further, we analyse the UV resonance line formation in the winds of rapidly rotating O stars. Methods: We consider both a (simplified) continuum model including scattering and thermal sources, and a UV resonance line transition approximated by a two-level-atom. We applied the short-characteristics (SC) method, using linear or monotonic B\'ezier interpolations, to solve the equation of radiative transfer on a non-uniform Cartesian grid. To calculate scattering dominated problems, our solution method is supplemented by an accelerated $\Lambda$-iteration scheme. Results: For the spherical test models, the mean relative error of the source function is on the $5-20\,\%$ level, depending on the applied interpolation technique and the complexity of the considered model. All calculated line profiles are in excellent agreement with corresponding 1D solutions. The predicted line profiles from fast rotating stars show a distinct behaviour as a function of rotational speed and inclination. This behaviour is tightly coupled to the wind structure and the description of gravity darkening and stellar surface distortion. Conclusions: Our SC methods are ready to be used for quantitative analyses of UV resonance line profiles. When calculating optically thick continua, both SC methods give reliable results, in contrast to the alternative finite-volume method.

(Abridged) The behaviour of mass loss across bi-stability jump is a key uncertainty in models of massive stars. While an increase in mass loss is theoretically predicted, this has so far not been observationally confirmed. However, radiation-driven winds of massive stars are known to exhibit clumpy structures triggered by the line-deshadowing instability (LDI). Wind clumping affects empirical mass-loss rates inferred from density square-dependent spectral diagnostics. If clumping properties differ significantly for O and B supergiants across the bi-stability jump, this may help alleviate discrepancies between theory and observations. We investigate with analytical and numerical tools how the onset of clumpy structures behaves in the winds of O supergiants (OSG) and B supergiants (BSG) across the bi-stability jump. We derive a scaling relation for the linear growth rate of the LDI for a single optically thick line and apply it in both regimes. We run 1D time-dependent line-driven instability simulations to study the non-linear evolution of the LDI in clumpy OSG and BSG winds. Linear perturbation analysis for a single line shows that the LDI linear growth rate scales strongly with stellar effective temperature and terminal wind speed. This implies significantly lower growth rates for (cooler, slower) BSG winds than for OSG winds. This is confirmed by the non-linear simulations, which show significant differences in OSG and BSG wind structure formation, with the latter characterized by significantly weaker clumping factors and lower velocity dispersions. This suggests that lower correction factors due to clumping should be employed when deriving empirical mass-loss rates for BSGs on the cool side of the bi-stability jump. Moreover, the non-linear simulations provide a theoretical background toward explaining the general lack of observed intrinsic X-ray emission in (single) B star winds.

Radiative feedback from luminous, massive stars during their formation is a key process in moderating accretion onto the stellar object. In the prior papers in this series, we showed that one form such feedback takes is UV line-driven disc ablation. Extending on this study, we now constrain the strength of this effect in the parameter range of star and disc properties appropriate to forming massive stars. Simulations show that ablation rate depends strongly on stellar parameters, but that this dependence can be parameterized as a nearly constant, fixed enhancement over the wind mass loss rate, allowing us to predict the rate of disc ablation for massive (proto)stars as a function of stellar mass and metallicity. By comparing this to predicted accretion rates, we conclude that ablation is a strong feedback effect for very massive (proto)stars which should be considered in future studies of massive star formation.

In studying the formation of massive stars, it is essential to consider the strong radiative feedback on the stars' natal environments from their high luminosities ($10^4 \sim 10^6 L_\odot$). Given that massive stars contract to main-sequence-like radii before accretion finishes, one form this feedback takes is UV line-acceleration, resulting in outflows much like those expected from main-sequence massive stars. As shown by the prior papers in this series, in addition to driving stellar winds, such line forces also ablate the surface layers off of circumstellar discs within a few stellar radii of the stellar photosphere. This removal of material from an accretion disc in turn results in a decreased accretion rate onto the forming star. Quantifying this, however, requires accounting for the continuum optical depth of the disc along the non-radial rays required for the three-dimensional line-acceleration prescription used in this paper series. We introduce the "thin disc approximation", allowing these continuum optical depths arising from an optically thick but geometrically thin disc to be dynamically treated in the context of radiation-hydrodynamics simulations. Using this approximation in full dynamical simulations, we show that such continuum optical depth effects only reduce the disc ablation by 30 percent or less relative to previous simulations that ignored continuum absorption.

Context: State of the art quantitative spectroscopy of OB-stars compares synthetic spectra (calculated by means of 1D, spherically symmetric computer codes) with observations. Certain stellar atmospheres, however, show strong deviations from spherical symmetry, and need to be treated in 3D. Aims: We present a newly developed 3D radiative transfer code, tailored to the solution of the radiation field in rapidly expanding stellar atmospheres. We apply our code to the continuum transfer in wind-ablation models, and to the UV resonance line formation in magnetic winds. Methods: We have used a 3D finite-volume method for the solution of the equation of radiative transfer, to study continuum- and line-scattering problems. Convergence has been accelerated by a non-local approximate Lambda-iteration scheme. Particular emphasis has been put on careful (spherically symmetric) test cases. Results: Typical errors of the source functions, when compared to 1D solutions, are of the order of 10-20 %, and increase for optically thick continua. In circumstellar discs, the radiation temperatures in the (optically thin) transition region from wind to disc are quite similar to corresponding values in the wind. For MHD simulations of dynamical magnetospheres, the line profiles, calculated with our 3D code, agree well with previous solutions using a 3D-SEI method. When compared with profiles resulting from the `analytic dynamical magnetosphere' (ADM) model, significant differences become apparent. Conclusions: Due to similar radiation temperatures in the wind and the transition region to the disc, the same line-strength distribution can be applied within radiation hydrodynamic calculations for circumstellar discs in `accreting high-mass stars'. To properly describe the UV line formation in dynamical magnetospheres, the ADM model needs to be further developed, at least in a large part of the outer wind.

Spectropolarimetric surveys reveal that 8-10\% of OBA stars harbor large-scale magnetic fields, but thus far no such fields have been detected in any classical Be stars. Motivated by this, we present here MHD simulations for how a pre-existing Keplerian disc -- like that inferred to form from decretion of material from rapidly rotating Be stars -- can be disrupted by a rotation-aligned stellar dipole field. For characteristic stellar and disc parameters of a near-critically rotating B2e star, we find that a polar surface field strength of just 10 G can significantly disrupt the disc, while a field of 100 G, near the observational upper limit inferred for most Be stars, completely destroys the disc over just a few days. Our parameter study shows that the efficacy of this magnetic disruption of a disc scales with the characteristic plasma beta {(defined as the ratio between thermal and magnetic pressure)} in the disc, but is surprisingly insensitive to other variations, e.g. in stellar rotation speed, or the mass loss rate of the star's radiatively driven wind. The disc disruption seen here for even a modest field strength suggests that the presumed formation of such Be discs by decretion of material from the star would likely be strongly inhibited by such fields; this provides an attractive explanation for why no large-scale fields are detected from such Be stars.

We extend our previous study of radiative ablation of circumstellar disks by line-scattering of the star's radiation, accounting now for the effect of multiple line resonances off the stellar limb. For an analytic, three-dimensional model of the velocity structure of an equatorial Keplerian disk bounded at higher latitudes by a radially accelerating stellar wind outflow, we use root-finding methods to identify multiple resonances from a near-disk circumstellar location along starward rays both on and off the stellar core. Compared to our previous study that accounted only for the effect of on-core resonances in reducing the radiative driving through the scattering of radiation away from a near-disk circumstellar location, including off-limb resonances leads to additional radiative driving from scattering toward this location. Instead of the up-to-50% reduction in line-acceleration previously inferred from multiple resonance effects, we now find a more modest 15-20% net reduction.

OB stars exhibit various types of spectral variability historically associated with wind structures, including the apparently ubiquitous discrete absorption components (DACs). These features have been proposed to be caused either by magnetic fields or non-radial pulsations. In this second paper of this series, we revisit the canonical phenomenological hydrodynamical modelling used to explain the formation of DACs by taking into account modern observations and more realistic theoretical predictions. Using constraints on putative bright spots located on the surface of the O giant $\xi$ Persei derived from high precision space-based broadband optical photometry obtained with the Microvariability and Oscillations of STars (MOST) space telescope, we generate two-dimensional hydrodynamical simulations of co-rotating interaction regions in its wind. We then compute synthetic ultraviolet (UV) resonance line profiles using Sobolev Exact Integration and compare them with historical timeseries obtained by the International Ultraviolet Explorer (IUE) to evaluate if the observed behaviour of $\xi$ Persei's DACs is reproduced. Testing three different models of spot size and strength, we find that the classical pattern of variability can be successfully reproduced for two of them: the model with the smallest spots yields absorption features that are incompatible with observations. Furthermore, we test the effect of the radial dependence of ionization levels on line driving, but cannot conclusively assess the importance of this factor. In conclusion, this study self-consistently links optical photometry and UV spectroscopy, paving the way to a better understanding of cyclical wind variability in massive stars in the context of the bright spot paradigm.

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.

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.

This work explores a Pulsationally Driven Orbital Mass Ejection (PDOME) model for the launching of Classical Be star disks. Under this model, a combination of rapid rotation and non-radial pulsation modes contribute to placing material into the circumstellar environment. Several varieties of non-radial pulsation modes, characterized by their propagation direction and the relative phase of their velocity and density perturbations, are considered. As well, the orbital stability of material launched by such a mechanism is investigated.

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.