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

J. M. Pittard

J. M. Pittard appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.HEastro-phastro-ph.COastro-ph.IM

Catalog footprint

What is connected

40works
6topics
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

40 published item(s)

preprint2022arXiv

Momentum and energy injection by a wind-blown bubble into an inhomogeneous interstellar medium

We investigate the effect of mass-loading from embedded clouds on the evolution of wind-blown bubbles. We use 1D hydrodynamical calculations and assume that the clouds are numerous enough that they can be treated in the continuous limit, and that rapid mixing occurs so that the injected mass quickly merges with the global flow. The destruction of embedded clouds adds mass into the bubble, increasing its density. Mass-loading increases the temperature of the unshocked stellar wind due to the frictional drag, and reduces the temperature of the hot shocked gas as the available thermal energy is shared between more particles. Mass-loading may increase or decrease the volume-averaged bubble pressure. Mass-loaded bubbles are smaller, have less retained energy and lower radial momentum, but in all cases examined are still able to do significant $PdV$ work on the swept-up gas. In this latter respect, the bubbles more closely resemble energy-conserving bubbles than the momentum-conserving-like behaviour of ``quenched'' bubbles.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shocking interactions of supernova remnants with atomic and molecular clouds -- the interplay between shocks, thermal instability and gravity in the large cloud regime

Using the adaptive mesh refinement code MG, we perform 3D hydrodynamic simulations of a supernova-cloud interaction in the "large cloud regime". The cloud is initially atomic and evolving due to the thermal instability (TI) and gravity. We study interactions in a "pre-TI" and "post-TI" stage when cold and dense clumps are present, and compare these results to idealised shock-cloud scenarios in the "small cloud regime", and a scenario without shocks. On aggregate, the supernova disruption is significantly weaker than that from an idealised shock due to the supernova impact being instantaneous, and not continuous. In both supernova-cloud interactions, we observe two shocks impact the cloud, followed by the development of a weak 10 km s$^{-1}$ upstream flow on the cloud interface, and a global ambient pressure drop. When the cloud is still atomic, it expands due to this drop. Additionally, the TI is triggered at the front of the cloud, causing the formation of a cap-like structure with clumps embedded inside. The upstream flow converges in this region, resulting in a lobe-like cloud morphology. When the cloud is molecular, the transmitted shock disrupts the inter-clump material and causes the clumps' outer envelopes to expand slightly and form tail-like morphologies. These effects are less pronounced than those in our shock-cloud scenarios, and more pronounced that those in our un-shocked scenario. After 3.5 Myrs, the effects from the supernova decay and the cloud returns to an almost indistinguishable state from an un-shocked cloud, in spite of the global ambient pressure drop. In neither supernova-cloud scenario do we see any local gravitational collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

How D-type HII region expansion depends on numerical resolution

We investigate the resolution dependence of HII regions expanding past their Strömgren spheres. We find that their structure and size, and the radial momentum that they attain at a given time, is in good agreement with analytical expectations if the Strömgren radius is resolved with $dr \leq 0.3\,R_{\rm st}$. If this is not satisfied, the radial momentum may be over- or under-estimated by factors up to 10 or more. Our work has significance for the amount of radial momentum that a HII region can impart to the ambient medium in numerical simulations, and thus on the relative importance of ionizing feedback from massive stars.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

COBRaS: The e-MERLIN 21 cm Legacy survey of Cygnus OB2

The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Colliding-wind binary systems: Diffusive shock acceleration and non-thermal emission

We present a model for the non-thermal emission from a colliding-wind binary. Relativistic protons and electrons are assumed to be accelerated through diffusive shock acceleration (DSA) at the global shocks bounding the wind-wind collision region. The non-linear effects of the back-reaction due to the cosmic ray pressure on the particle acceleration process and the cooling of the non-thermal particles as they flow downstream from the shocks are included. We explore how the non-thermal particle distribution and the keV-GeV emission changes with the stellar separation and the viewing angle of the system, and with the momentum ratio of the winds. We confirm earlier findings that DSA is very efficient when magnetic field amplification is not included, leading to significantly modified shocks. We also find that the non-thermal flux scales with the binary separation in a complicated way and that the anisotropic inverse Compton emission shows only a moderate variation with viewing angle due to the spatial extent of the wind-wind collision.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interactions of a shock with a molecular cloud at various stages of its evolution due to thermal instability and gravity

