Showing 1–50 of 105 results for author: Pittard, J M

First spatially resolved Na I and He I transitions towards an MYSO. Finding new tracers for the gaseous star/disc interface

Authors: Evgenia Koumpia, M. Koutoulaki, W. -J. de Wit, R. D. Oudmaijer, A. J. Frost, S. L. Lumsden, J. M. Pittard

Abstract: With steady observational advances, the formation of massive stars is being understood in more detail. Numerical models are converging on a scenario where accretion discs play a key role. Direct observational evidence of such discs at a few au scales is scarce, due to the rarity of such objects and the observational challenges, including the lack of adequate diagnostic lines in the near-IR. We pre… ▽ More With steady observational advances, the formation of massive stars is being understood in more detail. Numerical models are converging on a scenario where accretion discs play a key role. Direct observational evidence of such discs at a few au scales is scarce, due to the rarity of such objects and the observational challenges, including the lack of adequate diagnostic lines in the near-IR. We present the analysis of K-band spectro-interferometric observations toward the Massive Young Stellar Object IRAS 13481-6124, which is known to host an accreting dusty disc. Using GRAVITY on the VLTI, we trace the crucial au-scales of the warm inner interface between the star and the accretion dusty disc. We detect and spatially resolve the Na I doublet and He I transitions towards an object of this class for the first time. The new observations in combination with our geometric models allowed us to probe the smallest au-scales of accretion/ejection around an MYSO. We find that Na I originates in the disc at smaller radii than the dust disc and is more compact than any of the other spatially resolved diagnostics (Br$γがんま$, He I, and CO). Our findings suggest that Na I can be a new powerful diagnostic line in tracing the warm star/disc accreting interface of forming (massive) stars, while the similarities between He I and Br$γがんま$ point towards an accretion/ejection origin of He I △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS Letters; 5 pages, 3 Figures

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

Authors: J. M. Pittard

Abstract: 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, increasi… ▽ More 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. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 15 pages, 7 figures, accepted by MNRAS

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

Authors: M. M. Kupilas, J. M. Pittard, C. J. Wareing, S. A. E. G. Falle

Abstract: 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… ▽ More 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. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: 15 pages, 9 figures

Journal ref: MNRAS 513, 2022, 3345

Exploring dust growth in the episodic WCd system WR140

Authors: J. W. Eatson, J. M. Pittard, S. Van Loo

Abstract: The wind collision region (WCR) in a colliding wind binary (CWB) is a particularly violent place, as such, it is surprising that it is also a region where significant quantities of interstellar dust can form. In extreme cases, approximately 30% of the total mass loss rate of a system can be converted into dust. These regions are poorly understood, as observation and simulation of these systems are… ▽ More The wind collision region (WCR) in a colliding wind binary (CWB) is a particularly violent place, as such, it is surprising that it is also a region where significant quantities of interstellar dust can form. In extreme cases, approximately 30% of the total mass loss rate of a system can be converted into dust. These regions are poorly understood, as observation and simulation of these systems are difficult. In our previous paper we simulated dust growth in CWB systems using an advected scalar model and found our model to be suitable for qualitative study. For this paper we simulated the periodic dust forming CWB (WCd) system WR140 with our dust model, to determine how dust growth changes over the systems periastron passage. We found that dust production increases significantly at periastron passage, which is consistent with IR emission of the surrounding dusty shell. We also find that the dust production rate of the system decreases rapidly as the stars recede from each other, though the rate of decrease is significantly lower than the rate of increase during periastron passage. This was found to be due to strong cooling and its associated thermal instabilities, resulting in cool, high-density pockets of gas in the WCR where dust forms. The WCR also shows a degree of hysteresis, resulting in a radiative post-shock flow even when the stars are separated enough for the region to behave adiabatically. △ Less

Submitted 17 October, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: 9 pages, 9 figures, 5 tables

An exploration of dust grain growth within WCd systems using an advected scalar dust model

Authors: J. W. Eatson, J. M. Pittard, S. Van Loo

Abstract: Dust production is one of the more curious phenomena observed in massive binary systems with interacting winds. The high temperatures, UV photon flux and violent shocks should destroy any dust grains that condense. However, in some extreme cases dust production yields of approximately 30% of the total mass loss rate of the stellar winds have been observed. In order to better understand this phenom… ▽ More Dust production is one of the more curious phenomena observed in massive binary systems with interacting winds. The high temperatures, UV photon flux and violent shocks should destroy any dust grains that condense. However, in some extreme cases dust production yields of approximately 30% of the total mass loss rate of the stellar winds have been observed. In order to better understand this phenomenon a parameter space exploration was performed using a series of numerical models of dust producing carbon phase Wolf-Rayet (WCd) systems. These models incorporated a passive scalar dust model simulating dust growth, destruction and radiative cooling. We find that reasonable dust yields were produced by these simulations. Significant changes in the dust yield were caused by changing the mass loss rates of the stars, with a greater mass loss rate contributing to increased dust yields. Similarly, a close orbit between the stars also resulted in higher dust yields. Finally, a high velocity wind shear, which induces Kelvin-Helmholtz (KH) instabilities and wind mixing, drastically increases the dust yields. △ Less

Submitted 17 October, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Comments: 13 pages, 18 figures, 10 tables

How D-type HII region expansion depends on numerical resolution

Authors: J. M. Pittard, M. M. Kupilas, C. J. Wareing

Abstract: 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… ▽ More 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. △ Less

Submitted 16 December, 2021; originally announced December 2021.

Comments: 6 pages, 2 figures, accepted for publication in MNRAS

How to inflate a wind-blown bubble

Authors: J. M. Pittard, C. J. Wareing, M. M. Kupilas

Abstract: Stellar winds are one of several ways that massive stars can affect the star formation process on local and galactic scales. In this paper we investigate the numerical resolution needed to inflate an energy-driven stellar wind bubble in an external medium. We find that the radius of the wind injection region, $r_{\rm inj}$, must be below a maximum value, $r_{\rm inj,max}$, in order for a bubble to… ▽ More Stellar winds are one of several ways that massive stars can affect the star formation process on local and galactic scales. In this paper we investigate the numerical resolution needed to inflate an energy-driven stellar wind bubble in an external medium. We find that the radius of the wind injection region, $r_{\rm inj}$, must be below a maximum value, $r_{\rm inj,max}$, in order for a bubble to be produced, but must be significantly below this value if the bubble properties are to closely agree with analytical predictions. The final bubble momentum is within 25 per cent of the value from a higher resolution reference model if $χかい= r_{\rm inj}/r_{\rm inj,max}$ = 0.1. Our work has significance for the amount of radial momentum that a wind-blown bubble can impart to the ambient medium in simulations, and thus on the relative importance of stellar wind feedback. △ Less

Submitted 17 September, 2021; v1 submitted 30 July, 2021; originally announced July 2021.

Comments: 9 pages, 6 figures. Accepted by MNRAS

Particle acceleration and non-thermal emission in colliding-wind binary systems

Authors: J. M. Pittard, G. E. Romero, G. S. Vila

