Prof Steven Balbus FRS, FInstP

Abstract:

About ten per cent of 'massive' stars (those of more than 1.5 solar masses) have strong, large-scale surface magnetic fields1-3. It has been suggested that merging of main-sequence and pre-main-sequence stars could produce such strong fields4,5, and the predicted fraction of merged massive stars is also about ten per cent6,7. The merger hypothesis is further supported by a lack of magnetic stars in close binaries8,9, which is as expected if mergers produce magnetic stars. Here we report three-dimensional magnetohydrodynamical simulations of the coalescence of two massive stars and follow the evolution of the merged product. Strong magnetic fields are produced in the simulations, and the merged star rejuvenates such that it appears younger and bluer than other coeval stars. This can explain the properties of the magnetic 'blue straggler' star τ Sco in the Upper Scorpius association that has an observationally inferred, apparent age of less than five million years, which is less than half the age of its birth association10. Such massive blue straggler stars seem likely to be progenitors of magnetars, perhaps giving rise to some of the enigmatic fast radio bursts observed11, and their supernovae may be affected by their strong magnetic fields12.

Abstract:

We present solutions to the relativistic thin disc evolutionary equation using an α-model for the turbulent stress tensor. Solutions with a finite stress at the innermost stable circular orbit (ISCO) give rise to bolometric light curves with a shallow power-law time dependence, in good agreement with those observed in tidal disruption events. A self-similar model based on electron scattering opacity, for example, yields a power-law index of −11/14, as opposed to −19/16 for the case of zero ISCO stress. These solutions correspond to an extended period of relaxation of the evolving disc which, like the light curves they produce, is not sustainable indefinitely. Cumulative departures from the approximation of exact circular orbits cause the power-law index to evolve slowly with time, leading eventually to the steeper fall-off associated with traditional zero ISCO stress models. These modified solutions are discussed in detail in a companion paper.

Abstract:

We present solutions to the relativistic thin disc evolutionary equation using a modified description of the mean fluid flow within the disc. The model takes into account the effects of sub-circular velocities in the innermost disc regions, and resolves otherwise unsustainable behaviour present in simple finite innermost stable circular orbit (ISCO) stress disc models. We show that the behaviour of a relativistic thin disc evolving with a finite ISCO stress is comprised of three distinct stages which join the ordinarily distinct finite and vanishing ISCO stress solutions into a fully continuous model parametrization. The most important prediction of our model is the existence of an intermediate stage of ‘stalled accretion’, controlled by a single dimensionless parameter. The hallmarks of this evolutionary phase appear to have been seen in General Relativistic MHD simulations as well as in the late time X-ray observations of tidal disruption events, but dedicated simulations and extended observations are needed for a deeper understanding.

Abstract:

An encounter between a passing star and a massive black hole at the centre of a galaxy, a so-called tidal disruption event or TDE, may leave a debris disc that subsequently accretes onto the hole. We solve for the time evolution of such a TDE disc, making use of an evolutionary equation valid for both the Newtonian and Kerr regimes. The late time luminosity emergent from such a disc is of interest as a model diagnostic, as it tends to follow a power law decline. The original simple ballistic fallback model, with equal mass in equal energy intervals, produces a −5/3 power law, while standard viscous disc descriptions yield a somewhat more shallow decline, with an index closer to −1.2. Of four recent, well-observed tidal disruption event candidates however, all had fall-off power law indices smaller than 1 in magnitude. In this work, we revisit the problem of thin disc evolution, solving this reduced problem in full general relativity. Our solutions produce power law indices that are in much better accord with observations. The late time observational data from many TDEs are generally supportive, not only of disc accretion models, but of finite stress persisting down to the innermost stable circular orbit.

Abstract:

Understanding what determines the strength of MHD turbulence in accretion discs is a question of fundamental theoretical and observational importance. In this work we investigate whether the dependence of the turbulent accretion disc stress (α) on the magnetic Prandtl number (Pm) is sufficiently sensitive to induce thermal-viscous instability using 3D MHD simulations. We first investigate whether the α-Pm dependence, found by many previous authors, has a physical or numerical origin by conducting a suite of local shearing-box simulations. We find that a definite α-Pm dependence persists when simultaneously increasing numerical resolution and decreasing the absolute values of both the viscous and resistive dissipation coefficients. This points to a physical origin of the α-Pm dependence. Using a further set of simulations which include realistic turbulent heating and radiative cooling, and by giving Pm a realistic physical dependence on the plasma temperature and density, we demonstrate that the α-Pm dependence is sufficiently strong to lead to a local instability. We confirm that the instability manifests itself as an unstable limit cycle by mapping the local thermal-equilibrium curve of the disc. This is the first self-consistent MHD simulation demonstrating the Pm instability from first principles. This result is important because a physical Pm instability could lead to the global propagation of heating and cooling fronts and a transition between disc states on timescales compatible with the observed hard/soft state transitions in black hole binaries.