We report the discovery of a handful of optical hydrogen-poor knots in the central part of an extended infrared nebula centred on the [WO1] star WR 72, obtained by spectroscopic and imaging observations with the Southern African Large Telescope (SALT). Wide-field Infrared Survey Explorer (WISE) images of the nebula show that it is composed of an extended almost circular halo (of $\approx6$ arcmin or $\approx2.4$ pc in diameter) and an elongated and apparently bipolar inner shell (of a factor of six smaller size), within which the knots are concentrated. Our findings indicate that WR 72 is a new member of the rare group of hydrogen-poor planetary nebulae, which may be explained through a very late thermal pulse of a post-AGB star, or by a merger of two white dwarfs.

We report the discovery of a square axisymmetric circumstellar nebula around the emission-line star HD 93795 in archival Spitzer Space Telescope 24 micron data. We classify HD 93795 as an B9 Ia star using optical spectra obtained with the Southern African Large Telescope (SALT). A spectral analysis carried out with the stellar atmosphere code FASTWIND indicates that HD 93795 only recently left the main sequence and is evolving redward for the first time. We discuss possible scenarios for the origin of the nebula and suggest that HD 93795 was originally a binary system and that the nebula was formed because of merger of the binary components. We also discuss a discrepancy between distance estimates for HD 93795 based on the Gaia data and the possible membership of this star of the Car OB1 association, and conclude that HD 93795 could be at the same distance as Car OB1.

We report the discovery of a spiral-like nebula with the Wide-field Infrared Survey Explorer (WISE) and the results of optical spectroscopy of its associated star TYC 8606-2025-1 with the Southern African Large Telescope (SALT). We find that TYC 8606-2025-1 is a G8 III star of $\approx3$ M$_\odot$, showing a carbon depletion by a factor of two and a nitrogen enhancement by a factor of three. We also derived an excess of s-process elements, most strongly for barium, which is a factor of three overabundant, indicating that TYC 8606-2025-1 is a mild barium star. We thereby add a new member to the small group of barium stars with circumstellar nebulae. Our radial velocity measurements indicate that TYC 8606-2025-1 has an unseen binary companion. The advanced evolutionary stage of TYC 8606-2025-1, together with the presence of a circumstellar nebula, implies an initial mass of the companion of also about 3 M$_\odot$. We conclude that the infrared nebula, due to its spiral shape, and because it has no optical counterpart, was ejected by the companion as a consequence of a very late thermal pulse, during about one orbital rotation.

Gravitational wave emission can lead to the coalescence of close pairs of compact objects orbiting each other. For the case of neutron stars such mergers may yield masses above the Tolman-Oppenheimer-Volkoff limit, leading to the formation of black holes. For the case of white dwarfs the merger product may exceed the Chandrasekhar limit, leading either to a thermonuclear explosion as Type Ia supernova, or to a collapse forming a neutron star. If a Type Ia supernova explosion is avoided, the merger of two massive white dwarfs is expected to form a hydrogen- and helium-free circumstellar nebula with a hot and luminous, rapidly rotating and highly magnetized central star for several 10,000 yr before its final collapse. Here we report the discovery of a hot star with an emission line dominated spectrum in the centre of a circular mid-infrared nebula. Both the star and the nebula appear to be free of hydrogen and helium. Our tailored stellar atmosphere and wind models indicate a stellar surface temperature of about 200,000 K, and a record outflow velocity of 16,000 km/s. This extreme velocity, together with the derived mass outflow rate, imply rapid stellar rotation and a strong magnetic field aiding the wind acceleration. The Gaia distance of the star leads to a luminosity of 10^{4.5} Lsun, which matches models of the post-merger evolution of super-Chandrasekhar mass white dwarfs. The high stellar temperature and the nebular size argue for a short remaining lifetime of the star, which will produce a bright optical and high-energy transient upon collapse. Our observations indicate that super-Chandrasekhar mass white dwarf mergers can indeed avoid a thermonuclear explosion as Type Ia supernova, and provide empirical evidence for magnetic field generation in stellar mergers.

The Bubble Nebula (or NGC 7635) is a parsec-scale seemingly spherical wind-blown bubble around the relatively unevolved O star BD+60$^\circ$2522. The small dynamical age of the nebula and significant space velocity of the star suggest that the Bubble Nebula might be a bow shock. We have run 2D hydrodynamic simulations to model the interaction of the central star's wind with the interstellar medium (ISM). The models cover a range of possible ISM number densities of $n=50-200 {\rm cm}^{-3}$ and stellar velocities of $v_{\star}=20-40$ km s$^{-1}$. Synthetic H$\alpha$ and 24 $\mu$m emission maps predict the same apparent spherical bubble shape with quantitative properties similar to observations. The synthetic maps also predict a maximum brightness similar to that from the observations and agree that the maximum brightness is at the apex of the bow shock. The best-matching simulation had $v_{\star}\approx20$ km s$^{-1}$ into an ISM with $n\sim100 {\rm cm}^{-3}$, at an angle of 60$^\circ$ with respect to the line of sight. Synthetic maps of soft ($0.3-2$ keV) and hard ($2-10$ keV) X-ray emission show that the brightest region is in the wake behind the star and not at the bow shock itself. The unabsorbed soft X-rays have luminosity $\sim10^{32}-10^{33}$ erg s$^{-1}$. The hard X-rays are fainter, luminosity $\sim 10^{30} - 10^{31}$ erg s$^{-1}$, and may be too faint for current X-ray instruments to successfully observe. Our results imply that the O star creates a bow shock as it moves through the ISM and in turn creates an asymmetric bubble visible at optical and infrared wavelengths, and predicted to be visible in X-rays. The Bubble Nebula does not appear to be unique, it could be just a favourably oriented very dense bow shock. The dense ISM surrounding BD+60$^\circ$2522 and its strong wind suggest that it could be a good candidate for detecting non-thermal emission.

We report the discovery of a circular shell centred on the Be X-ray binary (BeXB) SXP 1323 in the Small Magellanic Cloud (SMC). The shell was detected in an Halpha image obtained with the Very Large Telescope (VLT). Follow-up spectroscopy with the Southern African Large Telescope (SALT) showed that the shell expands with a velocity of $\approx$ 100 km/s and that its emission is due to shock excitation. We suggest that this shell is the remnant of the supernova explosion that led to the formation of the SXP 1323's neutron star $\approx$ 40 000 yr ago. SXP 1323 represents the second known case of a BeXB associated with a supernova remnant (the first one is SXP 1062). Interestingly, both these BeXBs harbour long period pulsars and are located in a low-metallicity galaxy.

We report the discovery of two mid-infrared nebulae in the northern hemisphere with the Wide-field Infrared Survey Explorer and the results of optical spectroscopy of their central stars, BD+60 2668 (composed of two components, separated from each other by \approx 3 arcsec) and ALS 19653, with the Calar Alto 3.5-m telescope and the Southern African Large Telescope (SALT), respectively. We classify the components of BD+60 2668 as stars of spectral types B0.5 II and B1.5 III. ALS 19653 is indicated in the SIMBAD data base as a planetary nebula, while our observations show that it is a massive B0.5 Ib star, possibly in a binary system. Using the stellar atmosphere code FASTWIND, we derived fundamental parameters of the three stars as well as their surface element abundances, implying that all of them are either on the main sequence or only recently left it. This provides further evidence that massive stars can produce circumstellar nebulae while they are still relatively unevolved. We also report the detection of optical counterparts to the mid-infrared nebulae and a second, more extended optical nebula around ALS 19653, and present the results of SALT spectroscopy of both nebulae associated with this star. The possible origin of the nebulae is discussed.

Many high-mass X-ray binaries (HMXBs) are runaways. Stellar wind and radiation of donor stars in HMXBs along with outflows and jets from accretors interact with the local interstellar medium and produce curious circumstellar structures. Several such structures are presented and discussed in this contribution.

We show that Westerlund 1 is the likely birth cluster of the runaway bona fide LBV MN44.

We report the results of our study of the high-velocity (\approx160 km/s) runaway O star CPD-64 2731 and its associated horseshoe-shaped nebula discovered with the Wide-field Infrared Survey Explorer. Spectroscopic observations with the Southern African Large Telescope and spectral analysis indicate that CPD-64 2731 is a fast-rotating main-sequence O5.5 star with enhanced surface nitrogen abundance. We derive a projected rotational velocity of \approx300 km/s which is extremely high for this spectral type. Its kinematic age of \approx6 Myr, assuming it was born near the Galactic plane, exceeds its age derived from single star models by a factor of two. These properties suggest that CPD-64 2731 is a rejuvenated and spun-up binary product. The geometry of the nebula and the almost central location of the star within it argue against a pure bow shock interpretation for the nebula. Instead, we suggest that the binary interaction happened recently, thereby creating the nebula, with a cavity blown by the current fast stellar wind. This inference is supported by our results of 2D numerical hydrodynamic modelling.

We report the discovery of filamentary structures stretched behind the bow-shock-producing high-mass X-ray binary Vela X-1 using the SuperCOSMOS H-alpha Survey and present the results of optical spectroscopy of the bow shock carried out with the Southern African Large Telescope (SALT). The geometry of the detected structures suggests that Vela X-1 has encountered a wedge-like layer of enhanced density on its way and that the shocked material of the layer partially outlines a wake downstream of Vela X-1. To substantiate this suggestion, we carried out 3D magnetohydrodynamic simulations of interaction between Vela X-1 and the layer for three limiting cases. Namely, we run simulations in which (i) the stellar wind and the interstellar medium (ISM) were treated as pure hydrodynamic flows, (ii) a homogeneous magnetic field was added to the ISM, while the stellar wind was assumed to be unmagnetized, and (iii) the stellar wind was assumed to possess a helical magnetic field, while there was no magnetic field in the ISM. We found that although the first two simulations can provide a rough agreement with the observations, only the third one allowed us to reproduce not only the wake behind Vela X-1, but also the general geometry of the bow shock ahead of it.

