Source author record

Atsuo T. Okazaki

Atsuo T. Okazaki appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2022arXiv

Morpho-kinematic modeling of the expanding ejecta of the extremely slow nova V1280 Scorpii

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 $γ$-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.

preprint2021arXiv

ALMA Observations of PSR B1259-63/LS 2883 in an Inactive Period: Variable Circumstellar Disk?

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.

preprint2016arXiv

Be discs in binary systems I. Coplanar orbits

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.

preprint2016arXiv

Effect of Be Disk Evolution on Global One-Armed Oscillations

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.

preprint2016arXiv

Radiative transfer on decretion discs of Be binaries

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.

preprint2015arXiv

Probing the nature of the TeV gamma-ray binary HESS J0632+057 by monitoring Be disk variability

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.

preprint2015arXiv

Warping and tearing of misaligned circumbinary disks around eccentric SMBH binaries

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.

preprint2014arXiv

3D Dynamical Modeling of Wind Accretion in Cyg X-3

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.

preprint2014arXiv

3D Hydrodynamic & Radiative Transfer Models of X-ray Emission from Colliding Wind Binaries

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 $η$ Carinae and WR 140. The models reproduce the observations well away from periastron passage, but only $η$ Carinae's X-ray spectrum is reproduced at periastron; the WR 140 model produces too much flux during this more complicated phase.

preprint2014arXiv

Radiation-Driven Warping of Circumbinary Disks Around Eccentric Young Star Binaries

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.

preprint2014arXiv

Sudden End of X-Ray Outbursts around Periastron of Circinus X-1 Observed with MAXI

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.

preprint2014arXiv

Warped Circumbinary Disks in Active Galactic Nuclei

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.

preprint2013arXiv

Constraints on decreases in Eta Carinae's mass loss from 3D hydrodynamic simulations of its binary colliding winds

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.

preprint2013arXiv

Precessing Warped Be Disk Triggering the Giant Outbursts in 2009 and 2011 in A 0535+262/V725 Tau

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.

preprint2012arXiv

Origin of Two Types of X-Ray Outbursts in Be/X-Ray Binaries. I. Accretion Scenarios

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.

preprint2011arXiv

Constraining the Absolute Orientation of Eta Carinae's Binary Orbit: A 3-D Dynamical Model for the Broad [Fe III] Emission

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 θ 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, θ \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.

preprint2011arXiv

Drastic Spectroscopic Variability of the Be/X-ray Binary A0535+262/V725 Tau during and after the 2009 Giant Outburst

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.

preprint2011arXiv

Hydrodynamic Interaction between the Be Star and the Pulsar in the TeV Binary PSR B1259-63/LS 2883

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.

preprint2011arXiv

Resonant Excitation of Disk Oscillations in Deformed Disks IV: A New Formulation Studying Stability

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.

preprint2011arXiv

The first determination of the viscosity parameter in the circumstellar disk of a Be Star

Be stars possess gaseous circumstellar decretion disks, which are well described using standard $α$-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 $α=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 $α$ 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.

preprint2011arXiv

X-ray Modeling of η Carinae and WR140 from SPH Simulations

The colliding wind binary (CWB) systems η 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 β 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 η 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.

preprint2011arXiv

X-ray Modeling of η Carinae and WR140 from SPH Simulations

The colliding wind binary (CWB) systems η 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 η Car and WR140. This point-source model can also explain the early recovery of η Car's X-ray light curve from the 2009.0 minimum by a factor of 2-4 reduction in the mass loss rate of η 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 η 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.

preprint2010arXiv

Detection of high-velocity material from the wind-wind collision zone of Eta Carinae across the 2009.0 periastron passage

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.

preprint2010arXiv

Excitation of Trapped g-Mode Oscillations in Warped Disks around Black Holes

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.

preprint2010arXiv

The Dynamical Evolution of Be star disks

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.

preprint2008arXiv

Global Disk Oscillations in Binary Be Stars

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.

preprint2000arXiv

Viscous Transonic Decretion in Disks of Be Stars

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.