On the evolution of giant protoplanets forming in circumbinary discs

We present the results of hydrodynamic simulations of Jovian mass protoplanets that form in circumbinary discs. The simulations follow the orbital evolution of the binary plus protoplanet system acting under their mutual gravitational forces, and forces exerted by the viscous circumbinary disc. The evolution involves the clearing of the inner circumbinary disc initially, so that the binary plus protoplanet system orbits within a low density cavity. Continued interaction between disc and protoplanet causes inward migration of the planet towards the inner binary. Subsequent evolution can take three distinct paths: (i) the protoplanet enters the 4 : 1 mean motion resonance with the binary, but is gravitationally scattered through a close encounter with the secondary star; (ii) the protoplanet enters the 4 : 1 mean motion resonance, the resonance breaks, and the planet remains in a stable orbit just outside the resonance; (iii) when the binary has initial eccentricity   e _bin≥ 0.2  , the disc becomes eccentric, leading to a stalling of the planet migration, and the formation of a stable circumbinary planet.
These results have implications for a number of issues in the study of extrasolar planets. The ejection of protoplanets in close binary systems provides a source of 'free-floating planets', which have been discovered recently. The formation of a large, tidally truncated cavity may provide an observational signature of circumbinary planets during formation. The existence of protoplanets orbiting stably just outside a mean motion resonance (4 : 1) in the simulations indicate that such sites may harbour planets in binary star systems, and these could potentially be observed. Finally, the formation of stable circumbinary planets in eccentric binary systems indicates that circumbinary planets may not be uncommon.     

On the angular momentum transfer on to compact stars in binary systems

Results of three-dimensional numerical simulations of the gas transfer in close binary systems show that, in addition to the formation of a tidally induced spiral shock wave, it is also possible for accretion streams to be produced, having low specific angular momentum in a region close to the accreting star. These streams are mainly placed above the orbital disc but are also unevenly present in the equatorial plane. The relevance of such flows is related to formation of hot coronae or bulges in regions very close to the accretor centre. The eventual formation of such bulges and shock-heated flows is interesting in the context of advection-dominated solutions and for the explanation of spectral properties of the black hole candidates in binary systems.     

Supernova kicks and misaligned Be star binaries

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron star. They are found to have non-zero eccentricities and there is evidence that some systems have a misalignment between the spin axis of the star and the spin axis of the binary orbit. The eccentricities in these systems are caused by a kick to the neutron star during the supernova that formed it. Such kicks would also give rise to misalignments. In this paper, we investigate the extent to which the same kick distribution can give rise to both the observed eccentricity distribution and the observed misalignments. We find that a Maxwellian distribution of velocity kicks with a low velocity dispersion,  σ_k≈ 15 km s⁻¹  , is consistent with the observed eccentricity distribution but is hard to reconcile with the observed misalignments, typically   i ≥ 25°  . Alternatively, a higher velocity kick distribution,  σ_k= 265 km s⁻¹  , is consistent with the observed misalignments but not with the observed eccentricities, unless post-supernova circularization of the binary orbits has taken place. We discuss briefly how this might be achieved.     

Two- and three-dimensional numerical simulations of accretion discs in a close binary system

We perform 2D and 3D numerical simulations of an accretion disc in a close binary system using the simplified flux vector splitting (SFS) finite volume method. In our calculations, the gas is assumed to be ideal with γ =1.01, 1.05, 1.1 and 1.2 . The mass ratio of the mass-losing star to the mass-accreting star is unity. Our results show that spiral shocks are formed on the accretion disc in all cases. In 2D calculations we find that the smaller γ is, the more tightly the spiral winds. We observe this trend in 3D calculations as well in a somewhat weaker sense. Mach numbers in our discs are less than 10. These values are lower than the values in observed accretion discs in close binary systems.
Recently, Steeghs, Harlaftis & Horne found the first convincing evidence for spiral structure in the accretion disc of the eclipsing dwarf nova binary IP Pegasi, using the technique known as Doppler tomography. Although the Mach numbers in present calculations are rather low, we may claim that the spiral structure that we discovered in earlier numerical simulations is now found observationally.     

GRB060602B = Swift J1749.4−2807: an unusual transiently accreting neutron-star X-ray binary