We report the results of optical spectroscopy of the blue supergiant Sk-69 279 and its circular shell in the Large Magellanic Cloud (LMC) with the Southern African Large Telescope (SALT). We classify Sk-69 279 as an O9.2 Iaf star and analyse its spectrum by using the stellar atmosphere code CMFGEN, obtaining a stellar temperature of \approx 30 kK, a luminosity of log(L_*/Lsun)=5.54, a mass-loss rate of log(\dot{M}/(Msun/yr))=-5.26, and a wind velocity of 800 km/s. We found also that Sk-69 279 possesses an extended atmosphere with an effective temperature of \approx 24 kK and that its surface helium and nitrogen abundances are enhanced, respectively, by factors of \approx 2 and 20--30. This suggests that either Sk-69 279 was initially a (single) fast-rotating (\ga 400 km/s) star, which only recently evolved off the main sequence, or that it is a product of close binary evolution. The long-slit spectroscopy of the shell around Sk-69 279 revealed that its nitrogen abundance is enhanced by the same factor as the stellar atmosphere, which implies that the shell is composed mostly of the CNO processed material lost by the star. Our findings support previous propositions that some massive stars can produce compact circumstellar shells and, presumably, appear as luminous blue variables while they are still on the main sequence or have only recently left it.

High-resolution mid-infrared observations carried out by the Spitzer Space Telescope allowed one to resolve the fine structure of many astrospheres. In particular, they showed that the astrosphere around the B0.7 Ia star kappa Cas (HD 2905) has a clear-cut arc structure with numerous cirrus-like filaments beyond it. Previously, we suggested a physical mechanism for the formation of such filamentary structures. Namely, we showed theoretically that they might represent the non-monotonic spatial distribution of the interstellar dust in astrospheres (viewed as filaments) caused by interaction of the dust grains with the interstellar magnetic field disturbed in the astrosphere due to colliding of the stellar and interstellar winds. In this paper, we invoke this mechanism to explain the structure of the astrosphere around kappa Cas. We performed 3D magnetohydrodynamic modelling of the astrosphere for realistic parameters of the stellar wind and space velocity. The dust dynamics and the density distribution in the astrosphere were calculated in the framework of a kinetic model. It is found that the model results with the classical MRN size distribution of dust in the interstellar medium do not match the observations, and that the observed filamentary structure of the astrosphere can be reproduced only if the dust is composed mainly of big (micron-sized) grains. Comparison of the model results with observations allowed us to estimate parameters (number density and magnetic field strength) of the surrounding interstellar medium.

When a massive star in a binary system explodes as a supernova its companion star may be polluted with heavy elements from the supernova ejecta. Such a pollution had been detected in a handful of post-supernova binaries (Gonzalez Hernandez et al. 2011), but none of them is associated with a supernova remnant. We report the discovery of a solar-type star in a close, eccentric binary system with a neutron star within the young Galactic supernova remnant RCW 86. Our discovery implies that the supernova progenitor was a moving star, which exploded near the edge of its wind bubble and lost most of its initial mass due to common-envelope evolution shortly before core collapse. We find that the solar-type star is strongly polluted with calcium and other elements, which places the explosion within the class of calcium-rich supernovae -- faint and fast transients (Filippenko et al. 2003; Kasliwal et al. 2012}, whose origin is strongly debated (Kawabata et al. 2010; Waldman et al. 2011), and provides the first observational evidence that supernovae of this type can arise from core-collapse explosions (Kawabata et al. 2010; Moriya et al. 2010).

High-resolution mid-infrared observations of astrospheres show that many of them have filamentary (cirrus-like) structure. Using numerical models of dust dynamics in astrospheres, we suggest that their filamentary structure might be related to specific spatial distribution of the interstellar dust around the stars, caused by a gyrorotation of charged dust grains in the interstellar magnetic filed. Our numerical model describes the dust dynamics in astrospheres under an influence of the Lorentz force and assumption of a constant dust charge. Calculations are performed for the dust grains with different sizes separately. It is shown that non-monotonic spatial dust distribution (viewed as filaments) appears for dust grains with the period of gyromotion comparable with the characteristic time scale of the dust motion in the astrosphere. Numerical modelling demonstrates that number of filaments depends on charge-to-mass ratio of dust.

We report the results of spectroscopic observations and numerical modelling of the H II region IRAS 18153-1651. Our study was motivated by the discovery of an optical arc and two main-sequence stars of spectral type B1 and B3 near the centre of IRAS 18153-1651. We interpret the arc as the edge of the wind bubble (blown by the B1 star), whose brightness is enhanced by the interaction with a photoevaporation flow from a nearby molecular cloud. This interpretation implies that we deal with a unique case of a young massive star (the most massive member of a recently formed low-mass star cluster) caught just tens of thousands of years after its stellar wind has begun to blow a bubble into the surrounding dense medium. Our two-dimensional, radiation-hydrodynamics simulations of the wind bubble and the H II region around the B1 star provide a reasonable match to observations, both in terms of morphology and absolute brightness of the optical and mid-infrared emission, and verify the young age of IRAS 18153-1651. Taken together our results strongly suggest that we have revealed the first example of a wind bubble blown by a main-sequence B star.

Long-slit spectroscopy with the Southern African Large Telescope (SALT) of central stars of mid-infrared nebulae detected with the Spitzer Space Telescope and Wide-field Infrared Survey Explorer (WISE) led to the discovery of numerous candidate luminous blue variables (cLBVs) and other rare evolved massive stars. With the recent advent of the SALT fibre-fed high-resolution echelle spectrograph (HRS), a new perspective for the study of these interesting objects is appeared. Using the HRS we obtained spectra of a dozen newly identified massive stars. Some results on the recently identified cLBV Hen 3-729 are presented.

Searches for compact mid-IR nebulae with the Spitzer Space Telescope and the Wide-field Infrared Survey Explorer (WISE), accompanied by spectroscopic observations of central stars of these nebulae led to the discovery of many dozens of massive stars at different evolutionary stages, of which the most numerous are candidate luminous blue variables (LBVs). In this paper, we give a census of candidate and confirmed Galactic LBVs revealed with Spitzer and WISE, and present some new results of spectroscopic observations of central stars of mid-IR nebulae.

In this paper, we report the results of spectroscopic and photometric observations of the candidate evolved massive star MN48 disclosed via detection of a mid-infrared circular shell around it with the Spitzer Space Telescope. Follow-up optical spectroscopy of MN48 with the Southern African Large Telescope (SALT) carried out in 2011--2015 revealed significant changes in the spectrum of this star, which are typical of luminous blue variables (LBVs). The LBV status of MN48 was further supported by photometric monitoring which shows that in 2009--2011 this star has brightened by approx 0.9 and 1 mag in the V and I_c bands, respectively, then faded by approx 1.1 and 1.6 mag during the next four years, and apparently started to brighten again recently. The detected changes in the spectrum and brightness of MN48 make this star the 18th known Galactic bona fide LBV and increase the percentage of LBVs associated with circumstellar nebulae to more than 70 per cent. We discuss the possible birth place of MN48 and suggest that this star might have been ejected either from a putative star cluster embedded in the HII region IRAS 16455-4531 or the young massive star cluster Westerlund 1.

We present the results of optical spectroscopic observations of 54 candidate evolved massive stars revealed through the detection of mid-infrared nebulae of various shapes surrounding them with the {\it Spitzer Space Telescope} and {\it Wide-field Infrared Survey Explorer}. These observations, carried out with the Southern African Large Telescope (SALT) in 2010-2015, led to the discovery of about two dozens emission-line stars, of which 15 stars we classify as candidate luminous blue variables (cLBVs). Spectroscopic and photometric monitoring revealed significant changes in the spectra and brightness of four newly identified cLBVs, meaning that they are new members of the class of bona fide LBVs. We present an updated list of the Galactic bona fide LBVs. Currently, this list contains eighteen stars, of which more than 70 per cent are associated with circumstellar nebulae. We also discovered a very rare [WN] star - the central star of the planetary nebula Abell 48, and a WN3 star in a close, eccentric binary system with an O6V star in the Large Magellanic Cloud - the first-ever extragalactic massive star identified via detection of a circular shell around it. Most of the remaining targets are tentatively classified as OB, A and M stars.

Mid-infrared arcs of dust emission are often seen near ionizing stars within HII regions. A possible explanations for these arcs is that they could show the outer edges of asymmetric stellar wind bubbles. We use two-dimensional, radiation-hydrodynamics simulations of wind bubbles within HII regions around individual stars to predict the infrared emission properties of the dust within the HII region. We assume that dust and gas are dynamically well-coupled and that dust properties (composition, size distribution) are the same in the HII region as outside it, and that the wind bubble contains no dust. We post-process the simulations to make synthetic intensity maps at infrared wavebands using the TORUS code. We find that the outer edge of a wind bubble emits brightly at 24um through starlight absorbed by dust grains and re-radiated thermally in the infrared. This produces a bright arc of emission for slowly moving stars that have asymmetric wind bubbles, even for cases where there is no bow shock or any corresponding feature in tracers of gas emission. The 24um intensity decreases exponentially from the arc with increasing distance from the star because the dust temperature decreases with distance. The size distribution and composition of the dust grains has quantitative but not qualitative effects on our results. Despite the simplifications of our model, we find good qualitative agreement with observations of the HII region RCW120, and can provide physical explanations for any quantitative differences. Our model produces an infrared arc with the same shape and size as the arc around CD -38 11636 in RCW120, and with comparable brightness. This suggests that infrared arcs around O stars in HII regions may be revealing the extent of stellar wind bubbles, although we have not excluded other explanations.

