text:
# Theory of Turbulent Accretion Disks

## 1 Introduction

It is an observational fact that at least part of the mass in protostellar disks is accreted onto the central objet. In a Keplerian accretion disk, the specific angular momentum increases with radius. Therefore, a particle can be accreted only if its angular momentum is removed or transferred to particles located at lager radii. Whether angular momentum is removed from or redistributed within the disk depends on whether the disk is subject to external or internal torques. Possible external torques can either be magnetic (when an outflow is present) or tidal (in binary systems), whereas possible internal torques can either be gravitational (massive disks) or turbulent. These mechanisms have been discussed by Papaloizou & Lin (1995).

During the early stages of disk evolution, when the disk is still embedded (class 0/1 object) and has a significant mass compared to the central star, there may exist strong disk winds and bipolar outflows (e.g. Reipurth et al. 1997) with associated magnetic fields. During this stage, a hydromagnetic disk wind may be an important means of angular momentum removal for the system (Blandford & Payne 1982). Because of the action of magnetic torques, material ejected from the disk is able to carry away significantly more angular momentum than it originated with in the disk. Therefore, even a modest ejection rate can lead to a significant accretion rate through the disk. However, observations indicate that outflows may exist only in the early stages of disk evolution, so that this mechanism cannot account for angular momentum transport during the whole life of accretion disks. In addition, it may affect only the inner parts.

When the mass of the disk is significant compared to that of the star, gravitational instabilities may develop, leading to outward angular momentum transport (Papaloizou & Savonije 1991; Heemskerk et al. 1992; Laughlin & Bodenheimer 1994; Pickett et al. 1998) that results in additional mass growth of the central star. This redistribution of mass may occur on the dynamical timescale (a few orbits) of the outer part of the disk and thus may be quite rapid: on the order of $`10^5`$ yr for a disk radius of 500 AU. The parameter governing the importance of disk self–gravity is the Toomre parameter, $`QM_{}H/(M_dr)`$, with $`M_{}`$ being the central mass, $`M_d`$ being the disk mass contained within radius $`r`$ and $`H`$ being the disk semi–thickness. Typically $`H/r0.1`$ (Stapelfeldt et al. 1998) such that the condition for the importance of self–gravity, $`Q1`$, gives $`M_d0.1M_{}`$. During the period of global gravitational instability, it is reasonable to suppose that the disk mass is quickly redistributed and reduced by accretion onto the central object such that the effects of self–gravity become negligible.

If the disk surrounds a star which is in a pre–main sequence binary system, tidal torques transport angular momentum outward, from the disk rotation to the orbital motion of the binary. However, although tidal effects are important for truncating protostellar disks and for determining their size, it is unlikely that tidally–induced angular momentum transport plays a dominant role in the evolution of protostellar disks (see Terquem 2001 and references therein). In a non self–gravitating disk, the amount of transport provided by tidal waves is probably too small to account for the lifetime of protostellar disks. In addition, tidal effects tend to be localized in the disk outer regions.

When the disk mass is such that self–gravity can be ignored and the jet activity has significantly decreased, turbulent torques may become the most important way of redistributing angular momentum in the disk. Historically, the first angular momentum transport mechanism to be considered was through the action of viscosity (von Weizsäcker 1948). However, in order to result in evolution on astronomically interesting timescales, it is necessary to suppose that an anomalously large viscosity is produced through the action of some sort of turbulence.

In this paper we present the theory of turbulent accretion disks. We first calculate the molecular viscosity in a disk, and show this is far too small to account for the evolutionary timescale of disks. We then describe how turbulence may result in enhanced transport of (angular) momentum. About ten years ago, a linear magnetic instability was identified (Balbus & Hawley 1991) that results in MHD turbulence. We present this instability. We then describe both the basis and the structure of $`\alpha `$ disk models, which were introduced by Shakura & Sunyaev (1973).

## 2 Fundamental of hydrodynamics

### 2.1 Fluid mechanics vs. kinetic theory

The fundamental equations of fluids can be derived by considering them as either a collection of particles (kinetic theory) or as a smooth continuum. This latter approach is justified when the mean free path $`\lambda `$ of the particles is very small compared to the macroscopic length scale $`L`$ of interest in the fluid. This condition is not met in stellar systems (galaxies or star clusters), where $`\lambda L`$, nor in planetary rings, where $`\lambda L`$, but it is fulfilled in disks. The conservation equations are obtained more straightforwardly when considering the fluid as a continuum, which is the approach we will use here.

This does not lead to a precise expression for the transport coefficients, however, in contrast to kinetic theory. Transport of energy, mass and momentum occurs in a gas which is out of equilibrium (i.e. in which the distribution function is not a Maxwell–Boltzmann distribution) through molecular collisions. Most of the time, the departure from equilibrium is tiny, so that the distribution function is nearly maxwellian. Within the context of the kinetic theory, in which molecular collisions are explicitly calculated, the Chapman–Enskog procedure gives the transport coefficients by considering small variations of the distribution function around the Maxwell–Boltzmann distribution. Such a calculation is not possible when fluids are viewed as continua, as in this case molecular collisions are not explicitly calculated. It is possible however to get a phenomenological expression for the transport coefficients in this context, as we shall see below.

For more details on the kinetic theory of gases, the reader is referred to, e.g., Huang (1987) and Shu (1992). For the description of fluids as continua, which is developed below, a more complete presentation can be found, e.g., in Acheson (1990), Batchelor (1967), Landau & Lifchitz (1987) and Shu (1992).

### 2.2 Conservation laws

#### 2.2.1 Mass conservation

If there is no sink or source of matter, the time rate of mass variation within a fixed volume $`V`$ is equal to minus the flux of mass (mass advected by the flow) through the surface $`S`$ which encloses the volume:

$$\frac{}{t}_V\rho 𝑑V=_S\rho 𝐯𝐧𝑑S,$$
(1)
where $`\rho (𝐫,t)`$ and $`𝐯(𝐫,t)`$ denote the density and velocity vector of the fluid at location $`𝐫`$ and time $`t`$, and $`𝐧`$ is the unit normal to the surface oriented outward.

Since this has to be satisfied for any volume $`V`$, and using the divergence theorem, we get:

$$\frac{\rho }{t}+\mathbf{}\left(\rho 𝐯\right)=0.$$
(2)
The mass conservation equation derived above is also know as the continuity equation. It states that mass variation in a fixed volume is due to a continuous flow of mass through the surfaces which enclose the volume, and not to ’jumps’ from one location to another remote one.

Equation (2) has been derived by considering the variation of mass of a volume fixed in space. There are two contributions to the time rate of mass change in this volume: $`\rho \mathbf{}𝐯`$, which is due to the change of the volume itself, and $`𝐯\mathbf{}\rho `$, which is due to the mass flowing through the volume.

We can also consider the time rate of mass change of a volume of fluid that we follow along its trajectory. For that we define the Lagrangian derivative $`d/dt`$, which is the time rate of change in a frame moving with a particular fluid element along its trajectory. The density $`\rho `$ depends both on $`𝐫`$ and $`t`$ and, if $`𝐫`$ is the location of the moving fluid element, it depends also on $`t`$. Then

$$\frac{d\rho }{dt}=\frac{\rho }{t}+\frac{\rho }{x_i}\frac{dx_i}{dt}=\frac{\rho }{t}+𝐯\mathbf{}\rho ,$$
(3)
were the $`x_i`$ ($`i=1,2,3`$) denote the coordinates and we adopt the standard Einstein notation of summation over repeated indices.

Equation (2) can therefore be rewritten as:

$$\frac{d\rho }{dt}+\rho \mathbf{}𝐯=0.$$
(4)
This expresses the fact that the mass change in a volume moving with the fluid is due only to the change in the volume.

#### 2.2.2 Incompressibility

A fluid is said to be incompressible if a volume element is neither compressed nor dilated when moving with the flow. Let us consider a fixed volume $`V`$ of the fluid bounded by the surface $`S`$. The volume of fluid which leaves $`V`$ through an infinitesimal surface $`dS`$ per unit time is $`𝐯𝐧dS`$, where $`𝐧`$ is the unit normal to the surface oriented outward. Therefore

$$\frac{dV}{dt}=_S𝐯𝐧𝑑S=_V\left(\mathbf{}𝐯\right)𝑑V,$$
(5)
where we have used the divergence theorem to get the final expression. The condition for the fluid to be incompressible is therefore:

$$\mathbf{}𝐯=0,$$
(6)
or, equivalently, using the mass conservation equation (4):

$$\frac{d\rho }{dt}=0.$$
(7)
This expresses the fact that, in an incompressible fluid, the density of a fluid element is constant along its trajectory.

#### 2.2.3 Equation of motion and momentum conservation

The product of the acceleration of a fluid element moving with the flow, $`d𝐯/dt`$, by the mass density $`\rho `$ is equal to the net force per unit volume exerted onto this element. This has contribution from external forces (e.g., gravitational, magnetic…) and/or from momentum transport within the flow due to molecular collisions (pressure and viscous forces). The local form of the equation of motion is then:

$$\rho \frac{d𝐯}{dt}=\rho \left(\frac{𝐯}{t}+\left(𝐯\mathbf{}\right)𝐯\right)=\mathbf{}P+𝐅_{visc}+𝐅,$$
(8)
where $`P`$ is the pressure, $`𝐅_{visc}`$ is the viscous force per unit volume and $`𝐅`$ is any other force per unit volume which does not arise from pressure and viscosity. When $`𝐅_{visc}=\mathrm{𝟎}`$, this equation is known as Euler equation.

Using equation (2), we can rewrite the equation of motion under the form:

$$\frac{}{t}\left(\rho 𝐯\right)=\mathbf{}\left(\rho \mathrm{𝐯𝐯}\right)\mathbf{}P+𝐅_{visc}+𝐅,$$
(9)
where, according to our notation, the $`i`$ component of the vector $`\mathbf{}\left(\rho \mathrm{𝐯𝐯}\right)`$ is $`\mathbf{}\left(\rho v_i𝐯\right)`$.

