galactic black hole candidates ( bhcs ) are known to exhibit transitions between different x - ray spectral states .
typically , a transient outburst begins in the low / hard state at a low luminosity .
the x - ray spectrum in the low / hard state is roughly described by a hard power law with a high energy cutoff at @xmath6 .
as the luminosity increases , these systems undergo a transition to the high / soft state ( so - called a hard - to - soft transition ) .
the x - ray spectrum in the high / soft state is dominated by the disk emission of characteristic temperature @xmath7 .
recently , two distinct types of hard - to - soft transitions , the bright / slow transition and the dark / fast transition , are reported ( e.g. , * ? ? ?
* ; * ? ? ?
the bright / slow transition occurs at @xmath8 and takes more than @xmath9 days .
the system undergoes a transition from the low / hard state to the high / soft state via the `` bright / hard '' state and the very high / steep power law ( vh / spl ) state during the bright / slow transition .
the x - ray spectrum in the bright / hard state is described by a hard power - law with a ( moderately ) high energy cutoff at @xmath10 and the luminosity is `` brighter '' than that in the low / hard state .
the dark / fast transition occurs at less than @xmath11 and takes less than @xmath12 days .
the system immediately switches from the low / hard state to the high / soft state during the dark / fast transition .
here , @xmath13 is the eddington luminosity , @xmath14 is the black hole mass , and @xmath15 is the electron scattering opacity .
@xcite analyzed the results of _ rxte _ observations of the bhc gx 339 - 4 in the rising phases of the transient outbursts ( this object showed the bright / slow transition in the 2002/2003 outburst and the dark / fast transition in the 2004 outburst ) .
they found that the cutoff energy strongly anti - correlates with the luminosity and decreases from @xmath16 to @xmath17 in the bright / hard state , while the cutoff energy is roughly constant at @xmath16 in the low / hard state .
this suggests that the electron temperature of an accretion disk emitting hard x - rays decreases as the luminosity increases in the bright / hard state .
furthermore , the bright / hard state has been observed from @xmath18 up to @xmath8 .
they concluded that such anti - correlation is explained by the scenario that the heating rate from protons to electrons via the coulomb collision balances the radiative cooling rate of inverse compton scattering .
the main purpose of this paper is to present a model explaining the bright / hard state .    in the conventional theory of accretion disks , the concept of phenomenological @xmath19-viscosity
is introduced . in this framework ,
the @xmath20-component of the stress tensor , which appears in the angular momentum equation and the viscous heating term , is assumed to be proportional to the gas pressure , @xmath21 ( we ignore the radiation pressure in this paper because we focus on optically thin disks ) .
here @xmath22 is the viscosity parameter introduced by @xcite .
the magnetic field can be an origin of the @xmath19-viscosity because the maxwell stress generated by the magneto - rotational instability ( mri ) efficiently transports angular momentum in accretion disks ( e.g. , * ? ? ? * ) and the dissipation of magnetic energy contributes the disk heating ( e.g. , * ? ? ? * ) .
optically thin , hot accretion disks have been studied to explain hard x - rays from bhcs .
@xcite proposed that hard x - rays from cyg x-1 are produced in an inner optically thin hot disk .
@xcite studied the structure and stability of optically thin hot accretion disks .
we have to consider two - temperature plasma in such disks because the electron temperature is expected to become lower than the ion temperature in such a low density , high temperature region .
the energy equations for ions and electrons are written in the form @xmath23 @xmath24 where @xmath25 is the viscous heating rate , @xmath26 is the energy transfer rate from ions to electrons via coulomb collisions , @xmath27 is the radiative cooling rate , and @xmath28 is the fraction of heating to electrons .
the left hand - sides represent the heat advection terms for ions ( @xmath29 ) and electrons ( @xmath30 ) .
we note that early works on optically thin , hot , two - temperature accretion disks assumed that the viscous heating acts primarily on ions ( @xmath31 ) .
@xcite and ( * ? ? ?
* hereafter sle ) constructed a model for optically thin two - temperature accretion disks . in the sle solutions ,
the dissipated energy is transferred from ions to electrons via coulomb collisions ( @xmath32 ) and radiated away by electrons ( @xmath33 ) .
although the electron temperature is high enough to explain the x - ray spectrum in the low / hard state , the sle solutions are thermally unstable in the framework of the @xmath19-prescription of viscosity .
@xcite pointed out the importance of heat advection in hot , magnetized accretion flows , and obtained steady solutions of optically thin disks .
such geometrically thick , optically thin , advection - dominated accretion flows ( adafs ) or radiatively inefficient accretion flows ( riafs ) have been studied extensively by ( * ? ? ?
* @xcite ) and @xcite . in the adaf / riaf solutions , a substantial fraction of the dissipated energy
is stored in the gas as entropy and advected into the central object ( @xmath34 ) .
only a small fraction of the dissipated energy is transferred to electrons and radiated away .
the adaf / riaf solutions are thermally stable .
@xcite found that the maximum mass accretion rate of the adaf / riaf solutions is @xmath35 which corresponds to @xmath36 , where @xmath37 is the eddington mass accretion rate .
@xcite showed that the electron temperature in the adaf / riaf solutions weakly anti - correlates with the luminosity and decreases to @xmath38 .
these features are consistent with the facts that the energy cutoff weakly anti - correlates with the luminosity around @xmath16 in the low / hard state , and that these systems undergo a transition from the low / hard state to other x - ray spectral states ( i.e. , the high / soft state during the dark / fast transition , and the bright / hard state during the bright / slow transition ) about at this maximum luminosity of the adaf / riaf solutions
. however , the adaf / riaf solutions can not explain the strong anti - correlation in the range of @xmath39 and @xmath40 observed in the bright / hard state .
the heat advection works as an effective cooling in the adaf / riaf solutions .
@xcite presented a luminous hot accretion flow ( lhaf ) in which the heat advection for ions works as an effective heating . above the maximum mass accretion rate of the adaf / riaf solutions ,
the heat advection overwhelms the viscous heating and balances the energy transfer from ions to electrons ( @xmath41 ) .
the lhafs are thermally unstable .
however , @xcite concluded that the thermal instability will have no effect on the dynamics of the lhafs because the accretion timescale is shorter than the timescale of growth of the local perturbation at such high mass accretion rate .
the lhaf solutions also can not explain the bright / hard state because the electron temperature is high and roughly constant at @xmath42 .    in these models mentioned above ,
magnetic fields are not considered explicitly , and the ratio of the gas pressure to the magnetic pressure is assumed to be constant ( typically , @xmath43 ) .
@xcite suggested that an accretion disk evolves toward two types of disks , a high-@xmath1 disk and a low-@xmath1 disk , by carrying out two - dimensional magnetohydrodynamic ( mhd ) simulations of the buoyant escape of the magnetic flux owing to the parker instability @xcite .
in the high-@xmath1 disk , the magnetic flux escapes from the disk owing to the parker instability and @xmath1 inside the disk is maintained at a high value .
global three - dimensional mhd simulations of optically thin , radiatively inefficient accretion disks also indicate that the amplification of magnetic fields becomes saturated when @xmath44 in a quasi - steady state except in the plunging region very close to the black hole ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
on the other hand , once a disk is dominated by the magnetic pressure , it can stay in the low-@xmath1 state because the strong magnetic tension suppresses the growth of the parker instability .
@xcite demonstrated that an optically thin , radiatively inefficient , hot , high-@xmath1 ( adaf / riaf - like ) disk undergoes transition to an optically thin , radiatively efficient , cool , low-@xmath1 disk except in the plunging region when the mass accretion rate exceeds the threshold for the onset of a cooling instability . during this transition ,
the magnetic flux @xmath45 is almost conserved at each radius because the cooling timescale is shorter than that of the buoyant escape of the magnetic flux , where @xmath46 is the mean azimuthal magnetic field and @xmath47 is the half thickness of the disk . in this way , the magnetic pressure becomes dominant and supports the disk as the gas pressure decreases owing to the cooling instability .
eventually , the disk stays in a quasi - steady , cool , low-@xmath1 state .
because the mri is not yet stabilized in this quasi - steady state , the magnetic field still remains turbulent and dominated by the azimuthal component .
as a result , the heating owing to the dissipation of the turbulent magnetic field balances the radiative cooling .    @xcite performed local three - dimensional mhd simulations of strongly magnetized , vertically stratified accretion disks in a keplerian potential .
they showed that strongly magnetized state is maintained near the equatorial plane because the buoyantly escaping magnetic flux is replenished by stretching of a radial field .
the mri feeds off both vertical and azimuthal fields and drives turbulence .
the maxwell and reynolds stresses generated by the turbulence become significant .
therefore , they indicated that highly magnetized disks are astrophysically viable .
we note that such low-@xmath1 disks are quite different from magnetically dominated accretion flows ( mdafs ; * ? ? ?
* ) observed in the plunging region of optically thin accretion disks in global mhd and general relativistic mhd simulations ( e.g. , * ? ? ?
* ) in terms of their energy balance and configuration of magnetic fields .
outside the plunging region , the magnetic field becomes turbulent and dominated by the azimuthal component because the growth timescale of the mri is shorter than the inflow timescale . as a result
, the released gravitational energy is efficiently converted into the thermal energy via the dissipation of the turbulent magnetic field .
on the other hand in the plunging region , the ratio of the timescales is reversed because the inflow velocity increases with decreasing the radius in such disks . therefore ,
magnetic field lines are stretched out in the radial direction before turbulence is generated by the mri and is dissipated . as a result , a substantial fraction of the gravitational energy is converted into the radial infall kinetic energy . since the heating owing to the dissipation of turbulent magnetic fields becomes inefficient , the gas pressure becomes low , and the flow becomes magnetically dominated .
although both of the low-@xmath1 disk and the mdaf are cool and magnetic pressure dominant , they are essentially different .
we focus on the low-@xmath1 disk in this paper .
@xcite constructed a steady model of optically thin , one - temperature accretion disks incorporating magnetic fields on the basis of these results of three - dimensional mhd simulations of accretion disks .
@xcite extended it to the optically thick regime .
they assumed that the @xmath20-component of the stress tensor is proportional to the total pressure . in order to complete the set of basic equations
, they specified the advection rate of the azimuthal magnetic flux instead of @xmath1 .
they found a new thermally stable solution , a low-@xmath1 solution , which can explain the results by @xcite . in the low-@xmath1 solutions , the magnetic heating enhanced by the strong magnetic pressure balances the radiative cooling . the disk temperature is lower than that in the adaf / riaf solutions and strongly anti - correlates with the mass accretion rate
they also found that the low-@xmath1 solutions exist above the maximum mass accretion rate of the adaf / riaf solutions .
therefore , they concluded that the optically thin low-@xmath1 disk can qualitatively explain the bright / hard state .
however , they considered one - temperature plasma and bremsstrahlung emission as a radiative cooling mechanism .
it is expected that the electron temperature becomes lower than the ion temperature and that synchrotron cooling and/or inverse compton scattering become effective in such disks .    in this paper , we extend the model of optically thin , one - temperature disk to that of optically thin , two - temperature disks .
we consider synchrotron emission and inverse compton scattering as a cooling mechanism in addition to bremsstrahlung emission .
we obtained the thermal equilibrium curves and found that the optically thin low-@xmath1 solutions can quantitatively explain the bright / hard state .
the basic equations are presented in section [ basic_eq ] . in section
[ results ] , we present the thermal equilibrium curves . section [ discussion ] is devoted to a discussion .
we summarize the paper in section [ summary ] .
we extended the basic equations for one - dimensional steady , optically thin , two - temperature black hole accretion flows ( e.g. , * ? ? ? * ) incorporating magnetic fields .
we adopt cylindrical coordinates @xmath48 .
general relativistic effects are simulated using the pseudo - newtonian potential @xmath49 @xcite , where @xmath50 is the gravitational constant , @xmath14 is the black hole mass ( we assume @xmath51 in this paper ) , @xmath52 , and @xmath53 is the schwarzschild radius . for simplicity , the gas is assumed to consist of protons ( ions ) and electrons .
the number density of ions and electrons are equal by charge neutrality , @xmath54 .
we start with the resistive mhd equations @xmath55
@xmath56 =    - \rho \nabla \psi - \nabla p_{\rm gas } + \frac{\mbox{\boldmath{$j$ } } \times \mbox{\boldmath{$b$}}}{c } ~ ,   \label{eq : vec_mom}\end{aligned}\ ] ] @xmath57 - \left (   \mbox{\boldmath{$v$ } } \cdot \nabla \right ) p_{\rm i } \nonumber \\   = \left(1 -   \delta_{\rm heat } \right ) q^{+ } - q^{\rm ie } ~ ,   \label{eq : vec_ene_i}\end{aligned}\ ] ] @xmath58 - \left (   \mbox{\boldmath{$v$ } } \cdot \nabla \right ) p_{\rm e } \nonumber
\\   = \delta_{\rm   heat } q^{+ } +   q^{\rm ie } - q_{\rm rad}^{- } ~ ,   \label{eq : vec_ene_e}\end{aligned}\ ] ] @xmath59 where @xmath60 is the density , @xmath61 and @xmath62 are the ion and electron densities , @xmath63 and @xmath64 are the ion and electron masses , @xmath65 is the velocity , @xmath66 is the magnetic field , @xmath67 is the current density , @xmath68 is the gas pressure , @xmath69 and @xmath70 are the ion and electron gas pressure , @xmath71 and @xmath72 are the ion and electron temperature , @xmath73 is the boltzmann constant , @xmath74 and @xmath75 are the internal energy of ions and electrons .
here , @xmath76 and @xmath77 are the specific heat ratio for ions and electrons . in the energy equations for ions ( [ eq : vec_ene_i ] ) and electrons ( [ eq : vec_ene_e ] ) , @xmath25 is the heating rate , @xmath27 is the radiative cooling rate , and @xmath26 is the energy transfer rate from ions to electrons via coulomb collisions . here , @xmath28 represents the fraction of heating to electrons . in the induction equation ( [ eq : vec_ind ] ) , @xmath78 is the magnetic diffusivity , where @xmath79 is the electric conductivity .
three - dimensional global mhd and local radiation - mhd simulations of black hole accretion disks showed that magnetic fields inside the disk are turbulent and dominated by the azimuthal component in a quasi - steady state ( e.g. , * ? ? ?
* ; * ? ? ?
on the basis of results of the simulations , we decomposed the magnetic fields into the mean fields @xmath80 and fluctuating fields @xmath81 and also decomposed the velocity into the mean velocity @xmath82 and the fluctuating velocity @xmath83 .
we assume that the fluctuating components vanish when azimuthally averaged , @xmath84 , and that the radial and vertical components of the magnetic fields are negligible compared with that of the azimuthal component , @xmath85 , @xmath86 . here
@xmath87 denotes the azimuthal average .
let us derive the azimuthally averaged equations .
we assume that the disk is in a steady state and in hydrostatic balance in the vertical direction . by
azimuthally averaging equations ( [ eq : vec_con ] ) - ( [ eq : vec_ind ] ) and ignoring the second order terms of @xmath88 , @xmath89 , and @xmath90 , we obtain @xmath91 @xmath92 @xmath93 + \frac{\partial}{\partial z }   \left ( \frac{\langle b_{\varphi } b_{z } \rangle}{4 \pi}\right)~,\end{aligned}\ ] ] @xmath94 @xmath95   + \frac{v_{\varpi}}{\varpi } \left ( \rho_{\rm i } \epsilon_{\rm i } +   p_{\rm i } \right ) +   \frac{\partial}{\partial z } \left [   \left ( \rho_{\rm i } \epsilon_{\rm i } + p_{\rm i } \right ) v_{z}\right ]   \nonumber \\    - v_{\varpi } \frac{\partial}{\partial \varpi } p_{\rm i } - v_{z }    \frac{\partial}{\partial z } p_{\rm i } = \left ( 1 - \delta_{\rm heat } \right ) q^{+ } -    q^{\rm ie } ~,\end{aligned}\ ] ] @xmath96   + \frac{v_{\varpi}}{\varpi } \left ( \rho_{\rm e } \epsilon_{\rm e } +   p_{\rm e } \right ) +   \frac{\partial}{\partial z } \left [   \left ( \rho_{\rm e } \epsilon_{\rm e } + p_{\rm e } \right ) v_{z}\right ]   \nonumber \\    - v_{\varpi } \frac{\partial}{\partial \varpi } p_{\rm e } - v_{z }    \frac{\partial}{\partial z } p_{\rm e } = \delta_{\rm heat } q^{+ } + q^{\rm ie } -    q^{-}_{\rm rad } ~,\end{aligned}\ ] ] @xmath97 -\frac{\partial}{\partial \varpi } \left [ v_{\varpi } \langle   b_{\varphi } \rangle \right ]   \nonumber \\    + \ { \nabla \times \langle
\delta    \mbox{\boldmath{$v$ } } \times \delta \mbox{\boldmath{$b$ } } \rangle   \}_{\varphi } - \ { \eta_{\rm m } \nabla \times \left ( \nabla \times    \mbox{\boldmath{$\bar{b}$ } } \right ) \}_{\varphi } ~,\end{aligned}\ ] ] where @xmath98 is the total pressure and @xmath99 is the azimuthally averaged magnetic pressure .
the third and fourth terms on the right - hand side of equation ( [ eq : ind_phi ] ) represent the dynamo term and the magnetic diffusion term which we approximate later on the basis of the results of the numerical simulations .
we assume that the radial velocity @xmath100 , the specific angular momentum @xmath101 , and @xmath102 are independent of @xmath103 , and that the disks are isothermal in the vertical direction . under these assumptions ,
the surface density @xmath104 , the vertically integrated total pressure @xmath105 , and the half thickness of the disk @xmath47 are defined as @xmath106 @xmath107 @xmath108 where @xmath109 is the keplerian angular velocity . here
the subscript @xmath110 refers to quantities in the equatorial plane . using the equation of state for the ideal gas ,
the vertically integrated total pressure is expressed as @xmath111    now we integrate the other basic equations in the vertical direction .
we obtain @xmath112 @xmath113 @xmath114 @xmath115 @xmath116 @xmath117 d\varpi dz + \mbox{const.}\end{aligned}\ ] ] where @xmath118 is the mass accretion rate , @xmath119 is the keplerian angular momentum and @xmath120 is the specific angular momentum swallowed by the black hole . in the energy equations , @xmath121 and @xmath122
are the vertically integrated heat advection terms for ions and electrons , @xmath123 , @xmath124 , and @xmath125 are the vertically integrated heating rate , radiative cooling rate , and energy transfer rate from ions to electrons via coulomb collisions . in equation ( [ eq : ind_int ] ) ,
@xmath126 is the radial advection rate of the azimuthal magnetic flux ( hereafter we call it the magnetic flux advection rate ) .    in this paper
, we assume that @xmath127 instead of equation ( [ eq : mom_pi_int ] ) because we focus on local thermal equilibrium solutions .
global mhd simulations of radiatively inefficient , accretion flows ( e.g. , * ? ? ? * ; * ? ? ?
* ) showed that the ratio of the azimuthally averaged maxwell stress to the sum of the azimuthally averaged gas pressure and magnetic pressure is nearly constant ( @xmath128 ) except in the plunging region . following the simulation results , we assume that the azimuthally averaged @xmath129-component of the maxwell stress inside a disk is proportional to the total ( gas and magnetic ) pressure @xmath130 integrating in the vertical direction
, we obtain @xmath131 this is one of the key assumptions in this paper . when the magnetic pressure is high , the stress level can be high even though the gas pressure is low .
we can rewrite this relation in terms of the kinematic viscosity , @xmath132 , as @xmath133 where @xmath134 @xmath135 is the sound speed , and @xmath136 is the alfvn speed .
we complete the set of basic equations by specifying the radial distribution of the magnetic flux advection rate .
if we perform the integration in the second term of the induction equation ( [ eq : ind_int ] ) , we obtain @xmath137 + \mbox{const . } \nonumber\end{aligned}\ ] ] where @xmath138 is the mean azimuthal magnetic field in the equatorial plane .
according to the result of the global three - dimensional mhd simulation by @xcite , the magnetic flux advection rate at a radius is roughly unchanged before and after the transition from the adaf / riaf - like disk to the low-@xmath1 disk .
hence , we adopt the magnetic flux advection rate as the parameter in order to complete the set of the basic equations .
the magnetic flux advection rate depends on various mechanisms such as the escape of magnetic fluxes due to the magnetic buoyancy , the regeneration of azimuthal magnetic fields by the shear motion , the generation of magnetic turbulence through the mri , dissipation of magnetic fields due to the magnetic diffusivity , and magnetic reconnection .
if the sum of the dynamo term and the magnetic diffusion term is zero in the whole region , the magnetic flux advection rate is constant in the radial direction .
the global three - dimensional mhd simulation performed by @xcite indicated that the magnetic advection rate increases with decreasing radius , specifically , @xmath139 , in the quasi steady state as a result of magnetic dynamo and diffusivity processes . because it is hard to compute the dynamo term and the magnetic diffusion term explicitly from the local quantities , we parameterize the dependence of @xmath126 on @xmath140 by introducing a parameter , @xmath141 , as follows .
@xmath142 where @xmath143 is the magnetic flux advection rate at the outer boundary @xmath144 .
when @xmath145 , the magnetic flux advection rate is constant in the radial direction .
when @xmath146 , the magnetic flux advection rate increases with decreasing radius ( see also figure 1 in * ? ? ?
we determine the parameter @xmath143 by imposing the outer boundary condition , @xmath147 \left(\varpi_{\rm out } / r_{\rm s } \right)^{-1}$ ] and @xmath148 at @xmath149 .
this leads that @xmath150 . in the local model presented in this paper ,
the value of @xmath141 just means the amount of the magnetic flux advection rate at a radius ( e.g , @xmath151 , @xmath152 , @xmath153 for @xmath145 , @xmath154 , @xmath155 , respectively ) .    equation ( [ eq : phidot ] )
is the second key assumption in this paper .
specifying the magnetic flux advection rate enables the magnetic pressure to increase when the disk temperature decreases .
by contrast , if we specified the plasma @xmath1 at each radius instead of the magnetic flux advection rate , the decrease in the temperature results in a decrease in magnetic pressure . this is inconsistent with the results of three - dimensional mhd simulations ( e.g. , * ? ? ?
we address the similarity between the concepts of the mass accretion rate and the magnetic flux advection rate in order to facilitate understanding of the concept of the magnetic flux advection rate .
the right - hand side of equation ( [ eq : con_int ] ) represents the mass flux crossing at @xmath140 per unit time and we have denoted it by @xmath118 .
if there is no mass loss or gain ( e.g , due to inflows and/or outflows ) , @xmath118 is constant in the radial direction .
otherwise @xmath118 is a function of @xmath140 ( e.g. , @xmath156 , @xmath157 is a parameter ) .
the right - hand side of equation ( [ eq : ind_int_ignore ] ) represents the azimuthal magnetic flux crossing at @xmath140 per unit time and we have denoted it by @xmath126 .
if there is no loss or gain of azimuthal magnetic fields ( e.g , due to the dynamo and/or the magnetic diffusion ) , @xmath126 is constant in the radial direction .
otherwise @xmath126 is a function of @xmath140 and we prescribed it as equation ( [ eq : phidot ] ) .
that is , when @xmath158 , the magnetic flux is conserved in the radial direction , and when @xmath159 ( or @xmath160 ) , the magnetic flux increases ( or decreases ) with a decreasing radius ( see figure 1(b ) in * ? ? ?
we note that the azimuthal magnetic flux inside a disk can increase when the azimuthal flux of opposite polarity buoyantly escapes from the disk ( e.g. , * ? ? ?
* ) .        in the conventional theory ,
the viscous heating was expressed as @xmath161 and assumed to heat primarily ions ( @xmath162 ) , where @xmath163 is the @xmath129-component of the total stress and @xmath164 is the angular velocity .
the results of three - dimensional mhd simulations indicate that the dissipation of the magnetic energy dominates the total dissipative heating rate throughout a disk and is expressed as @xmath165 ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
hereafter , we refer to it as the magnetic heating rate . following these simulation results ,
we employ magnetic heating as the heating mechanism inside a disk , and set the vertically integrated heating rate as follows :    @xmath166 dz = - \alpha w_{\rm tot } \varpi
\frac{d \omega}{d \varpi } ~ , \end{aligned}\ ] ]    where we have used equation ( [ eq : al ] )
. we note that if the magnetic pressure is high , the heating rate can also be high even when the gas pressure is low .
the eventual expression of the heating term is not at all unusual ( e.g. , * ? ? ?
* ; * ? ? ?
* ) except that we consider much lower values of @xmath1 .
we will discuss in section [ discussion_main ] a lower limit of @xmath1 below which the mri is stabilized ( e.g. , * ? ? ?
* ) so that this expression is no longer valid .
we note that the magnetic heating does not always heat primarily ions .
@xcite suggested that @xmath167 is required to fit the spectrum of sgr @xmath168 with riaf models .
@xcite carried out local shearing box simulations of the nonlinear evolution of the mri in a collisionless plasma considering pressure anisotropy and showed that @xmath169 ( note that the definition of @xmath170 is different from that of @xmath171 in * ? ? ?
the exact value we choose for this parameter is not so important , in particular , for the low-@xmath1 solutions .
we find that the low-@xmath1 solutions presented in this paper practically unchanged for any value of @xmath172 .
therefore , we adopt @xmath173 as a fiducial value
.      if the ion temperature is higher than the electron temperature , coulomb collisions transfer energy from ions to electrons .
the energy transfer rate from ions to electrons per unit volume via coulomb collisions is given by @xcite @xmath174 , \end{aligned}\ ] ] where @xmath175 is the thomson scattering cross section and @xmath176 is the coulomb logarithm ( roughly @xmath177 ) .
@xmath178 are modified bessel function of the second kind of the order @xmath179 , respectively .
the dimensionless electron and ion temperatures are defined by @xmath180    for technical reason , we use the following formula which uses no special functions , and is accurate to within a factor of 2 when @xmath181 @xcite @xmath182 integrating in the vertical direction , we obtain @xmath183 & ( \theta_{\rm i } > 0.2)\\ \\
\displaystyle    ( k t_{\rm e } - k   t_{\rm i } ) \frac{(2 \pi)^{1/2 } +    ( \theta_{\rm e } +     \theta_{\rm i})^{1/2}}{(\theta_{\rm e } + \theta_{\rm i } ) } &    ( \theta_{\rm i } < 0.2 ) \nonumber \end{array}\right .   ~ .\end{aligned}\ ] ]    we note that @xmath184 in almost all solutions presented in this paper .
we assume that the radiative cooling occurs through electrons and consider bremsstrahlung , synchrotron , and compton cooling by bremsstrahlung and synchrotron photons as cooling processes .
the vertically integrated radiative cooling rate is expressed as @xmath185    following @xcite ( see also * ? ? ?
* ; * ? ? ?
* ) , bremsstrahlung cooling rate per unit volume is @xmath186~,\end{aligned}\ ] ] where the subscripts @xmath187 and @xmath188 denote the electron - ion and electron - electron bremsstrahlung cooling rates , @xmath189 is fine - structure constant , the function @xmath190 and @xmath191 have the approximate form @xmath192 & ( \theta_{\rm e } > 1 ) \\ \\
\displaystyle 4 \left ( \frac{2 \theta_{\rm e}}{\pi^3}\right)^{1/2 } \left [ 1    + 1.781 \theta_{\rm e}^{1.34}\right ] & ( \theta_{\rm e } < 1 ) \end{array } \right .
~,\end{aligned}\ ] ] @xmath193 @xmath194 and @xmath195 is euler s number .
integrating in the vertical direction , we obtain @xmath196 ~.\end{aligned}\ ] ]    the synchrotron emissivity of a relativistic maxwellian distribution of electrons in the optically thin limit is given by @xcite : @xmath197 where @xmath198 with the fitting function @xmath199 given by @xcite @xmath200 integrating in the vertical direction , we obtain @xmath201 we assumed that the emission below a critical frequency , @xmath202 , is completely self - absorbed so that the emissivity can be approximated by the blackbody emission from the surface of the disk .
following @xcite , we estimate @xmath202 as the frequency at which the synchrotron emission from the region @xmath203 is equal to the blackbody emission ( in the rayleigh - jeans limit ) from the upper and lower surfaces of the region .
this condition gives the equation @xmath204 we obtain @xmath202 by solving this equation numerically . integrating over frequency , we obtain @xmath205 ,   \end{aligned}\ ] ] where the parameters @xmath206 , @xmath207 , @xmath208 , and @xmath209 are defined as @xmath210 and @xmath211 is the incomplete gamma function @xmath212    we adopted the prescription for the compton energy enhancement factor @xmath213 described by @xcite , which is defined to be the average change in energy of a photon between injection and escape : @xmath214 \nonumber \\   & \equiv & 1 + \eta_1 -\eta_2 \left ( \frac{x}{\theta_{\rm                  e}}\right)^{\eta_3 } ~,\end{aligned}\ ] ] where @xmath215 here , @xmath216 is the probability that an escaping photon is scattered , @xmath217 is the mean amplification factor in the energy of a scattered photon when the scattering electrons have a maxwellian velocity distribution of temperature @xmath218 .
following @xcite , the vertically integrated compton cooling by bremsstrahlung and synchrotron photons are given respectively by @xmath219\right\ } ~,\end{aligned}\ ] ] @xmath220 ~ , \end{aligned}\ ] ] where @xmath221 .
the vertically integrated heat advection terms for ions and electrons are expressed as @xmath222 @xmath223 where @xmath224 @xmath225 \frac{\partial \ln t_{\rm    e}}{\partial \ln \varpi } \nonumber \\
+ \frac{\partial \ln \sigma}{\partial \ln    \varpi } - \frac{\partial \ln h}{\partial \ln \varpi}~ , \end{aligned}\ ] ] are dimensionless quantities of the order of unity ( hereafter we call them the entropy gradient parameter ) , @xmath226 , @xmath227 , respectively .
the positivity ( negativity ) of the entropy gradient parameter means that the heat advection term works as an effective cooling ( heating ) .
the entropy gradient parameter for ions , @xmath228 , has a positive value in adaf / riaf solutions , while it has a negative value in lhaf solutions .
the global three - dimensional mhd simulations ( e.g. , * ? ? ?
* ) and steady , vertically integrated , one - dimensional transonic solutions ( e.g. , * ? ? ?
* ; * ? ? ?
* ) of optically thin black hole accretion flows indicate that @xmath229 .
following these results , we adopt @xmath230 as a fiducial value in this paper .
we note that a value of @xmath228 is not important in sle and low-@xmath1 solutions because the heat advection term is negligible compared to the other terms .
according to @xcite , @xmath231 can have a positive and negative value ( @xmath232 ) .
however , they assumed that no viscous heat goes into electrons ( @xmath233 ) . in the energy equation for electrons , when @xmath172 , the dissipated magnetic energy term is typically greater than or comparable to the heat advection term in the inner region of the disk .
therefore the exact value we choose for this parameter is not so important , in particular , for sle and low-@xmath1 solutions in which @xmath234 .
we find that the results presented in this paper are practically unchanged for any value of @xmath231 between @xmath235 and @xmath236 when @xmath237 .
we show the results for @xmath238 in most part of this paper .
we solved the above basic equations at @xmath239 for given parameters @xmath118 , @xmath19 , @xmath141 , @xmath28 , @xmath228 , and @xmath231 .
we obtained new thermal equilibrium solutions , low-@xmath1 solutions , in addition to the adaf / riaf ( for positive @xmath228 ) , sle , lhaf ( for negative @xmath228 ) solutions in the optically thin regime .
figure [ al05ds3xp05sifoo ] shows the sequences of each thermal equilibrium solution in the @xmath104 versus @xmath240 , @xmath241(thin line ) , @xmath72(thick line ) , and @xmath1 plane .
the disk parameters we adopted are @xmath242 , @xmath230 , @xmath238 , @xmath173 , @xmath243 ( solid ) , @xmath154 ( dashed ) , and @xmath110 ( dotted ) , respectively . here
@xmath244 is the eddington mass accretion rate , @xmath245 is the energy conversion efficiency , and @xmath246 is the electron scattering opacity .
we obtain three types of solutions , adaf / riaf ( for @xmath247 at this radius ) , sle ( for @xmath248 ) , and low-@xmath1 solutions ( for @xmath249 ) .
we find that the low-@xmath1 solutions exist above the maximum mass accretion rate of the adaf / riaf solutions , @xmath250 .
this indicates that the disk initially staying in the adaf / riaf state undergoes transition to the low-@xmath1 state when the mass accretion rate exceeds @xmath251 .
furthermore , the electron temperature in the low-@xmath1 solutions is lower ( @xmath252 ) than that in the adaf / riaf solutions .
the energy balance for ions and electrons is illustrated in figure [ al05ds3xp05sieie ] .
the upper panel shows the ratio of the heat advection to the magnetic heating for ions , @xmath253 ( this quantity is referred to as the advection factor @xmath254 , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the lower panel shows the ratio of the heat advection to the total heating for electrons ( thin line ) and the fraction of the energy transfer via coulomb collisions to the total heating ( thick line ) .
the electrons receive the total amount of the magnetic heating in the low-@xmath1 solutions as well as in the sle solutions even though we introduced the parameter @xmath28 which represents the fraction of heating to electrons .
the fraction @xmath255 of the magnetic heating goes into ions and the fraction @xmath28 of the magnetic heating goes into electrons
. however , almost all the magnetic heating going into ions is transferred to electrons via coulomb collisions in the low-@xmath1 solutions ( @xmath256 ) . eventually , the total amount of the magnetic heating goes into electrons .
this means that the parameter @xmath28 does not appear practically in the energy balance for electrons .
the radiative cooling overwhelms the heat advection in the low-@xmath1 solutions .
therefore , the magnetic heating balances the radiative cooling in the low-@xmath1 solutions ( @xmath257 ) .    the energy balance is essentially the same in both the sle solutions and the low-@xmath1 solutions ( @xmath258 )
the difference is which pressure dominates the magnetic heating . the magnetic heating , which is proportional to the total pressure in our model , is dominated by the gas pressure in the sle solutions and the magnetic pressure in the low-@xmath1 solutions .
next we describe the main cooling mechanism in each solution .
figure [ al05ds3xp05sirr ] shows the vertically integrated bremsstrahlung ( solid ) , bremsstrahlung - compton ( dotted ) , synchrotron ( dashed ) , and synchrotron - compton ( dash - dotted ) cooling rate for @xmath145 ( bottom ) , @xmath154 ( middle ) , and @xmath155 ( top ) . when @xmath145 , the synchrotron - compton cooling is dominant in the low-@xmath1 solutions for lower mass accretion rate while the bremsstrahlung - compton cooling is dominant for high mass accretion rates ( @xmath259 , @xmath260 ) even though @xmath261 .
the synchrotron cooling is relatively ineffective for higher mass accretion rates because of the lower electron temperature . as @xmath141 increases ,
the electron temperature become high because the large magnetic flux enhances not only the synchrotron cooling but also the magnetic heating . as a result ,
the synchrotron and synchrotron - compton cooling become efficient . when @xmath243 , the synchrotron - compton is dominant in whole low-@xmath1 solutions .
now we show that @xmath262 , @xmath263 , @xmath264 on the low-@xmath1 branch ( see also * ? ? ?
* ) , where @xmath265 is the mean temperature .
these relations depend on the dependence of @xmath143 on @xmath118 ( our outer boundary condition leads @xmath266 ) . here
we introduce the parameter @xmath157 , @xmath267 , in order to leave this dependence explicitly .
first , we derive the relations between @xmath118 and @xmath104 .
since @xmath268 on the low-@xmath1 branch , equations ( [ eq : h2 ] ) , ( [ eq : mom_phi_int ] ) , and ( [ eq : al_int ] ) yield @xmath269 and @xmath270
. using equations ( [ eq : con_int ] ) , ( [ eq : ind_int_ignore ] ) , and ( [ eq : phidot ] ) , we find that @xmath271 , @xmath272 , and @xmath273 .
next , we derive the relation between @xmath274 and @xmath104 from the energy balance of the low-@xmath1 solutions ( @xmath257 ) . here
we roughly approximate the radiative cooling rate by @xmath275 for simplicity ( this is the same dependence as non - relativistic bremsstrahlung cooling for single temperature plasma ) .
since @xmath276 , we find that @xmath277 .
therefore @xmath278 .
we find that @xmath279 , @xmath263 , and @xmath280 when @xmath281 .
we investigate the relation between the mass accretion rate and the electron temperature in order to explain the anti - correlation between the luminosity and the cutoff energy observed in the bright / hard state . in the low-@xmath1 solutions
, we expect that the luminosity is roughly proportional to the mass accretion rate since @xmath282 , and that the electron temperature roughly represents the cutoff energy since the inverse - compton scattering is the dominant radiative cooling mechanism .
therefore , the relation between the electron temperature and the mass accretion rate is useful for comparison with the observational data .
figure [ al05ds3xp05temd ] shows the relations between the electron temperature and the mass accretion rate for the same parameters as in figure [ al05ds3xp05sifoo ] .
we find that the electron temperature is typically @xmath283 in the low-@xmath1 solutions while it is @xmath284 in the adaf / riaf solutions . furthermore , the electron temperature in the low-@xmath1 solutions strongly anti - correlates with the mass accretion rate .
since @xmath285 and @xmath286 in the low-@xmath1 solutions , we find that @xmath287 for @xmath281 .
according to @xcite , the mri is stabilized for @xmath288 ( we plot this critical point ( filled circle in figure [ al05ds3xp05temd ] ) at which @xmath289 ) .
therefore , the low-@xmath1 solutions may not exist under the condition that @xmath290 .
we discuss this issue in section [ discussion ] .
we also show the results for @xmath291 in figure [ al20ds3xp05sifoo ] and [ al20ds3xp05temd ] .
the maximum mass accretion rate of the adaf / riaf solution is around @xmath292 .
this maximum mass accretion rate is higher than that for @xmath293 .
here we investigate the dependence of @xmath294 on @xmath19 .
since @xmath295 , @xmath296 and @xmath297 at @xmath298 , the energy equations yield ( i ) @xmath299 , ( ii ) @xmath300 , and ( iii ) @xmath301 q^{+ } \sim q_{\rm rad}^{-}$ ] . here
we assume @xmath28 to be constant for simplicity .
we find from equation ( i ) that @xmath302 .
if @xmath124 can be written in a simple form of @xmath303 like @xmath304 , we find from equations ( ii ) and ( iii ) that @xmath72 at @xmath305 is independent of @xmath19 and @xmath251 is proportional to @xmath306 exactly .
although @xmath124 has a complicated form in our model , our numerical results indicate that the @xmath72 is roughly independent of @xmath19 ( @xmath307 ) . using this value of @xmath72 instead of solving equation ( iii )
, we find from equation ( ii ) that roughly @xmath308 .
this result is consistent with @xmath309 by @xcite .
we investigated the dependence of the results on @xmath28 because @xmath28 is a poorly constrained parameter .
we considered three cases : @xmath310 as an example of model in which ions and electrons receive equal amounts of the dissipated magnetic energy , @xmath311 ( @xmath312 ) as an example of a conventional model in which ions receive a substantial amount of the energy , and @xmath313 as an intermediate example .
we note that @xmath314 lies between @xmath315 and @xmath316 in the solutions presented in this paper .
figure [ al05xp05sifoo ] shows the thermal equilibrium curves for the same parameters as in figure [ al05ds3xp05sifoo ] but for @xmath317 ( solid ) , @xmath236 ( dashed ) , @xmath318 ( dotted ) , and @xmath319 ( dot - dashed ) , respectively .
the energy balance for ions and electrons is illustrated in figure [ al05xp05sieie ] .
the ion temperature decreases as @xmath28 increases because the magnetic heating for ions ( @xmath320 ) which is the only heating source for ions decreases .
@xmath305 also decreases because of the decrease in the magnetic heating for ions .
the electron temperature increases as @xmath28 increases in the adaf / riaf solutions .
by contrast , in the low-@xmath1 solutions , the electron temperature is almost independent of the value of @xmath28 because @xmath321 does not appear practically in the energy equation for electrons ( @xmath322 ) .
therefore , we find that the exact value of @xmath321 is not so important for the low-@xmath1 solutions .
figure [ al05ds3xp05z000sifoo ] shows the thermal equilibrium curves for the same parameters as in figure [ al05ds3xp05sifoo ] but for @xmath323 ( dashed ) .
we also plotted the thermal equilibrium curves for @xmath324 ( solid ) for comparison .
we obtained the lhaf solutions in the high mass accretion rate and high surface density region .
the heat advection works as a heating for ions and balances the energy transfer via the coulomb collisions in the lhaf solutions . above @xmath305 ,
the heat advection term overwhelms the magnetic heating .
hence , the ion temperature becomes higher than that in the adaf / riaf solutions .    in the lhaf solutions ,
the electrons receive a larger amount of heat from ions than that in the adaf / riaf solutions . nonetheless , the electron temperature becomes lower than that in the adaf / riaf solutions ( yet higher than that in the low-@xmath1 solutions ) because the radiative cooling becomes more effective in such higher surface density region .
we find that our results are practically unchanged for any value of @xmath231 between @xmath235 and @xmath236 because the heat advection term is negligible compared to the other terms in energy equation for electrons except when @xmath311 .
furthermore , even when @xmath325 , the other quantities except the electron temperature in the adaf / riaf solutions does not change practically .
we show the thermal equilibrium curves plotted on the @xmath104 - @xmath71 ( thin ) and @xmath72 ( thick ) plane for @xmath311 , @xmath238 ( solid ) , @xmath110 ( dashed ) , @xmath235 ( dotted ) , @xmath243 ( left panel ) , and @xmath110 ( right panel ) in figure [ al05dd3sitt ] .
we obtained thermal equilibrium solutions for an optically thin , two - temperature accretion disk incorporating magnetic fields .
we included bremsstrahlung emission , synchrotron emission , and inverse compton scattering as the radiative cooling mechanisms , and introduced the parameter @xmath28 which represents the fraction of the magnetic heating to electrons .
we prescribed the @xmath129-component of the azimuthally averaged maxwell stress is proportional to the total pressure . in order to complete the set of basic equations
, we specified the radial distribution of the magnetic flux advection rate by introducing a parameter @xmath141 .
we found a branch of low-@xmath1 solutions in addition to the usual adaf / riaf ( for @xmath326 ) , sle , and lhaf ( for @xmath327 ) solutions .
here we remark why we can obtain the low-@xmath1 solutions .
first , we prescribed that the @xmath129-component of the maxwell stress is proportional to the total pressure . therefore , if the magnetic pressure is high , we can obtain the magnetic heating rate balancing the radiative cooling rate even in high surface density and low temperature region .
second , we specified the magnetic flux advection rate , @xmath328 , in order to complete the set of the basic equations . in the conventional theory , @xmath1
is assumed to be constant ( typically , @xmath329 ) , which means that the magnetic pressure is proportional to the gas pressure .
this implies that @xmath330 is just multiplied by a constant ( @xmath331 ) . as a result ,
we can not obtain the sequence of the low-@xmath1 solutions . on the other hand in our model , a decrease in temperature results in an increase in magnetic pressure ( therefore an increase in the magnetic heating ) under the conservation of the magnetic flux in the vertical direction .
this is the reason why we can obtain the sequence of the low-@xmath1 solutions .
we note that the exact values of @xmath228 , @xmath231 , and @xmath28 are not so important in the low-@xmath1 solutions . the magnetic heating enhanced by the high magnetic pressure balances the radiative cooling in the low-@xmath1 solutions .
the heat advection terms including @xmath228 and @xmath332 are negligible for both ions and electrons .
furthermore , @xmath28 does not appear in the energy balance for electrons practically .
let us discuss the lower limit of @xmath1 .
since @xmath280 and @xmath333 under our outer boundary condition ( @xmath281 ) , we find that @xmath334 , that is , @xmath1 decreases as the mass accretion rate increases in the low-@xmath1 solutions .
it has been confirmed by global mhd simulations ( e.g. , * ? ? ?
* ) that the mri is not stabilized and hence the magnetic heating rate is expressed in the form of equation ( [ eq : qmag ] ) , at least , when @xmath335 .
local mhd simulations ( e.g. , * ? ? ?
* ) also indicated that the maxwell and reynolds stresses generated by magnetic turbulence are significant and yield an effective @xmath19-viscosity ( @xmath336 ) in highly magnetized disks ( @xmath329 ) .
however , the expression of the magnetic heating employed in our paper may no longer be valid for much lower values of @xmath1 because strong magnetic fields suppress the growth of the mri .
we implicitly assumed that the dissipation energy of the turbulent magnetic fields generated by the mri is converted into the thermal energy of the disk gas .
therefore , such heating mechanism becomes ineffective if the mri is stabilized .
@xcite studied the evolution of the mri in differentially rotating , magnetized flows beyond the weak - field limit , and showed that the mri is stabilized for toroidal alfvn speeds , @xmath337 , exceeding the geometric mean of the sound speed , @xmath338 , and the rotational speed , @xmath339 , ( i.e. , @xmath288 , or equivalently @xmath340 ) .
our results satisfy this condition when @xmath341 , @xmath342 , @xmath343 for @xmath242 , and @xmath344 , @xmath345 , and @xmath346 for @xmath347 , respectively .
these critical points are denoted by filled circles in figure [ al05ds3xp05temd ] and [ al20ds3xp05temd ] ( see also table [ tbl-1 ] ) .
when @xmath1 falls below this critical value , the low-@xmath1 solutions may not exist because there is no heating source balancing the radiative cooling . as a result ,
the disk may undergo a transition to other states ( e.g , an mdaf - like disk or an optically thick disk ) .
let us discuss the thermal stability in this subsection .
a general criterion concerning the thermal instability of disks can be expressed as ( see * ? ? ?
* ; * ? ? ?
* ) @xmath348_{\sigma } >    0 ~ .\end{aligned}\ ] ] in the low-@xmath1 solutions , we ignore the heat advection term because @xmath349 .
once again we approximate the radiative cooling rate by @xmath350 for simplicity .
we find from @xmath351 and @xmath352 ( see section [ results_main ] ) that the low-@xmath1 solutions do not satisfy the criterion ( [ eq : criterion ] ) , that is , are thermally stable .
we note that the thermal stability is independent of @xmath157 .
we also remark the thermal stability of the other solutions ( adaf / riaf , sle , and lhaf ) in which the gas pressure is dominant ( @xmath353 ) . equations ( [ eq : h2 ] ) and ( [ eq : con_int ] ) yield @xmath354 and @xmath355 . we find from @xmath356 , @xmath357 , @xmath358 that the adaf / riaf solutions ( @xmath359 ) are thermally stable but the sle ( @xmath258 ) and lhaf ( @xmath360 ) solutions are thermally unstable .
the main purpose of this paper is to explain the bright / hard state observed during the bright / slow transition in the rising phases of the transient outbursts of bhcs . in the low / hard state ,
the x - ray spectrum is described by a hard power law with a high energy cutoff at @xmath16 .
when their luminosity exceed @xmath361 , these systems undergo a transition from the low / hard state to the bright / hard state . in the bright / hard state ,
the cutoff energy decreases from @xmath16 to @xmath362 as the luminosity increases from @xmath361 to @xmath363 ( e.g. , * ? ? ?
* ) . beyond the bright / hard state ,
these systems undergo a transition to the high / soft state going through the vh / spl state .
the adaf / riaf solution explains the x - ray spectrum in the low / hard state because the electron temperature is high ( @xmath364 ) .
however , this solution does not exist at the high mass accretion rates and does not show the strong anti - correlation between the electron temperature and the mass accretion rate observed in the bright / hard state .
the low-@xmath1 solution extends to such high mass accretion rates .
therefore , the disk initially staying in the adaf / riaf state undergoes transition to the low-@xmath1 state when the mass accretion rate exceeds @xmath294 . on the low-@xmath1 branches ,
the electron temperature is low ( @xmath365 @xmath366 ) and strongly anti - correlates with the mass accretion rate .
these features are consistent with the anti - correlation between the luminosity and the energy cutoff observed in the bright / hard state .
therefore , the low-@xmath1 solution can explain the bright / hard state .    in section [ discussion_main ]
, we have discussed the lower limit of @xmath1 below which the mri is stabilized therefore the low-@xmath1 solutions may not exist . since @xmath367 in the low-@xmath1 solutions
, this means that the low-@xmath1 solution has a maximum mass accretion rate , @xmath368 .
the mass accretion rate and the electron temperature at the lower limit of @xmath1 are depicted by filled circles in figure [ al05ds3xp05temd ] and [ al20ds3xp05temd ] ( see also table [ tbl-1 ] ) .
we found that when @xmath369 ( or equivalently , @xmath370 ) @xmath371 .
this indicates that the disk staying in the adaf / riaf state undergoes transition to the low-@xmath1 disk when the mass accretion rate exceeds @xmath372 , after that , the disk undergoes transition to the optically thick disk when the mass accretion rate exceeds @xmath373 . in this case , the hard - to - soft transition occurs at @xmath373 ( not at @xmath251 ) .
this can correspond to the bright / slow transition .    here
we also remark on the dark / fast transition during which the system undergoes transition from the low / hard state to the high / soft state without going through the bright / hard state and the vh / spl state .
we found that when @xmath374 ( or equivalently , @xmath375 ) @xmath376 , that is , there is no optically thin , thermally stable solution above @xmath372 .
therefore , the disk staying in the adaf / riaf state immediately undergoes transition to an optically thick disk without through the low-@xmath1 disk .
this can correspond to the dark / fast transition .
accordingly , we conclude that the bright / slow transition occurs when @xmath126 has a large value and the dark / fast transition occurs when @xmath126 has a small value .
we have obtained the low-@xmath1 solutions for optically thin , two - temperature accretion disks incorporating the mean azimuthal magnetic fields , and concluded that the low-@xmath1 solutions explain the bright / hard state observed during the bright / slow transition of bhcs .
we assumed that the energy transfer from ions to electrons occurs through coulomb collisions , and considered bremsstrahlung emission , synchrotron emission , and inverse compton scattering as the radiative cooling processes .
we prescribed that the maxwell stress is proportional to the total ( gas and magnetic ) pressure . in order to complete the set of basic equations
, we specified the radial distribution of the magnetic flux advection rate .
accordingly , a decrease in temperature can result in an increase in magnetic pressure under the conservation of the magnetic flux in the vertical direction . in the low-@xmath1 solutions ,
the magnetic heating is enhanced by the high magnetic pressure .
the fraction @xmath377 of the magnetic heating goes into ions and is transferred to electrons via coulomb collisions .
the fraction @xmath28 of the magnetic heating goes into electrons . eventually , the total amount of the magnetic heating goes into electrons and balances the radiative cooling ( compton cooling by bremsstrahlung and/or synchrotron photons ) . the electron temperature is lower than that in the adaf / riaf solutions ( @xmath378 ) and strongly anti - correlates with the mass accretion rate in the low-@xmath1 solutions .
these features are consistent with the x - ray spectrum observed in the bright / hard state .    according to @xcite ,
the mri is stabilized for @xmath379 .
this indicates that the low-@xmath1 solutions disappear below this critical point .
when the magnetic flux advection rate is high ( @xmath369 ) , the critical mass accretion rate of the low-@xmath1 solutions , @xmath373 , is above the maximum mass accretion rate of the adaf / riaf solutions , @xmath294 .
therefore , the disk initially staying the adaf / riaf state undergoes transition to the low-@xmath1 state when the mass accretion rate exceeds @xmath251 , after that , the disk undergoes transition to the optically thick disk when the mass accretion rate exceeds @xmath373 .
this corresponds to the bright / slow transition . on the other hand , when the magnetic flux advection rate is low ( @xmath380 ) , @xmath381 is below @xmath294 .
therefore , the disk initially staying the adaf / riaf state immediately undergoes transition to the optically thick disk when the mass accretion rate exceeds @xmath251 .
this can correspond to the dark / fast transition .
we are grateful to r. narayan for valuable discussion and closely examining my draft .
this work is supported by the grant - in - aid for science research of the ministry of education , culture , sports , science and technology ( rm : 20340040 ) and grant - in - aid for jsps fellows ( 20.1842 ) .
llccc @xmath382 & @xmath155 & @xmath383 & @xmath384 & @xmath385 + & @xmath154 & @xmath386 & @xmath387 & @xmath388 + & @xmath389 & @xmath390 & @xmath391 & @xmath392 + @xmath318 & @xmath155 & @xmath393 & @xmath394 & @xmath395 + & @xmath154 & @xmath396 & @xmath397 & @xmath398 + & @xmath110 & @xmath399 & @xmath400 & @xmath392 +