Abstract: We present a model for the creation of non-thermal particles via diffusive shock acceleration in a colliding-wind binary. Our model accounts for the oblique nature of the global shocks bounding the wind-wind collision region and the finite velocity of the scattering centres to the gas. It also includes magnetic field amplification by the cosmic ray induced streaming instability and the dynamical b… ▽ More We present a model for the creation of non-thermal particles via diffusive shock acceleration in a colliding-wind binary. Our model accounts for the oblique nature of the global shocks bounding the wind-wind collision region and the finite velocity of the scattering centres to the gas. It also includes magnetic field amplification by the cosmic ray induced streaming instability and the dynamical back reaction of the amplified field. We assume that the injection of the ions and electrons is independent of the shock obliquity and that the scattering centres move relative to the fluid at the Alfvén velocity (resulting in steeper non-thermal particle distributions). We find that the Mach number, Alfvénic Mach number, and transverse field strength vary strongly along and between the shocks, resulting in significant and non-linear variations in the particle acceleration efficiency and shock nature (turbulent vs. non-turbulent). We find much reduced compression ratios at the oblique shocks in most of our models compared to our earlier work, though total gas compression ratios that exceed 20 can still be obtained in certain situations. We also investigate the dependence of the non-thermal emission on the stellar separation and determine when emission from secondary electrons becomes important. We finish by applying our model to WR 146, one of the brightest colliding wind binaries in the radio band. We are able to match the observed radio emission and find that roughly 30 per cent of the wind power at the shocks is channelled into non-thermal particles. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 22 pages, 19 figures, accepted for publication in MNRAS

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

Authors: M. M. Kupilas, C. J. Wareing, J. M. Pittard, S. A. E. G. Falle

Abstract: 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 in… ▽ More 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. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Striations, integrals, hourglasses and collapse -- thermal instability driven magnetic simulations of molecular clouds

Authors: C. J. Wareing, J. M. Pittard, S. A. E. G. Falle

Abstract: The MHD version of the adaptive mesh refinement (AMR) code, MG, has been employed to study the interaction of thermal instability, magnetic fields and gravity through 3D simulations of the formation of collapsing cold clumps on the scale of a few parsecs, inside a larger molecular cloud. The diffuse atomic initial condition consists of a stationary, thermally unstable, spherical cloud in pressure… ▽ More The MHD version of the adaptive mesh refinement (AMR) code, MG, has been employed to study the interaction of thermal instability, magnetic fields and gravity through 3D simulations of the formation of collapsing cold clumps on the scale of a few parsecs, inside a larger molecular cloud. The diffuse atomic initial condition consists of a stationary, thermally unstable, spherical cloud in pressure equilibrium with lower density surroundings and threaded by a uniform magnetic field. This cloud was seeded with 10% density perturbations at the finest initial grid level around n=1.1 cm^{-3} and evolved with self-gravity included from the outset. Several cloud diameters were considered (100 pc, 200 pc and 400 pc) equating to several cloud masses (17,000 Msun, 136,000 Msun and 1.1x10^6 Msun). Low-density magnetic-field-aligned striations were observed as the clouds collapse along the field lines into disc-like structures. The induced flow along field lines leads to oscillations of the sheet about the gravitational minimum and an integral-shaped appearance. When magnetically supercritical, the clouds then collapse and generate hourglass magnetic field configurations with strongly intensified magnetic fields, reproducing observational behaviour. Resimulation of a region of the highest mass cloud at higher resolution forms gravitationally-bound collapsing clumps within the sheet that contain clump-frame supersonic (M~5) and super-Alfvenic (M_A~4) velocities. Observationally realistic density and velocity power spectra of the cloud and densest clump are obtained. Future work will use these realistic initial conditions to study individual star and cluster feedback. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: Accepted for publication in MNRAS. 20 pages, 11 figures + appendix and figure. Raw data: doi.org/10.5518/897

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

Authors: E. Mossoux, J. M. Pittard, G. Rauw, Y. Nazé

Abstract: 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… ▽ More 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. △ Less

Submitted 20 March, 2020; originally announced March 2020.

Comments: 12 pages, 8 figures, A&A accepted

Journal ref: A&A 636, A109 (2020)

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

Authors: J. C. Morford, D. M. Fenech, R. K. Prinja, R. Blomme, J. A. Yates, J. J. Drake, S. P. S. Eyres, A. M. S. Richards, I. R. Stevens, N. J. Wright, J. S. Clark, S. Dougherty, J. M. Pittard, H. Smith, J. S. Vink

Abstract: 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… ▽ More 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. △ Less

Submitted 20 March, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 41 pages, 12 figures, accepted in A&A

Journal ref: A&A 637, A64 (2020)

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

Authors: J. M. Pittard, G. S. Vila, G. E. Romero

Abstract: 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 parti… ▽ More 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. △ Less

Submitted 15 April, 2020; v1 submitted 11 December, 2019; originally announced December 2019.

Comments: 18 pages, 24 figures, accepted by MNRAS

The isothermal evolution of a shock-filament interaction

Authors: K. J. A. Goldsmith, J. M. Pittard

Abstract: 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 an… ▽ More 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. △ Less

Submitted 11 December, 2019; originally announced December 2019.

Comments: 21 pages, 26 figures

Momentum and energy injection by a supernova remnant into an inhomogeneous medium

Authors: Julian M. Pittard

Abstract: We investigate the effect of mass-loading from embedded clouds on the evolution of supernova remnants and on the energy and momentum that they inject into an inhomogeneous interstellar medium. We use 1D hydrodynamical calculations and assume that the clouds are numerous enough that they can be treated in the continuous limit. The destruction of embedded clouds adds mass into the remnant, increasin… ▽ More We investigate the effect of mass-loading from embedded clouds on the evolution of supernova remnants and on the energy and momentum that they inject into an inhomogeneous interstellar medium. We use 1D hydrodynamical calculations and assume that the clouds are numerous enough that they can be treated in the continuous limit. The destruction of embedded clouds adds mass into the remnant, increasing its density and pressure, and decreasing its temperature. The remnant cools more quickly, is less able to do PdV work on the swept-up gas, and ultimately attains a lower final momentum (by up to a factor of two or more). We thus find that the injection of momentum is more sensitive to an inhomogeneous environment than previous work has suggested, and we provide fits to our results for the situation where the cloud mass is not limited. The behaviour of the remnant is more complex in situations where the cloud mass is finite and locally runs out. In the case of multiple supernovae in a clustered environment, later supernova explosions may encounter higher densities than previous explosions due to the prior liberation of mass from engulfed clouds. If the cloud mass is finite, later explosions may be able to create a sustained hot phase when earlier explosions have not been able to. △ Less

Submitted 8 July, 2019; originally announced July 2019.

Comments: 21 pages, 16 figures, accepted for publication in MNRAS

Non-thermal X-rays from Colliding Wind Shock Acceleration in the Massive Binary Eta Carinae

Authors: Kenji Hamaguchi, Michael F. Corcoran, Julian M. Pittard, Neetika Sharma, Hiromitsu Takahashi, Christopher M. P. Russell, Brian W. Grefenstette, Daniel R. Wik, Theodore R. Gull, Noel D. Richardson, Thomas I. Madura, Anthony F. J. Moffat

Abstract: Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study, we still do not have a complete census of acceleration mechanisms. The collision of strong stellar winds in massive binary systems creates powerful shocks, which have been expected to produce high-energy cosmic-rays through Fermi acceleration at the shock interface. The accelerated particles should c… ▽ More Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study, we still do not have a complete census of acceleration mechanisms. The collision of strong stellar winds in massive binary systems creates powerful shocks, which have been expected to produce high-energy cosmic-rays through Fermi acceleration at the shock interface. The accelerated particles should collide with stellar photons or ambient material, producing non-thermal emission observable in X-rays and gamma-rays. The supermassive binary star eta Carinae drives the strongest colliding wind shock in the solar neighborhood. Observations with non-focusing high-energy observatories indicate a high energy source near eta Carinae, but have been unable to conclusively identify eta Carinae as the source because of their relatively poor angular resolution. Here we present the first direct focussing observations of the non-thermal source in the extremely hard X-ray band, which is found to be spatially coincident with the star within several arc-seconds. These observations show that the source of non-thermal X-rays varies with the orbital phase of the binary, and that the photon index of the emission is similar to that derived through analysis of the gamma-ray spectrum. This is conclusive evidence that the high-energy emission indeed originates from non-thermal particles accelerated at colliding wind shocks. △ Less

