EZ Canis Majoris is a classical Wolf-Rayet star whose binary nature has been debated for decades. It was recently modeled as an eccentric binary with a periodic brightening at periastron of the emission originating in a shock heated zone near the companion. The focus of this paper is to further test the binary model and to constrain the nature of the unseen close companion by searching for emission arising in the shock-heated region. We analyze over 400 high resolution the International Ultraviolet Explorer spectra obtained between 1983 and 1995 and XMM-Newton observations obtained in 2010. The light curve and radial velocity (RV) variations were fit with the eccentric binary model and the orbital elements were constrained. We find RV variations in the primary emission lines with a semi-amplitude K$_1\sim$30 km/s in 1992 and 1995, and a second set of emissions with an anti-phase RV curve with K$_2\sim$150 km/s. The simultaneous model fit to the RVs and the light curve yields the orbital elements for each epoch. Adopting a Wolf-Rayet mass M$_1\sim$20 M$_\odot$ leads to M$_2\sim$3-5 M$_\odot$, which implies that the companion could be a late B-type star. The eccentric (e=0.1) binary model also explains the hard X-ray light curve obtained by XMM-Newton and the fit to these data indicates that the duration of maximum is shorter than the typical exposure times. The anti-phase RV variations of two emission components and the simultaneous fit to the RVs and the light curve are concrete evidence in favor of the binary nature of EZ Canis Majoris. The assumption that the emission from the shock-heated region closely traces the orbit of the companion is less certain, although it is feasible because the companion is significantly heated by the WR radiation field and impacted by the WR wind.

