# Untitled Document

The structure of the central disk of NGC 1068: a clumpy disk model

Pawan Kumar

Institute for Advanced Study, Princeton, NJ 08540

Abstract

NGC 1068 is one of the best studied Seyfert II galaxies, for which the blackhole mass has been determined from the Doppler velocities of water maser. We show that the standard $`\alpha `$-disk model of NGC 1068 gives disk mass between the radii of 0.65 pc and 1.1 pc (the region from which water maser emission is detected) to be about 7x10<sup>7</sup> M (for $`\alpha =0.1`$), more than four times the blackhole mass, and a Toomre Q-parameter for the disk is $``$0.001. This disk is therefore highly self-gravitating and is subject to large-amplitude density fluctuations. We conclude that the standard $`\alpha `$-viscosity description for the structure of the accretion disk is invalid for NGC 1068.

In this paper we develop a new model for the accretion disk. The disk is considered to be composed of gravitationally bound clumps; accretion in this clumped disk model arises because of gravitational interaction of clumps with each other and the dynamical frictional drag exerted on clumps from the stars in the central region of the galaxy. The clumped disk model provides a self-consistent description of the observations of NGC 1068. The computed temperature and density are within the allowed parameter range for water maser emission, and the rotational velocity in the disk falls off as $`r^{0.35}`$.

Subject heading: accretion disks – galaxies: individual (NGC 1068)

Alfred P. Sloan Fellow

1. Introduction

The discovery of emission from water masers in the central region of Seyfert II galaxy NGC 1068 provides an opportunity to study the properties of the associated accretion disk. The mass of the central black hole is estimated to be 1.5x10<sup>7</sup> M (Greenhill & Gwinn, 1997) from the measurement of the Doppler velocities of the masing spots.

The nucleus of NGC 1068 is heavily obscured, and the luminosity of the central source is determined from the observed flux by modeling the dust obscuration and the scattering of photons by ionized gas in the nuclear region. According to a careful analysis carried out by Pier et al. (1994) the bolometric luminosity of NGC 1068 is estimated to be about 8x10<sup>44</sup> erg s<sup>-1</sup>; the luminosity is perhaps uncertain by a factor of a few.

We use the blackhole mass and the bolometric luminosity to construct the standard viscous disk model for NGC 1068 (§2) and find that the disk is highly self-gravitating thereby rendering the $`\alpha `$-disk model inapplicable. The effect of the irradiation of the disk from the central source does not modify this conclusion (§2b).

In section 3 we present a new model for the disk in NGC 1068 composed of gas clumps. The accretion in this case arises as a result of gravitational interaction amongst clumps (§3.1); the structure of the clumped-disk model (CD model) of NGC 1068 is described in §3.2.

The velocity of the masing spots appears to be falling off with distance from the center as r<sup>-0.35</sup> (Greenhill & Gwinn, 1997), which is less rapid than the Keplerian power law of $`0.5`$. One possible reason for this could be that the disk of NGC 1068 is sufficiently massive so as to modify the rotation curve. However, we show in §3 that this is not so. Another possibility is that the flattening of the rotation curve is due to a central star cluster. The influence of the star cluster on the accretion rate and the disk structure is discussed in section 3.2.

2. Standard thin disk model for NGC 1068

The standard viscous accretion disk model for NGC 1068, including the irradiation of the disk from the central source, is presented below. Throughout this paper we take the mass of the blackhole at the center of NGC 1068 to be 1.5x10<sup>7</sup> M (Greenhill & Gwinn, 1997), and the bolometric luminosity to be 8x10<sup>44</sup> erg s<sup>-1</sup> (Pier et al., 1994).

The theory of thin accretion disk is well developed and is described in a number of review articles and monographs e.g. Frank, King & Raine (1992). When the flux intercepted by disk is not small compared to the local energy generation rate, such as that expected for the masing disk of NGC 1068, the incident flux must be included in determining the thermal structure of the disk.

The fractional luminosity intercepted by the disk depends on the disk geometry, the scattering of radiation by the coronal gas etc. For instance if the dominant source of radiation intercepted by the disk were the scattered radiation from an extended corona then we might expect the incident flux at the disk to be roughly uniform. On the other hand for a flaring or a warped disk the radiation intercepted from the source directly might dominate, and the incident radiation in this case depends on the inclination angle of the normal to the disk. We adopt the second model to analyze the disk structure; it is straightforward to consider other possibilities, but we do not pursue these since the main result of this section viz., the standard $`\alpha `$-disk model for NGC 1068 is highly unstable to gravitational perturbations turns out to be independent of whether we include irradiation and radiative forces in the calculation of disk structure.