Submitted 19 April, 2019; originally announced April 2019.

Comments: 36 pages, 7 figures, 2 tables

Journal ref: Nature Astronomy 2 (2018) 731-736

Sheets, filaments and clumps - high resolution simulations of how the thermal instability can form molecular clouds

Authors: C. J. Wareing, S. A. E. G. Falle, J. M. Pittard

Abstract: This paper describes 3D simulations of the formation of collapsing cold clumps via thermal instability inside a larger cloud complex. The initial condition was a diffuse atomic, stationary, thermally unstable, 200pc diameter spherical cloud in pressure equilibrium with low density surroundings. This was seeded with 10% density perturbations at the finest initial grid level (0.29pc) around n_H = 1.… ▽ More This paper describes 3D simulations of the formation of collapsing cold clumps via thermal instability inside a larger cloud complex. The initial condition was a diffuse atomic, stationary, thermally unstable, 200pc diameter spherical cloud in pressure equilibrium with low density surroundings. This was seeded with 10% density perturbations at the finest initial grid level (0.29pc) around n_H = 1.1cm^{-3} and evolved with self-gravity included. No magnetic field was imposed. Resimulations at a higher resolution of a region extracted from this simulation (down to 0.039pc), show that the thermal instability forms sheets, then filaments and finally clumps. The width of the filaments increases over time, in one particular case from 0.26 to 0.56pc. Thereafter clumps with sizes of around 5pc grow at the intersections of filaments. 21 distinct clumps, with properties similar to those observed in molecular clouds, are found by using the FellWalker algorithm to find minima in the gravitational potential. Not all of these are gravitationally bound, but the convergent nature of the flow and increasing central density suggest they are likely to form stars. Further simulation of the most massive clump shows the gravitational collapse to a density >10^6 cm^{-3}. These results provide realistic initial conditions that can be used to study feedback in individual clumps, interacting clumps and the entire molecular cloud complex. △ Less

Submitted 12 March, 2019; v1 submitted 21 December, 2018; originally announced December 2018.

Comments: 19 pages, 10 figures. MNRAS accepted. Moderate revision and extension

A comparison of shock-cloud and wind-cloud interactions: effect of increased cloud density contrast on cloud evolution

Authors: K. J. A. Goldsmith, J. M. Pittard

