Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: magnetic fields are ubiquitous in diverse astrophysical settings on various scales ranging from proto - stellar discs , young stellar objects ( ysos ) , microquasars , quasars , galaxies and grandiose astrophysical jets to clusters of galaxies . in some cases , magnetic fields are relatively weak so that dynamical processes are almost unaffected by their presence . however , there do exist numerous cases where magnetic fields are necessary and important for both dynamics and diagnostics , especially in spiral galaxies and accretion disc systems ( e.g. , sofue et al . 1986 ; beck et al . 1996 ; balbus & hawley 1998 ; tagger & pellat 1999 ; widrow 2002 ) . in general , it is challenging to model magnetic fields realistically because of their complexities ( sofue et al . 1986 ; kronberg 1994 ; caunt & tagger 2001 ; widrow 2002 ) . in the problems involving protostar formation and disc galaxies , shu & li ( 1997 ) presumed the so - called ` isopedic ' magnetic field configuration where the mass - to - magnetic flux ratio in a razor thin disc remains constant . using a singular isothermal disc ( sid ) model ( mestel 1963 ) that is isopedically magnetized , shu et al . ( 2000 ) studied stationary coplanar perturbations and proposed that these global modes are bifurcations to either secularly or dynamically unstable configurations . along a separate yet complementary line , lou ( 2002 ) studied an azimuthally magnetized singular isothermal disc ( msid ) and derived two different global stationary magnetohydrodynamic ( mhd ) perturbation modes . based on the msid model , lou ( 2002 ) explored the manifestation of interlaced optical and magnetic field spiral arms in the outer portion of a disc with a nearly flat rotation curve such as the case of the nearby spiral galaxy ngc 6946 ( beck & hoernes 1996 ; fan & lou 1996 ; lou & fan 1998a , 2002 ; frick et al . 2000 , 2001 ; lou et al . 2002 ) .    in contexts of spiral galaxies , it is natural and important to consider a composite disc system consisting of gravitationally coupled gaseous and stellar discs . this is because various physical processes on different scales occur in the gaseous disc yet large - scale gas dynamics and environment are significantly affected by large - scale structures in the stellar disc ( lou & fan 1998b ) . to construct global coplanar perturbation structures in a systematic manner , we started with a composite sid system of two fluid discs without magnetic field ( lou & shen 2003 ) to generalize the work of shu et al . ( 2000 ) for a single sid . in terms of the global axisymmetric stability for such a composite sid system , we provided a more straightforward @xmath12criterion ( shen & lou 2003 ) in contrast to the axisymmetric stability criteria of elmegreen ( 1995 ) and jog ( 1996 ) . to be more general than a background sid profile , we recently constructed global stationary perturbation patterns in a composite system of two fluid scale - free discs without magnetic field ( shen & lou 2004a , b ) to further generalize the work of syer & tremaine ( 1996 ) . with proper adaptations , these global stationary perturbation solutions may be utilized to model large - scale structures spiral galaxies . while a scale - free disc has no characteristic scales , one may impose boundary conditions to describe a modified scale - free disc , e.g. , by cutting a central hole in a disc ( e.g. , zang 1976 ; evans & read 1998a , b ) . with such boundary conditions , discrete eigenfunctions for perturbations may be constructed as growing normal modes . on the other hand , by specifying a phase relation for a postulated reflection of spiral waves from the origin @xmath13 , goodman & evans ( 1999 ) could also define discrete normal modes even for an unmodified gaseous sid . shu et al . ( 2000 ) speculated that the swing amplification process ( e.g. , goldreich & lynden - bell 1965 ; julian & toomre 1966 ; toomre 1977 ; fan & lou 1997 ) across corotation allows a continuum of normal modes ( lynden - bell & lemos 1993 ) while proper ` boundary conditions ' may select from this continuum a discrete spectrum of unstable normal modes ; they also suggested that zero - frequency ( stationary ) disturbances do signal the onset of instability . on the basis of the work of shu et al . ( 2000 ) for a single isopedically magnetized sid and of lou & shen ( 2003 ) for a composite sid system without magnetic field , lou & wu ( 2004 ) analyzed coplanar mhd perturbation structures in a composite sid system with one of the disc being isopedically magnetized ( shu & li 1997 ) . as x - ray emitting hot gases can freely flow along open magnetic field lines anchored at the gaseous disc , global stationary solutions constructed by lou & wu ( 2004 ) set an important preliminary stage for modelling spiral galactic mhd winds . in particular , if activities of star formation , starbursts and supernovae etc . are mainly responsible for producing x - ray emitting hot gases streaming out of the galactic plane from both sides , we would expect stronger winds from the circumnuclear starburst region ( e.g. lou et al . 2001 ) and along spiral arms than those winds from the remainder of the galactic disc plane . both sid and msid models offer valuable physical insight for disc dynamics and belong to a wider family of scale - free discs ( syer & tremaine 1996 ) where physical variables in the axisymmetric background equilibrium scale as powers of radius @xmath14 ( e.g. , the rotation curve @xmath15 ) . it turns out that coplanar perturbations in such scale - free discs in a proper range of @xmath8 can be treated globally and analytically ( lemos et al . 1991 ; syer & tremaine 1996 ) without invoking the usual wkbj or tight - winding approximation to solve the poisson equation for density wave perturbations ( lin & shu 1964 ; binney & tremaine 1987 ; bertin & lin 1996 ) . in other words , it becomes feasible to include azimuthal magnetic fields in a scale - free disc and to construct global coplanar mhd perturbation configurations that are stationary in an inertial frame of reference . such analytical solutions , still idealized and limited , are making important further steps and are extremely valuable for bench - marking numerical codes and for initializing numerical mhd simulations .    for an azimuthally magnetized scale - free disc with an infinitesimal thickness , the axisymmetric background equilibrium without radial flow is characterized by following features : a surface mass density @xmath0 , a rotation curve @xmath2 , an azimuthal magnetic field @xmath1 and a barotropic equation of state @xmath16 where @xmath5 is the vertically integrated gas pressure with constant @xmath17 for a warm disc and @xmath18 as the barotropic index . to maintain a radial force balance at all radii , there exists an explicit relationship among @xmath7 , @xmath8 , @xmath9 and @xmath6 such that only one parameter is independent . as will be seen in section 2 , it is convenient to use @xmath8 as this independent parameter ( syer & tremaine 1996 ; shen & lou 2004b ) . in order to satisfy relevant constraints , the prescribed range of @xmath8 turns out to be @xmath10 that includes the previous special case of @xmath19 ( i.e. , the msid case ) studied by lou ( 2002 ) , lou & zou ( 2004a ) and lou & wu ( 2004 ) . in section 2 , we introduce coplanar mhd perturbations and present the linearized perturbation equations . we construct and analyze global stationary aligned and logarithmic spiral perturbation configurations in section 3 and as an example of illustration , we discuss specifically the case of @xmath11 with the background azimuthal magnetic field being force - free . in section 4 , we derive phase relationships between perturbation enhancements of surface mass density and magnetic field for both aligned and logarithmic spiral cases , and discuss the problem of angular momentum conservation . finally , we summarize our results in section 5 . for the convenience of reference , relevant technical details are collected in the appendices . in the magnetofluid approximation , the disc is taken to be razor - thin ( i.e. , we use vertically integrated magnetofluid equations and neglect vertical derivatives of physical variables ) and large - scale stationary aligned and spiral coplanar disturbances develop in a background mhd rotational equilibrium of axisymmetry ( syer & tremaine 1996 ; lou 2002 ; lou & zou 2004a ; lou & wu 2004 ) . the background magnetic field is taken to be azimuthal to avoid the magnetic field winding dilemma ( lou & fan 1998a ) . for the sake of simplicity at this stage , non - ideal effects such as viscosity , resistivity and thermal diffusion etc . are ignored for large - scale perturbations . using cylindrical coordinates @xmath20 , we have the following two - dimensional nonlinear mhd equations for a razor - thin disc geometry : the mass conservation equation @xmath21 where @xmath22 is the surface mass density , @xmath23 is the radial bulk flow velocity , @xmath24 is the specific angular momentum in the vertical @xmath25-direction and @xmath26 is the azimuthal linear velocity ; the radial component of the momentum equation @xmath27\ , \end{split}\ ] ] where @xmath5 is the vertically integrated gas pressure ( sometimes referred to as two - dimensional pressure ) , @xmath28 is the total gravitational potential including that of an axisymmetric dark matter halo in contexts of disc galaxies , and @xmath29 and @xmath30 are the azimuthal and radial components of the coplanar magnetic field @xmath31 , respectively ; the azimuthal component of the momentum equation @xmath32\ ; \end{split}\ ] ] the poisson integral equation @xmath33^{1/2}}\ , \ ] ] where @xmath34 is defined as the ratio of the potential arising from the disc to that arising from the entire system including a massive dark matter halo that is presumed not to respond to coplanar perturbations in the disc plane ( syer & tremaine 1996 ; shu & li 1997 ; shu et al . 2000 ; lou & shen 2003 ; shen & lou 2004a , b ; lou & zou 2004a ) with @xmath35 for a full disc and @xmath36 for partial discs ; the divergence - free condition of magnetic field @xmath37 @xmath38 the radial component of the magnetic induction equation @xmath39 the azimuthal component of the magnetic field induction equation @xmath40 and the barotropic equation of state @xmath41 where constants @xmath17 and @xmath18 . among the three equations ( [ divergfree ] ) , ( [ radiaminduct ] ) and ( [ azimuminduct ] ) , we can freely choose two independent ones . by barotropic equation of state ( [ barotropic ] ) , the sound speed @xmath42 is defined by @xmath43 where subscript @xmath44 indicates the background equilibrium . an isothermal sound speed corresponds to a barotropic index @xmath45 . we now proceed to derive properties of an axisymmetric rotational mhd equilibrium characterized by physical variables associated with a subscript @xmath44 . in our notations , such a background equilibrium has the following power - law scalings : surface mass density of @xmath0 , a radial velocity profile @xmath46 and @xmath47 , a purely azimuthal magnetic field with @xmath48 and @xmath49 . substitutions of these power - law scalings into equations @xmath50 lead to the following condition for the radial force balance , namely @xmath51 where the alfvn speed @xmath52 in a thin disc is defined by @xmath53 and the poisson integral gives @xmath54 with the numerical factor @xmath55 explicitly defined by @xmath56 and @xmath57 being the standard gamma function . this expression ( [ p0 ] ) can also be included in a more general form of @xmath58 as in expression ( [ pm ] ) derived later with the limiting result of @xmath59 as @xmath60 .    to satisfy the radial force balance equation ( [ sfeqn ] ) at all radii ( i.e. , scale - free ) would necessarily require ) in general and for the special case of @xmath11 , we have @xmath61 for the equilibrium magnetic field to be force - free ( e.g. , low & lou 1990 ) . the non - force - free msid case studied by lou ( 2002 ) has @xmath62 , @xmath19 , @xmath63 and @xmath45 . ] @xmath64 which immediately leads to the explicit expressions of indices @xmath7 , @xmath9 and @xmath6 all in terms of @xmath8 @xmath65    in this formulation , @xmath66 remains finite only if @xmath67 which implies @xmath68 . for the total gravity force arising from the equilibrium surface mass densities to be finite , a larger @xmath69range @xmath70 is allowed . the @xmath69range is also constrained by @xmath10 where the first inequality is required for @xmath18 and the second inequality is imposed such that the central point mass will not diverge . by these considerations , the plausible range of @xmath8 falls in @xmath71 ( see syer & tremaine 1996 and shen & lou 2004a , b ) .    according to equations ( [ sfeqn ] ) and ( [ pararelation ] ) , we can explicitly write all physical variables of the background equilibrium state in power - law scalings of @xmath14 in terms of parameter @xmath8 , namely @xmath72 } { 2\pi g(2\beta{\cal p}_0)r^{1 + 2\beta}}\ , \ \ \qquad a^2=\frac{a^2}{(1 + 2\beta)r^{2\beta}}\ , \ \ \qquad c_a^2\equiv q^2a^2=\frac{q^2a^2}{(1 + 2\beta)r^{2\beta}}\ , \\ & v_0=\frac{ad}{r^{\beta}}\ , \ \ \qquad j_0\equiv rv_0=\frac{ad}{r^{\beta-1}}\ , \ \ \qquad \omega\equiv\frac{j_0}{r^2}=\frac{ad}{r^{1+\beta}}\ , \ \ \qquad\kappa\equiv [ ( 2\omega / r)d(r^2\omega)/dr]^{1/2}=[2(1-\beta)]^{1/2}\omega\ . \end{split}\ ] ] by introducing a reference radius @xmath73 , the constant @xmath74 is actually related to a scaled sound speed @xmath75 , the constant @xmath76 $ ] is a scaled rotational mach number , the constant @xmath77 is the ratio of the alfvn speed @xmath52 to the sound speed @xmath42 , @xmath78 and @xmath79 are the equilibrium angular speed and the epicyclic frequency , respectively ( lou 2002 ; shen & lou 2004a , b ) . the first expression of equation ( [ propequilib ] ) puts an upper limit on the magnetic field strength , namely , @xmath80 . alternatively , once a magnetic field parameter @xmath81 is chosen , the scaled rotational mach number @xmath82 must meet the following physical requirement @xmath83 for a positive equilibrium surface mass density @xmath84 . note that for @xmath85 , inequality ( [ restrt ] ) is automatically satisfied without actually restricting @xmath81 and the relevant physical requirement is simply @xmath86 . not all stationary solutions of @xmath87 can satisfy condition ( [ restrt ] ) . this then implies limits on allowed magnetic field strength for a known background rotational mhd equilibrium . it is also important to note that for the axisymmetric background rotational mhd equilibrium under consideration , the magnetic force arising from the azimuthal magnetic field is radially inward when @xmath88 , radially outward when @xmath89 and is zero when @xmath11 [ i.e. , the azimuthal magnetic field is force - free ( lou & fan 1998a ) and the equilibrium density distribution reduces to that of the single unmagnetized disc case ( lemos et al . 1991 ; syer & tremaine 1996 ) ] . these different possibilities can be physically understood in terms of the competition between the magnetic pressure and tension forces . in order to see this more specifically for the background , we can write @xmath90 $ ] referring to equation ( [ radialmomeneqn ] ) . the first term in the square brackets stands for the magnetic pressure force and the second term stands for the magnetic tension force for a purely azimuthal magnetic field @xmath91 in the cylindrical geometry . as @xmath92 scales @xmath93 , we have @xmath94 . therefore for @xmath95 or equivalently @xmath88 , the magnetic tension force dominates and the total magnetic force is radially inward ; for @xmath96 or equivalently @xmath89 , the magnetic pressure force dominates and the total magnetic force is radially outward . the two magnetic forces balance each other for @xmath61 ( or @xmath11 ) . on the basis of the full mhd equations @xmath97 , we readily derive the linearized coplanar mhd perturbation equations as @xmath98 @xmath99 -\frac{1}{\sigma_0}\int\frac{dzb_{\theta}}{4\pi r } \frac{\partial(rb_0)}{\partial r}\ , \end{split}\ ] ] @xmath100 @xmath101^{1/2}},\ ] ] @xmath102 @xmath103 @xmath104 where we use subscript @xmath105 to denote associations with small disturbances in physical variables and @xmath106 stands for the coplanar magnetic field perturbation . as the background rotational mhd equilibrium is stationary and axisymmetric , these perturbed physical variables can be decomposed in terms of fourier harmonics with the periodic dependence @xmath107 where @xmath108 is the angular frequency and @xmath109 is an integer for azimuthal variations . more specifically , we write @xmath110 where @xmath111 , @xmath112 , @xmath113 , @xmath114 , @xmath115 and @xmath116 are all radial variations of the corresponding perturbed physical variables and can be complex in general for possible radial oscillations . without loss of generality , we take @xmath117 in our analysis .    with fourier decomposition ( [ fourier ] ) , it is then straightforward to cast coplanar mhd perturbation equations @xmath118 into the following forms of @xmath119 @xmath120\ , \end{split}\ ] ] where @xmath121 is a short - hand notation , and @xmath122 @xmath123 @xmath124 @xmath125 @xmath126 where equation ( [ perturb07 ] ) can be derived by combining equations ( [ perturb05 ] ) and ( [ perturb06 ] ) . we now rearrange the time - dependent coplanar mhd perturbation equations by taking @xmath127 from equations ( [ perturb01])@xmath128([perturb07 ] ) ; and the special case of @xmath129 will be analyzed in details at the end of this subsection .    from equations ( [ perturb05 ] ) and ( [ perturb06 ] ) , we readily obtain @xmath130 a substitution of equation ( [ zr ] ) into the radial and azimuthal components of the momentum equation ( [ perturb02 ] ) and ( [ perturb03 ] ) yields @xmath131 \bigg(\frac{u}{\omega - m\omega}\bigg)\ , \end{split}\ ] ] and @xmath132 respectively .    in order to construct global solutions without the wkbj approximation , we are mainly interested in stationary configurations with zero pattern speed @xmath133 , also referred to as neutral @xmath134 modes ( e.g. , syer & tremaine 1996 ; shu et al . 2000 ; lou 2002 ; shen & lou 2004a , b ; lou & zou 2004a ) . by setting @xmath134 and taking @xmath127 , equations ( [ perturb01 ] ) , ( [ zr ] ) , ( [ radial ] ) and ( [ azimuthal ] ) can be readily reduced to @xmath135 @xmath136 @xmath137\frac{\hbox{i}u}{\omega}\ , \end{split}\ ] ] @xmath138    equations ( [ stationary1 ] ) through ( [ stationary4 ] ) together with equation ( [ perturb04 ] ) are the basic mhd perturbation equations for constructing stationary non - axisymmetric @xmath127 configurations of both aligned and unaligned logarithmic spiral cases . note that aligned and spiral stationary global solutions both involve propagations of fast and slow mhd density waves ( lou 2002 ) . we next present the basic coplanar mhd perturbation equations for stationary axisymmetric @xmath129 configurations with or without radial propagations . starting from equations ( [ perturb01 ] ) , ( [ perturb06 ] ) , ( [ perturb07 ] ) , ( [ radial ] ) and ( [ azimuthal ] ) by setting @xmath129 and @xmath139 , we obtain @xmath140\hbox{i}u\ , \\ & \omega j=\frac{r\kappa^2}{2\omega}\hbox{i}u\ , \\ & r=0\ , \qquad\qquad \omega z=\frac{\partial}{\partial r}(b_0\hbox{i}u)\ . \end{split}\ ] ] this set of equations is applicable as we proceed to investigate stationary axisymmetric coplanar mhd perturbations ( lou & zou 2004a , b ) . for coplanar hydrodynamic perturbations or coplanar mhd perturbations in an isopedically magnetized sid ( shu et al . 2000 ) , there is no essential difference for changing the order of limiting procedure for @xmath129 and @xmath141 . however , for msid with a coplanar azimuthal background magnetic field , the results would be different by changing the order of limiting procedure for @xmath129 and @xmath141 ( lou 2002 ; lou & fan 2002 ; lou & zou 2004a ) . to solve the poisson equation connecting the perturbed surface mass density and the perturbed gravitational potential , we may consistently assume disturbances to carry proper scale - free forms ( lynden - bell & lemos 1993 ) for aligned and logarithmic spiral perturbations . for aligned perturbations , @xmath111 contains only an amplitude variation in @xmath14 of perturbed surface mass density and does not involve phase variation in @xmath14 so that the maximum density perturbations at different radii line up in the azimuth . for spiral perturbations in comparison , in addition to an amplitude variation in @xmath14 , @xmath111 also involves a phase variation in @xmath14 such that a spiral pattern emerges .      for aligned global coplanar mhd perturbations , we select those perturbations that carry the same power - law variation in @xmath14 of the background equilibrium , namely @xmath142 where @xmath143 is a small constant coefficient and the numerical factor @xmath58 is given explicitly by @xmath144 where @xmath145 ( qian 1992 ; syer & tremaine 1996 ; shen & lou 2004a , b ) . in the limit of @xmath60 , we simply have @xmath146 consistent with earlier results on sids and msids ( shu et al . 2000 ; lou 2002 ; lou & shen 2003 ; lou & zou 2004a ; lou & wu 2004 ) . in fact , a more general class of self - consistent potential - density pairs satisfying the poisson integral ( [ perturb04 ] ) is also available ; see footnote 3 in shen & lou ( 2004b ) .      for aligned neutral modes of axisymmetry ( @xmath129 ) , we can start from equation ( [ m0 ] ) by setting @xmath134 ; this implies @xmath147 and gives no constraints on @xmath148 . as @xmath149 , the scale - free condition requires @xmath150 with @xmath151 and @xmath152 . it turns out that such stationary coplanar mhd perturbations are simply alternative states for the axisymmetric background equilibrium with a proper rescaling factor ( shu et al . 2000 ; lou 2002 ) . as this is somewhat trivial , we now turn our attention mainly to cases of non - axisymmetric @xmath127 stationary aligned coplanar mhd perturbations .      for the potential - density pair @xmath149 and @xmath153 , we have @xmath154 and @xmath152 according to equations @xmath155 . it then follows that @xmath156}{2m\omega r^2}\hbox{i}u\ , \\ & m\omega j+r(1-\beta)\omega \hbox{i}u=\frac{ma^2}{r\sigma_0}\bigg(\frac{r\sigma_0}{a^2}2\pi g { \cal p}_m-1\bigg)rs+\frac{c_a^2(1 - 4\beta)}{2\omega r}\hbox{i}u\ , \\ & z=-\bigg(\frac{1}{2}-2\beta\bigg)\frac{b_0\hbox{i}u}{m\omega r}\ , \qquad\qquad\qquad r=\frac{b_0u}{\omega r}\ . \end{split}\ ] ] rearranging the first three equations of ( [ alignedmge1 ] ) , we obtain a set of three algebraic equations , namely @xmath157 where coefficients @xmath158 , @xmath159 , @xmath160 ( @xmath161 ) are explicitly defined by @xmath162\ , \quad\qquad b_2\equiv m\omega+\frac{[(1 - 4\beta)^2 - 2m^2]c_a^2}{2m\omega r^2}\ , \quad\qquad c_2\equiv\frac{2\omega}{r}\ , \\ & a_3\equiv -mr\bigg[2\pi g{\cal p}_m(\beta ) -\frac{a^2}{\sigma_0 r}\bigg]\ , \qquad\qquad b_3\equiv r(1-\beta)\omega-\frac{(1 - 4\beta)c_a^2}{2\omega r}\ , \qquad\qquad c_3\equiv m\omega\ . \end{split}\ ] ] for nontrivial solutions of @xmath163 , the determinant of algebraic equations ( [ suj ] ) must vanish , namely @xmath164 using coefficient definitions ( [ abc ] ) in condition ( [ aligned0 ] ) , we obtain a more informative form of stationary dispersion relation @xmath165}{[m^2\omega^2r^2 -c_a^2(m^2 - 1/2 - 8\beta^2 + 4\beta ) ] [ ( 1 + 2\beta)\omega^2r^2-c_a^2(1/2 - 2\beta)]}\ \\ & \qquad\qquad\qquad -\frac{4\beta a^2 - 8\beta\pi g{\cal p}_m\sigma_0r + ( 1 - 4\beta)c_a^2}{2[m^2\omega^2r^2 -c_a^2(m^2 - 1/2 - 8\beta^2 + 4\beta)]}\ . \end{split}\ ] ] note that equation ( [ aligned1 ] ) takes a simpler form in the limit of @xmath60 ( i.e. , a flat rotation curve of a sid ) and reduces to equation ( 3.1.9 ) of lou ( 2002 ) as expected .    by using background mhd equilibrium variables ( [ propequilib ] ) in condition ( [ aligned1 ] ) , we obtain the stationary dispersion relation for global aligned @xmath127 coplanar mhd configurations in the form of a quadratic equation in terms of @xmath166 , namely @xmath167 where the three coefficients @xmath168 ( @xmath161 ) are defined by @xmath169q^2 -{\cal a}_m(1-f { \cal c}{\cal p}_m)\ , \\ & c_0\equiv -\frac{1 - 4\beta}{4(1 + 2\beta ) } \bigg[\frac{-2m^2+(1 - 4\beta)(1 - 2\beta)}{1 + 2\beta } f{\cal c}{\cal p}_m+1 - 4\beta\bigg]q^4 -\frac{-2m^2+(1 - 4\beta)(1 - 2\beta)}{2(1 + 2\beta ) } ( 1-f{\cal c}{\cal p}_m)q^2\ , \end{split}\ ] ] with auxiliary parameters explicitly defined by @xmath170\ . \end{split}\ ] ] which are exactly the same as those defined in shen & lou ( 2004b ) for aligned cases of global perturbations in a purely hydrodynamic composite disc system . here , the effect of coplanar magnetic field is associated with the @xmath81 parameter in definitions ( [ alignc2c1c0 ] ) . the determinant of quadratic equation ( [ aligned ] ) , @xmath171 , depends on various parameters and can become negative under some circumstances ( i.e. , no physical solutions for @xmath166 ) . for @xmath172 , the two real solutions of @xmath166 to equation ( [ aligned ] ) are @xmath173    as a first check of necessary consistency , we consider the limiting case of vanishing magnetic field , that is , @xmath174 . then , equation ( [ aligned ] ) simply reduces to @xmath175 which gives a non - trivial solution @xmath176 this @xmath87 solution ( [ alignednomsol ] ) is exactly the same result for the case of a single scale - free gas disc without magnetic field ( see syer & tremaine 1996 and subsection 3.1 of shen & lou 2004b ) .    as a second check of necessary consistency , we consider another limiting case of @xmath60 which corresponds to the msid case with a flat rotation curve studied by lou ( 2002 ) . one can readily show that equation ( [ aligned ] ) reduces to equation ( 3.2.11 ) of lou ( 2002 ) as expected . in contrast to the case of a single fluid disc studied by syer & tremaine ( 1996 ) without azimuthal magnetic fields , there are now two branches of @xmath87 solutions of algebraic equation ( [ aligned ] ) in general . a @xmath87 solution is considered to be physical when it satisfies both the condition @xmath86 and inequality ( [ restrt ] ) . when the magnetic field becomes sufficiently weak , among these two branches of @xmath87 solutions , @xmath177 is the counterpart of the single disc case and @xmath178 is additional due to the coupling between the surface mass density and magnetic field . we recall that in a composite system of one stellar disc and one gaseous disc coupled through the mutual gravitational interaction , there are also two branches of @xmath87 solutions ( lou & fan 1998b , 2000 ; lou & shen 2003 ; shen & lou 2004b ) . in a composite stellar - gaseous disc system , the two different classes of perturbation modes correspond to either in - phase or out - of - phase of surface mass density perturbation enhancements ( see also lou & wu 2004 ) . in a single azimuthally magnetized disc , the two coplanar mhd modes @xmath177 and @xmath178 will be distinguished by either in - phase or out - of - phase perturbation enhancements of the surface mass density and the azimuthal magnetic field in the wkbj or tight - winding approximation . in the wkbj regime , the @xmath177 and @xmath178 branches correspond to stationary fast and slow mhd density waves , respectively ( fan & lou 1996 ; lou & fan 1998a ) . although the phase relationships for perturbation enhancements of the surface mass density and the azimuthal magnetic field become @xmath179 for open spiral structures ( lou 2002 ; lou & fan 2002 ) , we shall still refer to @xmath177 and @xmath178 solutions as global stationary fast and slow mhd density waves , respectively . we shall come back to this discussion of ` phase relationship ' in more details in section 4 . as an example of illustration , we now focus on the specific case of @xmath11 with inequality ( [ restrt ] ) automatically satisfied . in such a magnetized disc system , the rotation curve scales as @xmath180 , the surface mass density scales as @xmath181 , the barotropic index is @xmath182 and the azimuthal magnetic field scales as @xmath183 ( i.e. , the background equilibrium magnetic field is force - free ) . for this specific case , the determinant @xmath184 of algebraic equation ( [ aligned ] ) remains always positive for both full and partial discs ( see proofs in appendix a ) . as seen from @xmath87 solutions ( [ solalign ] ) , the sign of @xmath185 determines whether @xmath177/@xmath178 is the upper / lower branch or reverse . for @xmath186 , we always have @xmath187 and therefore @xmath177 and @xmath178 remain to be the higher and lower branches , respectively . the situation becomes somewhat involved for the @xmath188 case when the sign of @xmath185 is dependent on the value of @xmath189 parameter . we show in fig . 1 the @xmath166 solution branches for @xmath188 and @xmath190 cases , respectively , versus the gravitational potential ratio @xmath189 .    for the @xmath188 case , there exists a diverging point at @xmath191 $ ] for @xmath192 solution branch where @xmath193 . the value of this critical @xmath194 is numerically determined to be @xmath195 . more specifically , we have @xmath196 for @xmath197 , while we have @xmath187 for @xmath198 . this sign variation of @xmath185 is indicated by the relative locations of @xmath177 and @xmath178 solution branches as shown in panel @xmath199 of fig . in other words , for @xmath200 , @xmath177 remains negative and the only physical solution is @xmath178 . on the other hand for @xmath198 , both @xmath177 and @xmath178 solution branches become physically plausible .    for @xmath11 , the explicit form of algebraic equation ( [ aligned ] ) becomes @xmath201d^4 -\bigg[\bigg(\frac{2{f}{\cal c}{\cal p}_m}{3 } + 1\bigg)m^2q^2+\bigg(m^2+\frac{3}{4}\bigg)(1-{f } { \cal c}{\cal p}_m)\bigg]d^2+\frac{2}{3}(1-{f } { \cal c}{\cal p}_m)m^2q^2=0\ , \ ] ] where auxiliary parameters @xmath202 , @xmath203 , @xmath204 , @xmath205 and @xmath206 are all evaluated at @xmath11 . it is fairly straightforward to show that @xmath207 for all @xmath127 when @xmath11 ( shen & lou 2004b ) . as @xmath208 , we then have @xmath209 .      in our formulation , stationary surface mass density perturbations are characterized by @xmath210 . for aligned coplanar mhd perturbations , we took @xmath211 where @xmath212 is a real exponent . in section 3.1 , we have chosen @xmath213 that carries the same power - law dependence as the equilibrium disc does . for a complex @xmath212 , a perturbation pattern would appear spiral , that is , a logarithmic spiral in the form of @xmath214 $ ] where @xmath215 and @xmath216 are real and imaginary parts of @xmath212 . to ensure the gravitational potential arising from this perturbed surface mass density being finite , we should require @xmath217 ( e.g. , qian 1992 ) . in fact , there exists a more general class of self - consistent potential - density pairs as indicated in footnote 7 of shen & lou ( 2004b ) . however , for the analysis of axisymmetric stability problem in section 3.2.1 , equation ( [ spirperturb ] ) is used such that dispersion relation ( [ m0u ] ) derived later on is real on both sides ( e.g. , lemos et al . 1991 ; syer & tremaine 1996 ; shu et al . 2000 ) . to be specific here , we consistently take the logarithmic spiral surface mass density perturbations and the resulting gravitational potentials as ( kalnajs 1971 ; lemos et al . 1991 ; shu et al . 2000 ; lou 2002 ; lou & shen 2003 ; shen & lou 2004b ; lou & zou 2004a ) @xmath218 where @xmath143 is a small constant , @xmath219 is a parameter related to the radial wavenumber and the kalnajs function @xmath220 is defined by @xmath221 ( kalnajs 1971 ) . as @xmath220 is an even function of @xmath219 , a consideration of @xmath222 would suffice . in our convention of notations , @xmath223 and @xmath224 correspond to leading and trailing logarithmic spiral waves , respectively , and @xmath225 relates the parameter @xmath219 to the radial wavenumber @xmath226 . we note that @xmath220 decreases monotonically with increasing @xmath223 , and @xmath227 for @xmath127 ; for @xmath129 , @xmath228 is positive and can be greater than 1 when @xmath219 becomes sufficiently small . the choice of unaligned perturbations is not unique , and the perturbation potential - density pairs may involve @xmath8 parameter such that the background plus coplanar mhd perturbations are altogether scale - free ( e.g. , syer & tremaine 1996 ) . as an example of illustration , we here presume coplanar logarithmic spiral perturbations defined by the density - potential pair ( [ spirperturb ] ) in the following analysis . for the specific @xmath11 case , the background surface mass density and the perturbed surface mass density bear the same radial scaling @xmath229 . we start from the axisymmetric coplanar mhd perturbation equation ( [ m0 ] ) with @xmath139 . for radial oscillations with a logarithmic potential - density pair , we presume a surface mass density @xmath230 , a gravitational potential @xmath231 and thus a radial speed @xmath232 according to the first of equation ( [ m0 ] ) . after proper rearrangements , the set of equations ( [ m0 ] ) can be combined to form a single equation in terms of @xmath112 , namely @xmath233u\\ & = \frac{a^2}{r^2}\bigg\{\bigg[1-{f}{\cal c } { \cal n}_0\bigg(1+d^2-\frac{1 - 4\beta}{2 + 4\beta}q^2\bigg)\bigg ] \bigg(\frac{1}{4}+\xi^2\bigg ) + q^2\bigg(\frac{1}{4}+\xi^2-\beta\bigg)\bigg\}u\ , \end{split}\ ] ] where @xmath234 , @xmath235 $ ] and @xmath236 is extremely small .    by inspection , equation ( [ m0u ] ) bears similar features of the dispersion relation in the classic wkbj regime involving magnetic field . as the right - hand side of equation ( [ m0u ] ) is real , axisymmetric instabilities first set via neutral @xmath134 modes . it should be emphasized that the scale - free condition is met only for neutral modes . with @xmath237 , the marginal stability curves are then given by @xmath238\}(\xi^2 + 1/4)+q^2(\xi^2-\beta+1/4)}{{f } { \cal c}{\cal n}_0(\xi^2 + 1/4)+(1 + 2\beta)(-2 + 2\beta)}\ .\ ] ] as a check of necessary consistency , the axisymmetric marginal stability curve for an msid when @xmath19 ( @xmath239 ) then gives @xmath240 ( \xi^2 + 1/4)}{{f}{\cal n}_0(\xi^2 + 1/4)-2}\ .\ ] ] one can readily show ) is automatically satisfied for @xmath85 and we need only to examine the range of @xmath88 . when the denominator of the right - hand side of equation ( [ msc ] ) is negative , the numerator must also be negative for @xmath86 . it then follows that @xmath241 and inequality ( [ restrt ] ) is met . when the denominator of the right - hand side of equation ( [ msc ] ) is positive , it is easy to show that inequality ( [ restrt ] ) holds . as solution ( [ mscbeta0 ] ) is a special case of solution ( [ msc ] ) , inequality ( [ restrt ] ) is thus met . ] that stationary dispersion relations ( [ msc ] ) and ( [ mscbeta0 ] ) automatically satisfy inequality ( [ restrt ] ) for physical solutions @xmath86 . in the msid case with different values of @xmath77 parameter for magnetic field strengths , one readily obtains the marginal stability curves that separate two distinct unstable regimes , namely , the collapse regime for large radial spatial scales and the ring fragmentation regime for relatively short radial wavelengths ( lemos et al . 1991 ; shen & lou 2003 , 2004a ) . there exists a critical @xmath242 where @xmath87 diverges and this @xmath243 is determined by @xmath244 . for @xmath245 , we have @xmath246 and thus an enhancement of @xmath77 tends to suppress the collapse boundary . for @xmath247 , we have @xmath248 and thus an enhancement of @xmath77 tends to raise the ring fragmentation boundary . in short , an enhancement of ring magnetic field tends to reduce dangers of both collapse and ring fragmentation instabilities in the msid case . while the enhancement of ring magnetic field tends to suppress ring fragmentation instabilities in general , it only suppresses collapse instabilities for @xmath88 and tends to aggravate collapse instabilities for @xmath89 ( see fig . 2 for details ) . it has been shown that the magnetic force is radially inward and outward for @xmath88 and @xmath89 , respectively . why does the trend of variations for the collapse - stability appear seemingly paradoxical ? this situation can be understood because the unperturbed background is in an mhd radial force balance and the surface mass density , the rotation speed and the magnetic field are coupled explicitly through condition ( [ propequilib ] ) . furthermore , for fast mhd density waves with @xmath129 ( lou & fan 1998a ) , the gas pressure and magnetic pressure together is associated with the radial wavenumber squared , while the surface mass density and self - gravity are associated linearly with the radial wavenumber ( see the first line of dispersion relation [ m0u ] with small @xmath236 ) . for @xmath88 , the surface mass density becomes smaller for stronger magnetic field strength , while for @xmath89 , the situation reverses . for the ring fragmentation instability at relatively large radial wavenumbers , the dominant magnetic pressure effect tends to stabilize the disc . for the jeans collapse instability at relatively small radial wavenumbers , the dominant self - gravity effect is proportional to the background surface mass density ; the roles of magnetic field for two different situations can then be readily understood . it is the coupling effect between the surface mass density and the magnetic field of the background that gives rise to this collapse feature .    the marginal stability curves for scale - free discs with @xmath249 and without magnetic fields are @xmath250 that are consistent with those for the single - disc case of shen & lou ( 2004b ; e.g. , see their subsection 3.2.1 ) . the marginal axisymmetric stability curves generally consist two branches , namely , the collapse branch and the ring - fragmentation branch ( lemos et al . 1991 ; shu et al . 2000 ; lou 2002 ; lou & shen 2003 ; shen & lou 2003 ) . the lowest value of @xmath87 for stability is the maximum of the collapse branch and the highest value of @xmath87 for stability is the minimum of the ring - fragmentation branch . we note that the maximum of the collapse branch is always located at @xmath251 for @xmath252 , while for @xmath253 the maximum of the collapse branch may locate at @xmath223 with @xmath251 corresponding to a local minimum of @xmath87 . it was proven in appendix c of shen & lou ( 2004b ) that @xmath251 is always a local extremum for @xmath87 . variations of the stable range for @xmath87 with parameters @xmath254 , @xmath8 and @xmath77 are shown in panels ( a ) and ( b ) of fig . finally we note that the axisymmetric destabilization involves only the fast mhd wave branch because the slow mhd wave branch is negative and thus unphysical . for non - axisymmetric ( @xmath255 ) perturbations , both fast and slow mhd density waves are possible and stationary global configurations may represent transitions from stable to unstable situations ( e.g. , shu et al . 2000 ) . here , the stability of our slow mhd density waves differs from those models involving finite disc thickness where either vertical or horizontal weak magnetic fields may induce velikhov - chandrasekhar - balbus - hawley instability and magneto - rotational instabilities ( mri ) through slow mhd modes ( chandrasekhar 1961 ; balbus & hawley 1998 ; kim & ostriker 2000 ; lou , yuan & fan 2001 ) .    it would be interesting to recall the axisymmetric stability analysis in the wkbj approximation . the local dispersion relation for fast mhd waves in an azimuthal magnetic field is ( fan & lou 1996 ) @xmath256 where @xmath226 is the radial wavenumber and other notations have the conventional meanings . equation ( [ drwkbj ] ) gives the marginal stability curve for axisymmetric @xmath129 perturbations as @xmath257 by using scale - free disc equilibrium condition ( [ propequilib ] ) and inserting @xmath258 for large @xmath226 and @xmath219 into dispersion relation ( [ stability ] ) , we readily obtain the corresponding marginal stability curve in the wkbj regime as @xmath259\xi } { f{\cal c}\xi+(1 + 2\beta)(-2 + 2\beta)}\ , \ ] ] which bears a strong resemblance to expression ( [ msc ] ) of our global analysis , especially in view of the asymptotic expression for @xmath220 , namely @xmath260 when @xmath261 . another equivalent way of this comparison is to note the parallel terms in the wkbj condition ( [ stability ] ) and in the first line of the corresponding global condition ( [ m0u ] ) , also in reference to the above asymptotic expression for @xmath220 . a specific illustrating example of such a comparison is shown in fig . 3 with parameters @xmath19 ( @xmath239 ; flat rotation curve ) , @xmath35 ( full disc ) and @xmath262 . in the short - wavelength limit ( i.e. , @xmath263 ) , the wkbj approximation is well justified and the results are qualitatively consistent , while in the long - wavelength of jeans collapse regime , the wkbj analysis differs from the global exact treatment significantly .    , @xmath35 and @xmath262 . the marginal stability curves for the long - wavelength jeans collapse and the short - wavelength ring fragmentation regimes are denoted . for @xmath263 ( i.e. , the short - wavelength limit ) , the wkbj approximation is well applicable , while for small @xmath219 the results in the wkbj approximation deviates from those of the exact global treatment considerably ( see shen & lou 2003 , 2004b ) . ]      in parallel with the aligned @xmath127 cases , we now construct unaligned logarithmic spiral solutions from coplanar mhd perturbation equations ( [ stationary1])@xmath128([stationary4 ] ) . for the density - potential pair @xmath230 and @xmath264 , we then have radial flow speed @xmath265 and @xmath266component specific angular momentum @xmath267 . coplanar mhd perturbation equations @xmath155 can therefore be written in the forms of @xmath268\hbox{i}u\ , \\ & m\omega j+r(1-\beta)\omega \hbox{i}u = m\bigg(2\pi g{\cal n}_m -\frac{a^2}{r\sigma_0}\bigg)rs+\frac{c_a^2(1 - 4\beta ) } { 2\omega r}\hbox{i}u\ , \\ & z=-\hbox{i}\xi\frac{b_0iu}{m\omega r}\ , \qquad\qquad\qquad\qquad r=\frac{b_0u}{\omega r}\ . \end{split}\ ] ] the first three equations of set ( [ spiralmge1 ] ) immediately lead to @xmath269 where coefficients @xmath158 , @xmath159 and @xmath160 ( @xmath161 ) are explicitly defined by @xmath270c_a^2}{m\omega r^2}\ , \qquad c_2\equiv \frac{2\omega}{r}\ , \\ & a_3\equiv -mr\bigg(2\pi g{\cal n}_m-\frac{a^2}{\sigma_0 r}\bigg)\ , \qquad\qquad b_3\equiv r(1-\beta)\omega-\frac{(1 - 4\beta)c_a^2}{2\omega r}\ , \qquad\qquad c_3\equiv m\omega\ . \end{split}\ ] ] for nontrivial solutions of the set of homogeneous algebraic equations ( [ sujspiral ] ) , the determinant must vanish @xmath271 by using the set of expressions ( [ abcspiral ] ) and applying equilibrium conditions ( [ propequilib ] ) , we obtain the final dispersion relation for stationary logarithmic spiral @xmath127 mhd configurations in the form of a quadratic equation in terms of @xmath166 , namely @xmath272 where the coefficients @xmath185 , @xmath273 and @xmath274 are explicitly defined by @xmath275q^2 -{\cal a}_m(1-{f}{\cal c}{\cal n}_m)\ , \\ & c_0\equiv -\frac{(1 - 4\beta ) } { 4(1 + 2\beta ) } \bigg[\frac{(-2m^2 - 2\xi^2 - 2\beta+1/2)}{(1 + 2\beta ) } { f}{\cal c}{\cal n}_m+1 - 4\beta\bigg]q^4 -\frac{(-2m^2 - 2\xi^2 - 2\beta+1/2)}{2(1 + 2\beta ) } ( 1-{f}{\cal c}{\cal n}_m)q^2\ , \end{split}\ ] ] with auxiliary parameters defined by @xmath276\ . \end{split}\ ] ] that are the same as those in shen & lou ( 2004b ) for unaligned logarithmic spiral cases . for a positive determinant @xmath277 , the two real solutions to quadratic equation ( [ spiral ] ) are given by @xmath278    similar to @xmath127 cases for aligned perturbations in the limit of vanishing magnetic field ( i.e. , @xmath279 ) , equation ( [ spiral ] ) reduces to @xmath280 that leads to a nontrivial solution @xmath281 solution ( [ spiralnomsol ] ) is exactly the same result for the single - disc case discussed in section 3.2 of shen & lou ( 2004b ) . another limiting regime of @xmath60 has been analyzed by lou ( 2002 ) for global stationary perturbation structures in a single msid .      for stationary spiral mhd configurations , we here analyze again the @xmath11 case as a specific example of illustration . the background equilibrium ring magnetic field is force - free and inequality ( [ restrt ] ) is satisfied . the determinant @xmath184 of quadratic equation ( [ spiral ] ) remains always positive as shown in appendix a and hence two real solutions of @xmath282 to equation ( [ spiral ] ) exist . in parallel with the aligned @xmath11 case , we solve quadratic equation ( [ spiral ] ) for given values of parameter @xmath283 and display in diagrams the two branches of solutions @xmath177 and @xmath178 ( @xmath166 ) versus @xmath219 variation . as @xmath284 for @xmath188 and when @xmath283 becomes sufficiently small , there is a diverging point for @xmath192 solution where @xmath193 . since @xmath285 increases monotonically with increasing @xmath219 for @xmath222 and attains its minimum value at @xmath251 ( @xmath285 is an even function of @xmath219 ; see appendix c of shen & lou 2004b ) , the existence of such a critical point @xmath243 for a diverging @xmath166 requires @xmath286 for the specific @xmath11 case , this inequality ( [ criticalf ] ) implies a critical @xmath283 value of @xmath287 $ ] below which there exists one divergent point of @xmath243 where @xmath177 diverges . numerically , we have determined @xmath288 [ note that this @xmath289 is of the same value as that of the aligned @xmath188 case because @xmath290 . for selected @xmath254 values in the range of @xmath291 , the corresponding critical @xmath243 at which @xmath292 can be numerically computed ; for example , we have @xmath293 for @xmath294 . it becomes much simpler for cases of @xmath186 because @xmath185 remains always positive and there is no diverging point for @xmath87 solutions . therefore , @xmath177 and @xmath178 remain always as the upper and lower branches , respectively . for @xmath188 and @xmath190 cases with different parameters of @xmath283 and @xmath77 , we show in four panels @xmath295 of fig . 4 relevant curves of @xmath166 versus @xmath219 variation to illustrate the basic features . to be more explicit for the case of @xmath11 in quadratic equation ( [ spiral ] ) , we write @xmath296d^4 -\bigg[(m^2+\xi^2)\bigg(\frac{2}{3 } { f}{\cal c}{\cal n}_m+1\bigg)q^2\\ & \hbox { } \qquad\qquad\qquad\qquad\qquad\qquad\qquad + \bigg(m^2+\xi^2+\frac{3}{4}\bigg ) ( 1-{f}{\cal c}{\cal n}_m)\bigg]d^2 + \frac{2}{3}(m^2+\xi^2)(1-{f}{\cal c}{\cal n}_m)q^2=0\ , \end{split}\ ] ] where auxiliary parameters @xmath297 , @xmath203 , @xmath204 , @xmath205 and @xmath206 are all evaluated for the special value of @xmath11 . for multi - wavelength observations of nearby spiral galaxies ( sofue et al . 1986 ; kronberg 1994 ; beck et al . 1996 ) , it is possible to identify large - scale spatial phase relationships of spiral patterns using various observational diagnostics ( e.g. , mathewson et al . 1972 ; visser 1980a , b ; neininger 1992 ; beck & hoernes 1996 ; frick et al . 2000 , 2001 ; ferrire 2001 ; lou et al . 2002 ) . for this reason , we derive below spatial phase relationships for various coplanar mhd perturbation variables , keeping in mind of potential applications to magnetized spiral galaxies ( e.g. , fan & lou 1996 ; lou & fan 1998a ) .      here we systematically examine the spatial phase relationships among enhancements of gas surface mass density , magnetic field and velocity perturbations by analysing the relationships among @xmath298 , @xmath148 , @xmath73 , @xmath299 and @xmath300 for aligned and logarithmic spiral cases , separately . starting from the set of homogeneous algebraic equations ( [ suj ] ) for the aligned cases , we obtain the relationships among @xmath298 , @xmath301 and @xmath299 in the following , namely @xmath302 where @xmath158 , @xmath159 and @xmath160 ( @xmath161 ) are defined by expressions ( [ abc ] ) . by using definitions ( [ abc ] ) and equations ( [ alignedmge1 ] ) , we obtain @xmath303 for the aligned cases , where the common dimensionless real factor @xmath304 is defined below by @xmath305    similarly , by using the set of homogeneous algebraic equations ( [ sujspiral ] ) , definitions ( [ abcspiral ] ) and equations ( [ spiralmge1 ] ) , we readily derive for the cases of global stationary logarithmic spiral configurations @xmath306\ , \end{split}\ ] ] where the common dimensionless complex factor @xmath307 is defined below by @xmath308    the two sets of algebraic equations ( [ relsurjz ] ) and ( [ relsurjzspiral ] ) are the most general phase relationships for aligned and logarithmic spiral cases , respectively . once a @xmath8 value is specified and a particular @xmath87 for a stationary solution inserted , we immediately obtain the phase relationship between any two perturbation variables . as an example of illustration , we examine again the special @xmath11 case and determine the relevant phase relationships in the next subsection .      when @xmath11 , the results of phase relations derived in the last subsection 4.1 can be simplified considerably as all terms involving the factor @xmath309 vanish . meanwhile , inequality ( [ restrt ] ) is automatically satisfied . the only requirement for a physical @xmath87 solution is simply @xmath86 . as already noted in subsections 3.1.3 and 3.2.3 , there are two branches of @xmath87 solutions . we insert relevant physical @xmath87 solutions in equations ( [ relsurjz ] ) and ( [ relsurjzspiral ] ) to determine phase relationships among perturbation variables . in particular , we are interested in the phase relationship between the perturbed surface mass density and the perturbed azimuthal magnetic field . we now analyze below this mass - magnetic field phase relationships for aligned and spiral cases , separately . some illustrative examples are shown in fig . 5 and fig . 6 to provide more direct visual impressions .      for aligned @xmath11 cases , we readily find @xmath310 which means the azimuthal magnetic field remains unchanged . from equation ( [ relsurjz ] ) , we derive @xmath311 meanwhile , equation ( [ align0d25 ] ) can be rewritten in the form of @xmath312 } = \frac{d^2\{d^2-(1-f{\cal c}{\cal p}_m)/ [ f{\cal c}{\cal p}_m+3(m^2 - 3/2)/(2m^2 + 3/2)]\ } } { d^2-(1-f{\cal c}{\cal p}_m)/(f{\cal c}{\cal p}_m+3/2 ) } \ , \ ] ] that contains two distinct @xmath87 solutions according to subsection 3.1.3 , namely , the plus - sign solution @xmath177 and the minus - sign solution @xmath178 in the solution form ( [ solalign ] ) , respectively . as shown in appendix b , we emphasize that the plus - sign solution @xmath177 ( if positive and thus physical ) makes @xmath313 and the minus - sign solution @xmath178 makes @xmath314 on the right - hand side of equation ( [ q2aligned ] ) , respectively . therefore , the two branches of @xmath87 solutions ( if positive and thus physical ) will give either @xmath315 for the plus - sign @xmath177 solution or @xmath316 for the minus - sign @xmath178 solution .      for spiral cases with @xmath11 , we readily obtain from equations ( [ relsurjzspiral ] ) and ( [ thetas ] ) that @xmath317(3 + 2\hbox{i}\xi)\ .\ ] ] following the same procedure , we can show that the two branches of @xmath87 solutions , namely , the plus - sign solution @xmath177 and the minus - sign solution @xmath178 given by equation ( [ solspiral ] ) , if positive and thus physical , correspond to the factor @xmath318 being positive and negative , respectively ( see appendix b for a detailed proof ) . we now examine the phase relationship between coplanar perturbation enhancements of surface mass density and magnetic field .    in the short - wavelength limit ( i.e. , @xmath263 with a logarithmic spiral pattern tightly wound ) , we have @xmath319 and @xmath320 where the plus- and minus - signs correspond to the @xmath177 and @xmath178 solutions , respectively . the phase relations between perturbation enhancements of surface mass density and magnetic field are in - phase for the @xmath177 solution ( if positive and physical ) and out - of - phase for the @xmath178 solution . for @xmath186 , @xmath177 remains always greater than @xmath178 . in contexts of magnetized spiral galaxies , these solution properties in the wkbj regime are consistent with the physical identifications of the @xmath177 mode with a stationary fast mhd density wave and the @xmath178 mode with a stationary slow mhd density wave , respectively ( fan & lou 1996 ; lou & fan 1998a ; lou & fan 2002 ) . it is now possible to model magnetized spiral galaxies in terms of coplanar stationary fast and slow mhd density waves with more general rotation curves including the case of a flat rotation curve for a magnetized singular isothermal disc ( shu et al . 2000 ; lou 2002 ; lou & fan 2002 ; lou & zou 2004a , b ; lou & wu 2004 ) .    in the long - wavelength limit ( i.e. , @xmath321 with a logarithmic spiral pattern widely open ) , we have @xmath322 and @xmath323 where the minus- and plus - signs correspond to the @xmath177 and @xmath178 solutions , respectively . therefore , the coplanar perturbation enhancement of magnetic field is either ahead of or lagging behind the perturbation enhancement of surface mass density by a phase difference of @xmath324 .    for intermediate radial wavelengths , the phase difference between perturbation enhancements of surface mass density and magnetic field in a logarithmic spiral pattern can be readily determined for given specific conditions . on the basis of our model analysis here , it appears that except for global stationary fast mhd density waves in the extreme wkbj regime , magnetic and gas spiral arms tend to be phase shifted relative to each other for a global stationary logarithmic spiral pattern for a much broader set of rotation curves including the special case of a flat rotation curve . one can readily show that for aligned global stationary mhd perturbation configurations , there is no net radial angular momentum flux as aligned stationary mhd density waves do not propagate radially ( shu et al . 2000 ; lou 2002 ) . we here briefly discuss the process of angular momentum transfer in association with the stationary logarithmic spiral mhd modes obtained in this paper . according to the analyses of lynden - bell & kalnajs ( 1972 ) and fan & lou ( 1999 ) , the total angular momentum flux carried by coplanar mhd density waves should contain three separate contributions , namely , the advective transport flux @xmath325 defined by @xmath326 the gravity torque flux @xmath327 defined by @xmath328\bigg[\frac{\partial\re(\phi_1 ) } { \partial r}\bigg]\\ ] ] and the magnetic torque flux @xmath329 defined by @xmath330    the perturbed three - dimensional gravitational potential @xmath331 associated with a fourier component of a coplanar logarithmic spiral perturbation in surface mass density @xmath332 is then @xmath333 where @xmath334 is the cylindrical bessel function of order @xmath109 with an argument @xmath23 ( e.g. , binney & tremaine 1987 ) . in fact , integral ( [ 3dphi1 ] ) gives a simple form of kalnajs function in the disc plane at @xmath335 with the real part of @xmath331 explicitly given by @xmath336 it is fairly straightforward to show that @xmath325 and @xmath329 are independent of @xmath14 where the phase relationships for stationary logarithmic spiral cases derived in section 4.1 have been use . while mathematically tedious , one can further show that @xmath327 is also independent of @xmath14 ( see appendix c ) . therefore , the total angular momentum flux @xmath337 remains constant at all radii , a natural result of angular momentum conservation given that perturbation surface mass density falls off as @xmath181 . therefore , in order to maintain a stationary logarithmic spiral pattern , there must be a source or a sink at the disc center . the perturbation energy flux , however , is zero at all radii for the stationary coplanar mhd configurations [ see equation ( 4.17 ) of fan & lou ( 1999 ) in the wkbj regime ] . the mass flux is also zero at all radii since the first - order mass flux @xmath338 , and the second - order mass flux @xmath339 [ where @xmath340 is the second - order perturbed radial velocity , this equation holds due to the second - order mass conservation ( fan & lou 1999 ) ] . therefore our model result differs from that of spruit ( 1987 ) for stationary logarithmic spiral shocks , where only the gravity of the central object is considered . in spruit s model , there is a net mass accretion through the disc and an effective viscosity parameter `` @xmath7 '' ( shakura & sunyaev 1973 ) may be estimated . physically , the central source of angular momentum required to sustain global stationary logarithmic spiral patterns may be plausibly identified with a fast spinning supermassive black hole ( smbh ) in the nucleus of a spiral galaxy . the kerr spacetime associated with such a spinning smbh excites and sustains mhd density waves in the surrounding disc through a long - range gravitational interaction . for trailing spiral patterns , a net flux of angular momentum is transport radially outward in the disc . conceptually , this process should cause the spinning smbh to gradually slow down , entirely analogous to the spin - down of a rotating neutron star by emitting electromagnetic waves that carry a net angular momentum away . in this paper , we have explored and analyzed stationary perturbation structures in coplanarly magnetized razor - thin scale - free discs in a more general manner . to be specific , we start from a rotational and magnetized background equilibrium of axisymmetry that is dynamically self - consistent with a surface mass density @xmath341 , a rotation curve @xmath2 , a purely azimuthal ( ring ) magnetic field @xmath1 with a vertically integrated barotropic equation of state in the form of @xmath16 . the radial force balance at all radii in an mhd disc implies several simple relationships among these power indices @xmath7 , @xmath8 , @xmath9 and @xmath6 , as stated explicitly in equation ( [ pararelation ] ) . without loss of generality , we can simply use @xmath8 as an independent parameter to specify properties of an important class of radial variations for rotational mhd background equilibria . the allowed range of @xmath8 falls within the interval @xmath71 ( syer & tremaine 1996 ; shen & lou 2004a , b ) with the special case of @xmath19 corresponding to a magnetized singular isothermal disc ( msid ) system ( shu et al . 2000 ; lou 2002 ; lou & fan 2002 ; lou & zou 2004a , b ; lou & wu 2004 ) . for clarity and convenience of our analysis , we introduced several dimensionless parameters . the first parameter is the partial disc parameter @xmath34 for the ratio of the gravitational potential ( @xmath342 ) arising from the disc to that ( @xmath28 ) arising from the entire system ( syer & tremaine 1996 ; shu et al . 2000 ; lou 2002 ) including an axisymmetric halo mass distribution that is involved in the background equilibrium but is presumed to be unresponsive to coplanar mhd disc disturbances ( a massive dark matter halo is an example in mind ) ; for a partial disc system , we have @xmath343 , while for a full disc system , we have @xmath344 by conventional definitions . the second parameter is an effective rotational mach number @xmath82 defined as the ratio of the disc rotation speed to the sound speed yet with an additional scaling factor @xmath345 ; @xmath82 is constant at all radii . the third parameter @xmath77 is a measure for the background ring magnetic field strength defined as the ratio of the azimuthal alfvn speed to the sound speed ; and @xmath77 remains constant at all radii . for @xmath89 , there is no constraint on @xmath77 and the net lorentz force of the background azimuthal magnetic field is radially outward ( the magnetic pressure force is stronger than the magnetic tension force ) to resist the self - gravitation ; for @xmath11 , the background azimuthal magnetic field is force - free ( the magnetic pressure and tension forces cancel each other exactly ) ; for @xmath88 , the net lorentz force of the background azimuthal magnetic field points radially inward ( the magnetic pressure force is weaker than the magnetic tension force ) to aggravate the self - gravity and thus a restriction [ i.e. , inequality ( [ restrt ] ) ] is imposed on the magnetic field strength in order to guarantee that the surface mass density be positive in equation ( [ propequilib ] ) . inequality ( [ restrt ] ) provides a necessary physical criterion for plausible @xmath87 solutions .    with such a rotational mhd background equilibrium chosen , we introduce coplanar mhd disturbances to construct global stationary perturbation structures as viewed in an inertial frame of reference . perturbation variables are expressed in terms of fourier components for either aligned or logarithmic spiral pattern forms ( syer & tremaine 1996 ; shu et al . 2000 ; lou 2002 ; lou & shen 2003 ; shen & lou 2003 ; lou & zou 2004a ) . we derive analytical solutions for stationary coplanar mhd perturbation configurations in both aligned and logarithmic spiral cases . we now summarize key results below . 1 .   cases of aligned coplanar mhd perturbations 0.2 cm + for aligned cases , we choose perturbation variables that bear the same power - law dependence in radius @xmath14 as the background mhd equilibrium does . for the special aligned case of @xmath129 , the stationary coplanar mhd perturbations actually describe alternative equilibria to the axisymmetric background equilibrium with proper rescaling factors ( shu et al . 2000 ; lou 2002 ; lou & zou 2004a ) and are therefore somewhat trivial . + for aligned cases of @xmath127 , we derive a quadratic equation ( [ aligned ] ) in terms of @xmath166 . the resulting positive @xmath87 is the rotational parameter required for sustaining a stationary nonaxisymmetric mhd perturbation configuration . other parameters involved in quadratic equation ( [ aligned ] ) include @xmath8 , @xmath109 , @xmath283 and @xmath77 . for a general combination of @xmath8 , @xmath109 , @xmath283 and @xmath77 , the determinant @xmath184 of quadratic equation ( [ aligned ] ) may not be always positive definite . on the other hand , for specific values of @xmath8 , equation ( [ aligned ] ) always gives two different real solutions ( see appendix a ) . once we prescribe parameters @xmath8 , @xmath109 , @xmath283 and @xmath77 , the two solutions of @xmath87 can be readily determined . for physical solutions , @xmath87 must be positive and satisfy inequality ( [ restrt ] ) simultaneously . + as a specific case study , we examined the case of @xmath11 . here the background azimuthal magnetic field remains force - free and the @xmath87 solution is free from constraint ( [ restrt ] ) . it has been further shown in appendix a that for this @xmath11 case , equation ( [ aligned ] ) always has two distinct real solutions , the plus - signed solution @xmath177 and the minus - signed solution @xmath178 , as contained in equation ( [ solalign ] ) . all these aspects contribute to a significant simplification of the relevant computational procedures . we derive two branches of @xmath87 solutions and need to require only @xmath346 . there is yet a qualitative difference between the @xmath188 case and the @xmath186 cases . for @xmath188 , there exists one critical value @xmath347 ; for @xmath283 below @xmath289 , @xmath177 and @xmath178 are negative and positive , respectively , while @xmath283 above @xmath289 , both branches of @xmath87 solutions are positive and thus physical . for @xmath186 in comparison , both branches of @xmath87 solutions remain always positive and thus physical . moreover , @xmath177 remains always to be the upper branch for azimuthal stationary fast mhd density waves and represents the hydrodynamic counterpart of the unmagnetized single disc case ; complementarily , @xmath178 remains always to be the lower branch for azimuthal stationary slow mhd density waves ( fan & lou 1996 ; lou & fan 1998a ) and is caused by the very presence of the azimuthal magnetic field ( detailed computations reveal that @xmath178 must be smaller than @xmath348 which seems to be too small for galactic applications ) . when a disc rotates sufficiently fast , it could only support the upper @xmath177 stationary mhd configurations . an increase of magnetic field parameter @xmath77 will raise both physical branches of @xmath87 solutions as displayed in fig . the variation trend of the upper @xmath177 branch appears to be more sensitive . therefore for a much stronger magnetic field , a sufficiently rapidly rotating disc may sustain stationary coplanar mhd configurations in terms of the upper @xmath349 branch . + 0.2 cm 2 . cases of logarithmic spiral configurations for coplanar mhd perturbations 0.2 cm + for cases of stationary spiral perturbation structures , we choose coplanar mhd perturbations in terms of kalnajs logarithmic spirals ( kalnajs 1971 ; lemos et al . 1991 ; shu et al . 2000 ; lou 2002 ; lou & fan 2002 ; lou & shen 2003 ; lou & zou 2004a ; lou & wu 2004 ) . for the special @xmath129 case with radial propagations , it happens that inequality ( [ restrt ] ) is automatically satisfied and the stationary coplanar mhd configurations with @xmath350 obtained here represent marginal stability curves ( syer & tremaine 1996 ; shu et al . 2000 ; shen & lou 2003 , 2004b ; lou & zou 2004b ) . only when the rotational parameter @xmath87 falls within a specific finite range can the magnetized disc be stable against axisymmetric disturbances at all wavelengths . a disc with too slow a rotation speed will succumb to jeans instability in the collapse regime corresponding to long wavelength perturbations , while a disc with too fast a rotation speed will suffer the ring - fragmentation instability corresponding to relatively short wavelength perturbations ( safronov 1960 ; toomre 1964 ; lemos et al . 1991 ; lou & fan 1998a , 2000 ; shu et al . 2000 ; shen & lou 2003 , 2004a ; lou & zou 2004b ) . in comparison to a hydrodynamics disc system , the enhancement of azimuthal magnetic field tends to suppress ring - fragmentation instabilities as expected from the perspective of the mhd @xmath351 parameter ( lou & fan 1998a ; lou 2002 ; lou & zou 2004a , b ) , while an enhancement of azimuthal magnetic field tends to suppress jeans collapse instabilities for @xmath88 , to aggravate jeans collapse instabilities for @xmath89 and to bear no effect on jeans collapse instabilities when @xmath11 . we have provided intuitive interpretations for the roles of magnetic field earlier . + in parallel with the aligned cases , we derive quadratic equation ( [ spiral ] ) in terms of @xmath166 for logarithmic spiral cases of @xmath127 . in addition to the four parameters @xmath8 , @xmath109 , @xmath283 and @xmath77 , we now have one more parameter @xmath219 for the radial wavenumber . while the determinant @xmath184 of quadratic equation ( [ spiral ] ) may not be always positive in general , we prove in appendix a that @xmath352 for the specific @xmath11 case . once the five parameters @xmath8 , @xmath109 , @xmath283 , @xmath77 and @xmath219 are specified , the @xmath87 solutions can be readily obtained . for physical solutions , @xmath87 must be non - negative and also satisfy inequality ( [ restrt ] ) . + as an example of illustration , we study the specific @xmath11 case inequality ( [ restrt ] ) being satisfied automatically . not surprisingly , the aligned cases and the logarithmic spiral cases parallel with each other very well , because of a mere change of wave propagation direction . for the special spiral case of @xmath188 , there again exists the same critical @xmath347 for @xmath283 value as in the aligned @xmath188 case ; for @xmath353 , there is one critical point of @xmath243 where @xmath354 diverges . more specifically , we have unphysical @xmath355 for @xmath356 , while we have physical @xmath357 for @xmath247 . there is no divergent point for @xmath87 solution branch and we always have @xmath358 ( @xmath359 when @xmath174 ) . the solution structures become simpler for @xmath186 cases when both @xmath87 branches of solutions are non - negative and there is no divergent point of @xmath354 . moreover , @xmath177 remains always to be the upper branch for stationary spiral fast mhd density waves and is the hydrodynamic counterpart for those in the unmagnetized single disc case , while @xmath178 remains always to be the lower branch for stationary spiral slow mhd density waves ( fan & lou 1996 ) and is caused by the very presence of the magnetic field ( detailed computations reveal @xmath178 must be smaller than @xmath348 ) . + 0.2 cm 3 . phase relationships among coplanar mhd perturbation variables 0.2 cm + so far , we have constructed global stationary coplanar mhd configurations for both aligned and logarithmic spiral cases in magnetized discs with a range of rotation curves . as expected , the inclusion of azimuthal magnetic field gives rise to one additional distinct physical solution for stationary slow mhd density waves to the problem ( fan & lou 1996 ; lou & fan 1998a ; lou 2002 ; lou & zou 2004a ) ; the alfvnic modes are excluded by our restricted consideration for coplanar mhd perturbations . the two @xmath87 branches of solutions are at least mathematically reasonable and may both have applications to astrophysical disc systems such as magnetized disc galaxies and so forth . these devised mhd disc problems are highly idealized with yet precious global analytical perturbation solutions , and conceptually still belong to a subclass of more general time - dependent mhd perturbation solutions . for applications to magnetized spiral galaxies , it would be of considerable interest to discuss pertinent physical aspects and identify relevant observational diagnostics . + for physical solutions of @xmath86 , the @xmath177 branch would be applicable to mhd discs with relatively fast rotation while the @xmath178 branch would be applicable to mhd discs with relatively slow rotation in general . for spatial phase relationships , the two branches of @xmath86 solutions give different results . we take the specific @xmath11 case as an example of illustration . for corresponding aligned cases , the plus - solution @xmath177 will lead to @xmath315 and the minus - solution @xmath178 will lead to @xmath316 ; there are no azimuthal magnetic field perturbations for aligned @xmath11 cases . for tightly wound logarithmic spiral cases , perturbations of both radial and azimuthal magnetic field are approximately in phase with the perturbation enhancement of the surface mass density for the @xmath177 branch , while they are approximately out of phase for the @xmath178 branch . for open logarithmic spirals , the radial and azimuthal magnetic field perturbations are either ahead of or lagging behind the enhancement of the surface mass by a significant phase difference . these phase relationships together with those associated with flow perturbations ( visser 1980a , b ) are valuable to interpret the spatial phase shifts between optical arms and magnetic arms in lopsided , barred or other spiral galaxies ( mathewson et al . 1972 ; neininger 1992 ; beck & hoernes 1996 ; beck et al . 1996 ; fan & lou 1996 ; lou & fan 1998 , 2000 , 2002 ; frick et al . 2000 , 2001 ; lou et al . + 0.2 cm 4 . angular momentum transfer in steady logarithmic spirals for coplanar mhd perturbation configurations 0.2 cm + stationary logarithmic spiral configurations of mhd density waves carry constant angular momentum flux either outward or inward associated with the advective transport , the gravity torque and the magnetic torque ( lynden - bell & kalnajs 1972 ; goldreich & tremaine 1978 ; fan & lou 1999 ) . therefore , the net angular momentum flux in the entire disc system is conserved and there must be either a source or a sink at the disc center . for the @xmath177 branch with @xmath86 , we find that the total angular momentum flux is inward for leading ( i.e. , @xmath223 ) spiral mhd density waves and is outward for trailing ( i.e. , @xmath224 ) spiral mhd density waves . for the @xmath178 branch in comparison , the total angular momentum flux is outward for leading ( i.e. , @xmath223 ) spiral mhd density waves and is inward for trailing ( i.e. , @xmath224 ) spiral mhd density waves . as a real mhd disc system typically rotates at a relatively large @xmath87 and inequality ( [ restrt ] ) can also rule out small @xmath87 for @xmath360 , the lower @xmath178 branch might be rarely useful in galactic applications . in either case , however , the net mass accretion rate is zero since if we focus on materials between two concentric circles , at any time the influent angular momentum from the inner ring is equal to the effluent angular momentum from the outer ring such that materials inside this belt neither gain nor lose any angular momentum and therefore no net mass accretion occurs . this research has been supported in part by the asci center for astrophysical thermonuclear flashes at the university of chicago under department of energy contract b341495 , by the special funds for major state basic science research projects of china , by the tsinghua center for astrophysics , by the collaborative research fund from the national natural science foundation of china ( nsfc ) for young outstanding overseas chinese scholars ( nsfc 10028306 ) at the national astronomical observatory , chinese academy of sciences , by nsfc grant 10373009 at the tsinghua university , and by the yangtze endowment from the ministry of education through the tsinghua university . affiliated institutions of y.q.l . share this contribution .    9    balbus s. a. , hawley j. f. , 1998 , rev . , 70 , 1 beck r. , brandenburg a. , moss d. , shukurov a. , sokoloff d. , 1996 , ara&a , 34 , 155    beck r. , hoernes p. , 1996 , nat , 379 , 47    binney j. , tremaine s. , 1987 , galactic dynamics . princeton university press , princeton , new jersey    caunt s. e. , tagger m. , 2001 , a&a , 367 , 1095 chandrasekhar s. , 1961 , hydrodynamic and hydromagnetic stability . dover , new york    elmegreen b. g. , 1995 , mnras , 275 , 944    evans n. w. , read j. c. a. , 1998a , mnras , 300 , 83    evans n. w. , read j. c. a. , 1998b , mnras , 300 , 106    fan z. h. , lou y .- q . , 1996 , nat , 383 , 800    fan z. h. , lou y .- q . , 1999 , mnras , 307 , 645 ferrire k. m. , 2001 , rev . mod . phys . , 73 , 1031 frick p. , beck r. , shukurov a. , sokoloff d. , ehle m. , kamphuis j. , 2000 , mnras , 318 , 925 frick p. , beck r. , berkhuijsen e. m. , patrickeyev i. , 2001 , mnras , 327 , 1145 goodman j. , evans n. w. , 1999 , mnras , 309 , 599 goldreich p. , lynden - bell d. , 1965 , mnras , 130 , 125    goldreich p. , tremaine s. , 1978 , icarus , 34 , 240 jog c. j. , 1996 , mnras , 278 , 209    julian w. h. , toomre a. , 1966 , apj , 146 , 810    kalnajs a. j. , 1971 , apj , 166 , 275    kim w .- , ostriker e. c. , 2000 , apj , 540 , 372 kronberg p. p. , 1994 , reports on progress in physics , 57 , 325 lemos j. p. s. , kalnajs a. j. , lynden - bell d. , 1991 , apj , 375 , 484 lin c. c. , shu f. h. , 1964 , apj , 140 , 646    lou y .- q . , 2002 , mnras , 337 , 225    lou y .- q . , fan z. h. , 1998a , apj , 493 , 102    lou y .- q . , fan z. h. , 1998b , mnras , 297 , 84 lou y .- q . , fan z. h. , 2000 , mnras , 315 , 646 lou y .- q . , fan z. h. , 2002 , mnras , 329 , l62 lou y.q . , walsh w. m. , han j. l. , fan z. h. , 2002 , apj , 567 , 289 lou y. q. , yuan c. , fan z. h. , 2001 , apj , 552 , 189 lou y. q. , yuan c. , fan z. h. , leon s. , 2001 , apj , 553 , l35 lou y .- q . , shen y. , 2003 , mnras , 343 , 750 ( astro - ph/0304270 )    lou y .- q . , wu y. , 2004 , mnras , submitted lou y .- q . , zou y. , 2004a , mnras , 350 , 1220 ( astro - ph/0312082 ) lou y .- q . , zou y. , 2004b , mnras , submitted low b. c. , lou y .- q . , 1990 , apj , 352 , 343 lynden - bell d. , kalnajs a. j. , 1972 , mnras , 157 , 1 lynden - bell d. , lemos j. p. s. , 1993 , preprint ( astro - ph/9907093 )    mathewson d. s. , van der kruit p. c. , brouw w. n. , 1972 , a&a , 17,468 mestel l. , 1963 , mnras , 126 , 553 neininger n. , 1992 , a&a , 263 , 30 qian e. , 1992 , mnras , 257 , 581 rdiger g. , kitchatinov l. l. , 2000 , astro . nachr . , 3 , 181 safronov v. s. , 1960 , ann . , 23 , 979    shakura n. i. , sunyaev r. a. , 1973 , a&a , 24 , 337 shen y. , lou y .- q . , 2003 , mnras , 345 , 1340 ( astro - ph/0308063 )    shen y. , lou y .- q . , 2004a , chjaa , in press ( astro - ph/0404190 )    shen y. , lou y .- q . , 2004b , mnras , 353 , 249 ( astro - ph/0405444 )    shu f. h. , laughlin g. , lizano s. , galli d. , 2000 , apj , 535 , 190 shu f. h. , li z .- y . , 1997 , apj , 475 , 251    sofue y. , fujimoto m. , wielebinski r. , 1986 , ara&a , 24 , 459 spruit h. c. , 1987 , a&a , 184 , 173 syer d. , tremaine s. , 1996 , mnras , 281 , 925 tagger m. , pellat r. , 1999 , a&a , 349 , 1003 toomre a. , 1964 , apj , 139 , 1217 toomre a. , 1977 , ara&a , 15 , 437    visser h. c. d. , 1980a , a&a , 88 , 149    visser h. c. d. , 1980b , a&a , 88 , 159    widrow l. m. , 2002 , rev . phys . , 74 , 775 zang t. a. , 1976 , phd thesis , mit , cambridge ma for @xmath127 , we first examine the solution criterion of quadratic equation ( [ aligned ] ) for the aligned case . the determinant @xmath184 of quadratic equation ( [ aligned ] ) is @xmath361 where the three coefficients @xmath362 ( @xmath363 ) are defined explicitly by @xmath364 ^ 2\\ & \qquad\qquad\qquad\qquad\qquad + \frac{(1 - 4\beta)}{(1 + 2\beta ) } \bigg[\frac{-2m^2+(1 - 4\beta)(1 - 2\beta)}{(1 + 2\beta)}{f } { \cal c}{\cal p}_m+1 - 4\beta\bigg]{\cal h}_m\ , \\ c_1\equiv\   & 2(1-{f}{\cal c}{\cal p}_m)\bigg\{{\cal a}_m \bigg[\frac{(3 - 4\beta)m^2-(1 - 4\beta)(1 - 6\beta+4\beta^2)}{2(1 + 2\beta ) } { f}{\cal c}{\cal p}_m+m^2-(1 - 4\beta)(2-\beta)\bigg]\\ & \qquad\qquad\qquad\qquad\qquad + \frac{-2m^2+(1 - 4\beta)(1 - 2\beta)}{(1 + 2\beta)}{\cal h}_m\bigg\}\ , \\ c_0\equiv\   & { \cal a}_m^2(1-{f}{\cal c}{\cal p}_m)^2\ , \end{split}\ ] ] respectively . for the quadratic expression ( [ a01 ] ) in terms of @xmath81 , we further introduce another determinant @xmath365 such that @xmath366 [ ( 1 + 8\beta)m^2-(1 - 4\beta)(1 + 4\beta-8\beta^2)]\ .\ ] ] the determinant @xmath184 of quadratic equation ( [ aligned ] ) in terms of @xmath87 can be either positive or negative in general for different parameter combinations of @xmath8 , @xmath109 , @xmath77 and @xmath283 . it is fairly straightforward to obtain @xmath87 solutions by assigning values of all relevant parameters . in particular , for one important case of @xmath19 ( i.e. , the msid case with a flat rotation curve and an isothermal equation of state ) , the condition that @xmath172 for @xmath127 and @xmath367 when @xmath188 holds for all values of @xmath81 in a full msid with @xmath344 and the condition that @xmath352 for @xmath127 holds in a partial msid with @xmath343 ( lou 2002 ) .    for the magnetic force - free case of @xmath11 with a declining rotation curve with increasing @xmath14 and a polytropic index @xmath182 , we have @xmath352 for @xmath127 and an arbitrary @xmath368 in both full and partial magnetized discs . we now proceed to prove this statement rigorously . for @xmath11 of the aligned case , we explicitly have @xmath369 ^ 2\ , \qquad c_1\equiv 2(1-{f}{\cal c}{\cal p}_m)m^2 [ { \cal a}_m(2{f}{\cal c}{\cal p}_m/3 + 1)-4{\cal h}_m/3]\ , \qquad c_0\equiv { \cal a}_m^2(1-{f}{\cal c}{\cal p}_m)^2\ , \ ] ] where @xmath209 and @xmath370 for @xmath11 and @xmath365 introduced above simply reduces to @xmath371 it is apparent that for @xmath372 , we have @xmath373 . as @xmath374 , it follows immediately that @xmath352 in definition ( [ a01 ] ) . on the other hand for @xmath375 , we have @xmath376 and therefore @xmath352 by definition ( [ a01 ] ) .    for @xmath11 in the logarithmic spiral case , we again have the determinant @xmath352 for quadratic equation ( [ spiral ] ) , where in parallel with definitions ( [ a01 ] ) and ( [ a02 ] ) of the aligned case , the relevant parameters in the spiral case are @xmath377 ^ 2\ , \ c_1=2(1-{f}{\cal c } { \cal n}_m)(m^2+\xi^2)\bigg[{\cal a}_m\bigg(\frac{2}{3}{f } { \cal c}{\cal n}_m+1\bigg)-\frac{4}{3}{\cal h}_m\bigg]\ , \ c_0={\cal a}_m^2(1-{f}{\cal c}{\cal n}_m)^2\ , \ ] ] with @xmath378 and @xmath370 for @xmath11 . and @xmath365 simply reduces to @xmath379 hence , as the aligned @xmath11 case , it is apparent that for @xmath372 , we have @xmath373 . as @xmath374 , it follows immediately that @xmath352 . on the other hand for @xmath380 , we have @xmath376 and therefore @xmath352 . in this appendix , we first provide a proof that for the aligned case of @xmath11 , the physical @xmath87 solutions of the plus @xmath357 and the minus @xmath381 determined from quadratic equation ( [ q2aligned ] ) correspond to the factor @xmath382 being positive and negative , respectively . according to appendix a , the determinant @xmath184 for both aligned and spiral cases remains always positive for the case of @xmath11 . we start from equation ( [ q2aligned ] ) by casting it in the following compact form of @xmath383 where @xmath166 and in reference to equation ( [ q2aligned ] ) , @xmath384 , @xmath385 and @xmath386 are constants defined explicitly by @xmath387 } \ .\ ] ] it is apparent that @xmath385 is always positive . while @xmath384 and @xmath386 remain always positive for @xmath186 , their signs depend on the value of @xmath283 for the special case of @xmath188 . in fact , both @xmath384 and @xmath386 diverge at a critical value of @xmath283 for the @xmath188 case as already discussed in subsection 3.1.3 . let us first analyze the simpler cases of @xmath186 . as @xmath388 and @xmath389 , we have two @xmath87 solutions @xmath390^{1/2}}{2 } & > \frac{\wp_2+\wp_3+|\wp_2-\wp_3|}{2}\ge \wp_2\ ; \\ y_2=\frac{\wp_1+\wp_3-[(\wp_1+\wp_3)^2 - 4\wp_2\wp_3]^{1/2}}{2 } & = \frac{\wp_1+\wp_3-[(\wp_1+\wp_3 - 2\wp_2)^2 + 4\wp_2(\wp_1-\wp_2)]^{1/2}}{2}\\ & < \frac{\wp_1+\wp_3-|\wp_1+\wp_3 - 2\wp_2|}{2}\le \wp_2\ . \end{split}\ ] ] we next examine the @xmath188 case . the situation becomes somewhat involved since @xmath391 , @xmath392 when @xmath291 where the critical @xmath283 value @xmath393\simeq 0.6842 $ ] is determined from the condition @xmath292 ( see subsection 3.1.3 ) . on the other hand , we have @xmath394 when @xmath395 in general . for the aligned case and @xmath178 are consistent with the general solution form ( [ solalign ] ) for the aligned case . ] and for @xmath291 , we therefore have two @xmath87 solutions satisfying inequalities @xmath396^{1/2}}{2}<0\ ; \\ 0<y_2&=\frac{\wp_1+\wp_3+[(\wp_1+\wp_3)^2 - 4\wp_2\wp_3]^{1/2}}{2 } \\ & = \frac{\wp_1+\wp_3+[(\wp_1+\wp_3 - 2\wp_2)^2   + 4\wp_2(\wp_1-\wp_2)]^{1/2}}{2}\\ & < \frac{\wp_1+\wp_3+|\wp_1+\wp_3 - 2\wp_2|}{2}=\wp_2\ , \end{split}\ ] ] and for @xmath395 the situation remains the same as ( [ b2 ] ) .    for the spiral case of @xmath11 in parallel , we can repeat this same procedure of analysis . as before , we cast quadratic equation ( [ beta0d25s ] ) in the form of @xmath397 } \\ & \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad = \frac{d^2\ { d^2-(1-{f}{\cal c}{\cal n}_m)/ [ { f}{\cal c}{\cal n}_m+3(m^2 - 3/2)/(2m^2 + 2\xi^2 + 3/2)]\ } } { d^2-(1-{f}{\cal c}{\cal n}_m)/ ( { f}{\cal c}{\cal n}_m+3/2)}\ , \end{split}\ ] ] which yields two distinct @xmath87 solutions according to the analysis in subsection 3.2.3 . the following treatment is identical to that for the aligned @xmath11 case described in the first part of this appendix . we here provide a proof to show that the angular momentum flux associated with stationary coplanar mhd perturbations due to the gravity torque @xmath327 , the advective transport @xmath325 and the magnetic torque @xmath329 remains constant at all radii for the entire range @xmath10 . we start from expression ( [ rmphi ] ) for @xmath398 to derive the following two expressions in a straightforward manner : @xmath399j_m(k^{\prime}r ) \int_0^{\infty}r_a^{\prime\ -1/2}\cos(m\theta-\xi\ln r_a^{\prime})j_m(k^{\prime}r_a^{\prime})dr_a^{\prime}\ . \end{split}\ ] ] then in integral ( [ gravangm ] ) , we first integrate the kernel over @xmath400 to yield @xmath401dkdk^{\prime}\ , \ ] ] where @xmath402dr_adr_a^{\prime}\ .\ ] ] we next integrate over the vertical coordinate @xmath403 to yield the final result @xmath404dkdk^{\prime}\ .\ ] ] by introducing the following integral transformation @xmath405 we reduce expression ( [ final ] ) to @xmath406dxdx^{\prime}\ , \ ] ] where @xmath407dr_bdr_b^{\prime}={\cal n}_mx^{-1/2}x^{\prime\ -1/2}\sin\bigg[\xi\ln\bigg(\frac{x}{x^{\prime}}\bigg)\bigg]\ , \ ] ] where @xmath297 is the kalnajs function ( kalnajs 1971 ) .    with these technical preparations , we finally determine the angular momentum flux associated with coplanar stationary mhd perturbations due to the gravity torque as @xmath408 where @xmath409dxdx^{\prime}\ , \ ] ] which is clearly independent of radius @xmath14 . in the tight - winding regime or the wkbj limit ( i.e. , @xmath410 ) , we may directly apply the result of fan & lou ( 1999 , see their eqn . 4.13 ) by noting @xmath411 such that @xmath412 where @xmath413 has been defined by equation ( [ thetas ] ) ; here , @xmath414 is inward for leading ( @xmath223 ) spiral waves and outward for trailing ( @xmath224 ) spiral waves . in general , @xmath415 can be calculated via numerical integrations and the specific result is displayed in fig . c1 . for fairly open mhd density waves ( @xmath416 ) with @xmath417 , the angular momentum flux associated with the gravity torque is therefore inward for leading ( @xmath223 ) spiral waves and outward for trailing ( @xmath224 ) spiral waves ; with increasing @xmath418 , @xmath415 falls off rapidly and oscillates around zero with the envelope profile closely fitting the wkbj result better and better .     versus @xmath219 , shown in heavy solid line . we plot @xmath419 in light dashed line for comparison . for increasing @xmath418 , @xmath415 oscillates around zero and its envelope profile fits the @xmath419 curve with an increasing accuracy . ] meanwhile , we compute the other two perturbation contributions to the angular momentum flux , namely , those from the advective transfer and the magnetic torque . for the part of advective transport , we start from equation ( [ adveangm ] ) and make use of phase relationship ( [ relsurjzspiral ] ) to explicitly derive @xmath420r^{-1/2}\cos(\xi\ln r - m\theta)\\ & \qquad\qquad\qquad\qquad\qquad\qquad\qquad\ -(\omega\sigma/\sigma_0)[(1/2+\beta)\im(\theta_s ) -\xi\re(\theta_s)]r^{-1/2}\sin(\xi\ln r - m\theta)\ . \end{split}\ ] ] it immediately follows that the advective torque is given by @xmath421=m\pi^2g\sigma^2\xi^a\ , \ ] ] where @xmath422 } { { f}{\cal c}[1+d^2-(1 - 4\beta)q^2/(2 + 4\beta)]}\ , \ ] ] that is independent of radius @xmath14 . by specifying parameters @xmath8 , @xmath109 , @xmath283 , @xmath77 and @xmath219 , one finds that for the @xmath177 branch @xmath423 for @xmath223 and @xmath424 for @xmath224 . therefore , the advective angular momentum flux is radially inward for leading stationary spiral density waves and radially outward for trailing stationary spiral density waves for the @xmath177 branch of @xmath87 solutions . for the @xmath178 branch in comparison , the situation switches and the advective angular momentum flux is outward for leading stationary spiral density waves and inward for trailing stationary spiral density waves . practically , to achieve the @xmath178 branch of stationary logarithmic spiral configurations requires a much lower @xmath87 which would be unusual for real magnetized disc systems . an example of illustration for @xmath190 , @xmath344 and @xmath11 is shown in panel ( a ) of fig . c2 .    for the momentum flux associated with the magnetic torque , we start from equation ( [ magnangm ] ) and make use of phase relationship ( [ relsurjzspiral ] ) to first derive @xmath425 \ r^{-3/2}\sin(\xi\ln r - m\theta ) + [ mb_0\sigma\im(\theta_s)/\sigma_0]\ r^{-3/2 } \cos(\xi\ln r - m\theta)\ , \\ & \re(b_\theta)=[\xi b_0\sigma\re(\theta_s ) /\sigma_0]\ r^{-3/2}\sin(\xi\ln r - m\theta ) + [ \xi b_0\sigma\im(\theta_s)/\sigma_0 ] \ r^{-3/2}\cos(\xi\ln r - m\theta)\ . ] it immediately follows that the momentum flux @xmath329 associated with the magnetic torque is given by @xmath426 where @xmath427 } \ , \ ] ] which is again independent of radius @xmath14 ; here , @xmath329 is inward for leading spiral density waves ( i.e. , @xmath223 ) and outward for trailing spiral density waves ( i.e. , @xmath224 ) . an example of illustration for @xmath190 , @xmath344 and @xmath11 is shown in panel ( b ) of fig . c2 .    by comparing the three contributions to the total angular momentum flux transport , we reach the following conclusions . first , the contribution from the advective transport part tends to be more dominant over those from the magnetic torque and the gravity torque . secondly , for the @xmath177 branch for stationary fast mhd density waves that has a counterpart in the unmagnetized case , leading spiral density waves ( @xmath223 ) transport angular momentum inward and trailing spiral density waves ( @xmath224 ) transport angular momentum outward . thirdly , for the @xmath178 branch for stationary slow mhd density waves caused by the very presence of the azimuthal magnetic field , leading spiral density waves ( @xmath223 ) carry angular momentum outward and trailing spiral density waves ( @xmath224 ) carry angular momentum inward . note that all these results are only valid for stationary logarithmic spiral mhd density waves . there might be some reasons to rule out the @xmath178 branch in galactic contexts . for example , the @xmath178 branch is of lower @xmath87 which can not be sustained in real galactic discs . furthermore , for rotation curve parameter @xmath360 , the @xmath178 branch may be ruled out by inequality ( [ restrt ] ) .
Abstract:

Answer: we construct analytically stationary global configurations for both aligned and logarithmic spiral coplanar magnetohydrodynamic ( mhd ) perturbations in an axisymmetric background mhd disc with a power - law surface mass density @xmath0 , a coplanar azimuthal magnetic field @xmath1 , a consistent self - gravity and a power - law rotation curve @xmath2 where @xmath3 is the linear azimuthal gas rotation speed . the barotropic equation of state @xmath4 is adopted for both mhd background equilibrium and coplanar mhd perturbations where @xmath5 is the vertically integrated pressure and @xmath6 is the barotropic index . for a scale - free background mhd equilibrium , a relation exists among @xmath7 , @xmath8 , @xmath9 and @xmath6 such that only one parameter ( e.g. , @xmath8 ) is independent . for a linear axisymmetric stability analysis , we provide global criteria in various parameter regimes . for nonaxisymmetric aligned and logarithmic spiral cases , two branches of perturbation modes ( i.e. , fast and slow mhd density waves ) can be derived once @xmath8 is specified . to complement the magnetized singular isothermal disc ( msid ) analysis of lou , we extend the analysis to a wider range of @xmath10 . as an example of illustration , we discuss specifically the @xmath11 case when the background magnetic field is force - free . angular momentum conservation for coplanar mhd perturbations and other relevant aspects of our approach are discussed . mhd  ism : magnetic fields  stars : formation  galaxies : kinematics and dynamics  galaxies : spiral  galaxies : structure .