We report the results of optical spectroscopy of the candidate evolved massive star MN44 revealed via detection of a circular shell with the Spitzer Space Telescope. First spectra taken in 2009 May--June showed the Balmer lines in emission as well as numerous emission lines of iron, which is typical of luminous blue variables (LBVs) near the visual maximum. New observations carried out in 2015 May--September detected significant changes in the spectrum, indicating that the star became hotter. We found that these changes are accompanied by significant brightness variability of MN44. In particular, the I_c-band brightness decreased by \approx 1.6 mag during the last six years and after reaching its minimum in 2015 June has started to increase. Using archival data, we also found that the I_c-band brightness increased by \approx 3 mag in \approx 30 yr preceding our observations. MN44 therefore represents the seventeenth known example of the Galactic bona fide LBVs. We detected a nitrogen-rich knot to the northwest of the star, which might represent an interstellar cloudlet interacting with the circumstellar shell. We discuss a possible association between MN44 and the INTEGRAL transient source of hard X-ray emission IGR J16327-4940, implying that MN44 might be either a colliding-wind binary or a high-mass X-ray binary.

We report the results of spectrophotometric observations of the massive star MN18 revealed via discovery of a bipolar nebula around it with the Spitzer Space Telescope. Using the optical spectrum obtained with the Southern African Large Telescope, we classify this star as B1 Ia. The evolved status of MN18 is supported by the detection of nitrogen overabundance in the nebula, which implies that it is composed of processed material ejected by the star. We analysed the spectrum of MN18 by using the code CMFGEN, obtaining a stellar effective temperature of \approx 21 kK. The star is highly reddened, E(B-V)\approx 2 mag. Adopting an absolute visual magnitude of M_V=-6.8\pm0.5 (typical of B1 supergiants), MN18 has a luminosity of log L/Lsun \approx 5.42\pm0.30, a mass-loss rate of \approx (2.8-4.5)\times10^{-7} Msun/yr, and resides at a distance of \approx 5.6^{+1.5} _{-1.2} kpc. We discuss the origin of the nebula around MN18 and compare it with similar nebulae produced by other blue supergiants in the Galaxy (Sher 25, HD 168625, [SBW2007] 1) and the Large Magellanic Cloud (Sk-69 202). The nitrogen abundances in these nebulae imply that blue supergiants can produce them from the main sequence stage up to the pre-supernova stage. We also present a K-band spectrum of the candidate luminous blue variable MN56 (encircled by a ring-like nebula) and report the discovery of an OB star at \approx 17 arcsec from MN18. The possible membership of MN18 and the OB star of the star cluster Lynga 3 is discussed.

In this Letter, we report the results of spectroscopic and photometric monitoring of the candidate luminous blue variable (LBV) WS1, which was discovered in 2011 through the detection of a mid-infrared circular shell and follow-up optical spectroscopy of its central star. Our monitoring showed that WS1 brightened in the B, V and I bands by more than 1 mag during the last three years, while its spectrum revealed dramatic changes during the same time period, indicating that the star became much cooler. The light curve of WS1 demonstrates that the brightness of this star has reached maximum in 2013 December and then starts to decline. These findings unambiguously proved the LBV nature of WS1 and added one more member to the class of Galactic bona fide LBVs, bringing their number to sixteen (an updated census of these objects is provided).

Interstellar bubbles around O stars are driven by a combination of the star's wind and ionizing radiation output. The wind contribution is uncertain because the boundary between the wind and interstellar medium is difficult to observe. Mid-infrared observations (e.g., of the H II region RCW 120) show arcs of dust emission around O stars, contained well within the H II region bubble. These arcs could indicate the edge of an asymmetric stellar wind bubble, distorted by density gradients and/or stellar motion. We present two-dimensional, radiation-hydrodynamics simulations investigating the evolution of wind bubbles and H II regions around massive stars moving through a dense (n=3000 cm^{-3}), uniform medium with velocities ranging from 4 to 16 km/s. The H II region morphology is strongly affected by stellar motion, as expected, but the wind bubble is also very aspherical from birth, even for the lowest space velocity considered. Wind bubbles do not fill their H II regions (we find filling factors of 10-20%), at least for a main sequence star with mass M~30 Msun. Furthermore, even for supersonic velocities the wind bow shock does not significantly trap the ionization front. X-ray emission from the wind bubble is soft, faint, and comes mainly from the turbulent mixing layer between the wind bubble and the H II region. The wind bubble radiates <1 per cent of its energy in X-rays; it loses most of its energy by turbulent mixing with cooler photoionized gas. Comparison of the simulations with the H II region RCW 120 shows that its dynamical age is <=0.4 Myr and that stellar motion <=4 km/s is allowed, implying that the ionizing source is unlikely to be a runaway star but more likely formed in situ. The region's youth, and apparent isolation from other O or B stars, makes it very interesting for studies of massive star formation and of initial mass functions.

We report the discovery of a circular mid-infrared shell around the emission-line star Wray 16-137 using archival data of the Spitzer Space Telescope. Follow-up optical spectroscopy of Wray 16-137 with the Southern African Large Telescope revealed a rich emission spectrum typical of the classical luminous blue variables (LBVs) like P Cygni. Subsequent spectroscopic and photometric observations showed drastic changes in the spectrum and brightness during the last three years, meaning that Wray 16-137 currently undergoes an S Dor-like outburst. Namely, we found that the star has brightened by \approx 1 mag in the V and I_c bands, while its spectrum became dominated by Fe ii lines. Taken together, our observations unambiguously show that Wray 16-137 is a new member of the family of Galactic bona fide LBVs.

At least 5 per cent of the massive stars are moving supersonically through the interstellar medium (ISM) and are expected to produce a stellar wind bow shock. We explore how the mass loss and space velocity of massive runaway stars affect the morphology of their bow shocks. We run two-dimensional axisymmetric hydrodynamical simulations following the evolution of the circumstellar medium of these stars in the Galactic plane from the main sequence to the red supergiant phase. We find that thermal conduction is an important process governing the shape, size and structure of the bow shocks around hot stars, and that they have an optical luminosity mainly produced by forbidden lines, e.g. [OIII]. The Ha emission of the bow shocks around hot stars originates from near their contact discontinuity. The H$\alpha$ emission of bow shocks around cool stars originates from their forward shock, and is too faint to be observed for the bow shocks that we simulate. The emission of optically-thin radiation mainly comes from the shocked ISM material. All bow shock models are brighter in the infrared, i.e. the infrared is the most appropriate waveband to search for bow shocks. Our study suggests that the infrared emission comes from near the contact discontinuity for bow shocks of hot stars and from the inner region of shocked wind for bow shocks around cool stars. We predict that, in the Galactic plane, the brightest, i.e. the most easily detectable bow shocks are produced by high-mass stars moving with small space velocities.

Betelgeuse, a nearby red supergiant, is a runaway star with a powerful stellar wind that drives a bow shock into its surroundings. This picture has been challenged by the discovery of a dense and almost static shell that is three times closer to the star than the bow shock and has been decelerated by some external force. The two physically distinct structures cannot both be formed by the hydrodynamic interaction of the wind with the interstellar medium. Here we report that a model in which Betelgeuse's wind is photoionized by radiation from external sources can explain the static shell without requiring a new understanding of the bow shock. Pressure from the photoionized wind generates a standing shock in the neutral part of the wind and forms an almost static, photoionization-confined shell. Other red supergiants should have significantly more massive shells than Betelgeuse, because the photoionization-confined shell traps up to 35 per cent of all mass lost during the red supergiant phase, confining this gas close to the star until it explodes. After the supernova explosion, massive shells dramatically affect the supernova lightcurve, providing a natural explanation for the many supernovae that have signatures of circumstellar interaction.

We discuss two important effects for the astrospheres of runaway stars: the propagation of ionizing photons far beyond the astropause, and the rapid evolution of massive stars (and their winds) near the end of their lives. Hot stars emit ionizing photons with associated photoheating that has a significant dynamical effect on their surroundings. 3D simulations show that HII regions around runaway O stars drive expanding conical shells and leave underdense wakes in the medium they pass through. For late O stars this feedback to the interstellar medium is more important than that from stellar winds. Late in life, O stars evolve to cool red supergiants more rapidly than their environment can react, producing transient circumstellar structures such as double bow shocks. This provides an explanation for the bow shock and linear bar-shaped structure observed around Betelgeuse.

The circumstellar medium around massive stars is strongly impacted by stellar winds, radiation, and explosions. We use numerical simulations of these interactions to constrain the current properties and evolutionary history of various stars by comparison with observed circumstellar structures. Two- and three-dimensional simulations of bow shocks around red supergiant stars have shown that Betelgeuse has probably only recently evolved from a blue supergiant to a red supergiant, and hence its bow shock is very young and has not yet reached a steady state. We have also for the first time investigated the magnetohydrodynamics of the photoionised H II region around the nearby runaway O star Zeta Oph. Finally, we have calculated a grid of models of bow shocks around main sequence and evolved massive stars that has general application to many observed bow shocks, and which forms the basis of future work to model the explosions of these stars into their pre-shaped circumstellar medium.

