Search

We study the effect of irradiation from two accretion disks (minidisks) around respective black holes of stellar to intermediate masses in a circular binary on the spectrum of a circumbinary disk (CBD) surrounding them. We assume the CBD to be a standard disk and adopt the orbit-averaged irradiation flux because the viscous timescale is much longer than the orbital period. We then solve the energy equation both analytically and numerically to compute the CBD temperature distribution and the corresponding disk spectrum. We find that the analytically calculated spectra are in good agreement with the numerical ones. The CBD spectrum is almost independent of the binary mass ratio. We also find that the combined spectra of two minidisks and the CBD have double peaks, one peak in the soft X-ray band and the other in the infrared (IR) band. The former peak comes from the two minidisks, while the latter peak from the CBD. The observed flux density increases with frequency as $\nu^{1/3}$ towards the soft X-ray peak, while it decreases with frequency away from the IR peak as $\nu^{-5/3}$. The latter feature is testable with near-IR observations with Subaru and JWST., Comment: 19 pages, 7 figures. Submitted to ApJ

We have performed the near-infrared photometric monitoring observations of two TeV gamma-ray binaries with O-stars (LS 5039 and 1FGL J1018.6-5856), using IRSF/SIRIUS at SAAO, in order to study the stellar parameters and their perturbations caused by the binary interactions. The whole orbital phase was observed multiple times and no significant variabilities including orbital modulations are detected for both targets. Assuming that the two systems are colliding wind binaries, we estimate the amplitude of flux variation caused by the difference in the optical depth of O-star wind at inferior conjunction, where the star is seen through the cavity created by pulsar wind, and other orbital phases without pulsar-wind intervention. The derived amplitude is <0.001 mag, which is about two orders of magnitude smaller than the observed upper limit. Also using the upper limits of the near-infrared variability, we for the first time obtain the upper limit of the dust formation rate resulting from wind-wind collision in O-star gamma-ray binaries., Comment: 7 figures, 2 tables, Accepted to publications in PASJ

Morphology of nova ejecta is essential for fully understanding the physical processes involved in nova eruptions. We studied the 3D morphology of the expanding ejecta of the extremely slow nova V1280 Sco with a unique light curve. Synthetic line profile spectra were compared to the observed [O III] 4959, 5007 and [N II] 5755 emission line profiles in order to find the best-fit morphology, inclination angle, and maximum expansion velocity of the ejected shell. We derive the best fitting expansion velocity, inclination, and squeeze as $V_{\rm exp} = 2100^{+100}_{-100}$ \,km\,s$^{-1}$, $i = 80^{+1}_{-3}$ deg, and $squ = 1.0^{+0.0}_{-0.1}$ using [O III] line profiles, and $V_{\rm exp} = 1600^{+100}_{-100}$ \,km\,s$^{-1}$, $i = 81^{+2}_{-4}$ deg, and $squ = 1.0^{+0.0}_{-0.1}$ using [N II] 5755 line profile. A high inclination angle is consistent with the observational results showing multiple absorption lines originating from clumpy gases which are produced in dense and slow equatorially focused outflows. Based on additional observational features such as optical flares near the maximum light and dust formation on V1280 Sco, a model of internal shock interaction between slow ejecta and fast wind proposed for the $\gamma$-ray emission detected in other novae seems to be applicable to this extremely slow and peculiar nova. Increasing the sample size of novae whose morphology is studied will be helpful in addressing long-standing mysteries in novae such as the dominant energy source to power the optical light at the maximum, optical flares near the maximum, clumpiness of the ejecta, and dust formation., Comment: 21 pages, 8 figures, accepted for publication in ApJ

We use a smoothed particle hydrodynamics (SPH) code to examine the effects of a binary companion on a Be star disk for a range of disk viscosities and misalignment angles, i.e. the angle between the orbital plane and the primary's spin axis. The density structures in the disk due to the tidal interaction with the binary companion are investigated. Expanding on our previous work, the shape and density structure of density enhancements due to the binary companion are analyzed and the changes in observed interferometric features due to these orbiting enhancements are also predicted. We find that larger misalignment angles and viscosity values result in more tightly wound spiral arms with densities that fall-off more slowly with radial distance from the central star. We show that the orbital phase has very little effect on the structure of the spiral density enhancements. We demonstrate that these spiral features can be detected with an interferometer in H$\alpha$ and K-band emission. We also show that the spiral features affect the axis ratios determined by interferometry depending on the orientation of these features and the observer. For example, our simulations show that the axis ratios can vary by 20% for our co-planar binary disk system depending on the location of the disk density enhancements., Comment: 13 pages, 14 figures

