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.     

The dynamics of eccentric accretion discs in superhump systems

We have applied an eccentric accretion disc theory in simplified form to the case of an accretion disc in a binary system, where the disc contains the 3:1 Lindblad resonance. This is relevant to the case of superhumps in SU Ursae Majoris cataclysmic variables and other systems, where it is thought that this resonance leads to growth of eccentricity and a modulation in the light curve due to the interaction of a precessing eccentric disc with tidal stresses. A single differential equation is formulated which describes the propagation, resonant excitation and viscous damping of eccentricity. The theory is first worked out in the simple case of a narrow ring and leads to the conclusion that the eccentricity distribution is locally suppressed by the presence of the resonance, creating a dip in the eccentricity at the resonant radius. Application of this theory to the superhump case confirms this conclusion and produces a more accurate expression for the precession rate of the disc than has been previously accomplished with simple dynamical estimates.     

Tilted accretion discs in cataclysmic variables: tidal instabilities and superhumps

We investigate the growth of tidal instabilities in accretion discs in a binary star potential, using three-dimensional numerical simulations. As expected from analytic work, the disc is prone to an eccentric instability provided that it is large enough to extend to the 3:1 resonance. The eccentric disc leads to positive superhumps in the light curve. It has been proposed that negative superhumps might arise from a tilted disc, but we find no evidence that the companion gravitational tilt instability can grow fast enough in a fluid disc to create a measurable inclination. The origin of negative superhumps in the light curves of cataclysmic variables remains a puzzle.     

A spectroscopic study of V603 Aquilae: stellar parameters and continuumline variations

A time-resolved spectroscopic study of V603 Aql (Nova Aquilae 1918) is presented. An orbital period of P _orb=01385±00002, consistent with previous results, and a radial velocity semi-amplitude of K =20±3 km s¹ are obtained from the radial velocity variations of the H emission line. Similar K values are also found in H , H , and He  i emission lines. Using the measured FWHM of the H line and assuming that the derived semi-amplitude is that of the white dwarf, we deduce a most likely mass ratio of q =0.24±0.05 and stellar masses of M ₂=0.29±0.04 M and M ₁=1.2±0.2 M for the secondary and primary (the white dwarf) star, respectively. The dynamical solution also indicates a very low orbital inclination, i =13°±2°. We find that the continuum and line variations are modulated with both the positive and the negative superhump periods, indicating that they arise from similar regions of the accretion disc. Moreover, we find, for the first time from spectroscopy, evidence of negative superhumps in addition to the positive superhumps. Positive superhumps are explained within the disc instability model as caused by an eccentric disc surrounding the white dwarf, which is precessing (apsidal advance) because of tidal instabilities, causing the observed positive superhumps. A nodal precession in the accretion disc is currently believed to be the cause of the observed negative superhumps. The low value of q is consistent with the expected value for systems that show superhumps, in accordance with the eccentric disc model. We find no evidence of periodicity associated with the spin period.     

Investigating the structure of the accretion disc in WZ Sge from multiwaveband time-resolved spectroscopic observations – I

We present the first of two papers describing an in-depth study of multiwaveband phase-resolved spectroscopy of the unusual dwarf nova WZ Sge. In this paper we present an extensive set of Doppler maps of WZ Sge covering optical and infrared emission lines, and describe a new technique for studying the accretion discs of cataclysmic variables using ratioed Doppler maps. Applying the ratioed Doppler map technique to our WZ Sge data shows that the radial temperature profile of the disc is unlike that predicted for a steady state α disc. Time-averaged spectra of the accretion disc line flux (with the bright spot contribution removed) show evidence in the shapes of the line profiles for the presence of shear broadening in a quiescent non-turbulent accretion disc. From the positions of the bright spots in the Doppler maps of different lines, we conclude that the bright spot region is elongated along the ballistic stream, and that the density of the outer disc is low. The velocity of the outer edge of the accretion disc measured from the H α line is found to be 723±23 km s⁻¹. Assuming that the accretion disc reaches to the 3:1 tidal resonance radius, we derive a value for the primary star mass of 0.82 M_⊙. We discuss the implications of our results on the present theories of WZ Sge type dwarf nova outbursts.     

Simulations of superhumps and superoutbursts