Context: GR290 (M33 V0532=Romano's star) is a luminous M33 object undergoing photometric variability typical for luminous blue variable (LBV) stars. It lies inside Wolf-Rayet region in the Hertzsprung-Russell diagram and possesses a WN8 type spectrum at the light minima. Analysis of Gran Telescopio Canarias (GTC) spectra obtained in 2016 led to the conclusion that it is surrounded by an unresolved HII region formed mostly of ejected material from the central star, and disclosed the presence of a second, more extended asymmetrical emission region. Aims: The aim of this paper is to further explore the structure of the nearby environment of GR290. Methods: Long-slit spectra of GR290 were obtained with three slit orientations in the visual and red spectral regions. The emission-line distribution for each slit was analyzed. Results: We confirm the presence of an asymmetric HII region that extends ~50 pc to the south; ~30pc to the north and southeast; ~20 pc to the east and northwest and ~10pc to the west. We also present the first spectrum to be acquired of a star belonging to the neighboring OB88 association, J013501.87+304157.3, which we classify as a B-type supergiant with a possible binary companion.

GR 290 (M 33 V0532 = Romano's Star) is a unique variable star in the M33 galaxy, which simultaneously displays variability typical for luminous blue variable (LBV) stars and physical parameters typical for nitrogen-rich Wolf-Rayet (WR) stars (WN). As of now, GR 290 is the first object which is confidently classified as a post-LBV star. In this paper, we outline the main results achieved from extensive photometric and spectroscopic observations of the star: the structure and chemical composition of its wind and its evolution over time, the systematic increase of the bolometric luminosity during the light maxima, the circumstellar environment. These results show that the current state of Romano's Star constitutes a fundamental link in the evolutionary path of very massive stars.

Context. The variability of the Wolf-Rayet star EZ CMa has been documented for close to half a century, and a clear periodicity of $\sim$3.7 days is established. However, all attempts to prove that it is a binary have failed because the photometric, spectroscopic, and polarimetric variations are not coherent over more than a few orbital cycles. Aims. In this letter we show that the lack of coherence in the variability can be explained with a very rapid apsidal motion in a binary orbit.} Methods. We measured the times of minima in a recently published exceptionally long photometric light curve obtained by the Toronto {\emph{BRITE}} satellite. The apsidal motion and the system eccentricity are determined from the length of the time intervals between these minima, which alternate in their duration, following a pattern that is clearly associated with apsidal motion. These minima are superposed on brightness enhancements of the emission from a shock zone, which occur at about the times of periastron phases. Results. We determine the orbital periodicity, $P_{a}=3.63\, $d, and the period of the apsidal motion, $U\simeq 100\, $d, which together yield an average sidereal period of $P_{s}=3.77\,$d. The eccentricity is found to be close to 0.1. The rate of periapsis retreat changes significantly over the period of observation and is determined to be $-16^\circ\,\mathrm{P}^{-1}_a$ at the beginning of the observing period and $-10^\circ\,\mathrm{P}^{-1}_a$ at the end. Conclusions. We demonstrate that by introducing a fast apsidal motion, the basic photometric variability is very well explained. The binary nature of EZ CMa is now established. This might imply that other apparently single Wolf-Rayet stars that emit hard X-rays, similar to EZ CMa, are also binaries.

Context: GR290 (M33/V532=Romano's Star) is a suspected post-LBV star located in M33 galaxy that shows a rare Wolf-Rayet spectrum during its minimum light phase. In spite of many studies, its atmospheric structure, its circumstellar environment and its place in the general context of massive stars evolution is poorly known. Aims: Detailed study of its wind and mass loss, and study of the circumstellar environment associated to the star. Methods: Long-slit spectra of GR290 were obtained during its present minimum luminosity phase with the GTC together with contemporaneous BVRI photometry. The data were compared with non-LTE model atmosphere synthetic spectra computed with CMFGEN and with CLOUDY models for ionized interstellar medium regions. Results: The current $m_V=18.8$ mag, is the faintest at which this source has ever been observed. The non-LTE models indicate effective temperature $T_{eff}$=27-30 kK at radius $R_{2/3}$=27-21 Rsun and mass loss rate $\dot{M}=1.5\times10^{-5}$ Msun yr$^{-1}$. The terminal wind speed $V_\infty$=620 ${\rm km~s^{-1}}$ is faster than ever before recorded while the current luminosity $L_*=(3.1-3.7)\times 10^5$ Lsun is the lowest ever deduced. It is overabundant in He and N and underabundant in C and O. It is surrounded by an unresolved compact HII region with dimensions $\leq$4 pc, from where H-Balmer, HeI lines and [OIII] and [NII] are detected. In addition, we find emission from a more extended interstellar medium (ISM) region which appears to be asymmetric, with a larger extent to the East (16-40 pc) than to the West. Conclusions: In the present long lasting visual minimum, GR290 is in a lower bolometric luminosity state with higher mass loss rate. The nearby nebular emission seems to suggest that the star has undergone significant mass loss over the past $10^4-10^5$ years and is nearing the end stages of its evolution.

$\eta$ Car is a massive, eccentric binary with a rich observational history. We obtained the first high-cadence, high-precision light curves with the BRITE-Constellation nanosatellites over 6 months in 2016 and 6 months in 2017. The light curve is contaminated by several sources including the Homunculus nebula and neighboring stars, including the eclipsing binary CPD$-$59$^\circ$2628. However, we found two coherent oscillations in the light curve. These may represent pulsations that are not yet understood but we postulate that they are related to tidally excited oscillations of $\eta$ Car's primary star, and would be similar to those detected in lower-mass eccentric binaries. In particular, one frequency was previously detected by van Genderen et al. and Sterken et al. through the time period of 1974 to 1995 through timing measurements of photometric maxima. Thus, this frequency seems to have been detected for nearly four decades, indicating that it has been stable in frequency over this time span. These pulsations could help provide the first direct constraints on the fundamental parameters of the primary star if confirmed and refined with future observations.

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

Abridged: Alpha Virginis is a binary system whose proximity and brightness allow detailed investigations of the internal structure and evolution of stars undergoing time-variable tidal interactions. Previous studies have led to the conclusion that the internal structure of Spica's primary star may be more centrally condensed than predicted by theoretical models of single stars, raising the possibility that the interactions could lead to effects that are currently neglected in structure and evolution calculations. The key parameters in confirming this result are the values of the orbital eccentricity $e$, the apsidal period $U$, and the primary star's radius, R_1. We analyze the impact that line profile variability has on the derivation of its orbital elements and R_1. We use high SNR observations obtained in 2000, 2008, and 2013 to derive the orbital elements from fits to the radial velocity curves. We produce synthetic line profiles using an ab initio tidal interaction model. Results: The variations in the line profiles can be understood in terms of the tidal flows, whose large-scale structure is relatively fixed in the rotating binary system reference frame. Fits to the radial velocity curves yield $e$=0.108$\pm$0.014. However, the analogous RV curves from theoretical line profiles indicate that the distortion in the lines causes the fitted value of $e$ to depend on the argument of periastron; i.e., on the epoch of observation. As a result, the actual value of $e$ may be as high as 0.125. We find that $U$=117.9$\pm$1.8, which is in agreement with previous determinations. Using the value $R_1=6.8 R_\odot$ derived by Palate et al. (2013) the value of the observational internal structure constant $k_{2,obs}$ is consistent with theory. We confirm the presence of variability in the line profiles of the secondary star.

Planetary systems are born in the disks of gas, dust and rocky fragments that surround newly formed stars. Solid content assembles into ever-larger rocky fragments that eventually become planetary embryos. These then continue their growth by accreting leftover material in the disc. Concurrently, tidal effects in the disc cause a radial drift in the embryo orbits, a process known as migration. Fast inward migration is predicted by theory for embryos smaller than three to five Earth masses. With only inward migration, these embryos can only rarely become giant planets located at Earth's distance from the Sun and beyond, in contrast with observations. Here we report that asymmetries in the temperature rise associated with accreting infalling material produce a force (which gives rise to an effect that we call "heating torque") that counteracts inward migration. This provides a channel for the formation of giant planets and also explains the strong planet-metallicity correlation found between the incidence of giant planets and the heavy-element abundance of the host stars.

New spectroscopic observations of the LBV/WR multiple system HD5980 in the Small Magellanic Cloud are used to address the question of the masses and evolutionary status of the two very luminous stars in the 19.3d eclipsing binary system. Two distinct components of the N V 4944 A line are detected in emission and their radial velocity variations are used to derive masses of 61 and 66 Mo, under the assumption that binary interaction effects on this atomic transition are negligible. We propose that this binary system is the product of quasi-chemically homogeneous evolution with little or no mass transfer. Thus, both of these binary stars may be candidates for gamma-ray burst progenitors or even pair instability supernovae. Analysis of the photospheric absorption lines belonging to the third-light object in the system confirm that it consists of an O-type star in a 96.56d eccentric orbit (e=0.82) around an unseen companion. The 5:1 period ratio and high eccentricities of the two binaries suggest that they may constitute a hierarchical quadruple system.

We discuss the potential role that tidal flows in asynchronous binary stars may play in transporting chemically enriched material from deep layers towards the surface and the corresponding observational consequences of these processes. We suggest that the time-dependent velocity field induced by asynchronous rotation may contribute significantly to the mixing, thus providing a channel for the formation of chemically enriched slowly rotating massive stars.

Context: The technique of matching synthetic spectra computed with theoretical stellar atmosphere models to the observations is widely used in deriving fundamental parameters of massive stars. When applied to binaries, however, these models generally neglect the interaction effects present in these systems Aims: The aim of this paper is to explore the uncertainties in binary stellar parameters that are derived from single-star models Methods: Synthetic spectra that include the tidal perturbations and irradiation effects are computed for the binary system alpha Virginis (Spica) using our recently-developed CoMBiSpeC model. The synthetic spectra are compared to S/N~2000 observations and optimum values of Teff and log(g) are derived. Results: The binary interactions have only a small effect on the strength of the photospheric absorption lines in Spica (<2% for the primary and <4% for the secondary). These differences are comparable to the uncertainties inherent to the process of matching synthetic spectra to the observations and thus the derived values of Teff and log(g) are unaffected by the binary perturbations. On the other hand, the interactions do produce significant phase-dependent line profile variations in the primary star, leading to systematic distortions in the shape of its radial velocity curve. Migrating sub-features ("bumps") are predicted by our model to be present in the same photospheric lines as observed, and their appearance does not require any a priori assumptions regarding non-radial pulsation modes. Matching the strength of lines in which the most prominent "bumps" occur requires synthetic spectra computed with larger "microturbulence" than that required by other lines.

Aims: We simulate the spectra of massive binaries at different phases of the orbital cycle, accounting for the gravitational influence of the companion star on the shape and physical properties of the stellar surface. Methods: We used the Roche potential modified to account for radiation pressure to compute the stellar surface of close circular systems and we used the TIDES code for surface computation of eccentric systems. In both cases, we accounted for gravity darkening and mutual heating generated by irradiation to compute the surface temperature. We then interpolated NLTE plane-parallel atmosphere model spectra in a grid to obtain the local spectrum at each surface point. We finally summed all contributions, accounting for the Doppler shift, limb-darkening, and visibility to obtain the total synthetic spectrum. We computed different orbital phases and sets of physical and orbital parameters. Results: Our models predict line strength variations through the orbital cycle, but fail to completely reproduce the Struve-Sahade effect. Including radiation pressure allows us to reproduce a surface temperature distribution that is consistent with observations of semi-detached binary systems. Conclusions: Radiation pressure effects on the stellar surface are weak in (over)contact binaries and well-detached systems but can become very significant in semi-detached systems. The classical von Zeipel theorem is sufficient for the spectral computation. Broad-band light curves derived from the spectral computation are different from those computed with a model in which the stellar surfaces are equipotentials of the Roche potential scaled by the instantaneous orbital separation. In many cases, the fit of two Gaussian/Lorentzian profiles fails to properly measure the equivalent width of the lines and leads to apparent variations that could explain some of the effects reported in the literature.

The mass of the neutron star in Vela X-1 has been found to be more massive than the canonical 1.5 Mo. This result relies on the assumption that the amplitude of the optical component's measured radial velocity curve is not seriously affected by the interactions in the system. In this paper we explore the effect on the radial velocity curve caused by surface motions excited by tidal interactions. We use a calculation from first principles that involves solving the equations of motion of a Lagrangian grid of surface elements. The velocities on the visible surface of the star are projected along the line-of-sight to the observer to obtain the absorption-line profile in the observer's reference frame. The centroid of the line-profiles for different orbital phases is then measured and a simulated RV curve constructed. Models are run for the "standard" (vsini=116 km/s) and "slow" (56 km/s) supergiant rotation velocities. We find that the surface velocity field is complex and includes fast, small-spatial scale structures. It leads to strong variability in the photospheric line profiles which, in turn, causes significant deviations from a Keplerian RV curve. The peak-to-peak amplitudes of model RV curves are in all cases larger than the amplitude of the orbital motion. Keplerian fits to RV curves obtained with the "standard" rotation velocity imply a neutron star >1.7 Mo. However, a similar analysis of the "slow" rotational velocity models allows for m_ns ~ 1.5 Mo. Thus, the stellar rotation plays an important role in determining the characteristics of the perturbed RV curve. Given the observational uncertainty in GP Vel's projected rotation velocity and the strong perturbations seen in the published and the model RV curves, we are unable to rule out a small (~1.5 Mo) mass for the neutron star companion.

Over the past 40 years, the massive LBV/WR system HD 5980 in the Small Magellanic Cloud has undergone a long-term S Doradus type variability cycle and two brief and violent eruptions in 1993 and 1994. In this paper we analyze a collection of UV and optical spectra obtained between 1979 and 2009 and perform CMFGEN model fits to spectra of 1994, 2000, 2002 and 2009. The results are as follows: a) The long term S Dor-type variability is associated with changes of the hydrostatic radius; b) The 1994 eruption involved changes in its bolometric luminosity and wind structure; c) the emission-line strength, the wind velocity and the continuum luminosity underwent correlated variations in the sense that a decreasing V$_\infty$ is associated with increasing emission line and continuum levels; and d) The spectrum of the third star in the system ({\it Star C}) is well-fit by a T$_{eff}$=32 K model atmosphere with SMC chemical abundances. For all epochs, the wind of the erupting star is optically thick at the sonic point and is thus driven mainly by the continuum opacity. We speculate that the wind switches between two stable regimes driven by the "hot" (during the eruption) and the "cool" (post-eruption) iron opacity bumps as defined by Lamers & Nugis (2002) and Gr\"afener and Hamann (2008), and thus the wind may undergo a bi-stability jump of a different nature from that which occurs in OB-stars.