We report the first-ever discovery of a Wolf-Rayet (WR) star in the Large Magellanic Cloud via detection of a circular shell with the Spitzer Space Telescope. Follow-up observations with Gemini-South resolved the central star of the shell into two components separated from each other by approx 2 arcsec (or approx 0.5 pc in projection). One of these components turns out to be a WN3 star with H and He lines both in emission and absorption (we named it BAT99 3a using the numbering system based on extending the Breysacher et al. catalogue). Spectroscopy of the second component showed that it is a B0 V star. Subsequent spectroscopic observations of BAT99 3a with the du Pont 2.5-m telescope and the Southern African Large Telescope revealed that it is a close, eccentric binary system, and that the absorption lines are associated with an O companion star. We analyzed the spectrum of the binary system using the non-LTE Potsdam Wolf-Rayet (PoWR) code, confirming that the WR component is a very hot (approx 90 kK) WN star. For this star, we derived a luminosity of log L/Lsun =5.45 and a mass-loss rate of 10^{-5.8} Msun/yr, and found that the stellar wind composition is dominated by helium with 20 per cent of hydrogen. Spectroscopy of the shell revealed an He iii region centred on BAT99 3a and having the same angular radius (approx 15 arcsec) as the shell. We thereby add a new example to a rare class of high-excitation nebulae photoionized by WR stars. Analysis of the nebular spectrum showed that the shell is composed of unprocessed material, implying that the shell was swept-up from the local interstellar medium. We discuss the physical relationship between the newly identified massive stars and their possible membership of a previously unrecognized star cluster.

In this Letter, we explore the hypothesis that the smooth appearance of bow shocks around some red supergiants (RSGs) might be caused by the ionization of their winds by external sources of radiation. Our numerical simulations of the bow shock generated by IRC-10414 (the first-ever RSG with an optically detected bow shock) show that the ionization of the wind results in its acceleration by a factor of two, which reduces the difference between the wind and space velocities of the star and makes the contact discontinuity of the bow shock stable for a range of stellar space velocities and mass-loss rates. Our best fit model reproduces the overall shape and surface brightness of the observed bow shock and suggests that the space velocity and mass-loss rate of IRC-10414 are $\approx$50 ${\rm km} \, {\rm s}^{-1}$ and $\approx$$10^{-6}$ $M_\odot \, {\rm yr}^{-1}$, respectively, and that the number density of the local ISM is $\approx$3 ${\rm cm}^{-3}$. It also shows that the bow shock emission comes mainly from the shocked stellar wind. This naturally explains the enhanced nitrogen abundance in the line-emitting material, derived from the spectroscopy of the bow shock. We found that photoionized bow shocks are $\approx$15$-$50 times brighter in optical line emission than their neutral counterparts, from which we conclude that the bow shock of IRC-10414 must be photoionized.

Massive stars rapidly change their masses through strong stellar winds and mass transfer in binary systems. We show that such mass changes leave characteristic signatures in stellar mass functions of young star clusters which can be used to infer their ages and to identify products of binary evolution. We model the observed present day mass functions of the young Galactic Arches and Quintuplet star clusters using our rapid binary evolution code. We find that shaping of the mass function by stellar wind mass loss allows us to determine the cluster ages to 3.5$\pm$0.7 Myr and 4.8$\pm$1.1 Myr, respectively. Exploiting the effects of binary mass exchange on the cluster mass function, we find that the most massive stars in both clusters are rejuvenated products of binary mass transfer, i.e. the massive counterpart of classical blue straggler stars. This resolves the problem of an apparent age spread among the most luminous stars exceeding the expected duration of star formation in these clusters. We perform Monte Carlo simulations to probe stochastic sampling, which support the idea of the most massive stars being rejuvenated binary products. We find that the most massive star is expected to be a binary product after 1.0$\pm$0.7 Myr in Arches and after 1.7$\pm$1.0 Myr in Quintuplet. Today, the most massive 9$\pm$3 stars in Arches and 8$\pm$3 in Quintuplet are expected to be such objects. Our findings have strong implications for the stellar upper mass limit and solve the discrepancy between the claimed 150 $\mathrm{M}_\odot$ limit and observations of fours stars with initial masses of 165-320 $\mathrm{M}_\odot$ in R136 and of SN 2007bi, which is thought to be a pair-instability supernova from an initial 250 $\mathrm{M}_\odot$ star. Using the stellar population of R136, we revise the upper mass limit to values in the range 200-500 $\mathrm{M}_\odot$.

We report the results of optical spectroscopy of a candidate evolved massive star in the Cygnus X region, TYC 3159-6-1, revealed via detection of its curious circumstellar nebula in archival data of the Spitzer Space Telecope. We classify TYC 3159-6-1 as an O9.5-O9.7 Ib star and derive its fundamental parameters by using the stellar atmosphere code FASTWIND. The He and CNO abundances in the photosphere of TYC 3159-6-1 are consistent with the solar abundances, suggesting that the star only recently evolved off the main sequence. Proper motion and radial velocity measurements for TYC 3159-6-1 show that it is a runaway star. We propose that Dolidze 7 is its parent cluster. We discuss the origin of the nebula around TYC 3159-6-1 and suggest that it might be produced in several successive episodes of enhanced mass-loss rate (outbursts) caused by rotation of the star near the critical, \Omega-limit.

Most runaway OB stars, like the majority of massive stars residing in their parent clusters, go through the red supergiant (RSG) phase during their lifetimes. Nonetheless, although many dozens of massive runaways were found to be associated with bow shocks, only two RSG bow-shock-producing stars, Betelgeuse and \mu Cep, are known to date. In this paper, we report the discovery of an arc-like nebula around the late M-type star IRC-10414 using the SuperCOSMOS H-alpha Survey. Our spectroscopic follow-up of IRC-10414 with the Southern African Large Telescope (SALT) showed that it is a M7 supergiant, which supports previous claims on the RSG nature of this star based on observations of its maser emission. This was reinforced by our new radio- and (sub)millimeter-wavelength molecular line observations made with the Atacama Pathfinder Experiment (APEX) 12 meter telescope and the Effelsberg 100 m radio telescope, respectively. The SALT spectrum of the nebula indicates that its emission is the result of shock excitation. This finding along with the arc-like shape of the nebula and an estimate of the space velocity of IRC-10414 (\approx 70\pm20 km/s) imply the bow shock interpretation for the nebula. Thus, IRC-10414 represents the third case of a bow-shock-producing RSG and the first one with a bow shock visible at optical wavelengths. We discuss the smooth appearance of the bow shocks around IRC-10414 and Betelgeuse and propose that one of the necessary conditions for stability of bow shocks generated by RSGs is the ionization of the stellar wind. Possible ionisation sources of the wind of IRC-10414 are proposed and discussed.

At least 25 per cent of massive stars are ejected from their parent cluster, becoming runaways or exiles, travelling with often-supersonic space velocities through the interstellar medium (ISM). Their overpressurised H II regions impart kinetic energy and momentum to the ISM, compress and/or evaporate dense clouds, and can constrain properties of both the star and the ISM. Here we present one-, two-, and (the first) three-dimensional simulations of the H II region around a massive star moving supersonically through a uniform, magnetised ISM, with properties appropriate for the nearby O star Zeta Oph. The H II region leaves an expanding overdense shell behind the star and, inside this, an underdense wake that should be filled with hot gas from the shocked stellar wind. The gas column density in the shell is strongly influenced by the ISM magnetic field strength and orientation. H-alpha emission maps show the H II region remains roughly circular, although the star is displaced somewhat from the centre of emission. For our model parameters, the kinetic energy feedback from the H II region is comparable to the mechanical luminosity of the stellar wind, and the momentum feedback rate is >100X larger than that from the wind and about 10X larger than the total momentum input rate available from radiation pressure. Compared to the star's eventual supernova explosion, the kinetic energy feedback from the H II region over the star's main sequence lifetime is >100X less, but the momentum feedback is up to 4X larger. H II region dynamics are found to have only a small effect on the ISM conditions that a bow shock close to the star would encounter.

We describe the results of the world-wide observing campaign of the highly eccentric Be binary system delta Scorpii 2011 periastron passage which involved professional and amateur astronomers. Our spectroscopic observations provided a precise measurement of the system orbital period at 10.8092+/- 0.0005 years. Fitting of the He II 4686A line radial velocity curve determined the periastron passage time on 2011 July 3, UT 9:20 with a 0.9--day uncertainty. Both these results are in a very good agreement with recent findings from interferometry. We also derived new evolutionary masses of the binary components (13 and 8.2 Msun) and a new distance of 136 pc from the Sun, consistent with the HIPPARCOS parallax. The radial velocity and profile variations observed in the H_alpha line near the 2011 periastron reflected the interaction of the secondary component and the circumstellar disk around the primary component. Using these data, we estimated a disk radius of 150 Rsun. Our analysis of the radial velocity variations measured during the periastron passage time in 2000 and 2011 along with those measured during the 20th century, the high eccentricity of the system, and the presence of a bow shock-like structure around it suggest that delta Sco might be a runaway triple system. The third component should be external to the known binary and move on an elliptical orbit that is tilted by at least 40 degree with respect to the binary orbital plane for such a system to be stable and responsible for the observed long-term radial velocity variations.

A considerable fraction of the central stars of planetary nebulae (CSPNe) are hydrogen-deficient. Almost all of these H-deficient central stars (CSs) display spectra with strong carbon and helium lines. Most of them exhibit emission line spectra resembling those of massive WC stars. Therefore these stars are classed as CSPNe of spectral type [WC]. Recently, quantitative spectral analysis of two emission-line CSs, PB 8 and IC 4663, revealed that these stars do not belong to the [WC] class. Instead PB 8 has been classified as [WN/WC] type and IC 4663 as [WN] type. In this work we report the spectroscopic identification of another rare [WN] star, the CS of Abell 48. We performed a spectral analysis of Abell 48 with the Potsdam Wolf-Rayet (PoWR) models for expanding atmospheres. We find that the expanding atmosphere of Abell 48 is mainly composed of helium (85 per cent by mass), hydrogen (10 per cent), and nitrogen (5 per cent). The residual hydrogen and the enhanced nitrogen abundance make this object different from the other [WN] star IC 4663. We discuss the possible origin of this atmospheric composition.