Let the luminosity of the central source be $`L_c`$. The flux incident at the disk, $`F_{in}`$, a distance $`r`$ from the central source is taken to be

$$F_{in}(r)=\frac{L_c}{4\pi r^2}\frac{dH}{dr}=\frac{\beta HL_c}{4\pi r^3},$$
$`(1)`$   
where $`H`$ is the vertical scale height, the constant factor $`\beta `$ is defined by $`dH/dr=\beta H/r`$, and determines the fraction of the flux that is intercepted by the disk. The upward energy flux due to mass accretion rate $`\dot{M}`$ is

$$F_{up}(r)=\frac{3}{8\pi }\mathrm{\Omega }^2\dot{M},$$
$`(2)`$   
and the ratio of these fluxes $`F_{in}/F_{up}ϵ\beta H/R_{sc}`$, where $`ϵ0.1`$ is the efficiency of the conversion of the rest mass to energy by the central source, and $`R_{sc}`$ is the Schwarzschild radius of the blackhole, $`\mathrm{\Omega }=(GM_t/r^3)^{1/2}`$ is the angular rotation speed, and $`M_t`$ is the mass contained inside the radius $`r`$. The effective temperature of the disk at $`r`$ is given by

$$\sigma T_{eff}^4(r)=F_{up}+F_{in}F_t.$$
$`(3)`$   
The temperature at the mid-plane of the disk ($`T_c`$) can be calculated by considering the first moment of the radiative transfer equation, and making use of the Eddington approximation, and is given by

$$T_c^4=\frac{3}{4}T_{eff}^4\left(\frac{F_{up}}{F_t}\right)\left[\tau _0+\frac{4}{3}\frac{F_t}{F_{up}}\frac{2}{3}\right],$$
$`(4)`$   
where $`\tau _0=\kappa \mathrm{\Sigma }`$ is the optical depth of the disk. The opacity $`\kappa `$, for masing disks, with $`H_2`$ density and temperature in the range 10<sup>8</sup>–10<sup>10</sup> cm<sup>-3</sup> and 100–1000 K respectively is dominated by metal grains and is given by (Bell & Lin 1994)

$$\kappa 0.1T^{1/2}\mathrm{cm}^2\mathrm{g}^1.$$
$`(5)`$   
The temperature in the disk mid-plane is not sensitive to the height where the energy is deposited so long as the optical depth of the disk to the incident radiation is much greater than one.

The solution to the hydrostatic equilibrium equation in the vertical direction, which includes the radiation pressure at the disk surface $`(2F_{in}+F_{up})/c`$, yields

$$\rho _cc_s^2\mathrm{\Omega }^2\mathrm{\Sigma }H+\frac{2F_t}{c}\left(1\frac{\tau _0F_{up}}{2F_t}\right).$$
$`(6)`$   
where $`\mathrm{\Sigma }`$ is the surface mass density, $`c_s`$ & $`\rho _c`$ are the sound speed and gas density in the mid plane of the disk. Taking the effective viscosity in the disk to be $`\alpha c_sH`$ the mass accretion rate $`\dot{M}`$ is given by

$$\dot{M}=\frac{2\pi \alpha \mathrm{\Sigma }c_s^2}{\mathrm{\Omega }}.$$
$`(7)`$   
Equations (4), (6) and (7) are three equations in three unknowns viz. $`\mathrm{\Sigma }`$, $`H`$ and $`T_c`$, which are solved numerically and the results are shown in fig. 1. Also shown in fig. 1, for comparison, are the solutions for $`\beta =0`$ i.e. when the irradience of the disk is neglected. The solution in the latter case can be obtained analytically and is given below

$$T_c=200\alpha ^{2/9}\dot{m}^{4/9}M_7^{7/9}r^1\mathrm{K},$$
$`(8)`$   
$$n=10^{10}\alpha ^{2/3}\dot{m}^{1/3}M_7^{5/6}r^{3/2}\mathrm{cm}^3,$$
$`(9)`$   
and the Toomre-Q parameter for the stability of the disk is