A number of binary systems present evidence of enhanced activity around periastron passage, suggesting a connection between tidal interactions and these periastron effects. The aim of this investigation is to study the time-dependent response of a star's surface as it is perturbed by a binary companion. We derive expressions for the rate of dissipation, $\dot{E}$, of the kinetic energy by the viscous flows driven by tidal interactions on the surface layer. The method is tested by comparing the results from a grid of model calculations with the analytical predictions of Hut (1981) and the synchronization timescales of Zahn (1977, 2008). Our results for the orbital cycle averaged energy dissipation on orbital separation are consistent with those of Hut for model binaries with orbital separations at periastron >8 stellar radii. The model also reproduces the predicted pseudo-synchronization angular velocity for moderate eccentricities and the same scaling of synchronization timescales for circular orbits with separation as given by Zahn. The computations gives the distribution of $\dot{E}$ over the stellar surface, and show that it is generally concentrated at the equatorial latitude, with maxima generally located around four clearly defined longitudes, corresponding to the fastest azimuthal velocity perturbations. Maximum amplitudes occur around periastron passage or slightly thereafter for supersynchronously rotating stars. In very eccentric binaries, the distribution of $\dot{E}$ over the surface changes significantly as a function of orbital phase, with small spatial structures appearing after periastron. An exploratory calculation for the highly eccentric binary system delta Sco suggests that the sudden and large amplitude variations in surface properties around periastron may contribute toward the activity observed around this orbital phase.