This is the momentum conservation equation, which states that the time rate of momentum variation within a fixed volume is equal to minus the momentum which is advected by the flow through the surface which encloses the volume plus the net force exerted onto the volume.

### 2.3 Viscous forces: Molecular transport of momentum

The viscous force in a fluid is due to an irreversible transport of momentum from regions where the velocity is higher to regions where it is lower. A variational calculation shows that uniform rotation is an extremum for the energy. Sometimes it is a minimum, in which case dissipation of energy at fixed (angular) momentum results in uniform rotation (e.g. synchronous binaries). This tends to be the case when pressure is dominant against gravity. Sometimes the extremum is a maximum, in which case all the matter tends to accumulate at the center, with all the angular momentum at infinity. This tends to be the case when rotation is dominant against gravity, and this is what happens in disks (for further details, see Lynden–Bell & Kalnajs 1972; Lynden–Bell & Pringle 1974). The dissipation of energy results from molecular collisions, or friction.

We denote by $`\sigma _{ij}`$ the flux of the $`i`$ component of the momentum in the $`j`$ direction. The tensor $`\left[\sigma \right]`$ is called the viscous stress tensor. The $`i`$ component of the viscous force per unit volume is then:

$$F_{visc,i}\frac{\sigma _{ij}}{x_j}.$$
(10)
Friction occurs only when different parts of the fluid have different velocities. Therefore $`\sigma _{ij}`$ should depend on the velocity gradients. If the velocity varies on a scale large compared to the mean free path, i.e. to the scale over which molecular transport arises, one can suppose that $`\sigma _{ij}`$ depends only on the first derivatives of the velocity with respect to the coordinates. Furthermore, we suppose that the dependence is linear, i.e. we limit ourselves to Newtonian fluids. The most general form of $`\sigma _{ij}`$ is then:

$$\sigma _{ij}\frac{v_i}{x_j}+A\frac{v_j}{x_i}+B\delta _{ij}\frac{v_k}{x_k},$$
(11)
where $`A`$ and $`B`$ are constants to be determined, and $`\delta _{ij}`$ is the Kronecker symbol. If the flow is uniformly rotating in the $`(xy)`$–plane for instance, we must have $`\sigma _{xy}=0`$. Since $`v_xy`$ and $`v_yx`$ in that case, that implies $`A=1`$. Therefore the tensor $`\left[\sigma \right]`$ is symmetrical. We note that:

$$Tr\left[\sigma \right]=\left(2+3B\right)\mathbf{}𝐯,$$

which shows that $`Tr\left[\sigma \right]`$ is a measure of the volume change of a fluid element. It is an experimental fact that the stresses which change the volume of a fluid element give different viscous forces than the stresses that preserve the volume. Therefore we rewrite $`\sigma _{ij}`$ under the form:

$$\sigma _{ij}\left(\frac{v_i}{x_j}+\frac{v_j}{x_i}\frac{2}{3}\delta _{ij}\mathbf{}𝐯\right)+b\delta _{ij}\mathbf{}𝐯,$$

where the term in bracket on the right hand side is trace free, i.e. does not modify the volume of a fluid element.

The shear and bulk viscosities, that we denote $`\eta `$ and $`\zeta `$ respectively, are then experimentally defined as:

$$\sigma _{ij}=\eta \left(\frac{v_i}{x_j}+\frac{v_j}{x_i}\frac{2}{3}\delta _{ij}\mathbf{}𝐯\right)+\zeta \delta _{ij}\mathbf{}𝐯.$$
(12)
The bulk viscosity is associated with internal degrees of freedom of the molecules in the fluid. It becomes negligible if the equipartition between these different degrees of freedom is reached over a timescale shorter than the timescale between two collisions. Furthermore, for a perfect monoatomic gas it can be shown that $`\zeta =0`$ (Huang 1987). Therefore, from now on we shall ignore it.

With this expression for $`\sigma _{ij}`$, we can rewrite the $`i`$ component of the momentum equation (9) under the form:

$$\frac{\left(\rho v_i\right)}{t}+\frac{}{x_j}\left(\rho v_iv_j+T_{ij}\right)=F_i,$$
(13)
where

$$T_{ij}=P\delta _{ij}\sigma _{ij}.$$
(14)
Equation (13) makes is clear that the flux of the $`i`$ component of the momentum in the $`j`$ direction, $`\rho v_iv_j+T_{ij}`$, has contributions from both advection and molecular transport (pressure and viscous forces).

For an incompressible fluid in which $`\eta `$ is constant, the corresponding equation of motion, derived from equation (8) with the above expression of the stress tensor, is called the Navier–Stokes equation.

We note that since friction converts mechanical energy into heat, the rate of change of the kinetic energy, $`dE_k/dt`$, must be negative if only viscous forces are acting. That can be shown from equation (13) to be equivalent to the condition $`\eta >0`$ (e.g., Landau & Lifchitz 1987, § II.16).

### 2.4 Expression of the shear viscosity

The expression (12) for the stress tensor need not be exact. It has been derived phenomenologically assuming that $`\sigma _{ij}`$ depends only on a linear combination of the first derivatives of the velocity with respect to the coordinates. However, the kinetic theory applied to dilute gases leads to the same expression for $`\sigma _{ij}`$. Here we are going to derive an approximate expression for the shear viscosity $`\eta `$ by using the kinetic theory of gases, and we will point out that it arises from correlations between the particle velocities.

#### 2.4.1 Kinetic theory of gases

For simplicity, let us assume that the velocity in cartesian coordinates is $`𝐯=v_x(y)𝐞_x`$, where $`𝐞_x`$ is a unit vector along the $`x`$–axis. The expression (12) for the stress tensor then gives:

$$\sigma _{xy}=\sigma _{yx}=\eta \frac{dv_x}{dy}.$$
(15)
This relation had already been proposed by Newton and Hooke, but it is Maxwell who gave the derivation of the viscosity $`\eta `$ that we develop now.

In average, a molecule has a collision with another molecule after it travels through a distance $`\lambda `$, which is the mean free path of the particles. We suppose that after the collision, the momentum of the molecule is the same as that of its new environment.

Let us consider the momentum which is transported during the time $`\delta t`$ across a surface element $`\delta S`$ perpendicular to the $`y`$–axis and with ordinate $`y`$. There are $`nu\delta t\delta S/6`$ particles crossing $`\delta S`$ from above during $`\delta t`$, where $`n`$ is the number density of particles, $`u`$ is their random (thermal) velocity relative to the mean flow and the factor 6 comes about because there are three possible directions for the particles, each with two orientations. Each of these particles travel through $`\lambda `$ before it suffers a collision below $`\delta S`$, which results in its momentum varying by:

$$m\left[v_x(y)v_x(y+\lambda )\right]m\lambda \frac{dv_x}{dy},$$

to first order in $`\lambda /L`$, where $`L`$ is the scale of variation of the velocity. In other words, each particle carries below $`\delta S`$ the excess of momentum $`m\lambda dv_x/dy`$. Here $`m`$ is the mass of a particle. On the other hand, each particle traveling upward carries above $`\delta S`$ the deficit of momentum $`m\lambda dv_x/dy`$. Therefore, the net $`x`$–component of the momentum which is carried downward during $`\delta t`$ by the particles crossing $`\delta S`$ is:

$$\delta ^2p_x=2\left(\frac{1}{6}nu\delta t\delta S\right)\left(m\lambda \frac{dv_x}{dt}\right)=\frac{1}{3}nmu\lambda \frac{dv_x}{dy}\delta S\delta t.$$
(16)
Let us now consider a box with horizontal faces at $`y`$ and $`y+\delta y`$ and surface area $`\delta S`$. The exchange of particles across the upper face during the time $`\delta t`$ results in the momentum $`\delta ^2p_x(y+\delta y)`$ being added to the volume, whereas the exchange across the lower surface results in the momentum $`\delta ^2p_x(y)`$ being removed from the volume. Therefore, the time rate of change of the momentum content of the box is:

$$\frac{1}{\delta t}\left[\delta ^2p_x(y+\delta y)\delta ^2p_x(y)\right]\frac{d}{dy}\left(\frac{1}{3}nmu\lambda \frac{dv_x}{dy}\right)\delta S\delta y$$
(17)
to first order in $`\lambda /L`$. This is also $`F_{visc,x}\delta S\delta y`$, where $`F_{visc,x}`$ is the viscous force per unit volume. Therefore:

$$F_{visc,x}\frac{d}{dy}\left(\frac{1}{3}nmu\lambda \frac{dv_x}{dy}\right).$$
(18)
Since $`F_{visc,x}=d\sigma _{xy}/dy`$ and using equation (15), this leads to:

$$\eta \frac{1}{3}\rho u\lambda ,$$
(19)
where $`\rho =nm`$ is the mass density of the particles. We also define the kinematic viscosity $`\nu \eta /\rho `$, i.e.

$$\nu \frac{1}{3}u\lambda .$$
(20)
The thermal velocity $`u`$ is on the order of the sound speed.

We denote by $``$ the Reynolds number, which measures the relative strength of the inertial and viscous forces. The inertial force is $`\rho v^2/L`$, where $`v`$ is the flow velocity, and the viscous force is $`\rho \nu v/L^2`$ (from eq.  and ). Therefore,

$$=Lv/\nu .$$
(21)
#### 2.4.2 Velocity correlation

Here again, for simplicity, we consider the case where the flow velocity is along the $`x`$–axis and the motion is in the $`(xy)`$–plane. The components of the instantaneous velocity of a particle in the fluid are $`(v_x+u_x,u_y)`$, where $`u_x`$ and $`u_y`$ are the components of the random (thermal) velocity relative to the mean flow. We have $`<u_x>=<u_y>=0`$, where the brackets denote a time average. The flux of the $`x`$–component of the momentum along the $`y`$–direction is:

$$nm\left(v_x+u_x\right)u_y.$$

Averaged over a large number of particles, or, equivalently, over time, this gives:

$$\rho u_xu_y,$$