Two optically obscured Wolf-Rayet (WR) stars have been recently discovered by means of their infrared (IR) circumstellar shells, which show signatures of interaction with each other. Following the systematics of the WR star catalogues, these stars obtain the names WR\,120bb and WR\,120bc. In this paper, we present and analyse new near-IR, $J$, $H$, and $K$-band, spectra using the Potsdam Wolf-Rayet (PoWR) model atmosphere code. For that purpose, the atomic data base of the code has been extended in order to include all significant lines in the near-IR bands. The spectra of both stars are classified as WN9h. As their spectra are very similar the parameters that we obtained by the spectral analyses hardly differ. Despite their late spectral subtype, we found relatively high stellar temperatures of 63 kK. The wind composition is dominated by helium, while hydrogen is depleted to 25 per cent by mass. Because of their location in the Scutum-Centaurus arm, WR\,120bb and WR\,120bc appear highly reddened, $A_{K_{\rm s}} \approx 2$ mag. We adopt a common distance of 5.8\,kpc to both stars, which complies with the typical absolute $K$-band magnitude for the WN9h subtype of -6.5 mag, is consistent with their observed extinction based on comparison with other massive stars in the region, and allows for the possibility that their shells are interacting with each other. This leads to luminosities of $\log(L/L_\odot) = 5.66$ and 5.54 for WR\,120bb and WR\,120bc, with large uncertainties due to the adopted distance. The values of the luminosities of WR\,120bb and WR\,120bc imply that the immediate precursors of both stars were red supergiants (RSG). This implies in turn that the circumstellar shells associated with WR\,120bb and WR\,120bc were formed by interaction between the WR wind and the dense material shed during the preceding RSG phase.

We report the discovery of a bow-shock-producing star in the vicinity of the young massive star cluster NGC 3603 using archival data of the Spitzer Space Telescope. Follow-up optical spectroscopy of this star with Gemini-South led to its classification as O6 V. The orientation of the bow shock and the distance to the star (based on its spectral type) suggest that the star was expelled from the cluster, while the young age of the cluster (~2 Myr) implies that the ejection was caused by a dynamical few-body encounter in the cluster's core. The relative position on the sky of the O6 V star and a recently discovered O2 If*/WN6 star (located on the opposite side of NGC 3603) allows us to propose that both objects were ejected from the cluster via the same dynamical event -- a three-body encounter between a single (O6 V) star and a massive binary (now the O2 If*/WN6 star). If our proposal is correct, then one can "weigh" the O2 If*/WN6 star using the conservation of the linear momentum. Given a mass of the O6 V star of 30 Msun, we found that at the moment of ejection the mass of the O2 If*/WN6 star was 175 Msun. Moreover, the observed X-ray luminosity of the O2 If*/WN6 star (typical of a single star) suggests that the components of this originally binary system have merged (e.g., because of encounter hardening).

Mass-loss rate, $\dot{M}$, is one of the key parameters affecting evolution and observational manifestations of massive stars, and their impact on the ambient medium. Despite its importance, there is a factor of ~100 discrepancy between empirical and theoretical $\dot{M}$ of late-type O dwarfs, the so-called weak-wind problem. In this Letter, we propose a simple novel method to constrain $\dot{M}$ of runaway massive stars through observation of their bow shocks and Stromgren spheres, which might be of decisive importance for resolving the weak-wind problem. Using this method, we found that $\dot{M}$ of the well-known runaway O9.5 V star zeta Oph is more than an order of magnitude higher than that derived from ultraviolet (UV) line-fitting (Marcolino et al. 2009) and is by a factor of 6 to 7 lower than those based on the theoretical recipe by Vink et al. (2000) and the Halpha line (Mokiem et al. 2005). The discrepancy between $\dot{M}$ derived by our method and that based on UV lines would be even more severe if the stellar wind is clumpy. At the same time, our estimate of $\dot{M}$ agrees with that predicted by the moving reversing layer theory by Lucy (2010).

A significant fraction of massive stars in the Milky Way and other galaxies are located far from star clusters. It is known that some of these stars are runaways and therefore most likely were formed in embedded clusters and then ejected into the field because of dynamical few-body interactions or binary-supernova explosions. However, there exists a group of field O stars whose runaway status is difficult to prove via direct proper motion measurements or whose low space velocities and/or young ages appear to be incompatible with their large separation from known star clusters. The existence of this group led some authors to believe that field O stars can form in situ. In this paper, we examine the runaway status of the best candidates for isolated formation of massive stars in the Milky Way and the Magellanic Clouds by searching for bow shocks around them, by using the new reduction of the Hipparcos data, and by searching for stellar systems from which they could originate within their lifetimes. We show that most of the known O stars thought to have formed in isolation are instead very likely runaways. We show also that the field must contain a population of O stars whose low space velocities and/or young ages are in apparent contradiction with the large separation of these stars from their parent clusters and/or the ages of these clusters. These stars (the descendants of runaway massive binaries) cannot be traced back to their parent clusters and therefore can be mistakenly considered as having formed in situ. We argue also that some field O stars could be detected in optical wavelengths only because they are runaways, while their cousins residing in the deeply embedded parent clusters might still remain totally obscured. The main conclusion of our study is that there is no significant evidence whatsoever in support of the in situ proposal on the origin of massive stars.

We report the discovery of a bipolar nebula around the peculiar emission-line star MWC 349A using archival Spitzer Space Telescope 24 um data. The nebula extends over several arcminutes (up to 5 pc) and has the same orientation and geometry as the well-known subarcsecond-scale (~400 times smaller) bipolar radio nebula associated with this star. We discuss the physical relationship between MWC 349A and the nearby B0 III star MWC 349B and propose that both stars were members of a hierarchical triple system, which was ejected from the core of the Cyg OB2 association several Myr ago and recently was dissolved into a binary system (now MWC 349A) and a single unbound star (MWC 349B). Our proposal implies that MWC 349A is an evolved massive star (likely a luminous blue variable) in a binary system with a low-mass star. A possible origin of the bipolar nebula around MWC 349A is discussed.

We report the discovery of two new Galactic candidate luminous blue variable (cLBV) stars via detection of circular shells (typical of known confirmed and cLBVs) and follow-up spectroscopy of their central stars. The shells were detected at 22 um in the archival data of the Mid-Infrared All Sky Survey carried out with the Wide-field Infrared Survey Explorer (WISE). Follow-up optical spectroscopy of the central stars of the shells conducted with the renewed Southern African Large Telescope (SALT) showed that their spectra are very similar to those of the well-known LBVs P Cygni and AG Car, and the recently discovered cLBV MN112, which implies the LBV classification for these stars as well. The LBV classification of both stars is supported by detection of their significant photometric variability: one of them brightened in the R- and I-bands by 0.68\pm0.10 mag and 0.61\pm0.04 mag, respectively, during the last 13-18 years, while the second one (known as Hen 3-1383) varies its B,V,R,I and K_s brightnesses by \simeq 0.5-0.9 mag on time-scales from 10 days to decades. We also found significant changes in the spectrum of Hen 3-1383 on a timescale of \simeq 3 months, which provides additional support for the LBV classification of this star. Further spectrophotometric monitoring of both stars is required to firmly prove their LBV status. We discuss a connection between the location of massive stars in the field and their fast rotation, and suggest that the LBV activity of the newly discovered cLBVs might be directly related to their possible runaway status.

We have undertaken a near-IR spectral survey of stars associated with compact nebulae recently revealed by the Spitzer and WISE imaging surveys. These circumstellar nebulae, produced by massive evolved stars, display a variety of symmetries and shapes and are often only evident at mid-IR wavelengths. Stars associated with ~50 of these nebulae have been observed. We also obtained recent spectra of previously confirmed (known) luminous blue variables (LBVs) and candidate LBVs (cLBVs). The spectral similarity of the stars observed when compared directly to known LBVs and Wolf-Rayet (WR) stars indicate many are newly identified cLBVs, with a few being newly discovered WR stars, mostly of WN8-9h spectral type. These results suggest that a large population of previously unidentified cLBVs and related transitional stars reside in the Galaxy and confirm that circumstellar nebulae are inherent to most (c)LBVs.

We have undertaken a near-infrared spectral survey of stars associated with compact mid-IR shells recently revealed by the MIPSGAL (24 micron) and GLIMPSE (8 micron) Spitzer surveys, whose morphologies are typical of circumstellar shells produced by massive evolved stars. Through spectral similarity with known Luminous Blue Variable (LBV) and Wolf-Rayet (WR) stars, a large population of candidate LBVs (cLBVs) and a smaller number of new WR stars are being discovered. This significantly increases the Galactic cLBV population and confirms that nebulae are inherent to most (if not all) objects of this class. Keywords - stars: emission-line, Be, stars: mass loss, stars: winds, outflows, stars: Wolf-Rayet

The question whether or not massive stars can form in isolation or only in star clusters is of great importance for the theory of (massive) star-formation as well as for the stellar initial mass function of whole galaxies (IGIMF-theory). While a seemingly easy question it is rather difficult to answer. Several physical processes (e.g. star-loss due to stellar dynamics or gas expulsion) and observational limitations (e.g. dust obscuration of young clusters, resolution) pose severe challenges to answer this question. In this contribution we will present the current arguments in favour and against the idea that all O stars form in clusters.

We report the first-ever discovery of an extragalactic Wolf-Rayet (WR)star with Spitzer. A new WR star in the Large Magellanic Cloud (LMC) was revealed via detection of its circumstellar shell using 24 {\mu}m images obtained in the framework of the Spitzer Survey of the Large Magellanic Cloud (SAGE-LMC). Subsequent spectroscopic bservations with the Gemini South resolved the central star in two components, one of which is a WN3b+abs star, while the second one is a B0V star. We consider the lopsided brightness distribution over the circumstellar shell as an indication that the WR star is a runaway and use this interpretation to identify a possible parent cluster of the star.