$$Q=\frac{\mathrm{\Omega }c_s}{\pi G\mathrm{\Sigma }}=2.4\times 10^3\alpha ^{2/3}\dot{m}^{1/3}M_7^{1/6}r^{3/2},$$
$`(10)`$   
where $`\dot{m}`$ is the accretion rate in terms of the Eddington rate, $`M_7=M/(10^7M_{})`$, and $`r`$ is the radial distance from the center in parsecs. For the NGC 1068 system $`M_71.5`$ and $`\dot{m}0.4`$. The inner and the outer radii of the masing disk are at 0.65pc and 1.1pc respectively.

For the numerical results shown in fig. 1 we took the value of $`\beta `$ such that the disk intercepts about 50% of the flux from the central source. In this case the disk temperature is dominated by the incident flux, and the temperature in the disk midplane is close to $`T_{eff}511`$K at $`r=1`$pc. The structure of the disk ($`\mathrm{\Sigma }`$, $`Q`$, $`H`$, $`T_c`$) in this case is almost independent of the opacity of the gas and consequently it is unaffected by any uncertainty in $`\kappa `$. The molecular mass of the masing disk, when irradiation is included, is $`7.0`$x10$`{}_{}{}^{7}M_{}^{}`$ (for $`\alpha =0.1`$), and the $`Q`$ is about 10<sup>-3</sup> (see fig. 1). For the irradiation dominated disk, the disk mass decreases with $`\alpha `$ as $`\alpha ^{0.9}`$, and the $`Q`$ increases as $`\alpha ^{0.9}`$. Since $`Q`$ should be greater than one for stability, we see that the masing disk of NGC 1068 is highly unstable to gravitational perturbation. A decrease in the irradiation flux makes the disk more unstable.

The main conclusion of this section is that the standard $`\alpha `$-disk model for the NGC 1068 masing disk is inconsistent. The disk according to this model is highly unstable, and should fragment into clumps. A consistent analysis in this case should include the effect of self-gravity, and clump interaction, which is described below.

3. A model for clumpy self gravitating disks

A disk with small value for ($`1Q`$) is likely to develop spiral structure which can transport angular momentum outward. However, when $`Q1`$, as in the case of NGC 1068 (see §2), the disk is likely to breakup into clumps, and the accretion rate is determined by the gravitational interaction among these clumps. Accretion in self gravitating disks was considered by Paczynski (1978), Lin & Pringle (1987), Shlosman & Begelman (1987), Shlosman et al. (1990), and has been investigated more recently by Kumar (1998) in some detail for clumpy disks. In §3.1 we derive the results that we need to construct model for clumpy disks (CDs), and its application to NGC 1068 is discussed in §3.2.

§3.1 Velocity dispersion and accretion rate in a clumpy disk

Consider a disk consisting of clumps of size $`l_c`$, and take the mean separation between clumps to be $`d_c`$. The tidal radius of a cloud, the distance to which the gravity of the cloud dominates over the gravity of the central mass, is $`d_tr(m_c/M_t)^{1/3}`$; where $`m_c\sigma _cl_c^2`$ is the cloud mass, $`\sigma _c`$ is the surface mass density of the cloud and $`M_t`$ is the total axisymmetrically distributed mass contained inside the radius $`r`$.

The change to the displacement amplitude of the epicyclic oscillation of a clump as a result of gravitational interaction with another clump, with impact parameter $`d\mathrm{}>d_t`$, can be shown to be $`m_cr^3/(M_td^2)`$. Gravitational encounters with $`d\mathrm{}<d_t`$ are almost adiabatic and these leave the epicyclic energy of clumps unchanged (two particles moving on circular orbits merely interchange their trajectories in such an encounter). Thus the dominant gravitational interactions for exciting epicyclic oscillation are those with impact parameter of $`d_t`$, and the change to the epicyclic amplitude in such an encounter is $`\delta rm_cr^3/(M_td_t^2)d_t`$, or

$$\delta rd_c^{2/3}r^{1/3}\left(\frac{M_r}{\pi M_t}\right)^{1/3},$$
$`(11)`$   
where

$$M_r=\frac{\pi r^2\sigma _cl_c^2}{d_c^2}\pi r^2\overline{\sigma }(r),$$
$`(12)`$   
and $`\overline{\sigma }(r)`$ is the mean surface mass density of the disk at radius $`r`$.

