Although ultra-luminous X-ray sources (ULX) are important for astrophysics due to their extreme apparent super-Eddington luminosities, their nature is still poorly known. Theoretical and observational studies suggest that ULXs could be a diversified
group of objects composed of low-mass X-ray binaries, high-mass X-ray binaries and marginally also systems containing intermediate-mass black holes, which is supported by their presence in a variety of environments. Observational data on the ULX donors could significantly boost our understanding of these systems, but only a few were detected. There are several candidates, mostly red supergiants (RSGs), but surveys are typically biased toward luminous near-infrared objects. Nevertheless, it is worth exploring if RSGs can be members of ULX binaries. In such systems matter accreted onto the compact body would have to be provided by the stellar wind of the companion, since a Roche-lobe overflow could be unstable for relevant mass-ratios. Here we present a comprehensive study of the evolution and population of wind-fed ULXs and provide a theoretical support for the link between RSGs and ULXs. Our estimated upper limit on contribution of wind-fed ULX to the overall ULX population is $sim75$--$96%$ for young ($<100$ Myr) star forming environments, $sim 49$--$87%$ for prolonged constant star formation (e.g., disk of Milky Way), and $lesssim1%$ for environments in which star formation ceased long time ($>2$ Gyr) ago. We show also that some wind-fed ULXs (up to $6%$) may evolve into merging double compact objects (DCOs), but typical systems are not viable progenitors of such binaries because of their large separations. We demonstrate that, the exclusion of wind-fed ULXs from population studies of ULXs, might have lead to systematical errors in their conclusions.
It is now widely accepted that most ultraluminous X-ray sources (ULXs) are binary systems whose large (above $10^{39}$ erg s$^{-1}$) apparent luminosities are explained by super-Eddington accretion onto a stellar-mass compact object. Many of the ULXs
, especially those containing magnetized neutron stars, are highly variable; some exhibit transient behaviour. Large luminosities might imply large accretion discs that could be therefore prone to the thermal-viscous instability known to drive outbursts of dwarf novae and low-mass X-ray binary transient sources. The aim of this paper is to extend and generalize the X-ray transient disc-instability model to the case of large (outer radius larger than $10^{12}$ cm) accretion discs and apply it to the description of systems with super-Eddington accretion rates at outburst and, in some cases, super-Eddington mass transfer rates. We have used our disc-instability-model code to calculate the time evolution of the accretion disc and the outburst properties. We show that, provided that self-irradiation of the accretion disc is efficient even when the accretion rate exceeds the Eddington value, possibly due to scattering back of the X-ray flux emitted by the central parts of the disc on the outer portions of the disc, heating fronts can reach the discs outer edge generating high accretion rates. We also provide analytical approximations for the observable properties of the outbursts. We have reproduced successfully the observed properties of galactic transients with large discs, such as V404 Cyg, as well as some ULXs such as M51 XT-1. Our model can reproduce the peak luminosity and decay time of ESO 243-39 HLX-1 outbursts if the accretor is a neutron star. Observational tests of our predicted relations between the outburst duration and decay time with peak luminosity would be most welcome.
We show that ultraluminous X-ray sources (ULXs) with coherent X-ray pulsing (PULXs) probably have neutron-star spin axes significantly misaligned from their central accretion discs. Scattering in the funnels collimating their emission and producing t
heir apparent super-Eddington luminosities is the most likely origin of the observed correlation between pulse fraction and X-ray photon energy. Pulsing is suppressed in systems with the neutron-star spin closely aligned to the inner disc, explaining why some ULXs show cyclotron features indicating strong magnetic fields, but do not pulse. We suggest that alignment (or conceivably, field suppression through accretion) generally occurs within a fairly short fraction of the ULX lifetime, so that most neutron-star ULXs become unpulsed. As a result we further suggest that almost all ULXs actually have neutron-star accretors, rather than black holes or white dwarfs, reflecting their progenitor high-mass X-ray binary and supersoft X-ray source populations.
We consider the current observed ensemble of pulsing ultraluminous X-ray sources (PULXs). We show that all of their observed properties (luminosity, spin period, and spinup rate) are consistent with emission from magnetic neutron stars with fields in
the usual range $10^{11} - 10^{13}, {rm G}$, which is collimated (`beamed) by the outflow from an accretion disc supplied with mass at a super-Eddington rate, but ejecting the excess, in the way familiar for other (non-pulsing) ULXs. The observed properties are inconsistent with magnetar-strength fields in all cases. We point out that all proposed pictures of magnetar formation suggest that they are unlikely to be members of binary systems, in agreement with the observation that all confirmed magnetars are single. The presence of magnetars in ULXs is therefore improbable, in line with our conclusions above.
We have analyzed how anisotropic emission of radiation affects the observed sample of ultraluminous X-ray sources (ULXs) by performing simulations of the evolution of stellar populations, employing recent developments in stellar and binary physics, a
nd by utilizing a geometrical beaming model motivated by theory and observation. Whilst ULXs harboring black hole accretors (BH ULXs) are typically emitting isotropically, the majority of ULXs with neutron star accretors (NS ULXs) are found to be beamed. These findings confirm previous assertions that a significant fraction of ULXs are hidden from view due to a substantial misalignment of the emission beam and the line-of-sight. We find the total number of NS ULXs in regions with constant star formation, solar metallicity, and ages above ~1 Gyr to be higher than the BH ULXs, although observationally both populations are comparable. For lower metallicities BH ULX dominate both the total and observed ULX populations. As far as burst star-formation is concerned, young ULX populations are dominated by BH ULXs, but this changes as the population ages and, post star-formation, NS ULXs dominate both the observed and total population of ULXs. We also compare our simulation output to a previous analytical prediction for the relative ratio of BH to NS ULXs in idealized flux-limited observations and find broad agreement for all but the lowest metallicities. In so doing we find that in such surveys the observed ULX population should be heavily dominated by black-hole systems rather than by systems containing neutron stars.
The disk instability model attributes the outbursts of dwarf novae to a thermal-viscous instability of their accretion disk, an instability to which nova-like stars are not subject. We aim to test the fundamental prediction of the disk instability mo
del: the separation of cataclysmic variables (CVs) into nova-likes and dwarf novae depending on orbital period and mass transfer rate from the companion. We analyse the lightcurves from a sample of ~130 CVs with a parallax distance in the Gaia DR2 catalogue to derive their average mass transfer rate. The method for converting optical magnitude to mass accretion rate is validated against theoretical lightcurves of dwarf novae. Dwarf novae (resp. nova-likes) are consistently placed in the unstable (resp. stable) region of the orbital period - mass transfer rate plane predicted by the disk instability model. None of the analyzed systems present a challenge to the model. These results are robust against the possible sources of error and bias that we investigated. Lightcurves from Kepler or, in the future, the LSST or Plato surveys, could alleviate a major source of uncertainty, the irregular sampling rate of the lightcurves, assuming good constraints can be set on the orbital parameters of the CVs that they happen to target. The disk instability model remains the solid base on which to construct the understanding of accretion processes in cataclysmic variables.
FO Aquarii, an asynchronous magnetic cataclysmic variable (intermediate polar) went into a low-state in 2016, from which it slowly and steadily recovered without showing dwarf nova outbursts. This requires explanation since in a low-state, the mass-t
ransfer rate is in principle too low for the disc to be fully ionized and the disc should be subject to the standard thermal and viscous instability observed in dwarf novae. We investigate the conditions under which an accretion disc in an intermediate polar could exhibit a luminosity drop of 2 magnitudes in the optical band without showing outbursts. We use our numerical code for the time evolution of accretion discs, including other light sources from the system (primary, secondary, hot spot). We show that although it is marginally possible for the accretion disc in the low-state to stay on the hot stable branch, the required mass-transfer rate in the normal state would then have to be extremely high, of the order of 10$^{19}$ gs$^{-1}$ or even larger. This would make the system so intrinsically bright that its distance should be much larger than allowed by all estimates. We show that observations of FO Aqr are well accounted for by the same mechanism that we have suggested as explaining the absence of outbursts during low states of VY Scl stars: during the decay, the magnetospheric radius exceeds the circularization radius, so that the disc disappears before it enters the instability strip for dwarf nova outbursts. Our results are unaffected, and even reinforced, if accretion proceeds both via the accretion disc and directly via the stream during some intermediate stages; the detailed process through which the disc disappears still needs investigations.
Recent observations of two unusual Z Cam systems, V513 Cas and IW And have shown light curves that seem to contradict the disc-instability model for dwarf novae: outbursts are appearing during standstills of the system when according to the model, th
e disc is supposed to be in a hot quasi-equilibrium state. We investigate what additional physical processes need to be included in the model to reconcile it with observations of such anomalous Z Cam systems. We used our code for modeling thermal-viscous outbursts of the accretion discs and determined what types of mass-transfer variations reproduce the observed light curves. Outbursts of mass transfer (with a duration of a few days, with a short rise time and an exponential decay) from the stellar companion will account for the observed properties of V513 Cas and IW And, provided they are followed by a short but significant mass-transfer dip. The total mass involved in outbursts is of the order of 10$^{23}$g. We studied the possible origins of these mass transfer outbursts and showed that they most probably result from a giant flare near the secondary star surface, possibly due to the absence of star spots in the $L_1$ region.
The study of outer disc radius variations in close binary systems is important for understanding the structure and evolution of accretion discs. These variations are predicted by models of both quasi steady and time-dependent discs, and these predict
ions can be confronted with observations. We consider theoretical and observational consequences of such variations in cataclysmic variables and low mass X-ray binaries. We find that the action of tidal torques, that determine the outer radius at which the disc is truncated, must be important also well inside the tidal radius. We conclude that it is doubtful that the tidal-thermal instability is responsible for the superoutburst/superhump phenomena in dwarf novae, and confirm that it cannot be the reason for the outbursts of soft X-ray transients. It is likely that tidal torques play a role during superoutbursts of very-low mass-ratio systems but they cannot be the main and only cause of superhumps.
