A&A 440, L13-L16 (2005)
DOI: 10.1051/0004-6361:200500157

## The effects of expansion line opacity in accretion disks

G. Shaviv1,2 - R. Wehrse1

1 - Institut f r Theoretische Astrophysik Tiergartenstrasse 15, 69121 Heidelberg, Germany
2 - On leave from Department of Physics and Asher Space Research Institute, Technion, 32000 Haifa, Israel

Recieved 10 June 2005 / Accepted 19 July 2005

Abstract
We estimate the effect of the shear flow in accretion disks on the radiation field. We apply the recent results of Wehrse et al. (2003, A&A, 401, 43) for the expansion line opacity in shear flows. We use an approximation and derive an approximate equation for the radiation field in an accretion disk in the diffusion approximation. We find that in the appropriate limit for the flow in disks:
- The effect depends on a product of three factors: the first depends solely atomic physics and hence on the composition of the accreted matter, the second depends on the density while the third depends solely on the shear flow.

- The radiation field in accretion disks around neutron stars, black holes and white dwarfs is strongly affected by the expansion line opacity.

- The radiation field is very non isotropic leading to significantly higher specific intensities in the vertical direction, a higher temperature gradient and increased susceptibility to convection and other instabilities.

Key words: disks, shear flow, expansion opacity

### 1 Introduction

Radiative transfer in moving media is affected by velocity gradients. The effect, dubbed "expansion opacity'', was first evaluated by Karp et al. (1977). In the simplest model of an accretion disks the flow is Keplerian and hence possess a shear in the radial direction. Hence, a photon propagating inside the disk sees an anisotrpic shear flow. The enhanced effective line opacity, called expansion opacity, is very non isotropic. The main purpose of this letter is to expose the importance of this effect to the structure of accretion disks (here after AD).

### 2 The expansion opacity

#### 2.1 Effective extinction in a given direction

Let be a unit vector in the direction of the ray and a velocity field. Define the scalar "velocity gradient'' as:

 (1)

In the general case of a moving media the radiative flux vector (to distinguish from the flux in a given direction ) is given by (Wehrse et al., Paper II):

 (2)

where is a generalized extinction tensor namely:

 (3)

is the wavelength and s0 the location. The static extinction coefficient is independent of direction and hence the static flux vector is:

 (4)

#### 2.2 Total effective opacity

Clearly, is a scalar function, the value of which depends on direction. Consequently, the effective opacity is a function of direction. On the other hand, the extinction  becomes a tensor. The deviation of the radiation vector, which is proportional to  , from  is important for induced circulation in the ADs and will be treated elsewhere. The total flux in direction  is obtained by integration over the solid angle and is:

 (5)

where

 (6)

and

 (7)

### 3 The effective opacity

We consider here a general case of a series of infinitely sharp lines with a distribution in amplitude. Let the line density be:

 (8)

where is the number density per unit logarithmic wavelength interval and f(A) is the normalized line strength distribution. In this simplified picture we assume that the line strengths and line density are independent. The mean line strength is given by:

 (9)

The total line extinction is then written as: and a new parameter  is defined as

 (10)

with

 (11)

where and . We note that the  dependence in  cancels out. is the Doppler spread over a mean free path  of a photon in the continuum and  is the line density  times the velocity change over one mean free path of a photon in an average spectral line. Once f(A) is given, the effective opacity can be calculated. The result under the present assumptions can be cast into the following general form:

 (12)

where the function depends on the assumption about the line distribution.

The function was calculated by Wehrse et al. (2003, Paper III) for several assumptions about the distribution of lines. For example, in the case of Poisson distributed infinitely sharp lines it has the form:

 (13)

Under these conditions, the effect is negligible for and reaches saturation for . The second parameter  may be as high as few  which means that the saturation level of the effect due to velocity gradient can be as high as 105 times the continuum.

For we find (Wehrse et al., Paper III) that

 (14)

where is given by (Wehrse et al., Paper II, Eq. (21)) by:

 (15)

We find it is advantageous to remove the  dependence and define . We write now

 (16)

and we write similarly for the diffusion in the moving media:

 (17)

where:

 (18)

and .

Define now an effective opacity as:

 (19)

where is the correction function. should be calculated from a given list of lines, if this is given, or else under the assumption of a Poisson Point Process for example. Note that  for spectral lines.

The function is given in Fig. 6 of Wehrse et al. (2003) in the case of Poisson distributed spectral lines and for a range of parameters for the spectral lines. The sensitivity to the parameters of the spectral lines (as is clear from the figure) is small and will be neglected here. The general behavior of  is as follows: for there is no significant effect, while for the function reaches saturation at about 102-104 times the value at (cf. Table 2 ibid.). In the limit of very large velocity gradients we get:

 (20)

for Poisson distributed lines of equal strength with Lorentz profile and Dirac  line profile.

