The formation, properties, and evolution of massive stars remain subject to considerable uncertainty; impacting on fields as diverse as galactic feedback and the nature of the progenitors of both electromagnetic and gravitational wave transients. The clusters many such stars reside within provide a unique laboratory for addressing these issues, and in this work we provide a comprehensive stellar census of Westerlund 1 to underpin such efforts. 69 new members were identified via I-band spectroscopy, yielding a total cluster population of 166 stars with initial masses of ~25Msun to ~50Msun, with more massive stars already lost to supernova. The stellar population follows a smooth and continuous morphological sequence from late-O giant through to OB supergiant. Subsequently, the progression bifurcates, with one branch yielding mid- to late-B hypergiants and cool super-/hypergiants, and the other massive blue stragglers, prior to a diverse population of H-depleted Wolf-Rayets. A substantial population of O-type stars with anomalously broad Paschen series lines are seen, a property which we attribute to binarity. Binary interaction is clearly required to yield the uniquely rich cohort of hypergiants, which includes both mass-stripped primaries and rejuvenated secondaries/stellar mergers. As a consequence future observations of Wd1 and similar stellar aggregates hold out the prospect of characterising both single- and binary- evolutionary channels for massive stars and determining their relative contributions. This in turn will permit the physical properties of such objects at the point of core-collapse to be predicted; of direct relevance for understanding the formation of relativistic remnants such as the magnetars associated with Wd1 and other young massive clusters (Abridged).

Despite the first detection of X-rays from massive stars occurring four decades ago, the physical dependence of the emission mechanism(s) on the underlying stellar and binary properties of the emitters remains uncertain. The young massive cluster Westerlund 1 provides an ideal testbed for understanding this phenomenon, with over 50 cluster members detected in historical X-ray observations. In the decade since these data were obtained, significant new multi-epoch observations of the cluster have been made, allowing a fundamental reappraisal of the nature of both X-ray bright and dark stars. A total of 45 X-ray sources within Wd1 now have precise spectral classifications. These observations identify 16 candidate and confirmed massive binaries; by comparison 22 stars have X-ray properties that imply a contribution from a wind collision zone. X-ray emission appears confined to O9-B0.5 supergiants, Wolf-Rayets and a small group of highly luminous interacting/post-interaction binaries. Despite their presence in large numbers, no emission is seen from earlier, less evolved O stars or later, cooler B super-/hypergiants. We suppose that the lack of X-ray emission from O giants is due to their comparatively low bolometric luminosities if, as expected, they follow the canonical L_X/L_bol relation for hot stars. The transition away from X-ray emission for OB supergiants occurs at the location of the bistability jump; we speculate that below this limit, stellar wind velocities are insufficient for internal, X-ray emitting shocks to form. Our results are consistent with recent findings that massive binaries are not uniformly brighter than single stars of comparable luminosity and spectral type, although it is noteworthy that the brightest and hardest stellar X-ray sources within Wd1 are all either confirmed or candidate massive, interacting/post-interaction binaries (Abridged).

Observational studies suggest that many OB stars are found within binary systems which may be expected to interact during their lifetimes. Significant mass transfer or merger of both components will modify evolutionary pathways, facilitating the production of exceptionally massive stars which will present as blue stragglers. Identification and characterisation of such objects is crucial if the outcomes of binary evolutionary channels are to be quantified. The massive cluster Westerlund 1 hosts a rich population of X-ray bright stars where the emission is thought to derive from the wind collision zones of massive binaries. Selected on this basis, we present the results of a multiwavelength analysis of the X-ray luminous O stars Wd1-27 and -30a. We find both to be early/mid-O hypergiants with luminosities, temperatures and masses significantly in excess of other early stars within Wd1, hence qualifying as massive blue stragglers. The nature of Wd1-27 remains uncertain but the detection of radial velocity changes and the X-ray properties of Wd1-30a suggest that it is a binary. Analysis of Gaia proper motion and parallactic data indicates that both are cluster members; we also provide a membership list for Wd1 based on this analysis. The presence of hypergiants of spectral types O to M within Wd1 cannot be understood via single-star evolution. We suppose that the early-B and mid-O hypergiants formed via binary-induced mass-stripping of the primary and mass-transfer to the secondary, respectively. This implies that for a subset of objects massive star-formation may be regarded as a two-stage process, with binary-driven mass-transfer or merger yielding stars with masses significantly in excess of their initial birth mass (Abridged).

