We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of the gamma-ray binary system containing the pulsar PSR B1259-63 orbiting around a massive star LS 2883 in an inactive period between the 2017 and 2021 periastron passages. We detected radio continuum emission from the binary system at 97 GHz (Band 3) and 343 GHz (Band 7). Compared with our previous ALMA observations performed soon after the 2017 periastron passage, the fluxes have decreased by an factor of six at 97 GHz and two at 343 GHz. The flux at 343 GHz is large relative to that at 97 GHz and appears to be thermal emission from the circumstellar disk around LS 2883. The decrease of the 343 GHz flux may indicate that the disk has expanded and become partially optically-thin since the disk is no longer affected by pulsar winds. The flux at 97 GHz is consistent with that expected from the pulsed emission from the pulsar, which indicates that the unpulsed emission that had been produced through pulsar-disk or pulsar-stellar wind interaction has disappeared. The image of the system is consistent with a point source and shows no sign of ejecta.

As the largest population of high mass X-ray binaries, Be/X-ray binaries provide an excellent laboratory to investigate the extreme physics of neutron stars. It is generally accepted that Be stars possess a circumstellar disc, providing an additional source of accretion to the stellar winds present around young hot stars. Interaction between the neutron star and the disc is often the dominant accretion mechanism. A large amount of work has gone into modelling the properties of these circumstellar discs, allowing for the explanation of a number of observable phenomena. In this paper, smoothed particle hydroynamics simulations are performed whilst varying the model parameters (orbital period, eccentricity, the mass ejection rate of the Be star and the viscosity and orientation of the disc). The relationships between the model parameters and the disc's characteristics (base gas density, the accretion rate of the neutron star and the disc's size) are presented. The observational evidence for a dependency of the size of the Be star's circumstellar disc on the orbital period (and semi-major axis) is supported by the simulations.

The tidal interaction of a Be star with a binary companion forms two spiral arms that cause orbital modulation of the Be disc structure. The aim of this work is to identify observables in which this modulation is apparent. The structure of a Be disc in a coplanar circular binary system is computed with a smoothed-particle hydrodynamics code, and a radiation transfer code calculates the spectral energy distribution. Line depolarisation was confirmed, with polarisation profiles nearly reverse to emission-line profiles. The continuum flux maximizes for pole-on discs, but photometric variability maximizes for edge-on discs. The linear polarisation exhibits one or two maxima per orbital cycle. While polarisation variability in visible passbands is important only at low inclinations, infrared bands may demonstrate high orbital variability even at large inclinations. More evident is the modulation in the polarisation angle (PA) for low inclinations. The latter can be used to track azimuthal asymmetries for pole-on discs, where the spectroscopic variability in the violet-to-red (V/R) emission-component ratio disappears. PA reversals coincide with phases where V/R=1, tracking lines of sight directed towards regions where the approaching and receding arms overlap. Continuum flux and polarisation are mostly in phase for neighbouring wavelength regions. It is suggested that studies of non-symmetric discs distorted by tidal forces from a secondary star may be used to study disc variabilities of other origins.

SXP 5.05 is a Be/X-ray binary with a neutron star companion located in the Small Magellanic Cloud. It was first detected in 2013, and later that year, SXP 5.05 underwent a massive optical and X-ray outburst. This outburst dwarfs any other optical event that has been observed for this system during the last 5 years. The large increase in optical brightness of the system implies an increase in the size and density of the Be star's circumstellar disc. The X-ray data show two occultations of the neutron star per orbit and is not consistent with a neutron star passing behind the Be star, and hence the disc is responsible for these occultations. In this paper, we model the outburst of Be/neutron star binary SXP 5.05 as being due to a large increase in mass ejection by the Be star. The neutron star passes directly through the growing disc, and it is shown that the resulting obscuration can qualitatively explain the observed X-ray behaviour of the system. We find the only way to reproduce the timescales of the observed optical behaviour of the system is to increase the mass ejection substantially for a short time (<2 orbits) and to decrease the viscosity during the event. The general behaviour of the observed X-ray and Ha line emission are also reproduced by the model. However, the inferred mass ejection and viscosity needed to produce a sufficiently rapid increase of disc size are both higher than suggested by previous works.

X-ray photometry and optical spectra are presented covering the periastron passage of the highly-eccentric, ~50 year binary system PSR J2032+4127 in November 2017. This system consists of a 143 ms pulsar in orbit around a massive OB star, MT 91-213. The data show dramatic changes during the encounter as the pulsar wind collided with the stellar wind. The X-ray flux rose on the approach to periastron, then underwent a major dip in the few days around periastron, and then gradually declined over the next few weeks. The optical spectroscopy revealed a steady decline in the H? line strength on the approach to periastron (from an Equivalent Width of -15A to -7A) implying a truncation of the OB star's circumstellar disk by the approaching neutron star. Smooth Particle Hydrodynamic (SPH) modelling is used here to model the system within the context of the observed behaviour and predict the geometrical configuration of the circumstellar disk with respect to the pulsar's orbit.

There are over 100 Be stars that are known to have neutron star companions but only one such system with a black hole. Previous theoretical work suggests this is not due to their formation but due to differences in X-ray luminosity. It has also been proposed that the truncation of the Be star's circumstellar disc is dependent on the mass of the compact object. Hence, Be star discs in black hole binaries are smaller. Since accretion onto the compact object from the Be star's disc is what powers the X-ray luminosity, a smaller disc in black hole systems leads to a lower luminosity. In this paper, simulations are performed with a range of eccentricities and compact object mass. The disc's size and density are shown to be dependent on both quantities. Mass capture and, in turn, X-ray luminosity are heavily dependent on the size and density of the disc. Be/black hole binaries are expected to be up to ~10 times fainter than Be/neutron star binaries when both systems have the same eccentricity and can be 100 times fainter when comparing systems with different eccentricity.

The first results of radiative transfer calculations on decretion discs of binary Be stars are presented. A smoothed particle hydrodynamics code computes the structure of Be discs in coplanar circular binary systems for a range of orbital and disc parameters. The resulting disc configuration consists of two spiral arms, and can be given as input into a Monte Carlo code, which calculates the radiative transfer along the line of sight for various observational coordinates. Making use of the property of steady disc structure in coplanar circular binaries, observables are computed as functions of the orbital phase. Orbital-phase series of line profiles are given for selected parameter sets under various viewing angles, to allow comparison with observations. Flat-topped profiles with and without superimposed multiple structures are reproduced, showing, for example, that triple-peaked profiles do not have to be necessarily associated with warped discs and misaligned binaries. It is demonstrated that binary tidal effects give rise to phase-locked variability of the violet-to-red (V/R) ratio of hydrogen emission lines. The V/R ratio exhibits two maxima per cycle; in certain cases those maxima are equal, leading to a clear new V/R cycle every half orbital period. This study opens a way in identifying binaries and in constraining the parameters of binary systems that exhibit phase-locked variations induced by tidal interaction with a companion star.

We report on the RXTE detection of a sudden increase in the absorption column density, $N_\mathrm{H}$, during the 2011 May outburst of GX 304-1. The $N_\mathrm{H}$ increased up to ${\sim}16\times 10^{22}$ atoms cm$^{-2}$, which is a factor of 3-4 larger than what is usually measured during the outbursts of GX 304-1 as covered by RXTE. Additionally, an increase in the variability of the hardness ratio as calculated from the energy resolved RXTE-PCA light curves is measured during this time range. We interpret these facts as an occultation event of the neutron star by material in the line of sight. Using a simple 3D model of an inclined and precessing Be disk around the Be type companion, we are able to qualitatively explain the $N_\mathrm{H}$ evolution over time. We are able to constrain the Be-disk density to be on the order of $10^{-11}$ g cm$^{-3}$. Our model strengthens the idea of inclined Be disks as origin of double-peaked outbursts as the derived geometry allows accretion twice per orbit under certain conditions.

We use a smoothed particle hydrodynamics (SPH) code to examine the effects of misaligned binary companions on Be star discs. We systematically vary the degree of misalignment between the disc and the binary orbit, as well as the disc viscosity and orbital period to study their effects on the density in the inner and outer parts of the disc. We find that varying the degree of misalignment, the viscosity, and the orbital period affects both the truncation radius and the density structure of the outer disc, while the inner disc remains mostly unaffected. We also investigate the tilting of the disc in the innermost part of the disc and find the tilt increases with radius until reaching a maximum around 5 stellar radii. The direction of the line of nodes, with respect to the equator of the central star, is found to be offset compared to the orbital line of nodes, and to vary periodically in time, with a period of half a orbital phase. We also compare the scale height of our discs with the analytical scale height of an isothermal disc, which increases with radius as $r^{1.5}$. We find that this formula reproduces the scale height well for both aligned and misaligned systems but underestimates the scale height in regions of the disc where density enhancements develop.

