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.     

Microlensing of close binary stars

The gravity due to a multiple-mass system has a remarkable gravitational effect: the extreme magnification of background light sources along extended so-called caustic lines. This property has been the channel for some remarkable astrophysical discoveries over the past decade, including the detection and characterization of extrasolar planets, the routine analysis of limb darkening, and, in one case, limits set on the apparent shape of a star several kiloparsec distant. In this paper, we investigate the properties of the microlensing of close binary star systems. We show that in some cases it is possible to detect flux from the Roche lobes of close binary stars. Such observations could constrain models of close binary stellar systems.     

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.     

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.     

On the evolution of the cosmic supernova rates

The Galactic radio-emitting X-ray binary Cygnus X-3 is known to be a source of large-scale radio jets associated with periods of intense radio flaring. These jets have been found to have an expansion velocity of ∼0.3 c and are believed (on kinematic grounds) to lie close to the plane of the sky. We present new observations of Cygnus X-3 using the VLBA at 15 GHz. These observations, which included the detection of two small flares, show an additional kind of behaviour with apparent superluminal expansion along both major and minor axes. Evidence for superluminal activity has been found in a number of X-ray binary systems such as GRS 1915+105 and GRO J1655−40 with their superluminal radio jets. Apparently similar morphologies of the Galactic and extragalactic jet sources have led to the X-ray binaries being described as 'micro-quasars'. The superluminal expansion seen in our results appears to be different in nature from these other two sources, and a number of mechanisms are presented and discussed.     

Predicting the properties of binary stellar systems: the evolution of accreting protobinary systems

We investigate the formation of binary stellar systems. We consider a model where a 'seed' protobinary system forms, via fragmentation, within a collapsing molecular cloud core and evolves to its final mass by accreting material from an infalling gaseous envelope. This accretion alters the mass ratio and orbit of the binary, and is largely responsible for forming the circumstellar and/or circumbinary discs.
Given this model for binary formation, we predict the properties of binary systems and how they depend on the initial conditions within the molecular cloud core. We predict that there should be a continuous trend such that closer binaries are more likely to have equal-mass components and are more likely to have circumbinary discs than wider systems. Comparing our results with observations, we find that the observed mass-ratio distributions of binaries and the frequency of circumbinary discs as a function of separation are most easily reproduced if the progenitor molecular cloud cores have radial density profiles between uniform and 1/ r (e.g., Gaussian) with near-uniform rotation. This is in good agreement with the observed properties of pre-stellar cores. Conversely, we find that the observed properties of binaries cannot be reproduced if the cloud cores are in solid-body rotation and have initial density profiles which are strongly centrally condensed. Finally, in agreement with the radial-velocity searches for extrasolar planets, we find that it is very difficult to form a brown dwarf companion to a solar-type star with a separation ≲10 au, but that the frequency of brown dwarf companions should increase with larger separations or lower mass primaries.     

Numerical simulations of protostellar encounters — I. Star–disc encounters

It appears that most stars are born in clusters, and that at birth most stars have circumstellar discs which are comparable in size to the separations between the stars. Interactions between neighbouring stars and discs are therefore likely to play a key role in determining disc lifetimes, stellar masses, and the separations and eccentricities of binary orbits. Such interactions may also cause fragmentation of the discs, thereby triggering the formation of additional stars.   We have carried out a series of simulations of star–disc interactions using an SPH code which treats self-gravity, hydrodynamic and viscous forces. We find that interactions between discs and stars provide a mechanism for removing energy from, or adding energy to, the orbits of the stars, and for truncating the discs. However, capture during such encounters is unlikely to be an important binary formation mechanism.   A more significant consequence of such encounters is that they can trigger fragmentation of the disc, via tidally and compressionally induced gravitational instabilities, leading to the formation of additional stars and substellar objects. When the disc spins and stellar orbits are randomly oriented, encounters lead to the formation of new companions to the original star in 20 per cent of encounters. If most encounters are prograde and coplanar, as suggested by simulations of dynamically triggered star formation, then new companions are formed in approximately 50 per cent of encounters.     

Doppler Tomography

