BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

E. R. Parkin

E. R. Parkin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEastro-ph.SRastro-ph.IMastro-ph.GAphysics.comp-ph

Catalog footprint

What is connected

18works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

18 published item(s)

preprint2016arXiv

The turbulent destruction of clouds - III. Three dimensional adiabatic shock-cloud simulations

We present 3D hydrodynamic simulations of the adiabatic interaction of a shock with a dense, spherical cloud. We compare how the nature of the interaction changes with the Mach number of the shock, $M$, and the density contrast of the cloud, $χ$. We examine the differences with 2D axisymmetric calculations, perform detailed resolution tests, and compare "inviscid" results to those obtained with the inclusion of a $k$-$ε$ subgrid turbulence model. Resolutions of 32-64 cells per cloud radius are the minimum necessary to capture the dominant dynamical processes in 3D simulations, while the 3D inviscid and $k$-$ε$ simulations typically show very good agreement. Clouds accelerate and mix up to 5 times faster when they are poorly resolved. The interaction proceeds very similarly in 2D and 3D - although non-azimuthal modes lead to different behaviour, there is very little effect on key global quantities such as the lifetime of the cloud and its acceleration. In particular, we do not find significant differences in the hollowing or "voiding" of the cloud between 2D and 3D simulations with $M=10$ and $χ=10$, which contradicts previous work in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the origin of γ-ray emission in η Carina

η Car is the only colliding-wind binary for which high-energy γ rays are detected. Although the physical conditions in the shock region change on timescales of hours to days, the variability seen at GeV energies is weak and on significantly longer timescales. The γ-ray spectrum exhibits two features that can be interpreted as emission from the shocks on either side of the contact discontinuity. Here we report on the first time-dependent modelling of the non-thermal emission in η Car. We find that emission from primary electrons is likely not responsible for the γ-ray emission, but accelerated protons interacting with the dense wind material can explain the observations. In our model, efficient acceleration is required at both shocks, with the primary side acting as a hadron calorimeter, whilst on the companion side acceleration is limited by the flow time out of the system, resulting in changing acceleration conditions. The system therefore represents a unique laboratory for the exploration of hadronic particle acceleration in non-relativistic shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Global simulations of magnetorotational turbulence III: influence of field configuration and mass injection

The stresses produced by magnetorotational turbulence can provide effective angular momentum transport in accretion disks. However, questions remain about the ability of simulated disks to reproduce observationally inferred stress-to-gas-pressure ratios. In this paper we present a set of high resolution global magnetohydrodynamic disk simulations which are initialised with different field configurations: purely toroidal, vertical field lines, and nested poloidal loops. A mass source term is included which allows the total disk mass to equilibrate in simulations with long run times, and also enables the impact of rapid mass injection to be explored. Notably different levels of angular momentum transport are observed during the early-time transient disk evolution. However, given sufficient time to relax, the different models evolve to a statistically similar quasi-steady state with a stress-to-gas-pressure ratio, $α\sim 0.032-0.036$. The indication from our results is that {\it steady, isolated} disks may be unable to maintain a large-scale magnetic field or produce values for the stress-to-gas-pressure ratio implied by some observations. Supplementary simulations exploring the influence of trapping magnetic field, injecting vertical field, and rapidly injecting additional mass into the disk show that large stresses ($α\sim 0.1-0.25$) can be induced by these mechanisms. The simulations highlight the common late-time evolution and characteristics of turbulent disks for which the magnetic field is allowed to evolve freely. If the boundaries of the disk, the rate of injection of magnetic field, or the rate of mass replenishment are modified to mimic astrophysical disks, markedly different disk evolution occurs.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The 2.35 year itch of Cygnus OB2#9 III. X-ray and radio emission analysis based on three dimensional hydrodynamical modelling