Massive stars play an important role in both cluster and galactic evolution and the rate at which they lose mass is a key driver of both their own evolution and their interaction with the environment up to and including their SNe explosions. Young massive clusters provide an ideal opportunity to study a co-eval population of massive stars. We performed 3mm continuum observations with the Atacama Large Millimetre/submillimetre Array of the Galactic cluster Westerlund 1, to study the constituent massive stars and determine mass-loss rates for the diverse post-main sequence population. We detected emission from 50 stars in Westerlund 1, comprising all 21 Wolf-Rayets within the field of view, eight cool and 21 OB super-/hypergiants. Emission nebulae were associated with a number of the cool hypergiants while, unexpectedly, a number of hot stars also appear spatially resolved. We measured the mass-loss rates for a unique population of massive post-main sequence stars at every stage of evolution, confirming a significant increase as stars transition from OB supergiant to WR states. The range of spectral types exhibited provides a critical test of radiatively driven wind theory and the reality of the bi-stability jump. The extreme mass-loss rate inferred for the interacting binary Wd1-9 in comparison to other cluster members confirmed the key role binarity plays in massive stellar evolution. The presence of compact nebulae around a number of OB and WR stars is unexpected; by analogy to the cool super-/hypergiants we attribute this to confinement and sculpting of the stellar wind via interaction with the intra-cluster medium/wind. Given the morphology of core collapse SNe depend on the nature of the pre-explosion circumstellar environment, if this hypothesis is correct then the properties of the explosion depend not just on the progenitor, but also the environment in which it is located.

The first soft gamma-ray repeater was discovered over three decades ago, and subsequently identified as a magnetar. However there is currently no consenus on the formation channel(s) of these objects. The presence of a magnetar in the starburst cluster Westerlund 1 implies a >40Msun progenitor, favouring formation in a binary that was disrupted at supernova. To test this hypothesis we searched for the putative pre-SN companion via a radial velocity survey to identify high-velocity runaways, with subsequent atmospheric analysis of the resultant candidate, Wd1-5. Wd1-5 appears overluminous for its spectroscopic mass and we find evidence of He- and N-enrichment, O-depletion, and critically C-enrichment, a combination of properties that is difficult to explain for a single star. We infer a pre-SN history for Wd1-5 which supposes an initial close binary comprising two stars of comparable (~41Msun+35Msun) masses. Efficient mass transfer from the initially more massive component leads to the mass-gainer evolving more rapidly, initiating luminous blue variable/common envelope evolution. Reverse, wind-driven mass transfer during its subsequent WC Wolf-Rayet phase leads to the carbon pollution of Wd1-5, before a type Ibc supernova disrupts the binary system. Under the assumption of a physical association between Wd1-5 and J1647-45, the secondary is identified as the magnetar progenitor; its common envelope evolutionary phase prevents spin-down of its core prior to SN and the seed magnetic field for the magnetar forms either in this phase or during the earlier episode of mass transfer in which it was spun-up. Our results suggest that binarity is a key ingredient in the formation of at least a subset of magnetars by preventing spin-down via core-coupling and potentially generating a seed magnetic field.

The location of the supergiant (sg)B[e] star Wd1-9 within the coeval starburst cluster Westerlund 1 allows for its placement in a detailed post-Main Sequence evolutionary scheme and hence we have utilised a comprehensive multiwavelength dataset to determine its physical properties and relation to other sgB[e] stars and the evolved stellar population of Wd1. Wd1-9 is found to exhibit the rich optical emission line spectrum that is characteristic of sgB[e] stars, while dust mass, composition and disc geometry determined from mid-IR data resemble those of other well known examples such as R126. Extreme historical and ongoing mass loss is inferred from radio observations, while the X-ray properties of Wd1-9 imply the presence of high temperature plasma within the system and are directly comparable to a number of confirmed short-period colliding wind binaries within Wd1. The most complete explanation for the observational properties of Wd1-9 is that it is a massive interacting binary currently undergoing, or recently exited from, rapid Roche-lobe overflow, supporting the hypothesis that binarity mediates the formation of (a subset of) sgB[e] stars. The mass loss rate of Wd1-9 is consistent with such an assertion, while viable progenitor and descendent systems are present within Wd1 and other examples of sgB[e] binaries have previously been found. Moreover, the rarity of sgB[e] stars - only two examples are identified from a census of ~68 young massive Galactic clusters and associations containing ~600 post-Main Sequence stars - is explicable given the rapidity (~10^4yr) expected for this phase of massive binary evolution.