In this work we explore the effect of binarity in the decretion disc of Be stars, in order to explain their variability. To this aim, we performed smoothed particle hydrodynamics (SPH) simulations on Be binary systems, following the matter ejected isotropically from the equator of the Be star towards the base of an isothermal decretion disc. We let the system evolve for time long enough to be considered at steady state, and focus on the effect of viscosity for coplanar prograde binary orbits. The disc structure is found to be locked to the orbital phase, exhibiting also a dependence on the azimuthal angle. Additionally, we present the first results from detailed non-local thermodynamic equilibrium (non-LTE) radiative transfer calculations of the disc structure computed with the SPH code. This is achieved by the use of the three-dimensional (3D) Monte Carlo code HDUST, which can produce predictions with respect to a series of observables.

A global disk oscillation implemented in the viscous decretion disk (VDD) model has been used to reproduce most of the observed properties of the well known Be star $\zeta$ Tau. 48 Librae shares several similarities with $\zeta$ Tau -- they are both early-type Be stars, they display shell characteristics in their spectra, and they exhibit cyclic $V/R$ variations -- but has some marked differences as well, such as a much denser and more extended disk, a much longer $V/R$ cycle, and the absence of the so-called triple-peak features. We aim to reproduce the photometric, polarimetric, and spectroscopic observables of 48 Librae with a self-consistent model, and to test the global oscillation scenario for this target. Our calculations are carried out with the three-dimensional NLTE radiative transfer code HDUST. We employ a rotationally deformed, gravity-darkened central star, surrounded by a disk whose unperturbed state is given by the VDD model. A two-dimensional global oscillation code is then used to calculate the disk perturbation, and superimpose it on the unperturbed disk. A very good, self-consistent fit to the time-averaged properties of the disk is obtained with the VDD. The calculated perturbation has a period $P = 12$ yr, which agrees with the observed period, and the behaviour of the $V/R$ cycle is well reproduced by the perturbed model. The perturbed model improves the fit to the photometric data and reproduces some features of the observed spectroscopic data. Some suggestions to improve the synthesized spectroscopy in a future work are given.

Be stars are surrounded by outflowing circumstellar matter structured in the form of decretion discs. They are often members of binary systems, where it is expected that the decretion disc interacts both radiatively and gravitationally with the companion. In this work we study how various orbital (period, mass ratio and eccentricity) and disc (viscosity) parameters affect the disc structure in coplanar systems. We simulate such binaries with the use of a smoothed particle hydrodynamics code. The main effects of the secondary on the disc are its truncation and the accumulation of material inwards of truncation. We find two cases with respect to the effects of eccentricity: (i) In circular or nearly circular prograde orbits, the disc maintains a rotating, constant in shape, configuration, which is locked to the orbital phase. The disc is smaller in size, more elongated and more massive for low viscosity parameter, small orbital separation and/or high mass ratio. (ii) Highly eccentric orbits are more complex, with the disc structure and total mass strongly dependent on the orbital phase and the distance to the secondary. We also study the effects of binarity in the disc continuum emission. Since the infrared and radio SED are sensitive to the disc size and density slope, the truncation and matter accumulation result in considerable modifications in the emergent spectrum. We conclude that binarity can serve as an explanation for the variability exhibited in observations of Be stars, and that our model can be used to detect invisible companions.

We study the effect of density distribution evolution on the global one-armed oscillation modes in low viscosity disks around isolated and binary Be stars. Observations show that some Be stars exhibit evidence of formation and dissipation of the equatorial disk. In this paper, we first calculate the density evolution in disks around isolated Be stars. To model the formation stage of the disk, we inject mass at a radius just outside the star at a constant rate for 30-50 years. As the disk develops, the density distribution approaches the form of the steady disk solution. Then, we turn off the mass injection to model the disk dissipation stage. The innermost part of the disk starts accretion, and a gap forms between the star and the disk. Next, we calculate the one-armed modes at several epochs. We neglect the effect of viscosity because the time-scale of oscillations is much shorter than the disk evolution time-scale for low viscosity. In the disk formation stage, the eigenfrequency increases with time toward the value for the steady state disk. On the other hand, one-armed eigenmodes in dissipating Be disks have significantly higher eigenfrequencies and narrower propagation regions. Observationally, such a change of mode characteristics can be taken as an evidence for gap opening around the star. In binary Be stars, the characteristics of the disk evolution and the eigenmodes are qualitatively the same as in isolated Be stars, but quantitatively they have shorter evolution time-scales and higher eigenfrequencies, which is in agreement with the observed trend.

We have used the radiative transfer code HDUST to analyze and interpret the long-term photometric behavior of the Be star omega CMa, considering four complete cycles of disk formation and dissipation. This is the first time in which a full lightcurve of a Be star was investigated and modeled including both disk build-up and dissipation phases. Based on the quite good fit of the observed data we were able to derive the history of stellar mass decretion rates (including long- and short-term changes) during the disk formation and dissipation phases in all four cycles.

We report on monitoring observations of the TeV gamma-ray binary HESS J0632+057, which were carried out to constrain the interaction between the Be circumstellar disk and the compact object of unknown nature, and provide for the first time high-dispersion (R > 50000) optical spectra in the second half of the orbital cycle, from apastron through periastron. The Halpha, Hbeta, and Hgamma line profiles are found to exhibit remarkable short-term variability for ~1 month after the apastron (phase 0.6--0.7), whereas they show little variation near the periastron. These emission lines show "S-shaped" variations with timescale of ~150 days, which is about twice that reported previously. In contrast to the Balmer lines, no profile variability is seen in any FeII emission line. We estimate the radii of emitting regions of the Halpha, Hbeta, Hgamma, and FeII emission lines to be ~30, 11, 7, and 2 stellar radii (R_*), respectively. The amplitudes of the line profile variations in different lines indicate that the interaction with the compact object affects the Be disk down to, at least, the radius of 7 R_* after the apastron. This fact, together with little profile variability near the periastron, rules out the tidal force as the major cause of disk variability. Although this leaves the pulsar wind as the most likely candidate mechanism for disk variations, understanding the details of the interaction, particularly the mechanism for causing a large-scale disk disturbance after the apastron, remains an open question.

Context. Various Be stars exhibit intensity variations of the violet and red emission peaks in their HI lines observed in emission. This so-called $V/R$ phenomenon is usually explained by the precession of a one-armed spiral density perturbation in the circumstellar disk. That global-disk oscillation scenario was confirmed, both observationally and theoretically, in the previous series of two papers analyzing the Be shell star {\zeta} Tauri. The vertically averaged (2D) global-disk oscillation model used at the time was able to reproduce the $V/R$ variations observed in H{\alpha}, as well as the spatially resolved interferometric data from AMBER/VLTI. Unfortunately, that model failed to reproduce the $V/R$ phase of Br15 and the amplitude of the polarization variation, suggesting that the inner disk structure predicted by the model was incorrect. Aims. The first aim of the present paper is to quantify the temporal variations of the shell-line characteristics of {\zeta} Tauri. The second aim is to better understand the physics underlying the $V/R$ phenomenon by modeling the shell-line variations together with the $V/R$ and polarimetric variations. The third aim is to test a new 2.5D disk oscillation model, which solves the set of equations that describe the 3D perturbed disk structure but keeps only the equatorial (i.e., 2D) component of the solution. This approximation was adopted to allow comparisons with the previous 2D model, and as a first step toward a future 3D model. Results. The new 2.5D formalism improves the agreement with the observed $V/R$ variations of H{\alpha} and Br15, under the proviso that a large value of the viscosity parameter, {\alpha} = 0.8, be adopted. Nonetheless, it remains challenging for the models to reproduce consistently the amplitude and the average level of the polarization data, whatever formalism is adopted.

We study the warping and tearing of a geometrically thin, non-self-gravitating disk surrounding binary supermassive black holes on an eccentric orbit. The circumbinary disk is significantly misaligned with the binary orbital plane, and is subject to the time-dependent tidal torques. In principle, such a disk is warped and precesses, and is torn into mutually misaligned rings in the region, where the tidal precession torques are stronger than the local viscous torques. We derive the tidal-warp and tearing radii of the misaligned circumbinary disks around eccentric SMBH binaries. We find that in disks with the viscosity parameter, alpha, larger than a critical value depending on the disk aspect ratio, the disk warping appears outside the tearing radius. This condition is expressed as alpha > sqrt{H/3r} for H/r ~<0.1, where H is the disk scale height. If alpha < sqrt{H/3r}, only the disk tearing occurs because the tidal warp radius is inside the tearing radius, where most of disk material is likely to rapidly accrete onto SMBHs. In warped and torn disks, both the tidal-warp and the tearing radii most strongly depend on the binary semi-major axis, although they also mildly depend on the other orbital and disk parameters. This strong dependence enables us to estimate the semi-major axis, once the tidal warp or tearing radius is determined observationally: For the tidal warp radius of 0.1 pc, the semi-major axis is estimated to be ~10^{-2} pc for 10^7 Msun black hole with typical orbital and disk parameters. We also briefly discuss the possibility that central objects of observed warped maser disks in active galactic nuclei are supermassive black hole binaries.