Dynamical few-body encounters in the dense cores of young massive star clusters are responsible for the loss of a significant fraction of their massive stellar content. Some of the escaping (runaway) stars move through the ambient medium supersonically and can be revealed via detection of their bow shocks (visible in the infrared, optical or radio). In this paper, which is the second of a series of papers devoted to the search for OB stars running away from young (several Myr) Galactic clusters and OB associations, we present the results of the search for bow shocks around the star-forming region NGC 6357. Using the archival data of the Midcourse Space Experiment (MSX) satellite and the Spitzer Space Telescope, and the preliminary data release of the Wide-Field Infrared Survey Explorer (WISE), we discovered seven bow shocks, whose geometry is consistent with the possibility that they are generated by stars expelled from the young star clusters, Pismis 24 and AH03 J1725-34.4, associated with NGC 6357. Two of the seven bow shocks are driven by the already known O stars. Follow-up spectroscopy of three other bow shock-producing stars showed that they are O-type stars as well, while the 2MASS photometry of the remaining two stars suggests that they could be B0 V stars, provided that both are located at the same distance as NGC 6357. Detection of numerous massive stars ejected from the very young clusters is consistent with the theoretical expectation that star clusters can effectively lose massive stars at the very beginning of their dynamical evolution and lends strong support to the idea that probably all field OB stars have been dynamically ejected from their birth clusters. A by-product of our search for bow shocks around NGC 6357 is the detection of three circular shells typical of luminous blue variable and late WN-type Wolf-Rayet stars.

We report the discovery of a bow shock around the high-mass X-ray binary (HMXB) 4U 1907+09 using the Spitzer Space Telescope 24 $\mu$m data (after Vela X-1 the second example of bow shocks associated with HMXBs). The detection of the bow shock implies that 4U 1907+09 is moving through the space with a high (supersonic) peculiar velocity. To confirm the runaway nature of 4U 1907+09, we measured its proper motion, which for an adopted distance to the system of 4 kpc corresponds to a peculiar transverse velocity of $\simeq 160 \pm 115$ km/s, meaning that 4U 1907+09 is indeed a runaway system and supporting the general belief that most of HMXBs possess high space velocities. The direction of motion of 4U 1907+09 inferred from the proper motion measurement is consistent with the orientation of the symmetry axis of the bow shock, and shows that the HMXB is running away from the Galactic plane. We also present the Spitzer images of the bow shock around Vela X-1 (a system similar to 4U 1907+09) and compare it with the bow shock generated by 4U 1907+09.

Using archival Spitzer Space Telescope data, we identified for the first time a dozen runaway OB stars in the Small Magellanic Cloud (SMC) through the detection of their bow shocks. The geometry of detected bow shocks allows us to infer the direction of motion of the associated stars and to determine their possible parent clusters and associations. One of the identified runaway stars, AzV 471, was already known as a high-velocity star on the basis of its high peculiar radial velocity, which is offset by ~40 km/s from the local systemic velocity. We discuss implications of our findings for the problem of the origin of field OB stars. Several of the bow shock-producing stars are found in the confines of associations, suggesting that these may be "alien" stars contributing to the age spread observed for some young stellar systems. We also report the discovery of a kidney-shaped nebula attached to the early WN-type star SMC-WR3 (AzV 60a). We interpreted this nebula as an interstellar structure created owing to the interaction between the stellar wind and the ambient interstellar medium.

Very massive stars preferentially reside in the cores of their parent clusters and form binary or multiple systems. We study the role of tight very massive binaries in the origin of the field population of very massive stars. We performed numerical simulations of dynamical encounters between single (massive) stars and a very massive binary with parameters similar to those of the most massive known Galactic binaries, WR 20a and NGC 3603-A1. We found that these three-body encounters could be responsible for the origin of high peculiar velocities ($\geq$ 70 km/s) observed for some very massive ($\geq$ 60-70 Msun) runaway stars in the Milky Way and the Large Magellanic Cloud (e.g., $\lambda$ Cep, BD+43 3654, Sk-67 22, BI 237, 30 Dor 016), which can hardly be explained within the framework of the binary-supernova scenario. The production of high-velocity massive stars via three-body encounters is accompanied by the recoil of the binary in the opposite direction to the ejected star. We show that the relative position of the very massive binary R145 and the runaway early B-type star Sk-69 206 on the sky is consistent with the possibility that both objects were ejected from the central cluster, R136, of the star-forming region 30 Doradus via the same dynamical event -- a three-body encounter.

The origin of massive field stars in the Large Magellanic Cloud (LMC) has long been an enigma. The recent measurements of large offsets (~100 km/s) between the heliocentric radial velocities of some very massive (O2-type) field stars and the systemic LMC velocity provides a possible explanation of this enigma and suggests that the field stars are runaway stars ejected from their birth places at the very beginning of their parent cluster's dynamical evolution. A straightforward way to prove this explanation is to measure the proper motions of the field stars and to show that they are moving away from one of the nearby star clusters or OB associations. This approach however is complicated by the large distance to the LMC, which makes accurate proper motion measurements difficult. We use an alternative approach for solving the problem, based on the search for bow shocks produced by runaway stars. The geometry of detected bow shocks would allow us to infer the direction of stellar motion and thereby to determine their possible parent clusters. In this paper we present the results of a search for bow shocks around six massive field stars which were suggested in the literature as candidate runaway stars. Using archival (Spitzer Space Telescope) data, we found a bow shock associated with one of our program stars, the O2 V((f*)) star BI 237, which is the first-ever detection of bow shocks in the LMC. Orientation of the bow shock suggests that BI 237 was ejected from the OB association LH 82 (located at ~120 pc in projection from the star). A by-product of our search is the detection of bow shocks generated by four OB stars in the field of the LMC and an arc-like structure attached to the candidate luminous blue variable R81 (HD 269128). The geometry of two of these bow shocks is consistent with the possibility that their associated stars were ejected from the 30 Doradus star forming complex.

Massive evolved stars loss a large fraction of their mass via copious stellar wind or instant outbursts and during certain evolutionary phases they can be identified via the presence of their circumstellar nebulae. In this paper, we present the results of search for compact nebulae (reminiscent of circumstellar nebulae around evolved massive stars) using archival 24 $\mu$m data obtained with the Multiband Imaging Photometer for Spitzer. We discovered 115 nebulae, most of which bear a striking resemblance to the circumstellar nebulae associated with Luminous Blue Variables (LBVs) and late WN-type (WNL) Wolf-Rayet (WR) stars in the Milky Way and the Large Magellanic Cloud (LMC). We interpret this similarity as an indication that the central stars of detected nebulae are either LBVs or related evolved massive stars. Our interpretation is supported by follow-up spectroscopy of two dozens of these central stars, most of which turns out to be either candidate LBVs (cLBVs), blue supergiants or WNL stars. We expect that the forthcoming spectroscopy of the remaining objects from our list, accompanied by the spectrophotometric monitoring of the already discovered cLBVs, will further increase the known population of Galactic LBVs, which in turn would have profound consequences for better understanding the LBV phenomenon and its role in the transition between hydrogen burning O stars and helium burning WR stars. We also report the detection of an arc-like structure attached to the cLBV HD326823 and an arc associated with the LBV R99 (HD269445) in the LMC.