Luminous blue variables (LBVs) are expected to play an important role in massive stellar evolution as well as being the progenitors of some of the most luminous supernovae known. In this paper we provide a multiwavelength study of the population of (candidate) LBVs identified within M33. New spectra provide an observational baseline of >4yr with respect to published data, which is well suited to identifying LBV outbursts. Multi-epoch optical and mid-IR surveys of M33 further constrain the variability of the sample and permit a search for dusty circumstellar ejecta. Spectroscopic and photometric variability appears common amongst the sample, although in many cases further observations will be required to determine its physical origin. Nevertheless, we report a new outburst of M33 Var C, while the transition of the WNLh star B517 to a cooler B supergiant phase between 1993-2010 confirms an LBV classification. Proof-of-concept quantitative analysis is provided for Romano's star; the results being consistent with the finding that its bolometric luminosity varies during its LBV excursions. The combination of the temperature and luminosity of two stars, the B hypergiant [HS80] 110A and the cool hypergiant B324, appears to be in violation of the empirical Humphreys-Davidson limit. Mid-IR observations demonstrate that a number of candidates appear associated with hot circumstellar dust, although no objects as extreme as Eta Carinae are identified. The combined multiwavelength dataset suggests that the population of LBVs studied is contaminated by stars demonstrating the B[e] phenomenon. Of these, a subset of optically faint, low luminosity stars associated with hot dust are of particular interest since they appear similar to the likely progenitors of SN 2008S and the 2008 NGC300 transient, albeit suffering less intrinsic extinction. [ABRIDGED]

Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs); Cyg OB2 #12, zeta Sco, HD190603 and BP Cru. These are combined with an analysis of their long-term spectroscopic and photometric behaviour in order to determine their evolutionary status. The long-term datasets revealed that they are remarkably stable over long periods (>40yr), with the possible exception of zeta Sco prior to the 20th century, in contrast to the typical excursions that characterise luminous blue variables (LBVs). Zeta Sco, HD190603 and BP Cru possess physical properties intermediate between B supergiants and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such stars (for initial masses in the range ~30-60Msun). In contrast, while the wind properties of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (~110Msun) and comparatively low temperature means it cannot be accommodated in such a scheme. Likewise, despite its co-location with several LBVs above the Humphreys-Davidson (HD) limit, the lack of long term variability and its unevolved chemistry apparently excludes such an identification. Since such massive stars are not expected to evolve to such cool temperatures, the properties of Cyg OB2 #12 are difficult to understand under current evolutionary paradigms. [ABRIDGED]