Abstract: The similarities, or otherwise, of a shock or wind interacting with a cloud of density contras t$ χかい= 10$ were explored in a previous paper. Here, we investigate such interactions with clouds of higher density contrast. We compare the adiabatic hydrodynamic interaction of a Mach 10 shock with a spherical cloud of $χかい= 10^{3}$ with that of a cloud embedded in a wind with identical parameters to the p… ▽ More The similarities, or otherwise, of a shock or wind interacting with a cloud of density contras t$ χかい= 10$ were explored in a previous paper. Here, we investigate such interactions with clouds of higher density contrast. We compare the adiabatic hydrodynamic interaction of a Mach 10 shock with a spherical cloud of $χかい= 10^{3}$ with that of a cloud embedded in a wind with identical parameters to the post-shock flow. We find that initially there are only minor morphological differences between the shock-cloud and wind-cloud interactions, compared to when $χかい= 10$. However, once the transmitted shock exits the cloud, the development of a turbulent wake and fragmentation of the cloud differs between the two simulations. On increasing the wind Mach number, we note the development of a thin, smooth tail of cloud material, which is then disrupted by the fragmentation of the cloud core and subsequent `mass-loading' of the flow. We find that the normalized cloud mixing time ($t_{mix}$) is shorter at higher $χかい$. However, a strong Mach number dependence on tmix and the normalized cloud drag time, $t'_{drag}$, is not observed. Mach-number-dependent values of $t_{mix}$ and $t'_{drag}$ from comparable shock-cloud interactions converge towards the Mach-number-independent time-scales of the wind-cloud simulations. We find that high $χかい$ clouds can be accelerated up to 80-90 per cent of the wind velocity and travel large distances before being significantly mixed. However, complete mixing is not achieved in our simulations and at late times the flow remains perturbed. △ Less

Submitted 12 May, 2018; originally announced May 2018.

Comments: 11 pages, 5 figures

Journal ref: MNRAS 476 (2018) 2209

Colliding Stellar Winds Structure and X-ray Emission

Authors: J. M. Pittard, B. Dawson

Abstract: We investigate the structure and X-ray emission from the colliding stellar winds in massive star binaries. We find that the opening angle of the contact discontinuity (CD) is overestimated by several formulae in the literature at very small values of the wind momentum ratio. We find also that the shocks in the primary (dominant) and secondary winds flare by approx 20 degrees compared to the CD, an… ▽ More We investigate the structure and X-ray emission from the colliding stellar winds in massive star binaries. We find that the opening angle of the contact discontinuity (CD) is overestimated by several formulae in the literature at very small values of the wind momentum ratio. We find also that the shocks in the primary (dominant) and secondary winds flare by approx 20 degrees compared to the CD, and that the entire secondary wind is shocked when the wind momentum ratio < 0.02. Analytical expressions for the opening angles of the shocks, and the fraction of each wind that is shocked, are provided. We find that the X-ray luminosity scales with the wind momentum ratio, and that the spectrum softens slightly as the wind momentum ratio decreases. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: 6 pages, 3 figures, accepted by MNRAS

A new mechanical stellar wind feedback model for the Rosette Nebula

Authors: C. J. Wareing, J. M. Pittard, N. J. Wright, S. A. E. G. Falle

Abstract: The famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2-6 million-year-old co-distant and co-moving central star cluster NGC 2244. However, with upper age estimates of less than 110,000 years, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data r… ▽ More The famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2-6 million-year-old co-distant and co-moving central star cluster NGC 2244. However, with upper age estimates of less than 110,000 years, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data reveals the ejection of the most massive star in the Rosette, HD46223, from NGC 2244 occurred 1.73 (+0.34,-0.25)Myr (1$σしぐま$ uncertainty) in the past. Assuming this ejection was at the birth of the most massive stars in NGC 2244, including the dominant centrally positioned HD46150, the age is set for the famous ionised region at more than ten times that derived for the cavity. Here, we are able to reproduce the structure of the Rosette Nebula, through simulation of mechanical stellar feedback from a 40M$_{\odot}$ star in a thin sheet-like molecular cloud. We form the 135,000M$_{\odot}$ cloud from thermally-unstable diffuse interstellar medium under the influence of a realistic background magnetic field with thermal/magnetic pressure equilibrium. Properties derived from a snapshot of the simulation at 1.5Myr, including cavity size, stellar age, magnetic field and resulting inclination to the line of sight, match those derived from observations. An elegant explanation is thus provided for the stark contrast in age estimates based on realistic diffuse ISM properties, molecular cloud formation and stellar wind feedback. △ Less

Submitted 13 February, 2018; originally announced February 2018.

Comments: This is the Version of Record, published RAS Open Access CC BY online in MNRAS, 2018 Feb 13. 15 Pages, 9 figures

Journal ref: MNRAS, Volume 475, Issue 3, Pages 3598-3612, 2018

A changing wind collision

Authors: Yael Naze, Gloria Koenigsberger, Julian M. Pittard, Elliot Ross Parkin, Gregor Rauw, Michael F. Corcoran, D. John Hillier

Abstract: We report on the first detection of a global change in the X-ray emitting properties of a wind-wind collision, thanks to XMM-Newton observations of the massive SMC system HD5980. While its lightcurve had remained unchanged between 2000 and 2005, the X-ray flux has now increased by a factor of ~2.5, and slightly hardened. The new observations also extend the observational coverage over the entire o… ▽ More We report on the first detection of a global change in the X-ray emitting properties of a wind-wind collision, thanks to XMM-Newton observations of the massive SMC system HD5980. While its lightcurve had remained unchanged between 2000 and 2005, the X-ray flux has now increased by a factor of ~2.5, and slightly hardened. The new observations also extend the observational coverage over the entire orbit, pinpointing the lightcurve shape. It has not varied much despite the large overall brightening, and a tight correlation of fluxes with orbital separation is found, without any hysteresis effect. Moreover, the absence of eclipses and of absorption effects related to orientation suggests a large size for the X-ray emitting region. Simple analytical models of the wind-wind collision, considering the varying wind properties of the eruptive component in HD5980, are able to reproduce the recent hardening and the flux-separation relationship, at least qualitatively, but they predict a hardening at apastron and little change in mean flux, contrary to observations. The brightness change could then possibly be related to a recently theorized phenomenon linked to the varying strength of thin-shell instabilities in shocked wind regions. △ Less

Submitted 15 December, 2017; originally announced December 2017.

Comments: accepted by ApJ

Journal ref: ApJ, 853, 164 (11p, 2018)

Variable millimetre radiation from the colliding-wind binary Cyg OB2 #8A

Authors: R. Blomme, D. M. Fenech, R. K. Prinja, J. M. Pittard, J. C. Morford

Abstract: In the colliding-wind region of massive binaries, non-thermal radio emission occurs. This non-thermal radio emission (due to synchrotron radiation) has so far been observed at centimetre wavelengths. At millimetre wavelengths, the stellar winds and the colliding-wind region emit more thermal free-free radiation, and it is expected that any non-thermal contribution will be difficult or impossible t… ▽ More In the colliding-wind region of massive binaries, non-thermal radio emission occurs. This non-thermal radio emission (due to synchrotron radiation) has so far been observed at centimetre wavelengths. At millimetre wavelengths, the stellar winds and the colliding-wind region emit more thermal free-free radiation, and it is expected that any non-thermal contribution will be difficult or impossible to detect. We aim to determine if the material in the colliding-wind region contributes substantially to the observed millimetre fluxes of a colliding-wind binary. We also try to distinguish the synchrotron emission from the free-free emission. We monitored the massive binary Cyg OB2 #8A at 3 mm with the NOrthern Extended Millimeter Array (NOEMA) interferometer of the Institut de Radioastronomie Millimetrique (IRAM). The data were collected in 14 separate observing runs (in 2014 and 2016), and provide good coverage of the orbital period. The observed millimetre fluxes range between 1.1 and 2.3 mJy, and show phase-locked variability, clearly indicating that a large part of the emission is due to the colliding-wind region. A simple synchrotron model gives fluxes with the correct order of magnitude, but with a maximum that is phase-shifted with respect to the observations. Qualitatively this phase shift can be explained by our neglect of orbital motion on the shape of the colliding-wind region. A model using only free-free emission results in only a slightly worse explanation of the observations. Additionally, on the map of our observations we also detect the O6.5 III star Cyg OB2 #8B, for which we determine a 3 mm flux of 0.21 +- 0.033 mJy. The question of whether synchrotron radiation or free-free emission dominates the millimetre fluxes of Cyg OB2 #8A remains open. More detailed modelling of this system, based on solving the hydrodynamical equations, is required to give a definite answer. △ Less

Submitted 16 August, 2017; originally announced August 2017.

Comments: 9 pages, 4 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 608, A69 (2017)

A comparison of shock-cloud and wind-cloud interactions: The longer survival of clouds in winds

Authors: K. J. A. Goldsmith, J. M. Pittard

Abstract: The interaction of a hot, high-velocity wind with a cold, dense molecular cloud has often been assumed to resemble the evolution of a cloud embedded in a post-shock flow. However, no direct comparative study of these two processes currently exists in the literature. We present 2D adiabatic hydrodynamical simulations of the interaction of a Mach 10 shock with a cloud of density contrast $χかい= 10$ and… ▽ More The interaction of a hot, high-velocity wind with a cold, dense molecular cloud has often been assumed to resemble the evolution of a cloud embedded in a post-shock flow. However, no direct comparative study of these two processes currently exists in the literature. We present 2D adiabatic hydrodynamical simulations of the interaction of a Mach 10 shock with a cloud of density contrast $χかい= 10$ and compare our results with those of a commensurate wind-cloud simulation. We then investigate the effect of varying the wind velocity, effectively altering the wind Mach number $M_{wind}$, on the cloud's evolution. We find that there are significant differences between the two processes: 1) the transmitted shock is much flatter in the shock-cloud interaction; 2) a low-pressure region in the wind-cloud case deflects the flow around the edge of the cloud in a different manner to the shock-cloud case; 3) there is far more axial compression of the cloud in the case of the shock. As $M_{wind}$ increases, the normalised rate of mixing is reduced. Clouds in winds with higher $M_{wind}$ also do not experience a transmitted shock through the cloud's rear and are more compressed axially. In contrast with shock-cloud simulations, the cloud mixing time normalised by the cloud-crushing time-scale $t_{cc}$ increases for increasing $M_{wind}$ until it plateaus (at $t_{mix} \simeq 25 \, t_{cc}$) at high $M_{wind}$, thus demonstrating the expected Mach scaling. In addition, clouds in high Mach number winds are able to survive for long durations and are capable of being moved considerable distances. △ Less

Submitted 12 June, 2017; originally announced June 2017.

Comments: 13 pages, 5 figures

Journal ref: MNRAS 470 (2017) 2427

Hydrodynamic simulations of mechanical stellar feedback in a molecular cloud formed by thermal instability

Authors: Christopher J. Wareing, Julian M. Pittard, Samuel A. E. G. Falle

Abstract: We have used the AMR hydrodynamic code, MG, to perform 3D hydrodynamic simulations with self-gravity of stellar feedback in a spherical clumpy molecular cloud formed through the action of thermal instability. We simulate the interaction of the mechanical energy input from 15 Msun, 40 Msun, 60 Msun and 120 Msun stars into a 100 pc-diameter 16,500 Msun cloud with a roughly spherical morphology with… ▽ More We have used the AMR hydrodynamic code, MG, to perform 3D hydrodynamic simulations with self-gravity of stellar feedback in a spherical clumpy molecular cloud formed through the action of thermal instability. We simulate the interaction of the mechanical energy input from 15 Msun, 40 Msun, 60 Msun and 120 Msun stars into a 100 pc-diameter 16,500 Msun cloud with a roughly spherical morphology with randomly distributed high density condensations. The stellar winds are introduced using appropriate non-rotating Geneva stellar evolution models. In the 15 Msun star case, the wind has very little effect, spreading around a few neighbouring clumps before becoming overwhelmed by the cloud collapse. In contrast, in the 40 Msun, 60 Msun and 120 Msun star cases, the more powerful stellar winds create large cavities and carve channels through the cloud, breaking out into the surrounding tenuous medium during the wind phase and considerably altering the cloud structure. After 4.97 Myrs, 3.97 Myrs and 3.01 Myrs respectively, the massive stars explode as supernovae (SNe). The wind-sculpted surroundings considerably affect the evolution of these SN events as they both escape the cloud along wind-carved channels and sweep up remaining clumps of cloud/wind material. The `cloud' as a coherent structure does not survive the SN from any of these stars, but only in the 120 Msun case is the cold molecular material completely destabilised and returned to the unstable thermal phase. In the 40 Msun and 60 Msun cases, coherent clumps of cold material are ejected from the cloud by the SN, potentially capable of further star formation. △ Less