We present an analysis of the Swift Burst Alert Telescope (BAT) and X-ray telescope (XRT) data of GRB060602B, which is most likely an accreting neutron star in a binary system and not a gamma-ray burst. Our analysis shows that the BAT burst spectrum is consistent with a thermonuclear flash (type I X-ray burst) from the surface of an accreting neutron star in a binary system. The X-ray binary nature is further confirmed by the report of a detection of a faint point source at the position of the XRT counterpart of the burst in archival XMM–Newton data approximately six year before the burst and in more recent XMM–Newton data obtained at the end of 2006 September (nearly four months after the burst). Since the source is very likely not a gamma-ray burst, we rename the source Swift J1749.4−2807, based on the Swift /BAT discovery coordinates. Using the BAT data of the type I X-ray burst, we determined that the source is at most at a distance of  6.7 ± 1.3 kpc  . For a transiently accreting X-ray binary, its soft X-ray behaviour is atypical: its 2–10 keV X-ray luminosity (as measured using the Swift /XRT data) decreased by nearly three orders of magnitude in about 1 day, much faster than what is usually seen for X-ray transients. If the earlier phases of the outburst also evolved this rapidly, then many similar systems might remain undiscovered because the X-rays are difficult to detect and the type I X-ray bursts might be missed by all the sky surveying instruments. This source might be part of a class of very fast transient low-mass X-ray binary systems of which there may be a significant population in our Galaxy.     

On the precession of accretion discs in X-ray binaries

In this Letter, recent results on the nodal precession of accretion discs in close binaries are applied to the discs in some X-ray binary systems. The ratio between the tidally forced precession period and the binary orbital period is given, as well as the condition required for the rigid precession of gaseous Keplerian discs. Hence the minimum precessional period that may be supported by a fluid Keplerian disc is determined. It is concluded that near-rigid body precession of tilted accretion discs can occur and generally reproduce observationally inferred precession periods, for reasonable system parameters. In particular, long periods in SS 433, Her X-1, LMC X-4 and SMC X-1 can be fitted by the tidal model. It is also found that the precession period that has been tentatively put forward for Cyg X-2 cannot be accommodated by a tidally precessing disc model for any realistic choice of system parameters.     

Stellar-mass black hole binaries as ultraluminous X-ray sources

Ultraluminous X-ray sources (ULXs) with   L _x > 10³⁹ erg s⁻¹  have been discovered in great numbers in external galaxies with ROSAT , Chandra and XMM-Newton . The central question regarding this important class of sources is whether they represent an extension in the luminosity function of binary X-ray sources containing neutron stars and stellar-mass black holes (BHs), or a new class of objects, e.g. systems containing intermediate-mass BHs  (100–1000 M_⊙)  . We have carried out a theoretical study to test whether a large fraction of the ULXs, especially those in galaxies with recent star formation activity, can be explained with binary systems containing stellar-mass BHs. To this end, we have applied a unique set of binary evolution models for BH X-ray binaries, coupled to a binary population synthesis code, to model the ULXs observed in external galaxies. We find that for donor stars with initial masses  ≳10 M_⊙  the mass transfer driven by the normal nuclear evolution of the donor star is sufficient to potentially power most ULXs. This is the case during core hydrogen burning and, to an even more pronounced degree, while the donor star ascends the giant branch, although the latter phases last only ∼5 per cent of the main-sequence phase. We show that with only a modest violation of the Eddington limit, e.g. a factor of ∼10, both the numbers and properties of the majority of the ULXs can be reproduced. One of our conclusions is that if stellar-mass BH binaries account for a significant fraction of ULXs in star-forming galaxies, then the rate of formation of such systems is  ∼3 × 10⁻⁷ yr⁻¹  normalized to a core-collapse supernova rate of 0.01 yr⁻¹.     

Alignment time-scale of the microquasar GRO J1655−40

The microquasar GRO J1655−40 has a black hole with spin angular momentum apparently misaligned to the orbital plane of its companion star. We analytically model the system with a steady-state disc warped by Lense–Thirring precession and find the time-scale for the alignment of the black hole with the binary orbit. We make detailed stellar evolution models so as to estimate the accretion rate and the lifetime of the system in this state. The secondary can be evolving at the end of the main sequence or across the Hertzsprung gap. The mass-transfer rate is typically 50 times higher in the latter case but we find that, in both the cases, the lifetime of the mass-transfer state is at most a few times the alignment time-scale. The fact that the black hole has not yet aligned with the orbital plane is therefore consistent with either model. We conclude that the system may or may not have been counter aligned after its supernova kick but that it is most likely to be close to alignment rather than counter alignment now.     

Neutron star high‐mass binaries as the origin of SGR/AXP

A close high‐mass binary system consisting of a neutron star (NS) and a massive OB supergiant companion is expected to lead to a Thorne‐Żytkow object (TZO) structure, which consists of a NS core and a stellar envelope. We use the scenario machine program to calculate the formation tracks of TZOs in close high‐mass NS binaries and their subsequent evolution. We propose and demonstrate that the explosion and instant contraction of a TZO structure leave its stellar remnant as a soft gamma‐ray repeater and an anomalous X‐ray pulsar respectively. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Close encounters in young stellar clusters: implications for planetary systems in the solar neighbourhood

