Search

Radial structure of accretion discs around compact objects is often described using analytic approximations which are derived from averaging or integrating vertical structure equations. For non-solar chemical composition, partial ionization, or for supermassive black holes, this approach is not accurate. Additionally, radial extension of 'analytically-described' disc zones is not evident in many cases. We calculate vertical structure of accretion discs around compact objects, with and without external irradiation, with radiative and convective energy transport taken into account. For this, we introduce a new open Python code, allowing different equations of state and opacity laws, including tabular values. As a result, radial structure and stability 'S-curves' are calculated for specific disc parameters and chemical composition. In particular, based on more accurate power-law approximations for opacity in the disc, we supply new analytic formulas for the farthest regions of the hot disc around stellar-mass object. On calculating vertical structure of a self-irradiated disc, we calculate a self-consistent value of the irradiation parameter C irr for stationary α-disc. We find that, for a fixed shape of the X-ray spectrum, C irr depends weakly on the accretion rate but changes with radius, and the dependence is driven by the conditions in the photosphere and disc opening angle. The hot zone extent depends on the ratio between irradiating and intrinsic flux: corresponding relation for |$T_{\rm irr,\, crit}$| is obtained. [ABSTRACT FROM AUTHOR]

We present a model of a viscously evolving accretion disc around a magnetized neutron star. The model features the varying outer radius of the hot ionized part of the disc due to cooling and the varying inner radius of the disc due to interaction with the magnetosphere. It also includes hindering of accretion on the neutron star because of the centrifugal barrier and irradiation of the outer disc and companion star by X-rays from the neutron star and disc. When setting inner boundary conditions, we take into account that processes at the inner disc occur on a time-scale much less than the viscous time-scale of the whole disc. We consider three types of outflow from the disc inner edge: zero outflow, one based on MHD calculations, and a very efficient propeller mechanism. The light curves of an X-ray transient after the outburst peak can be calculated by a corresponding, publicly available code. We compare observed light curves of the 2013 burst of Aql X-1 in X-ray and optical bands with modelled ones. We find that the fast drop of the 0.3–10 keV flux can be solely explained by a radial shrinking of the hot disc. At the same time, models with the neutron star magnetic field >108 G have better fits because the accretion efficiency behaviour emphasizes the 'knee' on the light curve. We also find that a plato emission can be produced by a 'disc-reservoir' with stalled accretion. [ABSTRACT FROM AUTHOR]