X-ray and radio data recently acquired as part of a project to study Cyg OB2#9 are used to constrain physical models of the binary system, providing in-depth knowledge about the wind-wind collision and the thermal, and non-thermal, emission arising from the shocks. We use a three-dimensional, adaptive mesh refinement simulation (including wind acceleration, radiative cooling, and the orbital motion of the stars) to model the gas dynamics of the wind-wind collision. The simulation output is used as the basis for radiative transfer calculations considering the thermal X-ray emission and the thermal/non-thermal radio emission. To obtain good agreement with the X-ray observations, our initial mass-loss rate estimates require a down-shift by a factor of roughly 7.7 to $6.5\times10^{-7}$ and $7.5\times10^{-7}$ solar mass per year for the primary and secondary star, respectively. Furthermore, the low gas densities and high shock velocities in Cyg OB2#9 are suggestive of unequal electron and ion temperatures, and the X-ray analysis indicates that an (immediately post-shock) electron-ion temperature ratio of $\simeq 0.1$ is also required. The radio emission is dominated by (non-thermal) synchrotron emission. A parameter space exploration provides evidence against models assuming equipartition between magnetic and relativistic energy densities. However, fits of comparable quality can be attained with models having stark contrasts in the ratio of magnetic-to-relativistic energy densities. The radio models also reveal a subtle effect whereby inverse Compton cooling leads to an increase in emissivity as a result of the synchrotron characteristic frequency being significantly reduced. Finally, using the results of the radio analysis, we estimate the surface magnetic field strengths to be $\approx 0.3-52\;$G. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global simulations of magnetorotational turbulence I: convergence and the quasi-steady state

Magnetorotational turbulence provides a viable mechanism for angular momentum transport in accretion disks. We present global, three dimensional (3D), MHD accretion disk simulations that investigate the dependence of the turbulent stresses on resolution. Convergence in the time-and-volume-averaged stress-to-gas-pressure ratio, at a value of $\sim0.04$, is found for a model with radial, vertical, and azimuthal resolution of 12-51, 27, and 12.5 cells per scale-height (the simulation mesh is such that cells per scale-height varies in the radial direction). A control volume analysis is performed on the main body of the disk (|z|<2H) to examine the production and removal of magnetic energy. Maxwell stresses in combination with the mean disk rotation are mainly responsible for magnetic energy production, whereas turbulent dissipation (facilitated by numerical resistivity) predominantly removes magnetic energy from the disk. Re-casting the magnetic energy equation in terms of the power injected by Maxwell stresses on the boundaries of, and by Lorentz forces within, the control volume highlights the importance of the boundary conditions (of the control volume). The different convergence properties of shearing-box and global accretion disk simulations can be readily understood on the basis of choice of boundary conditions and the magnetic field configuration. Periodic boundary conditions restrict the establishment of large-scale gradients in the magnetic field, limiting the power that can be delivered to the disk by Lorentz forces and by stresses at the surfaces. The factor of three lower resolution required for convergence in turbulent stresses for our global disk models compared to stratified shearing-boxes is explained by this finding. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global simulations of magnetorotational turbulence II: turbulent energetics

Magnetorotational turbulence draws its energy from gravity and ultimately releases it via dissipation. However, the quantitative details of this energy flow have not been assessed for global disk models. In this work we examine the energetics of a well-resolved, three-dimensional, global magnetohydrodynamic accretion disk simulation by evaluating statistically-averaged mean-field equations for magnetic, kinetic, and internal energy using simulation data. The results reveal that turbulent magnetic (kinetic) energy is primarily injected by the correlation between Maxwell (Reynolds) stresses and shear in the (almost Keplerian) mean flow, and removed by dissipation. This finding differs from previous work using local (shearing-box) models, which indicated that turbulent kinetic energy was primarily sourced from the magnetic energy reservoir. Lorentz forces provide the bridge between the magnetic and kinetic energy reservoirs, converting ~ 1/5 of the total turbulent magnetic power input into turbulent kinetic energy. The turbulent energies (both magnetic and kinetic) are mainly driven by terms associated with the turbulent fields, with only a minor influence from mean magnetic fields. The interaction between mean and turbulent fields is most evident in the induction equation, with the mean radial magnetic field being strongly influenced by the turbulent electromotive force (EMF). During the quasi-steady turbulent state roughly 2/3 of the Poynting flux travels into the corona, with the remainder transporting magnetic energy in the radial direction. In contrast to previous studies, the stress-related part of the Poynting flux is found to dominate, which may have important implications for "reflection" models of Seyfert galaxy coronae that typically invoke a picture of buoyant rising of magnetic flux tubes via advection.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Long time-scale variability of X-ray binaries with late type giant companions