The binary system PSR B1259-63/LS 2883 is well sampled in radio, X-rays, and TeV gamma-rays, and shows orbital phase-dependent variability in these frequencies. The first detection of GeV gamma-rays from the system was made around the 2010 periastron passage. In this Letter, we present an analysis of X-ray and gamma-ray data obtained by the Swift/XRT, NuSTAR/FPM, and Fermi/LAT, through the recent periastron passage which occurred on 2014 May 4. While PSR B1259-63/LS 2883 was not detected by the LAT before and during this passage, we show that the GeV flares occurred at a similar orbital phase as in early 2011, thus establishing the repetitive nature of the post-periastron GeV flares. Multiple flares each lasting for a few days have been observed and short-term variability is seen as well. We also found X-ray flux variation contemporaneous with the GeV flare for the first time. A strong evidence of the keV-to-GeV connection came from the broadband high-energy spectra, which we interpret as synchrotron radiation from the shocked pulsar wind.

We study a warping instability of a geometrically thin, non-self-gravitating, circumbinary disk around young binary stars on an eccentric orbit. Such a disk is subject to both the tidal torques due to a time-dependent binary potential and the radiative torques due to radiation emitted from each star. The tilt angle between the circumbinary disk plane and the binary orbital plane is assumed to be very small. We find that there is a radius within/beyond which the circumbinary disk is unstable to radiation-driven warping, depending on the disk density and temperature gradient indices. This marginally stable warping radius is very sensitive to viscosity parameters, a fiducial disk radius and the temperature measured there, the stellar luminosity, and the disk surface density at a radius where the disk changes from the optically thick to thin for the irradiation from the central stars. On the other hand, it is insensitive to the orbital eccentricity and binary irradiation parameter, which is a function of the binary mass ratio and luminosity of each star. Since the tidal torques can suppress the warping in the inner part of the circumbinary disk, the disk starts to be warped in the outer part. While the circumbinary disks are most likely to be subject to the radiation-driven warping on a AU to kilo-AU scale for binaries with young massive stars more luminous than 10^4Lsun, the radiation driven warping does not work for those around young binaries with the luminosity comparable to the solar luminosity.

We study a warping instability of a geometrically thin, non-self-gravitating disk surrounding binary supermassive black holes on a circular orbit. Such a circumbinary disk is subject to not only tidal torques due to the binary gravitational potential but also radiative torques due to radiation emitted from an accretion disk around each black hole. We find that a circumbinary disk initially aligned with the binary orbital plane is unstable to radiation-driven warping beyond the marginally stable warping radius, which is sensitive to both the ratio of vertical to horizontal shear viscosities and the mass-to-energy conversion efficiency. As expected, the tidal torques give no contribution to the growth of warping modes but tend to align the circumbinary disk with the orbital plane. Since the tidal torques can suppress the warping modes in the inner part of circumbinary disk, the circumbinary disk starts to be warped at radii larger than the marginally stable warping radius. If the warping radius is of the order of 0.1 pc, a resultant semi-major axis is estimated to be of the order of 10^-2 pc to 10^-4 pc for 10^7 Msun black hole. We also discuss the possibility that the central objects of observed warped maser disks in active galactic nuclei are binary supermassive black holes with a triple disk: two accretion disks around the individual black holes and one circumbinary disk surrounding them.

MAXI/GSC observed 21 outbursts from Circinus X-1 between 2009 August and 2013 December. Although 14 outbursts showed ordinary gradual decays, in 7 outbursts we found sudden luminosity decrease in a time scale of a few hours around the periastron, and then the outbursts terminated. These sudden decreases started at the estimated luminosity of a few times $10^{37}$ erg s$^{-1}$ and reached to $\lesssim3\times10^{36}$ erg s$^{-1}$. We propose three interpretations for the sudden luminosity decrease: (1) the end of the outburst during the dip, (2) the propeller effect, and (3) the stripping effect by the stellar wind of the companion star. It is difficult to explain the phenomenon with any of these interpretations alone. The interpretation of (1) is possible for only two outbursts assuming rapid decay. The propeller effect (2) is expected to occur at a constant luminosity, which is incompatible with the observed facts. In wind stripping effect (3), the ram pressure of a typical stellar wind is not sufficient to blow out most of the accretion disk. In this paper, we discuss a possibility of a modified effect of (3) assuming other additional conditions such as wind clumping and disk instability.

Cyg X-3 is a high mass X-ray binary consisting of a Wolf-Rayet star and a compact object in a very short orbital period of 4.8h. The only confirmed microquasar with high energy gamma-ray emission, Cyg X-3 provides a unique opportunity to study the relationship between the accretion power and the power in high energy emission. Because of a compact orbit and a slow Wolf-Rayet wind, the flow structure around the compact object is thought to be strongly affected by the orbital motion, details of which can be obtained only by numerical simulations. In this paper, we report on the results from 3D hydrodynamic simulations of the wind accretion in Cyg X-3. For simplicity we adopt an anti-gravity-like force that emulates the radiative acceleration consistent with the beta-velocity wind. Due to the rapid orbital motion, the flow around the compact object has large density gradients. As a result, the accretion rate onto the compact object is significantly lower than that of the Bondi-Hoyle-Lyttleton rate. We also calculate the model X-ray light curve. Although it roughly agrees with the observed light curve, more detailed modeling is needed for detailed comparison.

Colliding wind binaries (CWBs) are unique laboratories for X-ray astrophysics. The massive stars in these systems possess powerful stellar winds with speeds up to $\sim$3000 km s$^{-1}$, and their collision leads to hot plasma (up to $\sim10^8$K) that emit thermal X-rays (up to $\sim$10 keV). Many X-ray telescopes have observed CWBs, including Suzaku, and our work aims to model these X-ray observations. We use 3D smoothed particle hydrodynamics (SPH) to model the wind-wind interaction, and then perform 3D radiative transfer to compute the emergent X-ray flux, which is folded through X-ray telescopes' response functions to compare directly with observations. In these proceedings, we present our models of Suzaku observations of the multi-year-period, highly eccentric systems $\eta$ Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only $\eta$ Carinae's X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.

Recent work suggests that the mass-loss rate of the primary star (Eta A) in the massive colliding wind binary Eta Carinae dropped by a factor of 2-3 between 1999 and 2010. We present results from large- (r=1545au) and small- (r=155au) domain, 3D smoothed particle hydrodynamic (SPH) simulations of Eta Car's colliding winds for 3 Eta A mass-loss rates (2.4, 4.8, and 8.5 x 10^-4 M_sun/yr), investigating the effects on the dynamics of the binary wind-wind collision (WWC). These simulations include orbital motion, optically thin radiative cooling, and radiative forces. We find that Eta A's mass-loss rate greatly affects the time-dependent hydrodynamics at all spatial scales investigated. The simulations also show that the post-shock wind of the companion star (Eta B) switches from the adiabatic to the radiative-cooling regime during periastron passage. The SPH simulations together with 1D radiative transfer models of Eta A's spectra reveal that a factor of 2 or more drop in Eta A's mass-loss rate should lead to substantial changes in numerous multiwavelength observables. Recent observations are not fully consistent with the model predictions, indicating that any drop in Eta A's mass-loss rate was likely by a factor < 2 and occurred after 2004. We speculate that most of the recent observed changes in Eta Car are due to a small increase in the WWC opening angle that produces significant effects because our line-of-sight to the system lies close to the dense walls of the WWC zone. A modest decrease in Eta A's mass-loss rate may be responsible, but changes in the wind/stellar parameters of Eta B cannot yet be fully ruled out. We suggest observations during Eta Car's next periastron in 2014 to further test for decreases in Eta A's mass-loss rate. If Eta A's mass-loss rate is declining and continues to do so, the 2014 X-ray minimum should be even shorter than that of 2009.