since $`<v_xu_y>=v_x<u_y>=0`$. By construction, this quantity is a component of the stress tensor. Therefore,

$$\sigma _{xy}=\rho u_xu_y,$$
(22)
where the minus sign comes about because of the definition of the viscous stress tensor (see § 2.3). For a Maxwell–Boltzmann distribution function (or any $`xy`$ symmetric function), $`<u_xu_y>=0`$ and there is no transport. Equations (15) and (22) lead to:

$$u_xu_y=\nu \frac{dv_x}{dy}.$$
(23)
Note that in general the time average of the fluctuation velocity may not be zero. In an accretion disk for instance, there is a net radial drift of mass, i.e. the mean value of the radial component of the fluctuation velocity is finite. Here, if $`<u_y>`$ were non zero for instance, the flux of the $`x`$–component of the momentum along the $`y`$–direction would have the extra contribution $`\rho v_x<u_y>`$. This represents the transport of mean momentum by the fluctuations. Formally however, what we call the stress tensor would still be given by equation (22), as this is the only contribution from friction between particles.

We denote by $`𝒞`$ the correlation coefficient between the velocities $`u_x`$ and $`u_y`$:

$$𝒞\frac{u_xu_y}{u^2},$$

where $`u`$ is the characteristic velocity of the particles in the fluid relative to the mean flow. We see from (22) that $`𝒞`$ gives a measure of the stress tensor. If $`c_s`$ is the sound speed, then $`uc_s`$. Therefore

$$𝒞\frac{\nu }{c_s^2}\frac{dv_x}{dy}\frac{\lambda }{L}\frac{v_x}{c_s}=\frac{\lambda }{L},$$

where we have used equation (20) and $`dv_x/dyv_x/L`$. Here $`=v_x/c_s`$ is the Mach number. Using equations (20) and (21), we can write $`=L/\lambda `$, and therefore $`𝒞=^2/`$.

In an accretion disk, $`10`$–20 and the Reynolds number corresponding to the molecular viscosity is larger than $`10^{14}`$, so that $`𝒞10^{12}1`$!!

This is a consequence of the small value of the ratio of the mean free path to the scale of the mean flow ($`10^{12}`$) and it means the state of the gas is not affected at all by the molecular transport of angular momentum. In other words, the flow and the molecular transport are completely decoupled.

## 3 Angular momentum transport by turbulence

The realization that molecular transport of angular momentum is so inefficient led the theorists to look for another mechanism of transport in accretion disks. Because Reynolds numbers are so high, it was thought that probably accretion disks would be subject to the same hydrodynamical nonlinear instabilities that shear flows experience in laboratory. The resulting turbulence would then transport angular momentum efficiently. Although today much doubt has been cast on hydrodynamical instabilities in disks, turbulence is still a strong candidate for transport since it has been shown relatively recently that a linear magnetohydrodynamical instability can develop in disks (see below). Therefore, we turn now to turbulent transport, and contrast it with molecular transport. Much of this section is based on Tennekes & Lumley (1972), which the reader is referred to for more details (see also Balbus & Hawley 1998).

We shall restrict ourselves here to the study of incompressible flows, as this simplifies the discussion. For our argument, it is sufficient to take into account only pressure and viscous forces, but in principle any other (external or inertial) force could be added. The equations describing the fluid are:

$$\frac{\stackrel{~}{v}_i}{t}+\stackrel{~}{v}_j\frac{\stackrel{~}{v}_i}{x_j}=\frac{1}{\rho }\frac{}{x_j}\stackrel{~}{\sigma }_{ij},$$
(24)
and

$$\frac{\stackrel{~}{v}_i}{x_i}=0,$$
(25)
where the tilde symbol above a variable means we consider the instantaneous value of the variable at the location $`x_i`$ and time $`t`$. Here

$$\stackrel{~}{\sigma }_{ij}=\stackrel{~}{p}\delta _{ij}+\eta \stackrel{~}{s}_{ij},$$

with

$$\stackrel{~}{s}_{ij}=\frac{\stackrel{~}{v}_i}{x_j}+\frac{\stackrel{~}{v}_j}{x_i}$$

and $`\stackrel{~}{p}`$ is the pressure. We suppose $`\eta `$ is a constant.

We now use the so–called Reynolds decomposition, in which an instantaneous value is written as the sum of a mean value (denoted by a capital letter) plus a fluctuation (denoted by a small letter):

$$\stackrel{~}{v}_i=V_i+v_i.$$

This fluctuation is characteristic of the turbulence. This decomposition is meaningful only if the timescale on which the fluctuations vary and the evolution timescale of the flow are well separated. The mean values are then taken over a timescale large compared to the turbulence timescale but short compared to that of the flow evolution. As here we are not interested in the long term evolution of the flow, we neglect the derivative with respect to time of the mean values (i.e. we consider a quasi–steady state). To simplify the discussion, we suppose that the average of $`v_i`$ over time is zero: $`<v_i>=0`$.

Equation (25) averaged over time then leads to $`V_i/x_i=0,`$ i.e. the mean flow is incompressible. Equation (25) thus implies $`v_i/x_i=0,`$ i.e. the fluctuations are also incompressible.

Using

$$\stackrel{~}{\sigma }_{ij}=\mathrm{\Sigma }_{ij}+\sigma _{ij}$$

with $`\sigma _{ij}=0`$, the equation of motion averaged over time gives:

$$\frac{}{x_j}V_jV_i+\frac{}{x_j}v_jv_i=\frac{1}{\rho }\frac{}{x_j}\mathrm{\Sigma }_{ij},$$
(26)
where we have used $`V_i/t=0`$ and the incompressibility of the mean flow and the fluctuations. The term $`<v_jv_i>`$ represents the averaged transport of the fluctuations of the momentum by the fluctuations of the velocities. This is the turbulent transport. It is non zero only if the turbulent velocities in the different directions are correlated. It is an experimental fact that this is in general the case for shear turbulence (see below). Equation (26) shows that momentum is transferred between the fluctuations and the mean flow through the term $`<v_jv_i>/x_j`$.

We can rewrite equation (26) under the form:

$$\frac{}{x_j}V_jV_i=\frac{1}{\rho }\frac{}{x_j}\left(\mathrm{\Sigma }_{ij}+\tau _{ij}\right),$$
(27)
where $`\tau _{ij}=\rho v_jv_i`$ is called the Reynolds stress tensor, or turbulent stress. This is the contribution from the fluctuations to the averaged total stress tensor $`\mathrm{\Sigma }_{ij}+\tau _{ij}`$. Note that if we had not supposed $`<v_i>=0`$, we would also have had terms like $`V_j<v_i>`$, representing transport (or advection) of mean momentum by the fluctuations. This is what Lynden–Bell & Kalnajs (1972) call lorry–transport, because it is a direct ”shipment” by the equilibrium flow. Formally, this term is not part of what we call the turbulent stress however.

As we have no expression for the components of $`\tau _{ij}`$ (six of which are independent), the problem has more unknowns than equations. This is the well–known closure problem of turbulence. Since $`\tau _{ij}`$ appears in the same way as $`\mathrm{\Sigma }_{ij}`$ in equation (27), it is tempting to express $`\tau _{ij}`$ by analogy with $`\mathrm{\Sigma }_{ij}`$, which depends on the molecular motion as we have seen in the previous section. This is the basis for the mixing length theory, in which $`\tau _{ij}`$ is written exactly like the tensor deriving from molecular motions using a so–called turbulent viscosity $`\nu _T`$. By analogy with expression (20) for the molecular viscosity, it is supposed that $`\nu _Tv_T\mathrm{\Lambda }`$, where $`v_T`$ is a characteristic velocity of the turbulent eddies and $`\mathrm{\Lambda }`$ is the so–called mixing length, which is the “mean free path” of the eddies, i.e. the distance they travel through before they mix with their environment.

The same analogy is used in accretion disk theory through the $`\alpha `$ model that we shall describe in details below.

It is important to note that the basis for this analogy is very weak. For a thorough discussion of the differences between molecular and turbulent transport, we refer the reader to Tennekes & Lumley (1972). Here we recall the main points.

There is a fundamental difference between molecular and turbulent transport, in that turbulence is part of the flow whereas molecules are part of the fluid. We saw above that molecular transport is decoupled from the state of the flow. This is not true for turbulent transport, which depends completely on the flow for its existence. Turbulence can be sustained only if the eddies are able to tap the energy of the mean flow. Observations of laboratory shear flows for instance indicate that the eddies which are most efficient in tapping the energy of the mean flow are also most efficient in preserving a good correlation between the different components of the turbulent velocities, thus giving a large stress tensor. This is very different than for the molecules, whose velocities do not depend on non local properties of the mean flow.

And as a matter of fact, transport of momentum is observed to be a characteristic of turbulent shear flows. In the air for instance, the correlation coefficient between the velocities is on the order of $`10^6`$ for the molecules and $`0.4`$ for turbulent motions. This is a characteristic of most shear flows.

In addition, we saw when we derived the expression of the viscosity (19) that this requires $`\lambda L`$. For a turbulent flow, that would mean that the scale of the turbulent eddies is very small compared to the scale over which the characteristics of the flow, like the mean velocity, vary. This is not necessarily satisfied in an accretion disk, where both may be on the order of the disk semi–thickness.

The conclusion is that shear turbulence is indeed a very efficient way of transporting momentum, but great care should be taken when expressing the resulting stress tensor by analogy with molecular transport.

In particular, while representing gross transport by an effective viscosity can often be useful, doing a detailed stability analysis on a viscous fluid model for turbulent flow is generally not self–consistent, and can be very misleading (e.g. Hawley, Balbus, & Stone 2001).

Note that the above discussion applies to shear turbulence only, i.e. flows where the fluctuations get their energy from the mean velocity gradients. Transport of momentum may not be nearly as efficient when the energy source for the fluctuations is a gradient of temperature or magnetic field for instance. As a matter of fact, there are strong indications that the transport of momentum associated with thermal convection is orders of magnitude weaker than that associated with shear turbulence (Balbus & Hawley 1998 and references therein).