In this paper we propose and examine a physical mechanism which can lead to the generation of noise in the mass accretion rate of low mass X-ray binaries on time-scales comparable to the orbital period of the system. We consider modulations of mass captured by the compact object from the companion star's stellar wind in binaries with late type giants, systems which usually have long orbital periods. We show that a hydrodynamical interaction of the wind matter within a binary system even without eccentricity results in variability of the mass accretion rate with characteristic time-scales close to the orbital period. The cause of the variability is an undeveloped turbulent motion (perturbed motion without significant vorticity) of wind matter near the compact object. Our conclusions are supported by 3D simulations with two different hydrodynamic codes based on Lagrangian and Eulerian approaches -- the SPH code GADGET and the Eulerian code PLUTO. In this work we assume that the wind mass loss rate of the secondary is at the level of $(0.5-1)\times10^{-7} M_\odot$/year, required to produce observable variations of the mass accretion rate on the primary. This value is higher than that, estimated for single giant stars of this type, but examples of even higher mass loss rate of late type giants in binaries do exist. Our simulations show that the stellar wind matter intercepted by the compact object might create observational appearances similar to that of an accretion disc corona/wind and could be detected via high energy resolution observations of X-ray absorption lines, in particular, highly ionized ions of heavy elements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Self Regulated Shocks in Massive Star Binary Systems

In an early-type, massive star binary system, X-ray bright shocks result from the powerful collision of stellar winds driven by radiation pressure on spectral line transitions. We examine the influence of the X-rays from the wind-wind collision shocks on the radiative driving of the stellar winds using steady state models that include a parameterized line force with X-ray ionization dependence. Our primary result is that X-ray radiation from the shocks inhibits wind acceleration and can lead to a lower pre-shock velocity, and a correspondingly lower shocked plasma temperature, yet the intrinsic X-ray luminosity of the shocks, LX remains largely unaltered, with the exception of a modest increase at small binary separations. Due to the feedback loop between the ionizing X-rays from the shocks and the wind-driving, we term this scenario as self regulated shocks. This effect is found to greatly increase the range of binary separations at which a wind-photosphere collision is likely to occur in systems where the momenta of the two winds are significantly different. Furthermore, the excessive levels of X-ray ionization close to the shocks completely suppresses the line force, and we suggest that this may render radiative braking less effective. Comparisons of model results against observations reveals reasonable agreement in terms of log(LX/Lbol). The inclusion of self regulated shocks improves the match for kT values in roughly equal wind momenta systems, but there is a systematic offset for systems with unequal wind momenta (if considered to be a wind-photosphere collision).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The 2.35 year itch of Cyg OB2 #9 - I. Optical and X-ray monitoring