Submitted 6 June, 2017; originally announced June 2017.

Comments: 33 pages, 26 figures, accepted to MNRAS 2017 June 6. Raw data available from http://doi.org/10.5518/201

Hydrodynamic ablation of protoplanetary disks via supernovae

Authors: J. L. Close, J. M. Pittard

Abstract: We present three-dimensional simulations of a protoplanetary disk subject to the effect of a nearby (0.3pc distant) supernova, using a time-dependent flow from a one dimensional numerical model of the supernova remnant (SNR), in addition to constant peak ram pressure simulations. Simulations are performed for a variety of disk masses and inclination angles. We find disk mass-loss rates that are ty… ▽ More We present three-dimensional simulations of a protoplanetary disk subject to the effect of a nearby (0.3pc distant) supernova, using a time-dependent flow from a one dimensional numerical model of the supernova remnant (SNR), in addition to constant peak ram pressure simulations. Simulations are performed for a variety of disk masses and inclination angles. We find disk mass-loss rates that are typically 1e-7 to 1e-6 Msol/yr (but peak near 1e-5 Msol/yr during the "instantaneous" stripping phase) and are sustained for around 200 yr. Inclination angle has little effect on the mass loss unless the disk is close to edge-on. Inclined disks also strip asymmetrically with the trailing edge ablating more easily. Since the interaction lasts less than one outer rotation period, there is not enough time for the disk to restore its symmetry, leaving the disk asymmetrical after the flow has passed. Of the low-mass disks considered, only the edge-on disk is able to survive interaction with the SNR (with 50% of its initial mass remaining). At the end of the simulations, disks that survive contain fractional masses of SN material up to 5e-6. This is too low to explain the abundance of short-lived radionuclides in the early solar system, but a larger disk and the inclusion of radiative cooling might allow the disk to capture a higher fraction of SN material. △ Less

Submitted 20 April, 2017; originally announced April 2017.

Comments: 16 pages, 16 figures, 2 tables, accepted by MNRAS

The interaction of a magnetohydrodynamical shock with a filament

Authors: K. J. A. Goldsmith, J. M. Pittard

Abstract: 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 ef… ▽ More 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. △ Less

Submitted 13 July, 2016; originally announced July 2016.

Comments: 30 pages, 42 figures, accepted by MNRAS

Journal ref: MNRAS 461 (2016) 578

Magnetohydrodynamic simulations of mechanical stellar feedback in a sheet-like molecular cloud

Authors: C. J. Wareing, J. M. Pittard, S. A. E. G. Falle

Abstract: We have used the AMR hydrodynamic code, MG, to perform 3D magnetohydrodynamic simulations with self-gravity of stellar feedback in a sheet-like molecular cloud formed through the action of the thermal instability. We simulate the interaction of the mechanical energy input from a 15 solar mass star and a 40 solar mass star into a 100 pc-diameter 17000 solar mass cloud with a corrugated sheet morpho… ▽ More We have used the AMR hydrodynamic code, MG, to perform 3D magnetohydrodynamic simulations with self-gravity of stellar feedback in a sheet-like molecular cloud formed through the action of the thermal instability. We simulate the interaction of the mechanical energy input from a 15 solar mass star and a 40 solar mass star into a 100 pc-diameter 17000 solar mass cloud with a corrugated sheet morphology that in projection appears filamentary. The stellar winds are introduced using appropriate Geneva stellar evolution models. In the 15 solar mass star case, the wind forms a narrow bipolar cavity with minimal effect on the parent cloud. In the 40 solar mass star case, the more powerful stellar wind creates a large cylindrical cavity through the centre of the cloud. After 12.5 Myrs and 4.97 Myrs respectively, the massive stars explode as supernovae (SNe). In the 15 solar mass star case, the SN material and energy is primarily deposited into the molecular cloud surroundings over ~10^5 years before the SN remnant escapes the cloud. In the 40 solar mass star case, a significant fraction of the SN material and energy rapidly escapes the molecular cloud along the wind cavity in a few tens of kiloyears. Both SN events compress the molecular cloud material around them to higher densities (so may trigger further star formation), and strengthen the magnetic field, typically by factors of 2-3 but up to a factor of 10. Our simulations are relevant to observations of bubbles in flattened ring-like molecular clouds and bipolar HII regions. △ Less

Submitted 16 November, 2016; v1 submitted 16 May, 2016; originally announced May 2016.

Comments: 29 pages, 20 figures, MNRAS accepted. Additions after review include an extended analysis in a new Section 4, clarification of numerous minor points throughout and quantitative data figures in a new Appendix

Eta Carinae's Thermal X-ray Tail Measured with XMM-Newton and NuSTAR

Authors: Kenji Hamaguchi, Michael F. Corcoran, Theodore R. Gull, Hiromitsu Takahashi, Brian W. Grefenstette, Takayuki Yuasa, Martin Stuhlinger, Christopher M. P. Russell, Anthony F. J. Moffat, Neetika Sharma, Thomas I. Madura, Noel D. Richardson, Jose Groh, Julian M. Pittard, Stanley Owocki

Abstract: The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 k… ▽ More The evolved, massive highly eccentric binary system, eta Carinae, underwent a periastron passage in the summer of 2014. We obtained two coordinated X-ray observations with XMM-Newton and NuSTAR during the elevated X-ray flux state and just before the X-ray minimum flux state around this passage. These NuSTAR observations clearly detected X-ray emission associated with eta Car extending up to ~50 keV for the first time. The NuSTAR spectrum above 10 keV can be fit with the bremsstrahlung tail from a kT ~6 keV plasma. This temperature is Delta kT ~2 keV higher than those measured from the iron K emission line complex, if the shocked gas is in collisional ionization equilibrium. This result may suggest that the companion star's pre-shock wind velocity is underestimated. The NuSTAR observation near the X-ray minimum state showed a gradual decline in the X-ray emission by 40% at energies above 5 keV in a day, the largest rate of change of the X-ray flux yet observed in individual eta Car observations. The column density to the hardest emission component, NH ~1e24 cm-2, marked one of the highest values ever observed for eta Car, strongly suggesting increased obscuration of the wind-wind colliding X-ray emission by the thick primary stellar wind prior to superior conjunction. Neither observation detected the power-law component in the extremely hard band that INTEGRAL and Suzaku observed prior to 2011. If the non-detection by NuSTAR is caused by absorption, the power-law source must be small and located very near the WWC apex. Alternatively, it may be that the power-law source is not related to either eta Car or the GeV gamma-ray source. △ Less

Submitted 2 February, 2016; originally announced February 2016.

Comments: 24 pages, 6 figures, 2 tables

Journal ref: The Astrophysical Journal, 817:23, 2016 January 20

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

Authors: C. J. Wareing, J. M. Pittard, S. A. E. G. Falle, S. Van Loo

Abstract: 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 o… ▽ More 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. △ Less

Submitted 8 March, 2016; v1 submitted 17 January, 2016; originally announced January 2016.

Comments: 17 pages, 10 figures. Accepted to MNRAS, minor changes and additions to text, primarily in Section 6

Numerical Simulations of a Shock-Filament Interaction

Authors: J. M. Pittard, K. J. A. Goldsmith

Abstract: 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 differenc… ▽ More 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. △ Less

Submitted 11 February, 2016; v1 submitted 18 December, 2015; originally announced December 2015.

Comments: 29 pages, 46 figures, accepted by MNRAS

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

Authors: J. M. Pittard, E. R. Parkin

Abstract: 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 th… ▽ More 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. △ Less