We carried out optical high-dispersion spectroscopic monitoring of the Be disk in a Be/X-ray binary A 0535+262/V725 Tau from 2009 to 2012, covering two giant outbursts and several normal outbursts. This monitoring was performed in order to investigate variabilities of the Be disk due to the interaction with the neutron star in recent X-ray active phase from 2008 to 2011. Such variabilities give a clue to uncleared detailed mechanism for very bright X-ray outbursts, which are unique to some Be/X-ray binaries with relatively wide and eccentric orbit. In the previous letter (Moritani et al. 2011), a brief overview of line profile variabilities around the 2009 giant outburst was given and the possibility of a warped Be disk was discussed. In this paper, a full analysis of the Halpha line profiles as well as other line profiles is carried out. A bright blue component, or blue "shoulder", showing up after periastron indicates the presence of a dense gas stream toward the neutron star, which is associated with observed outbursts. We re-analyze the Halpha line profiles before 2009 (down to 2005) in order to investigate the variability of the the disk structure in the innermost region, which seems to have detached from the Be star surface by 2008. A redshifted enhanced component is remarkable in all emission lines observed around the 2009 giant outburst, occasionally forming a triple peak. These features indicate that the Be disk was warped in X-ray active phase. We estimate the position of the warped region from fitting the radial velocity of the redshifted enhanced component of Halpha, and find that it was very close to the periastron when two giant outbursts in 2009 and 2011 and a bright normal outburst in 2010 March occurred. These facts strongly suggest that the warped Be disk triggered these giant outbursts.

We present the results of the optical spectropolarimetry of Earthshine on the Moon for Earth phase angles ranging from 49 to 96 degrees. The observations were conducted on 2011 March 9-13 (UT) using the spectropolarimeter HBS installed on the 1.88 m telescope at Okayama Astrophysical Observatory. The wavelength coverage was 450-850 nm with a resolution of 6 nm. The observed Earthshine polarization degree spectra exhibit decreasing polarization degree with increasing wavelength at any phase. The overall degree of polarization increases as the Earth approaches a quadrature phase. The phase dependence differs with the wavelengths; the maximum polarization for the V band occurs at a phase angle of ~90 degrees, whereas that for longer wavelengths is reached at larger phase angles. This is interpreted as indicating that Earthshine polarization at shorter wavelengths is dominated by atmospheric Rayleigh scattering, whereas that at longer wavelengths has an increasingly effective contribution from the Earth surface reflection. The wavelength dependence of the phase angle of the maximum polarization appears to be unique among the terrestrial planetary bodies in the Solar System. Therefore this might constitute important evidence pointing toward a distinctive characteristic of the Earth: the planet has a scattering but transparent atmosphere above its surface.

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.

We propose the new scenario for X-ray outbursts in Be/X-ray binaries that normal and giant outbursts are respectively caused by radiatively inefficient accretion flows (RIAFs) and Bondi-Hoyle-Lyttleton (BHL) accretion of the material transferred from the outermost part of a Be disk misaligned with the binary orbital plane. Based on simulated mass-transfer rates from misaligned Be disks, together with simplified accretion flow models, we show that mass-accretion rates estimated from the luminosity of the normal X-ray outbursts are consistent with those obtained with advection-dominated accretion flows, not with the standard, radiative-cooling dominated, accretion. Our RIAF scenario for normal X-ray outbursts resolves problems that have challenged the standard disk picture for these outbursts. When a misaligned Be disk crosses the orbit of the neutron star, e.g., by warping, the neutron star can capture a large amount of mass via BHL-type accretion during the disk transit event. We numerically show that such a process can reproduce the X-ray luminosity of giant X-ray outbursts. In the case of very high Be disk density, the accretion flow associated with the disk transit becomes supercritical, giving rise to the luminosity higher than the Eddington luminosity.

Be stars possess gaseous circumstellar disks that modify in many ways the spectrum of the central B star. Furthermore, they exhibit variability at several timescales and for a large number of observables. Putting the pieces together of this dynamical behavior is not an easy task and requires a detailed understanding of the physical processes that control the temporal evolution of the observables. There is an increasing body of evidence that suggests that Be disks are well described by standard $\alpha$-disk theory. This paper is the first of a series that aims at studying the possibility of inferring several disk and stellar parameters through the follow-up of various observables. Here we study the temporal evolution of the disk density for different dynamical scenarios, including the disk build-up as a result of a long and steady mass injection from the star, the disk dissipation that occurs after mass injection is turned off, as well as scenarios in which active periods are followed by periods of quiescence. For those scenarios, we investigate the temporal evolution of continuum photometric observables using a 3-D non-LTE radiative transfer code. We show that lightcurves for different wavelengths are specific of a mass loss history, inclination angle and $\alpha$ viscosity parameter. The diagnostic potential of those lightcurves is also discussed.

Temporal changes of X-ray to very-high-energy gamma-ray emissions from the pulsar-Be star binary PSR B1259-63/LS 2883 are studied based on 3-D SPH simulations of pulsar wind interaction with Be-disk and wind. We focus on the periastron passage of the binary and calculate the variation of the synchrotron and inverse-Compton emissions using the simulated shock geometry and pressure distribution of the pulsar wind. The characteristic double-peaked X-ray light curve from observations is reproduced by our simulation under a dense Be disk condition (base density ~10^{-9} g cm^{-3}). We interpret the pre- and post-periastron peaks as being due to a significant increase in the conversion efficiency from pulsar spin down power to the shock-accelerated particle energy at orbital phases when the pulsar crosses the disk before periastron passage, and when the pulsar wind creates a cavity in the disk gas after periastron passage, respectively. On the contrary, in the model TeV light curve, which also shows a double peak feature, the first peak appears around the periastron phase. The possible effects of cooling processes on the TeV light curve are briefly discussed.

Be stars possess gaseous circumstellar decretion disks, which are well described using standard $\alpha$-disk theory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous decretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter $\alpha=1.0\pm0.2$, corresponding to a mass injection rate $\dot{M}=(3.5\pm 1.3) \times 10^{-8}\ M_\sun\,\mathrm{yr}^{-1}$. Such a large value of $\alpha$ suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass loss rate inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, but some other mechanism.

We present a three-dimensional (3-D) dynamical model for the broad [Fe III] emission observed in Eta Carinae using the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS). This model is based on full 3-D Smoothed Particle Hydrodynamics (SPH) simulations of Eta Car's binary colliding winds. Radiative transfer codes are used to generate synthetic spectro-images of [Fe III] emission line structures at various observed orbital phases and STIS slit position angles (PAs). Through a parameter study that varies the orbital inclination i, the PA {\theta} that the orbital plane projection of the line-of-sight makes with the apastron side of the semi-major axis, and the PA on the sky of the orbital axis, we are able, for the first time, to tightly constrain the absolute 3-D orientation of the binary orbit. To simultaneously reproduce the blue-shifted emission arcs observed at orbital phase 0.976, STIS slit PA = +38 degrees, and the temporal variations in emission seen at negative slit PAs, the binary needs to have an i \approx 130 to 145 degrees, {\theta} \approx -15 to +30 degrees, and an orbital axis projected on the sky at a PA \approx 302 to 327 degrees east of north. This represents a system with an orbital axis that is closely aligned with the inferred polar axis of the Homunculus nebula, in 3-D. The companion star, Eta B, thus orbits clockwise on the sky and is on the observer's side of the system at apastron. This orientation has important implications for theories for the formation of the Homunculus and helps lay the groundwork for orbital modeling to determine the stellar masses.

Eta Carinae is considered to be a massive colliding wind binary system with a highly eccentric (e \sim 0.9), 5.54-yr orbit. However, the companion star continues to evade direct detection as the primary dwarfs its emission at most wavelengths. Using three-dimensional (3-D) SPH simulations of Eta Car's colliding winds and radiative transfer codes, we are able to compute synthetic observables across multiple wavebands for comparison to the observations. The models show that the presence of a companion star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line profiles, as well as the ground-based photometric monitoring. Here, we focus on the Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in the dense primary wind, allowing increased escape of radiation from the hotter/deeper layers of the primary's extended wind photosphere. The results have important implications for interpretations of Eta Car's observables at multiple wavelengths.

The colliding wind binary (CWB) systems \eta\ Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these X-ray light curves and spectra, we apply 3D hydrodynamic simulations of the wind-wind collision using smoothed particle hydrodynamics (SPH), with the recent improvements of radiative cooling and the acceleration of the stellar winds according to a \beta\ law. For both systems, the 2-10 keV RXTE light curves are well-reproduced in absolute units for most phases, but the light curve dips associated with the periastron passages are not well matched. In WR140, the dip is too weak, and in \eta\ Carinae, the large difference in wind speeds of the two stars leads to a hot, post-periastron bubble that produces excess emission toward the end of the X-ray minimum.

The colliding wind binary (CWB) systems \eta\ Carinae and WR140 provide unique laboratories for X-ray astrophysics. Their wind-wind collisions produce hard X-rays that have been monitored extensively by several X-ray telescopes, including RXTE. To interpret these RXTE X-ray light curves, we model the wind-wind collision using 3D smoothed particle hydrodynamics (SPH) simulations. Adiabatic simulations that account for the absorption of X-rays from an assumed point source at the apex of the wind-collision shock cone by the distorted winds can closely match the observed 2-10keV RXTE light curves of both \eta\ Car and WR140. This point-source model can also explain the early recovery of \eta\ Car's X-ray light curve from the 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of \eta\ Car. Our more recent models relax the point-source approximation and account for the spatially extended emission along the wind-wind interaction shock front. For WR140, the computed X-ray light curve again matches the RXTE observations quite well. But for \eta\ Car, a hot, post-periastron bubble leads to an emission level that does not match the extended X-ray minimum observed by RXTE. Initial results from incorporating radiative cooling and radiatively-driven wind acceleration via a new anti-gravity approach into the SPH code are also discussed.