We numerically study the tidal instability of accretion discs in close binary systems using a two-dimensional SPH code. We find that the precession rate of tidally unstable, eccentric discs does not only depend upon the binary mass ratio q . Although the (prograde) disc precession rate increases with the strength of the tidal potential, we find that increasing the shear viscosity ν also has a significant prograde effect. Increasing the disc temperature has a retrograde impact upon the precession rate.   We find that motion relative to the binary potential results in superhump-like, periodic luminosity variations in the outer reaches of an eccentric disc. The nature and location of the luminosity modulation are functions of ν. Light curves most similar to observations are obtained for ν values appropriate for a dwarf nova in outburst.   We investigate the thermal–tidal instability model for superoutburst. A dwarf nova outburst is simulated by instantaneously increasing ν, which causes a rapid radial expansion of the disc. Should the disc encounter the 3: 1 eccentric inner Lindblad resonance and become tidally unstable, then tidal torques become much more efficient at removing angular momentum from the disc. The disc then shrinks and M _d increases. The resulting increase in disc luminosity is found to be consistent with the excess luminosity of a superoutburst.     

Reprocessing of X-rays in the outer accretion disc of the black hole binary XTE J1817−330

We build a simple model of the optical/ultraviolet (UV) emission from irradiation of the outer disc by the inner disc and coronal emission in black hole binaries. We apply this to the broad-band Swift data from the outburst of the black hole binary XTE J1817−330 to confirm previous results that the optical/UV emission in the soft state is consistent with a reprocessing a constant fraction of the bolometric X-ray luminosity. However, this is very surprising as the disc temperature drops by more than a factor of 3 in the soft state, which should produce a marked change in the reprocessing efficiency. The easiest way to match the observed constant reprocessed fraction is for the disc skin to be highly ionized (as suggested 30 yr ago by van Paradijs), so that the bulk of the disc flux is reflected and only the hardest X-rays heat the disc. The constant reprocessed fraction also favours direct illumination of the disc over a scattering origin as the optical depth/solid angle of any scattering material (wind/corona) over the disc should decrease as the source luminosity declines. By contrast, the reprocessed fraction increases very significantly (by a factor of ∼6) as the source enters the hard state. This dramatic change is not evident from X-ray/UV flux correlations as it is masked by bandpass effects. However, it does not necessarily signal a change in emission, for example, the emergence of the jet dominating the optical/UV flux as the reflection albedo must change with the dramatic change in spectral shape.     

Disc precession in the M31 dipping X-ray binary Bo 158?

We present results from three XMM–Newton observations of the M31 low mass X-ray binary (LMXB) XMMU J004314.4+410726.3 (Bo 158), spaced over 3 d in 2004 July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-h period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anticorrelated with source intensity. In light of the 2004 XMM–Newton observations of Bo 158, we suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio ≲0.3 (period ≲4 h), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of  ∼11 P _orb  , and prograde disc precession on a period of  29 ± 1 P _orb  . This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on an  81 ± 3 h  cycle.     

Superhumps in V348 Pup

The eclipsing nova-like cataclysmic variable star V348 Pup exhibits a persistent luminosity modulation with a period 6 per cent longer than its 2.44-h orbital period ( P _orb). This has been interpreted as a 'positive superhump' resulting from a slowly precessing non-axisymmetric accretion disc gravitationally interacting with the secondary. We find a clear modulation of mid-eclipse times on the superhump period, which agrees well with the predictions of a simple precessing eccentric disc model. Our modelling shows that the disc light centre is on the far side of the disc from the donor star when the superhump reaches maximum light. This phasing suggests a link between superhumps in V348 Pup and late superhumps in SU UMa systems. Modelling of the full light curve and maximum entropy eclipse mapping both show that the disc emission is concentrated closer to the white dwarf at superhump maximum than at superhump minimum. We detect additional signals consistent with the beat periods between the implied disc precession period and both and     

Superhumps in low-mass X-ray binaries

We propose a mechanism for the superhump modulations observed in optical photometry of at least two black-hole X-ray transients (SXTs). As in extreme mass-ratio cataclysmic variables (CVs), superhumps are assumed to result from the presence of the 3:1 orbital resonance in the accretion disc. This causes the disc to become non-axisymmetric and precess. However, the mechanism for superhump luminosity variations in low-mass X-ray binaries (LMXBs) must differ from that in CVs, where it is attributed to a tidally-driven modulation of the disc's viscous dissipation, varying on the beat between the orbital and disc precession period. By contrast in LMXBs, tidal dissipation in the outer accretion disc is negligible: the optical emission is overwhelmingly dominated by reprocessing of intercepted central X-rays. Thus a different origin for the superhump modulation is required. Recent observations and numerical simulations indicate that in an extreme mass-ratio system the disc area changes on the superhump period. We deduce that the superhumps observed in SXTs arise from a modulation of the reprocessed flux by the changing area. Therefore, unlike the situation in CVs, where the superhump amplitude is inclination-independent, superhumps should be best seen in low-inclination LMXBs, whereas an orbital modulation from the heated face of the secondary star should be more prominent at high inclinations. Modulation at the disc precession period (10 s of days) may indicate disc asymmetries such as warping. We comment on the orbital period determinations of LMXBs, and the possibility and significance of possible permanent superhump LMXBs.     

