A previous spectroscopic study identified the very massive O2 III star VFTS 16 in the Tarantula Nebula as a runaway star based on its peculiar line-of-sight velocity. We use the Gaia DR2 catalog to measure the relative proper motion of VFTS 16 and nearby bright stars to test if this star might have been ejected from the central cluster, R136, via dynamical ejection. We find that the position angle and magnitude of the relative proper motion (0.338 +/- 0.046 mas/yr, or approximately 80 +\- 11 km/s) of VFTS 16 are consistent with ejection from R136 approximately 1.5 +/- 0.2 Myr ago, very soon after the cluster was formed. There is some tension with the presumed age of VFTS 16 that, from published stellar parameters, cannot be greater than 0.9 +0.3/-0.2 Myr. Older ages for this star would appear to be prohibited due to the absence of He I lines in its optical spectrum, since this sets a firm lower limit on its effective temperature. The dynamical constraints may imply an unusual evolutionary history for this object, perhaps indicating it is a merger product. Gaia DR2 also confirms that another very massive star in the Tarantula Nebula, VFTS 72 (alias BI253; O2 III-V(n)((f*)), is also a runaway on the basis of its proper motion as measured by Gaia. While its tangential proper motion (0.392 +/-0.062 mas/yr or 93 +/-15 km/s) would be consistent with dynamical ejection from R136 approximately 1 Myr ago, its position angle is discrepant with this direction at the 2$\sigma$ level. From their Gaia DR2 proper motions we conclude that the two ~100 solar mass O2 stars, VFTS 16 and VFTS72, are fast runaway stars, with space velocities of around 100 km/s relative to R136 and the local massive star population. The dynamics of VFTS16 are consistent with it having been ejected from R136, and this star therefore sets a robust lower limit on the age of the central cluster of ~1.3 Myr.

We search for runaway stars in the Large Magellanic Cloud (LMC) by computing the space velocities of the visually brightest stars in the LMC that are included in the Gaia TGAS proper motion catalog. We compare with predictions from stellar dynamical models to obtain (peculiar) velocities relative to their local stellar environment. Two of the 31 stars have unusually high proper motions. Of the remaining 29 stars we find that most objects in this sample have velocities in very good agreement with model predictions of a circularly rotating disk model. Indeed the excellent fit to the model implies that the TGAS uncertainty estimates are likely overestimated. The fastest outliers in this subsample contain the LBV R71 and a few other well known emission line objects though in no case do we derive velocities consistent with fast (~100 km/s) runaways. Our results imply that R 71 in particular has a moderate deviation from the local stellar velocity field (40 km/s) lending support to the proposition that this object cannot have evolved as a normal single star since it lies too far from massive star forming complexes to have arrived at its current position during its lifetime. Our findings therefore strengthen the case for this LBV being the result of binary evolution. Of the two stars with unusually high proper motions we find that one, the isolated B1.5 Ia+ supergiant Sk-67 2 (HIP 22237), is a candidate hypervelocity star, the TGAS proper motion implying a very large peculiar transverse velocity (~360 km/s) directed radially away from the LMC centre. If confirmed, for example by Gaia Data Release 2, it would imply that this massive supergiant, on the periphery of the LMC, is leaving the galaxy where it will explode as a supernova.

We present chemical and radiative transfer models for the many far-IR ortho- and para-H2O lines that were observed from the Orion-KL region in high resolution Fabry-Perot (FP) mode by the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO). The chemistry of the region was first studied by simulating the conditions in the different known components of Orion-KL: chemical models for a hot core, a plateau and a ridge were coupled with an accelerated Lambda$-iteration (ALI) radiative transfer model to predict H2O line fluxes and profiles. Our models include the first 45 energy levels of ortho- and para-H2O. We find that lines arising from energy levels below 560 K were best reproduced by a gas of density 3x10^5 cm^-3 at a temperature of 70-90 K, expanding at a velocity of 30 km s^-1 and with a H2O/H2 abundance ratio of the order of 2 - 3 x 10^-5, similar to the abundance derived by Cernicharo et al. (2006). However, the model that best reproduces the fluxes and profiles of H2O lines arising from energy levels above 560 K has a significantly higher H2O/H2 abundance, 1 - 5 x 10^-4, originating from gas of similar density, in the Plateau region, that has been heated to 300 K, relaxing to 90-100 K. We conclude that the observed water lines do not originate from high temperature shocks.