NGC346-013 is a peculiar double-lined eclipsing binary in the Small Magellanic Cloud discovered by the VLT-FLAMES survey of massive stars. Spectra obtained with VLT-FLAMES are used to construct a radial velocity curve and photometry obtained with the Faulkes Telescope South is then used to derive orbital parameters, while spectra of the secondary are compared with synthetic spectra from TLUSTY model atmospheres. The orbital period is found to be 4.20381(12) days, with masses of 19.1+/-1.0 and 11.9+/-0.6 Msun. The primary is a rapidly rotating late-O dwarf while the secondary, an early-B giant, displays near-synchronous rotation and has filled its Roche lobe, implying that it was originally the more massive component with recent mass transfer `spinning up' the primary to near-critical rotation. Comparison with synthetic spectra finds temperatures of 34.5kK and 24.5kK for the primary and secondary respectively, with the nitrogen abundance of the secondary enhanced compared to baseline values for the SMC, consistent with the predictions of models of interacting binaries. NGC346-013 likely evolved via non-conservative mass transfer in a system with initial masses ~22+15Msun, with the well-constrained orbital solution and atmospheric parameters making it an excellent candidate for tailored modelling with binary evolution codes. This system will form a cornerstone in constraining the physics of thermal timescale mass transfer, and the associated mass transfer efficiency, in massive close binary systems.

There is growing evidence that a treatment of binarity amongst OB stars is essential for a full theory of stellar evolution. However the binary properties of massive stars - frequency, mass ratio and orbital separation - are still poorly constrained. In order to address this shortcoming we have undertaken a multiepoch spectroscopic study of the stellar population of the young massive cluster Westerlund 1. In this paper we present an investigation into the nature of the dusty Wolf-Rayet star and candidate binary W239. We determined an orbital period of only 5.05 days - making it one of the most compact WC binaries yet identified - while we interpret a significant near-IR flare between 2004-6 as evidence for a third massive stellar component in the system in a wide eccentric orbit. The presence of a near-IR excess characteristic of hot dust at every epoch is consistent with the expectation that the subset of persistent dust forming WC stars are short period binaries. In conjunction with previous results, the simultaneous presence of both short period Wolf-Rayet binaries and cool hypergiants within Wd 1 provides compelling evidence for a bifurcation in the post-Main Sequence evolution of massive stars due to binarity. Short period O+OB binaries will evolve directly to the Wolf-Rayet phase, either due to an episode of binary mediated mass loss or via chemically homogenous evolution. Conversely, long period binaries and single stars will instead undergo a red loop across the HR diagram via a cool hypergiant phase. Future analysis of the full spectroscopic dataset for Wd 1 will constrain the proportion of massive stars experiencing each pathway; hence quantifying the importance of binarity in massive stellar evolution up to and beyond supernova and the resultant production of relativistic remnants and X-ray binaries. [ABRIDGED]

Recent IR surveys of the Galactic plane have revealed a large number of candidate Luminous Blue Variables. In order to verify these classifications we have been undertaking a long term spectroscopic and photometric monitoring campaign supplemented with tailored non-LTE model atmosphere analysis. Here we present a brief overview of selected aspects of this program, highlighting the prospects for identification, classification and quantitative analysis of LBVs in the near-IR spectral window.

We present initial results from a long-baseline radial velocity survey for massive binaries in the cluster Westerlund 1. Four systems are examined: the dust-producing WC binary W239, the double-lined eclipsing binary W13, and the single-lined B0 supergiants W43a and W3003. Finally, the evolutionary implications for the population of massive stars in Westerlund 1 are discussed.

Westerlund 1 is a young, massive Galactic starburst cluster that contains a rich coeval population of Wolf-Rayet stars, hot- and cool-phase transitional supergiants, and a magnetar. We use spectroscopic and photometric observations of the eclipsing double-lined binary W13 to derive dynamical masses for the two components, in order to determine limits for the progenitor masses of the magnetar CXOU J164710.2-455216 and the population of evolved stars in Wd1. W13 has an orbital period of 9.2709+/-0.0015 days and near-contact configuration. The shallow photometric eclipse rules out an inclination greater than 65 degrees, leading to lower limits for the masses of the emission-line optical primary and supergiant optical secondary of 21.4+/-2.6Msun and 32.8+/-4.0Msun respectively, rising to 23.2 +3.3/-3.0Msun and 35.4 +5.0/-4.6 Msun for our best-fit inclination 62 +3/-4 degrees. Comparison with theoretical models of Wolf-Rayet binary evolution suggest the emission-line object had an initial mass in excess of 35Msun, with the most likely model featuring highly non-conservative late-Case-A/Case-B mass transfer and an initial mass in excess of 40Msun. This confirms the high magnetar progenitor mass inferred from its membership in Wd1, and represents the first dynamical constraint on the progenitor mass of any magnetar. The red supergiants in Wd1 must have similar progenitor masses to W13 and are therefore amongst the most massive stars to undergo a red supergiant phase, representing a challenge for population models that suggest stars in this mass range end their redwards evolution as yellow hypergiants. [ABRIDGED]

After leaving the main sequence, massive stars undergo complex evolution, still poorly understood. With a population of 100s OB stars, the starburst cluster Westerlund~1 offers an unparallelled environment to study their evolutionary tracks. We used the VLT/FORS2 to obtain intermediate-resolution spectroscopy over the range 5800-9000A of ~ 100 likely members. We developed criteria for their spectral classification using only spectral features in the range observed. We discuss these criteria, useful for spectral classification of early-type stars in the GAIA spectral region, in the appendix. Using these criteria, we obtain spectral classifications, probably accurate to one subtype, for 57 objects, most of which had no previous classification or a generic classification. We identify more than 50 objects as OB supergiants. We find almost 30 luminous early-B supergiants and a number of less luminous late-O supergiants. In addition, we find a few mid B supergiants with very high luminosity, some of them displaying signs of heavy mass loss. All these stars form a sequence compatible with theoretical evolutionary tracks. In addition, two early B supergiants also show indication of heavy mass loss and may represent the evolutionary phase immediately prior to the Wolf-Rayet stage. We investigate cluster properties using the spectral types and existing photometry. We find that the reddening law to the cluster does not deviate strongly from standard, even though extinction is quite variable, with an average value A_v=10.8. Though evolutionary tracks for high-mass stars are subject to large uncertainties, our data support an age of >~5Myr and a distance ~5kpc for Westerlund 1. The spectral types observed are compatible with a single burst of star formation (the age range is very unlikely to be >1Myr).[ABRIDGED]

Massive stars exhibit spectroscopic and photometric variability over a wide range of timescales. However the physical mechanisms driving this behaviour remain poorly understood. Westerlund 1 presents an ideal laboratory for studying these processes in a rich, coeval population of post-main sequence stars and we present a pathfinding study aimed at characterising their variability utilising the large body of data that has accumulated since the 1960s. Spectroscopic variability attributable to both wind asphericity and photospheric pulsations was present amongst both the hot and cool stellar populations. Given the limitations imposed by the data, we were unable to determine the physical origin of the wind structure inferred for the OB supergiants, although the inhomogineities in the winds of the Wolf Rayets are likely driven by binary interactions and, conversely, by pulsations in at least one of the cool hypergiants. Photospheric pulsations were found for stars ranging from spectral types as early as O9 I through to the mid F Ia+ yellow hypergiants - with a possible dependence on the luminosity class amongst the OB stars. The spectroscopically variable red supergiants (M2-5 Ia) are also potential pulsators but require further observations to confirm this hypothesis. Determination of the amplitude and periodicity of these pulsations as a function of temperature, luminosity and evolutionary state holds open the prospect of identifying the physical mechanisms driving the instabilities that constrain and define stellar evolution in the upper reaches of the HR diagram, while the presence of highly luminous yellow hypergiants and red supergiants within Wd1 place strong constraints on post-main sequence evolutionary pathways, apparently contradicting current theoretical predictions for >25Msun stars at solar metallicites. [ABRIDGED]

The massive post-Main Sequence star W243 in the galactic starburst cluster Westerlund 1 has undergone a spectral transformation from a B2Ia supergiant devoid of emission features in 1981 to an A-type supergiant with a rich emission-line spectrum by 2002/03. We used VLT/UVES and VLT/FLAMES to obtain high-resolution spectra on six epochs in 2003/04 (UVES) and ten epochs in 2008/09 (FLAMES). These spectra are used alongside other VLT/FLAMES and NTT/EMMI spectra to follow the evolution of W243 from 2002 to 2009. W243 displays a complex, time-varying spectrum with emission lines of Hydrogen, Helium and Lyman-pumped metals, forbidden lines of Nitrogen and Iron, and a large number of absorption lines from neutral and singly-ionized metals. Many lines are complex emission/absorption blends. The LBV has a current temperature of ~8500K (spectral type A3Ia+), and displays signs of photospheric pulsations and weak episodic mass loss. Nitrogen is highly overabundant, with Carbon and Oxygen depleted, indicative of surface CNO-processed material and considerable previous mass-loss, although current time-averaged mass-loss rates are low. The emission-line spectrum forms at large radii, when material lost by the LBV in a previous mass-loss event is ionized by an unseen hot companion. Monitoring of the near-infrared spectrum suggests that the star has not changed significantly since it finished evolving to the cool state, close to the Humphreys-Davidson limit, in early 2003. [ABRIDGED]

Multiwavelength observations of the young massive cluster Westerlund 1 have revealed evidence for a large number of OB supergiant and Wolf-Rayet binaries. However, in most cases these findings are based on secondary binary characteristics such as hard X-ray emission and/or non-thermal radio spectra and hence provide little information on binary properties such as mass ratio and orbital period. To overcome this shortcoming we have initiated a long temporal baseline, VLT/FLAMES+GIRAFFE multi-epoch radial velocity survey that will provide the first direct constraints on these parameters. This study presents first results from twenty of the most luminous supergiant stars in Wd1. Statistically significant radial velocity changes are detected in ~60% of targets. W43a is identified as a short-period binary, while W234 and the newly-identified cluster member W3003 are probable binaries and W2a is a strong binary candidate. The cool hypergiants W243 and W265 display photospheric pulsations, while a number of early-mid B supergiants display radial velocity changes that we cannot distinguish between orbital or photospheric motion in our initial short-baseline survey. When combined with existing observations, we find 30% of our sample to be binary (6/20) while additional candidate binaries support a binary fraction amongst Wd1 supergiants in excess of ~40% [ABRIDGED]

We characterise the variability in the physical properties of the luminous blue variable AFGL2298 between 1989-2008. In conjunction with published data from 1989-2001, we have undertaken a long term (2001-2008) near-IR spectroscopic and photometric observational campaign for this star and utilise a non-LTE model atmosphere code to interpret these data. We find AFGL2298 to have been highly variable during the two decades covered by the observational datasets. Photometric variations of >1.6 mag have been observed in the JHK wavebands; however, these are not accompanied by correlated changes in near-IR colour. Non-LTE model atmosphere analysis of 4 epochs of K band spectroscopy obtained between 2001-7 suggests that the photometric changes were driven by expansion and contraction of the stellar photosphere accompanied by comparatively small changes in the stellar temperature. Unclumped mass loss rates throughout this period were modest and directly comparable to those of other highly luminous LBVs. However, the bolometric luminosity of AFGL2298 appears to have varied by at least a factor of ~2 between 1989-2008, with it being one of the most luminous stars in the Galaxy during maximum. Comparison to other LBVs that have undergone non bolometric luminosity conserving `eruptions' shows such events to be heterogeneous, with AFGL2298 the least extreme example. These results - and the diverse nature of both the quiescent LBVs and associated ejecta - may offer support to the suggestion that more than one physical mechanism is responsible for such behaviour. [ABRIDGED]