I review the method of Doppler tomography which is widely used to understand the complex emission line profile variations displayed by accreting binary stars. Doppler tomography uses the information encoded in line profiles as a function of orbital phase to calculate the strength of emission as a function of velocity, using a process closely related to medical X-ray tomographic imaging. I review applications which have revealed spiral structures in accretion discs, the accretion flows within magnetically-dominated systems and irradiation induced emission in X-ray binary stars. I also cover some of the more recent extensions to the method which variously allow for Roche geometry, modulation of the fluxes and motion at angles to the orbital plane.     

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.     

Observational implications of precessing protostellar discs and jets

We consider the dynamics of a protostellar disc in a binary system where the disc is misaligned with the orbital plane of the binary, with the aim of determining the observational consequences for such systems. The disc wobbles with a period approximately equal to half the orbital period of the binary and precesses on a longer time-scale. We determine the characteristic time-scale for realignment of the disc with the orbital plane as a result of dissipation. If the dissipation is determined by a simple isotropic viscosity then we find, in line with previous studies, that the alignment time-scale is of the order of the viscous evolution time-scale. However, for typical protostellar disc parameters, if the disc tilt exceeds the opening angle of the disc, then tidally induced shearing within the disc is transonic. In general, hydrodynamic instabilities associated with the internally driven shear result in extra dissipation that is expected to drastically reduce the alignment time-scale. For large disc tilts the alignment time-scale is then comparable with the precession time-scale, while for smaller tilt angles δ , the alignment time-scale varies as (sin δ )⁻¹. We discuss the consequences of the wobbling, precession and rapid realignment for observations of protostellar jets and the implications for binary star formation mechanisms.     

Accretion and dynamical interactions in small-N star-forming clusters: N= 5

We present results from high-resolution hydrodynamical simulations that explore the effects of small-scale clustering in star-forming regions. A large ensemble of small- N clusters with five stellar seeds have been modelled and the resulting properties of stars and brown dwarfs statistically derived and compared with observational data.
Close dynamical interactions between the protostars and competitive accretion driven by the cloud collapse are shown to produce a distribution of final masses that is bimodal, with most of the mass residing in the binary components. When convolved with a suitable core mass function, the final distribution of masses resembles the observed initial mass function, in both the stellar and substellar regimes. Binaries and single stars are found to constitute two kinematically distinct populations, with about half of the singles attaining velocities ≥2 km s⁻¹, which might deprive low-mass star-forming regions of their lightest members in a few crossing times. The eccentricity distribution of binaries and multiples is found to follow a distribution similar to that of observed long-period (uncircularized) binaries.
The results obtained support a mechanism in which a significant fraction of brown dwarfs form under similar circumstances as those of normal stars but are ejected from the common envelope of unstable multiple systems before their masses exceed the hydrogen burning limit. We predict that many close binary stars should have wide brown dwarf companions. Brown dwarfs, and, in general, very low-mass stars, would be rare as pure binary companions. The binary fraction should be a decreasing function of primary mass, with low-mass or substellar primaries being scarce. Where such binaries exist, they are expected either to be close enough (semimajor axis ∼10 au) to survive strong interactions with more massive binaries or to be born in very small molecular cloud cores.     

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⁻¹.     

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)     

A SAURON look at galaxy bulges