Although turbulence has been considered as a way of transporting angular momentum in accretion disks for more than fifty years, it is only relatively recently that an instability able to extract the energy of the shear and put it in the fluctuations has been elucidated. This instability, which requires a magnetic field, is described in the next section.

## 4 Magnetohydrodynamic instabilities

We consider now the stability of a disk containing a magnetic field. Balbus & Hawley (1991, 1998 and references therein) have shown that in under very general conditions a magnetized disk is subject to a linear instability which nonlinear development was subsequently found to result in turbulence (Hawley, Gammie & Balbus 1995; Brandenburg et al. 1995). Following the presentation by Balbus & Hawley (1998), we first briefly describe the waves which propagate in a fluid containing a magnetic field (more details can be found in, e.g., the textbooks by Jackson 1975, Shu 1992 or Sturrock 1994), and we then go on to describe the linear instability.

### 4.1 Waves propagating in a magnetized fluid

Here we take into account only pressure and magnetic forces, as these are the only important forces for our argument. The governing equations are the induction equation:

$$\frac{𝐁}{t}=\mathbf{}\mathbf{\times }\left(𝐯\mathbf{\times }𝐁\right),$$
(28)
which ensures that $`\mathbf{}𝐁=0`$, the equation of continuity:

$$\frac{\rho }{t}+\mathbf{}\left(\rho 𝐯\right)=0,$$
(29)
and the equation of motion:

$$\frac{𝐯}{t}+\left(𝐯\mathbf{}\right)𝐯=\frac{1}{\rho }\mathbf{}P+\frac{1}{\mu _0\rho }\left(\mathbf{}\mathbf{\times }𝐁\right)\mathbf{\times }𝐁,$$
(30)
where the last term on the right hand side is the Lorentz force per unit volume in SI units ($`\mu _0`$ is the permeability of vacuum). Note that it can also be written:

$$\frac{1}{\mu _0}\left(\mathbf{}\mathbf{\times }𝐁\right)\mathbf{\times }𝐁=\mathbf{}\left(\frac{B^2}{2\mu _0}\right)+\left(\frac{𝐁}{\mu _0}\mathbf{}\right)𝐁$$
(31)
which shows that the magnetic force if the sum of a magnetic pressure and a magnetic tension (first term and second term on the right hand side, respectively).

We first consider a nonrotating fluid and we look for small perturbations to an equilibrium state characterized by the quantities $`\rho _0`$, $`𝐯_0=\mathrm{𝟎}`$, $`P_0`$ and $`𝐁_0`$. We denote the Eulerian perturbations by a prime, so that $`\rho =\rho _0+\rho ^{}`$ etc, and we suppose they vary on a scale much smaller than the scale of variation of the equilibrium quantities. The linearized equations (28), (29) and (30) are then:

$`{\displaystyle \frac{𝐁^{}}{t}}`$ $`=`$ $`\mathbf{}\mathbf{\times }\left(𝐯^{}\mathbf{\times }𝐁_\mathrm{𝟎}\right),`$ (32)
$`{\displaystyle \frac{\rho ^{}}{t}}`$ $`+`$ $`\rho _0\mathbf{}𝐯^{}=0,`$ (33)
$`{\displaystyle \frac{𝐯^{}}{t}}`$ $`=`$ $`{\displaystyle \frac{1}{\rho _0}}\mathbf{}P^{}+{\displaystyle \frac{1}{\mu _0\rho _0}}\left(\mathbf{}\mathbf{\times }𝐁^{}\right)\mathbf{\times }𝐁_\mathrm{𝟎}.`$ (34)
Since the coefficients of the above equations are constant, we perform a Fourier transform for the variables $`𝐫`$ and $`t`$ and look for solutions under the form exp$`\left[i\left(\omega t𝐤𝐫\right)\right]`$ with $`\omega `$ and $`|𝐤|`$ real constants. The above equations can then be rewritten by replacing $`/t`$ by $`i\omega `$ and $`\mathbf{}`$ by $`i𝐤`$.

To close the system of equations, we need an equation of state for the fluid. We consider adiabatic perturbations (i.e. we suppose there is no dissipation of heat), so that:

$$\frac{P^{}}{P}=\gamma \frac{\rho ^{}}{\rho },$$
(35)
where $`\gamma `$ is the ratio of specific heats ($`\gamma =5/3`$ for a perfect gas).

The system of equations (32)–(35) gives a dispersion relation of order six. It can be simplified by noting that the solutions associated with motions perpendicular to the $`(𝐤,𝐁_0)`$ plane and those associated with motions in this plane can be decoupled.

The first type of solutions is an Alfven wave, and its dispersion relation is:

$$\omega ^2=k^2v_A^2\mathrm{cos}^2\psi ,$$
(36)
where $`\psi `$ is the angle between $`𝐤`$ and $`𝐁_0`$ and

$$𝐯_A=\frac{𝐁_0}{\sqrt{\mu _0\rho _0}}$$

is the Alfven speed. These are transverse waves (with no motion in the $`(𝐤,𝐁_0)`$ plane) which propagate along the field lines. They do not involve any compression across the field lines and their restoring force is the magnetic tension. These waves are analogous to the waves which propagate along a stretched string. Within the factor $`\mathrm{cos}\psi `$, we see that the Alfven speed is the square root of the magnetic tension divided by a mass density, just like the phase speed in the case of a string.

The second type of solutions is associated with motions only in the $`(𝐤,𝐁_0)`$ plane. Its dispersion relation is:

$$\omega ^2=\frac{k^2}{2}\left\{\left(v_A^2+c_s^2\right)\pm \left[\left(v_A^2+c_s^2\right)^24v_A^2c_s^2\mathrm{cos}^2\psi \right]^{1/2}\right\},$$
(37)
where $`c_s^2=\gamma P_0/\rho _0`$ is the adiabatic sound speed. The upper sign, which gives a larger phase speed $`\omega /k`$, corresponds to the fast MHD wave, whereas the lower sign corresponds to the slow MHD wave. The fast mode is also referred to as magneto–acoustic wave. If either $`c_sv_A`$, $`v_Ac_s`$ or $`\mathrm{cos}\psi 1`$, then

$$\omega _+^2=k^2\left(v_A^2+c_s^2\right)$$

and

$$\omega _{}^2=\frac{k^2v_A^2c_s^2\mathrm{cos}^2\psi }{v_A^2+c_s^2},$$

where the subscripts + and - denote the fast and slow waves, respectively. For the fast mode, the magnetic and thermal pressures act in phase. If $`v_Ac_s`$, the slow mode is an acoustic wave propagating along the field lines, whereas if $`v_Ac_s`$ it degenerates into an Alfven mode in its dispersion properties (the eigenvector is distinct from that of the Alfven mode however, as the motions are not in the same plane). For the fast mode it is the opposite.

In the absence of rotation, these modes are stable. However, Balbus & Hawley have shown that when rotation is introduced, the slow mode can become unstable, and this what we describe now.

### 4.2 Linear magnetorotational instability

We consider the simplest system in which the instability can develop. This is an axisymmetric disk with a vertical uniform magnetic field. For the original presentation, which includes the case of a radial field, see Balbus & Hawley (1991), and for the stability of a toroidal field see Balbus & Hawley (1992b), Fogglizo & Tagger (1995), Terquem & Papaloizou (1996) and Ogilvie & Pringle (1996).

Since it is the slow mode which is destabilized, one can consider an incompressible fluid. That amounts to taking the limit $`c_s\mathrm{}`$ in the above equations, so that the fast mode disappears (its frequency becomes infinite).

The system of equations describing the fluid is then

$$\frac{𝐁}{t}=\mathbf{}\mathbf{\times }\left(𝐯\mathbf{\times }𝐁\right),$$
(38)
$$\mathbf{}𝐯=0,$$
(39)
$$\frac{𝐯}{t}+\left(𝐯\mathbf{}\right)𝐯=\frac{1}{\rho }\mathbf{}P+\frac{1}{\mu _0\rho }\left(\mathbf{}\mathbf{\times }𝐁\right)\mathbf{\times }𝐁\mathbf{}\mathrm{\Psi },$$
(40)
where $`\mathrm{\Psi }`$ is the gravitational potential due to the central star.

We use the cylindrical coordinates $`(r,\varphi ,z)`$ and we denote the unit vectors in the three directions by $`𝐞_r`$, $`𝐞_\varphi `$ and $`𝐞_z`$. The components of the vector $`\left(𝐯\mathbf{}\right)𝐯`$ which appears in the equation of motion are:

$$v_r\frac{v_r}{r}+\frac{v_\varphi }{r}\frac{v_r}{\varphi }+v_z\frac{v_r}{z}\frac{v_\varphi ^2}{r},$$

$$v_r\frac{v_\varphi }{r}+\frac{v_\varphi }{r}\frac{v_\varphi }{\varphi }+v_z\frac{v_\varphi }{z}+\frac{v_rv_\varphi }{r},$$

$$v_r\frac{v_z}{r}+\frac{v_\varphi }{r}\frac{v_z}{\varphi }+v_z\frac{v_z}{z}.$$

The equilibrium quantities, that we suppose uniform, are $`\rho _0`$, $`P_0`$, and $`𝐁_0=B_0𝐞_z`$. The velocity is $`𝐯_0=r\mathrm{\Omega }(r)𝐞_\varphi `$. We now suppose that this equilibrium is slightly perturbed, so that $`\rho =\rho _0+\rho ^{}`$ etc, where the prime denotes the Eulerian perturbations, and we look for solutions proportional to exp$`\left[i\left(\omega tk_zzk_rr\right)\right]`$. Here we consider axisymmetric perturbations, but more general solutions can be obtained (see the above mentioned papers).

The linearized system of equations is then:

$`i\omega B_r^{}`$ $`=`$ $`ik_zBv_r^{},`$ (41)
$`i\omega B_\varphi ^{}`$ $`=`$ $`ik_zBv_\varphi ^{}ik_zr\mathrm{\Omega }B_z^{}+\left({\displaystyle \frac{d\left(r\mathrm{\Omega }\right)}{dr}}ik_rr\mathrm{\Omega }\right)B_r^{},`$ (42)
$`i\omega B_z^{}`$ $`=`$ $`{\displaystyle \frac{B}{r}}v_r^{}+ik_rBv_r^{},`$ (43)
$`\left({\displaystyle \frac{1}{r}}ik_r\right)v_r^{}ikv_z^{}`$ $`=`$ $`0,`$ (44)
$`i\omega v_r^{}2\mathrm{\Omega }v_\varphi ^{}`$ $`=`$ $`{\displaystyle \frac{ik_rP^{}}{\rho }}{\displaystyle \frac{ik_zBB_r^{}}{\mu _0\rho }}+{\displaystyle \frac{ik_rBB_z^{}}{\mu _0\rho }},`$ (45)
$`i\omega v_\varphi ^{}+\left({\displaystyle \frac{d\left(r\mathrm{\Omega }\right)}{dr}}+\mathrm{\Omega }\right)v_r^{}`$ $`=`$ $`{\displaystyle \frac{ik_zBB_\varphi ^{}}{\mu _0\rho }},`$ (46)
$`i\omega v_z^{}`$ $`=`$ $`{\displaystyle \frac{ik_zP^{}}{\rho }},`$ (47)
where, for simplicity, we have dropped the subscript ’0’ for the equilibrium quantities. We now consider large vertical wavenumbers, such that $`|k_z||k_r|`$ and $`|k_z|1/r`$. Then, the above system of equations leads to the following dispersion relation:

$$\omega ^4\omega ^2\left(2k^2v_A^2+\frac{d\mathrm{\Omega }^2}{d\mathrm{ln}r}+4\mathrm{\Omega }^2\right)+k^2v_A^2\left(k^2v_A^2+\frac{d\mathrm{\Omega }^2}{d\mathrm{ln}r}\right)=0.$$
(48)
There is instability if $`\omega `$ has a negative imaginary part (the perturbation then grows exponentially with time). Equation (48) is a quartic for $`\omega ^2`$ which solutions are real. Therefore there is instability if $`\omega ^2<0`$, which requires:

$$k^2v_A^2<\frac{d\mathrm{\Omega }^2}{d\mathrm{ln}r}.$$
(49)
This criterion has a very simple physical explanation. It states that there is an instability when the magnetic tension that acts on a segment of a field line is smaller than the net tidal force (i.e. centrifugal force minus gravitational force) acting on it.

For a given equilibrium field $`B`$, and therefore a given Alfven speed $`v_A`$, there will always be a wavenumber $`k`$ such that this inequality is satisfied provided the right hand side is positive. Therefore, all the disks with

$$\frac{d\mathrm{\Omega }^2}{d\mathrm{ln}r}<0$$
(50)
are unstable, and this is the criterion for instability.

The heart of the instability resides in the fact that a perturbed fluid element tends to conserve its angular velocity when a magnetic field is present. This is to be contrasted with a non magnetized disk, in which a perturbed element tends to conserve its specific angular momentum. When displaced inward therefore, it has too much angular momentum for its new location (as the angular momentum increases outward in an accretion disk), and it moves back to its initial unperturbed position. When a magnetic field is present, the magnetic tension along the field line tends to enforce isorotation of the elements to which it is connected. A fluid element displaced inward has therefore a lower angular velocity than the elements at its new location, and thus not enough angular momentum for its new position. As a result it sinks further in. On the opposite, a fluid element connected to the same field line and displaced outward will tend to move further out. Angular momentum is transferred via magnetic torques from the inner fluid element to the outer fluid element. Note that the source of free energy for the instability is not in the magnetic field, but in the disk differential rotation. The magnetic field just provides a path to extract the energy.

From equation (48), we can write the negative values of $`\omega ^2`$ as a function of $`k^2v_A^2`$, and show that the maximum growth rate is:

$$\left|\omega _{\mathrm{max}}\right|=\frac{1}{2}\left|\frac{d\mathrm{\Omega }}{d\mathrm{ln}r}\right|.$$
(51)
For a Keplerian disk,

$$\left|\omega _{\mathrm{max}}\right|=\frac{3}{4}\mathrm{\Omega }$$
(52)
and is attained for

$$kv_A=\frac{\sqrt{15}}{4}\mathrm{\Omega }.$$
(53)
This holds even if the field has radial and azimuthal components and if the perturbed quantities are allowed to vary with $`r`$ and $`\varphi `$ provided we then replace $`kv_A`$ by $`𝐤𝐯_A`$. Note that the non–axisymmetric case is more subtle though, as in that case plane waves cannot be sustained, being sheared out by the differential rotation. If we write $`𝐤=k_r𝐞_r+(m/r)𝐞_\varphi +k_z𝐞_z`$, with $`m`$ being the azimuthal wavenumber, then $`k_r`$ is time dependent and

$$k_r(t)=k_0mt\frac{d\mathrm{\Omega }}{dr}$$

(Goldreich & Lynden–Bell 1965), which means that a disturbance always becomes trailing in a disk where the angular velocity decreases outward. If $`k_r`$ is initially large and positive (leading disturbance), then the mode is stable as indicated by (49). But as time goes on, $`k_r`$ decreases, so that $`𝐤`$ enters a region of instability. As $`k_r`$ becomes negative and keeps decreasing though, the mode becomes stable again. Formally, the instability is therefore not purely exponential. The important question however is whether the mode can be amplified significantly before its wavelength becomes small enough to be affected by ohmic resistivity. This, of course, depends on the magnetic Reynolds number.

We have seen above that the so–called magnetorotational or Balbus–Hawley instability can develop in any disk in which the angular velocity decreases outward. In principle there is no other condition. However, it may be that the scale of the modes which are unstable according to (49) do not fit into the disk, i.e. they have a wavelength larger than the disk semithickness $`H`$. This is the case if $`v_A/\mathrm{\Omega }>H`$. Since in a thin disk $`\mathrm{\Omega }Hc_s`$ (see § 5.3.1), the disk will be stable if $`v_A>c_s`$.

Another condition which is implicit in the above presentation is, of course, that the magnetic field be coupled to the fluid. This may not be the case everywhere in protostellar disks, which are rather cold and dense (Gammie 1996; Fromang, Terquem & Balbus 2001).

Ohmic resistivity can prevent the development of the instability if the scale on which it acts is comparable to the wavelength of the unstable modes. The time it takes for the magnetic field to diffuse over a scale $`1/k`$ due to the effect of an Ohmic resistivity $`\eta _B`$ is $`1/(\eta _Bk^2)`$. Since in the absence of resistivity it would grow on a timescale $`1/\mathrm{\Omega }`$, we see that Ohmic dissipation prevents the instability if $`\eta _B\mathrm{\Omega }/k^2`$, i.e. if the magnetic Reynolds number is of order unity. For more details, the reader is referred to Balbus & Blaes (1994), Jin (1996) and Papaloizou & Terquem (1997). Note that Fleming, Stone & Hawley (2000) have shown that even though the linear instability is not affected at higher values of the magnetic Reynolds number, the nonlinear instability cannot be sustained for values as high as $`10^4`$ in some circumstances, which depend on the field topology for instance. However, the issue is not yet settled as other effects like Hall electromotive forces may resurrect the instability in regions where Ohmic resistivity acting alone would inhibit it (Wardle 1999, Balbus & Terquem 2001).

### 4.3 Angular momentum transport in magnetized disks

Numerical simulations have shown that this instability puts the energy it extracts from the disk differential rotation into fluctuations which transport angular momentum outward (see Balbus & Hawley 1998 and references therein). As noted above, transport of momentum is not a property of all types of turbulence. Furthermore, when it happens, it is not always against the velocity gradient, i.e. outward in an accretion disk. That has to be the case in an isolated dissipative system, but not if the energy driving the turbulence has an external source. We have already mentioned that the transport associated with thermal convection appears to be very weak compared to that associated with shear turbulence. In addition, there are strong indications that convection driven by external heating transports angular momentum inward in accretion disks (Ryu & Goodman 1992, Cabot 1996, Stone & Balbus 1996). In the case of the magnetorotational instability though, it can be shown that in the linear regime the transport is outward (Balbus & Hawley 1992a). This still holds in the nonlinear regime, as expected since the energy source for the turbulence is internal. Numerical simulations also show that most of the transport is done by the (magnetic) Maxwell stress, which dominates over the (hydrodynamic) Reynolds stress. Furthermore, magnetic fields are regenerated through a dynamo action so that the mechanism can sustain itself in an isolated system.

It is important to realize that this instability and the resulting turbulence do not tend to make the disk angular velocity uniform. This is because the angular velocity is imposed by the gravitational potential of the object around which the disk rotates. The turbulence can smear the shear out only on scales smaller than the disk semithickness $`H`$. To smear the shear out on a scale $`l`$, the turbulent velocity $`v_t`$ has indeed to be on the order of the shear over $`l`$, i.e. $`l|rd\mathrm{\Omega }/dr|l\mathrm{\Omega }`$, where $`\mathrm{\Omega }`$ is the angular velocity in the disk. If $`lH`$, then $`v_tc_s`$, as in a thin disk $`c_sH\mathrm{\Omega }`$ (see § 5.3.1). Therefore, only supersonic turbulence could modify the shear on scales larger than $`H`$. However, such a turbulence could not be maintained as the supersonic fluctuations would dissipate into shocks. The state toward which the disk evolves as a result of the turbulence is not a state of uniform rotation but one where all the mass is at the center and all the angular momentum at infinity (Lynden–Bell & Pringle 1974).

Long before a mechanism for producing turbulence in accretion disks had been identified, Shakura & Sunyaev (1973) proposed a prescription for modeling turbulent disks. We now describe this prescription, and discuss its validity in the context of magnetic turbulence. We then develop disk models based on this so–called $`\alpha `$ prescription.

## 5 $`\alpha `$ disk models

### 5.1 Evolution of turbulent disks