We aim to characterise the properties of a third massive, red supergiant dominated galactic cluster. To accomplish this we utilised a combination of near/mid-IR photometry and spectroscopy to identify and classify the properties of cluster members, and statistical arguments to determine the mass of the cluster. We found a total of 16 strong candidates for cluster membership, for which formal classification of a subset yields spectral types from K3-M4 Ia and luminosities between log(L/L_sun)~4.5-4.8 for an adopted distance of 6+/-1 kpc. For an age in the range of 16-20 Myr, the implied mass is 2-4x10^4 M_sun, making it one of the most massive young clusters in the Galaxy. This discovery supports the hypothesis that a significant burst of star formation occurred at the base of Scutum-Crux arm between 10-20 Myr ago, yielding a stellar complex comprising at least ~10^5M_sun of stars (noting that since the cluster identification criteria rely on the presence of RSGs, we suspect that the true stellar yield will be significantly higher). We highlight the apparent absence of X-ray binaries within the star formation complex and finally, given the physical association of at least two pulsars with this region, discuss the implications of this finding for stellar evolution and the production and properties of neutron stars.

We present the results of high-resolution AP3M+SPH simulations of merging clusters of galaxies. We find that the compression and shocking of the core gas during a merger can lead to large increases in bolometric X-ray luminosities and emission-weighted temperatures of clusters. Cooling flows are completely disrupted during equal-mass mergers, with the mass deposition rate dropping to zero as the cores of the clusters collide. The large increase in the cooling time of the core gas strongly suggests that cooling flows will not recover from such a merger within a Hubble time. Mergers with subclumps having 1/8th of the mass of the main cluster are also found to disrupt a cooling flow if the merger is head-on. However, in this case the entropy injected into the core gas is rapidly radiated away and the cooling flow restarts within a few Gyr of the merger. Mergers in which the subcluster has an impact parameter of 500kpc do not disrupt the cooling flow, although the mass deposition rate is reduced by \~30%. Finally, we find that equal mass, off-centre mergers can effectively mix gas in the cores of clusters, while head on mergers lead to very little mixing. Gas stripped from the outer layers of subclumps results in parts of the outer layers of the main cluster being well mixed, although they have little effect on the gas in the core of the cluster. None of the mergers examined here resulted in the ICM being well mixed globally.

We adapt the Smoothed-Particle Hydrodynamics technique to allow a multiphase fluid in which SPH particles of widely differing density may be freely intermixed. Applications include modelling of galaxy formation and cooling flows.