Submitted 6 January, 2016; v1 submitted 19 October, 2015; originally announced October 2015.

Comments: 33 pages, 40 figures, accepted to MNRAS

Suzaku Monitoring of Hard X-ray Emission from Eta Carinae over a Single Binary Orbital Cycle

Authors: Kenji Hamaguchi, Michael F. Corcoran, Hiromitsu Takahashi, Takayuki Yuasa, Manabu Ishida, Theodore R. Gull, Julian M. Pittard, Christopher M. P. Russell, Thomas I. Madura

Abstract: The Suzaku X-ray observatory monitored the supermassive binary system Eta Carinae 10 times during the whole 5.5 year orbital cycle between 2005-2011. This series of observations presents the first long-term monitoring of this enigmatic system in the extremely hard X-ray band between 15-40 keV. During most of the orbit, the 15-25 keV emission varied similarly to the 2-10 keV emission, indicating an… ▽ More The Suzaku X-ray observatory monitored the supermassive binary system Eta Carinae 10 times during the whole 5.5 year orbital cycle between 2005-2011. This series of observations presents the first long-term monitoring of this enigmatic system in the extremely hard X-ray band between 15-40 keV. During most of the orbit, the 15-25 keV emission varied similarly to the 2-10 keV emission, indicating an origin in the hard energy tail of the kT ~4 keV wind-wind collision (WWC) plasma. However, the 15-25 keV emission declined only by a factor of 3 around periastron when the 2-10 keV emission dropped by two orders of magnitude due probably to an eclipse of the WWC plasma. The observed minimum in the 15-25 keV emission occurred after the 2-10 keV flux had already recovered by a factor of ~3. This may mean that the WWC activity was strong, but hidden behind the thick primary stellar wind during the eclipse. The 25-40 keV flux was rather constant through the orbital cycle, at the level measured with INTEGRAL in 2004. This result may suggest a connection of this flux component to the gamma-ray source detected in this field. The Helium-like Fe Kalpha line complex at ~6.7 keV became strongly distorted toward periastron as seen in the previous cycle. The 5-9 keV spectra can be reproduced well with a two-component spectral model, which includes plasma in collision equilibrium (CE) and a plasma in non-equilibrium ionization (NEI) with tau ~1e11 cm-3 s-1. The NEI plasma increases in importance toward periastron. △ Less

Submitted 22 October, 2014; originally announced October 2014.

Comments: 19 pages, 6 figures

Journal ref: Astrophysical Journal, 2014, 795, 119

V444 Cyg X-ray and polarimetric variability: Radiative and Coriolis forces shape the wind collision region

Authors: Jamie R. Lomax, Yael Naze, Jennifer L. Hoffman, Christopher M. P. Russell, Michael De Becker, Michael F. Corcoran, James W. Davidson, Hilding R. Neilson, Stan Owocki, Julian M. Pittard, Andy M. T. Pollock

Abstract: We present results from a study of the eclipsing, colliding-wind binary V444 Cyg that uses a combination of X-ray and optical spectropolarimetric methods to describe the 3-D nature of the shock and wind structure within the system. We have created the most complete X-ray light curve of V444 Cyg to date using 40 ksec of new data from Swift, and 200 ksec of new and archived XMM-Newton observations.… ▽ More We present results from a study of the eclipsing, colliding-wind binary V444 Cyg that uses a combination of X-ray and optical spectropolarimetric methods to describe the 3-D nature of the shock and wind structure within the system. We have created the most complete X-ray light curve of V444 Cyg to date using 40 ksec of new data from Swift, and 200 ksec of new and archived XMM-Newton observations. In addition, we have characterized the intrinsic, polarimetric phase-dependent behavior of the strongest optical emission lines using data obtained with the University of Wisconsin's Half-Wave Spectropolarimeter. We have detected evidence of the Coriolis distortion of the wind-wind collision in the X-ray regime, which manifests itself through asymmetric behavior around the eclipses in the system's X-ray light curves. The large opening angle of the X-ray emitting region, as well as its location (i.e. the WN wind does not collide with the O star, but rather its wind) are evidence of radiative braking/inhibition occurring within the system. Additionally, the polarimetric results show evidence of the cavity the wind-wind collision region carves out of the Wolf-Rayet star's wind. △ Less

Submitted 22 October, 2014; originally announced October 2014.

Comments: 19 pages, 14 figures, accepted A&A

Journal ref: A&A 573, A43 (2015)

Numerical Simulations of a Shock Interacting with Multiple Magnetized Clouds

Authors: R. Alūzas, J. M. Pittard, S. A. E. G. Falle, T. W. Hartquist

Abstract: 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 o… ▽ More 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. △ Less

Submitted 24 July, 2014; originally announced July 2014.

Comments: 24 pages, 30 figures. Accepted for publication in MNRAS

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

Authors: E. R. Parkin, J. M. Pittard, Y. Nazé, R. Blomme

Abstract: 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 motio… ▽ More 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) △ Less

Submitted 22 June, 2014; originally announced June 2014.

Comments: 20 pages, 18 figures. Accepted for publication in A&A

Journal ref: A&A 570, A10 (2014)

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

Authors: H. Rogers, J. M. Pittard

Abstract: 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 f… ▽ More 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. △ Less

Submitted 20 May, 2014; originally announced May 2014.

X-ray Emission from Eta Carinae near Periastron in 2009 I: A Two State Solution

Authors: Kenji Hamaguchi, Michael F. Corcoran, Christopher Russell, Andrew M. T. Pollock, Theodore R. Gull, Mairan Teodoro, Thomas I. Madura, Augusto Damineli, Julian M. Pittard

Abstract: X-ray emission from the supermassive binary system Eta Carinae declines sharply around periastron. This X-ray minimum has two distinct phases - the lowest flux phase in the first ~3 weeks and a brighter phase thereafter. In 2009, the Chandra X-ray Observatory monitored the first phase five times and found the lowest observed flux at ~1.9e-12 ergs cm-2 s-1 (3-8 keV). The spectral shape changed such… ▽ More X-ray emission from the supermassive binary system Eta Carinae declines sharply around periastron. This X-ray minimum has two distinct phases - the lowest flux phase in the first ~3 weeks and a brighter phase thereafter. In 2009, the Chandra X-ray Observatory monitored the first phase five times and found the lowest observed flux at ~1.9e-12 ergs cm-2 s-1 (3-8 keV). The spectral shape changed such that the hard band above ~4 keV dropped quickly at the beginning and the soft band flux gradually decreased to its lowest observed value in ~2 weeks. The hard band spectrum had begun to recover by that time. This spectral variation suggests that the shocked gas producing the hottest X-ray gas near the apex of the wind-wind collision (WWC) is blocked behind the dense inner wind of the primary star, which later occults slightly cooler gas downstream. Shocked gas previously produced by the system at earlier orbital phases is suggested to produce the faint residual X-ray emission seen when the emission near the apex is completely blocked by the primary wind. The brighter phase is probably caused by the re-appearance of the WWC plasma, whose emissivity significantly declined during the occultation. We interpret this to mean that the X-ray minimum is produced by a hybrid mechanism of an occultation and a decline in emissivity of the WWC shock. △ Less

Submitted 22 January, 2014; originally announced January 2014.

Comments: 10 pages, 3 figures, accepted for publication in ApJ

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

Authors: J. L. Close, J. M. Pittard, T. W. Hartquist, S. A. E. G. Falle

Abstract: 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 stri… ▽ More 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 △ Less

Submitted 20 September, 2013; originally announced September 2013.

Comments: 10 pages, 8 figures, 3 tables; accepted for publication in MNRAS

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

Authors: H. Rogers, J. M. Pittard