We report on high-dispersion optical spectroscopic observations of the Be/X-ray binary A0535+262/V725 Tau during the giant outburst in November/December 2009 and after it. The observed emission line profiles, reflecting the structure of the geometrically thin circumstellar envelope of the Be star (Be disk), show drastic variabilities and indicate the existence of a warped component. The enhanced blue shoulder seen after periastron passage implies the gas stream from a dense part of the Be disk to the neutron star.

We study the interaction between the Be star and the pulsar in the TeV binary PSR B1259-63/LS 2883, using 3-D SPH simulations of the tidal and wind interactions in this Be-pulsar system. We first run a simulation without pulsar wind nor Be wind, taking into account only the gravitational effect of the pulsar on the Be disk. In this simulation, the gas particles are ejected at a constant rate from the equatorial surface of the Be star, which is tilted in a direction consistent with multi-waveband observations. We run the simulation until the Be disk is fully developed and starts to repeat a regular tidal interaction with the pulsar. Then, we turn on the pulsar wind and the Be wind. We run two simulations with different wind mass-loss rates for the Be star, one for a B2V type and the other for a significantly earlier spectral type. Although the global shape of the interaction surface between the pulsar wind and the Be wind agrees with the analytical solution, the effect of the pulsar wind on the Be disk is profound. The pulsar wind strips off an outer part of the Be disk, truncating the disk at a radius significantly smaller than the pulsar orbit. Our results, therefore, rule out the idea that the pulsar passes through the Be disk around periastron, which has been assumed in the previous studies. It also turns out that the location of the contact discontinuity can be significantly different between phases when the pulsar wind directly hits the Be disk and those when the pulsar wind collides with the Be wind. It is thus important to adequately take into account the circumstellar environment of the Be star, in order to construct a satisfactory model for this prototypical TeV binary.

We analyze the intrinsic polarization of two classical Be stars in the process of losing their circumstellar disks via a Be to normal B star transition originally reported by Wisniewski et al. During each of five polarimetric outbursts which interrupt these disk-loss events, we find that the ratio of the polarization across the Balmer jump (BJ+/BJ-) versus the V-band polarization traces a distinct loop structure as a function of time. Since the polarization change across the Balmer jump is a tracer of the innermost disk density whereas the V-band polarization is a tracer of the total scattering mass of the disk, we suggest such correlated loop structures in Balmer jump-V band polarization diagrams (BJV diagrams) provide a unique diagnostic of the radial distribution of mass within Be disks. We use the 3-D Monte Carlo radiation transfer code HDUST to reproduce the observed clockwise loops simply by turning "on/off" the mass decretion from the disk. We speculate that counter-clockwise loop structures we observe in BJV diagrams might be caused by the mass decretion rate changing between subsequent "on/off" sequences. Applying this new diagnostic to a larger sample of Be disk systems will provide insight into the time-dependent nature of each system's stellar decretion rate.

The possibility has been suggested that high-frequency quasi-periodic oscillations observed in low-mass X-ray binaries are resonantly excited disk oscillations in deformed (warped or eccentric) relativistic disks (Kato 2004). In this paper we examine this wave excitation process from a viewpoint somewhat different from that of previous studies. We study how amplitudes of a set of normal mode oscillations change secularly with time by their mutual couplings through disk deformation. As a first step, we consider the case where the number of oscillation modes contributing to the resonance coupling is two. The results show that two prograde oscillations interacting through disk deformation can grow if their wave energies have opposite signs.

We present a novel theoretical tool to analyze the dynamical behaviour of a Be disk fed by non-constant decretion rates. It is mainly based on the computer code HDUST, a fully three-dimensional radiative transfer code that has been successfully applied to study several Be systems so far, and the SINGLEBE code that solves the 1D viscous diffusion problem. We have computed models of the temporal evolution of different types of Be star disks for different dynamical scenarios. By showing the behaviour of a large number of observables (interferometry, polarization, photometry and spectral line profiles), we show how it is possible to infer from observations some key dynamical parameters of the disk.

Perhaps the most extreme examples of "Active OB stars" are the subset of high-mass X-ray binaries -- consisting of an OB star plus compact companion -- that have recently been observed by Fermi and ground-based Cerenkov telescopes like HESS to be sources of very high energy (VHE; up to 30 TeV) gamma-rays. This paper focuses on the prominent gamma-ray source, LS5039, which consists of a massive O6.5V star in a 3.9-day-period, mildly elliptical (e = 0.24) orbit with its companion, assumed here to be a black-hole or unmagnetized neutron star. Using 3-D SPH simulations of the Bondi-Hoyle accretion of the O-star wind onto the companion, we find that the orbital phase variation of the accretion follows very closely the simple Bondi-Hoyle-Lyttleton (BHL) rate for the local radius and wind speed. Moreover, a simple model, wherein intrinsic emission of gamma-rays is assumed to track this accretion rate, reproduces quite well Fermi observations of the phase variation of gamma-rays in the energy range 0.1-10 GeV. However for the VHE (0.1-30 TeV) radiation observed by the HESS Cerenkov telescope, it is important to account also for photon-photon interactions between the gamma-rays and the stellar optical/UV radiation, which effectively attenuates much of the strong emission near periastron. When this is included, we find that this simple BHL accretion model also quite naturally fits the HESS light curve, thus making it a strong alternative to the pulsar-wind-shock models commonly invoked to explain such VHE gamma-ray emission in massive-star binaries.

In order to study the origin of high-frequency quasi-periodic oscillations observed in X-ray binaries, Kato (2004) suggested a resonant excitation mechanism of disk oscillations in deformed disks. In this paper, we study numerically, following his formulation, whether trapped g-mode oscillations in a warped disk, where the warp amplitude varies with radius, can be excited by this mechanism. For simplicity, we adopt Newtonian hydrodynamic equations with relativistic expressions for the characteristic frequencies of disks. We also assume that the accretion disk is isothermal. We find that the fundamental modes of trapped g-mode oscillations with eigenfrequencies close to the maximum of epycyclic frequency are excited. The intermediate oscillations found are isolated in a narrow region around the resonance radius. After varying some parameters, we find that the growth rate increases as the warp amplitude or the black hole spin parameter increases, while it decreases as the sound speed increases.

We report near-IR spectroscopic observations of the Eta Carinae massive binary system during 2008-2009 using VLT/CRIRES. We detect a strong, broad absorption wing in He I 10833 extending up to -1900 km/s across the 2009.0 spectroscopic event. Archival HST/STIS ultraviolet and optical data shows a similar high-velocity absorption (up to -2100 km/s) in the UV resonance lines of Si IV 1394, 1403 across the 2003.5 event. UV lines from low-ionization species, such as Si II 1527, 1533 and C II 1334, 1335, show absorption up to -1200 km/s, indicating that the absorption with v from -1200 to -2100 km/s originates in a region markedly faster and more ionized than the nominal wind of the primary star. Observations obtained at the OPD/LNA during the last 4 spectroscopic cycles (1989-2009) also display high-velocity absorption in He I 10833 during periastron. Based on the OPD/LNA dataset, we determine that material with v < -900 km/s is present in the phase range 0.976 < phi < 1.023 of the spectroscopic cycle, but absent in spectra taken at phi < 0.947 and phi > 1.049. Therefore, we constrain the duration of the high-velocity absorption to be 95 to 206 days (or 0.047 to 0.102 in phase). We suggest that the high-velocity absorption originates from shocked gas in the wind-wind collision zone, at distances of 15 to 45 AU in the line-of-sight to the primary star. Using 3-D hydrodynamical simulations of the wind-wind collision zone, we find that the dense high-velocity gas is in the line-of-sight to the primary star only if the binary system is oriented in the sky so that the companion is behind the primary star during periastron, corresponding to a longitude of periastron of omega ~ 240 to 270 degrees. We study a possible tilt of the orbital plane relative to the Homunculus equatorial plane and conclude that our data are broadly consistent with orbital inclinations in the range i=40 to 60 degrees.