The mean time interval for a clump to undergoes strong gravitational interaction with another clump i.e. with impact parameter $`d_t`$ is $`t_{int}\mathrm{\Omega }^1(d_c/d_t)^2`$, and so long as the cloud size ($`l_c`$) is not much smaller than $`d_c`$ the timescale for physical collision between clumps is also of the same order i.e. clouds undergo physical collision after having undergone one strong gravitational encounter on the average, and so the cloud velocity dispersion is $`\mathrm{\Omega }\delta r`$ or

$$\overline{v}_r(r\mathrm{\Omega })\left[\frac{M_r}{\pi M_t}\right]^{1/3}\left[\frac{d_c}{r}\right]^{2/3},$$
$`(13)`$   
The expression for $`\overline{v}_r`$ is similar to that given by Gammie et al. (1991) for the velocity dispersion of molecular clouds in the Galaxy.

We assume that the kinetic energy of the relative motion of colliding clouds is dissipated and their orbit is circularized. Colliding clumps might coalesce so long as the size of the resulting object does not exceed the maximum length for gravitational instability, $`l_{max}(rM_r)/M_t`$; clouds of larger size are susceptible to fragmentation due to the rotational shear of the disk which limits their size.

The characteristic time for a clump to fall to the center can be estimated using equation (11) and is given by

$$t_r\mathrm{\Omega }^1\left(\frac{d_c}{r}\right)^2\left(\frac{M_t}{m_c}\right)^{4/3},$$
$`(14)`$   
and the associated average mass accretion rate is

$$\dot{M}\mathrm{\Omega }(r)M_r\left(\frac{M_r}{\pi M_t}\right)^{4/3}\left(\frac{d_c}{r}\right)^{2/3}.$$
$`(15)`$   
This result is applicable so long as the density contrast between the clump and the inter-clump medium is about a factor of two or larger, and $`l_c`$ is of order $`d_c`$.

For $`l_cd_c`$ clouds undergo several gravitational encounters with other clouds before undergoing a physical collision, and during this time the random velocity of clouds continues to increase. The rate of increase of velocity dispersion in two body gravitational interaction is given by

$$\frac{dv^2}{dt}\frac{G^2m_cM_r\mathrm{\Omega }}{r^2v^2}.$$
$`(16)`$   
We assume that the scattering gives rise to isotropic velocity dispersion, and so the scale height for the vertical distribution of clumps is $`v/\mathrm{\Omega }`$. The time for a clump to undergo physical collision, assumed to be completely inelastic, is $`t_{col}\mathrm{\Omega }^1(d_c/l_c)^2`$, and therefore the velocity dispersion is

$$\overline{v}^2\frac{Gm_c}{l_c}.$$
$`(17)`$   
This corresponds to the Safronov number being equal to one. Since the cloud collision time is much greater than $`\mathrm{\Omega }^1`$ the effective viscosity is suppressed by a factor $`(\mathrm{\Omega }t_{col})^2`$ compared to the case where collision frequency is greater than $`\mathrm{\Omega }`$, (cf. Goldreich and Tremaine, 1978), i.e. $`\nu _e\overline{v}_r^2(t_{col}^1/\mathrm{\Omega }^2)`$, and thus the mass accretion rate is given by

$$\dot{M}M_r\mathrm{\Omega }\left(\frac{M_r}{M_t}\right)\left(\frac{l_c}{r}\right)\mathrm{\Omega }M_r\left(\frac{M_r}{M_t}\right)^{4/3}\left(\frac{d_c}{r}\right)^{2/3}\left(\frac{l_c}{d_t}\right).$$
$`(18)`$   
This accretion rate is smaller than given by equation (15) by a factor of $`(d_t/l_c)`$.

We can parameterize the effect of unknown size distribution and separation between clumps on the accretion rate and the velocity dispersion by a dimensionless parameter $`\eta `$, and rewrite equations for $`\dot{M}`$ and $`\overline{v}_r`$ in the following form for future use:

$$\dot{M}=\eta \mathrm{\Omega }(r)M_r(r)\left(\frac{M_r}{\pi M_t}\right)^2,$$
$`(19)`$   
and