Using the adaptive mesh refinement code MG, we perform hydrodynamic simulations of the interaction of a shock with a molecular cloud evolving due to thermal instability and gravity. To explore the relative importance of these processes, three case studies are presented. The first follows the formation of a molecular cloud out of an initially quiescent atomic medium due to the effects of thermal instability and gravity. The second case introduces a shock whilst the cloud is still in the warm atomic phase, and the third scenario introduces a shock once the molecular cloud has formed. The shocks accelerate the global collapse of the clouds with both experiencing local gravitational collapse prior to this. When the cloud is still atomic, the evolution is shock dominated and structures form due to dynamical instabilities within a radiatively cooled shell. While the transmitted shock can potentially trigger the thermal instability, this is prevented as material is shocked multiple times on the order of a cloud crushing time-scale. When the cloud is molecular, the post-shock flow is directed via the pre-existing structure through low-density regions in the inter-clump medium. The clumps are accelerated and deformed as the flow induces clump-clump collisions and mergers that collapse under gravity. For a limited period, both shocked cases show a mixture of Kolmogorov and Burgers turbulence-like velocity and logarithmic density power spectra, and strongly varying density spectra. The clouds presented in this work provide realistic conditions that will be used in future feedback studies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Search for non-thermal X-ray emission in the colliding wind binary Cyg OB2 #8A

Cyg OB2 #8a is a massive O-type binary displaying strong non-thermal radio emission. Owing to the compactness of this binary, emission of non-thermal X-ray photons via inverse Compton scattering is expected. We first revised the orbital solution for Cyg OB2 #8a using new optical spectra. We then reduced and analysed X-ray spectra obtained with XMM-Newton, Swift, INTEGRAL, and NuSTAR. The analysis of the XMM-Newton and Swift data allows us to better characterise the X-ray emission from the stellar winds and colliding winds region at energies below 10 keV. We confirm the variation of the broad-band light curve of Cyg OB2 #8a along the orbit with, for the first time, the observation of the maximum emission around phase 0.8. The minimum ratio of the X-ray to bolometric flux of Cyg OB2 #8a remains well above the level expected for single O-type stars, indicating that the colliding wind region is not disrupted during the periastron passage. The analysis of the full set of publicly available INTEGRAL observations allows us to refine the upper limit on the non-thermal X-ray flux of the Cyg OB2 region between 20 and 200 keV. Two NuSTAR observations (phases 0.028 and 0.085) allow us to study the Cyg OB2 #8a spectrum up to 30 keV. These data do not provide evidence of the presence of non-thermal X-rays, but bring more stringent constraints on the flux of a putative non-thermal component. Finally, we computed, thanks to a new dedicated model, the anisotropic inverse Compton emission generated in the wind shock region. The theoretical non-thermal emission appears to be compatible with observational limits and the kinetic luminosity computed from these models is in good agreement with the unabsorbed flux observed below 10 keV.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The isothermal evolution of a shock-filament interaction

Studies of filamentary structures that are prevalent throughout the interstellar medium are of great significance to a number of astrophysical fields. Here, we present 3D hydrodynamic simulations of shock-filament interactions where the equation of state has been softened to become almost isothermal. We investigate the effect of such an isothermal regime on the interaction (where both the shock and filament are isothermal), and we examine how the nature of the interaction changes when the orientation of the filament, the shock Mach number, and the filament density contrast are varied. We find that only sideways-oriented filaments with a density contrast of $10^2$ form a three-rolled structure, dissimilar to the results of a previous study. Moreover, the angle of orientation of the filament plays a large role in the evolution of the filament morphology: the greater the angle of orientation, the longer and less turbulent the wake. Turbulent stripping of filament material leading to fragmentation of the core occurs in most filaments; however, filaments orientated at an angle of $85^{\circ}$ to the shock front do not fragment and are longer-lived. In addition, values of the drag time are influenced by the filament length, with longer filaments being accelerated faster than shorter ones. Furthermore, filaments in an isothermal regime exhibit faster acceleration than those struck by an adiabatic shock. Finally, we find that the drag and mixing times of the filament increase as the angle of orientation of the filament is increased.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