Context: Nonthermal radio emission in massive stars is expected to arise in wind-wind collisions occurring inside a binary system. One such case, the O-type star Cyg OB2 #9, was proven to be a binary only four years ago, but the orbital parameters remained uncertain. The periastron passage of 2011 was the first one to be observable under good conditions since the discovery of binarity. Aims: In this context, we have organized a large monitoring campaign to refine the orbital solution and to study the wind-wind collision. Methods: This paper presents the analysis of optical spectroscopic data, as well as of a dedicated X-ray monitoring performed with Swift and XMM. Results: In light of our refined orbital solution, Cyg OB2 #9 appears as a massive O+O binary with a long period and high eccentricity; its components (O5-5.5I for the primary and O3-4III for the secondary) have similar masses and similar luminosities. The new data also provide the first evidence that a wind-wind collision is present in the system. In the optical domain, the broad Ha line varies, displaying enhanced absorption and emission components at periastron. X-ray observations yield the unambiguous signature of an adiabatic collision because, as the stars approach periastron, the X-ray luminosity closely follows the 1/D variation expected in that case. The X-ray spectrum appears, however, slightly softer at periastron, which is probably related to winds colliding at slightly lower speeds at that time. Conclusions: It is the first time that such a variation has been detected in O+O systems, and the first case where the wind-wind collision is found to remain adiabatic even at periastron passage.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

He II 4686 in eta Carinae: collapse of the wind-wind collision region during periastron passage

The periodic spectroscopic events in eta Carinae are now well established and occur near the periastron passage of two massive stars in a very eccentric orbit. Several mechanisms have been proposed to explain the variations of different spectral features, such as an eclipse by the wind-wind collision boundary, a shell ejection from the primary star or accretion of its wind onto the secondary. All of them have problems explaining all the observed phenomena. To better understand the nature of the cyclic events, we performed a dense monitoring of eta Carinae with 5 Southern telescopes during the 2009 low excitation event, resulting in a set of data of unprecedented quality and sampling. The intrinsic luminosity of the He II 4686 emission line (L~310 Lsun) just before periastron reveals the presence of a very luminous transient source of extreme UV radiation emitted in the wind-wind collision (WWC) region. Clumps in the primary's wind probably explain the flare-like behavior of both the X-ray and He II 4686 light-curves. After a short-lived minimum, He II 4686 emission rises again to a ne

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Investigating the X-ray emission from the massive WR+O binary WR 22 using 3D hydrodynamical models

We examine the dependence of the wind-wind collision and subsequent X-ray emission from the massive WR+O star binary WR~22 on the acceleration of the stellar winds, radiative cooling, and orbital motion. Simulations were performed with instantaneously accelerated and radiatively driven stellar winds. Radiative transfer calculations were performed on the simulation output to generate synthetic X-ray data, which are used to conduct a detailed comparison against observations. When instantaneously accelerated stellar winds are adopted in the simulation, a stable wind-wind collision region (WCR) is established at all orbital phases. In contrast, when the stellar winds are radiatively driven, and thus the acceleration regions of the winds are accounted for, the WCR is far more unstable. As the stars approach periastron, the ram pressure of the WR's wind overwhelms the O star's and, following a significant disruption of the shocks by non-linear thin-shell instabilities (NTSIs), the WCR collapses onto the O star. X-ray calculations reveal that when a stable WCR exists the models over-predict the observed X-ray flux by more than two orders of magnitude. The collapse of the WCR onto the O star substantially reduces the discrepancy in the $2-10\;$keV flux to a factor of $\simeq 6$ at $ϕ=0.994$. However, the observed spectrum is not well matched by the models. We conclude that the agreement between the models and observations could be improved by increasing the ratio of the mass-loss rates in favour of the WR star to the extent that a normal wind ram pressure balance does not occur at any orbital phase, potentially leading to a sustained collapse of the WCR onto the O star. Radiative braking may then play a significant rôle for the WCR dynamics and resulting X-ray emission.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Adaptive image ray-tracing for astrophysical simulations

A technique is presented for producing synthetic images from numerical simulations whereby the image resolution is adapted around prominent features. In so doing, adaptive image ray-tracing (AIR) improves the efficiency of a calculation by focusing computational effort where it is needed most. The results of test calculations show that a factor of >~ 4 speed-up, and a commensurate reduction in the number of pixels required in the final image, can be achieved compared to an equivalent calculation with a fixed resolution image.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerical heat conduction in hydrodynamical models of colliding hypersonic flows