$$\overline{v}_r=\eta r\mathrm{\Omega }(r)\left(\frac{M_r}{\pi M_t}\right).$$
$`(20)`$   
For $`l_cd_cl_{max}`$, we see from equations (16) and (18) that $`\eta 1`$. The effective kinematic viscosity in this case is approximately $`ł_{max}^2\mathrm{\Omega }Q^2H^2\mathrm{\Omega }`$ (where $`Ql_{max}/H`$ is the Toomre $`Q`$-parameter, $`H\overline{v}_r/\mathrm{\Omega }`$), same as in the ansatz suggested by Lin & Pringle (1987) for self-gravitating-disk. We note that the relative velocity of collision between clouds is smaller than their orbital speed by a factor of the ratio of the total mass to the molecular mass, and as long as the cooling time for post-collision gas is less than the time between collisions the clouds are not smeared away due to heating and the disk remains clumpy.

§3.2 Application to NGC 1068

The equation (19) can be recast in the following form

$$\frac{M_r^2}{M_r+M_c}=\psi ,$$
$`(21)`$   
where

$$\psi \frac{\pi ^2}{\eta ^{2/3}}\frac{r\dot{M}^{2/3}}{G^{1/3}}.$$
$`(22)`$   
The solution to this equation in the case where the central mass dominates is

$$M_r\psi ^{1/2}\left(M_c^{1/2}+\frac{\psi ^{1/2}}{2}\right),$$
$`(23)`$   
and from equation (20) we find the velocity dispersion of clouds to be

$$\overline{v}_r\eta ^{2/3}[G\dot{M}]^{1/3}.$$
$`(24)`$   
For NGC 1068, the luminosity $``$ 8x10<sup>44</sup> erg s<sup>-1</sup>, $`M_c`$ 1.5x10<sup>7</sup> M, and $`\dot{M}8`$x10<sup>24</sup> g s<sup>-1</sup>. Substituting these numbers in the above equations we find

$$M_r1.3\times 10^6M_{}r^{1/2}\eta ^{1/3},$$
$`(25)`$   
and

$$\overline{v}_r8\mathrm{km}\mathrm{s}^1\eta ^{2/3},$$
$`(26)`$   
where $`r`$ is measured in parsecs. Note that the disk mass in this model is more than an order of magnitude smaller than in the standard $`\alpha `$-disk model discussed in §2. The velocity dispersion of clumps is larger than the sound speed, and is of order the observed scatter in the velocity of masing spots in NGC 1068.

The number density of molecules and the gas temperature can be obtained by solving the hydrostatic and the thermal equilibrium equations for blobs in the vertical direction in the presence of incident flux (see §2). We find the thermal temperature of blobs to be about 510 K at $`r=1`$pc (see fig. 1), and the density scale height $`H1.5\times 10^{16}`$ cm. The mean number density of $`H_2`$ molecules $`nM_r/(4\pi r^2H\mu )=5\times 10^8`$ cm<sup>-3</sup>. Thus both the number density and the temperature are within the allowed range for water maser emission.<sup>1</sup> The allowed range for number density for water maser emission is 10<sup>8</sup>—10<sup>11</sup> cm<sup>-3</sup> and the temperature range is 200—1000 K. The size of clumps is $`\mathrm{}<r(M_r/M_t)0.1`$ pc and their mass is $`10^3`$ M. Since the Jeans mass is $`10`$ M the clumps could have some star formation activity. However, the efficiency of star formation is usually quite low, of order a few percent even for the giant molecular clouds, so we don’t expect the gas in the clumps to be converted into stars in the short lifetime of order 10<sup>6</sup> years or less for the clumps in the masing disk of NGC 1068.

Note that the velocity dispersion of blobs is independent of $`r`$ (eq. ), so the scale height for the vertical distribution of blobs increases as $`r^{3/2}`$ provided that the velocity distribution of clouds is nearly isotropic; this increase of disk thickness with $`r`$ is more rapid than in the standard $`\alpha `$-disk model considered in the last section. At a distance of 2 pc from the center of NGC 1068 the scale height is about 0.3 pc (for $`\eta =5`$), so the flaring disk blocks a significant fraction of the radiation from the central source.

Clouds colliding at a relative speed of 8 km s<sup>-1</sup> raise the gas temperature to approximately 10<sup>4</sup> K which is too hot for maser emission. However, the cooling time of gas at the density of $`10^9`$ cm<sup>-3</sup>, is of order 100 years, which is short compared to both the orbital time and the collision time of order $`10^4`$ years<sup>2</sup> The density enhancement of the shocked gas depends on the strength of magnetic field in the clumps and is order unity for equipartition magnetic field. ensuring that the disk remains cold and that a steady state solution exists.

