By exploring more than sixty thousand quasars from the Sloan Digital Sky Survey Data Release 5, Steinhardt & Elvis discovered a sub-Eddington boundary and a redshift-dependent drop-off at higher black hole mass, possible clues to the growth history of massive black holes. Our contribution to this special issue of Universe amounts to an application of a model for black hole accretion and jet formation to these observations. For illustrative purposes, we include ~100,000 data points from the Sloan Digital Sky Survey Data Release 7 where the sub-Eddington boundary is also visible and propose a theoretical picture that explains these features. By appealing to thin disk theory and both the lower accretion efficiency and the time evolution of jetted quasars compared to non-jetted quasars in our “gap paradigm”, we explain two features of the sub-Eddington boundary. First, we show that a drop-off on the quasar mass-luminosity plane for larger black hole mass occurs at all redshifts. But the fraction of jetted quasars is directly related to the merger function in this paradigm, which means the jetted quasar fraction drops with decrease in redshift, which allows us to explain a second feature of the sub-Eddington boundary, namely a redshift dependence of the slope of the quasar mass–luminosity boundary at high black hole mass stemming from a change in radiative efficiency with time. We are able to reproduce the mass dependence of, as well as the oscillating behavior in, the slope of the sub-Eddington boundary as a function of time. The basic physical idea involves retrograde accretion occurring only for a subset of the more massive black holes, which implies that most spinning black holes in our model are prograde accretors. In short, this paper amounts to a qualitative overview of how a sub-Eddington boundary naturally emerges in the gap paradigm.

By exploring more than sixty thousand quasars from the Sloan Digital Sky Survey Data Release 5, Steinhardt & Elvis discovered a sub-Eddington boundary and a redshift-dependent drop-off at higher black hole mass, possible clues to the growth history of massive black holes. Our contribution to this special issue of Universe amounts to an application of a model for black hole accretion and jet formation to these observations. For illustrative purposes we include about 100,000 data points from the Sloan Digital Sky Survey Data Release 7 where the sub-Eddington boundary is also visible, and propose a theoretical picture that explains these features. By appealing to thin disk theory and both the lower accretion efficiency and the time evolution of jetted quasars compared to non-jetted quasars in our 'gap paradigm', we explain two features of the sub-Eddington boundary. First, we show that a drop-off on the quasar mass-luminosity plane for larger black hole mass occurs at all redshifts. But the fraction of jetted quasars is directly related to the merger function in this paradigm, which means the jetted quasar fraction drops with decrease in redshift, which allows us to explain a second feature of the sub-Eddington boundary, namely a redshift dependence of the slope of the quasar mass-luminosity boundary at high black hole mass stemming from a change in radiative efficiency with time. We are able to reproduce the mass dependence of, as well as the oscillating behavior in, the slope of the sub-Eddington boundary as a function of time. The basic physical idea involves retrograde accretion occurring only for a subset of the more massive black holes which implies that most spinning black holes in our model are prograde accretors. In short, this paper amounts to a qualitative overview of how a sub-Eddington boundary naturally emerges in the gap paradigm.

Approximating differential operators defined on two-dimensional surfaces is an important problem that arises in many areas of science and engineering. Over the past ten years, localized meshfree methods based on generalized moving least squares (GMLS) and radial basis function finite differences (RBF-FD) have been shown to be effective for this task as they can give high orders of accuracy at low computational cost, and they can be applied to surfaces defined only by point clouds. However, there have yet to be any studies that perform a direct comparison of these methods for approximating surface differential operators (SDOs). The first purpose of this work is to fill that gap. For this comparison, we focus on an RBF-FD method based on polyharmonic spline kernels and polynomials (PHS+Poly) since they are most closely related to the GMLS method. Additionally, we use a relatively new technique for approximating SDOs with RBF-FD called the tangent plane method since it is simpler than previous techniques and natural to use with PHS+Poly RBF-FD. The second purpose of this work is to relate the tangent plane formulation of SDOs to the local coordinate formulation used in GMLS and to show that they are equivalent when the tangent space to the surface is known exactly. The final purpose is to use ideas from the GMLS SDO formulation to derive a new RBF-FD method for approximating the tangent space for a point cloud surface when it is unknown. For the numerical comparisons of the methods, we examine their convergence rates for approximating the surface gradient, divergence, and Laplacian as the point clouds are refined for various parameter choices. We also compare their efficiency in terms of accuracy per computational cost, both when including and excluding setup costs.