Hydrodynamical models of colliding hypersonic flows are presented which explore the dependence of the resulting dynamics and the characteristics of the derived X-ray emission on numerical conduction and viscosity. For the purpose of our investigation we present models of colliding flow with plane-parallel and cylindrical divergence. Numerical conduction causes erroneous heating of gas across the contact discontinuity which has implications for the rate at which the gas cools. We find that the dynamics of the shocked gas and the resulting X-ray emission are strongly dependent on the contrast in the density and temperature either side of the contact discontinuity, these effects being strongest where the postshock gas of one flow behaves quasi-adiabatically while the postshock gas of the other flow is strongly radiative. Introducing additional numerical viscosity into the simulations has the effect of damping the growth of instabilities, which in some cases act to increase the volume of shocked gas and can re-heat gas via sub-shocks as it flows downstream. The resulting reduction in the surface area between adjacent flows, and therefore of the amount of numerical conduction, leads to a commensurate reduction in spurious X-ray emission, though the dynamics of the collision are compromised. The simulation resolution also affects the degree of numerical conduction. A finer resolution better resolves the interfaces of high density and temperature contrast and although numerical conduction still exists the volume of affected gas is considerably reduced. However, since it is not always practical to increase the resolution, it is imperative that the degree of numerical conduction is understood so that inaccurate interpretations can be avoided. This work has implications for the dynamics and emission from astrophysical phenomena which involve high Mach number shocks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spiralling out of control: 3D hydrodynamical modelling of the colliding winds in $η\thinspace$Carinae

Three dimensional (3D) adaptive-mesh refinement (AMR) hydrodynamical simulations of the wind-wind collision between the enigmatic super-massive star \etacar and its mysterious companion star are presented which include radiative driving of the stellar winds, gravity, optically-thin radiative cooling, and orbital motion. Simulations with static stars with a periastron passage separation reveal that the preshock companion star's wind speed is sufficiently reduced that radiative cooling in the postshock gas becomes important, permitting the runaway growth of non-linear thin shell (NTSI) instabilities which massively distort the WCR. However, large-scale simulations which include the orbital motion of the stars, show that orbital motion reduces the impact of radiative inhibition, and thus increases the acquired preshock velocities. As such, the postshock gas temperature and cooling time see a commensurate increase, and sufficient gas pressure is preserved to stabilize the WCR against catastrophic instability growth. We then compute synthetic X-ray spectra and lightcurves and find that, compared to previous models, the X-ray spectra agree much better with {\it XMM-Newton} observations just prior to periastron. The narrow width of the 2009 X-ray minimum can also be reproduced. However, the models fail to reproduce the extended X-ray mimimum from previous cycles. We conclude that the key to explaining the extended X-ray minimum is the rate of cooling of the companion star's postshock wind. If cooling is rapid then powerful NTSIs will heavily disrupt the WCR. Radiative inhibition of the companion star's preshock wind, albeit with a stronger radiation-wind coupling than explored in this work, could be an effective trigger.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The interactions of winds from massive young stellar objects

The supersonic stellar and disk winds possessed by massive young stellar objects will produce shocks when they collide against the interior of a pre-existing bipolar cavity (resulting from an earlier phase of jet activity). The shock heated gas emits thermal X-rays which may be observable by spaceborne observa- tories such as the Chandra X-ray Observatory. Hydrodynamical models are used to explore the wind-cavity interaction. Radiative transfer calculations are performed on the simulation output to produce synthetic X-ray observations, allowing constraints to be placed on model parameters through comparisons with observations. The model reveals an intricate interplay between the inflowing and outflowing material and is successful in reproducing the observed X-ray count rates from massive young stellar objects.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

3D modelling of the colliding winds in Eta Carinae - evidence for radiative inhibition