We present the results of high precision, high resolution (R~68000) optical observations of the short-period (4d) eccentric binary system Alpha Virginis (Spica) showing the photospheric line-profile variability that in this system can be attributed to non-radial pulsations driven by tidal effects. Although scant in orbital phase coverage, the data provide S/N>2000 line profiles at full spectral resolution in the wavelength range delta-lambda = 4000--8500 Angstroms, allowing a detailed study of the night-to-night variability as well as changes that occur on ~2 hr timescale. Using an ab initio theoretical calculation, we show that the line-profile variability can arise as a natural consequence of surface flows that are induced by the tidal interaction.

Binary stars in eccentric orbits are frequently reported to present increasing levels of activity around periastron passage. In this paper we present results of a calculation from first principles of the velocity field on the surface of a star that is perturbed by a binary companion. This allows us to follow the orbital phase-dependence of the amount of kinetic energy that may be dissipated through the viscous shear, dot-E, driven by tidal interactions. For stars with relatively small stellar radii compared with the orbital separation at periastron (R/r_per< 0.14), a clear maximum occurs before periastron passage for sub-synchronous rotation and after periastron for super-synchronous rotation. For larger stellar radii however, dot-E oscillates over the orbital cycle and periastron passage does not cause a particularly greater enhancement in energy dissipation rates than some of the other orbital phases. Finally, we perform exploratory calculations for a WR/LBV binary system that in 1993-1994 underwent an LBV-like eruption, HD 5980. Our dot-E calculations reproduce the oscillations that appear around periastron passage in HD5980's recent visual light curve. We suggest that the energy dissipation from tidal flows in asynchronously rotating binary stars may provide a mechanism by which evolving stars may be driven into an active state. Given the nature of the tidal perturbations, the resulting mass-loss distribution is expected to be non-uniform over the stellar surface and highly time-dependent.