We report on optical high-dispersion spectroscopic monitoring observations of the Be/X-ray binary A0535+26/V725 Tau, carried out from November 2005 to March 2009. The main aim of these monitoring observations is to study spectral variabilities in the Be disc, on both the short (a week or so) and long (more than hundreds of days) timescales, by taking long-term frequent observations. Our four-year spectroscopic observations indicate that the V/R ratio, i.e., the relative intensity of the violet (V) peak to the red (R) one, of the double-peaked H-alpha line profile varies with a period of 500 days. The H-beta line profile also varies in phase with the H-alpha profile. With these observations covering two full cycles of the V/R variability, we reconstruct the 2-D structure of the Be disc by applying the Doppler tomography method to the H-alpha and H-beta emission line profiles, using a rigidly rotating frame with the V/R variability period. The resulting disc structure reveals non-axisymmetric features, which can be explained by a one-armed perturbation in the Be disc. It is the first time that an eccentric disc structure is directly detected by using a method other than the interferometric one.

We report on H-alpha spectroscopy of the 2009.0 spectroscopic event of eta Carinae collected via SMARTS observations using the CTIO 1.5 m telescope and echelle spectrograph. Our observations were made almost every night over a two month interval around the predicted minimum of eta Car. We observed a significant fading of the line emission that reached a minimum seven days after the X-ray minimum. About 17 d prior to the H-alpha flux minimum, the H-alpha profile exhibited the emergence of a broad, P Cygni type, absorption component (near a Doppler shift of -500 km/s) and a narrow absorption component (near -144 km/s and probably associated with intervening gas from the Little Homunculus Nebula). All these features were observed during the last event in 2003.5 and are probably related to the close periastron passage of the companion. We argue that these variations are consistent with qualitative expectations about changes in the primary star's stellar wind that result from the wind-wind collision with a massive binary companion and from atmospheric eclipses of the companion.

Emission lines formed in decretion disks of Be stars often undergo long-term cyclic variations, especially in the violet-to-red (V/R) ratio of their primary components. From observations of the bright Be-shell star zeta Tau, the possibly broadest and longest data set illustrating the prototype of this behaviour was compiled from our own and archival observations. It comprises optical and infrared spectra, broad-band polarimetry, and interferometric observations. From 3 V/R cycles between 1997 and 2008, a mean cycle length in H alpha of 1400-1430 days was derived. After each minimum in V/R, the shell absorption weakens and splits into two components, leading to 3 emission peaks. This phase makes the strongest contribution to the variability in cycle length. V/R curves of different lines are shifted in phase. Lines formed on average closer to the central star are ahead of the others. The shell absorption lines fall into 2 categories differing in line width, ionization/excitation potential, and variability of the equivalent width. The interferometry has resolved the continuum and the line emission in Br gamma and He I 2.06. The phasing of the Br gamma emission shows that the photocenter of the line-emitting region lies within the plane of the disk but is offset from the continuum source. The plane of the disk is constant throughout the observed V/R cycles. The self-consistent, one-armed, disk-oscillation model is developed in Paper II.

Aims. In this paper we model, in a self-consistent way, polarimetric, photometric, spectrophotometric and interferometric observations of the classical Be star $\zeta$ Tauri. Our primary goal is to conduct a critical quantitative test of the global oscillation scenario. Methods. We have carried out detailed three-dimensional, NLTE radiative transfer calculations using the radiative transfer code HDUST. For the input for the code we have used the most up-to-date research on Be stars to include a physically realistic description for the central star and the circumstellar disc. We adopt a rotationally deformed, gravity darkened central star, surrounded by a disc whose unperturbed state is given by a steady-state viscous decretion disc model. We further assume that disc is in vertical hydrostatic equilibrium. Results. By adopting a viscous decretion disc model for $\zeta$ Tauri and a rigorous solution of the radiative transfer, we have obtained a very good fit of the time-average properties of the disc. This provides strong theoretical evidence that the viscous decretion disc model is the mechanism responsible for disc formation. With the global oscillation model we have successfully fitted spatially resolved VLTI/AMBER observations and the temporal V/R variations of the H$\alpha$ and Br$\gamma$ lines. This result convincingly demonstrates that the oscillation pattern in the disc is a one-armed spiral. Possible model shortcomings, as well as suggestions for future improvements, are also discussed.

Circumbinary disks are considered to exist in a wide variety of astrophysical objects, e.g., young binary stars, protoplanetary systems, and massive binary black hole systems in active galactic nuclei (AGNs). However, there is no definite evidence for the circumbinary disk except for some in a few young binary star systems. In this Letter, we study possible oscillation modes in circumbinary disks around eccentric and circular binaries. We find that progarde, nonaxisymmetric waves are induced in the inner part of the circumbinary disk by the tidal potential of the binary. Such waves would cause variabilities in emission line profiles from circumbinary disks. Because of prograde precession of the waves, the distance between each component of the binary and the inner edge of the circumbinary disk varies with the beat period between the precession period of the wave and the binary orbital period. As a result, light curves from the circumbinary disks are also expected to vary with the same period. The current study thus provides a new method to detect circumbinary disks in various astrophysical systems.

LS 5039 is a TeV gamma-ray binary with extended radio emission. It consists of a compact object in the mildly eccentric (e=0.35), 3.9-day orbit around a massive O star. The nature of the compact object is not yet established. In this paper, assuming that the compact object is a black hole, we study the accretion of O-star wind by the black hole, by performing three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations. In order to roughly emulate the effect of the stellar radiation effectively canceling the stellar gravity, we assume that the O star's gravity does not exert on the wind. The wind particles are ejected with half the observed terminal velocity in a narrow range of azimuthal and vertical angles toward the black hole, in order to emulate the wind significantly slower than the terminal speed, and optimize the resolution and computational efficiency of simulations. We find that the mass-accretion rate closely follows the classical Bondi-Hoyle-Littleton accretion rate, which is of the order of 10^{16}g/s around periastron. The accretion rate at this level would provide jets enough power to produce the gamma-rays detected by HESS. Since the accretion peak occurs near the periastron passage, we need a strong gamma-ray absorption around periastron in order for the microquasar scenario to be consistent with the observed orbital modulation of the TeV gamma-ray flux.

We study the effects of the tidal interaction with the companion, via orbital separation and binary mass ratio, on the global one-armed oscillation modes in disks around binary Be stars. Our model takes into account the three-dimensional effect that contributes to the mode confinement, which was recently found by Ogilvie(2008). We find that the one-armed oscillations are well confined in systems with disks larger than a few tens of stellar radii. In such systems, the oscillation period depends little on the binary parameters. On the other hand, in systems with smaller disks, where the mode confinement is incomplete, the oscillation period increases with increasing orbital separation and/or decreasing binary mass ratio. The eigenmode is insensitive to the spectral type of the central star. Our results suggest that the dependence of V/R oscillation period on the orbital separation and binary mass ratio should be observed only in short period binary systems, and that, for systems with a similar orbital period, those with higher mass ratios will show shorter V/R variations.

The very massive star system $\eta$ Carinae exhibits regular 5.54-year (2024-day) period disruptive events in wavebands ranging from the radio to X-ray. There is a growing consensus that these events likely stem from periastron passage of an (as yet) unseen companion in a highly eccentric ($\epsilon \sim 0.9$) orbit. This paper presents three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulations of the orbital variation of the binary wind-wind collision, and applies these to modeling the X-ray light curve observed by the Rossi X-ray Timing Explorer (RXTE). By providing a global 3-D model of the phase variation of the density of the interacting winds, the simulations allow computation of the associated variation in X-ray absorption, presumed here to originate from near the apex of the wind-wind interaction cone. We find that the observed RXTE light curve can be readily fit if the observer's line of sight is within this cone along the general direction of apastron. Specifically, the data are well fit by an assumed inclination $i = 45^{\circ}$ for the orbit's polar axis, which is thus consistent with orbital angular momentum being along the inferred polar axis of the Homunculus nebula. The fits also constrain the position angle $\phi$ that an orbital-plane projection makes with the apastron side of the semi-major axis, strongly excluding positions $\phi < 9^{\circ}$ along or to the retrograde side of the axis, with the best fit position given by $\phi = 27^{\circ}$. Overall the results demonstrate the utility of a fully 3-D dynamical model for constraining the geometric and physical properties of this complex colliding-wind binary system.

We study colliding winds in the superluminous binary eta Carinae by performing three-dimensional, Smoothed Particle Hydrodynamics (SPH) simulations. For simplicity, we assume both winds to be isothermal. We also assume that wind particles coast without any net external forces. We find that the lower density, faster wind from the secondary carves out a spiral cavity in the higher density, slower wind from the primary. Because of the phase-dependent orbital motion, the cavity is very thin on the periastron side, whereas it occupies a large volume on the apastron side. The model X-ray light curve using the simulated density structure fits very well with the observed light curve for a viewing angle of i=54 degrees and phi=36 degrees, where i is the inclination angle and phi is the azimuth from apastron.