Here we consider an axisymmetric disk rotating around a central object. We suppose the motion is in the plane of the disk, or, equivalently, we use the vertically averaged equations of mass conservation and motion. The velocity is $`𝐯(u_r,r\mathrm{\Omega }+u_\varphi )`$. The term $`r\mathrm{\Omega }`$ is the circular velocity around the central mass, and $`(u_r,u_\varphi )`$ are the components of the fluctuation velocity. Note that as the disk accretes, there is a net radial drift and the mean value of $`u_r`$ is non zero. Here, the equation of mass conservation (2) and the azimuthal component of the equation of motion (8) are:

$$\frac{\mathrm{\Sigma }}{t}+\frac{1}{r}\frac{}{r}\left(r\mathrm{\Sigma }v_r\right)=0,$$
(54)
$$\mathrm{\Sigma }\left(\frac{v_\varphi }{t}+v_r\frac{v_\varphi }{r}+\frac{v_rv_\varphi }{r}\right)=0,$$
(55)
where $`\mathrm{\Sigma }`$ is the surface mass density, i.e. the mass density integrated over the disk thickness. Here, according to the discussion in § 2.4, we have neglected the viscous force arising from molecular viscosity. Multiplying equation (55) by $`r`$ and using (54), we obtain the angular momentum equation:

$$\frac{}{t}\left(r\mathrm{\Sigma }\left(r\mathrm{\Omega }+u_\varphi \right)\right)+\frac{1}{r}\frac{}{r}\left(r^2\mathrm{\Sigma }\left(r\mathrm{\Omega }+u_\varphi \right)u_r\right)=0.$$
(56)
As pointed out by Balbus & Papaloizou (1999), to get a diffusion equation describing the disk evolution we need to smooth out the fluctuations over radius. To do so, we average the above equation over a scale large compared to that of the fluctuations, but small compared to the characteristic scale of the flow (i.e. $`r`$). Equation (56) then yields:

$$\frac{}{t}\left(\mathrm{\Sigma }r^2\mathrm{\Omega }\right)+\frac{1}{r}\frac{}{r}\left(\mathrm{\Sigma }r^3\mathrm{\Omega }<u_r>+\mathrm{\Sigma }r^2<u_ru_\varphi >\right)=0,$$
(57)
where the brackets denote the radial average and we have neglected $`<u_\varphi >`$ compared to $`r\mathrm{\Omega }`$ in the time derivative. This is justified because $`|<u_\varphi >|r\mathrm{\Omega }`$ and the systematic, evolutionary time derivative of $`<u_\varphi >`$ is limited. In the radial divergence term however, both $`<u_r>`$ and $`<u_ru_\varphi >`$ are second order, and therefore we retain all the terms. This equation is the same as that describing a viscous flow with $`𝐯(u_r,r\mathrm{\Omega })`$ and a stress tensor $`\sigma _{r\varphi }\mathrm{\Sigma }<u_ru_\varphi >`$ as given by equation (22).

There are two contributions to the flux of angular momentum: the term $`\mathrm{\Sigma }r^2\mathrm{\Omega }<u_r>`$ is the mean angular momentum advected through the disk by the velocity fluctuations (because of the accretion of mass), whereas the term $`\mathrm{\Sigma }r<u_ru_\varphi >`$ represents the angular momentum fluctuations transported by the velocity fluctuations.

Using equation (54) averaged over radius, we can rewrite equation (57) under the form:

$$r\mathrm{\Sigma }<u_r>=\left[\frac{d}{dr}\left(r^2\mathrm{\Omega }\right)\right]^1\frac{}{r}\left(\mathrm{\Sigma }r^2<u_ru_\varphi >\right).$$

Using equation (54) again to eliminate $`<u_r>`$, this leads to the diffusion equation:

$$\frac{\mathrm{\Sigma }}{t}=\frac{1}{r}\frac{}{r}\left(\left[\frac{d}{dr}\left(r^2\mathrm{\Omega }\right)\right]^1\frac{}{r}\left(\mathrm{\Sigma }r^2<u_ru_\varphi >\right)\right).$$
(58)
In a steady state, equation (54) gives $`r\mathrm{\Sigma }<u_r>=\mathrm{constant}`$. Then the accretion rate

$$\dot{M}2\pi r\mathrm{\Sigma }<u_r>$$
(59)
is constant through the disk. Integration of the angular momentum equation (57) then yields:

$$\mathrm{\Sigma }<u_ru_\varphi >=\frac{\dot{M}}{2\pi }\mathrm{\Omega }\left[1\left(\frac{R_i}{r}\right)^{1/2}\right],$$
(60)
where $`R_i`$ is the disk inner boundary. We have assumed here that the turbulent stress $`<u_ru_\varphi >`$ vanishes at the disk inner edge (i.e. there is no torque at the boundary) and that the disk is Keplerian, so that $`\mathrm{\Omega }r^{3/2}`$. The above equation shows that for the mass to be accreted inward (i.e. $`<u_r>`$ negative), the flux of angular momentum due to the fluctuations has to be positive, i.e. the fluctuations have to transport angular momentum outward.

### 5.2 The $`\alpha `$ prescription

We pointed out that the angular momentum equation (57) is analogous to that describing a viscous flow with $`𝐯(u_r,r\mathrm{\Omega })`$ and a stress tensor $`\sigma _{r\varphi }\mathrm{\Sigma }<u_ru_\varphi >`$. Therefore, it is tempting to push the analogy further and express the turbulent stress $`\mathrm{\Sigma }<u_ru_\varphi >`$ as if it derived from an enhanced ’turbulent viscosity’ $`\nu `$, defined such that (equation ):

$$\mathrm{\Sigma }<u_ru_\varphi >\mathrm{\Sigma }\nu r\frac{d\mathrm{\Omega }}{dr}.$$
(61)
In a Keplerian disk, this gives:

$$<u_ru_\varphi >\nu \mathrm{\Omega }.$$
(62)
The equations presented in the previous section have been derived by Lynden–Bell & Pringle (1974; see also Pringle 1981) assuming a viscous flow with this expression of the stress tensor.

The prescription proposed by Shakura & Sunyaev (1973) consists in writing $`\nu =v_tH`$, where $`H`$ is the disk thickness, assumed to be the maximum scale of the turbulent cells, and $`v_t`$ is the turbulent velocity. They further define

$$\alpha \frac{v_t}{c_s},$$

where $`c_s`$ is the sound velocity. Note that $`\alpha <1`$, otherwise the fluctuations would dissipate into shocks in such a way as to restore $`v_t<c_s`$. Equation (62) can then be rewritten under the form:

$$<u_ru_\varphi >\alpha c_s^2.$$
(63)
Here we have used the fact that in a thin disk $`Hc_s/\mathrm{\Omega }`$ (see § 5.3.1).

So far we have focussed on non magnetized disks. In these, there are strong indications that $`<u_ru_\varphi >`$ is either zero or negative. Magnetism is needed to correlate the velocities. The above discussion does apply to magnetized disks provided we replace $`<u_ru_\varphi >`$ by $`<u_ru_\varphi u_{Ar}u_{A\varphi }>`$, where $`(u_{Ar},u_{A\varphi })`$ are the components of the fluctuations of the Alfven velocity (Shakura & Sunyaev 1973, Balbus & Hawley 1998). The extra term represents the Maxwell stress.

The validity of the $`\alpha `$ prescription in the context of magnetic turbulence was discussed by Balbus & Papaloizou (1999). They first pointed out that, as long as $`<u_ru_\varphi >`$ (or $`<u_ru_\varphi u_{Ar}u_{A\varphi }>`$) is positive, the disk dynamics is the same as if it were evolving under the action of a viscosity. In that case indeed, the diffusion coefficient in equation (58) is positive. We can then always define an $`\alpha `$ parameter according to (63), although $`\alpha `$ may not be constant through the disk.

Note however that this implicitly assumes it is possible to average the equations over radius in the way described in § 5.1. If the scale of the fluctuations and the characteristic disk scale are not well separated, such an average cannot be performed. Since the maximum scale of the fluctuations is of order the disk thickness $`H`$, this procedure requires that there is a scale large compared to $`H`$ and small compared to $`r`$. This condition may be only marginally fulfilled in protostellar disks, in which $`H`$ may be up to 0.1–$`0.2r`$.

Balbus & Papaloizou (1999) further noted that for the $`\alpha `$ prescription to apply, the disk had to behave viscously not only in its dynamics but also in its energetics. The key point here is that a viscous disk dissipates locally the energy it extracts from the shear, whether in a steady state or not. This may not be the case in a turbulent disk where, if the turbulent cascade is not efficient, part of the energy may be advected with the flow. As we have not addressed the energetics of viscous disks above, we will not go into the details of the discussion here. These can be found in Balbus & Papaloizou (1999), who showed that in disks subject to MHD turbulence the energy extracted from the shear is indeed dissipated locally (through the turbulent cascade) whether the disk is evolving or not. Note that this is in general not the case when the turbulence is due to self–gravitating instabilities. In that case indeed, part of the energy is transported by waves through the disk.

Most theoretical protostellar disk models have relied on the $`\alpha `$ parametrization. We have commented that MHD turbulence does lend itself to this prescription in thin disks. However, because MHD instabilities develop only in an adequately ionized fluid, they may not operate everywhere in protostellar disks (Gammie 1996). Therefore it is likely that the parameter $`\alpha `$ is not constant through these disks. It may even be that only parts of these disks can be described using this $`\alpha `$ prescription. However, we may still learn about disks from these models in the same way as we learned about stars from simple polytropic models. Therefore, we now describe these models.

### 5.3 Vertical structure

In a thin disk, the radial gradients of temperature and pressure are much smaller than the vertical gradients. The vertical and radial structures therefore decouple, and can be calculated separately. Here we present the calculation of the vertical structure, i.e. we determine pressure, temperature and density as a function of height $`z`$. Once this is done, given an initial surface density distribution, the radial structure can be computed using the diffusion equation (58). The section below is based on Papaloizou & Terquem (1999), which the reader is referred to for more details.

#### 5.3.1 Basic equations