In the limit of small velocity gradients we can write .

### 4 The diffusion flux

Clearly, the diffusion flux is very non-isotropic. Expanding the flux in the direction  in terms of w yields:

 (21)

where index n implies flux in the direction and

 (22)

Note the particular form of the expansion, namely functions of  times functions of w. To obtain the flux vector one has to integrate Eq. (21) over the solid angle to get (cf. Wehrse et al., Paper IV):

 (23)

with the matrices:

 (24)

 (25)

and

 (26)

namely, only the symmetric part in the Jacobian contributes to the effect.

### 5 Application to disk shear flow

The zeroth order disk structure is a Keplerian shear flow, namely

 (27)

so that for k=1,2 can be calculated for the AD. is the Schwarzchid radius. First we note that there are 4 main directions in which a moving photon does not experience any effect: up and down, radially outward and inward. In all other directions the photon sees a velocity shear. Since the effect depends on w2 there is no difference between the outgoing or ingoing photon. Once the velocity field is given, we can determine whether w is small or large, namely what order of effect we expect. Next, for Keplerian shear flow  and the term with f1 vanishes upon integration over wavelength (it exists only in the monochromatic case). For ADs,

 (28)

Returning to Eq. (14) we have approximately

 (29)

The estimate for depends on the assumption on the shape of the lines. Typical values of  are shown in Fig. 1. The  was calculated for solar composition material and various number of lines per and Doppler widths (which is equivalent to the temperature). We find extremely large values of   over a large range of parameters, namely it is very frequent that   is extremely large. As the temperature decreases and with it the Doppler width the effect increases. Clearly, the effect increases with the number of lines but for any given Doppler width, with sufficiently large number of lines, the lines start to overlap each other and decreases the effect.

The effect contains three factors, the density, the composition and the geometry. For given density and composition the effect is largest near the surface of the accreting star or the Schwarzchild radius in the case of a black hole. The effect is important for

 (30)

where x=r* is the radius of the object in units of the Schwarzschild radius for a  object. For a  compact Neutron star . Hence the condition becomes:

 (31)

We find that for most parameter combinations the effect of the expansion opacity reaches saturation, namely it reaches the maximum (cf. Eq. (20) is the total line absorption  which can be as high as 103 times the Thomson continuum.

 Figure 1: The dependence of on the number of lines per and on the Doppler width. The logarithm of the Doppler width is shown on the left side of the corresponding curve. Open with DEXTER

 Figure 2: The critical density for saturation of the expansion opacity as a function of the location in the disk for an accreting object of  . Open with DEXTER

 Figure 3: The dependence of the shear as a function of direction. Open with DEXTER

### 6 The effect on the structure of the disk

The enhancement of the opacity is very anisotropic. In Fig. 3 we show the geometrical factor  as a function of  at the location (1,0,0). In reality, this term is multiplied by which accentuates the unisotropy.

Because of the extremely large , saturation is reached very quickly and radiation escapes in the directions of smallest opacity. Namely, where is the solid angle in the direction of unchanged opacity, F0 is the would be radiative flux if no expansion opacity effects existed and  is the specific intensity.

The extremely large plus the very anisotropic leads to a situation where the effect enters already at very small shear. In other words, already in a very small deviation from the primary directions where the shear vanishes, the effect is at its maximum value so that  is extremely small.

As a consequence, the temperature gradient increases and leads right away to convective instability. This instability can be the source of macroscopic viscosity in the above ADs. The temperature gradient is given by

 (32)

where the indices w,0 refer to the case with velocity shear and the case without velocity shear respectively.

Clearly, the problem of an ideal ADs with only continuum opacity remains. The shear in the velocity has no effect when spectral lines are absent.

In the particular case of ADs we have to substitute the divergence of the total flux vector into the energy equation. Under the present approximation we return to Eq. (23) which transforms to

 (33)

where Q(r,z) is the energy source and f(r,z) a known function of position. Since the total energy radiated away is essentially fixed (for fixed ) the increase in increases but does not increase the location of the surface and the overall averaged  . (Saturation prevents the singularity.)

### 7 Conclusions

While the velocity gradients keep the effect in the limit of low w, the large  drives the effect into saturation.

We assumed that averaging over angle is plausible (as first order). The strong anisotropy calls for an improved treatment of the radiation field. This extension is underway and will be reported in due course along with the stability analysis.

The limb darkening effect of the disk must be very peculiar with significant consequences to the estimated total luminosity of disks observed at various angles. The expected structure of the atmosphere with non isotropic effective opacity differs significantly from hitherto assumed structure (work in progress).

Magnetic fields, which were ignored here, add to the velocity fields and hence broaden the lines sooner without affecting the global Keplerian flow and general effect.

Acknowledgements
G.S. acknowledges partial support by the Israel Science Foundation and by the DFG (Sonderforschungsberiech 359/C2)