We investigate the tidal interactions of a red giant with a main sequence in the dense stellar core of globular clusters by Smoothed Particle Hydrodynamics method. Two models of $0.8 \msun$ red giant with the surface radii 20 and $85 R_\sun$ are used with 0.6 or $0.8M_\sun$ main sequence star treated as a point mass. We demonstrate that even for the wide encounters that two stars fly apart, the angular momentum of orbital motion can be deposited into the red giant envelope to such an extent as to trigger rotational mixing and to explain the fast rotation observed for the horizontal branch stars, and also that sufficient mass can be accreted on the main sequence stars to disguise their surface convective zone with the matter from the red giant envelope. On the basis of the present results, we discuss the parameter dependence of these transfer characteristics with non-linear effects taken into account, and derive fitting formulae to give the amounts of energy and angular momentum deposited into the red giant and of mass accreted onto the perturber as functions of stellar parameters and the impact parameter of encounter. These formulae are applicable to the encounters not only of the red giants but also of the main sequence stars, and useful in the study of the evolution of stellar systems with the star-star interactions taken into account.

LS I +61 303 is a puzzling object detected from radio up to high-energy gamma-rays. Variability has recently been observed in its high-energy emission. The object is a binary system, with a compact object and a Be star as primary. The nature of the secondary and the origin of the gamma-ray emission are not clearly established at present. Recent VLBA radio data have been used to claim that the system is a Be/neutron star colliding wind binary, instead of a microquasar. We review the main views on the nature of LS I +61 303 and present results of 3D SPH simulations that can shed some light on the nature of the system. Our results support an accretion powered source, compatible with a microquasar interpretation.

LS I +61 303 is a puzzling Be/X-ray binary with variable gamma-ray emission at up TeV energies. The nature of the compact object and the origin of the high-energy emission are unclear. One family of models invokes particle acceleration in shocks from the collision between the B-star wind and a relativistic pulsar wind, while another centers on a relativistic jet powered by accretion. Recent high-resolution radio observations showing a putative "cometary tail" pointing away from the Be star near periastron have been cited as support for the pulsar-wind model. We wish here to carry out a quantitative assessment of these competing models for this extraordinary source. We apply a 3D SPH code for dynamical simulations of both the pulsar-wind-interaction and accretion-jet models. The former yields a description of the shape of the wind-wind interaction surface. The latter provides an estimation of the accretion rate. The results allow critical evaluation of how the two distinct models confront the data in various wavebands under a range of conditions. When one accounts for the 3D dynamical wind interaction under realistic constraints for the relative strength of the B-star and pulsar winds, the resulting form of the interaction front does not match the putative "cometary tail" claimed from radio observations. On the other hand, dynamical simulations of the accretion-jet model indicate that the orbital phase variation of accretion power includes a secondary broad peak well away from periastron, thus providing a plausible way to explain the observed TeV gamma ray emission toward apastron. We conclude that the colliding-wind model is not clearly established for LS I +61 303, while the accretion-jet model can reproduce many key characteristics of the observed TeV gamma-ray emission.

We present optical linear spectropolarimetry of the microquasar LS I +61$^{\circ}$ 303. The continuum emission is mildly polarized (up to 1.3 %) and shows almost no temporal change. We find a distinct change of polarization across the H$\alpha$ emission line, indicating the existence of polarization component intrinsic to the microquasar. We estimate the interstellar polarization (ISP) component from polarization of the H$\alpha$ line and derive the intrinsic polarization component. The wavelength dependence of the intrinsic component is well explained by Thomson scattering in equatorial disk of the Be-type mass donor. The position angle (PA) of the intrinsic polarization $\sim 25^{\circ}$ represents the rotational axis of the Be disk. This PA is nearly perpendicular to the PA of the radio jet found during quiescent phases. Assuming an orthogonal disk-jet geometry around the compact star, the rotational axis of the accretion disk is almost perpendicular to that of the Be disk. Moreover, according to the orbital parameters of the microquasar, the compact star is likely to get across the Be disk around their periastron passage. We discuss the peculiar circumstellar structure of this microquasar inferred from our observation and possible connection with its high-energy activities.

We numerically study the long-term evolution of the accretion disc around the neutron star in a coplanar Be/X-ray binary with a short period and a moderate eccentricity. From three dimensional Smoothed Particle Hydrodynamics simulations, we find that the disc evolves through three distinct phases, each characterized by different mass accretion patterns. In the first "developing phase", the disc is formed and develops towards a nearly Keplerian disc. It has a relatively large, double-peaked mass-accretion rate with the higher peak by the direct accretion at periastron, which is followed by the lower peak by the accretion induced by a one-armed spiral wave. In the second "transition phase", the disc is approximately Keplerian and grows with time. The mass-accretion rate increases as the disc grows. In the second phase, there is a transition in the mass accretion rate from a double peaked to a single peaked pattern. In the final quasi-steady state, the mass-accretion rate is on average balanced with the mass-transfer rate from the Be disc and exhibits a regular orbital modulation. In the quasi-steady state, the mass-accretion rate has a single peak by the wave-induced accretion as in a later stage of the transition phase. The orbital modulation of X-ray maxima could provide not only a circumstantial evidence for the persistent disc but also an observational diagnosis of the disc evolutionary state.