The slope of the rotation curve in the disk follows from the use of equation (25)

$$\frac{d\mathrm{ln}V_{rot}}{d\mathrm{ln}r}=0.5+\frac{1}{2}\frac{d\mathrm{ln}M_t}{d\mathrm{ln}r}=0.5+\frac{\psi ^{1/2}}{4M_c^{1/2}}0.5+\frac{M_r}{4M_c}.$$
$`(27)`$   
The disk mass for the NGC 1068 system, in our model, is small ($`M_r/M_t0.1`$), and so the slope of the rotation curve is very close to $`0.5`$.

A slower fall off of the rotation curve requires a more rapid increase of the total mass ($`M_t`$) with $`r`$ than in the model discussed above. This could arise for instance if there is a star cluster at the center. The rotation curve falls off as $`r^{0.35}`$ when the mass of the star cluster within a parsec of the center of NGC 1068 is about 8x10<sup>6</sup> M. Thatte et al. (1997) have carried out near infrared speckle imaging of the central 1” of NGC 1068, and conclude that about 6% of the near-infrared light from this region is contributed by a star cluster. They estimate that the stellar mass within 1” ($`50`$pc) of the nucleus of NGC 1068 is about 6x10<sup>8</sup> M; if the density in this cluster falls off as $`r^2`$, as in a singular isothermal sphere, then the expected stellar mass inside 1 pc is 1.2x10<sup>7</sup> M. We note that the stellar mass within 100 pc of our own galactic center is estimated to be in excess of 5x10<sup>8</sup> M (Kormendy & Richstone 1995, Genzel et al. 1994), and the mass enclosed within radius $`r`$ is seen to increase almost linearly with $`r`$ for $`r\mathrm{}>2`$pc. Thus the possibility of a similar stellar mass cluster at the center of NGC 1068 is not surprising.

We discuss below the effect a star cluster has on the structure of the accretion disk.

§3.2.1 Effect of a star cluster on accretion disk

Let us consider that the stellar mass within a radius $`r`$ of the center is $`M_s(r)`$. The disk mass, as before, is taken to be $`M_r(r)`$, and the central mass is $`M_c`$. The use of equation (23), which still applies with $`M_c`$ replaced by $`(M_c+M_s)`$, implies that the disk structure is not much affected at small radii where $`M_s(r)M_c`$. At larger radii, where the stellar mass becomes comparable to or exceeds $`M_c`$, the mass of the gas disk ($`M_r`$) increases with $`r`$ as $`r^{1/2}(M_c+M_s)^{1/2}`$, which is somewhat more rapid than the case of $`M_s=0`$ considered in §3.1. However, for the masing disk of NGC 1068 $`M_s\mathrm{}<M_c`$ for $`r\mathrm{}<1.6`$ pc, and the effect of the stellar cluster on the disk structure, in particular the number density of molecules and the gas temperature, is small.

The velocity dispersion of clouds due to gravitational encounters is proportional to the local surface mass density $`\overline{\sigma }`$ of gas and is unaffected by the stellar cluster (see eq. ). Therefore, the disk thickness (H) at first increases with radius as $`r^{3/2}`$ and then flattens out when $`M_s`$ starts to dominate the total mass:

$$H\frac{\overline{v}_r}{\mathrm{\Omega }(r)}\frac{\eta r^3\overline{\sigma }(r)}{M_t}.$$
$`(28)`$   
For $`\overline{\sigma }r^{3/2}`$, expected of the CD model with constant $`\dot{M}`$, $`H/rr^{1/2}`$, however a more rapid decrease of $`\overline{\sigma }`$ with $`r`$ leads to a corresponding decrease of $`H/r`$.<sup>3</sup> A decrease of $`\overline{\sigma }`$ with $`r`$ that is more rapid than $`r^{3/2}`$ leads to a drop in $`\dot{M}`$. However, the assumption of steady state accretion breaks down beyond some radius where the accretion time, $`\mathrm{\Omega }^1(M_c+M_s)^2/M_r^2`$, is longer compared with the evolution time scale of the disk.