Abstract: 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… ▽ More 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. △ Less

Submitted 11 February, 2013; originally announced February 2013.

Comments: 16 pages, 18 figures, accepted for publication in MNRAS

The 2.35 year itch of Cyg OB2 #9. II. Radio monitoring

Authors: R. Blomme, Y. Naze, D. Volpi, M. De Becker, R. K. Prinja, J. M. Pittard, E. R. Parkin, O. Absil

Abstract: Cyg OB2 #9 is one of a small set of non-thermal radio emitting massive O-star binaries. The non-thermal radiation is due to synchrotron emission in the colliding-wind region. Cyg OB2 #9 was only recently discovered to be a binary system and a multi-wavelength campaign was organized to study its 2011 periastron passage. We report here on the results of the radio observations obtained in this monito… ▽ More Cyg OB2 #9 is one of a small set of non-thermal radio emitting massive O-star binaries. The non-thermal radiation is due to synchrotron emission in the colliding-wind region. Cyg OB2 #9 was only recently discovered to be a binary system and a multi-wavelength campaign was organized to study its 2011 periastron passage. We report here on the results of the radio observations obtained in this monitoring campaign. We used the Expanded Very Large Array (EVLA) radio interferometer to obtain 6 and 20 cm continuum fluxes. The observed radio light curve shows a steep drop in flux sometime before periastron. The fluxes drop to a level that is comparable to the expected free-free emission from the stellar winds, suggesting that the non-thermal emitting region is completely hidden at that time. After periastron passage, the fluxes slowly increase. We introduce a simple model to solve the radiative transfer in the stellar winds and the colliding-wind region, and thus determine the expected behaviour of the radio light curve. From the asymmetry of the light curve, we show that the primary has the stronger wind. This is somewhat unexpected if we use the astrophysical parameters based on theoretical calibrations. But it becomes entirely feasible if we take into account that a given spectral type - luminosity class combination covers a range of astrophysical parameters. The colliding-wind region also contributes to the free-free emission, which can help to explain the high values of the spectral index seen after periastron passage. Combining our data with older Very Large Array (VLA) data allows us to derive a period P = 860.0 +- 3.7 days for this system. With this period, we update the orbital parameters that were derived in the first paper of this series. △ Less

Submitted 11 December, 2012; originally announced December 2012.

Comments: 10 pages, 4 figures, accepted for publication in A&A

Journal ref: A&A, 550, A90 (10p, 2013)

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

Authors: Y. Naze, L. Mahy, Y. Damerdji, H. A. Kobulnicky, J. M. Pittard, E. R. Parkin, O. Absil, R. Blomme

Abstract: 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 th… ▽ More 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. △ Less

Submitted 28 September, 2012; v1 submitted 25 September, 2012; originally announced September 2012.

Comments: 16p incl. online material, accepted by A&A

Journal ref: A&A, 546, A37 (14p, 2012)

Numerical simulations of shocks encountering clumpy regions

Authors: R. Alūzas, J. M. Pittard, T. W. Hartquist, S. A. E. G. Falle, R. Langton

Abstract: 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 m… ▽ More 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. △ Less

Submitted 26 June, 2012; originally announced June 2012.

Comments: 17 pages, 19 figures, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 425 (2012) 2212-2227

Cosmic Rays and Radiative Instabilities

Authors: T. W. Hartquist, A. Y. Wagner, S. A. E. G. Falle, J. M. Pittard, S. Van Loo

Abstract: 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 th… ▽ More 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. △ Less

Submitted 5 August, 2011; originally announced August 2011.

Comments: 5 pages, to appear in the Proceedings of Cosmic Rays & Their Interstellar Medium Environment CRISM-2011 (Montpellier, 26th June - 1st July 2011), to be published in an issue of MEMORIE della Società Astronomica Italiana, eds. A. Marcowith, A.Bykov, K.Ferriere, Y.Gallant, P.Hennebelle, T.Montmerle

The Integrated Diffuse X-ray Emission of the Carina Nebula Compared to Other Massive Star-forming Regions

Authors: Leisa K. Townsley, Patrick S. Broos, You-Hua Chu, Robert A. Gruendl, M. S. Oey, Julian M. Pittard

Abstract: The Chandra Carina Complex Project (CCCP) has shown that the Carina Nebula displays bright, spatially-complex soft diffuse X-ray emission. Here we `sum up' the CCCP diffuse emission work by comparing the global morphology and spectrum of Carina's diffuse X-ray emission to other famous sites of massive star formation with pronounced diffuse X-ray emission: M17, NGC 3576, NGC 3603, and 30 Doradus. A… ▽ More The Chandra Carina Complex Project (CCCP) has shown that the Carina Nebula displays bright, spatially-complex soft diffuse X-ray emission. Here we `sum up' the CCCP diffuse emission work by comparing the global morphology and spectrum of Carina's diffuse X-ray emission to other famous sites of massive star formation with pronounced diffuse X-ray emission: M17, NGC 3576, NGC 3603, and 30 Doradus. All spectral models require at least two diffuse thermal plasma components to achieve adequate spectral fits, a softer component with kT = 0.2--0.6 keV and a harder component with kT = 0.5--0.9 keV. In several cases these hot plasmas appear to be in a state of non-equilibrium ionization that may indicate recent and current strong shocks. A cavity north of the embedded giant HII region NGC 3576 is the only region studied here that exhibits hard diffuse X-ray emission; this emission appears to be nonthermal and is likely due to a recent cavity supernova, as evidenced by a previously-known pulsar and newly-discovered pulsar wind nebula also seen in this cavity. All of these targets exhibit X-ray emission lines that are not well-modeled by variable-abundance thermal plasmas and that might be attributed to charge exchange at the shock between the hot, tenuous, X-ray-emitting plasma and cold, dense molecular material; this is likely evidence for dust destruction at the many hot/cold interfaces that characterize massive star-forming regions. △ Less

Submitted 8 March, 2011; originally announced March 2011.

Comments: Accepted for the ApJS Special Issue on the Chandra Carina Complex Project (CCCP), scheduled for publication in May 2011. All 16 CCCP Special Issue papers are available at http://cochise.astro.psu.edu/Carina_public/special_issue.html through 2011 at least. 27 pages, 7 figures

X-ray emission from the double-binary OB-star system QZ Car (HD 93206)

Authors: E. R. Parkin, P. S. Broos, L. K. Townsley, J. M. Pittard, A. F. J. Moffat, Y. Naze, G. Rauw, L. M. Oskinova, W. L. Waldron

Abstract: X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra X-ray Observatory over a period of roughly 2 years are presented. The orbit of systems A (O9.7 I+b2 v, PA = 21 d) and B (O8 III+o9 v, PB = 6 d) are reasonably well sampled by the observations, allowing the origin of the X-ray emission to be examined in detail. The X-ray spectra can be well fitted by… ▽ More X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra X-ray Observatory over a period of roughly 2 years are presented. The orbit of systems A (O9.7 I+b2 v, PA = 21 d) and B (O8 III+o9 v, PB = 6 d) are reasonably well sampled by the observations, allowing the origin of the X-ray emission to be examined in detail. The X-ray spectra can be well fitted by an attenuated three temperature thermal plasma model, characterised by cool, moderate, and hot plasma components at kT ~ 0.2, 0.7, and 2 keV, respectively, and a circumstellar absorption of ~ 0.2 x 10^22 cm-2. Although the hot plasma component could be indicating the presence of wind-wind collision shocks in the system, the model fluxes calculated from spectral fits, with an average value of ~ 7 x 10^-13 erg s-1 cm-2, do not show a clear correlation with the orbits of the two constituent binaries. A semi-analytical model of QZ Car reveals that a stable momentum balance may not be established in either system A or B. Yet, despite this, system B is expected to produce an observed X-ray flux well in excess of the observations. If one considers the wind of the O8 III star to be disrupted by mass transfer the model and observations are in far better agreement, which lends support to the previous suggestion of mass-transfer in the O8 III + o9 v binary. We conclude that the X-ray emission from QZ Car can be reasonably well accounted for by a combination of contributions mainly from the single stars and the mutual wind-wind collision between systems A and B. △ Less