We present low-resolution ($R\sim 90$) and medium-resolution ($R\sim 2500$) spectropolarimetry of Nova V475 Sct with the HBS instrument, mounted on the 0.91-m telescope at the Okayama Astrophysical Observatory, and with FOCAS, mounted on the 8.2-m Subaru telescope. We estimated the interstellar polarization toward the nova from the steady continuum polarization components and H$\alpha$ line emission components. After subtracting the interstellar polarization component from the observations, we found that the H$\alpha$ emission seen on 2003 October 7 was clearly polarized. In the polarized flux spectrum, the H$\alpha$ emission had a distinct red wing extending to $\sim +4900$ km s$^{-1}$ and a shoulder around $+3500$ km s$^{-1}$, showing a constant position angle of linear polarization $\theta_{\rm *}\simeq 155\arcdeg\pm 15\arcdeg$. This suggests that the nova had an asymmetric outflow with a velocity of $v_{\rm wind}\simeq 3500$ km s$^{-1}$ or more, which is six times higher than the expansion velocity of the ionized shell at the same epoch. Such a high-velocity component has not previously been reported for a nova in the `moderately fast' speed class. Our observations suggest the occurrence of violent mass-loss activity in the nova binary system even during the common-envelope phase. The position angle of the polarization in the H$\alpha$ wing is in good agreement with that of the continuum polarization found on 2003 September 26 ($p_{\rm *}\simeq 0.4$--0.6 %), which disappeared within the following 2 d. The uniformity of the PA between the continuum polarization and the wing polarization on October 7 suggests that the axis of the circumstellar asymmetry remained nearly constant during the period of our observations.

We compare the accretion flow onto the neutron star induced by Roche lobe overflow with that by the overflow from the Be disk, in a zero eccentricity, short period binary with the same mass transfer rate, performing three-dimensional Smoothed Particle Hydrodynamics simulations. We find that a persistent accretion disk is formed around the neutron star in both cases. The circularization radius of the material transferred via Roche lobe overflow is larger than that of the material transfered from the Be disk. Thus, the growth of the accretion disk in the former case becomes significantly slower than in the latter case. In both cases, the mass accretion rate is very small and varies little with orbital phase, which is consistent with the observed X-ray behaviour of Be/X-ray binaries with circular orbits (e.g. XTE J1543-568).

we study the long-term evolution of the accretion disk around the neutron star in Be/X-ray binaries. We confirm the earlier result by Hayasaki & Okazaki (2004) that the disk evolves via a two-stage process, which consists of the initial developing stage and the later developed stage. The peak mass-accretion rate is distributed around apastron after the disk is fully developed. This indicates that the modulation of the mass accretion rate is essentially caused by an inward propagation of the one-armed spiral wave. The X-ray luminosity peak around the apastron could provide circumstatial evidence for an persistent disk around the neutron star in Be/X-ray binaries.

A probable binary period has been detected in the optical counterpart to the X-ray source CXOU J005455.6-724510 = RX J0054.9-7245 = AXJ0054.8-7244 = SXP504 in the Small Magellanic Cloud. This source was detected by Chandra on 04 Jul 2002 and subsequently observed by XMM-Newton on 18 Dec 2003. The source is coincident with an Optical Gravitational Lensing (OGLE) object in the lightcurves of which several optical outburst peaks are visible at ~ 268 day intervals. Timing analysis shows a period of 268.6 +/- 0.1 days at > 99% significance. Archival Rossi X-ray Timing Explorer (RXTE) data for the 504s pulse-period has revealed detections which correspond closely with predicted or actual peaks in the optical data. The relationship between this orbital period and the pulse period of 504s is within the normal variance found in the Corbet diagram.

We study non-axisymmetric structure of accretion discs in Be/X-ray binaries, performing three dimensional Smoothed Particle Hydrodynamics simulations for a coplanar system with a short period and a moderate eccentricity. We find that ram pressure due to the phase-dependent mass transfer from the Be-star disc excites a one-armed, trailing spiral structure in the accretion disc around the neutron star. The spiral wave is transient; it is excited around the periastron passage, when the material is transferred from the Be disc, and is gradually damped afterwards. The disc changes its morphology from circular to eccentric with the development of the spiral wave, and then from eccentric to circular with the decay of the wave during one orbital period. It turns out that the inward propagation of the spiral wave significantly enhances the mass-accretion rate onto the neutron star. Thus, the detection of an X-ray luminosity peak corresponding to the peak in enhanced mass-accretion rate provides circumstantial evidence that an accretion disc is present in Be/X-ray binaries.

The non-axisymmetric structure of accretion disks around the neutron star in Be/X-ray binaries is studied by analyzing the results from three dimensional (3D) Smoothed Particle Hydrodynamics (SPH) simulations. It is found that ram pressure due to the phase-dependent mass transfer from the Be-star disk excites a one-armed, trailing spiral structure in the accretion disk around the neutron star. The spiral wave has a transient nature; it is excited around the periastron, when the material is transferred from the Be disk, and is gradually damped afterwards. It is also found that the orbital phase-dependence of the mass-accretion rate is mainly caused by the inward propagation of the spiral wave excited in the accretion disk.

We study the accretion on to the neutron star in Be/X-ray binaries, using a 3D SPH code and the data imported from a simulation by \citet{oka2} for a coplanar system with a short period ($P_{\rm orb}=24.3 \rm{d}$) and a moderate eccentricity $(e=0.34)$, which targeted the Be/X-ray binary 4U 0115+63. For simplicity, we adopt the polytropic equation of state. We find that a time-dependent accretion disc is formed around the neutron star regardless of the simulation parameters. In the long term, the disc evolves via a two-stage process, which consists of the initial developing stage and the later developed stage. The developed disc is nearly Keplerian. In the short term, the disc structure modulates with the orbital phase. The disc shrinks at the periastron passage of the Be star and restores its radius afterwards. The accretion rate on to the neutron star is also phase dependent, but its peak is broader and much lower than that of the mass-transfer rate from the Be disc, unless the polytropic exponent is as large as 5/3. Our simulations show that the truncated Be disk model for Be/X-ray binaries is consistent with the observed X-ray behaviour of 4U 0115+63.

We study accretion onto the neutron star in Be/X-ray binaries, using a 3D SPH code and the data imported from a high resolution simulation by Okazaki et al.\ (2002) for a coplanar system with a short period ($P_{\rm orb}=24.3 \rm d$) and moderate eccentricity $(\rm{e=0.34})$. We find that a time-dependent accretion disk is formed around the neutron star in Be/X-ray binaries. The disk shrinks after the periastron passage of the Be star and restores its radius afterwards. Our simulations show that the truncated Be disk model for Be/X-ray binaries is consistent with the observed X-ray behavior.

We study the viscous effects on the interaction between the coplanar Be-star disc and the neutron star in Be/X-ray binaries, using a three-dimensional, smoothed particle hydrodynamics code. For simplicity, we assume the Be disc to be isothermal at the temperature of half the stellar effective temperature. In order to mimic the gas ejection process from the Be star, we inject particles with the Keplerian rotation velocity at a radius just outside the star. Both Be star and neutron star are treated as point masses. We find that the Be-star disc is effectively truncated if the Shakura-Sunyaev viscosity parameter alpha_SS >> 1, which confirms the previous semi-analytical result. In the truncated disc, the material decreted from the Be star accumulates, so that the disc becomes denser more rapidly than if around an isolated Be star. The resonant truncation of the Be disc results in a significant reduction of the amount of gas captured by the neutron star and a strong dependence of the mass capture rate on the orbital phase. We also find that an eccentric mode is excited in the Be disc through direct driving due to a one-armed bar potential of the binary. The strength of the mode becomes greater in the case of a smaller viscosity. In a high-resolution simulation with alpha_SS=0.1, the eccentric mode is found to precess in a prograde sense. The mass capture rate by the neutron star modulates as the mode precesses.

When applied to Be/X-ray binaries, the viscous decretion disc model for Be stars naturally predicts the truncation of the circumstellar disc. The distance at which the circumstellar disc is truncated depends mainly on the orbital parameters and the viscosity. In systems with low eccentricity, the disc is expected to be truncated at the 3:1 resonance radius, for which the gap between the disc outer radius and the critical lobe radius of the Be star is so wide that, under normal conditions, the neutron star cannot accrete enough gas at periastron passage to show Type I X-ray outbursts. These systems will display only Type II X-ray outbursts. On the other hand, in systems with high orbital eccentricity, the disc truncation occurs at a much higher resonance radius, which is very close to or slightly beyond the critical lobe radius at periastron unless the viscosity is very low. In these systems, disc truncation cannot be efficient, allowing the neutron star to capture gas from the disc at every periastron passage and display Type I outbursts regularly. In contrast to the rather robust results for systems with low eccentricity and high eccentricity, the result for systems with moderate eccentricity depends on rather subtle details. Systems in which the disc is truncated in the vicinity of the critical lobe will regularly display Type I outbursts, whereas those with the disc significantly smaller than the critical lobe will show only Type II outbursts under normal conditions and temporary Type I outbursts when the disc is strongly disturbed.

We present multi-wavelength long-term monitoring observations of V635 Cas, the optical counterpart to the transient X-ray pulsar 4U0115+63. The evolution of emission lines and photometric magnitudes indicates that the Be star undergoes relatively fast (~ 3-5 yr) quasi-cyclic activity, losing and reforming its circumstellar disc. We show that the general optical, infrared and X-ray behaviour can be explained by the dynamical evolution of the viscous circumstellar disc around the Be star. After each disc-loss episode, the disc starts reforming and grows until it reaches the radius at which the resonant interaction of the neutron star truncates it. At some point, the disc becomes unstable to (presumably radiative) warping and then tilts and starts precessing. The tilting is very large and disc precession leads to a succession of single-peaked and shell profiles in the emission lines. Type II X-ray outbursts take place after the disc has been strongly disturbed and we speculate that the distortion of the disc leads to interaction with the orbiting neutron star. We discuss the implications of these correlated optical/X-ray variations for the different models proposed to explain the occurrence of X-ray outbursts in Be/X-ray binaries. We show that the hypothesis of mass ejection events as the cause of the spectacular variability and X-ray outbursts is unlikely to be meaningful for any Be/X-ray binary.

We present photometry and high SNR spectroscopy in the classification region of V635 Cas, the optical counterpart to the transient X-ray pulsator 4U 0115+63, taken at a time when the circumstellar envelope had disappeared. V635 Cas is classified as a B0.2Ve star at a distance of 7-8 kpc. We use the physical parameters derived from these observations and the orbit derived from X-ray observations to elaborate a model of the system based on the theory of decretion discs around Be stars. We show that the disc surrounding the Be star must be truncated by the tidal/resonant interaction with the neutron star and cannot be in a steady state. This explains many of the observed properties of 4U 0115+63. In particular, because of this effect, under normal circumstances, the neutron star cannot accrete from the disc, which explains the lack of regular Type I outbursts from the source.

We present a new scenario for the behaviour of Be/X-ray binaries based on long-term multiwavelength monitoring and the decretion disc model. The circumstellar discs of the primaries are truncated because of the tidal and resonant effect of the neutron star. The geometry of the systems and the value of viscosity determine the presence or absence of Type I X-ray outbursts. The interaction of a strongly disturbed disc with the neutron star originates Type II X-ray and optical outbursts.

We study the characteristics of the outflow in disks of Be stars, based on the viscous decretion disk scenario. In this scenario, the matter ejected from the equatorial surface of the star drifts outward because of the effect of viscosity, and forms the disk. For simplicity, we adopt the alpha-prescription for the viscous stress, and assume the disk to be isothermal. Solving the resulting wind equations, we find that a transonic solution exists for any value of alpha. The sonic point is located at r>100R for plausible values of parameters, where R is the stellar radius. The sonic radius is smaller for higher temperature and/or larger radiative force. We also find that the topology of the sonic point is nodal for alpha>0.95, while it is of saddle type for alpha<0.9. We expect that the sonic point in the former case is unstable, whereas that in the latter case is stable. The outflow is highly subsonic in the inner part of the disk. Roughly, the outflow velocity increases linearly with r and the surface density decreases as r^{-2}. Interestingly, the disk is near Keplerian in the inner subsonic region, while it is angular momentum conserving in the outer subsonic region and in the supersonic region. Our results, together with the observed range of the base density for Be star disks, suggest that the mass loss rate in the equatorial region is at most comparable with that in the polar region.