The accretion rate calculated in §3.1 is modified due to the dynamical friction suffered by clumps in the disk by the star cluster; the orbits of stars in the cluster are gravitationally perturbed by the clumps so that on the average the density of stars behind a clump is greater than the density in the front. We estimate below the accretion rate that arises as a result of the frictional drag exerted by stars. We assume that clumps are not stretched out by the tidal force of stars, which is valid so long as the cloud mass density is greater than the mean mass density associated with the star cluster i.e. $`n_{H_2}\mathrm{}>M_s/(m_{H_2}r^3)`$ or $`n_{H_2}\mathrm{}>2\times 10^8\mathrm{cm}^3(M_s/10^7M_{})(r/1pc)^3`$; this condition is satisfied for NGC 1068.

The dynamical friction timescale for a clump to fall to the center is (cf. Binney and Tremaine, 1987)

$$t_{df}\frac{r^2M_tv_{rot}}{Gm_cM_s\mathrm{ln}\mathrm{\Lambda }},$$
$`(29)`$   
where $`m_c`$ is clump mass, $`\mathrm{\Lambda }rv_{rot}^2/(Gm_c)`$, and $`v_{rot}^2=GM_t/r`$. Taking the clump size to be the largest length scale for gravitational instability in a shearing disk i.e. $`l_c\pi ^2G\overline{\sigma }/\mathrm{\Omega }^2`$, we find $`m_c\pi M_r^3/M_t^2`$. Thus, $`\mathrm{ln}\mathrm{\Lambda }3\mathrm{ln}(M_t/M_r)`$ is about 7 for the NGC 1068 system. Substituting these in the above equation we find

$$t_{df}\frac{\mathrm{\Omega }^1}{10}\left(\frac{M_t^4}{M_sM_r^3}\right).$$
$`(30)`$   
For $`M_sM_t/3`$ and $`M_r/M_t0.1`$, values applicable to the NGC 1068 system, the dynamical friction time is of the same order as the timescale for clumps to fall to the center due to gravitational encounter with other clumps (see eq. ). Thus we see that a star cluster of modest mass at the center of NGC 1068 can both give rise to a sub-Keplerian rotation curve, as perhaps observed by the water maser, and also because of its dynamical friction on clumps remove angular momentum of the gaseous disk resulting in accretion rate that is of the same order as needed to account for the observed luminosity (the total $`\dot{M}`$ is the sum of the accretion rate due to gravitational interaction between clumps, and the dynamical friction drag exerted by stars).

It was pointed out by Ostriker (1983) that even a smooth disk is subject to frictional drag from a star cluster. He showed that the characteristic drag time for removing angular momentum of a gaseous disk is of the same order as the relaxation time for the star cluster. The accretion rate due to this process, in a form applicable to the NGC 1068 system, is given by

$$\dot{M}M_r\mathrm{\Omega }I_0\left(\frac{r_{}}{r}\right)^2\left(\frac{M_s}{m_{}}\right),$$
$`(31)`$   
where $`r_{}`$ is the radius and $`m_{}`$ is the mass of stars in the cluster, and $`I_010`$ is a dimensionless quantity that depends on the ratio of the escape velocity at the surface of a star to the stellar velocity dispersion. We see that the accretion rate resulting from friction drag from stars on a smooth disk is much smaller than the rate resulting from the frictional drag on clumps calculated earlier.

4. Conclusion

We find that the structure of masing disk of NGC 1068, as determined using the standard $`\alpha `$-viscosity prescription, is both inconsistent with the observations and is self contradictory. In particular, the disk mass according to this model is 7x10<sup>7</sup> M, which is much greater than the blackhole mass, and the Toomre $`Q`$ parameter has an extremely low value of $`10^3`$ (for $`\alpha =1`$). This makes the disk highly unstable to gravitational perturbations, and suggests that it breaks up into clumps, a conclusion that is inconsistent with the smooth disk assumption of standard $`\alpha `$-disks.

We have described in this paper a different model for the disk. We consider the disk to consist of gas clumps which undergo gravitational interactions with one another leading to an inward accretion of gas. The mass of the masing disk of NGC 1068 according to this model is about 10<sup>6</sup> M, and the velocity dispersion of clumps is about 10 km s<sup>-1</sup> which is in rough agreement with the observed velocity dispersion of masing spots. The temperature and the density of clumps in this model are approximately 510 K and 5x10<sup>8</sup> cm<sup>-3</sup> respectively, which are within the allowed parameter range for maser emission.