As part of the first high resolution far-IR spectral survey of the Orion KL region (Lerate et al. 2006), we observed 20 CO emission lines with Jup=16 to Jup=39 (upper levels from approx 752 K to 4294 K above the ground state). Observations were taken using the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO), in its high resolution Fabry-Perot (FP) mode (approx 33 km s$^{-1}$). We present here an analysis of the final calibrated CO data, performed with a more sophisticated modelling technique than hitherto, including a detailed analysis of the chemistry, and discuss similarities and differences with previous results. The inclusion of chemical modelling implies that atomic and molecular abundances are time-predicted by the chemistry. This provides one of the main differences with previous studies in which chemical abundances needed to be assumed as initial condition. The chemistry of the region is studied by simulating the conditions of the different known components of the KL region: chemical models for a hot core, a plateau and a ridge are coupled with an accelerated Lambda-iteration (ALI)radiative transfer model to predict line fluxes and profiles. We conclude that the CO transitions with 18<Jup<25 mainly arise from a hot core of diameter 0.02 pc and a density of 10$^{7}$ cm$^{-3}$ rather from the plateau as previous studies had indicated.

We present the detection and modeling of more than 70 far-IR pure rotational lines of water vapor, including the 18O and 17O isotopologues, towards Orion KL. Observations were performed with the Long Wavelength Spectrometer Fabry-Perot (LWS/FP; R~6800-9700) on board the Infrared Space Observatory (ISO) between ~43 and ~197 um. The water line profiles evolve from P-Cygni type profiles (even for the H2O18 lines) to pure emission at wavelengths above ~100 um. We find that most of the water emission/absorption arises from an extended flow of gas expanding at 25+-5 kms^-1. Non-local radiative transfer models show that much of the water excitation and line profile formation is driven by the dust continuum emission. The derived beam averaged water abundance is 2-3x10^-5. The inferred gas temperature Tk=80-100 K suggests that: (i) water could have been formed in the "plateau" by gas phase neutral-neutral reactions with activation barriers if the gas was previously heated (e.g. by shocks) to >500 K and/or (ii) H2O formation in the outflow is dominated by in-situ evaporation of grain water-ice mantles and/or (iii) H2O was formed in the innermost and warmer regions (e.g. the hot core) and was swept up in ~1000 yr, the dynamical timescale of the outflow.

We have carried out a high spectral resolution line survey towards the Orion Kleinmann-Low (KL) cluster from 44-188 um. The observations were taken with the Long Wavelength Spectrometer (LWS) in Fabry-Perot mode, on board the Infrared Space Observatory (ISO). A total of 152 lines are clearly detected and a further 34 features are present as possible detections. The spectrum is dominated by the molecular species H2O, OH and CO, along with [OI] and [CII] lines from PDR or shocked gas and [OIII], [NIII] lines from the foreground M42 HII region. Several isotopic species, as well as NH3, are also detected. HDO and H3O+ are tentatively detected for the first time in the far-infrared range towards Orion-KL. A basic analysis of the line observations is carried out, by comparing with previous measurements and published models and deriving rotational temperatures and column densities in the case of the molecular species. The complexity of the region requires more sophisticated models for the interpretation of all the line observations.

As part of the first far-IR line survey towards Orion KL, we present the detection of seven new rotationally excited OH Lambda-doublets (at 48, 65, 71, 79, 98 and 115 um). Observations were performed with the Long Wavelength Spectrometer (LWS) Fabry-Perots on board the Infrared Space Observatory (ISO). In total, more than 20 resolved OH rotational lines, with upper energy levels up to 620 K, have been detected at an angular and velocity resolutions of 80$'' and 33 km s^-1 respectively. OH line profiles show a complex behavior evolving from pure absorption, P-Cygni type to pure emission. We also present a large scale 6' declination raster in the OH ^2\Pi_3/2 J=5/2^+-3/2^- and ^2\Pi_3/2 J=7/2^-5/2^+ lines (at 119.441 and 84.597 um) revealing the decrease of excitation outside the core of the cloud. From the observed profiles, mean intrinsic line widths and velocity offsets between emission and absorption line peaks we conclude that most of the excited OH arises from Orion outflow(s), i.e. the ``plateau'' component. We determine an averaged OH abundance relative to H_2 of X(OH)=(0.5-1.0)x10^-6, a kinetic temperature of 100 K and a density of n(H_2)=5x10^5 cm^-3. Even with these conditions, the OH excitation is heavily coupled with the strong dust continuum emission from the inner hot core regions and from the expanding flow itself.