MHD simulation of the formation of clumps and filaments in quiescent diffuse medium by thermal instability

We have used the AMR hydrodynamic code, MG, to perform idealised 3D MHD simulations of the formation of clumpy and filamentary structure in a thermally unstable medium without turbulence. A stationary thermally unstable spherical diffuse atomic cloud with uniform density in pressure equilibrium with low density surroundings was seeded with random density variations and allowed to evolve. A range of magnetic field strengths threading the cloud have been explored, from beta=0.1 to beta=1.0 to the zero magnetic field case (beta=infinity), where beta is the ratio of thermal pressure to magnetic pressure. Once the density inhomogeneities had developed to the point where gravity started to become important, self-gravity was introduced to the simulation. With no magnetic field, clouds and clumps form within the cloud with aspect ratios of around unity, whereas in the presence of a relatively strong field (beta=0.1) these become filaments, then evolve into interconnected corrugated sheets that are predominantly perpendicular to the magnetic field. With magnetic and thermal pressure equality (beta=1.0), filaments, clouds and clumps are formed. At any particular instant, the projection of the 3D structure onto a plane parallel to the magnetic field, i.e. a line of sight perpendicular to the magnetic field, resembles the appearance of filamentary molecular clouds. The filament densities, widths, velocity dispersions and temperatures resemble those observed in molecular clouds. In contrast, in the strong field case beta=0.1, projection of the 3D structure along a line of sight parallel to the magnetic field reveals a remarkably uniform structure.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical Simulations of a Shock-Filament Interaction

We present 3D hydrodynamic adiabatic simulations of a shock interacting with a dense, elongated cloud. We compare how the nature of the interaction changes with the filament's length and its orientation to the shock, and with the shock Mach number and the density contrast of the filament. We then examine the differences with respect to 3D spherical-cloud calculations. We find significant differences in the morphology of the interaction when M=10 and chi=100: in many cases 3 parallel rolls are formed, and spread further apart with time, and periodic vortex shedding can occur off the ends of oblique filaments. Sideways-on filaments are accelerated more quickly, and initially lose mass more quickly than spherical clouds due to their greater surface area to volume ratio. However, at late stages they lose mass more slowly, due to the reduced relative speed between the filament and the postshock flow. The acceleration and mixing timescales can vary by a factor of 2 as the filament orientation changes. Oblique filaments can achieve transverse velocities up to 10% of the shock speed. Some aspects of our simulations are compared against experimental and numerical work on rigid cylinders.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The interaction of a magnetohydrodynamical shock with a filament

We present 3D magnetohydrodynamic numerical simulations of the adiabatic interaction of a shock with a dense, filamentary cloud. We investigate the effects of various filament lengths and orientations on the interaction using different orientations of the magnetic field, and vary the Mach number of the shock, the density contrast of the filament, and the plasma beta, in order to determine their effect on the evolution and lifetime of the filament. We find that in a parallel magnetic field filaments have longer lifetimes if they are orientated more 'broadside' to the shock front, and that an increase in the density contrast hastens the destruction of the cloud, in terms of the modified cloud-crushing time-scale, tcs. The combination of a mild shock and a perpendicular or oblique field provides the best condition for extending the life of the filament, with some filaments able to survive almost indefinitely since they are cocooned by the magnetic field. A high value for the density contrast does not initiate large turbulent instabilities in either the perpendicular or oblique field cases but rather draws the filament out into long tendrils which may eventually fragment. In addition, flux ropes are only formed in parallel magnetic fields. The length of the filament is, however, not as important for the evolution and destruction of a filament.

0 citations0 reviewsNo code linkedOpen access
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
preprint2014arXiv