Submitted 8 March, 2011; v1 submitted 3 March, 2011; originally announced March 2011.

Comments: 11 pages, 7 figures. Accepted for the ApJS Special Issue on the Chandra Carina Complex Project (CCCP), scheduled for publication in May 2011. All 16 CCCP Special Issue papers are available at http://cochise.astro.psu.edu/Carina_public/special_issue.html through 2011 at least

Journal ref: ApJS, 194, 8 (11p (2011)

The Chandra Carina Complex Project: Deciphering the Enigma of Carina's Diffuse X-ray Emission

Authors: Leisa K. Townsley, Patrick S. Broos, You-Hua Chu, Marc Gagné, Gordon P. Garmire, Robert A. Gruendl, Kenji Hamaguchi, Mordecai-Mark Mac Low, Thierry Montmerle, Yaël Nazé, M. S. Oey, Sangwook Park, Robert Petre, Julian M. Pittard

Abstract: We present a 1.42 square degree mosaic of diffuse X-ray emission in the Great Nebula in Carina from the Chandra X-ray Observatory Advanced CCD Imaging Spectrometer camera. After removing >14,000 X-ray point sources from the field, we smooth the remaining unresolved emission, tessellate it into segments of similar apparent surface brightness, and perform X-ray spectral fitting on those tessellates… ▽ More We present a 1.42 square degree mosaic of diffuse X-ray emission in the Great Nebula in Carina from the Chandra X-ray Observatory Advanced CCD Imaging Spectrometer camera. After removing >14,000 X-ray point sources from the field, we smooth the remaining unresolved emission, tessellate it into segments of similar apparent surface brightness, and perform X-ray spectral fitting on those tessellates to infer the intrinsic properties of the X-ray-emitting plasma. By modeling faint resolved point sources, we estimate the contribution to the extended X-ray emission from unresolved point sources and show that the vast majority of Carina's unresolved X-ray emission is truly diffuse. Line-like correlated residuals in the X-ray spectral fits suggest that substantial X-ray emission is generated by charge exchange at the interfaces between Carina's hot, rarefied plasma and its many cold neutral pillars, ridges, and clumps. △ Less

Submitted 3 March, 2011; originally announced March 2011.

Comments: Accepted for the ApJS Special Issue on the Chandra Carina Complex Project (CCCP), scheduled for publication in May 2011. All 16 CCCP Special Issue papers are available at http://cochise.astro.psu.edu/Carina_public/special_issue.html through 2011 at least. 32 pages, 14 figures

Journal ref: ApJS, 194, 15 (23p, 2011)

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

Authors: Y. Naze, P. S. Broos, L. Oskinova, L. K. Townsley, D. Cohen, M. F. Corcoran, N. R. Evans, M. Gagne, A. F. J. Moffat, J. M. Pittard, G. Rauw, A. ud-Doula, N. R. Walborn

Abstract: 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 t… ▽ More 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. △ Less

Submitted 1 March, 2011; originally announced March 2011.

Comments: 30 pages, 9 figures. Accepted for the ApJS Special Issue on the Chandra Carina Complex Project (CCCP), scheduled for publication in May 2011. All 16 CCCP Special Issue papers are available at http://cochise.astro.psu.edu/Carina_public/special_issue.html through 2011 at least

Journal ref: ApJS, 194, 7 (2011)

An Introduction to the Chandra Carina Complex Project

Authors: Leisa K. Townsley, Patrick S. Broos, Michael F. Corcoran, Eric D. Feigelson, Marc Gagné, Thierry Montmerle, M. S. Oey, Nathan Smith, Gordon P. Garmire, Konstantin V. Getman, Matthew S. Povich, Nancy Remage Evans, Yaël Nazé, E. R. Parkin, Thomas Preibisch, Junfeng Wang, Scott J. Wolk, You-Hua Chu, David H. Cohen, Robert A. Gruendl, Kenji Hamaguchi, Robert R. King, Mordecai-Mark Mac Low, Mark J. McCaughrean, Anthony F. J. Moffat , et al. (10 additional authors not shown)

Abstract: The Great Nebula in Carina provides an exceptional view into the violent massive star formation and feedback that typifies giant HII regions and starburst galaxies. We have mapped the Carina star-forming complex in X-rays, using archival Chandra data and a mosaic of 20 new 60ks pointings using the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer, as a testbed for understanding recent… ▽ More The Great Nebula in Carina provides an exceptional view into the violent massive star formation and feedback that typifies giant HII regions and starburst galaxies. We have mapped the Carina star-forming complex in X-rays, using archival Chandra data and a mosaic of 20 new 60ks pointings using the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer, as a testbed for understanding recent and ongoing star formation and to probe Carina's regions of bright diffuse X-ray emission. This study has yielded a catalog of properties of >14,000 X-ray point sources; >9800 of them have multiwavelength counterparts. Using Chandra's unsurpassed X-ray spatial resolution, we have separated these point sources from the extensive, spatially-complex diffuse emission that pervades the region; X-ray properties of this diffuse emission suggest that it traces feedback from Carina's massive stars. In this introductory paper, we motivate the survey design, describe the Chandra observations, and present some simple results, providing a foundation for the 15 papers that follow in this Special Issue and that present detailed catalogs, methods, and science results. △ Less

Submitted 23 February, 2011; originally announced February 2011.

Comments: Accepted for the ApJS Special Issue on the Chandra Carina Complex Project (CCCP), scheduled for publication in May 2011. All 16 CCCP Special Issue papers are available at http://cochise.astro.psu.edu/Carina_public/special_issue.html through 2011 at least. 43 pages; 18 figures

Journal ref: ApJS, 194, 1 (28p, 2011)

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

Authors: R. Fahed, A. F. J. Moffat, J. Zorec, T. Eversberg, A. N. Chené, F. Alves, W. Arnold, T. Bergmann, L. F. Gouveia Carreira, F. Marques Dias, A. Fernando, J. Sanchez Gallego, T. Hunger, J. H. Knapen, R. Leadbeater, T. Morel, G. Rauw, N. Reinecke, J. Ribeiro, N. Romeo, E. M. dos Santos, L. Schanne, O. Stahl, Ba. Stober, Be. Stober , et al. (7 additional authors not shown)

Abstract: 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-los… ▽ More 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. △ Less

Submitted 7 January, 2011; originally announced January 2011.

Comments: 5 pages, 5 figures, Published in proceedings of "Stellar Winds in Interaction", editors T. Eversberg and J.H. Knapen. Full proceedings volume is available on http://www.stsci.de/pdf/arrabida.pdf

Theory and Models of Colliding Stellar Winds

Authors: J. M. Pittard

Abstract: I discuss some of the important aspects of the phenomena of colliding winds in massive binary star systems with a particular focus on WR140. I discuss some of the important aspects of the phenomena of colliding winds in massive binary star systems with a particular focus on WR140. △ Less

Submitted 7 January, 2011; originally announced January 2011.

Comments: 11 pages, 5 figures, Published in proceedings of "Stellar Winds in Interaction", editors T. Eversberg and J.H. Knapen. Full proceedings volume is available on http://www.stsci.de/pdf/arrabida.pdf