The stars that populate the solar neighbourhood were formed in stellar clusters. Through N -body simulations of these clusters, we measure the rate of close encounters between stars. By monitoring the interaction histories of each star, we investigate the singleton fraction in the solar neighbourhood. A singleton is a star which formed as a single star, has never experienced any close encounters with other stars or binaries, or undergone an exchange encounter with a binary. We find that, of the stars which formed as single stars, a significant fraction is not singletons once the clusters have dispersed. If some of these stars had planetary systems, with properties similar to those of the Solar System, the planets' orbits may have been perturbed by the effects of close encounters with other stars or the effects of a companion star within a binary. Such perturbations can lead to strong planet–planet interactions which eject several planets, leaving the remaining planets on eccentric orbits. Some of the single stars exchange into binaries. Most of these binaries are broken up via subsequent interactions within the cluster, but some remain intact beyond the lifetime of the cluster. The properties of these binaries are similar to those of the observed binary systems containing extrasolar planets. Thus, dynamical processes in young stellar clusters will alter significantly any population of Solar System-like planetary systems. In addition, beginning with a population of planetary systems exactly resembling the Solar System around single stars, dynamical encounters in young stellar clusters may produce at least some of the extrasolar planetary systems observed in the solar neighbourhood.     

Accretion disc in the massive V448 Cygni system

The UBV light curves of the early-type eclipsing binary V448 Cygni, obtained at the Abastumani Astrophysical Observatory from 1964 to 1967, are re-analysed here. The analysis was made assuming the presence of an accretion disc in the system, as inferred from the light-curve shape and spectroscopic characteristics of the system. The Roche model of a binary was used, containing a geometrically and optically thick accretion disc around the hotter and more massive star. By solving the inverse problem, the orbital elements and the physical parameters of the system components and of the accretion disc were estimated. This result is important for understanding the star formation and evolution processes in the systems with massive components.     

A possible detection mechanism for pre-low-mass X-ray binaries

This work presents a possible detection mechanism for close, detached, neutron star–red dwarf binaries, which are expected to be the evolutionary precursors of low-mass X-ray binaries (LMXBs). Although this pre-low-mass X-ray binary (pre-LMXB) phase of evolution is predicted theoretically, as yet no such systems have been identified observationally. The calculations presented here suggest that the X-ray luminosity of neutron star wind accretion in a pre-LMXB system can be expected to exceed the intrinsic X-ray luminosity of the red dwarf secondary star. Furthermore, the temperature of the radiation emitted from the neutron star wind accretion process is expected, within the confines of a reasonable set of conditions, to lie within the detection range of X-ray satellites. Sources with X-ray luminosities greater than that expected for a red dwarf star, but the positions of which coincide with that of a red dwarf star, are then candidate pre-LMXB systems. These candidate systems should be surveyed for the radial velocity shifts that would occur as a result of the orbital motion of a red dwarf star within a close binary system containing a high-mass compact object.     

Influence of the r-mode instability on hypercritically accreting neutron stars

We have investigated the influence of the r-mode instability on hypercritically accreting neutron stars in close binary systems during their common envelope phases, based on the scenario proposed by Brown et al. On the one hand, neutron stars are heated by the accreted matter at the stellar surface, but on the other hand they are also cooled down by the neutrino radiation. At the same time, the accreted matter transports its angular momentum and mass to the star. We have studied the evolution of the stellar mass, temperature and rotational frequency.
The gravitational-wave-driven instability of the r-mode oscillation strongly suppresses spinning up of the star, the final rotational frequency of which is well below the mass-shedding limit, in fact typically as low as 10 per cent of that of the mass-shedding state. On a very short time-scale the rotational frequency tends to approach a certain constant value and saturates there, as long as the amount of accreted mass does not exceed a certain limit to collapse to a black hole. This implies that a similar mechanism of gravitational radiation to that in the so-called 'Wagoner star' may work in this process. The star is spun up by accretion until the angular momentum loss by gravitational radiation balances the accretion torque. The time-integrated dimensionless strain of the radiated gravitational wave may be large enough to be detectable by gravitational wave detectors such as LIGO II.     

Supersonic wind accretion in massive close binary systems: a 3D SPH simulation of Cen X-3

In this paper we study the 3D SPH structure and dynamics of an accretion disc generated in a close binary system by supersonic wind accretion from a massive secondary on to a compact primary. The stellar masses and separation between the two components are characteristic of the Cen X-3 system: the secondary is a 19.1-M⊙ star not filling completely its Roche lobe, while the primary is a white dwarf or a neutron star of 1.4 M⊙.
An interesting result of our simulation is that, in a quasi-stationary state attained after ≃4 orbital periods, only about three-quarters of the particles released by the secondary penetrate the primary Roche lobe. The disc is remarkably elongated and thick, and consistent deviations from the 'standard model' of specific angular momentum and radial temperature distributions have been found. However, the most interesting result is that the azimuthal distribution of the radial Mach number shows oblique structures (spiral shocks), which persist from the outer edge to the inner regions, in contrast to the case of discs formed by the usual L1 accretion and even by wind accretion in much less massive systems.     