Feedback from Winds and Supernovae in Massive Stellar Clusters - II. X-Ray Emission

The X-ray emission from a simulated massive stellar cluster is investigated. The emission is calculated from a 3D hydrodynamical model which incorporates the mechanical feedback from the stellar winds of 3 O-stars embedded in a giant molecular cloud clump containing 3240 M$_{\odot}$ of molecular material within a 4 pc radius. A simple prescription for the evolution of the stars is used, with the first supernova explosion at t=4.4 Myrs. We find that the presence of the GMC clump causes short-lived attenuation effects on the X-ray emission of the cluster. However, once most of the material has been ablated away by the winds the remaining dense clumps do not have a noticable effect on the attenuation compared with the assumed interstellar medium column. We determine the evolution of the cluster X-ray luminosity, L$_X$, and spectra, and generate synthetic images. The intrinsic X-ray luminosity drops from nearly 10$^{34}$ ergs s$^{-1}$ while the winds are `bottled up', to a near constant value of 1.7$\times 10^{32}\rm ergs s^{-1}$ between t=1-4 Myrs. L$_X$ reduces slightly during each star's red supergiant stage due to the depressurization of the hot gas. However, L$_X$ increases to $\approx 10^{34}\rm\,ergs s^{-1}$ during each star's Wolf-Rayet stage. The X-ray luminosity is enhanced by 2-3 orders of magnitude to $\sim 10^{37}\rm ergs s^{-1}$ for at least 4600 yrs after each supernova, at which time the blast wave leaves the grid and the X-ray luminosity drops. The X-ray luminosity of our simulation is generally considerably fainter than predicted from spherically-symmetric bubble models, due to the leakage of hot gas material through gaps in the outer shell. This process reduces the pressure within our simulation and thus the X-ray emission. However, the X-ray luminosities and temperatures which we obtain are comparable to similarly powerful massive young clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Simulations of a Shock Interacting with Multiple Magnetized Clouds

We present 2D adiabatic magnetohydrodynamic (MHD) simulations of a shock interacting with groups of two or three cylindrical clouds. We study how the presence of a nearby cloud influences the dynamics of this interaction, and explore the resulting differences and similarities in the evolution of each cloud. The understanding gained from this small-scale study will help to interpret the behaviour of systems with many 10's or 100's of clouds. We observe a wide variety of behaviour in the interactions studied, which is dependent on the initial positions of the clouds and the orientation and strength of the magnetic field. We find: i) some clouds are stretched along their field-lines, whereas others are confined by their field-lines; ii) upstream clouds may accelerate past downstream clouds (though magnetic tension can prevent this); iii) clouds may also change their relative positions transverse to the direction of shock propagation as they "slingshot" past each other; iv) downstream clouds may be offered some protection from the oncoming flow as a result of being in the lee of an upstream cloud; v) the cycle of cloud compression and re-expansion is generally weaker when there are nearby neighbouring clouds; vi) the plasma $β$ in cloud material can vary rapidly as clouds collide with one another, but low values of $β$ are always transitory. This work is relevant to studies of multi-phase regions, where fast, low-density gas interacts with dense clouds, such as in circumstellar bubbles, supernova remnants, superbubbles and galactic winds.

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

Feedback from Winds and Supernovae in Massive Stellar Clusters. I: Hydrodynamics

We use 3D hydrodynamical models to investigate the effects of massive star feedback from winds and supernovae on inhomogeneous molecular material left over from the formation of a massive stellar cluster. We simulate the interaction of the mechanical energy input from a cluster with 3 O-stars into a giant molecular cloud (GMC) clump containing 3240 solar masses of molecular material within a 4 pc radius. The cluster wind blows out of the molecular clump along low-density channels, into which denser clump material is entrained. We find that the densest molecular regions are surprisingly resistant to ablation by the cluster wind, in part due to shielding by other dense regions closer to the cluster. Nonetheless, molecular material is gradually removed by the cluster wind during which mass-loading factors in excess of several 100 are obtained. Because the clump is very porous, 60-75 per cent of the injected wind energy escapes the simulation domain, with the difference being radiated. After 4.4 Myr, the massive stars in our simulation begin to explode as supernovae. The highly structured environment into which the SN energy is released allows even weaker coupling to the remaining dense material and practically all of the SN energy reaches the wider environment. The molecular material is almost completely dispersed and destroyed after 6 Myr. The escape fraction of ionizing radiation is estimated to be about 50 per cent during the first 4 Myr of the cluster's life. A similar model with a larger and more massive GMC clump reveals the same general picture, though more time is needed for it to be destroyed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Ram pressure stripping of the hot gaseous haloes of galaxies using the k-epsilon sub-grid turbulence model

We perform three dimensional hydrodynamic simulations of the ram pressure stripping of the hot extended gaseous halo of a massive galaxy using the k-epsilon sub-grid turbulence model at Mach numbers 0.9, 1.1 and 1.9. The k-epsilon model is used to simulate high Reynolds number flows by increasing the transport coefficients in regions of high turbulence. We find that the initial, instantaneous stripping is the same whether or not the k-epsilon model is implemented and is in agreement with the results of other studies. However the use of the k-epsilon model leads to five times less gas remaining after stripping by a supersonic flow has proceeded for 10 Gyr, which is more consistent with what simple analytic calculations indicate. Hence the continual Kelvin-Helmholtz stripping plays a significant role in the ram pressure stripping of the haloes of massive galaxies. To properly account for this, simulations of galaxy clusters will require the use of sub-grid turbulence models

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical simulations of shocks encountering clumpy regions

We present numerical simulations of the adiabatic interaction of a shock with a clumpy region containing many individual clouds. Our work incorporates a sub-grid turbulence model which for the first time makes this investigation feasible. We vary the Mach number of the shock, the density contrast of the clouds, and the ratio of total cloud mass to inter-cloud mass within the clumpy region. Cloud material becomes incorporated into the flow. This "mass-loading" reduces the Mach number of the shock, and leads to the formation of a dense shell. In cases in which the mass-loading is sufficient the flow slows enough that the shock degenerates into a wave. The interaction evolves through up to four stages: initially the shock decelerates; then its speed is nearly constant; next the shock accelerates as it leaves the clumpy region; finally it moves at a constant speed close to its initial speed. Turbulence is generated in the post-shock flow as the shock sweeps through the clumpy region. Clouds exposed to turbulence can be destroyed more rapidly than a similar cloud in an "isolated" environment. The lifetime of a downstream cloud decreases with increasing cloud-to-intercloud mass ratio. We briefly discuss the significance of these results for starburst superwinds and galaxy evolution.

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

Cosmic Rays and Radiative Instabilities

In the absence of magnetic fields and cosmic rays, radiative cooling laws with a range of dependences on temperature affect the stability of interstellar gas. For about four and a half decades, astrophysicists have recognised the importance of the thermal instablity for the formation of clouds in the interstellar medium. Even in the past several years, many papers have concerned the role of the thermal instability in the production of molecular clouds. About three and a half decades ago, astrophysicists investigating radiative shocks noticed that for many cooling laws such shocks are unstable. Attempts to address the effects of cosmic rays on the stablity of radiative media that are initially uniform or that have just passed through shocks have been made. The simplest approach to such studies involves the assumption that the cosmic rays behave as a fluid. Work based on such an approach is described. Cosmic rays have no effect on the stability of initially uniform, static media with respect to isobaric perturbations, though they do affect the stability of such media with respect to isentropic perturbations. The effect of cosmic rays on the stability of radiative shocked media depends greatly on the efficiency of the conversion of energy in accelerated cosmic rays into thermal energy in the thermalized fluid. If that efficiency is low, radiative cooling makes weak shocks propagating into upstream media with low cosmic-ray pressures more likely to be cosmic-ray dominated than adiabatic shocks of comparable strength. The cosmic-ray dominated shocks do not display radiative overstability. Highly efficient conversion of cosmic-ray energy into thermal energy leads shocked media to behave as they do when cosmic rays are absent.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Global X-ray properties of the O and B stars in Carina

The key empirical property of the X-ray emission from O stars is a strong correlation between the bolometric and X-ray luminosities. In the framework of the Chandra Carina Complex Project, 129 O and B stars have been detected as X-ray sources; 78 of those, all with spectral type earlier than B3, have enough counts for at least a rough X-ray spectral characterization. This leads to an estimate of the Lx/Lbol ratio for an exceptional number of 60 O stars belonging to the same region and triples the number of Carina massive stars studied spectroscopically in X-rays. The derived log(Lx/Lbol) is -7.26 for single objects, with a dispersion of only 0.21dex. Using the properties of hot massive stars listed in the literature, we compare the X-ray luminosities of different types of objects. In the case of O stars, the Lx/Lbol ratios are similar for bright and faint objects, as well as for stars of different luminosity classes or spectral types. Binaries appear only slightly harder and slightly more luminous in X-rays than single objects; the differences are not formally significant (at the 1% level), except for the Lx/Lbol ratio in the medium (1.0--2.5keV) energy band. Weak-wind objects have similar X-ray luminosities but they display slightly softer spectra compared to "normal" O stars with the same bolometric luminosity. Discarding three overluminous objects, we find a very shallow trend of harder emission in brighter objects. The properties of the few B stars bright enough to yield some spectral information appear to be different overall (constant X-ray luminosities, harder spectra), hinting that another mechanism for producing X-rays, besides wind shocks, might be at work. However, it must be stressed that the earliest and X-ray brightest amongst these few detected objects are similar to the latest O stars, suggesting a possibly smooth transition between the two processes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The WR 140 periastron passage 2009: first results from MONS and other optical sources

We present the results from the spectroscopic follow-up of WR140 (WC7 + O4-5) during its last periastron passage in January 2009. This object is known as the archetype of colliding wind binaries and has a relatively large period (~ 8 years) and eccentricity (~ 0.89). We provide updated values for the orbital parameters, new estimates for the WR and O star masses and new constraints on the mass-loss rates.

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

Tails of the Unexpected: The Interaction of an Isothermal Shell with a Cloud

A new mechanism for the formation of cometary tails behind dense clouds or globules is discussed. Numerical hydrodynamical models show that when a dense shell of swept-up matter overruns a cloud, material in the shell is focussed behind the cloud to form a tail. This mode of tail formation is completely distinct from other methods, which involve either the removal of material from the cloud, or shadowing from a strong, nearby source of ionization. This mechanism is relevant to the cometary tails seen in planetary nebulae and to the interaction of superbubble shells with dense clouds.

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
preprint2010arXiv

The Turbulent Destruction of Clouds - II. Mach Number Dependence, Mass-loss Rates, and Tail Formation

The turbulent destruction of a cloud subject to the passage of an adiabatic shock is studied. We find large discrepancies between the lifetime of the cloud and the analytical result of Hartquist et al. (1986). These differences appear to be due to the assumption in Hartquist et al. that mass-loss occurs largely as a result of lower pressure regions on the surface of the cloud away from the stagnation point, whereas in reality Kelvin-Helmholtz (KH) instabilities play a dominant role in the cloud destruction. We find that the true lifetime of the cloud (defined as when all of the material from the core of the cloud is well mixed with the intercloud material in the hydrodynamic cells) is about 6 times t_KHD, where t_KHD is the growth timescale for the most disruptive, long-wavelength, KH instabilities. These findings have wide implications for diffuse sources where there is transfer of material between hot and cool phases. The properties of the interaction as a function of Mach number and cloud density contrast are also studied. The interaction is milder at lower Mach numbers with the most marked differences occuring at low shock Mach numbers when the postshock gas is subsonic with respect to the cloud (i.e. M < 2.76). Material stripped off the cloud only forms a long "tail-like" feature if the density contrast of the cloud to the ambient medium, chi > 1e3.

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 - II. Thermal radio to sub-mm emission

In this work the thermal emission over cm to sub-mm wavelengths from the winds in short-period O+O-star binaries is investigated (potential non-thermal emission is presently ignored). The calculations are based on three-dimensional hydrodynamical models which incorporate gravity, the driving of the winds, orbital motion of the stars, and radiative cooling of the shocked plasma. The thermal emission arises from the stellar winds and the region where they collide. We investigate the flux and spectrum from a variety of models as a function of orbital phase and orientation of the observer, and compare to the single star case. The emission from the wind-wind collision region (WCR) is strongly dependent on its density and temperature, being optically thick in radiative systems, and optically thin in adiabatic systems. The flux from systems where the WCR is highly radiative, as investigated for the first time in this work, can be over an order of magnitude greater than the combined flux from identically typed single stars. This excess occurs over a broad range of wavelengths from cm to sub-mm. In contrast, when the WCR is largely adiabatic, a significant excess in the thermal flux occurs only below 100 GHz. Eccentric systems may show order of magnitude or greater flux variability, especially if the plasma in the WCR transitions from an adiabatic to a radiative state and vice-versa - in such cases the flux can display significant hysteresis with stellar separation. We further demonstrate that clumping can affect the variability of radio lightcurves. We investigate the spectral index of the emission, and often find indices steeper than +0.6. Our predictions are of interest to future observations with the next generation of radio and sub-mm telescopes (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

3D Models of Radiatively Driven Colliding Winds In Massive O+O Star Binaries: I. Hydrodynamics

The dynamics of the wind-wind collision in massive stellar binaries is investigated using three-dimensional hydrodynamical models which incorporate gravity, the driving of the winds, the orbital motion of the stars, and radiative cooling of the shocked plasma. In this first paper we restrict our study to main-sequence O+O binaries. The nature of the wind-wind collision region is highly dependent on the degree of cooling of the shocked plasma, and the ratio of the flow timescale of the shocked plasma to the orbital timescale. The pre-shock wind speeds are lower in close systems as the winds collide prior to their acceleration to terminal speeds. Radiative inhibition may also reduce the pre-shock wind speeds. Together, these effects can lead to rapid cooling of the post-shock gas. Radiative inhibition is less important in wider systems, where the winds are accelerated to higher speeds before they collide, and the resulting collision region can be largely adiabatic. In systems with eccentric orbits, cold gas formed during periastron passage can persist even at apastron, before being ablated and mixed into its surroundings and/or accelerated out of the system.

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
preprint2008arXiv

Chandra X-ray Grating Spectrometry of Eta Carinae near X-ray Minimum: I. Variability of the Sulfur and Silicon Emission Lines

We report on variations in important X-ray emission lines in a series of Chandra grating spectra of the supermassive colliding wind binary star Eta Carinae, including key phases around the X-ray minimum/periastron passage in 2003.5. The X-rays arise from the collision of the slow, dense wind of Eta Car with the fast, low-density wind of an otherwise hidden companion star. The X-ray emission lines provide the only direct measure of the flow dynamics of the companion's wind along the wind-wind collision zone. We concentrate here on the silicon and sulfur lines, which are the strongest and best resolved lines in the X-ray spectra. Most of the line profiles can be adequately fit with symmetric Gaussians with little significant skewness. Both the silicon and sulfur lines show significant velocity shifts and correlated increases in line widths through the observations. The R = forbidden-to-intercombination ratio from the Si XIII and S XV triplets is near or above the low-density limit in all observations, suggesting that the line-forming region is >1.6 stellar radii from the companion star. We show that simple geometrical models cannot simultaneously fit both the observed centroid variations and changes in line width as a function of phase. We show that the observed profiles can be fitted with synthetic profiles with a reasonable model of the emissivity along the wind-wind collision boundary. We use this analysis to help constrain the line formation region as a function of orbital phase, and the orbital geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Mass-Loaded Flows

A key process within astronomy is the exchange of mass, momentum, and energy between diffuse plasmas in many types of astronomical sources (including planetary nebulae, wind-blown bubbles, supernova remnants, starburst superwinds, and the intracluster medium) and dense, embedded clouds or clumps. This transfer affects the large scale flows of the diffuse plasmas as well as the evolution of the clumps. I review our current understanding of mass-injection processes, and examine intermediate-scale structure and the global effect of mass-loading on a flow. I then discuss mass-loading in a variety of diffuse sources.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Non-thermal X-ray and Gamma-ray Emission from the Colliding Wind Binary WR140

WR140 is the archetype long-period colliding wind binary (CWB) system, and is well known for dramatic variations in its synchrotron emission during its 7.9-yr, highly eccentric orbit. This emission is thought to arise from relativistic electrons accelerated at the global shocks bounding the wind-collision region (WCR). The presence of non-thermal electrons and ions should also give rise to X-ray and gamma-ray emission from several separate mechanisms, including inverse-Compton cooling, relativistic bremsstrahlung, and pion decay. We describe new calculations of this emission and make some preliminary predictions for the new generation of gamma-ray observatories. We determine that WR140 will likely require several Megaseconds of observation before detection with INTEGRAL, but should be a reasonably strong source for GLAST.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

An XMM-Newton observation of the massive binary HD159176

We report the analysis of an XMM-Newton observation of the close binary HD159176 (O7V + O7V). The observed L_X/L_bol ratio reveals an X-ray luminosity exceeding by a factor ~ 7 the expected value for X-ray emission from single O-stars, therefore suggesting a wind-wind interaction scenario. EPIC and RGS spectra are fitted consistently with a two temperature mekal optically thin thermal plasma model, with temperatures ranging from ~ 2 to 6 10^6 K. At first sight, these rather low temperatures are consistent with the expectations for a close binary system where the winds collide well before reaching their terminal velocities. We also investigate the variability of the X-ray light curve of HD159176 on various short time scales. No significant variability is found and we conclude that if hydrodynamical instabilities exist in the wind interaction region of HD159176, they are not sufficient to produce an observable signature in the X-ray emission. Hydrodynamic simulations using wind parameters from the literature reveal some puzzling discrepancies. The most striking one concerns the predicted X-ray luminosity which is one or more orders of magnitude larger than the observed one. A significant reduction of the mass loss rate of the components compared to the values quoted in the literature alleviates the discrepancy but is not sufficient to fully account for the observed luminosity. Because hydrodynamical models are best for the adiabatic case whereas the colliding winds in HD159176 are most likely highly radiative, a totally new approach has been envisaged... (see paper for complete abstract)

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

The Evolution of Mass Loaded Supernova Remnants. II. Temperature Dependent Mass Injection Rates

We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~10 million K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Mass-Loss Rate Determination for the Massive Binary V444 Cyg using 3-D Monte-Carlo Simulations of Line and Polarization Variability

A newly developed 3-D Monte Carlo model is used, in conjunction with a multi-line non-LTE radiative transfer model, to determine the mass-loss rate of the Wolf-Rayet (W-R) star in the massive binary \object{V444 Cyg} (WN5+O6). This independent estimate of mass-loss rate is attained by fitting the observed \HeI (5876) Åand \HeII (5412) Åline profiles, and the continuum light curves of three Stokes parameters ((I, Q, U)) in the (V) band simultaneously. The high accuracy of our determination arises from the use of many observational constraints, and the sensitivity of the continuum polarization to the mass-loss rate. Our best fit model suggests that the mass-loss rate of the system is (\dot{M}_{\WR}=0.6(\pm 0.2) \times 10^{-5} M_{\sun} \mathrm{yr}^{-1} ), and is independent of the assumed distance to \object{V444 Cyg}. The fits did not allow a unique value for the radius of the W-R star to be derived. The range of the volume filling factor for the W-R star atmosphere is estimated to be in the range of 0.050 (for $R_{\WR}=5.0 R_{\sun}$) to 0.075 (for $R_{\WR}=2.5 R_{\sun}$). We also found that the blue-side of \HeI (5876 ) Åand \HeII (5412) Ålines at phase 0.8 is relatively unaffected by the emission from the wind-wind interaction zone and the absorption by the O-star atmosphere; hence, the profiles at this phase are suitable for spectral line fittings using a spherical radiative transfer model.

0 citations0 reviewsNo code linkedOpen access