We report the discovery of a new Galactic candidate Luminous Blue Variable (cLBV) via detection of an infrared circular nebula and follow-up spectroscopy of its central star. The nebula, MN112, is one of many dozens of circular nebulae detected at $24 \mu$m in the {\it Spitzer Space Telescope} archival data, whose morphology is similar to that of nebulae associated with known (c)LBVs and related evolved massive stars. Specifically, the core-halo morphology of MN112 bears a striking resemblance to the circumstellar nebula associated with the Galactic cLBV GAL 079.29+00.46, which suggests that both nebulae might have a similar origin and that the central star of MN112 is a LBV. The spectroscopy of the central star showed that its spectrum is almost identical to that of the bona fide LBV P Cygni, which also supports the LBV classification of the object. To further constrain the nature of MN112, we searched for signatures of possible high-amplitude ($\ga 1$ mag) photometric variability of the central star using archival and newly obtained photometric data covering a 45 year period. We found that the B magnitude of the star was constant ($\simeq$ 17.1$\pm$0.3 mag) over this period, while in the I band the star brightened by $\simeq 0.4$ mag during the last 17 years. Although the non-detection of large photometric variability leads us to use the prefix `candidate' in the classification of MN112, we remind that the long-term photometric stability is not unusual for genuine LBVs and that the brightness of P Cygni remains relatively stable during the last three centuries.

We report the discovery of a new Wolf-Rayet star in Aquila via detection of its circumstellar nebula (reminiscent of ring nebulae associated with late WN stars) using the Spitzer Space Telescope archival data. Our spectroscopic follow-up of the central point source associated with the nebula showed that it is a WN7h star (we named it WR121b). We analyzed the spectrum of WR121b by using the Potsdam Wolf-Rayet (PoWR) model atmospheres, obtaining a stellar temperature of ~ 50 kK. The stellar wind composition is dominated by helium with ~ 20 per cent of hydrogen. The stellar spectrum is highly reddened (E_{B-V} = 2.85 mag). Adopting an absolute magnitude of M_v = -5.7, the star has a luminosity of log L/Lsun = 5.75 and a mass-loss rate of 10^{-4.7} Msun/yr, and resides in a distance of 6.3 kpc. We searched for a possible parent cluster of WR121b and found that this star is located at ~ 1 degree from the young star cluster embedded in the giant HII region W43 (containing a WN7+a/OB? star -- WR121a). We also discovered a bow shock around the O9.5III star ALS9956, located at ~ 0.5 degree from the cluster. We discuss the possibility that WR121b and ALS9956 are runaway stars ejected from the cluster in W43.

We performed numerical simulations of dynamical encounters between hard massive binaries and a very massive star (VMS; formed through runaway mergers of ordinary stars in the dense core of a young massive star cluster), in order to explore the hypothesis that this dynamical process could be responsible for the origin of high-velocity (\geq 200-400 km/s) early or late B-type stars. We estimated the typical velocities produced in encounters between very tight massive binaries and VMSs (of mass of \geq 200 Msun) and found that about 3-4 per cent of all encounters produce velocities of \geq 400 km/s, while in about 2 per cent of encounters the escapers attain velocities exceeding the Milky Ways's escape velocity. We therefore argue that the origin of high-velocity (\geq 200-400 km/s) runaway stars and at least some so-called hypervelocity stars could be associated with dynamical encounters between the tightest massive binaries and VMSs formed in the cores of star clusters. We also simulated dynamical encounters between tight massive binaries and single ordinary 50-100 Msun stars. We found that from 1 to \simeq 4 per cent of these encounters can produce runaway stars with velocities of \geq 300-400 km/s (typical of the bound population of high-velocity halo B-type stars) and occasionally (in less than 1 per cent of encounters) produce hypervelocity (\geq 700 km/s) late B-type escapers.

We discuss the origin of the runaway early B-type star HD271791 and show that its extremely high velocity (\simeq 530-920 km/s) cannot be explained within the framework of the binary-supernova ejection scenario. Instead, we suggest that HD271791 attained its peculiar velocity in the course of a strong dynamical encounter between two hard massive binaries or via an exchange encounter between a hard massive binary and a very massive star, formed through runaway mergers of ordinary massive stars in the dense core of a young massive star cluster.

We report the serendipitous discovery of a ring nebula around a candidate Wolf-Rayet (WR) star, HBHA 4202-22, in Cygnus using the Spitzer Space Telescope archival data. Our spectroscopic follow-up observations confirmed the WR nature of this star (we named it WR138a) and showed that it belongs to the WN8-9h subtype. We thereby add a new example to the known sample of late WN stars with circumstellar nebulae. We analyzed the spectrum of WR138a by using the Potsdam Wolf-Rayet (PoWR) model atmospheres, obtaining a stellar temperature of 40 kK. The stellar wind composition is dominated by helium with 20 per cent of hydrogen. The stellar spectrum is highly reddened and absorbed (E_{B-V} = 2.4 mag, A_V = 7.4 mag). Adopting a stellar luminosity of log L/Lsun = 5.3, the star has a mass-loss rate of 10^{-4.7} Msun yr^{-1}, and resides in a distance of 4.2 kpc. We measured the proper motion for WR138a and found that it is a runaway star with a peculiar velocity of ~50 km/s. Implications of the runaway nature of WR138a for constraining the mass of its progenitor star and understanding the origin of its ring nebula are discussed.

We explore the hypothesis that some high-velocity runaway stars attain their peculiar velocities in the course of exchange encounters between hard massive binaries and a very massive star (either an ordinary 50-100 Msun star or a more massive one, formed through runaway mergers of ordinary stars in the core of a young massive star cluster). In this process, one of the binary components becomes gravitationally bound to the very massive star, while the second one is ejected, sometimes with a high speed. We performed three-body scattering experiments and found that early B-type stars (the progenitors of the majority of neutron stars) can be ejected with velocities of $\ga$ 200-400 km/s (typical of pulsars), while 3-4 Msun stars can attain velocities of $\ga$ 300-400 km/s (typical of the bound population of halo late B-type stars). We also found that the ejected stars can occasionally attain velocities exceeding the Milky Ways's escape velocity.

The atmosphere of the extremely high-velocity (530-920 km/s) early B-type star HD271791 is enriched in $\alpha$-process elements, which suggests that this star is a former secondary component of a massive tight binary system and that its surface was polluted by the nucleosynthetic products after the primary star exploded in a supernova. It was proposed that the (asymmetric) supernova explosion unbind the system and that the secondary star (HD271791) was released at its orbital velocity in the direction of Galactic rotation. In this Letter we show that to explain the Galactic rest-frame velocity of HD271791 within the framework of the binary-supernova scenario, the stellar remnant of the supernova explosion (a $\la$ 10 Msun black hole) should receive an unrealistically large kick velocity of $\geq$ 750-1200 km/s$. We therefore consider the binary-supernova scenario as highly unlikely and instead propose that HD271791 attained its peculiar velocity in the course of a strong dynamical three- or four-body encounter in the dense core of the parent star cluster. Our proposal implies that by the moment of encounter HD271791 was a member of a massive post-supernova binary.

N-body simulations have shown that the dynamical decay of the young (~1 Myr) Orion Nebula cluster could be responsible for the loss of at least half of its initial content of OB stars. This result suggests that other young stellar systems could also lose a significant fraction of their massive stars at the very beginning of their evolution. To confirm this expectation, we used the Mid-Infrared Galactic Plane Survey (completed by the Midcourse Space Experiment satellite) to search for bow shocks around a number of young ($\la$ several Myr) clusters and OB associations. We discovered dozens of bow shocks generated by OB stars running away from these stellar systems, supporting the idea of significant dynamical loss of OB stars. In this paper, we report the discovery of three bow shocks produced by O-type stars ejected from the open cluster NGC6611 (M16). One of the bow shocks is associated with the O9.5Iab star HD165319, which was suggested to be one of "the best examples for isolated Galactic high-mass star formation" (de Wit et al. 2005, A&A, 437, 247). Possible implications of our results for the origin of field OB stars are discussed.

The astrometric data on the runaway star BD+43 3654 are consistent with the origin of this O4If star in the center of the Cyg OB2 association, while BD+43 3654 is younger than the association. To reconcile this discrepancy, we suggest that BD+43 3654 is a blue straggler formed via a close encounter between two tight massive binaries in the core of Cyg OB2. A possible implication of this suggestion is that the very massive (and therefore apparently very young) stars in Cyg OB2 could be blue stragglers as well. We also suggest that the binary-binary encounter producing BD+43 3654 might be responsible for ejection of two high-velocity stars (the stripped helium cores of massive stars) -- the progenitors of the pulsars B2020+28 and B2021+51.

Recent proper motion and parallax measurements for the pulsar PSR B1508+55 indicate a transverse velocity of ~1100 km/s, which exceeds earlier measurements for any neutron star. The spin-down characteristics of PSR B1508+55 are typical for a non-recycled pulsar, which implies that the velocity of the pulsar cannot have originated from the second supernova disruption of a massive binary system. The high velocity of PSR B1508+55 can be accounted for by assuming that it received a kick at birth or that the neutron star was accelerated after its formation in the supernova explosion. We propose an explanation for the origin of hyperfast neutron stars based on the hypothesis that they could be the remnants of a symmetric supernova explosion of a high-velocity massive star which attained its peculiar velocity (similar to that of the pulsar) in the course of a strong dynamical three- or four-body encounter in the core of dense young star cluster. To check this hypothesis we investigated three dynamical processes involving close encounters between: (i) two hard massive binaries, (ii) a hard binary and an intermediate-mass black hole, and (iii) a single star and a hard binary intermediate-mass black hole. We find that main-sequence O-type stars cannot be ejected from young massive star clusters with peculiar velocities high enough to explain the origin of hyperfast neutron stars, but lower mass main-sequence stars or the stripped helium cores of massive stars could be accelerated to hypervelocities. Our explanation for the origin of hyperfast pulsars requires a very dense stellar environment of the order of 10^6 -10^7 stars pc^{-3}. Although such high densities may exist during the core collapse of young massive star clusters, we caution that they have never been observed.

We propose an explanation for the origin of hyperfast neutron stars (e.g. PSR B1508+55, PSR B2224+65, RX J0822-4300) based on the hypothesis that they could be the remnants of a symmetric supernova explosion of a high-velocity massive star (or its helium core) which attained its peculiar velocity (similar to that of the neutron star) in the course of a strong three- or four-body dynamical encounter in the core of a young massive star cluster. This hypothesis implies that the dense cores of star clusters (located either in the Galactic disk or near the Galactic centre) could also produce the so-called hypervelocity stars -- the ordinary stars moving with a speed of ~1000 km/s.

The proper motion of massive stars could cause them to explode far from the geometric centers of their wind-driven bubbles and thereby could affect the symmetry of the resulting diffuse supernova remnants. We use this fact to explain the origin of SNRs consisting of two partially overlapping shells (e.g. Cygnus Loop, 3C 400.2, etc.).

Astrometric data on the pulsars B2020+28 and B2021+51 suggest that they originated within several parsecs of each other in the direction of the Cyg OB2 association. It was proposed that the pulsars share their origin in a common massive binary and were separated at the birth of the second pulsar following the asymmetric supernova explosion. We consider a different scenario for the origin of the pulsar pair based on a possibility that the pulsars were separated before their birth and that they are the remnants of runaway stars ejected (with velocities similar to those of the pulsars) from the core of Cyg OB2 due to strong three- or four-body dynamical encounters. Our scenario does not require any asymmetry in supernova explosions.

The supernova remnant S147 harbors the pulsar PSR J0538+2817 whose characteristic age is more than an order of magnitude greater than the kinematic age of the system (inferred from the angular offset of the pulsar from the geometric center of the supernova remnant and the pulsar proper motion). To reconcile this discrepancy we propose that PSR J0538+2817 could be the stellar remnant of the first supernova explosion in a massive binary system and therefore could be as old as its characteristic age. Our proposal implies that S147 is the diffuse remnant of the second supernova explosion (that disrupted the binary system) and that a much younger second neutron star (not necessarily manifesting itself as a radio pulsar) should be associated with S147. We use the existing observational data on the system to suggest that the progenitor of the supernova that formed S147 was a Wolf-Rayet star (so that the supernova explosion occurred within a wind bubble surrounded by a massive shell) and to constrain the parameters of the binary system. We also restrict the magnitude and direction of the kick velocity received by the young neutron star at birth and find that the kick vector should not strongly deviate from the orbital plane of the binary system.

It is shown that the gravitational potential outside an inhomogeneous ellipsoid of revolution (spheroid) whose isodensity surfaces are confocal spheroids is identical to the gravitational potential of a homogeneous spheroid of the same mass.

A scenario for the origin of the system PSR B1757-24/supernova remnant (SNR) G5.4-1.2 is proposed. It is suggested that both objects are the remnants of a supernova (SN) that exploded within a pre-existing bubble blown-up by a runaway massive star (the SN progenitor) during the final (Wolf-Rayet) phase of its evolution. This suggestion implies that (a) the SN blast centre was significantly offset from the geometric centre of the wind-blown bubble (i.e. from the centre of the future SNR), (b) the bubble was surrounded by a massive wind-driven shell, and (c) the SN blast wave was drastically decelerated by the interaction with the shell. Therefore, one can understand how the relatively young and low-velocity pulsar PSR B1757-24 was able to escape from the associated SNR G5.4-1.2 and why the inferred vector of pulsar transverse velocity does not point away from the geometric centre of the SNR. A possible origin of the radio source G5.27-0.9 (located between PSR B1757-24 and the SNR G5.4-1.2) is proposed. It is suggested that G5.27-0.9 is a lobe of a low Mach number (\simeq 1.7) jet of gas outflowing from the interior of G5.4-1.2 through the hole bored in the SNR's shell by the escaping pulsar. It is also suggested that the non-thermal emission of the comet-shaped pulsar wind nebula originates in the vicinity of the termination shock and in the cylindric region of subsonically moving shocked pulsar wind. The role of magnetized wind-driven shells (swept-up during the Wolf-Rayet phase from the ambient interstellar medium with the regular magnetic field) in formation of elongated axisymmetric SNRs is discussed.

The association of PSR B1757-24 and the supernova remnant (SNR) G5.4-1.2 was recently questioned by Thorsett et al. (2002) on the basis of proper motion measurements of the pulsar and the "incorrect" orientation of the vector of pulsar transverse velocity [inferred from the orientation of the cometary-shaped pulsar wind nebula (PWN)]. We showed, however, that the association could be real if both objects are the remnants of an off-centred cavity supernova (SN) explosion.

We discuss a concept of off-centred cavity supernova explosion as applied to neutron star/supernova remnant associations and show how this concept could be used to preclude the anti-humane decapitating the Duck (G5.4-1.2 + G5.27-0.9) and dismembering the Swan (Cygnus Loop), as well as to search for a stellar remnant associated with the supernova remnant RCW86.

We report the results of a search for a point X-ray source (stellar remnant) in the southwest protrusion of the supernova remnant G 315.4-2.30 (MSH 14-63, RCW 86) using the archival data of the Chandra X-ray Observatory. The search was motivated by a hypothesis that G 315.4-2.30 is the result of an off-centered cavity supernova explosion of a moving massive star, which ended its evolution just near the edge of the main-sequence wind-driven bubble. This hypothesis implies that the southwest protrusion in G 315.4-2.30 is the remainder of a pre-existing bow shock-like structure created by the interaction of the supernova progenitor's wind with the interstellar medium and that the actual location of the supernova blast center is near the center of this hemispherical structure. We have discovered two point X-ray sources in the "proper" place. One of the sources has an optical counterpart with the photographic magnitude $13.38\pm0.40$, while the spectrum of the source can be fitted with an optically thin plasma model. We interpret this source as a foreground active star of late spectral type. The second source has no optical counterpart to a limiting magnitude $\sim 21$. The spectrum of this source can be fitted almost equally well with several simple models (power law: photon index $=1.87$; two-temperature blackbody: $kT_1 =0.11$ keV, $R_1 =2.34 $ km and $kT_2 =0.71$ keV, $R_2 =0.06$ km; blackbody plus power law: $kT =0.07$ keV, photon index $=2.3$). We interpret this source as a candidate stellar remnant (neutron star), while the photon index and non-thermal luminosity of the source (almost the same as those of the Vela pulsar and the recently discovered pulsar PSR J 0205+6449 in the supernova remnant 3C 58) suggest that it can be a young "ordinary" pulsar.

The possible association of the supernova remnant (SNR) G 343.1-2.3 with the pulsar PSR B 1706-44 (superposed on the arclike "shell" of the SNR) has been questioned by some authors on the basis of an inconsistency between the implied and measured (scintillation) transverse velocities of the pulsar, the absence of any apparent interaction between the pulsar and the SNR's "shell", and some other indirect arguments. We suggest, however, that this association could be real if both objects are the remnants of a supernova (SN) which exploded within a mushroom-like cavity (created by the SN progenitor wind breaking out of the parent molecular cloud). This suggestion implies that the actual shape of the SNR's shell is similar to that of the well-known SNR VRO 42.05.01 and that the observed bright arc corresponds to the "half" of the SNR located inside the cloud. We report the discovery in archival radio data of an extended ragged radio arc to the southeast of the bright arc which we interpret as the "half" of the SN blast wave expanding in the intercloud medium.

We propose a new approach for studying the neutron star/supernova remnant associations, based on the idea that the supernova remnants (SNRs) can be products of an off-centered supernova (SN) explosion in a preexisting bubble created by the wind of a moving massive star. A cavity SN explosion of a moving star results in a considerable offset of the neutron star (NS) birth-place from the geometrical center of the SNR. Therefore: a) the high transverse velocities inferred for a number of NSs through their association with SNRs can be reduced; b) the proper motion vector of a NS should not necessarily point away from the geometrical center of the associated SNR. Taking into account these two facts allow us to enlarge the circle of possible NS/SNR associations, and could significantly affect the results of previous studies of associations. The possibilities of our approach are illustrated with some examples. We also show that the concept of an off-centered cavity SN explosion could be used to explain the peculiar structures of a number of SNRs and for searches for stellar remnants possibly associated with them.

We present some arguments in support of the association of the pulsar PSR B1706-44 with the supernova remnant G343.1-2.3, based on the idea that these objects could be the result of a supernova explosion within a mushroom-like cavity (created by the supernova progenitor wind breaking out of the parent molecular cloud). We suggest that in addition to the known bright "half" of G343.1-2.3 there should exist a more extended and weaker component, such that the actual shape of G343.1-2.3 is similar to that of the well-known SNR VRO 42.05.01. We have found such a component in archival radio data.

It is generally accepted that the PSR B1509-58 is associated with the supernova remnant (SNR) MSH 15-52 (G 320.4-01.2). The spin-down age of the pulsar is \simeq 1700 years, while the size and the general appearance of the SNR suggest that this system is much older. A few possible explanations of this discrepancy have been put forward. We offer an alternative one and suggest that the high spin-down rate of the pulsar characterizes only a relatively short period of its (present) spin history, and that the enhanced braking torque is connected with the interaction between the pulsar's magnetosphere and the dense matter of a circumstellar clump (created during the late evolutionary stages of the supernova (SN) progenitor star). Our suggestion implies that the "true" age of PSR B1509-58 could be much larger than the spin-down age, and therefore the SNR MSH 15-52 is a middle-aged remnant similar to the Vela SNR (G 263.9-3.3). We also suggest that the dense (neutral) gas of the circumstellar clump could be responsible for the enhanced neutral hydrogen absorption towards PSR B1509-58, and that the optical emission of an optical counterpart for PSR B1509-58 should rather be attributed to a bow shock around this pulsar than to the pulsar itself.

It is pointed out that a cavity supernova (SN) explosion of a moving massive star could result in a significant offset of the neutron star (NS) birth-place from the geometrical centre of the supernova remnant (SNR). Therefore: a) the high implied transverse velocities of a number of NSs (e.g. PSR B1610-50, PSR B1757-24, SGR0525-66) could be reduced; b) the proper motion vector of a NS should not necessarily point away from the geometrical centre of the associated SNR; c) the circle of possible NS/SNR associations could be enlarged. An observational test is discussed, which could allow to find the true birth-places of NSs associated with middle-aged SNRs, and thereby to get more reliable estimates of their transverse velocities.

It is shown that if the shell of the Vela supernova remnant is responsible for nearly all the scattering of the Vela pulsar, then the scintillation and proper motion velocities of the pulsar can only be reconciled with each other in the case of nonzero transverse velocity of the scattering material. A possible origin of large-scale transverse motions in the shell of the Vela supernova remnant is discussed.

Most of middle-aged supernova remnants (SNRs) have a distorted and complicated appearance which cannot be explained in the framework of the Sedov-Taylor model. We consider three typical examples of such SNRs (Vela SNR, MSH 15-52, G 309.2-00.6) and show that their structure could be explained as a result of interaction of a supernova (SN) blast wave with the ambient medium preprocessed by the action of the SN progenitor's wind and ionized emission.

The nature of the Vela X-ray "jet", recently discovered by Markwardt & \"Ogelman (1995), is examined. It is suggested that the "jet" arises along the interface of domelike deformations of the Rayleigh-Taylor unstable shell of the Vela supernova remnant; thereby the "jet" is interpreted as a part of the general shell of the remnant. The origin of deformations as well as the general structure of the remnant are discussed in the framework of a model based on a cavity explosion of a supernova star. It is suggested that the shell deformations viewed at various angles appear as filamentary structures visible throughout the Vela supernova remnant at radio, optical, and X-ray wavelengths. A possible origin of the nebula of hard X-ray emission detected by Willmore et al. (1992) around the Vela pulsar is proposed.