Asynchronous rotation in binary stars produces non-radial oscillations that are known to cause observable variability on orbital timescales. The horizontal perturbations of the surface velocity fields are referred to as "tidal flows". In this paper we illustrate the manner in which tidal flows perturb the surface velocity field from that of uniform rotation, using a one-layer stellar model for the calculations. We justify the validity of this simplified model by the striking similarity between the photospheric absorption line-profiles it predicts and observational data of the binary system alpha Virginis. The velocity perturbations are used to compute the mechanical energy dissipation rates, dot-E, due to the shearing flows for the case of a massive (50+28 Mo) binary system having a moderately eccentric (e=0.3) orbit. The largest value of dot-E around periastron phases is found on the hemisphere facing the companion. However, at other orbital phases the maximum dot-E may migrate towards the poles. Assuming that dot-E plays a role in the mass-loss characteristics of massive binary systems, this suggests that peculiar binaries such as HD 5980 and eta Carinae may have a highly non-spherically symmetric mass-loss distribution which, in addition, is time-variable.

We present the results of optical wavelength observations of the unusual SMC eclipsing binary system HD 5980 obtained in 1999 and 2004--2005. Radial velocity curves for the erupting LBV/WR object (star A) and its close WR-like companion (star B) are obtained by deblending the variable emission-line profiles of N IV and N V lines under the simplistic assumption that these lines originate primarily in the winds of star A and star B. The derived masses M_A=58--79 Mo and M_B=51--67 Mo, are more consistent with the stars' location near the top of the HRD than previous estimates. The presence of a wind-wind interaction region is inferred from the orbital phase-dependent behavior of He I P Cygni absorption components. The emission-line intensities continued with the declining trend previously seen in UV spectra. The behavior of the photospheric absorption lines is consistent with the results of Schweickhardt (2002) who concludes that the third object in the combined spectrum, star C, is also a binary system with P(starC)~96.5 days, e=0.83. The data used in this paper will be made publicly available for further analysis.