XMM–Newton observations of UW CrB: detection of X-ray bursts and evidence for accretion disc evolution

UW CrB (MS 1603+2600) is a peculiar short-period X-ray binary that exhibits extraordinary optical behaviour. The shape of the optical light curve of the system changes drastically from night to night, without any changes in overall brightness. Here we report X-ray observations of UW CrB obtained with XMM–Newton . We find evidence for several X-ray bursts, confirming a neutron star primary. This considerably strengthens the case that UW CrB is an accretion disc corona system located at a distance of at least 5–7 kpc (3–5 kpc above the Galactic plane). The X-ray and Optical Monitor (ultraviolet–optical) light curves show remarkable shape variation from one observing run to another, which we suggest are due to large-scale variations in the accretion disc shape resulting from a warp that periodically obscures the optical and soft X-ray emission. This is also supported by the changes in phase-resolved X-ray spectra.     

Super-orbital period in the high-mass X-ray binary 2S 0114+650

We report the detection of a stable super-orbital period in the high-mass X-ray binary 2S 0114+650. Analyses of data from the Rossi X-ray Timing Explorer All-Sky Monitor from 1996 January 5 to 2004 August 25 reveal a super-orbital period of 30.7±0.1 d, in addition to confirming the previously reported neutron star spin period of 2.7 h and the binary orbital period of 11.6 d. It is unclear if the super-orbital period can be ascribed to the precession of a warped accretion disc in the system.     

The impact of a supernova explosion in a very massive binary

We consider the effect of a supernova (SN) explosion in a very massive binary that is expected to form in a portion of Population III stars with the mass higher than  100 M_⊙  . In a Population III binary system, a more massive star can result in the formation of a black hole (BH) and a surrounding accretion disc. Such BH accretion could be a significant source of the cosmic reionization in the early Universe. However, a less massive companion star evolves belatedly and eventually undergoes a SN explosion, so that the accretion disc around a BH might be blown off in a lifetime of companion star. In this paper, we explore the dynamical impact of a SN explosion on an accretion disc around a massive BH, and elucidate whether the BH accretion disc is totally demolished or not. For the purpose, we perform three-dimensional hydrodynamic simulations of a very massive binary system, where we assume a BH of  10³ M_⊙  that results from a direct collapse of a very massive star and a companion star of  100 M_⊙  that undergoes a SN explosion. We calculate the remaining mass of a BH accretion disc as a function of time. As a result, it is found that a significant portion of gas disc can survive through three-dimensional geometrical effects even after the SN explosion of a companion star. Even if the SN explosion energy is higher by two orders of magnitude than the binding energy of gas disc, about a half of disc can be left over. The results imply that the Population III BH accretion disc can be a long-lived luminous source, and therefore could be an important ionizing source in the early Universe.     

Optical spectroscopy of the quiescent counterpart to EXO 0748−676

We present phase resolved optical spectroscopy and X-ray timing of the neutron star X-ray binary EXO 0748−676 after the source returned to quiescence in the autumn of 2008. The X-ray light curve displays eclipses consistent in orbital period, orbital phase and duration with the predictions and measurements before the return to quiescence. Hα and He  i emission lines are present in the optical spectra and show the signature of the orbit of the binary companion, placing a lower limit on the radial velocity semi-amplitude of   K ₂ > 405 km s⁻¹  . Both the flux in the continuum and the emission lines show orbital modulations, indicating that we observe the hemisphere of the binary companion that is being irradiated by the neutron star. Effects due to this irradiation preclude a direct measurement of the radial velocity semi-amplitude of the binary companion; in fact, no stellar absorption lines are seen in the spectrum. Nevertheless, our observations place a stringent lower limit on the neutron star mass of   M ₁ > 1.27 M_⊙  . For the canonical neutron star mass of   M ₁= 1.4 M_⊙  , the mass ratio is constrained to  0.075 < q < 0.105  .     

Accretion disc dynamics in extreme mass ratio compact binaries

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios   q = M ₂/ M ₁ < 0.1  , applicable to ultracompact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For   q ∼ 0.1  , the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q . At   q = 0.01  , the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com, which has   q = 0.02  . The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultracompact binaries.     

A model of AW UMa

The contact binary AW UMa has an extreme mass ratio, with the more-massive component (the current primary) close to the main sequence, while the low-mass star at   q ≈ 0.1  (the current secondary) has a much larger radius than a main-sequence star of a comparable mass. We propose that the current secondary has almost exhausted hydrogen in its centre and is much more advanced in its evolution, as suggested by Stpień. Presumably the current secondary lost most of its mass during its evolution with part of it transferred to the current primary. After losing a large fraction of its angular momentum, the binary may evolve into a system of FK Com type.