Magnetically moderated outbursts of WZ Sagittae

We argue that the quiescent value of the viscosity parameter of the accretion disc in WZ Sge may be  α_cold∼ 0.01  , in agreement with estimates of α_cold for other dwarf novae. Assuming the white dwarf in WZ Sge to be magnetic, we show that, in quiescence, material close to the white dwarf can be propelled to larger radii, depleting the inner accretion disc. The propeller therefore has the effect of stabilizing the inner disc and allowing the outer disc to accumulate mass. The outbursts of WZ Sge are then regulated by the (magnetically determined) evolution of the surface density of the outer disc at a radius close to the tidal limit. Numerical models confirm that the recurrence time can be significantly extended in this way. The outbursts are expected to be superoutbursts since the outer disc radius is forced to exceed the tidal (3:1 resonance) radius. The large, quiescent disc is expected to be massive, and to be able to supply the observed mass accretion rate during outburst. We predict that the long-term spin evolution of the white dwarf spin will involve a long cycle of spin-up and spin-down phases.     

Launching of conical winds and axial jets from the disc–magnetosphere boundary: axisymmetric and 3D simulations

We investigate the launching of outflows from the disc–magnetosphere boundary of slowly and rapidly rotating magnetized stars using axisymmetric and exploratory 3D magnetohydrodynamic simulations. We find long-lasting outflows in the following cases. (1) In the case of slowly rotating stars , a new type of outflow, a conical wind , is found and studied in simulations. The conical winds appear in cases where the magnetic flux of the star is bunched up by the disc into an X-type configuration. The winds have the shape of a thin conical shell with a half-opening angle  θ∼ 30°–40°  . About 10–30 per cent of the disc matter flows from the inner disc into the conical winds. The conical winds may be responsible for episodic as well as long-lasting outflows in different types of stars. There is also a low-density, higher velocity component (a jet) in the region inside the conical wind. (2) In the case of rapidly rotating stars (the 'propeller regime'), a two-component outflow is observed. One component is similar to the conical winds. A significant fraction of the disc matter may be ejected into the winds. The second component is a high-velocity, low-density magnetically dominated axial jet where matter flows along the opened polar field lines of the star. The jet has a mass flux of about 10 per cent of that of the conical wind, but its energy flux (dominantly magnetic) can be larger than the energy flux of the conical wind. The jet's angular momentum flux (also dominantly magnetic) causes the star to spin down rapidly. Propeller-driven outflows may be responsible for the jets in protostars and for their rapid spin-down. The jet is collimated by the magnetic force while the conical winds are only weakly collimated in the simulation region. Exploratory 3D simulations show that conical winds are axisymmetric about the rotational axis (of the star and the disc), even when the dipole field of the star is significantly misaligned.     

Mass transfer in tidally unstable compact binaries

The 2001 outburst of WZ Sagittae has shown the most compelling evidence yet for an enhancement of the mass-transfer rate from the donor star during a dwarf nova outburst in the form of hotspot brightening. I show that, even in this extreme case, the brightening can be attributed to tidal heating near the interaction point of an accretion stream with the expanding edge of an eccentric accretion disc, with no need at all for an increase in the mass-transfer rate. Furthermore, I confirm previous suggestions that an increase in mass-transfer rate through the stream damps any eccentricity in an accretion disc and suppresses the appearance of superhumps, in contradiction to observations. Tidal heating is expected to be most significant in systems with small mass ratios. It follows that systems like WZ Sagittae – which has a tiny mass ratio – are those most likely to show a brightening in the hotspot region.     

Substellar companions and isolated planetary-mass objects from protostellar disc fragmentation

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time-scale that is short compared with the orbital period. We use a combination of hydrodynamic simulations and N -body orbit integrations to study the long-term evolution of a fragmenting disc with an initial mass ratio to the star of   M _disc/ M _*= 0.1  . For a disc that is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar-mass star) up to  ≈0.01 M_⊙  . Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 – in which the companions are close to or beyond the deuterium burning limit – appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ∼1 au should differ from that of stars with lower mass planetary companions.     