The equations describing the disk vertical structure are the equation of vertical hydrostatic equilibrium:

$$\frac{1}{\rho }\frac{P}{z}=\frac{GM_{}z}{\left(r^2+z^2\right)^{3/2}},$$
(64)
and the energy equation, which states that the rate of energy removal by radiation is locally balanced by the rate of energy production by viscous dissipation:

$$\frac{}{z}=\frac{9}{4}\rho \nu \mathrm{\Omega }^2,$$
(65)
where $``$ is the radiative flux of energy through a surface of constant $`z`$ which is given by:

$$=\frac{16\sigma T^3}{3\kappa \rho }\frac{T}{z}.$$
(66)
Here $`G`$ is the gravitational constant, $`M_{}`$ is the central mass, $`P`$ is the pressure, $`T`$ is the temperature, $`\kappa `$ is the opacity, which in general depends on both $`\rho `$ and $`T`$, and $`\sigma `$ is the Stefan–Boltzmann constant.

To close the system of equations, we relate $`P`$, $`\rho `$ and $`T`$ through the equation of state of an ideal gas:

$$P=\frac{\rho kT}{\mu m_H},$$
(67)
where $`k`$ is the Boltzmann constant, $`\mu `$ is the mean molecular weight and $`m_H`$ is the mass of the hydrogen atom. If we limit our calculations to temperatures lower than 4,000 K (typical of protostellar disks), then, at the densities of interest, hydrogen is in molecular form. Since the main component of protostellar disks is hydrogen, we then have $`\mu =2`$. We note that for the temperatures we consider transport of energy by convection can be neglected.

We denote the isothermal sound speed by $`c_s`$ ($`c_s^2=P/\rho `$). We adopt the $`\alpha `$–parametrization of Shakura & Sunyaev (1973), so that the kinematic viscosity is written $`\nu =\alpha c_s^2/\mathrm{\Omega }`$ (equations  and ). In general, $`\alpha `$ is a function of both $`r`$ and $`z`$. With this formalism, equation (65) becomes:

$$\frac{}{z}=\frac{9}{4}\alpha \mathrm{\Omega }P.$$
(68)
Note that, in a thin disk, $`zr`$. Then, if we set $`|P/z|P/H`$ and $`zH`$ in equation (64), we get $`c_s\mathrm{\Omega }H`$, where $`\mathrm{\Omega }`$ is the Keplerian velocity ($`\mathrm{\Omega }^2=GM_{}/r^3`$) and we have used $`c_s^2=P/\rho `$. The thin disk approximation then requires $`c_sr\mathrm{\Omega }`$, i.e. the angular velocity to be highly supersonic.

#### 5.3.2 Boundary conditions

We have to solve three first order ordinary differential equations for the three variables $``$, $`P`$ (or equivalently $`\rho `$), and $`T`$ as a function of $`z`$ at a given radius $`r.`$ Accordingly, we need three boundary conditions at each $`r`$. We shall denote with a subscript $`s`$ values at the disk surface.

The first boundary condition is obtained by integrating equation (65) over $`z`$ between $`H`$ and $`H`$. Since by symmetry $`(z=0)=0`$, this gives:

$$_s=\frac{3}{8\pi }\dot{M}_{st}\mathrm{\Omega }^2,$$
(69)
where we have defined $`\dot{M}_{st}3\pi \nu \mathrm{\Sigma }`$, with $`<\nu >`$ being the vertically averaged viscosity. If the disk were in a steady state, $`\dot{M}_{st}`$ would not vary with $`r`$ and would be the constant accretion rate through the disk (eq.  and  with $`R_ir`$). In general however, this quantity does depend on $`r`$.

Another boundary condition is obtained by integrating equation (64) over $`z`$ between $`H`$ and infinity. A detailed derivation of this condition is presented in the Appendix A of Papaloizou & Terquem (1999). Here we just give the result:

$$P_s=\frac{\mathrm{\Omega }^2H\tau _{ab}}{\kappa _s},$$
(70)
where $`\tau _{ab}`$ is the optical depth above the disk. This condition is familiar in stellar structure, where $`\mathrm{\Omega }^2H`$ would be replaced by the acceleration of gravity at the stellar surface. Since we have defined the disk surface such that the atmosphere above the disk is isothermal, we have to take $`\tau _{ab}1`$. Providing this is satisfied, the results do not depend on the value of $`\tau _{ab}`$ we choose.

A third and final boundary condition is given by the expression of the surface temperature (see Appendix A of Papaloizou & Terquem 1999 for a detailed derivation of this expression):

$$2\sigma \left(T_s^4T_b^4\right)\frac{9\alpha kT_s\mathrm{\Omega }}{8\mu m_H\kappa _s}\frac{3}{8\pi }\dot{M}_{st}\mathrm{\Omega }^2=0.$$
(71)
Here the disk is assumed immersed in a medium with background temperature $`T_b`$. The surface opacity $`\kappa _s`$ in general depends on both $`T_s`$ and $`\rho _s`$ and we have used $`c_s^2=kT/(\mu m_H)`$. The boundary condition (71) is the same as that used by Levermore & Pomraning (1981) in the Eddington approximation (their eq.  with $`\gamma =1/2`$). In the simple case when $`T_b=0`$ and the surface dissipation term involving $`\alpha `$ is set to zero, with $`\dot{M}_{st}`$ being retained, it simply relates the disk surface temperature to the emergent radiation flux.

#### 5.3.3  Disk models

At a given radius $`r`$ and for a given value of the parameters $`\dot{M}_{st}`$ and $`\alpha `$ (which uniquely determine the disk model), the disk vertical structure is obtained by solving equations (64), (66) and (68) with the boundary conditions (69), (70) and (71).

Numerical results are published, e.g., by Papaloizou & Terquem (1999) (see also Lin & Papaloizou 1980, 1985; Bell & Lin 1994; Bell et al. 1997), in which the opacity is taken from Bell & Lin (1994). This has contributions from dust grains, molecules, atoms and ions. It is written in the form $`\kappa =\kappa _i\rho ^aT^b`$ where $`\kappa _i`$, $`a`$ and $`b`$ vary with temperature. For a specified $`\dot{M}_{st},`$ they calculate the value of $`H,`$ the vertical height of the disk surface, iteratively. Starting from an estimated value of $`H,`$ after satisfying the surface boundary conditions, the equations are integrated down to the mid–plane $`z=0.`$ The condition that $`=0`$ at $`z=0`$ is not in general satisfied. An iteration procedure is then used to adjust value of $`H`$ until $`=0`$ at $`z=0`$ to a specified accuracy.

## Acknowledgments

It is a pleasure to thank Steven Balbus and John Papaloizou for sharing with me their knowledge and thoughts about accretion disk theory, and for their comments on an early draft of this paper which led to much improvements.

## References

Acheson D. J. 1990, Elementary Fluid Dynamics, Clarendon Press, Oxford.

Batchelor G. K. 1967, An introduction to fluid dynamics, Cambridge University Press.

Blaes O. M., Balbus S. A. 1994, ApJ, 421, 163
Local shear instabilities in weakly ionized, weakly magnetized disks.

Balbus S. A., Hawley J. F. 1991, ApJ, 376, 214
A powerful local shear instability in weakly magnetized disks. I – Linear analysis.

Balbus S. A., Hawley J. F. 1992a, ApJ, 392, 662
Is the Oort A–value a universal growth rate limit for accretion disk shear instabilities?

Balbus S. A., Hawley J. F. 1992b, ApJ, 400, 610
A Powerful Local Shear Instability in Weakly Magnetized Disks. IV – Nonaxisymmetric Perturbations.

Balbus S. A., Hawley J. F. 1998, Rev. Mod. Phys., 70, 1
Instability, turbulence, and enhanced transport in accretion disks.

Balbus S. A., Papaloizou J. C. B. 1999, ApJ, 521, 650
On the Dynamical Foundations of alpha Disks.

Balbus S. A., Terquem C. 2001, ApJ, 552, 235
Linear Analysis of the Hall Effect in Protostellar Disks.

Bell K. R., Cassen P. M., Klahr H. H., Henning Th. 1997, ApJ, 486, 372
The Structure and Appearance of Protostellar Accretion Disks: Limits on Disk Flaring.

Bell K. R., Lin D. N. C. 1994, ApJ, 427, 987
Using FU Orionis outbursts to constrain self–regulated protostellar disk models.

Blandford R. D., Payne D. G. 1982, MNRAS, 199, 883
Hydromagnetic flows from accretion discs and the production of radio jets.

Brandenburg A., Nordlund A., Stein R. F., Torkelsson U. 1995, ApJ, 446, 741
Dynamo–generated Turbulence and Large–Scale Magnetic Fields in a Keplerian Shear Flow.

Cabot W., 1996, ApJ, 465, 874
Numerical Simulations of Circumstellar Disk Convection.

Fleming T. P., Stone J. M., Hawley J. F. 2000, ApJ, 530, 464
The Effect of Resistivity on the Nonlinear Stage of the Magnetorotational Instability in Accretion Disks.

Foglizzo T., Tagger M. 1995, A&A, 301, 293
The Parker–shearing instability in azimuthally magnetized discs.

Fromang S., Terquem C., Balbus S. A., 2001, MNRAS, submitted
The ionization fraction in $`\alpha `$–models of protoplanetary disks.

Gammie C. F. 1996, ApJ, 457, 355
Layered Accretion in T Tauri Disks.

Goldreich P., Lynden–Bell D. 1965, MNRAS, 130, 125
Spiral arms as sheared gravitational instabilities.

Hawley J. F., Balbus S. A., Stone J. M. 2001, ApJ, in press (astro–ph/0103522)
A Magnetohydrodynamic Nonradiative Accretion Flow in Three Dimensions.

Hawley J. F., Gammie C. F., Balbus S. A. 1995, ApJ, 440, 742
Local Three–dimensional Magnetohydrodynamic Simulations of Accretion Disks.

Heemskerk M. H. M., Papaloizou J. C. B., Savonije G. J. 1992, A&A, 260, 161
Non–linear development of m=1 instabilities in a self–gravitating gaseous disc.

