Pre-calculated turning points of Stability-curve for stationary Shakura-Sunyaev accretion disc

The pre-calculated turning points for surface density \(\Sigma\), effective temperature \(T_{\rm eff}\) and accretion rate \(\dot{M}\) of S-shaped stability curve for stationary Shakura-Sunyaev viscous accretion disc. Calculations performed for 20 linearly scaled values of central source mass from \(1\, M_{\odot}\) to \(20\, M_{\odot}\), 20 logarithmically scaled values of viscous parameter \(\alpha\) from \(3\cdot10^{-4}\) to 0.7, 20 logarithmically scaled values of radius from \(7\cdot10^7\) cm to \(5\cdot10^{11}\) cm.

The new open-source Python3 code is used to calculate the vertical structure of the disc with tabular opacity (solar chemical composition) and both convective and radiative energy transport taking into account. As a result, one can obtain S-shaped stability curve — a graphically depicted sequence of solutions of the vertical-structure equations, obtained at a single disc radius, in the coordinates of accretion rate or effective temperature versus the surface density. S-curve turn points allow to investigate the stability of disc with respect to temperature and surface density fluctuations. 

Reading of the data (Python3):

import numpy as np
data = np.genfromtxt('Turn_points_array.dat', names=True)