However, the rotational velocity of clumps, in this model, are close to the Keplerian value, whereas the observations of masing spots indicate a slower fall off of their rotational velocity (Greenhill & Gwinn, 1997). One obvious way these results can be reconciled is if the mass in stars within 1 pc of the center of NGC 1068 is of order the blackhole mass. Speckle observations of the central 1” region of NGC 1068 in near infrared in fact suggests that the stellar mass contained within 50 pc of the center is about 5x10<sup>8</sup> M (Thatte et al. 1997), which is sufficient to explain the flattening of the observed rotation curve in the masing disk. The dynamical friction exerted by this star cluster on the clumps in the disk removes angular momentum at a rate that is of the same order as needed for the nearly Eddington accretion rate for the system. Thus the clump-clump gravitational interaction and the dynamical friction drag force on clumps together determine the disk structure, which we find to be consistent with observations.

The clumpy-disk (CD) model considered in this paper for the disk of NGC 1068 should also apply to any AGN at a distance of about a pc or more from the center, depending on the accretion rate and the luminosity, where the disk becomes self-gravitating and the standard $`\alpha `$ model is too inefficient to account for the accretion rate (Shlosman & Begelman, 1989).

Acknowledgment: I am indebted to Ramesh Narayan for sharing his work and insights on the accretion disk of NGC 1068, and for his detailed comments on this work. I am grateful to Scott Tremaine for very helpful discussion about disk dynamics. I thank John Bahcall, Mitch Begelman, Kathryn Johnston, Douglas Richstone, Phil Maloney and an anonymous referee for helpful comments.

REFERENCES

Bell, K.R. & Lin, D.N.C., 1994, ApJ 427, 987

Binney, J.J. & Tremaine, S.D., 1987, Galactic Dynamics, Princeton University Press, Princeton

Frank, J., King, I.R., & Raine, D., 1992, Accretion Power in Astrophysics, (Cambridge: Cambridge univ. press)

Gammie, C., Ostriker, J.P., and Jog, C., 1991, ApJ 378, 565

Genzel, R., Hollenbach, D. and Townes, C.H., 1994, Rep. Prog. Phys. 57, 417

Goldreich, P., and Tremaine, S., 1978, ICARUS 34, 227

Greenhill, L.J., & Gwinn, C.R., 1997, CfA preprint no. 4508

Kormendy, J. & Richstone, D., 1995, ARA & A 33, 581

Lin, D.N.C., and Pringle, J.E., 1987, MNRAS 225, 607

Ostriker, J.P., 1983, ApJ 273, 99

Paczynski, B., 1978, Acta Astr. 28, 91

Pier, E.A., Antonucci, R., Hurt, T., Kriss, G., and Krolik, J., 1994, ApJ 428, 124

Shlosman, I., and Begelman, M., 1987, Nature 329, 810

Shlosman, I., and Begelman, M., 1989, ApJ 341, 685

Shlosman, I., Begelman, M., and Frank, J., 1990, Nature 345, 679

Thatte, N., Quirrenbach, A., Genzel, R., Maiolino, R. & Tecza, M., 1997, ApJ 490, 238

Figure Caption

FIG. 1.— The panel on the left shows the gas temperature ($`T_c`$) and the number density of $`H_2`$ molecules ($`n_c`$) in the mid plane of the disk for the standard $`\alpha `$-disk model without any irradiation from the central source (thin continuous line, and thin dashed line respectively). The thick continuous and dashed lines are $`T_c`$ and $`n_c`$ respectively for an $`\alpha `$-disk model that includes irradiation from the central source. The irradiation is taken to be proportional to $`H/r`$, where $`H`$ is the scale height of the disk, and the proportionality constant ($`\beta `$) is chosen so that the flux intercepted by the disk at distance $`r`$ is about 50% of the flux from the central source at $`r`$ (such a large value for the flux intercepted by the disk corresponds to $`\beta =50`$ which might arise if the disk were extremely warped). The panel on the right shows the Toomre Q-parameter for the two disk models; thin continuous curve is for the standard $`\alpha `$-disk with no irradiation, and the thick curve is for a disk that intercepts about 50% of the flux. Note that $`Q`$ is much less than one in both of these models which corresponds to the disk being highly unstable to gravitational perturbations. In both of these models we chose $`\alpha =0.1`$ ($`Q`$ scales as roughly $`\alpha ^{0.9}`$ when irradiation dominates), the blackhole mass $`M=1.5\times 10^7M_{}`$, and the mass accretion rate ($`\dot{m}`$) in units of the Eddington rate was taken to be 0.4.