Huang K. 1987, Statistical Mechanics, 2nd edition, Wiley.

Jackson J. D. 1975, Classical Electrodynamics, 2nd edition, Wiley.

Jin L., 1996, ApJ, 457, 798
Damping of the Shear Instability in Magnetized Disks by Ohmic Diffusion.

Landau L. D., Lifchitz E. M. 1987, Fluid Mechanics, 2nd edition, Pergamon, Oxford.

Laughlin G., Bodenheimer P. 1994, ApJ, 436, 335
Nonaxisymmetric evolution in protostellar disks.

Levermore C. D., Pomraning G. C 1981, ApJ, 248, 321
A flux–limited diffusion theory.

Lin D. N. C., Papaloizou J. C. B. 1980, MNRAS, 191, 37
On the structure and evolution of the primordial solar nebula.

Lin D. N. C., Papaloizou J. C. B. 1985, in Protostars and planets II, eds D. C. Black & M. S. Mathews (Tucson: University of Arizona Press)
On the dynamical origin of the solar system.

Lynden–Bell D., Kalnajs A. J. 1972, MNRAS, 157, 1
On the generating mechanism of spiral structure.

Lynden–Bell D., Pringle J. E. 1974, MNRAS, 168, 60
The evolution of viscous discs and the origin of the nebular variables.

Ogilvie G. I., Pringle J. E. 1996, MNRAS, 279, 152
The non–axisymmetric instability of a cylindrical shear flow containing an azimuthal magnetic field.

Papaloizou J. C. B. , Lin D. N. C. 1995, ARA&A, 33, 505
Theory Of Accretion Disks. I – Angular Momentum Transport Processes.

Papaloizou J. C., Savonije G. J. 1991, MNRAS, 248, 35
Instabilities in self–gravitating gaseous discs.

Papaloizou J. C. B., Terquem C. 1997, MNRAS, 287, 771
On the stability of an accretion disc containing a toroidal magnetic field: the effect of resistivity.

Papaloizou J. C. B., Terquem C. 1999, ApJ, 521, 823
Critical Protoplanetary Core Masses in Protoplanetary Disks and the Formation of Short–Period Giant Planets.

Pickett B. K., Cassen P., Durisen R. H., Link R. 1998, ApJ, 504, 468
The Effects of Thermal Energetics on Three–dimensional Hydrodynamic Instabilities in Massive Protostellar Disks.

Pringle J. E., 1981, ARA&A, 19, 137
Accretion discs in astrophysics.

Reipurth B., Baly J., Devine D. 1997, AJ, 114, 2708
Giant Herbig–Haro Flows.

Ryu D., Goodman J. 1992, ApJ, 388, 438
Convective instability in differentially rotating disks.

Shakura N. I., Sunyaev R. A. 1973, A&A, 24, 337
Black holes in binary systems. Observational appearance.

Shu F. H. 1992, The Physics of Astrophysics. II – Gas Dynamics, University Science Books.

Stapelfeldt K. R., Krist J. E., Menard F., Bouvier J., Padgett D. L., Burrows C. J. 1998, ApJ, 502, L65
An Edge–On Circumstellar Disk in the Young Binary System HK Tauri.

Stone J. M., Balbus S. A. 1996, ApJ, 464, 364
Angular Momentum Transport in Accretion Disks via Convection.

Sturrock P. A. 1994, Plasma Physics, Cambridge University Press.

Tennekes H., Lumley J. L. 1972, A first course in turbulence, The MIT press.

Terquem C. E. J. M. L. J., in The Formation of Binary Stars, eds. R.D. Mathieu and H. Zinnecker (ASP Conference Series), in press (astro–ph/0008514)
Tidally–induced angular momentum transport in disks.

Terquem C., Papaloizou J. C. B. 1996, MNRAS, 279, 767
On the stability of an accretion disc containing a toroidal magnetic field.

von Weizsäcker C. F. 1948, Z. Naturforsch, 3a, 524

Wardle M. 1999, MNRAS, 307, 849
The Balbus–Hawley instability in weakly ionized discs.


format: text

text:
# Transverse voltage in high-Tc superconductors in zero magnetic fields

## 1 Introduction

High Tc superconductors represent an interesting group of substances for both theoretical and experimental investigations. A lot of new phenomena has been observed for the first time on these materials. For practical applications the most important properties of interest are those belonging to the group of dynamic properties, i.e. physical effects in electric field. From the point of view of fundamental research the well characterized single crystalline samples are asked for. But one should have in mind that polycrystalline ceramic materials rather than single crystals will be the main components of superconducting devices. It is thus necessary to study the properties of such polycrystalline systems as well. In this paper we deal with transport properties of highly textured Bi- based (2223) polycrystalline bulk superconductors in zero magnetic field. The temperature dependence of longitudinal and transversal resistivity has been studied in close vicinity of the transition temperature.

## 2 Experimental

Two different methods of transport measurement have been used. For the sample in the form of thin stripe standard six point contact method was used with contacts for the measurement of longitudinal and transverse voltage, respectively. The misalignment of transverse contacts was corrected for by measuring the transverse voltage in the regime where no Hall voltage should appear i.e. well above transition temperature where the sample was in the normal state.

Several samples were measured by the van der Pauw method . In this case the samples were either in the form of a square or disk with the contacts equally spaced on the perimeter of the sample. The thickness of the samples was in all cases bellow 250 $`\mu `$m. The dependence of the observed nonzero transverse voltage on the current was also measured by both methods. The reason for using standard six point method as a supplement to the van der Pauw method was to verify if the non-homogeneous current distribution in the latter method is not responsible for the observed effect. The results for all type of the samples and methods were qualitatively the same. Results of the van der Pauw measurement are shown in Figs. 1 and 2 for the samples of disk shape of 15mm diameter with thickness of the disk 120$`\mu `$m . From this graph one can clearly see that in the close vicinity of the Tc nonzero transverse voltage appears which is absent both well bellow and above Tc. One can also see the non symmetric shape of the transverse voltage peak. As to the current dependence of the resistance peak its height increases with the increasing current passed through the sample. The current not exceeding 50 mA was used in all measurement. For such currents we did not observe any heating effects. This conclusion has been confirmed by two facts: no change of the measured voltage was detected upon varying the equilibration time ( i.e. the period between switching-on the current and measuring the voltage) and no shift was seen in the temperature where longitudinal voltage vents effectively to zero value.

## 3 Discussion

Glazman  proposed a model for the explanation of the observed effect. Magnetic field produced by the current going through sample can penetrate into the sample in the form of vortices of different sign. This sign is determined by different direction of magnetic field on the opposite edges of the sample resulting in vortices penetration from one side and antivortices penetration from the opposite one. Vortices and antivortices move in opposite direction under the influence of Lorentz force and can annihilate if the attractive interaction between them overcomes the Lorentz force. This means that the path of vortex and antivortex will be distorted and transverse voltage appears according to the Josephson relation. With increasing current the Lorentz force is stronger and the probability of annihilation decreases i.e. for high enough current the trajectories of vortices and antivortices are not influenced by their interaction. They can therefore move perpendicularly to the edges until they reach the opposite side of the sample.

According to this theory the transverse voltage value should increase for low transport current and again decrease for high enough current. We have changed the transport current by two orders of magnitude in our experiments and did not observe any decrease of transverse voltage. Similar effect was observed by Francavilla et al.  on thin sputtered YBaCuO films of unknown thickness using currents up to 180 mA. Moreover the calculated magnetic field on the surface of our sample is in the range of units of $`\mu `$T. At such small magnetic fields the concentration of vortices and antivortices is probably too small to initiate the observable voltage.

In the following we will propose another explanation of the above mentioned effect. Since Bi based materials are strongly anisotropic their electrical behavior should reflect the two-dimensional-like nature. Near the percolation transition when the phase of the wave function becomes coherent along a privileged path in the sample and weak links among grains become irrelevant , the behavior should be similar to that occurring at Kosterlitz-Thoules temperature TKT . Bellow TKT vortex-antivortex pairs are thermally activated . At TKT these pairs spontaneously dissociate into free vortices. These free vortices and antivortices can move under the influence of external fields and cause in this way dissipation of the energy above TKT. There are two possibilities how to determine TKT. First one, most simple, is to determine this temperature as that at which the resistance of the sample reaches zero value at zero magnetic field. It is commonly believed that the beginning of the dissipation is connected with the creation of free vortices and with their movement . According to the Coulomb gas model  there are no free vortices present below TKT and no flux- flow resistance . However free vortices should , according to this model, be generated below TKT provided a finite current is imposed across the superconducting sample. The Coulomb gas model prediction for the flux-flow resistance generated below TKT in this way is equivalent to a nonlinear I-V characteristic of the form $`V=I^{\alpha (T)}`$. Precisely at TKT the prediction for the exponent is $`\alpha =3`$ which was confirmed by many of experiments. The exponent $`\alpha `$(T) for our measurement is shown in Fig. 3. One can see that the value $`\alpha =3`$ corresponds to the temperature 104 K i.e. TKT = 104 K. This value roughly coincides with the temperature where resistance in zero magnetic field effectively reaches zero (see Fig. 1). It is shown in Fig. 2 that the nonzero transverse resistance starts at the same temperature. Set of free vortices and antivortices created above TKT is forced to move by Lorentz and Magnus forces and dissipation can start and according to the Josephson relation the nonzero transverse resistance appears. This suggestion is also supported by sudden increase of the resistivity from the side of low temperatures.

## Acknowledgments

This work has been partially supported by GACR projects Nos. 104/99/1440 and 106/99/1441 and by GAASCR under contract No. 11116.

## Figure captions

Fig.1: Temperature dependence of longitudinal resistance Rxx for different current in zero external magnetic field

Fig.2: Temperature dependence of transverse resistance Rxy for different current in zero external magnetic field

Fig.3: Temperature dependence of the exponent $`\alpha `$ for nonlinear I-V characteristics ($`V=I^\alpha `$) in zero external magnetic field