The binary of the pulsar PSRB1259$-$63 and the Be star LS 2883 has been observed at the 2010 and 2014 periastron passages in the near-infrared (NIR) bands using the IRSF/SIRIUS and SIRPOL. The light curves in the J-,H-, and Ks-bands are almost identical in these periastron passages. A flare starts no later than 10 days before periastron and the maximum brightening of about 0.1 magnitude is observed 12--17 days after periastron. The rising part of the light curve is steeper and reaches a peak slightly earlier in the Ks-band than in the other bands, thus a characteristic track appears on the NIR color-magnitude diagram. The time lag between the NIR light curves indicates that the variation in the Be circumstellar disk first occurs in an outer region. We propose that the initial rapid contraction followed by the gradual expansion of the disk is evoked by the rapidly changing tidal torque around periastron and the resultant change of the optically thick area causes the observed NIR light curves., Comment: 13 pages, 4 figures. Accepted for publication in PASJ

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., Comment: Accepted for publication in PASJ

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., Comment: 9 pages, 10 figures. Accepted by MNRAS

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., Comment: 20 pages, 19 figures; accepted for publication in MNRAS

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., Comment: 10 pages, 10 figures. Accepted by MNRAS

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., Comment: 13 pages 14 figures accepted for publication MNRAS

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., Comment: 8 pages, 8 figures

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., Comment: 12 pages, 8 figures, 3 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

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., Comment: 11 pages, 10 figures, to be published in MNRAS

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., Comment: 12 pages, 16 figures; accepted for publication by MNRAS

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., Comment: 6 pages, 3 figures; to appear in the proceedings of the conference "The B[e] phenomenon: Forty years of studies" held in Prague, 2016

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., Comment: 15 pages, 30 figures

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., Comment: 11 pages, 11 figures, accepted for publication in MNRAS

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., Comment: 4 pages, 4 figures

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., Comment: 20 pages, 4 figures, accepted to ApJ Letter

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., Comment: 14 pages, 9 figures, 3 tables

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., Comment: 17 pages, 6 figures, submitted to JCAP. arXiv admin note: substantial text overlap with arXiv:1406.2317

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., Comment: 16 pages, 3 figures, 2 tables, ApJ Letters, in press

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., Comment: 28 pages, 5 figures, 1 table, accepted for publication in ApJ. arXiv admin note: substantial text overlap with arXiv:1406.2317

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., Comment: 23 pages, 4 figures, 2 tables, accepted for publication in ApJ

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., Comment: 13 pages, 7 figures

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., Comment: 2 pages, 2 figures, Proceedings from "Suzaku-Maxi 2014: Expanding the Frontiers of the X-ray Universe"

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., Comment: 2 pages, 3 figures, Proceedings from "Suzaku-MAXI 2014: Expanding the Frontiers of the X-ray Universe"

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., Comment: 38 pages, 25 figures, 1 table. Accepted for publication in MNRAS

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., Comment: accepted to PASJ

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., Comment: 35 pages, 6 figures, 3 tables, accepted to ApJ

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., Comment: 16 pages, 10 figures, accepted for publication in PASJ

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., Comment: Accepted in ApJ

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., Comment: 32 pages, 6 figues. Accepted for publication in ApJ

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., Comment: 4 pages, 3 figures, accepted for publication in ApJ Letters

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., Comment: 23 pages, 12 color figures, plus 2 online-only appendices (available in the /anc folder of the Source directory). Accepted for publication in MNRAS

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., Comment: 3 pages, 3 figures, to appear in "Four Decades of Research on Massive Stars: A Scientific Meeting in Honour of Anthony F. J. Moffatt", eds. L. Drissen, C. Robert, and N. St-Louis, ASP Conference Series

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., Comment: accepted to PASJ

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., Comment: 10 pages, 5 figures, accepted for publication in PASJ

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., Comment: 14 pages, 1 figure, to be published in PASJ 63, No.2 (2011)

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., Comment: IAU 272 symposium proceedings, 3 pages, 1 figure

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., Comment: 20 pages, 4 figures, accepted for publication in PASJ

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., Comment: 18 pages, 15 figures, accepted for publication in A&A; high-resolution PDF version available also at http://www.mpifr.de/staff/jgroh/etacar.html

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., Comment: (10 pages, 9 figures, accepted to MNRAS)

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., Comment: 21 pages, 6 figures, 1 table, accepted to AJ

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., Comment: 1 pdf file, 17 pages. accepted to Astronomy and Astrophysics

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., Comment: 14 pages, 9 figures, accepted to A&A

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., Comment: 7 pages, 3 figures, submitted to the proceedings of the "7th INTEGRAL Workshop - An INTEGRAL View of Compact Objects", Copenhagen, 8-11 September 2008

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., Comment: 9 pages, 5 figures, accepted for publication in PASJ

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., Comment: Accepted for publication in ApJ Letters (3/Dec/2008), 4pages, 1figure