We build an evolution model of the central black hole that depends on the processes of gas accretion, the capture of stars, mergers as well as electromagnetic torque. In case of gas accretion in the presence of cooling sources, the flow is momentum-driven, after which the black hole reaches a saturated mass; subsequently, it grows only by stellar capture and mergers. We model the evolution of the mass and spin with the initial seed mass and spin in $\Lambda$CDM cosmology. For stellar capture, we have assumed a power-law density profile for the stellar cusp in a framework of relativistic loss cone theory that include the effects of black hole spin, Carter's constant, loss cone angular momentum, and capture radius. Based on this, the predicted capture rates of $10^{-5}$--$10^{-6}$ yr$^{-1}$ are closer to the observed range. We have considered the merger activity to be effective for $z \lesssim 4$, and we self-consistently include the Blandford-Znajek torque. We calculate these effects on the black hole growth individually and in combination, for deriving the evolution. Before saturation, accretion dominates the black hole growth ($\sim 95\%$ of the final mass), and subsequently, stellar capture and mergers take over with roughly equal contribution. The simulations of the evolution of the $M_{\bullet} - \sigma$ relation using these effects are consistent with available observations. We run our model backward in time and retrodict the parameters at formation. Our model will provide useful inputs for building demographics of the black holes and in formation scenarios involving stellar capture.

Black holes at the centers of the galaxies grow mainly by the processes of accretion, mergers, and consumption of stars. In the case of gas accretion with cooling sources, the flow is momentum driven, after which the black hole reaches a saturated mass, and subsequently, it grows only by consumption of stars. In addition, we include the effect of mergers on the growth of black hole spin and mass and study its evolution as a function of redshift in a $\Lambda$CDM cosmology using an initial seed mass and spin distribution functions that we have derived. For the stellar ingestion, we have assumed a power-law density profile for the galaxy in our framework of a new relativistic loss cone theory that includes the effect of the black hole spin. We predict the impact of the evolution on the $M_{\bullet} - \sigma$ relation and compare it with available observations.

To investigate the $M_\bullet -\sigma$ relation, we consider realistic elliptical galaxy profiles that are taken to follow a single power law density profile given by $\rho(r) = \rho_{0}(r/ r_{0})^{-\gamma}$ or the Nuker intensity profile. We calculate the density using Abel's formula in the latter case by employing the derived stellar potential in both cases, we derive the distribution function $f(E)$ of the stars in presence of the supermassive black hole (SMBH) at the center and hence compute the line of sight (LOS) velocity dispersion as a function of radius. For the typical range of values for masses of SMBH, we obtain $M_{\bullet} \propto \sigma^{p}$ for different profiles. An analytical relation $p = (2\gamma + 6)/(2 + \gamma)$ is found which is in reasonable agreement with observations (for $\gamma$ = 0.75 - 1.4, $p$ = 3.6 - 5.3). Assuming that a proportionality relation holds between the black hole mass and bulge mass, $ M_{\bullet} =f M_b$, and applying this to several galaxies we find the individual best fit values of $p$ as a function of $f$; also by minimizing $\chi^{2}$, we find the best fit global $p$ and $f$. For Nuker profiles we find that $p$ = $3.81 \pm 0.004$ and $f$ = $(1.23 \pm 0.09)\times 10^{-3}$ which are consistent with the observed ranges.