Kinematic and population studies show that bulges are generally rotationally flattened systems similar to lowluminosity ellipticals. However, observations with state‐of‐the‐art integral field spectrographs, such as SAURON, indicate that the situation is much more complex, and allow us to investigate phenomena such as triaxiality, kinematic decoupling and population substructure, and to study their connection to current formation and evolution scenarios for bulges of early‐type galaxies. We present the examples of two S0 bulges from galaxies in our sample of nearby galaxies: one that shows all the properties expected from classical bulges (NGC5866), and another case that presents kinematic features appropriate for barred disk galaxies (NGC7332). (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Three-dimensional smoothed particle hydrodynamics simulations of radiation-driven warped accretion discs

We present three-dimensional smoothed particle hydrodynamics calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc, resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and, in general, the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios, using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios   q ∼ 0.1  a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (cf. X 1916−053). For less extreme mass ratios, the entire disc tilts out of the orbital plane (cf. Her X–1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction, we obtain both positive and negative superhumps simultaneously in the dissipation light curve (cf. V603 Aql).     

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.     

General relativistic MHD simulations of black hole accretion disks and jets

Accretion disks orbiting black holes power high-energy systems such as X-ray binaries and Active Galactic Nuclei. Observations are providing increasingly detailed quantitative information about such systems. This data has been interpreted using standard toy-models that rely on simplifying assumptions such as regular flow geometry and a parameterized stress. Global numerical simulations offer a way to investigate the basic physical dynamics of accretion flows without these assumptions and, in principle, lead to a genuinely predictive theory. In recent years we have developed a fully three-dimensional general relativistic magnetohydrodynamic simulation code that evolves time-dependent inflows into Kerr black holes. Although the resulting global simulations of black hole accretion are still somewhat simplified, they have brought to light a number of interesting results. These include the formation of electro-magnetically dominated jets powered by the black hole’s rotation, and the presence of strong stresses in the plunging region of the accretion flow. The observational consequences of these features are gradually being examined. Increasing computer power and increasingly sophisticated algorithms promise a bright future for the computational approach to black hole accretion.     

Magnetized accretion flows: effects of gas pressure

We study how axisymmetric magnetohydrodynamic (MHD) accretion flows depend on γ adiabatic index in the polytropic equation of state. This work is an extension of Mościbrodzka & Proga, where we investigated the γ dependence of two-dimensional Bondi-like accretion flows in the hydrodynamical (HD) limit. Our main goal is to study if simulations for various γ can give us insights into the problem of various modes of accretion observed in several types of accretion systems, such as black hole binaries (BHBs), active galactic nuclei (AGN) and gamma-ray bursts. We find that for  γ≳ 4/3  , the fast-rotating flow forms a thick torus that is supported by rotation and gas pressure. As shown before for  γ= 5/3  , such a torus produces a strong, persistent bipolar outflow that can significantly reduce the polar funnel accretion of a slowly rotating flow. For low γ, close to 1, the torus is thin and is supported by rotation. The thin torus produces an unsteady outflow which is too weak to propagate throughout the polar funnel inflow. Compared to their HD counterparts, the MHD simulations show that the magnetized torus can produce an outflow and does not exhibit regular oscillations. Generally, our simulations demonstrate how the torus thickness affects the outflow production. They also support the notion that the geometrical thickness of the torus correlates with the power of the torus outflow. Our results, applied to observations, suggest that the torus ability to radiatively cool and become thin can correspond to a suppression of a jet as observed in the BHBs during a transition from a hard/low to soft/high spectral state and a transition from a quiescent to hard/low state in AGN.     

Mass loss from the inner regions of accretion discs due to centrifugally driven magnetic wind flows

Wind flows and collimated jets are believed to be a feature of a range of disc accreting systems. These include active galactic nuclei, T Tauri stars, X-ray binaries and cataclysmic variables. The observed collimation implies large-scale magnetic fields and it is known that dipole-symmetry fields of sufficient strength can channel wind flows emanating from the surfaces of a disc. The disc inflow leads to the bending of the poloidal magnetic field lines, and centrifugally driven magnetic winds can be launched when the bending exceeds a critical value. Such winds can result in angular momentum transport at least as effective as turbulent viscosity, and hence they can play a major part in driving the disc inflow.
It is shown here that if the standard boundary condition of vanishing viscous stress close to the stellar surface is applied, together with the standard connection between viscosity and magnetic diffusivity, then poloidal magnetic field bending increases as the star is approached with a corresponding increase in the wind mass loss rate. A significant amount of material can be lost from the system via the enhanced wind from a narrow region close to the stellar surface. This occurs for a Keplerian angular velocity distribution and for a modified form of angular velocity, which allows for matching of the disc and stellar rotation rates through a boundary layer above the stellar surface. The enhanced mass loss is significantly affected by the behaviour of the disc angular velocity as the stellar surface is approached, and hence by the stellar rotation rate. Such a mechanism may be related to the production of jets from the inner regions of disc accreting systems.