The X-ray emission from the super-massive star Eta Carinae is simulated using a three dimensional model of the wind-wind collision. In the model the intrinsic X-ray emission is spatially extended and energy dependent. Absorption due to the unshocked stellar winds and the cooled postshock material from the primary LBV star is calculated as the intrinsic emission is ray-traced along multiple sightlines through the 3D spiral structure of the circumstellar environment. The observable emission is then compared to available X-ray data, including the lightcurve observed by the Rossi X-ray Timing Explorer (RXTE) and spectra observed by XMM-Newton. The orientation and eccentricity of the orbit are explored, as are the wind parameters of the stars and the nature and physics of their close approach. Our modelling supports a viewing angle with an inclination of ~ 42 degrees, consistent with the polar axis of the Homunculus nebula (Smith 2006), and the projection of the observer's line-of-sight onto the orbital plane has an angle of ~ 0 - 30 degrees in the prograde direction on the apastron side of the semi-major axis. However, there are significant discrepancies between the observed and model lightcurves and spectra through the X-ray minimum. In particular, the hard flux in our synthetic spectra is an order of magnitude greater than observed. Further calculations reveal that radiative inhibition significantly reduces the preshock velocity of the companion wind. As a consequence the hard X-ray emission is quenched, but it is unclear whether the long duration of the minimum is due solely to this mechanism alone. Models incorporating a collapse/disruption of the WCR and/or reduced preshock companion wind velocities bring the predicted emission and the observations into much better agreement (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

3D models of radiatively driven colliding winds in massive O+O star binaries - III. Thermal X-ray emission

The X-ray emission from the wind-wind collision in short-period massive O+O-star binaries is investigated. The emission is calculated from three-dimensional hydrodynamical models which incorporate gravity, the driving of the winds, orbital motion of the stars, and radiative cooling of the shocked plasma. Changes in the amount of stellar occultation and circumstellar attenuation introduce phase-dependent X-ray variability in systems with circular orbits, while strong variations in the intrinsic emission also occur in systems with eccentric orbits. The X-ray emission in eccentric systems can display strong hysteresis, with the emission softer after periastron than at corresponding orbital phases prior to periastron, reflecting the physical state of the shocked plasma at these times. Furthermore, the rise of the luminosity to maximum does not necessarily follow a 1/D law. Our models further demonstrate that the effective circumstellar column can be highly energy dependent. We simulate Chandra and Suzaku observations, and fit these using standard XSPEC models. We find that the recovered temperatures from two or three-temperature mekal fits are comparable to those from fits to the emission from real systems with similar stellar and orbital parameters/nature. We also find that when the global abundance is thawed in the spectral fits, sub-solar values are exclusively returned, despite the calculations using solar values as input (abridged).

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The interactions of winds from massive young stellar objects: X-ray emission, dynamics, and cavity evolution

2D axis-symmetric hydrodynamical simulations are presented which explore the interaction of stellar and disk winds with surrounding infalling cloud material. The star, and its accompanying disk, blow winds inside a cavity cleared out by an earlier jet. The collision of the winds with their surroundings generates shock heated plasma which reaches temperatures up to ~10^8 K. Attenuated X-ray spectra are calculated from solving the equation of radiative transfer along lines-of-sight. This process is repeated at various epochs throughout the simulations to examine the evolution of the intrinsic and attenuated flux. We find that the dynamic nature of the wind-cavity interaction fuels intrinsic variability in the observed emission on timescales of several hundred years. This is principally due to variations in the position of the reverse shock which is influenced by changes in the shape of the cavity wall. The collision of the winds with the cavity wall can cause clumps of cloud material to be stripped away. Mixing of these clumps into the winds mass-loads the flow and enhances the X-ray emission measure. The position and shape of the reverse shock plays a key role in determining the strength and hardness of the X-ray emission. In some models the reverse shock is oblique to much of the stellar and disk outflows, whereas in others it is closely normal over a wide range of polar angles. For reasonable stellar and disk wind parameters the integrated count rate and spatial extent of the intensity peak for X-ray emission agree with \textit{Chandra} observations of the deeply embedded MYSOs S106 IRS4, Mon R2 IRS3 A, and AFGL 2591.(abridged)

0 citations0 reviewsNo code linkedOpen access