The LBV/WR system HD 5980 contains a short-period, eccentric binary system with interacting stellar winds. In this paper we summarize results from model calculations of the tidal flows on the LBV component showing that energy dissipation rates, E-dot, associated with turbulent viscosity are orbital-phase dependent as well as variable over the stellar surface. We speculate that if E-dot contributes towards driving mass-loss, the strongest wind-wind interaction effects may occur after periastron passage. In addition, the model suggests the presence of stronger outflows localized at polar angles theta~30 --50 degrees during part of the orbital cycle. Thus, the analysis of wind-wind interactions in this system requires that models be revised to incorporate non-stationary and asymmetric wind structures.

Stellar oscillations are excited in non-synchronously rotating stars in binary systems due to the tidal forces. Tangential components of the tides can drive a shear flow which behaves as a differentially forced rotating structure in a stratified outer medium. In this paper we show that our single-layer approximation for the calculation of the forced oscillations yields results that are consistent with the predictions for the synchronization timescales in circular orbits. In addition, calibrating our model results to fit Zahn's relationship between synchronization timescales and orbital separation, we are able to constrain the value of the kinematical viscosity parameter. For the 4Mo+5Mo binary modeled in this paper, the kinematical viscosity is in the range 0.0015 -- 0.0043 Ro^2/d for orbital periods in the range 2.5 -- 25 d. Furthermore, the energy dissipation rate due to the shear flow is found to decrease by ~2 orders of magnitude as synchronization is approached, implying that binary systems may approach synchronization relatively quickly but that it takes a much longer timescale to actually attain this condition. Our model can be used to calculate the energy dissipation rates throughout the orbital cycle for arbitrary values of eccentricity and stellar rotational velocity. We suggest that the asymmetric distribution of energy dissipation over the stellar surface may lead to the appearance of localized regions of enhanced surface activity.

The X-ray source 2S0114+650=LSI+65 010 is a binary system containing a B-type primary and a low mass companion believed to be a neutron star. The system has three reported periodicities: the orbital period, P{orb}~11.6 d, X-ray flaring with P{flare}~2.7 hr, and a "superorbital" X-ray periodicity P{super}~30.7 d. The objective of this paper is to show that the puzzling periodicities in the system may be explained in the context of scenarios in which tidal interactions drive oscillations in the B-supergiant star. We calculate the solution of the equations of motion for one layer of small surface elements distributed along the equator of the star, as they respond to the forces due to gas pressure, centrifugal, coriolis, viscous forces, and the gravitational forces of both stars. This calculation provides variability timescales that can be compared with the observations. In addition, we use observational data obtained at the Observatorio Astron\'omico Nacional en San Pedro M\'artir (OAN/SPM) between 1993-2004 to determine which of the periodicities may be present in the optical region. We suggest that the tidal oscillations lead to a structured stellar wind which, when fed to the neutron star, produces the X-ray modulations. The connection between the stellar oscillations and the modulation of the mass ejection may lie in the shear energy dissipation generated by the tangential motions that are produced by the tidal interaction, particularly in the tidal bulge region. The tidal oscillation scenario weakens the case for 2S0114+650 containing a magnetar descendent.

HD 5980 is a unique system containing one massive star (star A) that is apparently entering the luminous blue variable phase, and an eclipsing companion (star B) that may have already evolved beyond this phase to become a Wolf-Rayet star. In this paper we present the results from FUSE observations obtained in 1999, 2000, and 2002 and one far-UV observation obtained by ORFEUS/BEFS in 1993 shortly before the first eruption of HD 5980. The eight phase-resolved spectra obtained by FUSE in 2002 are analyzed in the context of a wind-eclipse model. This analysis shows that the wind of the eruptor obeyed a very fast velocity law in 2002, which is consistent with the line-driving mechanism. Large amplitude line-profile variations on the orbital period are shown to be due to the eclipse of star B by the wind of star A, although the eclipse due to gas flowing in the direction of star B is absent. This can only be explained if the wind of star A is not spherically symmetric, or if the eclipsed line radiation is "filled-in" by emission originating from somewhere else in the system, e.g., in the wind-wind collision region. Except for a slightly lower wind speed, the ORFEUS/BEFS spectrum is very similar to the spectrum obtained by FUSE at the same orbital phase: there is no indication of the impending eruption. However, the trend for decreasing wind velocity suggests the occurrence of the "bi-stability" mechanism, which in turn implies that the restructuring of the circumbinary environment caused by the transition from "fast, rarefied wind" to "slow, dense wind" was observed as the eruptive event. The underlying mechanism responsible for the long-term decrease in wind velocity that precipitated this change remains an open issue.