Optical photometry and spectroscopy of the accretion-powered millisecond pulsar HETE J1900.1 − 2455

We present phase resolved optical photometry and spectroscopy of the accreting millisecond pulsar HETE J1900.1−2455. Our R -band light curves exhibit a sinusoidal modulation, at close to the orbital period, which we initially attributed to X-ray heating of the irradiated face of the secondary star. However, further analysis reveals that the source of the modulation is more likely due to superhumps caused by a precessing accretion disc. Doppler tomography of a broad Hα emission line reveals an emission ring, consistent with that expected from an accretion disc. Using the velocity of the emission ring as an estimate for the projected outer disc velocity, we constrain the maximum projected velocity of the secondary to be 200 km s⁻¹, placing a lower limit of  0.05 M_⊙  on the secondary mass. For a  1.4 M_⊙  primary, this implies that the orbital inclination is low, ≲20°. Utilizing the observed relationship between the secondary mass and the orbital period in short-period cataclysmic variables, we estimate the secondary mass to be ∼0.085  M_⊙  , which implies an upper limit of  ∼2.4 M_⊙  for the primary mass.     

Black hole mergers: can gas discs solve the 'final parsec' problem?

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre-existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework.
The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high-photon-energy variability on the relatively long orbital time-scale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary.
We also find that the disc can drive the binary to merger on a reasonable time-scale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e.   a ₀≲ 0.05  pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is ≲0.001. Gas discs thus appear not to be an effective solution to the 'last parsec' problem for a significant class of mergers.     

Optical spectroscopy and Doppler tomography of Cygnus X-2

We present phase resolved optical spectroscopy and Doppler tomography of V1341 Cygni, the optical counterpart to the neutron star low-mass X-ray binary (LMXB) Cygnus X-2 (Cyg X-2). We derive a radial velocity (RV) curve for the secondary star, finding a projected RV semi-amplitude of   K ₂= 79 ± 3 km s⁻¹  , leading to a mass function of  0.51 ± 0.06 M_⊙, ∼30  per cent lower than the previous estimate. We tentatively attribute the lower value of K ₂ (compared to that obtained by other authors) to variations in the X-ray irradiation of the secondary star at different epochs of observations. The limited phase coverage and/or longer timebase of previous observations may also contribute to the difference in K ₂. Our value for the mass function implies a primary mass of  1.5 ± 0.3 M_⊙  , somewhat lower than previous dynamical estimates, but consistent with the value found by analysis of type-I X-ray bursts from this system. Our Doppler tomography of the broad He  ii λ4686 line reveals that most of the emission from this line is produced on the irradiated face of the donor star, with little emission from the accretion disc. In contrast, the Doppler tomogram of the N  iii λ4640.64 Bowen blend line shows bright emission from near the gas stream/accretion disc impact region, with fainter emission from the gas stream and secondary star. This is the first LMXB for which the Bowen blend is dominated by emission from the gas stream/accretion disc impact region, without comparable emission from the secondary star. This has implications for the interpretation of Bowen blend Doppler tomograms of other LMXBs for which the ephemeris may not be accurately known.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Superhumps: confronting theory with observation

We review the theory and observations related to the 'superhump' precession of eccentric accretion discs in close binary systems. We agree with earlier work, although for different reasons, that the discrepancy between observation and dynamical theory implies that the effect of pressure in the disc cannot be neglected. We extend earlier work that investigates this effect to include the correct expression for the radius at which resonant orbits occur. Using analytic expressions for the accretion disc structure, we derive a relationship between the period excess and mass ratio with the pressure effects included. This is compared to the observed data, recently derived results for detailed integration of the disc equations and the equivalent empirically derived relations and used to predict values for the mass ratio based on measured values of the period excess for 88 systems.     

Eccentric instability of a two-dimensional fluid disc with a central massive object

In this paper we investigate, by linear modal analysis, the one-armed dynamical instability of a two-dimensional fluid disc that has a massive object at its centre. The model of the disc is chosen to avoid the artificial instabilities that originate from the unrealistic disc configurations that have been adopted in previous studies. We find a one-armed instability for which the central massive object is displaced from the centre, which is generally called the 'eccentric instability'. However, to excite the eccentric instability, the mass of the central object should be appreciably smaller than that of the disc, and this mass ratio is far smaller than what was originally proposed. The instability shown in this paper is likely to be excited in a stellar system with a central massive object, e.g. a galactic nucleus harbouring a massive black hole, and further studies are desirable via techniques such as numerical simulations.