dynamics of self - gravitating stellar systems and plasma fluids are governed by the collisionless boltzmann equation ( cbe ) @xcite .
finding a general solution of the cbe has been a challenging problem in various disciplines of physical sciences . due to existing difficulties of the general problem , finding a solution to
the linearized cbe became the center of attraction in the twentieth century when @xcite and van kampen ( 1955 ) discovered the normal modes of collisionless ensembles . later in 1970
s , @xcite developed a matrix theory that was capable of computing normal modes of stellar systems through solving a nonlinear eigenvalue problem .
his theory remained as the only analytical perturbation theory used by the community of galactic dynamicists over the past three decades .
@xcite assumed an exponential form @xmath0 for the time - dependent part of physical quantities where @xmath1 , and solved the linearized cbe for the perturbed distribution function ( df ) @xmath2 in terms of the perturbed potential @xmath3 . after expanding the potential and density functions in terms of bi - orthogonal basis sets in the configuration space , he used the weighted residual form of the fundamental equation @xmath4 to obtain a nonlinear eigenvalue problem for the complex eigenfrequency @xmath5 .
for self - consistent perturbations the surface density @xmath6 is related to @xmath3 through poisson s integral , and the symbols @xmath7 and @xmath8 denote the elements of velocity and position vectors .
@xcite used kalnajs s theory to compute the modes of the isothermal disk @xcite , which has the astrophysically important property of a flat rotation curve .
his analysis was then extended by @xcite to general scale - free disks with arbitrary cusp slopes .
application of kalnajs s theory to soft - centered models of stellar disks has been mainly focused on the isochrone and kuzmin - toomre disks @xcite . a disk with exponential light profile and an approximately flat rotation curve was also investigated by @xcite .
more recently jalali & hunter ( 2005a , hereafter jh ) gave new results for soft - centered models of stellar disks .
they showed the importance of a boundary integral in the modal properties of unidirectional disks and computed a fundamental bar mode and a secondary spiral mode for the isochrone , kuzmin - toomre and a newly introduced family of cored exponential disks .
jh also extended kalnajs s first order perturbation theory to the second order , and illustrated energy and angular momentum content of different fourier components .
their bar charts showed that only a few number of expansion terms in the radial angle govern the perturbed dynamics .
implementation of kalnajs s ( 1977 ) theory , however , has some technical problems due to the nonlinear dependency of his matrix equations on @xmath5 .
most computational methods that deal with nonlinear eigenvalue problems are iterative .
they start with an initial guess of the solution and continue with a search scheme in the frequency space .
newton s method is perhaps the most efficient technique that guarantees a quadratic convergence should the initial guess be close enough to the solution .
the key issue in success of any iterative scheme is the attracting or repelling nature of an eigenvalue .
it is obvious that only attracting eigenvalues can be captured by iterative methods while we have no priori knowledge of their basins of attraction in order to make our initial guess .
the mentioned computational difficulties make it a formidable task to explore all normal modes of a stellar system , which include growing modes as well as stationary van kampen modes .
moreover , it is not easy to develop a general nonlinear theory based on kalnajs s method for studying the interaction of modes .
@xcite introduced an alternative method for the normal mode calculation of stellar disks whose outcome was a linear eigenvalue problem for @xmath5 .
his method is capable of finding all eigenmodes of a stellar disk should one use fine grids in the action space .
polyachenko s method is somehow costly because it results in a large linear system of equations to assure point - wise convergence in the action space and mean convergence of fourier expansions in the space of the radial angle .
extension of his method to nonlinear regime is another challenging problem yet to be investigated .
@xcite has also followed an approach similar to @xcite and studied the instability of stellar disks surrounding massive objects .
his eigenvalue equations involve action variables , and practically , need to be solved over a discretized grid in the space of actions .
recent developments in fluid mechanics @xcite inspired me to formulate the dynamics of stellar systems in a new framework , which is capable of solving the cbe not only in the linear regime , but also in its full nonlinear form .
the method systematically searches for smooth solutions of the cbe by expanding the perturbed df using fourier series of angle variables and an appropriate set of trial functions in the space of actions .
coefficients of expansion are unknown time - dependent amplitude functions whose evolution equations are obtained by the petrov - galerkin projection @xcite of the cbe .
that is indeed taking the weighted residual form of the cbe by integration over the action - angle space and deriving a system of nonlinear ordinary differential equations ( odes ) for the amplitude functions .
the associated first order system of odes leads to a linear eigenvalue problem , which is solved using standard numerical methods .    in this paper
i present my new method and apply it to explore the modal properties of a model galaxy . in a second paper , i will address the nonlinear evolution of modes and wave interactions .
the paper is organized as follows . in sections [ sec::dynamical - theory ] and [ sec::linear - theory ] ,
i use the petrov - galerkin method to project the cbe to a system of odes in the time domain and derive a system of linear eigenvalue equations .
section [ sec : modes - exp - disk ] presents the eigenfrequency spectra and their corresponding mode shapes of the cored exponential disk of jh .
the stars of this model move in the field of a soft - centered logarithmic potential .
i also investigate the effect of physical parameters of the equilibrium model on the modal content . in section [ sec : final - discussions ] , i discuss on the nature of a self - gravitating mode and compare the performance of my method with other theories . some fundamental achievements of this work are summarized in section [ sec : conclusions ] .
i use the usual polar coordinates @xmath9 and assume that the temporal evolution of the df and gravitational potential starts from an axisymmetric equilibrium state described by @xmath10 and @xmath11 so that @xmath12 here @xmath13 is the momentum vector conjugate to @xmath14 .
motion of stars in the equilibrium state is governed by the zeroth order hamiltonian @xmath15 for bounded orbits , @xmath16 and @xmath17 become librating and rotating , respectively .
one can therefore describe the dynamics using the action variables @xmath18 , @xmath19 and their conjugate angles @xmath20 .
a transformation @xmath21 leaves the hamiltonian @xmath22 as a function of actions only , @xmath23 , and therefore , the phase space flows of the equilibrium state lie on a two dimensional torus @xmath24 with @xmath25 being a constant 2-vector .
an action - angle transformation can locally be found for any bounded regular orbit , but it is a global transformation if only one orbit family occupies the phase space .
the axisymmetric potential @xmath11 supports only rosette orbits .
radial and circular orbits are the limiting cases of rosette orbits with @xmath26 and @xmath27 , respectively . by representing @xmath28 in terms of the action - angle variables ,
the cbe reads @xmath29=0 , \label{eq : cbe - action - angle}\ ] ] where @xmath30 $ ] denotes a poisson bracket taken over the action - angle space .
according to jeans theorem @xcite @xmath31 depends on the phase space coordinates through the integrals of motion , which are the actions in the present formulation , and one obtains @xmath32=0 $ ] . subsequently , equation ( [ eq : cbe - action - angle ] )
may be rewritten as @xmath33- \left [ f_0,{\cal h}_1 \right ] -\left [ f_1,{\cal h}_1 \right ] , \label{eq : cbe - perturbed}\ ] ] where @xmath34 is the perturbed hamiltonian .
a dark matter halo contributes both to @xmath22 and to @xmath34 if it is live , i.e. , if it exchanges momentum / energy with the luminous stellar component . in this paper
i confine myself to a rigid halo that only contributes to @xmath22 through @xmath35 and assume that @xmath36 is the perturbed potential due to self - gravity .
let me expand @xmath2 and @xmath3 in fourier series of angle variables and write @xmath37 where @xmath38 and @xmath39 are some trial functions in the space of action variables , and @xmath40 and @xmath41 are time - dependent amplitude functions . on the other hand , one can expand @xmath3 and its corresponding surface density @xmath42 in the configuration space as @xmath43 here @xmath44 and @xmath45 are bi - orthogonal potential  surface density pairs that satisfy the relation @xmath46 where @xmath47 is the kronecker delta and @xmath48 are some constants .
it is remarked that the real parts of @xmath42 , @xmath3 and @xmath2 describe physical solutions . in this work
i utilize the @xcite functions @xmath49 that yield @xcite @xmath50 @xmath51 are associated legendre functions with @xmath52 .
clutton - brock functions have a length scale @xmath53 , which makes them suitable for reproducing the potential and surface density of soft - centered models .
the choice of this parameter is an important step in the calculation of normal modes .
i will discuss on this issue later in
[ sec : modes - exp - disk ] .    equating ( [ eq : expansion - v1 ] ) and ( [ eq : expansion - v1-config ] ) , multiplying both sides of the resulting equation by @xmath54 $ ] and integrating over the @xmath55-space , lead to ( see also kalnajs 1977 and tremaine & weinberg 1984 ) @xmath56 d \theta _
r ,   \label{eq : fourier - coeffs}\end{aligned}\ ] ] where @xmath57 are the fourier coefficients of the basis potential functions in the configuration space .
the trial functions used in the expansion of @xmath3 in the action - angle space are not necessarily identical to @xmath57 .
however , subsequent mathematical derivations are greatly simplified by setting @xmath58 , which implies @xmath59 . to build a relation between @xmath40 and @xmath60
, i use the fundamental equation @xmath61 where @xmath62 is the volume of an infinitesimal phase space element . on substituting ( [ eq : expansion - f1 ] ) and ( [ eq : expansion - sigma1-config ] ) in ( [ eq : fundamental - equation ] )
, multiplying both sides of the resulting equation by @xmath63 and integrating , one obtains @xmath64 which is inserted in ( [ eq : expansion - v1 ] ) to represent @xmath3 in terms of the amplitude functions @xmath40 as @xmath65 it would be computationally favorable to collect @xmath40 in a single vector @xmath66 by defining a map @xmath67 . in practice
the infinite sums in ( [ eq : expansion - f1 ] ) are truncated and approximated by finite sums so that @xmath68 , @xmath69 and @xmath70 . for @xmath71 ,
a simple map between indices will be @xmath72 one can now use ( [ eq : expansion - f1 ] ) and ( [ eq : expansion - v1-vs - dmlj ] ) in ( [ eq : cbe - perturbed ] ) and apply the petrov - galerkin method to construct the weighted residual form of the cbe .
that is to multiply ( [ eq : cbe - perturbed ] ) by some weighting functions @xmath73 and to integrate the identity over the action - angle space .
the outcome is the following system of nonlinear odes @xmath74 for the amplitude functions @xmath75 .
the elements of @xmath76 and @xmath77 have been determined in appendix a. each equation in ( [ eq : nonlinear - ode ] ) is the projection of the cbe on a subspace spanned by a weighting function . therefore , the left hand side of ( [ eq : nonlinear - ode ] ) is the projection of @xmath78 , the summation over first order terms is the projection of @xmath79- \left [ f_0,{\cal h}_1 \right ] $ ] , and the second order terms are the projections of @xmath80 $ ] .
the second order terms of amplitude functions , characterized by @xmath77 , show the interaction of modes in both the radial and azimuthal directions .
distribution of angular momentum between different fourier components provides useful information of the disk dynamics .
i compute the rate of change of the total angular momentum @xmath81 using ( see appendix b in jh ) @xmath82 where a bar denotes complex conjugate .
substituting ( [ eq : expansion - f1 ] ) and ( [ eq : expansion - v1-vs - dmlj ] ) in ( [ eq : rate - change - l - one ] ) and evaluating the integrals , yield @xmath83 \lambda^{ml}_{jp } d^{ml}_{p}(t ) \nonumber \\ & { } & - m\left [ a^{m}_{j}(t ) + \overline{a^{(-m)}_{j}}(t ) \right ] \lambda^{ml}_{jp } \overline{d^{ml}_{p}}(t ) \biggr \}. \label{eq : rate - change - l - two}\end{aligned}\ ] ] define @xmath84 with @xmath85 and @xmath86 being real functions of time .
according to identity ( [ eq : a - versus - d ] ) , one may further simplify equation ( [ eq : rate - change - l - two ] ) to @xmath87 .
\label{eq : rate - change - l - four}\end{aligned}\ ] ] the share of the @xmath88th mode from @xmath89 is thus determined by @xmath90 .
as one could anticipate for an isolated stellar disk , @xmath89 vanishes and the total angular momentum remains constant because the terms @xmath90 and @xmath91 cancel each other in ( [ eq : rate - change - l - three ] ) , and @xmath92 is annulled by the factor @xmath88 in ( [ eq : rate - change - l - four ] ) .
choosing the trial functions @xmath38 is the most delicate step in the reduction of the cbe to a system of odes .
one possible way is to set @xmath93 in ( [ eq : cbe - perturbed ] ) and solve the first order equation @xmath94+\left [ f_0,{\cal h}_1 \right ] = 0 , \label{eq : cbe - linear - equilibrium}\ ] ] for @xmath2 .
substituting ( [ eq : expansion - f1 ] ) and ( [ eq : expansion - v1 ] ) in ( [ eq : cbe - linear - equilibrium ] ) gives @xmath95 where @xmath96 equation ( [ eq : trial - functions - one ] ) suggests to choose @xmath97 as the trial functions ( in the space of actions ) for an unsteady @xmath98 .
these functions have integrable singularities for resonant orbits with @xmath99 . for unidirectional disks with only prograde orbits
, they also include a term with the dirac delta function @xmath100 ( see jh ) .
one should therefore avoid the partial derivatives of @xmath101 with respect to the actions by evaluating the weighted residual form of @xmath102 $ ] through integration by parts ( appendix a ) .    for deriving the relation between @xmath103 and @xmath104 in ( [ eq : a - versus - d ] ) , the fundamental equation ( [ eq : fundamental - equation ] )
was multiplied by the complex conjugates of the basis functions used in the expansion of @xmath105 .
one may follow a similar approach for obtaining the weighted residual form of the cbe and set @xmath106 which are the complex conjugates of the basis functions used in the expansion of @xmath107 in the action - angle space .
the trial and weighting functions introduced as above , are not orthogonal but they result in a simple form for the linear part of the reduced cbe as i explain in  [ sec::linear - theory ] .
in a first order perturbation analysis , the second order terms of the amplitude functions are ignored .
the evolution of modes is then governed by the linear parts of ( [ eq : ode - for - z - one ] ) as @xmath108 a general solution of ( [ eq : ode - for - z - linear ] ) has the form @xmath109 , which leads to the following linear eigenvalue problem @xmath110 for a prescribed azimuthal wavenumber @xmath88 .
operating @xmath111 on ( [ eq : eigenvalue - problem - one ] ) yields @xmath112 where @xmath113 is a general non - symmetric matrix .
a reduction to hessenberg form followed by the qr algorithm @xcite gives all real and complex eigenvalues .
real eigenvalues correspond to van kampen modes and complex eigenvalues , which occur in conjugate pairs , give growing / damping modes .
i utilize the method of singular value decomposition for finding the eigenvectors and perform the decomposition @xmath114 where @xmath115 is the identity matrix and the diagonal matrix @xmath116 is composed of the singular values @xmath117 ( @xmath118 ) .
the column of @xmath119 that corresponds to the smallest @xmath117 is the eigenvector associated with @xmath5 .    calculation of @xmath120 and @xmath121 involves evaluation of some definite integrals in the action space .
there will be two types of such integrals ( instead of three ) if one uses the trial functions defined in ( [ eq : trial - functions - two ] ) .
let me introduce the auxiliary integral @xmath122 and apply the trial functions @xmath123 in ( [ eq : a - tensor ] ) .
the elements of @xmath121 and @xmath120 are thus computed from @xmath124 { \cal i}^{ml}_{jk}\lambda^{ml'}_{kj'}. \label{eq : bar - a - matrix}\end{aligned}\ ] ] both @xmath125 and @xmath126 consist of boundary integrals when the unperturbed stellar disk is unidirectional with the df @xmath127 . here
@xmath128 is the heaviside function .
the boundary terms are @xmath129_{j_{\phi}=0 } , \\ % { \tilde { \cal i}}^{ml}_{jk } & = & \int_{0}^{\infty } d j_r\left [ mf^p_0(\textbf{\textit{j } } ) \psi^{ml}_{j}(\textbf{\textit{j}})\psi^{ml}_{k}(\textbf{\textit{j } } ) \right ] _ { j_{\phi}=0}.\end{aligned}\ ] ] dynamics of modes with different azimuthal wavenumbers are decoupled in the linear regime and the matrix @xmath113 is an odd function of the wavenumber @xmath88 .
i.e. , @xmath130 .
an immediate result of this property is @xmath131 .
consequently , @xmath90 becomes equal to zero for all @xmath132 and each mode individually conserves the total angular momentum .
the present theory has three major advantages over kalnajs s formulation .
firstly , all eigenmodes relevant to a prescribed azimuthal wavenumber are obtained at once with classical linear algebraic algorithms .
this makes it possible to explore and classify all families of growing modes beside pure oscillatory van kampen modes .
secondly , the constituting integrals of the elements of @xmath133 , @xmath134 and @xmath135 ( appendix [ app::petrov - galerkin ] ) are regular at exact resonances when the condition @xmath136 holds .
finally , nonlinear interaction of modes , and the mass and angular momentum exchange between them , can be readily monitored by integrating the system of nonlinear odes given in ( [ eq : nonlinear - ode ] ) . in the proceeding section
i will be concerned with the calculation and classification of modes in the linear regime .
jh calculated barred and spiral modes of certain stellar disks for the wavenumber @xmath137 . among the models studied in jh ,
the cored exponential disk with the surface density profile @xmath138 and embedded in the field of the soft - centered logarithmic potential @xmath139 is a viable model that resembles most features of realistic spirals . here
@xmath140 is the core radius , @xmath141 is the length scale of the exponential decay , and @xmath142 is a density scaling factor .
the velocity of circular orbits in this model rises from zero at the galactic center and approaches to the constant value @xmath143 in outer regions where the light profile falls off exponentially .
@xcite have derived the gravitational potential corresponding to @xmath144 .
i denote this potential by @xmath145 . the gradient @xmath146,\ ] ] will give the gravitational force of a spherical dark matter component , computed inside the galactic disk .
the density profile of the dark component , @xmath147 , can then be determined using @xmath148 .
the positiveness of @xmath147 imposes some restrictions on the physical values of @xmath149 and @xmath150 as figure 5 in jh shows . for a given @xmath149
, @xmath151 can not exceed a critical value @xmath152 .
the parameter @xmath149 determines the shape of the dark matter density profile . a model with @xmath153 and @xmath154
is maximal in the region where the rotation curve is rising .
i.e. , there is no dark matter in that region .
models with @xmath155 and @xmath154 are still maximal but only in the vicinity of the center for @xmath156 . in such models
the rotational velocity of stars due to dark matter ( @xmath157 ) has a monotonically rising profile . for @xmath158 ,
dark matter penetrates into the galactic center and its density profile becomes cuspy in the limit of @xmath159 .
the role of the parameter @xmath151 is to control the fraction of dark to luminous matter .
models with @xmath160 are dominated by dark matter .
jh introduced a family of equilibrium dfs that reproduces @xmath144 and depends on an integer constant @xmath161 .
this parameter controls the population of near - circular orbits and the disk temperature : the parameter @xmath162 of @xcite decreases by increasing @xmath161 .
the dfs of jh have an isotropic part that determines the fraction of radial orbits .
that isotropic part , which reconstructs the central density of the equilibrium state , shrinks to central regions of the galaxy as @xmath161 increases .
i apply my new method to the cored exponential disks of jh and calculate the spectrum of @xmath163 .
subsequently , the eigenvector @xmath164 is calculated from ( [ eq : eigenvalue - problem - two ] ) and it is used in ( [ eq : a - versus - d ] ) to compute @xmath165 and the perturbed density @xmath166 , \label{eq : mode - shape}\ ] ] which is the real part of ( [ eq : expansion - sigma1-config ] ) .
@xmath167 and @xmath168 are the amplitude and phase functions of an @xmath88-fold circumferential wave that travels with the angular velocity @xmath169 .
the factor @xmath170 shows the exponential growth / decay of the wave amplitude .
i normalize all length , velocity , and time variables to @xmath140 , @xmath143 and @xmath171 , respectively , and set @xmath172 .
i begin my case studies in  [ sec : maximal - disk ] with a near maximal disk of @xmath173=@xmath174 and compute its eigenfrequency spectra for the wavenumbers @xmath175 .
i then classify unstable @xmath137 modes of this model and investigate their evolution as the parameter @xmath151 is varied . in
[ sec : variations - of - lambda ] and  [ sec : variation - of - n ] , i study the behavior of unstable @xmath137 waves as the parameters @xmath149 and @xmath161 are changed .
the eigenfrequency spectrum of a model with an inner cutout is also computed and discussed in
[ sec : cutout - models ] .
i pick up the first model from table 4 of jh with @xmath176 and start solving the eigensystem ( [ eq : eigenvalue - problem - two ] ) with @xmath177 and increase these limits until complex eigenfrequencies converge . for an error threshold of @xmath178 the program terminates when @xmath179 , which gives a size of @xmath180 for the matrix @xmath113 . in such a circumstance , out of @xmath181 eigenfrequencies of @xmath113 ( for each wavenumber @xmath88 ) , less than 15 pair have non - zero growth rates ( @xmath182 ) . further increasing of @xmath183 and @xmath184
does not alter the number and location of complex eigenfrequencies in the @xmath5-plane .
this shows that unstable modes do not constitute a continuous family .
figure [ pic : full - spectra ] displays the eigenfrequency spectra for the azimuthal wavenumbers @xmath175 .
eigenfrequencies on the real axis are oscillatory van kampen modes .
their calculation requires evaluation of cauchy s principal value @xcite if one uses kalnajs s first order theory . in the present formalism
, van kampen modes are found together with growing modes without any special treatment .
more van kampen modes are obtainable by increasing the truncation limit @xmath183 of fourier terms in the @xmath185-direction .
toomre s @xmath162 is marginally greater than 1 for the model ( figure 7_b _ in jh ) , and therefore , one could expect that the disk is stable for @xmath186 excitations ( see top - left panel in figure [ pic : full - spectra ] ) .
the model is highly unstable for @xmath187 excitations although the average growth rate of unstable modes decreases for larger wavenumbers .
it is evident that either unstable modes are isolated or they are grouped in _
discrete families_. depending on the wavenumber , there may be one or more discrete families .
the most prominent family bifurcates from van kampen modes .
members of this family have spiral patterns with multiple peaks in their @xmath167 functions .
the ( global ) fastest growing mode belongs to the spectrum of @xmath137 .
that is the bar mode of a two - member unstable family .
the length scale of clutton - brock functions has been set to @xmath188 for @xmath137 and @xmath189 for other wavenumbers .
changing this length scale slightly displaces the eigenfrequencies although the spectrum maintains its global pattern .
large values of @xmath53 lead to a better computation accuracy of extensive modes ( with smaller pattern speeds ) , while compact bar modes show a rapid convergence for small values of @xmath53 .
moreover , the suitable value of @xmath53 differs from one azimuthal wavenumber to another .
finding an optimum length scale that gives the best results for all modes and wavenumbers is an open problem yet to be investigated precisely .
for the cored exponential disks with @xmath190 , working in the range @xmath191 gives reasonable results .
for @xmath137 , i have zoomed out and plotted in figure [ pic : spectrum - m2-unstable]_a _ the portion of the spectrum that contains growing modes .
the first and second modes reported in table 4 of jh have been shown by circles in the same figure
. the most unstable mode ( labeled as b1 ) is a compact , rapidly rotating bar .
the majority of unstable modes belong to a discrete spiral family that bifurcates from a van kampen mode with @xmath192 .
i have labeled these modes by s1,@xmath193,s6 .
the number of density peaks along the spiral arms is proportional to the integer number in the mode name .
both b2 and g are double peaked spirals but i have classified b2 as a bar mode , and collected it with b1 in a two member family , for it takes a bar - like structure when it is stabilized by decreasing @xmath151 .
i classify mode g as an isolated mode because it does not behave similar to either of s- or b - modes as the model parameters vary .
there is another isolated mode in the spectrum , c1 , which exhibits a spiral pattern . by decreasing @xmath149 ,
mode c1 joins a new family of spiral modes , which are accumulated near the galactic center ( see  [ sec : variations - of - lambda ] ) .    reducing @xmath151 increases the abundance of dark matter and according to @xcite and jh the growth rate of modes should decrease .
my calculations show that by reducing @xmath151 , spiral modes are affected sooner and more effective than the bar mode , and they join to the stationary modes , one by one from the location of the bifurcation point until the whole s - family disappears .
this is a generic scenario for all @xmath194 models regardless of the disk temperature controlled by @xmath161 .
solid lines in figure [ pic : spectrum - m2-unstable]_b _ show the _ eigenfrequency loci _ of a model with @xmath195 as @xmath151 increases from @xmath196 to @xmath197 .
it is evident that mode b1 is destabilized through a pitchfork bifurcation while the loci of modes b2 and c1 , and s - modes exhibit a tangent bifurcation .
all modes except mode g are stable for @xmath198 .
surprisingly , mode g resists against stabilization even for very small values of @xmath151 .
this indicates that mode g is not characterized by the fraction of dark to luminous matter . in
[ sec : variations - of - lambda ] and  [ sec : variation - of - n ] , i will show that this mode is highly sensitive to the variations of @xmath149 and @xmath161 . according to my computations ( e.g. , figure [ pic : spectrum - m2-unstable]_b _ ) , by increasing @xmath151 all s - modes
are born at the same bifurcation frequency @xmath199 , but mode b2 comes out from a van kampen mode with @xmath200 .
this result completely rules out any skepticism that mode b2 is a member of s - family .
figure [ pic : mode - shape - m2 ] displays the wave patterns and amplitude functions of modes b1 , b2 , g , s1 , s2 and s3 .
it is seen that the patterns of s - modes rotate slower and become more extensive as the mode number increases .
mode s6 , which is at the bifurcation point of the spiral family , has the largest extent .
eight wave packets of this mode are distributed by a phase shift of 90 degrees along major spiral arms .
mode g has at most two density peaks on its major spiral arms but the magnitude of its second peak increases as the disk is cooled .
modes b1 , b2 , s1 , s2 and s3 are , respectively , analogous to modes a , b , c , e and f of a gaussian disk explored in @xcite .
there are three low - speed modes in figure 11 of @xcite that have not been labeled , but they are analogs of modes s4 , s5 and s6 .
mode g and toomre s mode d also have some similarities but are of different origins ( see  [ sec : final - discussions ] ) .
none of them can be stabilized only by increasing the fraction of dark matter .
figure [ pic : spectrum - m2-unstable]_a _ shows that the fundamental mode obtained by jh coincides with mode b1 .
the wave pattern of mode b1 ( displayed in figure [ pic : mode - shape - m2 ] ) is identical to the mode shape computed using kalnajs s theory and demonstrated in figure 8 of jh .
jh found a secondary mode which lies between modes b2 and g. that mode is also a double - peaked spiral and it is not easy to identify its true nature unless we investigate its evolution as the model parameters vary . by comparing figure 10_a _ of jh with figure [ pic : spectrum - m2-unstable]_b _
, one can see that both mode b2 and the secondary mode of jh are destabilized through a tangent bifurcation while mode g has a different nature . the bifurcation frequency of mode b2 that i find ( @xmath200 ) matches very well with the frequency of the stabilized secondary mode of jh ( see figure 10_a _ in jh but note that their vertical axis indicates @xmath201 ) .
i conclude that the secondary mode of jh is indeed mode b2 although it seems to be closer to mode g. the existing discrepancy is due to different length scale of clutton - brock functions that jh have used for finding the secondary mode . by adjusting
@xmath53 one can improve the location of b2 .
however , this is an unnecessary attempt given the fact that mode b2 has already been identified , and the computation accuracy of other eigenfrequencies has an impressive level .
the parameter @xmath149 controls the density profile of the dark matter component , specifically near the galactic center .
the fraction of dark to luminous matter has its minimum value in marginal models with @xmath202 .
i choose a marginal @xmath203 model with @xmath204 , which has also been investigated by jh .
figure [ pic : spectrum - m2-vary - lambda]_a _ shows the portion of the spectrum that contains complex eigenfrequencies of this model .
the spectrum has been computed for @xmath188 .
although @xmath137 bar and spiral modes survive in this model , their pattern speeds and growth rates drop considerably .
mode g has been wiped out of existence by dark matter penetration into the center , and four unstable modes ( c1 , c2 , c3 and c4 ) have emerged that constitute a new family of spiral modes . they populate the central regions of the disk in most of @xmath203 models .
again , the location of eigenfrequencies obtained by jh have been marked by circles .
the agreement between the results of jh and the present study is very good and the variance is less than @xmath205 .
the population of spiral modes is changed by varying @xmath151 , and the eigenfrequencies of unstable modes are altered significantly .
solid lines in figure [ pic : spectrum - m2-vary - lambda]_a _ show the loci of growing modes as @xmath151 increases from @xmath196 to @xmath206 .
similar to the previous @xmath153 model , s - modes and mode b1 are destabilized through tangent and pitchfork bifurcations , respectively .
all c - modes are born by a pitchfork bifurcation although some minor modes of the same nature come and go as @xmath151 varies .
the loci of modes b1 and s1 ( in figure [ pic : spectrum - m2-vary - lambda]_a _ ) are in harmony with the results of kalnajs s method plotted in figure 10_b _ of jh .
it is noted that the locus of mode b1 steeply joins the real axis , well before stabilizing the s - modes .
this is how slowly growing spirals may dominate a stellar disk .
there are no new families of growing modes in @xmath207 models .
dark matter in these models induces a rising rotation curve on the disk stars ( see figure 6 in jh ) and the population of s - modes declines .
the growth rate of mode b2 increases proportional to @xmath149 , but that of mode g falls off although mode g is still robust against the variations of @xmath151 .
mode g has its maximum growth rate in @xmath153 models , which suggests that it must be a self - gravitating response of the luminous matter that involves only the potential of the disk , @xmath208 .
the function @xmath209 of mode c1 loses its minor peaks and becomes smoother as @xmath149 increases .
figure [ pic : spectrum - m2-vary - lambda]_b _ shows the eigenfrequency loci of models with @xmath210 as @xmath151 increases from @xmath211 to @xmath212 .
the eigenfrequency loci of modes c1 and b2 , and s - modes ( as @xmath151 varies ) are similar to @xmath153 models , but the locus of mode b1 loses its steepness and stretches towards small pattern speeds in an approximately linear form until it joins the real @xmath213-axis .
the bifurcation frequency of mode b2 differs from s - modes and it is larger .    figure [ pic : mode - shape - m2-sr34 ] shows the wave patterns of modes b1 , b2 , s1 , c1 , c2 and c3 for the model with @xmath204 . the ( isolated ) mode b1 is still a single - peaked bar although its edge is more extensive as the flat part of its @xmath209 plot indicates
. mode s1 is a triple - peaked spiral ( as before ) and mode b2 is being stabilized ( @xmath214 ) .
it is seen that mode b2 has a bar - like structure , which justifies its classification as the secondary bar mode .
the pattern of s1 and its @xmath209 plot can be compared with figure 9 in jh .
the agreement is quite satisfactory .
there is a remarkable difference between the patterns of c- and s - modes although both families have spiral structures .
in contrast to s - modes that become more extensive as their growth rate decays , c - modes are shrunk to central regions because their pattern speed increases .
c - modes are also a bifurcating family , but their bifurcation point lies on large pattern speeds associated with the azimuthal frequency ( @xmath215 ) of central stars .    the parameters @xmath151 and @xmath149 are essentially controlling the fraction and density profile of dark matter component , respectively . however , the radial velocity dispersion @xmath216 of the equilibrium state is also playing an important role in the perturbed dynamics .
@xmath216 is an indicator of the initial temperature of the disk .
@xcite had already pointed out that the pitch angle of spiral patterns decreases when @xmath217 ( see their figure 7 ) . apart from this morphological implication , can the variation in the disk temperature affect the modal content ? in the following subsection ,
i trace the evolution of growing waves by changing the disk temperature and show that the population of s - modes is larger in rotationally supported , cold disks .
the parameter @xmath161 of the dfs of jh controls the disk temperature by adjusting the size of the isotropic core and the population of near circular orbits . as @xmath161 increases
, the streaming velocity @xmath218 approaches the rotational velocity of circular orbits and the stellar disk is cooled .
figure [ pic : spectrum - vary - n ] displays the eigenfrequencies of previous @xmath219 and @xmath220 models for @xmath221 and @xmath222 .
the spectra for the intermediate value of @xmath223 have already been shown in figures [ pic : spectrum - m2-unstable]_a _ and [ pic : spectrum - m2-vary - lambda]_a_.    increasing @xmath161 gives birth to more s - modes while the bifurcation point of the family is preserved . as a new member
is born at the bifurcation point , other members including mode s1 , are pushed away from the real axis on a curved path .
this behavior is observed in both models but the branch of s - family in the model with @xmath220 stays closer to the real axis than the other model .
the growth rates of c - modes increase remarkably as the disk is cooled . despite mode b1 which rotates and grows faster in cold disks ,
mode b2 grows faster in warmer disks .
variation in the disk temperature changes the eigenfrequency of mode g more effective than what @xmath151 could , but nothing is more influential than the role of @xmath149 .
another consequence of cooling the stellar disk is that @xmath186 waves are no longer stable .
the parameter @xmath162 of toomre was marginally greater than unity for @xmath223 models . for @xmath222 ,
i find @xmath224 over an annular region because the plot of @xmath162 versus @xmath16 exhibits a minimum at some finite radius ( e.g. , figure 7 of jh ) .
for instance , i find three growing @xmath186 modes for the model @xmath225 .
they correspond to pure complex eigenfrequencies @xmath226 , @xmath227 and @xmath228 .
mode shapes have ( obviously ) ringed structures but the number of rings , which is identical to the number of peaks of @xmath229 , depends on the growth rate .
the modes associated with @xmath230 , @xmath231 and @xmath232 have three , four and five rings , respectively .
ring modes are very sensitive to the variations of model parameters and they are suppressed by decreasing @xmath149 and @xmath151 .      in order to simulate an immobile bulge , which does not respond to density perturbations , jh utilized an inner cutout function of the form @xmath233 where @xmath234 is an angular momentum scale . multiplying @xmath235 by the self - consistent df of the equilibrium state ,
prohibits the stars with @xmath236 from participating in the perturbed dynamics .
consequently , incoming waves are reflected at some finite radius and the innermost wave packets of multiple - peaked modes are diminished .
my calculations show that all s - modes survive in cutout models , mode g disappears , and the growth rate of mode b2 increases .
the pattern speed of mode b1 is boosted so that the corotation resonance is destroyed , but its growth rate drops drastically .
figure [ pic : spectrum - cutout ] shows the eigenfrequency spectrum of a model with @xmath237 and @xmath176 .
circles show the eigenfrequencies found by jh .
again , the agreement between the results of jh and the present work is very good .
the reason that i have identified mode b2 as the second member of b - family , and not the most unstable s - mode , is that its @xmath209 function has an evolved double - peaked structure ( see figure 11 in jh ) and its locus versus @xmath151 does not emerge from the same bifurcation frequency of s - modes .
disappearance of mode g in cutout models confirms my earlier note that it is a self - gravitating mode .
there are similarities between mode g of this study and toomre s ( 1981 ) mode d. both of these modes resist against stabilization by increasing the fraction of dark to luminous matter and they have at most double peaks on their spiral arms
. nonetheless , these modes are not the same because mode g is amplified through a feedback from the galactic center but toomre s mode d has been identified as an edge mode .
a question remains to be answered : why toomre ( 1981 ) did not detect mode g and i do not find an edge mode ?
the most convincing explanation is that to excite a self - gravitating wave inside the core of the stellar component , the governing potential in that region should mainly come from the self - gravity of stars . this requirement is fulfilled in my @xmath238 models for @xmath156 .
however , the completely flat rotation curve imposed by toomre ( 1981 ) nowhere follows the rotational velocity induced by the self - gravity of stars and it prohibits the gaussian disk from developing a g - like mode . on the other hand , i do nt find an edge mode because the density profile of the cored exponential disk does not decay as steep as the gaussian disk to create an outer boundary at some finite radius for reflecting the outgoing waves .
similar to the first order analysis of  [ sec::linear - theory ] , polyachenko s ( 2005 ) approach results in the full spectrum of eigenfrequencies for a given azimuthal wavenumber .
there are some differences between his method and the present formulation .
polyachenko directly uses poisson s integral to establish a point - wise relation in the action space between the fourier components of the perturbed df and its self - consistent potential .
combination of equations ( 5 ) and ( 9 ) in his paper is analogous to equation ( [ eq : expansion - v1-vs - dmlj ] ) in this paper .
the main departure of the two theories is in the way that the linearized cbe is treated .
polyachenko forces a point - wise fulfillment of the cbe in the action space while the present method works with a weighted residual form of the cbe .
a point - wise formulation poses a challenge for the numerical calculation of the eigenvalues and their conjugate eigenvectors . according to the bar charts of jh , at least ten fourier components ( @xmath239 )
are needed in the @xmath185-direction to assure a credible convergence of @xmath2 in a typical soft - centered galaxy model .
therefore , if one chooses a grid of @xmath240 in the action space , polyachenko s eigenvector @xmath241 will have a dimension of @xmath242
. therefore , for a _ very coarse _ grid with @xmath243 that polyachenko uses , the unknown eigenvector will have a dimension of 4410 .
this number must be compared with the dimension of @xmath164 in equation ( [ eq : eigenvalue - problem - two ] ) .
that is indeed @xmath244 for the most accurate calculations carried out by setting @xmath179 which means that @xmath245 fourier components in the @xmath185-direction and @xmath246 expansion terms in the @xmath16-direction have been taken into account . noting that the definite integrals @xmath126 and @xmath125 are independently evaluated over the action space with any desired accuracy , the present theory proves to be more efficient for eigenmode calculation ( in the linear regime ) than other existing alternatives .
the agreement between the results of this work and those of jh , who have used kalnajs s method , is impressive .
there is only a discrepancy in the results for a double - peaked spiral mode of @xmath153 models .
in fact , these models have two double - peaked modes , modes b2 and g , and jh find mode b2 .
the origin of discrepancies was attributed to the length scale of clutton - brock functions , @xmath53 , which is a fixed number for the whole spectrum of a given azimuthal wavenumber .
provided that jh optimized @xmath53 for each growing mode that they calculated ( see also  [ sec : maximal - disk ] ) , some minor deviations from the results of this paper are reasonable . in most cases
the algorithm used by jh converges to mode b1 and the fastest rotating s - mode .
they capture mode b2 only if its growth rate is large enough .
other modes remain unexplored because newton s method needs an initial guess of @xmath5 , which has a little chance to be in the basin of attraction of the other members of s - family .
the separation of eigenfrequencies near the bifurcation point of s - modes is very small and one could anticipate complex boundaries for the basins of attraction of these eigenfrequencies . thus , there is no guarantee that successive newton s iterations keep an estimated eigenfrequency on the same basin that it was initially .
nevertheless , in kalnajs s formulation , a systematic search for all growing modes is possible by introducing the mathematical eigenvalue @xcite and investigating its loci as the pattern speed and growth rate vary .
after three decades of kalnajs s ( 1977 ) publication , it was not known exactly whether growing modes of stellar systems appear as distinct roots in the eigenfrequency space or they belong to continuous families as van kampen modes do . in this paper ,
i attempted to answer this question using the galerkin projection of the cbe and unveiled the full eigenfrequency spectrum of a stellar disk .
i showed that similar to gaseous disks @xcite , majority of growing modes emerge as _ discrete families _ through a bifurcation from stationary modes .
there are some exceptions for this rule , the most important of which are the isolated bar and g modes .
the model that i used to test my method allows for dark matter presence as a spherical component , whose potential inside the galactic disk contributes to the rotational velocity of stars . by varying the parameters of the model , and investigating the eigenfrequencies and their associated mode shapes
, i showed that it is not the fraction of dark to luminous matter that controls the variety of growing modes .
what determines that variety is indeed the shape of the dark matter density profile controlled by the parameter @xmath247 .
my survey in the parameter space revealed that the concentration of dark matter in the galactic center ( @xmath203 ) destroys mode g and weakens the growth of b - modes substantially .
emergence of spiral c - modes that accumulate near the galactic center is another remarkable consequence of dark matter presence in central regions of a cored stellar disk .
although the solution of the galerkin system showed a credible convergence of the series expansions , the existence of strong solutions for the cbe , in its full nonlinear form , is still an open problem .
it has been known for years that van kampen modes make a complete set @xcite , and therefore , they may be used for a series representation of stationary oscillations .
but there is not a mathematical proof for the completeness of the discrete families of growing modes . in other words ,
whether an observed galaxy can be assembled using the modes of a linear eigensystem , requires further analysis .    in the second part of this study
, i will investigate the mechanisms of wave interactions in the nonlinear regime and will probe the mass and angular momentum transfer between waves of different fourier numbers .
i am indebted to chris hunter for his instructive and valuable comments since the beginning of this work .
i also thank the referee for helpful suggestions that improved the presentation of the results .
this work was partially supported by the research vice - presidency at sharif university of technology .
let me define a nonlinear operator @xmath248 and denote @xmath249 as the @xmath250th prolongation @xcite of the physical quantity @xmath251 in the domain of independent variables .
assume a ( nonlinear ) partial differential equation @xmath252 and its associated initial and boundary conditions that govern the evolution of @xmath253 in the domain of the spatial variable @xmath254 and the time @xmath255 . a weighted residual method @xcite attempts to find an approximate solution of the form @xmath256 through determining the time - dependent functions @xmath257 for a given set of trial ( basis ) functions @xmath258 .
the trial functions should satisfy the boundary conditions and be linearly independent . using ( [ eq : expansion - for - u - general ] ) and taking the inner product of ( [ eq : general - operator - prolonged ] ) by some weighting functions @xmath259 , yield the determining equations of @xmath257 as @xmath260
there are several procedures for choosing @xmath259 and each procedure has its own name . the method with @xmath261
is called the bubnov - galerkin , or simply the galerkin method .
the petrov - galerkin method is associated with @xmath262 .
the well - known collocation method uses dirac s delta functions for the weighting purpose .
there is an alternative interpretation for the inner product @xmath263 .
that is projecting the equation @xmath264 on a subspace spanned by the weighting function @xmath265
. therefore , equation ( [ eq : weighted - residual - form - general ] ) is often called the _ galerkin projection _ of ( [ eq : general - operator - prolonged ] ) . in what follows ,
i use the petrov - galerkin method and construct the weighted residual form of the cbe .
assume the functions @xmath266 and @xmath267 , and define their inner product over the action - angle space as @xmath268 taking the inner product of the perturbed cbe by the weighting functions @xmath269 gives @xmath270,w^{ml}_{j}\right ) -\left
( \left [ f_0,{\cal h}_1 \right ] , w^{ml}_{j } \right ) - \left ( \left [ f_1,{\cal h}_1 \right ] , w^{ml}_{j}\right )
. \label{eq : cbe - weighted - residual-1}\ ] ] note that the cbe is the governing equation of the perturbed df whose trial functions are @xmath38 . with my choice of the weighting function ( as above )
i am following the petrov - galerkin method . on substituting ( [ eq : expansion - f1 ] ) and ( [ eq : expansion - v1-vs - dmlj ] ) in ( [ eq : cbe - weighted - residual-1 ] ) and after some rearrangements of summations , one obtains @xmath271 \lambda^{m'l'}_{kj ' } \int   d \textbf{\textit{j } } \left ( l{\partial f_0\over \partial j_r } + m'{\partial f_0\over \partial j_{\phi } } \right ) \psi^{m'l}_{j}(\textbf{\textit{j}})\psi^{m'l}_k(\textbf{\textit{j } } ) \nonumber \\ % & + & \sum_{m',l'}^{}\sum_{j'}^ { }      \sum_{m'',l''}^{}\sum_{j''}^ { } \delta_{m'',(m - m ' ) } d^{m'l'}_{j'}(t ) d^{m''l''}_{j''}(t ) \sum_{k}^ { } \left [ { 4\pi^2\over d_k(m '' ) } \right ] \lambda ^{m''l''}_{kj '' } \nonumber \\ % & \times & \biggl [ \int   d \textbf{\textit{j } } \psi^{ml}_{j}(\textbf{\textit{j}})\phi^{m'l'}_{j'}(\textbf{\textit{j } } ) \left ( l'{\partial \over \partial j_r}+          m'{\partial \over \partial j_\phi } \right ) \psi^{m''(l - l')}_{k}(\textbf{\textit{j } } ) \nonumber \\ % & { } & \qquad -\int   d \textbf{\textit{j } } \psi^{ml}_{j}(\textbf{\textit{j } } ) \psi^{m''(l - l')}_{k}(\textbf{\textit{j } } ) \left ( ( l - l'){\partial \over \partial j_r } +          m''{\partial \over \partial j_\phi } \right ) \phi^{m'l'}_{j'}(\textbf{\textit{j } } ) \biggr ] , \label{eq : weighted - residual - expan } % \end{aligned}\ ] ] where @xmath272 , @xmath273 and @xmath274 . using equation ( [ eq : map - mlj - to - i ] ) and carrying out the index mappings @xmath275 , @xmath276 and @xmath277 one
may introduce the arrays @xmath278 \lambda^{m'l'}_{kj ' } \int   d \textbf{\textit{j } } \left ( l{\partial f_0\over \partial j_r } + m'{\partial f_0\over \partial j_{\phi } } \right ) \psi^{m'l}_{j}(\textbf{\textit{j } } ) \psi^{m'l}_k(\textbf{\textit{j}})\biggr ] , \label{eq : a - tensor } \\
% k_{pqr}&= & \delta_{m',(m - m '' ) } \delta_{l',(l - l '' ) } \sum_{k=0}^{j_{\rm max } } \left [ { 4\pi^2\over d_k(m '' ) } \right ] \lambda ^{m''l''}_{kj '' } \nonumber \\ & { } & \times \biggl [ \int   d \textbf{\textit{j } } \psi^{ml}_{j}(\textbf{\textit{j } } ) \phi^{m'l'}_{j'}(\textbf{\textit{j } } ) \left ( l'{\partial \over \partial j_r}+          m'{\partial \over \partial j_\phi } \right ) \psi^{m''l''}_{k}(\textbf{\textit{j } } ) \nonumber \\ % & { } & \qquad -\int   d \textbf{\textit{j } } \psi^{ml}_{j}(\textbf{\textit{j } } ) \psi^{m''l''}_{k}(\textbf{\textit{j } } ) \left ( l''{\partial \over \partial j_r } +          m''{\partial \over \partial j_\phi } \right ) \phi^{m'l'}_{j'}(\textbf{\textit{j } } ) \biggr ] .\label{eq : b - tensor } % \end{aligned}\ ] ] consequently , equation ( [ eq : weighted - residual - expan ] ) takes the following matrix form @xmath279 let the matrix @xmath280 $ ] be the inverse of @xmath281 $ ] and left - multiply ( [ eq : ode - for - z - one ] ) by @xmath111 to get @xmath282 evaluation of the integrands in ( [ eq : b - tensor ] ) will be considerably simplified if one avoids the partial derivatives of @xmath38 through integrating ( [ eq : b - tensor ] ) by parts . that gives @xmath283 \lambda ^{m''l''}_{kj '' } \nonumber
\times \biggl \ { & { } & \int { \rm d}\textbf{\textit{j } } \psi^{ml}_{j}(\textbf{\textit{j } } ) \phi^{m'l'}_{j'}(\textbf{\textit{j } } ) \left ( l{\partial \over \partial j_r}+          m{\partial \over \partial j_\phi } \right ) \psi^{m''l''}_{k}(\textbf{\textit{j } } ) \nonumber \\ % & + & \int   d \textbf{\textit{j } } \phi^{m'l'}_{j'}(\textbf{\textit{j } } ) \psi^{m''l''}_{k}(\textbf{\textit{j } } ) \left ( l''{\partial \over \partial j_r } +          m''{\partial \over \partial j_\phi }
\right ) \psi^{ml}_{j}(\textbf{\textit{j } } ) \biggr \}.\label{eq : b - tensor - by - parts}\end{aligned}\ ] ] when all stars move on prograde orbits , the equilibrium df takes the form @xmath127 where @xmath128 is the heaviside function . upon using ( [ eq : trial - functions - two ] )
, this contributes a term including dirac s delta function @xmath100 to the trial functions .
thus , the following boundary terms @xmath284 \lambda ^{m''l''}_{kj '' } \nonumber \\ % \times \biggl \ { \int_{0}^{\infty } \!\ ! &   d j_r & \left [ { m ' f^p_0(\textbf{\textit{j}})\psi^{ml}_{j}(\textbf{\textit{j } } )                            \psi^{m'l'}_{j'}(\textbf{\textit{j } } ) \over           l'\omega_r(\textbf{\textit{j}})+m'\omega_\phi(\textbf{\textit{j } } )          } \left ( l{\partial \over \partial j_r}+          m{\partial \over \partial j_\phi } \right ) \psi^{m''l''}_{k}(\textbf{\textit{j } } ) \right ] _ { j_{\phi}=0 } \nonumber \\ % + \int_{0}^{\infty } \!\ ! &
d j_r & \left [          { m ' f^p_0(\textbf{\textit{j } } ) \psi^{m''l''}_{k}(\textbf{\textit{j } } )                               \psi^{m'l'}_{j'}(\textbf{\textit{j } } )     \over l'\omega_r(\textbf{\textit{j}})+m'\omega_\phi(\textbf{\textit{j } } )          } \left ( l''{\partial \over \partial j_r } +          m''{\partial \over \partial j_\phi } \right ) \psi^{ml}_{j}(\textbf{\textit{j } } ) \right ] _ { j_\phi=0 }   \biggr \ } , \label{eq : b - tensor - boundary}\end{aligned}\ ] ] must be added to @xmath135 when the equilibrium disk is unidirectional .
the partial derivatives of @xmath39 needed for equations ( [ eq : b - tensor - by - parts ] ) and ( [ eq : b - tensor - boundary ] ) are calculated by differentiating equation ( [ eq : fourier - coeffs ] ) partially with respect to an action : @xmath285 \nonumber \\ & { } & \qquad -m\psi^{|m|}_j(r){\partial \over \partial j_\nu }       \left(\theta_{\phi}-\phi\right )       \sin [ l\theta _ r+m(\theta _ \phi -\phi ) ] \biggr \ }   d \theta _ r,~~\nu \equiv r,\phi .
\label{eq::deriv - fourier - coeff}\end{aligned}\ ] ] jalali & hunter ( 2005b ) encountered these partial derivatives in their second order perturbation theory devised for computing the energy of eigenmodes .
i adopt their technique for calculating the quantities @xmath286 and @xmath287 .
the variables @xmath16 , @xmath288 , and @xmath289 are regarded as functions of @xmath290 because the action - angle transformation @xmath291 is defined in the phase space of an axisymmetric state . from @xmath292 one
may write @xmath293= { \partial^2 r \over \partial \theta_r \partial j_\nu}\frac{d\theta_r}{dt } = \omega_r\frac{\partial}{\partial j_\nu } \left({\partial r \over \partial \theta_r}\right ) = \omega_r\frac{\partial}{\partial j_\nu}\left(\frac{v_r}{\omega_r}\right ) = { \partial v_r \over \partial j_\nu}-\frac{v_r}{\omega_r } { \partial \omega_r \over \partial j_\nu}.\ ] ] similarly , one obtains @xmath294 & = & { \partial \omega_{\phi } \over \partial j_\nu } -\frac{\delta_{\nu,\phi}}{r^2 } + { 2j_{\phi } \over r^3}{\partial r \over \partial j_\nu } -\frac{1}{\omega_r}\left[\omega_{\phi}-\frac{j_{\phi}}{r^2}\right ] { \partial \omega_r \over \partial j_\nu } , \label{eq::dthetadj } \\ \frac { d } { d t}\left[{\partial v_r \over \partial j_\nu}\right ] & = & { 2j_{\phi } \over r^3}\delta_{\nu,\phi } -\left[{3j_{\phi}^2 \over r^4}+v_0^{\prime\prime}(r)\right ] { \partial r \over \partial j_\nu } -\frac{1}{\omega_r}\left[{j_{\phi}^2 \over r^3}-v_0^{\prime}(r)\right ] { \partial \omega_r \over \partial j_\nu}. \label{eq::dvdj}\end{aligned}\ ] ] the set of three equations ( [ eq::drdj ] ) through ( [ eq::dvdj ] ) can be integrated along an orbit , and they provide the additional values needed to evaluate the partial derivatives ( [ eq::deriv - fourier - coeff ] ) .
initial values are @xmath295 at @xmath296 where @xmath297 because @xmath298 for all orbits .
however the initial @xmath299 values change with the actions , and initial values for the derivatives of @xmath16 with respect to the actions are @xmath300_{r = r_{\rm min}}= \frac{r^3_{\rm min}\omega_r } { r^3_{\rm min}v_0^{\prime}(r_{\rm min})-j^2_{\phi } } , ~~ \left [ { \partial r \over \partial j_{\phi}}\right]_{r = r_{\rm min}}= \frac{r_{\rm min}(r^2_{\rm min}\omega_{\phi}-j_{\phi } ) } { r^3_{\rm min}v_0^{\prime}(r_{\rm min})-j^2_{\phi}}.\ ] ] they are obtained by differentiating the zeroth order energy equation .