We present high resolution spectroscopic observations of the massive X-ray binary system LS I+65 010=2S 0114+650 in the optical wavelength region. A correlation between equivalent width and radial velocity of photospheric absorption lines is found. The systemic velocity, inferred from the weaker lines is v_helio = -31 +/- 5 km s^-1, which, if attributed solely to the Galactic rotation curve, implies that LS I+65 010 lies within 3 kpc from the Sun. The ISM Na I D lines display 2 resolved high velocity components at v_helio = -70, -48 km s^-1, possibly associated with gas surrounding the binary system, in addition to the -24, -8 km s^-1 ISM features due to the Orion and the Perseus arm regions. Strong photospheric line profile variability is present on a night to night timescale, with He I 5875 A displaying an additional blue-shifted absorption in some of the spectra, similar to what is observed in the optical counterpart of Vela X-1. A connection between the extended blue wing and X-ray maximum is suggested. Short timescale variations in line profiles are detected on only two nights, but the evidence that these variations occur on the 2.78 hour X-ray flaring period is marginal.

Numerical simulations of the interaction between supernova ejecta and a stellar wind are presented. We follow the temporal evolution of the shock fronts that are formed through such an interaction and determine the velocities, temperatures and densities. We model the X-ray emission from the SNR-stellar wind collision region and we compare it with recent results from X-ray observations carried out with the Chandra satellite of the SMC supernova remnant SNR 0057-7226 which could be interacting with the wind of the Wolf-Rayet system HD 5980. The simulations predict the presence of shell-like regions of enhanced X-ray emission which are consistent with the presence of X-ray emitting arcs in the Chandra image. Also the observed X-ray luminosity is comparable to the X-ray luminosities we obtain from the simulations for a supernova with an initial energy in the (1-5)E50 erg range.

We present results from a Chandra observation of the NGC346 cluster. This cluster contains numerous massive stars and is responsible for the ionization of N66, the most luminous HII region and the largest star formation region in the SMC. In this first paper, we will focus on the characteristics of the main objects of the field. The NGC346 cluster itself shows only relatively faint X-ray emission (with L_X^{unabs} ~ 1.5 times 10^{34} erg s^{-1}), tightly correlated with the core of the cluster. In the field also lies HD5980, a LBV star in a binary (or possibly a triple system) that is detected for the first time at X-ray energies. The star is X-ray bright, with an unabsorbed luminosity of L_X^{unabs} ~ 1.7 times 10^{34} erg s^{-1}, but needs to be monitored further to investigate its X-ray variability over a complete 19d orbital cycle. The high X-ray luminosity may be associated either with colliding winds in the binary system or with the 1994 eruption. HD5980 is surrounded by a region of diffuse X-ray emission, which is a supernova remnant. While it may be only a chance alignment with HD5980, such a spatial coincidence may indicate that the remnant is indeed related to this peculiar massive star.

We present the results of an intensive campaign of spectroscopic and photometric monitoring of the peculiar Wolf-Rayet star WR 134 from 1989 to 1997. This unprecedentedly large data set allows us to confirm unambiguously the existence of a coherent 2.25 +/- 0.05 day periodicity in the line-profile changes of He II 4686, although the global pattern of variability is different from one epoch to another. This period is only marginally detected in the photometric data set. Assuming the 2.25 day periodic variability to be induced by orbital motion of a collapsed companion, we develop a simple model aiming at investigating (i) the effect of this strongly ionizing, accreting companion on the Wolf-Rayet wind structure, and (ii) the expected emergent X-ray luminosity. We argue that the predicted and observed X-ray fluxes can only be matched if the accretion on the collapsed star is significantly inhibited. Additionally, we performed simulations of line-profile variations caused by the orbital revolution of a localized, strongly ionized wind cavity surrounding the X-ray source. A reasonable fit is achieved between the observed and modeled phase-dependent line profiles of He II 4686. However, the derived size of the photoionized zone substantially exceeds our expectations, given the observed low-level X-ray flux. Alternatively, we explore rotational modulation of a persistent, largely anisotropic outflow as the origin of the observed cyclical variability. Although qualitative, this hypothesis leads to greater consistency with the observations.