single particle tracking microscopy provides the position time series @xmath0 of individual particle trajectories in a medium @xcite . the information garnered by single particle tracking yields insights into the mechanisms and forces , that drive or constrain the motion of the particle .
an early example of systematic single particle tracking is given by the work of jean perrin on diffusive motion @xcite .
due to the relatively short individual trajectories , perrin used an ensemble average over many trajectories to obtain meaningful statistics .
a few years later , nordlund conceived a method to record much longer time series @xcite , allowing him to evaluate individual trajectories in terms of the time average and thus to avoid averages over not perfectly identical particles . today
, single particle tracking has become a standard tool to characterise the microscopic rheological properties of a medium @xcite , or to probe the active motion of biomolecular motors @xcite . particularly in biological cells and complex fluids
single particle trajectory methods have become instrumental in uncovering deviations from normal brownian motion of passively moving particles .    ) as function of the lag time @xmath1 . left : normal diffusion growing like @xmath2 ( 1 ) , restricted ( confined ) diffusion with a turnover from @xmath3 to @xmath4 ( 2 ) , drift diffusion @xmath5 ( 3 ) .
right : ergodic subdiffusion @xmath6 ( 1 ) , restricted ergodic subdiffusion turning over from @xmath6 to @xmath4 ( 2 ) , non - ergodic subdiffusion @xmath3 ( 1  ) , restricted non - ergodic subdiffusion turning over from @xmath3 to @xmath7 ( 2  ) , superdiffusion @xmath8 ( 3 ) . here , @xmath9 .
note the double - logarithmic scale.,width=332 ]    classical diffusion patterns are sketched in the left panel of fig .
[ diffusion_modes ] .
accordingly , one may observe free diffusion , leading to a linear growth with time of the second moment [ line 1 in fig .
[ diffusion_modes ] ] .
brownian motion may also be restricted ( corralled , confined ) .
confinement in a cell , for instance , could be due to the cell walls . in that case
the second moment initially grows linearly with time and eventually saturates to its thermal value equalling the second moment of the corresponding boltzmann distribution [ line 2 ] . in the presence of a drift
the second moment grows with the square of time [ line 3 ] .
such results are typical for simple fluids . in more complex environments , different patterns
may be observed , as displayed on the right of fig .
[ diffusion_modes ] . here ,
subdiffusion may occur , for which the second moment grows slower than linearly with time [ line 1 ] .
restricted subdiffusion would depart from this behaviour to reach a plateau [ line 2 ] .
driven motion may lead to a superdiffusive power - law form of the second moment with an exponent between 1 and 2 [ line 3 ] .
however , as we demonstrate below , subdiffusion may also be non - ergodic , and the associated time averaged second moment may grow _ linearly _ with time [ line 1  ] .
similarly strange behaviour may be observed for restricted non - ergodic subdiffusion , which exhibits a _
power - law _
growth , not a saturation to a plateau [ line 2  ] .
non - ergodic processes come along with a significant scatter between individual trajectories .
this is an effect of the ageing nature of the process that persists for long measurement times . in the following we discuss in detail the behaviour of passive subdiffusive motion in terms of time and ensemble averages and
address the peculiarities , that may arise for non - ergodic systems .
non - ergodic behaviour of the above sense is indeed observed experimentally .
[ christine ] shows the time averaged mean squared displacement for lipid granules in a living fission yeast cell .
the motion is recorded by indirect tracking in an optical tweezers setup .
initially the granule is located in the bottom of the laser trap potential such that the granule moves freely .
eventually the granule ventures away from the centre of the trap and experiences the hookean trap force .
as demonstrated in a detailed analysis the granule motion indeed exhibits weak ergodicity breaking , giving rise to the characteristic turnover from an initially linear scaling @xmath2 with the lag time @xmath1 , to the power - law regime @xmath10 @xcite .
moreover a pronounced trajectory - to - trajectory scatter is observed , again typical for systems with weak ergodicity breaking .
the initially linear scaling @xmath11 turns over to the power - law regime @xmath12 , induced by the restoring force exerted on the granules by the laser trap @xcite .
note also the characteristic scatter between individual trajectories .
inset : average behaviour of the shown trajectories . in both graphs
the unit of @xmath13 is volt@xmath14 , i.e. , the direct output voltage of the quadrant photodiode .
the voltage is directly proportional to the distance from the centre of the laser trap.,width=332 ]    free diffusion is typically quantified in terms of the second moment .
the mean squared displacement @xmath15 is obtained as the spatial average over the probability density function @xmath16 to find the particle at position @xmath17 at time @xmath18 .
the quantity ( [ eamsd ] ) therefore corresponds to the ensemble averaged second moment of the particle position , denoted by angular brackets , @xmath19 .
in particular , the time @xmath18 enters into eq .
( [ eamsd ] ) only as a parameter .
conversely , single particle trajectories @xmath0 are usually evaluated in terms of the time averaged mean squared displacement defined as @xmath20 where we use an overline @xmath21 to symbolise the time average . here
@xmath1 is the so - called lag time constituting a time window swept along the time series , and @xmath22 is the overall measurement time . the time averaged mean squared displacement thus compares the particle positions along the trajectory as separated by the time difference @xmath1 .    in an ergodic system the time average of a certain quantity obtained from sufficiently long time series is equal to the corresponding ensemble mean @xcite .
for instance , for the mean squared displacement ergodicity would imply @xmath23 brownian motion is ergodic , as well as certain stationary processes leading to anomalous diffusion , such as fractional brownian motion considered below .
there exist , however , non - ergodic processes , which are intimately connected to ageing properties . in what follows we discuss two prominent models for anomalous diffusion , the stationary fractional brownian motion and the ageing continuous time random walk , and analyse in detail the features of the associated time averaged mean squared displacement
.    generally , anomalous diffusion denotes deviations from the classical linear dependence of the mean squared displacement , @xmath24 .
such anomalies include ultraslow diffusion of the form @xmath25 @xcite .
in contrast , anomalous diffusion processes may become faster than ballistic , for instance for systems with correlated jump lengths or in systems governed by generalised langevin equations @xcite .
here we are interested in anomalous diffusion with power - law dependence on time @xcite , @xmath26 for which the anomalous diffusion exponent belongs to the subdiffusive range @xmath9 , such that the limit @xmath27 corresponds to brownian motion .
the proportionality factor @xmath28 in eq .
( [ aneamsd ] ) is the anomalous diffusion coefficient of physical dimension @xmath29 .
the embedding spatial dimension is @xmath30 , and eq .
( [ aneamsd ] ) includes the complete gamma function @xmath31 .
subdiffusion of the form ( [ aneamsd ] ) with @xmath9 occurs in the following biologically relevant systems .
fluorescently labelled mrna in _
e.coli _ bacteria cells was observed to follow @xmath32 with @xmath33 @xcite .
this result is consistent with more recent findings according to which free rna tracers in living cells exhibit @xmath34 , while dna loci show @xmath35 @xcite .
telomeres in the nucleus of mammalian cells were reported to follow anomalous diffusion with @xmath36 at shorter times and value @xmath37 at intermediate times , before a turnover to normal diffusion occurs @xcite .
also larger tracer particles show anomalous diffusion , such as adeno - associated viruses of radius @xmath38 nm in a cell with @xmath39 @xcite and endogenous lipid granules , of typical size of few hundred nm with @xmath40 @xcite .
it should be noted that this subdiffusion observed by single particle tracking microscopy is consistent with results from other techniques , such as fluorescence correlation spectroscopy @xcite or dynamic light scattering @xcite .
subdiffusion of biopolymers larger than some 10 kd in living cells is due to molecular crowding , the excluded volume effect in the superdense cellular environment @xcite .
larger tracer particles also experience subdiffusion due to interaction with the semiflexible cytoskeleton @xcite .
knowledge of the time or ensemble averaged mean squared displacement of an anomalous diffusion process is insufficient to fully characterise the underlying stochastic mechanism , as the associated probability density @xmath16 is no longer necessarily gaussian , and therefore no longer specified by the first and second moments , only @xcite .
this property is in contrast to the universal gaussian nature of brownian motion which is effected by the central limit theorem . at the same time the very nature of the anomalous diffusion process
may result in decisively different behaviours for diffusional mixing and diffusion - limited reactions @xcite .
in biological cells this would imply significant differences for signalling and regulatory processes . for a better understanding of the dynamics in biological cells and other complex fluids knowledge of the underlying stochastic mechanism
is therefore imperative .
here we discuss the properties of continuous time random walk ( ctrw ) processes with diverging characteristic waiting time with respect to their time averaged behaviour , expanding on our earlier work @xcite .
we show that for free motion the lag time dependence of the time averaged mean squared displacement , @xmath41 , is insensitive to the anomalous diffusion exponent @xmath42 for ctrw processes .
in contrast , for confined ctrw subdiffusion a universal scaling behaviour emerges , @xmath43 , with dynamic exponent @xmath44 [ curves ( 1  ) and ( 2  ) in fig .
[ diffusion_modes ] ] .
subdiffusion governed by fractional brownian motion ( fbm ) leads to the scaling @xmath45 of the time averaged mean squared displacement , turning over to a saturation plateau under confinement , @xmath46 [ curves ( 1 ) and ( 2 ) in fig .
[ diffusion_modes ] ] .
we particularly emphasise the irreproducible nature of the time averaged quantities and their associated scatter around the ensemble mean for ctrw subdiffusion processes .
this randomness of the time averages is captured by the distribution function of the amplitude of the time average .
we show that even for relatively short trajectories this distribution is a good characteristic for the underlying process .
in contrast for fbm processes the scatter typical of many single particle experiments is not found in the long time limit .    in the remainder of this work ,
for simplicity we restrict the discussion to the one - dimensional case ( @xmath47 ) .
generalisation to higher dimensions is straightforward .
the article is structured as follows .
we start with a brief introduction to ctrw and fbm .
we then consider the cases of unbounded motion and confined anomalous diffusion in the subsequent two sections .
the distribution of the time averages will be presented thereafter .
finally , we discuss the velocity autocorrelation functions for subdiffusive ctrw and fbm processes , before presenting a concluding discussion .
although both ctrw and fbm give rise to an ensemble averaged mean squared displacement of the form ( [ eamsd ] ) , they are fundamentally different processes , as outlined here .
we note in passing that also in the random motion on a fractal support subdiffusion arises @xcite .
we will not pursue this type of anomalous diffusion in the following .
ctrw theory dates back to the work of montroll and weiss @xcite , and was championed in the analysis of charge carrier motion in amorphous semiconductors by scher and montroll @xcite .
ctrw has become a standard statistical tool to describe processes ranging from particle motion in actin networks @xcite to the tracer motion in groundwater @xcite .
each jump of a ctrw process is characterised by a random jump length and a random waiting time elapsing before the subsequent jump . at each jump
the jump length and waiting time are chosen independently . for a subdiffusive process
we assume that the variance of the jump lengths is given by a finite value @xmath48 , and we consider the unbiased case @xmath49 . on a lattice of spacing @xmath50 , we would have @xmath51 .
in contrast , the waiting times @xmath52 are drawn from the probability density @xmath53 for large @xmath52 , with @xmath9 .
this form of @xmath54 is scale - free , that is , the average waiting time @xmath55 diverges , causing effects such as ageing @xcite and weak ergodicity breaking @xcite . in eq .
( [ wtd ] ) , the quantity @xmath56 is a scaling factor . the anomalous diffusion constant in this case becomes @xcite @xmath57    the scale - freeness of @xmath58 allows individual waiting times @xmath52 to become quite large .
no matter how long a given time series is chosen , single @xmath52 values may become of the order of the length of the entire time span covered in the trajectory .
an example is shown in fig .
[ timeseries ] : the stalling events with large waiting times @xmath52 are quite distinct . for values of the anomalous diffusion exponent @xmath42 that are closer to 1
the stalling is less pronounced .
physically , the power - law form of the waiting time distribution @xmath58 may be related to comb models @xcite or random energy landscapes @xcite ; for more details , compare refs .  .     for brownian motion ( top ) , ctrw ( middle ) , and fbm ( bottom )
the anomalous diffusion exponent for fbm and ctrw is @xmath59 . for ctrw the stalling events
are outstanding , while in the case of fbm strong antipersistence occurs .
note that the stalling events and the antipersistence are less pronounced for larger values of @xmath42 .
typical experimental data include additional noise , such that the appearance of measured data will not display the ideal behaviour shown here.,width=332 ]      fbm is a gaussian process with stationary increments .
its position is defined in terms of the langevin equation @xmath60 or , alternatively , @xmath61 the motion is driven by stationary , fractional gaussian noise @xmath62 with zero mean @xmath63 and long - ranged noise correlation @xcite @xmath64 contrasting the uncorrelated noise for normal diffusion @xmath27 : @xmath65 . in eq .
( [ fbnoise ] ) the anomalous diffusion exponent is connected to the traditionally used hurst exponent by @xmath66 , and we introduce the abbreviation @xmath67 for consistency with standard fbm notation . for subdiffusion the fractional gaussian noise
is anticorrelated , decaying like @xmath68 .
this implies that a given step is likely to go into the direction opposite to the previous step .
the corresponding oscillating behaviour is seen in fig .
[ timeseries ] .
the position autocorrelation function of fbm becomes @xmath69 so that at equal times @xmath70 we recover the mean squared displacement ( [ aneamsd ] ) .    if the fractional gaussian noise is not considered external , but the validity of the fluctuation - dissipation theorem is imposed , one obtains the generalised langevin equation ( gle ) @xcite , @xmath71 where we write @xmath72 instead of the exponent @xmath42 in eq .
( [ fbnoise ] ) . in the gle ( [ gle ] )
, we defined the coupling coefficient @xmath73}$ ] according to the fluctuation - dissipation theorem , where @xmath74 is the boltzmann constant and @xmath75 the absolute temperature . in what follows
we only consider the overdamped limit , in which the inertia term @xmath76 can be neglected .
the gle then gives rise to the form @xmath77 of the mean squared displacement .
in contrast to fbm , that is , the gle leads to subdiffusion for persistent noise with @xmath78 , while @xmath79 yields superdiffusion .
fbm and the related gle are used to describe processes such as long term storage capacity of water reservoirs @xcite , climate fluctuations @xcite , economical market dynamics @xcite , single file diffusion @xcite , and elastic models @xcite .
motion of this type has also been associated with the relative motion of aminoacids in proteins @xcite , and the free diffusion of biopolymers under molecular crowding conditions @xcite .
let us begin with considering anomalous diffusion on an infinite domain and without drift . to find an analytical expression for the time averaged mean squared displacement ( [ tamsd ] ) we note that even a brownian process recorded over a finite time span @xmath22 will show fluctuations in the number of jumps performed during @xmath22 . to average out these trajectory - to - trajectory fluctuations we introduce the ensemble mean , @xmath80 we can then express the integrand in terms of the variance of the jump lengths , @xmath48 , and the average number of jumps @xmath81 in the time interval @xmath82 , as follows : @xmath83 for a regular random walk on average every jump occurs after the waiting time @xmath55 .
thus @xmath84 , and @xmath85 where we defined the diffusion constant @xmath86 $ ] . in the
brownian limit the time averaged mean squared displacement ( [ btamsd ] ) in terms of the lag time @xmath1 takes on exactly the same form as the ensemble averaged mean squared displacement ( [ eamsd ] ) as function of time @xmath18 .
this is not surprising , as brownian motion is an ergodic process . for sufficient duration @xmath22 of the time records any time average converges to the corresponding ensemble average , and thus the ensemble average in expression ( [ eatamsd ] ) is no longer necessary .
the number of jumps of a ctrw process with waiting time distribution of the form ( [ wtd ] ) on average grows sublinearly with time , @xmath87 $ ] @xcite .
this time evolution translates into the mean squared displacement ( [ eamsd ] ) with the anomalous diffusion coefficient ( [ an_coff ] ) . combining the time dependence of @xmath88 for ctrw subdiffusion with the definition ( [ eatamsd ] ) we obtain the following result for the time averaged mean squared displacement , @xmath89 in the limit @xmath90 @xcite .
this follows from expansion of the relation @xmath91 .
the noteworthy feature of this result is that the linear lag time dependence of normal diffusion remains completely unaffected by the anomalous nature of the stochastic process .
only the dependence on the overall measurement time @xmath22 witnesses the underlying subdiffusion .
( [ ctrwtamsd ] ) can be understood if we notice that for normal diffusion we have @xmath92 [ since according to einstein @xmath93 where @xmath94 is the mean time between jump events .
now , for the subdiffusive case @xmath55 diverges and must be replaced by @xmath95 , which explains the term @xmath96 .
[ ctrw_tamsd ] shows 20 trajectories of a ctrw process with @xmath97 , displaying the general trend @xmath98 .    ,
shown for 20 individual trajectories .
the overall measurement time is @xmath99 ( a.u . ) .
note the local changes of the slope in the trajectories as well as the complete stalling in the lowest curve , all bearing witness to the scale - free nature of the underlying waiting time distribution @xmath100.,width=332 ]    rewriting the result ( [ ctrwtamsd ] ) in the form @xmath101 we see that the effective diffusion coefficient @xmath102 decays as function of the measurement time @xmath22 .
the longer the system evolves after its initial preparation , the less mobile it appears , consistent with the ageing property of ctrw subdiffusion .
for instance , in the picture of the trapping events , in the course of time deeper and deeper traps may be encountered by the diffusing particle , such that it gets more and more stuck .
this increasing immobility is mirrored in the form @xmath102 .
note that eqs .
( [ ctrwtamsd ] ) and ( [ rew ] ) suggest that from measuring the @xmath1 dependence of an anomalous diffusion process one might draw the erroneous conclusion that normal diffusion were observed .
the disparity between the time averaged mean squared displacement ( [ ctrwtamsd ] ) and its ensemble averaged counterpart ( [ aneamsd ] ) demonstrates the weak ergodicity breaking characteristic of a process with diverging time scale @xcite .
single or few , long waiting times are also responsible for the pronounced deviations between individual realisations shown in fig .
[ ctrw_tamsd ] .
this apparent irreproducibility of the time averaged mean squared displacement is again intimately coupled to the weak ergodicity breaking nature of the ctrw subdiffusion process .
we will discuss this feature more quantitatively in terms of the distribution @xmath103 as function of the relative deviation @xmath104 of the time averaged mean squared displacement around its ensemble mean in section [ scatter ] .
in contrast to the above behaviour of ctrw subdiffusion processes , fbm is ergodic . indeed , by help of the position autocorrelation ( [ fbm_auto ] ) the time averaged mean squared displacement ( [ eatamsd ] ) becomes @xcite @xmath105 exactly matching the ensemble averaged mean squared displacement , @xmath106 .
however , ergodicity is reached algebraically slowly @xcite , see also below . in fig .
[ fbm_tamsd ] the comparatively minute scatter between individual trajectories supports the ergodic behaviour of fbm processes .
the somewhat larger deviations at longer lag times are due to worsening statistics when @xmath107 .    .
the scatter between individual trajectories is minute , mirroring the ergodicity of this process.,width=332 ]
since the motion of tracer particles is typically confined , for example , by the cell walls or internal membranes , we now consider the important case of anomalous diffusion in a bounded domain .    in a finite interval
@xmath108 $ ] the mean squared displacement of a brownian particle initially released well away from the boundaries will grow linearly in time and eventually turn over to a stationary plateau of magnitude @xmath109 . similarly , if the particle evolves in the confinement of an external harmonic potential of the form @xmath110 , the thermal value @xmath111 $ ] will eventually be attained .
how is this behaviour modified for a ctrw subdiffusion process ? under the influence of an arbitrary external potential @xmath112 defining the force @xmath113 , ctrw subdiffusion can be described by the fractional fokker - planck equation @xcite , or , equivalently , in terms of the following coupled langevin equations @xcite :    @xmath114    here the position @xmath115 is expressed in terms of the parameter @xmath116 ( the internal time ) , and driven by the white gaussian noise @xmath117 .
thus , eq .  ( [ le1 ] ) defines standard brownian motion @xmath118 , where @xmath119 is the diffusivity for the normal diffusion process in internal time @xmath116 .
laboratory time @xmath18 is introduced by the so - called subordination through the process @xmath120 , given by the probability density function @xcite @xmath121^{1/\alpha } \frac{t}{s^{1 + 1/\alpha}}\right),\ ] ] where @xmath122 is a one - sided lvy stable probability density with laplace transform @xmath123 .
thus , eq .  ( [ le2 ] ) transforms the brownian process @xmath118 with diffusivity @xmath119 into the subdiffusive motion with generalised diffusivity @xmath28 . on the basis of this scheme , the ensemble averaged position - position correlation in an arbitrary confining potential @xmath124 becomes @xmath125 valid in the limit @xmath126)^{1/\alpha}$ ] , where @xmath127 is the smallest non - zero eigenvalue of the corresponding fokker - planck operator @xcite .
( [ ctrw_conf_ea_auto ] ) demonstrates that , despite the confinement , the process is non - stationary . in result ( [ ctrw_conf_ea_auto ] )
we used the incomplete beta function @xmath128 and defined the first and second boltzmann moments , whose general definition is @xmath129 the partition function reads @xmath130 at large time separation , @xmath131 , the position autocorrelation decays algebraically , @xmath132 towards the value @xmath133 .
the corresponding limiting behaviour of the incomplete beta function reads @xmath134    inserting expression ( [ ctrw_conf_ea_auto ] ) into the definition of the time averaged mean squared displacement ( [ eatamsd ] ) we obtain @xcite @xmath135dt.\end{aligned}\ ] ] then , with relation ( [ bet_lim ] ) we arrive at the scaling behaviour @xmath136 valid in the limits @xmath137 and @xmath138 ) ^{1/\alpha}$ ] .
result ( [ conf_tamsd ] ) is quite remarkable : instead of the naively assumed saturation toward the stationary plateau a power - law growth @xmath139 is observed .
only when the lag time @xmath1 approaches the overall measurement time @xmath22 the singularity in expression ( [ eatamsd ] ) causes a dip toward the plateau of the ensemble averaged mean squared displacement . additionally eqs .
( [ ctrw_conf_ea_auto ] ) and ( [ conf_tamsd ] ) are universal in the sense that the exact form of the confining potential solely enters into the prefactor through the first two moments of the boltzmann distribution corresponding to the confining potential @xmath124 .
we note that the asymptotic scaling @xmath140 is consistent with the numerical analysis presented in ref .
@xcite .
[ confined_ctrw ] depicts the ensemble mean of the time averaged mean squared displacement ( [ eatamsd ] ) for two types of confinement , an harmonic potential and a box potential .
the particle is initially placed at the bottom of the potential well and in the middle of the box potential , respectively .
this is why at short lag times the particle exhibits the linear scaling @xmath141 typical for unconfined motion , before turning over to the confinement - induced scaling @xmath7 .
note that this plot is fit - free , i.e. , the theoretically calculated asymptotic behaviours nicely fall on top of the simulations results .
for an harmonic binding potential @xmath142 ( @xmath143 ) and a particle in a box of length @xmath144 ( @xmath145 ) , with @xmath146 , measurement time @xmath147 ( a.u . ) , @xmath148 , and @xmath149 . without fitting ,
the lines show the analytic results for the transition from @xmath150 for short lag times according to eq .
( [ ctrwtamsd ] ) ( @xmath151 and @xmath152 ) to @xmath7 for long lag times , eq .
( [ conf_tamsd ] ) (  ) . for long @xmath1 ,
@xmath153 exhibits universal behaviour , in the sense that the curve does not depend on the external field.,width=332 ]    as shown in fig .
[ conf_scatter ] individual trajectories still exhibit the pronounced scatter typical for ctrw subdiffusion .
visually the scatter does not change between the unbiased initial motion and the confinement - dominated part of the process .
this is due to the fact that the scatter is caused by the scale - freeness of the waiting times . in our process waiting times and jump lengths are decoupled , corresponding to the subordination property of ctrw subdiffusion @xcite .    .
the extent of the deviations between the trajectories does not change qualitatively between the initial unbiased motion and the later confinement - dominated regime .
the squares ( @xmath143 ) represent simulations results for the ensemble average @xmath154 .
same parameters as in fig .
[ confined_ctrw].,width=332 ]      what happens if we introduce a cutoff in the power - law of the waiting time distribution of the form @xmath155,\ ] ] in which a characteristic time scale @xmath156 is introduced , terminating the power - law scaling ? for the waiting time density ( [ wtd_cutoff ] ) the characteristic waiting time @xmath157 becomes finite . at sufficiently short times
one would expect this process to still exhibit the features of ctrw subdiffusion , while at times @xmath158 the process should converge to regular brownian motion with a gaussian propagator . in fig .
[ ctrw_cutoff_conf ] we demonstrate that for a suitable choice of the cutoff time @xmath156 the behaviour of the time averaged mean squared displacement @xmath13 under confinement preserves the characteristic non - ergodic features of ctrw subdiffusion , i.e. , the turnover from the initial scaling @xmath3 to the confinement - dominated scaling @xmath159 .
below we will show that at the same time the distribution of the time average around its ensemble mean is significantly altered .    ) , in the box @xmath160 $ ] with reflecting boundary conditions .
the generic behaviour of initial and final scaling , @xmath161 and @xmath7 remains unaltered .
we chose @xmath162 , @xmath163 , @xmath164 and overall measurement time @xmath165 ( a.u.).,width=332 ]      ergodicity remains unaffected for fbm , that is confined to an interval of size @xmath108 $ ] .
namely , convergence of @xmath13 is observed toward the stationary value @xmath109 @xcite . for fbm under the influence of an harmonic potential @xmath110
the position - position correlation can be calculated exactly @xcite .
inserting into the ensemble mean ( [ eatamsd ] ) of the time averaged mean squared displacement one can show that the initial behaviour @xmath166 turns over to the stationary value @xmath167 $ ] . in general , the time averaged mean squared displacement for confined fbm converges to a constant : @xmath168
the deviations of the time averaged mean squared displacement @xmath169 between individual trajectories can be quantified in terms of the probability density function @xmath170 of the dimensionless ratio @xmath171 of the time averaged mean squared displacement over its ensemble mean .
for an ergodic process this distribution is necessarily sharp , @xmath172 for long measurement time @xmath22 .
deviations from this form are expected for non - ergodic processes such as ctrw subdiffusion , but also for relatively short trajectories .
we here address both effects and demonstrate that the distribution @xmath170 is a quite reliable means to distinguish different stochastic processes even when the recorded time series are fairly short .      for ctrw subdiffusion with power - law waiting time distribution ( [ wtd ] )
the distribution assumes the form @xcite @xmath173 for @xmath174 , where @xmath122 is again a one - sided lvy stable distribution , whose laplace transform is @xmath175 .
special cases include the brownian limit ( [ phi1 ] ) for @xmath176 and the gaussian shape @xmath177 for @xmath146 .
in fact , we observe an exponentially fast decay of @xmath170 with large @xmath178 for all @xmath9 , compare appendix [ app1 ] .    for short trajectories
the basic shape of the distribution @xmath170 is surprisingly well preserved @xcite .
this is demonstrated in fig .
[ ctrw_scatter ] , in which the characteristic asymmetric shape with respect to the ergodic value @xmath179 for the case @xmath97 is reproduced for a process with @xmath180 , even for a lag time as long as @xmath181 .
no significant dependence on the confinement is observed .
also for larger values of @xmath42 the finite value for small @xmath178 is similarly reproduced for short trajectories @xcite , pointing at the quite remarkably reliability on the shape of the scatter distribution @xmath170 for ctrw subdiffusion .
as we show below , the distribution for fbm processes with its zero value at @xmath182 can be clearly distinguished from the ctrw form .     of the time
averaged mean squared displacement , with @xmath183 for a ctrw process with @xmath97 , @xmath163 , and @xmath180 ( a.u . ) .
the filled and open symbols , respectively , represent the unconfined case and confined motion on the interval @xmath184 $ ] .
we see quite good agreement with the expected gaussian limiting distribution of the scatter , eq .
( [ phi_gauss ] ) , centred around @xmath182 . in the inset we show the case for the largest measured lag time , @xmath185 , also in good agreement with the predicted shape of the distribution.,width=332 ]    a useful parameter to quantify the violation of ergodicity is the ergodicity breaking parameter @xcite @xmath186 eb varies from @xmath187 for @xmath188 monotonically to @xmath189 in the brownian limit @xmath27 . for the special case @xmath146
one finds @xmath190 , while for @xmath191 , @xmath192 .
the scatter distribution is sensitive to few extreme events .
when these are lacking , the distribution should be significantly different from the form discussed for ctrw subdiffusion with diverging characteristic time scale .
indeed , as shown in fig .
[ cutoff_scatter ] for the cutoff waiting time distribution defined through eq .
( [ wtd_cutoff ] ) , @xmath170 drops down to zero around @xmath182 and assumes an almost gaussian shape around the ergodic value @xmath179 .
note that the simulations were performed with the same parameters as for fig .
[ ctrw_cutoff_conf ] , in which the breaking of ergodicity still persists , despite the cutoff .
this demonstrates that different quantities have different sensitivity to the cutoff .     for ctrw subdiffusion with cutoff , as defined in eq .
( [ wtd_cutoff ] ) .
the observed form is clearly different from the much more asymmetric shape in fig .
[ ctrw_scatter ] .
the value of @xmath170 reaches zero for small values of @xmath178 . an almost gaussian shape centred around the ergodic value @xmath179
is observed .
same parameters as in fig .
[ ctrw_cutoff_conf].,width=332 ]      fbm is based on long - ranged correlations .
however , we can obtain an approximate expression for the scatter distribution @xcite @xmath193 where @xmath194 is an intrinsic time scale .
expression ( [ fbmscatter ] ) is valid for sufficiently small lag times @xmath1 , for which the correlations are neglected @xcite .
note that this distribution is independent of @xmath42 and centred around the ergodic value @xmath179 .
in particular , in the long measurement time limit @xmath174 , the distribution converges to the sharp form @xmath195 .
also this behaviour is surprisingly well preserved for short trajectories , as demonstrated in fig .
[ fbm_scatter ] .     of the time averaged mean squared displacement , with @xmath196 for fbm subdiffusion with @xmath97 and @xmath180 ( a.u . ) .
the filled and open symbols , respectively , represent the cases for unconfined and confined motion . for small values of the lag time
@xmath1 we see quite good agreement with the expected gaussian limiting distribution of the scatter , eq .
( [ fbmscatter ] ) , that is centred around @xmath179 . for larger values of @xmath1 deviations
are observed , however , even for @xmath185 the value around @xmath182 is consistently zero .
a similar behaviour is found for larger values of @xmath42 @xcite .
note that when @xmath174 the process is ergodic and @xmath197 approaches a @xmath198 function.,width=332 ]    for fbm one can define a quantity similar to the ergodicity breaking parameter ( [ eb ] ) .
namely , without taking the long measurement time limit , we obtain the normalised variance of the time averaged mean squared displacement @xmath199 it turns out that the convergence to ergodicity is algebraically slow , and for subdiffusion @xmath9 one obtains a decay , which is inversely proportional to @xmath22 @xcite : @xmath200 the coefficient @xmath201 here is defined as @xcite @xmath202 it increases continuously from zero at @xmath203 to @xmath204 .    the quantity @xmath205 is connected to the probability density ( [ fbmscatter ] ) by @xmath206 , obtained from approximating that the values of @xmath178 may range in the interval @xmath207 , which appears reasonable given the sharp decay of @xmath197 to 0 at @xmath182 for fbm .
apart from the fact that @xmath205 explicitly depends on @xmath42 while the approximation leading to eq .
( [ fbmscatter ] ) loses the @xmath42 dependence , both quantities decay @xmath208 , and the internal time scale @xmath194 takes on a role similar to the lag time @xmath1 .
a typical quantity accessible from experimental data is the velocity autocorrelation function , which is defined through @xmath209 here the velocity is defined as the difference quotient @xmath210 $ ] .
the velocity autocorrelation ( [ autov00 ] ) was suggested as a tool to distinguish between different subdiffusion models @xcite .
we show here that for confined processes the shape of the velocity autocorrelation function does not allow for a significant distinction between subdiffusive ctrw and fbm .
note that the calculation of @xmath211 amounts to determine four two - point correlation functions of the type @xmath212 .      for unbounded ctrw process starting with initial condition @xmath213 , the position correlation function
is given by the following expression @xcite , @xmath214^\alpha .
\label{autov03}\ ] ] this result for free ctrw is due to the fact that the jump lengths in the interval @xmath215 for @xmath216 have zero mean . with eq .
( [ autov03 ] ) we find @xmath217 0 & \tau\geq\epsilon \end{array}\right . \label{autov04}\ ] ] the velocity autocorrelation function for an unbounded ctrw process does not yield negative values due to the absence of correlations for different jumps . in this case
the velocity autocorrelations can be easily distinguished from the behaviour of fbm processes , see below .
note that @xmath211 in eq .
( [ autov04 ] ) is non - analytic at @xmath218 , an observation that is confirmed in simulations .
the behaviour of the velocity autocorrelations for ctrw subdiffusion are displayed in fig .
[ velcorr1 ] .    for a confined ctrw process
the situation is quite different . to explore the behaviour of the corresponding velocity autocorrelation function we use the general result for the correlation function of confined ctrw , eq .
( [ ctrw_conf_ea_auto ] ) .
for simplicity we assume a symmetric potential such that @xmath219 . with the initial condition
@xmath213 we obtain @xmath220 figs .
[ velcor2 ] and [ velcor3 ] ( for the absolute value ) display excellent agreement of eq .
( [ autov05 ] ) with simulations of ctrw subdiffusion for a lattice of size 10 with reflecting boundary conditions , and @xmath146 .
we observe that for confined motion the ctrw velocity autocorrelation function indeed attains negative values and has a minimum on @xmath218 , as the confinement effectively induces correlations .
for long @xmath52 the velocity autocorrelation function decays to zero ( from the negative side ) as the power - law @xmath221 valid for @xmath222 .      for free fbm
the position correlation function ( [ fbm_auto ] ) together with the definition ( [ autov00 ] ) produce the result @xmath223 this function yields negative values for sufficiently long @xmath52 , and its minimum value @xmath224 is assumed at @xmath218 . for long @xmath52 it decays toward zero from the negative side in the power - law form
@xmath225 valid for @xmath222 .     for confined fbm with @xmath97 .
the simulations were performed on the interval @xmath184 $ ] , and @xmath226 ( a.u . ) . for comparison ,
the velocity autocorrelation for free fbm is drawn based on eq .
( [ autov01]).,width=332 ]    we see that both the velocity autocorrelation function of the confined ctrw and the one for free fbm acquire negative values and decay as power - laws for large @xmath52 , see below .
further both obtain a sharp minimum for @xmath218 .
confined fbm behaves similarly to free fbm , see fig .
[ fbm_box ] . from our analysis
it becomes clear that the shape of the velocity autocorrelation function may not be a good diagnosis tool to distinguish between subdiffusive ctrw and fbm processes .
we note that for ctrw subdiffusion the time averaged correlation functions will be random variables , exactly like the time averaged mean squared displacement .
the physical mechanisms leading to subdiffusion in biological cells or other complex systems are varied . the anomalous diffusion of inert biopolymers larger than some 10kd in living cells is due to molecular crowding , these are excluded volume effects in the viscoelastic , superdense cellular environment @xcite .
potentially similar effects occur in single file diffusion : a tracer particle diffusing in a one dimensional channel and interacting with other brownian particles will slow down its random motion , and subdiffusion with @xmath146 is observed @xcite . what is the resulting motion like ? consider a particle , that gets repeatedly trapped during its random walk .
such trapping may give rise to a broad distribution of waiting times , as embodied in the ctrw model .
such broad distribution of trapping times could be due to chemical binding of the tracer particle to its environment on varied time scales , or due to active gelation / de - gelation of the crowding particles .
long - tailed trapping time distributions were observed for @xmath227m - sized particles due to interaction with a semiflexible actin mesh @xcite .
there , the anomalous diffusion exponent @xmath42 depends on the size of the tracer particles versus the typical actin mesh size .
alternatively , motion patterns of fbm type may be due to coupling to the viscoelastic crowding environment @xcite . yet other approaches are based on polymer dynamics .
de gennes reptation model of a polymer diffusing in a tube formed by foreign chains in a melt yields subdiffusion with @xmath228 , and was used as a possible explanation of the short time dynamics of telomere motion @xcite .    one may speculate whether the widely observed subdiffusion of biomolecules and passive tracers is a mere coincidence of nature that should be attributed to the dense environment in the cell and/or to the wide distribution of obstacles and reactions in the cell . or
, maybe , subdiffusion is by itself a goal which was obtained via evolution ?
there exist claims that sub - diffusion is useful in certain cellular search strategies @xcite .
one may also view subdiffusion to be a sufficiently slow process to help maintaining the organisation of the genome in the nucleus of the cell without the need for physical compartments @xcite . since reactions in many cases are controlled by diffusion , the emergence of slower - than - normal diffusion has far - reaching implications to signalling and regulatory processes in the cell .
the standard theories of diffusion - controlled reactions under such circumstances must be replaced by subdiffusion - controlled reaction models . while these issues are clearly important , we here focused more on the characterisation of trajectories of single particles in the cell , in particular , as a quantitative diagnosis method to probe the nature of the stochastic motion in cells and other complex media .
single particle tracking microscopy is a widely used technique .
it allows one to locally probe complex systems in the liquid phase _ in situ_. in particular it has become one of the standard tools in biophysical , colloidal , polymeric , and gel - like environments .
the motion of tracers in these systems is often anomalous .
typically the experimentally recorded time series are evaluated in terms of the time averaged mean squared displacement @xmath13 .
here we collected the behaviour of @xmath229 for the two most prominent anomalous stochastic processes , the continuous time random walk and fractional brownian motion .
for ctrw subdiffusion , connected with ageing and weak ergodicity breaking , in an unconfined system the time averaged mean squared displacement scales linearly in time , @xmath13 [ eq .
( [ ctrwtamsd ] ) ] , renouncing the anomalous nature of the process .
only the dependence on the overall measurement time @xmath22 pays tribute to the underlying subdiffusion .
this behaviour contrasts the anomalous scaling @xmath230 of the ensemble averaged mean squared displacement . under confinement
no plateau value is observed as in the ensemble average , instead , a power - law of the form @xmath43 , eq .
( [ conf_tamsd ] ) , is found .
interestingly , this characteristic behaviour is approximately preserved when an appropriate cutoff in the waiting time is introduced .
conversely the distribution of the time averaged mean squared displacement around its ensemble mean becomes almost gaussian in presence of the cutoff , while it has an exponential decay and a finite value at @xmath182 when the system is non - ergodic and exhibits ageing .
fbm , in contrast , is ergodic , although ergodicity is reached algebraically slowly .
free fbm subdiffusion shows @xmath231 , eq .
( [ fbm_free ] ) , while under confinement the plateau value of the ensemble average is reached , eq .
( [ fbm_cnf ] ) . for finite trajectories
the distribution of the time averaged mean squared displacement is approximately gaussian for short lag times .
our analysis demonstrates how different the two stochastic processes are , despite sharing the same form of the ensemble averaged mean squared displacement . at the same time
the velocity autocorrelation of confined ctrw subdiffusion is hardly distinguishable from that of fbm subdiffusion .
given a recorded time series of anomalous diffusion from experiment it is important to know more precisely which stochastic process is responsible for the observed behaviour , in particular , with respect to diffusion - limited reactions , general transport behaviour , and related processes such as gene regulation . the analysis presented here , along with complementary tools discussed in refs .
@xcite , will be instrumental in the classification of anomalous diffusion behaviour .
we thank c. selhuber - unkel and l. oddershede for providing the data used in fig .
[ christine ] , and for many useful discussions .
partial financial support from the deutsche forschungsgemeinschaft and the israel science foundation is gratefully acknowledged .
using the theory of fox @xmath232-functions we obtain the exact form for the distribution @xmath170 of the dimensionless variable @xmath178 @xcite : @xmath233.\ ] ] this function has the series expansion @xmath234^n}{n!\gamma(-\alpha n)},\ ] ] and the asymptotic form @xmath235^{1/(1-\alpha)}\right)\end{aligned}\ ] ] valid at @xmath236 .
s. b. yuste , g. oshanin , k. lindenberg , o. bnichou , and j. klafter , phys .
e * 78 * , 021105 ( 2008 ) ; e. abad .
s. b. yuste , and k. lindenberg , phys .
e * 81 * , 031115 ( 2010 ) ; i. m. sokolov , s. b. yuste , j. j. ruiz - lorenzo , and k. lindenberg , phys .
e * 79 * , 051113 ( 2009 ) ; d. froemberg and i. m. sokolov , phys .
lett . * 100 * , 108304 ( 2008 ) .
bouchaud , j. phys .
( paris ) i , * 2 * , 1705 ( 1992 ) ; g. bel and e. barkai , phys . rev
. lett . * 94 * , 240602 ( 2005 ) ; a. rebenshtok and e. barkai , _ ibid . _ * 99 * , 210601 ( 2007 ) ; m. a. lomholt , i. m. zaid , and r. metzler , _ ibid . _ , 98 , 200603 ( 2007 ) ; i. m. zaid , m. a. lomholt , and r. metzler , biophys .
j. * 97 * , 710 ( 2009 ) ; f. d. stefani , j. p. hoogenboom , and e. barkai , physics today * 62*(2 ) , 34 ( 2009 ) ; g. margolin and e. barkai , phys .
94 * , 080601 ( 2005 ) ; g. aquino , p. grigolini , and b. j. west , epl * 80 * , 10002 ( 2007 ) .
r. metzler , e. barkai , and j. klafter , phys .
lett . * 82 * , 3563 ( 1999 ) ; europhys . lett . * 46 * , 431 ( 1999 ) .
e. barkai , r. metzler , and j. klafter , phys .
e * 61 * , 132 ( 2000 ) ; r. metzler , j. klafter , and i. m. sokolov , _ ibid . _ * 58 * , 1621 ( 1998 ) .
a. n. kolmogorov , dokl .
26 * , 115 ( 1940 ) ; b. b. mandelbrot and j. w. van ness , siam rev . * 1 * , 422 ( 1968 ) ; h. qian , fractional brownian motion and fractional gaussian noise . in g. rangarajan and m.z . ding ( eds ) , _ processes with long - range correlations _ ( springer , lecture notes in physics , vol.621 ) , pp.22 - 33 .
t. e. harris , j. appl .
* 2 * , 323 ( 1965 ) ; l. lizana and t. ambjrnsson , phys .
. lett . * 100 * , 200601 ( 2008 ) ; l. lizana , t. ambjrnsson , a. taloni , e. barkai , and m. a. lomholt , phys .
e * 81 * , 051118 ( 2010 ) ; e. barkai and r. silbey , phys .
* 102 * , 050602 ( 2009 ) .
m. magdziarz , a. weron , k. burnecki , and j. klafter , phys .
* 103 * , 180602 ( 2009 ) .
compare also m. magdziarz and j. klafter , phys .
e , at press ; k. burnecki and a. weron , phys .
e * 82 * , 021130 ( 2010 ) .
h. c. fogedby , phys .
e * 50 * , 1657 ( 1994 ) ; a. baule and r. friedrich , phys .
e * 71 * , 026101 ( 2005 ) ; m. magdziarz , a. weron , and j. klafter , phys . rev .
* 101 * , 210601 ( 2008 ) ; m. magdziarz , a. weron , and k. weron , phys .
e * 75 * , 016708 ( 2007 ) ; e. heinsalu , m. patriarca , i. goychuk , g. schmid , and p. hnggi , phys . rev .
e * 73 * , 046133 ( 2006 ) .
a. m. mathai and r. k. saxena , the h - function with applications in statistics and other disciplines ( wiley eastern ltd . , new delhi , 1978 )
; a. m. mathai , r. k. saxena , and h. j. haubold , the h - function , theory and applications ( springer , new york , 2010 ) .


in previous papers ( kato 2004 , 2008a ) , we proposed a resonant excitation mechanism of disk oscillations in deformed disks .
the purpose is to propose an excitation mechanism of quasi - periodic oscillations ( qpos ) observed in neutron - star and black - hole x - ray binaries ( e.g. , kato and fukue 2007 , kato 2008b ) . in this resonant excitation model of disk oscillations ,
a deformation of disks from an axially - symmetric equilibrium state is essential .
the deformation to be considered is a warp or an eccentric deformation of disks in the equatorial plane .
an outline of the model is as follows .
a non - linear coupling between a disk oscillation ( hereafter we call it the original oscillation ) and a deformed part of the disks ( warp or eccentric deformation ) brings about some forced disk oscillations ( we call them intermediate oscillations ) .
the intermediate oscillations make a resonant coupling with the unperturbed disk at particular radii of the disk . after this resonant coupling , the intermediate oscillations feedback to the original oscillation by a nonlinear coupling with the deformed part of the disk . since this nonlinear feedback process involves a resonance , the original oscillation is amplified or dampened .
kato ( 2004 , 2008a ) examined this nonlinear feedback processes , and derived an excitation criterion and growth ( damping ) rate of the resonant oscillations .
it is of importance to note that in keplerian disks this resonant excitation process works only when the disks are general relativistic .
that is , a non - monotonic radial distribution of epicyclic frequency is necessary for appearance of resonance and for trapping of oscillations .    in these studies
the intermediate disk oscillations were assumed to be local near to the resonant radius in the sense that their radial wavelength is shorter than the characteristic radial length of disks .
that is , the radial derivative , @xmath0 , operated to the wave quantities was taken to be @xmath1 .
further , @xmath2 was assumed to be constant .
this treatment of @xmath3 const .
, however , was inadequate around the resonant point .
if this inadequate treatment is properly corrected , we find that unlike in the previous papers , the resonant point is not the place where a local dispersion relation of the intermediate oscillations is satisfied with a constant @xmath2 , but the radii of lindblad resonance for intermediate oscillations ( when resonance occurs by horizontal motions ) or the radii where the intermediate oscillations are trapped in the vertical direction ( when resonance occurs by vertical motions ) .    in this paper , we modify the analyses in our previous papers so that the above - mentioned inadequate treatment is corrected .
we apply the mathematical techniques used by meyer - vernet and sicardy ( 1987 ) in their study of resonant disk - satellite interaction .
the results show that the stability criterion derived in the previous papers are unchanged .
the expression for growth rate obtained in the previous papers is found to be still applicable , but in this paper the expression is changed in a form more suitable to numerical calculations .
we first summarize in section 2 the basic equations and relations necessary to the study of the present resonant excitation problem , although they are given in the previous papers . in section 3 ,
the stability criterion and growth rate are derived to the case of pressure - less disks as it is instructive , and to the case of disks with pressure in section 4 .
the final section is devoted to a brief discussion on a meaning of the instability criterion .
we briefly summarize here basic equations and relations [ see kato ( 2008a ) for details ] that are necessary in the subsequent sections .      in the present problem ,
the general relativity is essential . for simplicity ,
however , formulation in this paper is done in the framework of a pseudo - newtonian using the gravitational potential introduced by paczy ' nski and wiita ( 1980 ) , assuming that the central object has no rotation .
we adopt a lagrangian formulation by lynden - bell and ostriker ( 1967 ) .
the unperturbed disk is in a steady state with a steady flow @xmath4 . by using a displacement vector , @xmath5 , a weakly nonlinear hydrodynamical equation describing adiabatic , non - self - gravitating perturbations
is written as , after generalizing the linear equation derived by lynden - bell and ostriker ( 1967 ) , @xmath6 where @xmath7 is a linear hermitian operator with respect to @xmath5 ( lynden - bell and ostriker 1967 ) and is @xmath8                             \nonumber
\\         -p_0\nabla({\rm div}\mbox{\boldmath $ \xi$ } )         -\nabla[(\mbox{\boldmath $ \xi$}\cdot\nabla)p_0 ]         + ( \mbox{\boldmath $ \xi$}\cdot\nabla)(\nabla p_0 ) , \label{2.4}\end{aligned}\ ] ] and @xmath9 and @xmath10 are the density and pressure in the unperturbed state , and @xmath11 is the barotropic index specifying the linear part of the relation between lagrangian variations @xmath12 and @xmath13 , i.e. , @xmath14 . since , the self - gravity of the disk gas has been neglected , the gravitational potential , @xmath15 , is a given function and there is no eulerian perturbation of it . at the stage of equation ( [ 2.3 ] ) , there is no restriction on the form of @xmath4 , but in the main part of this paper , @xmath4 is taken to be a pure cylindrical rotation .    the right - hand side of wave equation ( [ 2.3 ] ) represents the weakly nonlinear terms .
derivation of a detailed expression for @xmath16 is boring , because of the presence of many terms .
however , after lengthy manipulations we can summarize them into a relatively simple form in the case of @xmath17 : @xmath18 using cartesican cordinates .
when @xmath19 , some terms should be added on the right - hand side , which are shown in equation ( 82 ) of kato ( 2008a ) .
an important characteristics of @xmath16 is that we have commutative relations ( kato 2008a ) for an arbitrary set of @xmath20 , @xmath21 , and @xmath22 : @xmath23 as shown later , the presence of these commutative relations is the reason why we have a simple stability criterion . we suppose that this commutative relation is a general property of conservative systems .
we now assume that the disks are deformed from an axisymmetric steady state by some external or internal cause .
the deformation is assumed , for simplicity , to be time - independent , and has the azimuthal wavenumber @xmath24 .
the lagrangian displacement associated with the deformation , @xmath25 , is denoted by , @xmath5 [ see equation ( [ wave ] ) ] , and @xmath26 [ see equation ( [ intermediate ] ) ] are expressed in complex forms .
when their nonlinear couplings are calculated , we must be careful since the real displacement associated with deformations , for example , is not @xmath27 itself , but the real part of @xmath27 given here . ] @xmath28 as the deformation , a warp will be the most probable , but it is not only one candidate of possible deformations .
a plane - symmetric one - armed spiral deformation is one of another possible candidates that can excite disk oscillations . in both cases we have @xmath29 .
our purpose here is to examine how the behavior of disk oscillations are affected by the disk deformation .
as shown below , some of oscillation modes are resonantly excited on deformed disks through nonlinear coupling with disk deformation .
the nonlinear coupling processes are schematically shown in figure 1 of kato ( 2004 ) .    the displacement vector associated with a disk oscillation , @xmath5 ,
is assumed to have frequency @xmath30 and azimuthal wavenumber @xmath31 . separating the time and azimuthal dependences from @xmath32 , we introduce @xmath33 as @xmath34        \hat { \mbox{\boldmath $ \xi$}}(r , z ) .
\label{wave}\ ] ]    the first step of the nonlinear interaction between the disk oscillation characterized by @xmath35 and the deformation characterized by @xmath36 introduces two kinds of intermediate oscillations with azimuthal wavenumber being @xmath37 and @xmath38 .
let us denote the displacement vector associated with these intermediate oscillations by @xmath39           { \hat{\mbox{\boldmath $ \xi$}}}_\pm^{\rm int}(r , z ) , \label{intermediate}\ ] ] where @xmath40 or @xmath38 .
here , @xmath41 represents the intermediate oscillations resulting from the coupling between @xmath42 and @xmath43 , while @xmath44 does those resulting from the coupling between @xmath42 and @xmath45 , where the asterisk represents the complex conjugate .
that is , @xmath41 and @xmath44 are described , respectively , by @xmath46 , \label{2.7}\ ] ] @xmath47 .
\label{2.8}\ ] ]    next , the second stage of the nonlinear coupling is considered , which is a feedback process returning to the original oscillation , @xmath48 , by the intermediate oscillations @xmath49 interacting with @xmath50 . in the case where the coupling occurs through @xmath51 , the feedback is described by @xmath52 .
\label{2.9}\ ] ] on the other hand , in the case where the feedback occurs through @xmath44 , the equation corresponding to equation ( [ 2.9 ] ) is @xmath53 .
\label{2.10}\ ] ]    an important point to be noted here is that as a result of this feedback process , the original disk oscillation is amplified or dampened , i.e. , the frequency , @xmath30 , can be no longer real , since in the feedback processes a resonance is involved .
how much is the imaginary part of @xmath30 ?
this can be examined from equations ( [ 2.9 ] ) and ( [ 2.10 ] ) using the fact that the operators @xmath54 and @xmath55 are hermitian ( lynden - bell and ostriker 1967 ) . after some calculations the imaginary part of frequency , @xmath56 , can be formally expressed as ( e.g. , kato 2008a ) @xmath57 where @xmath58 denotes the cases of the coupling through @xmath59 and @xmath60 , respectively , and @xmath61 dv , \label{2.16}\ ] ] @xmath62 dv , \label{2.18}\ ] ] and @xmath63                          { \hat{\mbox{\boldmath $ \xi$}}}dv , \label{2.12}\ ] ] where @xmath64 is the frequency of the original oscillation before the mode couplings .
the above expressions for @xmath65 and @xmath66 have physical meanings that @xmath65 is the rate at which work is done on the original oscillations by the nonlinear resonant processes , and @xmath66 is the wave energy of the original oscillations [ see equation ( 93 ) by kato ( 2001 ) ] .
it is important to change the above expressions for @xmath65 by using the commutative relation ( [ commutative ] ) as @xmath67dv , \label{work+}\ ] ] @xmath68dv .
\label{work-}\ ] ]      hereafter , we restrict our attention to oscillations in geometrically thin disks .
the steady unperturbed disks are axially - symmetric and have no motion except for rotation . here , cylindrical coordinates ( @xmath69 , @xmath70 , @xmath71 ) are employed , in which the @xmath71-axis is perpendicular to the disk plane and the origin of the coordinates is at the disk center . the unperturbed flow is then described as @xmath72 , where @xmath73 is the angular velocity of disk rotation .
we further assume , for simplicity , that the disk is isothermal in the vertical direction . in the vertically isothermal disks , the unperturbed density in disks , @xmath74 ,
is stratified as ( e.g. , kato et al .
1998 ) @xmath75 , \label{rho}\ ] ] where @xmath76 is the density on the equatorial plane , and @xmath77 is the half - thickness of the disk and is related to the vertical epicyclic frequency , @xmath78 , by @xmath79 where @xmath80 is the isothermal acoustic speed .
the vertical epicyclic frequency , @xmath78 , is equal to the angular velocity of the keplerian rotation , @xmath81 , in the case of the central object being non - rotating , and practically equal to @xmath82 , since the disk is assumed to be geometrically thin .
hereafter , however , we use @xmath78 without using @xmath81 or @xmath82 , so that we can trace back the effects of @xmath78 on the final results .
concerning oscillations superposed on such disks , we consider oscillations whose radial wavelength is moderately short so that the radial variations of physical quantities in the unperturbed state [ including the radial variation of @xmath83 are neglected compared with the radial variation of wave quantities .
furthermore , we assume that the oscillations also occur isothermally .
then , we can neglect the term of @xmath84 $ ] in the operator @xmath55 in equation ( [ 2.4 ] ) .
if the above simplifications and approximations are adopted , the @xmath69- and @xmath71- dependences of @xmath85 are approximately separated .
that is , the @xmath71-dependence of @xmath85 is obtained by solving an eigen - value problem , and the eigen - functions are found to be an orthogonal set of hermite polynomials ( okazaki et al .
1987 ) . restricting our attention to one of such oscillations , we write @xmath85 as @xmath86 @xmath87 @xmath88 where @xmath89 is the hermite polynomial of argument @xmath90 , and
@xmath91 characterizes the number of node(s ) of oscillations in the vertical direction .
it is noted that the number of node(s ) of @xmath92 in the vertical direction is smaller than those of @xmath93 and @xmath94 by one , as shown in equation ( [ xi - z ] ) .
however , the subscript @xmath95 ( not @xmath96 ) is attached to @xmath97 as @xmath98 in order to emphasize that @xmath99 , @xmath100 , and @xmath101 are a set of solutions .    under these approximations ,
we express the @xmath69- , @xmath70- , and @xmath71- components of the homogeneous parts of wave equation ( [ 2.3 ] ) as @xmath102\breve{\xi}_{r , n }           -i2\omega(\omega - m\omega)\breve{\xi}_{\varphi , n }           + \omega^2_\bot
h{d \breve{\xi}_{z , n}\over dr }                  = 0 , \label{4.1}\ ] ] @xmath103 @xmath104\breve{\xi}_{z , n }              -n\omega^2_\bot h{d \breve{\xi}_{r , n}\over dr }          = 0 . \label{4.3}\ ] ]    nonlinear coupling between @xmath105 and @xmath106 introduce intermediate oscillations @xmath107 .
the azimuthal wavenumber of the intermediate oscillations are @xmath37 or @xmath38 .
their @xmath71-dependence are characterized by @xmath108 or @xmath96 when @xmath109 since @xmath110 and by @xmath95 when @xmath111 . to consider these various coupling cases separately , we write @xmath112 in the forms of @xmath113 @xmath114 @xmath115 where the subscript @xmath116 represents @xmath108 or @xmath96 or @xmath95 .
nonlinear coupling terms are also separated into terms proportional to exp@xmath117 $ ] and @xmath118 .
that is , in the case of coupling through @xmath119 , we write the coupling terms as @xmath120_r         = \rho_0\sum_{{\tilde m},\tilde{n } } \breve{a}_{r,+,\tilde{n}}(r )                 { \rm exp}[i(\omega t-{\tilde m}\varphi ) ]                  { \cal h}_{\tilde { n}}(z / h)+\cdot\cdot\cdot \label{5.7}\ ] ] @xmath120_\varphi         = \rho_0\sum_{{\tilde m},\tilde{n } } \breve{a}_{\varphi,+,\tilde{n}}(r )                  { \rm exp}[i(\omega t-{\tilde m}\varphi ) ]                  { \cal h}_{\tilde { n}}(z / h)+\cdot\cdot\cdot \label{5.7'}\ ] ] @xmath120_z         = \rho_0\sum_{{\tilde m},\tilde{n } } \breve{a}_{z,+,\tilde{n}}(r )                  { \rm exp}[i(\omega t-{\tilde m}\varphi ) ]                  { \cal h}_{\tilde { n}-1}(z / h)+\cdot\cdot\cdot , \label{5.7''}\ ] ] where @xmath121 denotes terms orthogonal both to @xmath122 and @xmath123 , and the subscript @xmath124 is added to @xmath125 s in order to emphasize that they are related to the @xmath70-dependence of @xmath126 $ ] . in a similar way , in the case of coupling through the intermediate oscillations of @xmath127 , @xmath128 $ ] can be expressed in forms similar to equations ( [ 5.7 ] )  ( [ 5.7  ] ) , introducing @xmath129 , @xmath130 , @xmath131 .
the equations describing intermediate oscillations are then written as @xmath132               \breve{\xi}^{\rm int}_{r,\pm,{\tilde n } }           -i2\omega(\omega-{\tilde m}\omega )               \breve{\xi}^{\rm int}_{\varphi,\pm,{\tilde
n } }   \nonumber
\\           + \omega^2_\bot h{\partial \breve{\xi}^{\rm int}_{z,\pm,{\tilde n } }                 \over\partial r }                  = \breve{a}_{r,\pm,\tilde{n } } , \label{4.1'}\end{aligned}\ ] ] @xmath133 @xmath134              \breve{\xi}^{\rm int}_{z,\pm , { \tilde n } }              -{\tilde n}\omega^2_\bot h{d               \breve{\xi}^{\rm int}_{r,\pm,{\tilde n}}\over dr }          = \breve{a}_{z , \pm,\tilde{n } } ,   \label{4.3'}\ ] ] where @xmath135 represents @xmath37 or @xmath38 , and @xmath136 does @xmath108 or @xmath96 when @xmath109 and @xmath137 when @xmath111 . from equations ( [ 4.1 ] ) and ( [ 4.2 ] ) we can eliminate @xmath138 to give @xmath139\breve{\xi}^{\rm int}_{r,\pm , { \tilde n } }         -c_{\rm s}^2{d^2 \breve{\xi}^{\rm int}_{r,\pm , { \tilde n } }                \over dr^2 }         + \omega^2_\bot h { d\breve{\xi}^{\rm int}_{z,\pm , { \tilde n } }              \over dr }          = \breve{a}_{r,\pm,{\tilde n } }           -i{2\omega\over\omega-{\tilde m}\omega}\breve{a}_{\varphi,\pm,{\tilde n}}. \label{4.1''}\ ] ] hereafter , we use equations ( [ 4.1 ] )  ( [ 4.1  ] ) as the basic equations describing the intermediate oscillations .    finally , we should notice that in the case of isothermal disks described above , @xmath140 and @xmath141 given by equations ( [ work+ ] ) and ( [ work- ] ) are written in the forms of @xmath142          \cdot \biggr[{\breve \xi}_{r,+}^{\rm int}{\breve a}_{r,+}^ *                + { \breve \xi}_{\varphi,+}^{\rm int}{\breve a}_{\varphi,+}^ *                + \frac{1}{{\tilde n}}{\breve \xi}_{z , + } { \breve a}_{z,+}^*\biggr]dx , \label{w+}\ ] ] and @xmath143          \cdot \biggr[{\breve \xi}_{r,-}^{\rm int}{\breve a}_{r,-}^ *                + { \breve \xi}_{\varphi,-}^{\rm int}{\breve a}_{\varphi,-}^
*                + \frac{1}{{\tilde n}}{\breve \xi}_{z , -}{\breve a}_{z,-}^*\biggr]dx , \label{w-}\ ] ] where @xmath144 $ ] comes from the part of volume integration in the @xmath70- and @xmath71-directions .
furthermore , the wave energy , @xmath66 , given by equation ( [ 2.12 ] ) can be expressed as @xmath145 where @xmath146 is a dimensionless quantity given by is different from that given in the previous paper ( kato 2008a ) , since the previous one is not suitable to represent the sign of wave energy . ] @xmath147 and the subscript c denotes the values at resonant radius , which is defined by @xmath148 or @xmath149 ( see below ) .
it is noted that in the case of @xmath150 , @xmath151 is zero .
before examining general cases of disks with pressure , we study here the limiting case of pressure - less disks , i.e. , @xmath152 and @xmath153 . equations ( [ 4.1  ] ) and ( [ 4.3 ] ) show that @xmath154 ( and @xmath155 ) and @xmath156 respond resonantly , respectively , at the radii where @xmath157 holds and at the radii where @xmath158 holds .
we define here @xmath159 and @xmath160 by @xmath161 we call the resonance at @xmath148 horizontal resonances ( lindblad resonances ) and the resonance at @xmath149 vertical resonances .
we consider these two resonances separately .    as a typical example in infinitesimally thin disks , we consider here the case where all motions including those associated with disk deformations are in the equatorial plane , i.e. , @xmath150 , @xmath111 , and @xmath162 .
, @xmath163 , and @xmath136 in general without specifying to particular values . ] at the radii of @xmath148 , @xmath154 and @xmath155 resonantly respond to @xmath125 s as [ see equations ( [ 4.1  ] ) and ( [ 4.2 ] ) ] @xmath164 @xmath165 .
\label{3.2}\ ] ] the quantity @xmath156 can be taken to be zero , since it does not respond resonantly and has no contribution to the final results .    to examine the resonant interaction at the radii of @xmath148
, we introduce a small imaginary part of @xmath30 , i.e. , @xmath56 , as @xmath166 , where the imaginary part , @xmath56 , is tentatively assumed to be negative so that causality is satisfied .
later , the results obtained by @xmath167 is extended analytically to the whole region of @xmath56 as usually done in stability analyses .
hereafter , the real part , @xmath168 , is written by @xmath30 or @xmath64 , without confusion .
near to the resonance , @xmath169 is expanded as @xmath170 , \label{3.4}\ ] ] where the subscript c represents the values at the resonant radii , and @xmath171 is @xmath172_{\rm c}. \label{3.5}\ ] ] let us now consider an integration of @xmath173 , @xmath174 being an arbitrary smooth real function , along the radial direction including the resonant region .
the imaginary part of the integration comes from the path near to the pole .
thus , we have @xmath175 this can be summarized as @xmath176 where sign@xmath177 is the sign of @xmath178 at the resonant radii .
let us now estimate the rate of work done by the resonance on the original oscillations by using equation ( [ 3.6 ] ) . in the case of coupling through @xmath179 ,
the work is given by equations ( [ w+ ] ) and ( [ w- ] ) .
since @xmath180 and @xmath181 are given by equations ( [ 3.1 ] ) and ( [ 3.2 ] ) , respectively , we have from equation ( [ w+ ] ) and ( [ w- ] ) @xmath182 where we have used @xmath183 , and subscript @xmath77 is attached to @xmath56 in order to emphasize that this is a case of horizontal resonances .
hence , the growth rate , @xmath184 [ see equation ( [ growth ] ) ] , can be expressed as @xmath185 it is noted that the above expression for growth rate is valid even when @xmath186
. equations ( [ 3.7 ] )  ( [ 5.30 ] ) formally hold even in the case with pressure , as shown in the next section
.    the vertical resonance might be unrealistic in pressure - less disks , but we briefly summarize the formal results here , since the expression for growth rate is applicable even in the case of disks with pressure , as shown in the next section .
we now expand @xmath187 around @xmath149 as @xmath188 , \label{5.32}\ ] ] where @xmath189_{\rm c } , \label{5.33}\ ] ] and @xmath190 is the resonant radius of @xmath149 . in
the above two cases of the horizontal and vertical resonances , the resonant radii are different each other .
however , the same notation , @xmath190 , is used here and hereafter without confusion .
after this preparation , we calculate the rate of work done on oscillations by the same procedures as those in the case of horizontal resonance .
the results show that the growth rate is given by @xmath191
equations describing the intermediate oscillations [ equations ( [ 4.1 ] )  ( [ 4.3 ] ) or equations ( [ 4.1  ] ) and ( [ 4.3 ] ) ] are now solved without assuming @xmath152 and @xmath153 .
since we are interested in behaviors of @xmath192 in a region close to the resonant point of @xmath148 or @xmath149 , all the coefficients of @xmath192 in equations ( [ 4.1 ] )  ( [ 4.3 ] ) are assumed to be constant except for @xmath169 or @xmath187 .
that is , in the case of the horizontal resonance ( lindblad resonance ) , @xmath193 , and @xmath169 is taken as [ see equations ( [ 3.4 ] ) and ( [ 3.5 ] ) ] @xmath194 where @xmath195 is the radial distance from the resonant radius , @xmath190 , i.e. , @xmath196 . in the case of the vertical resonance , on the other hand , @xmath197 , and @xmath187
is taken as [ see equations ( [ 5.32 ] ) and ( [ 5.33 ] ) ] @xmath198 where @xmath196 .
the resonant radii in the horizontal and vertical resonances are different each other , but the same notation , @xmath190 , has been adopted here without confusion .    in disks with pressure , different from the case of pressure - less disks ,
the resonant region is not narrow .
that is , the resonant region is widened by the effects of pressure ( see below ) , and thus the coupling terms , @xmath199 , @xmath200 , and @xmath201 , can be no longer regarded to be spatially constant in general .
considering these situations , we introduce fourier transforms of @xmath202 and @xmath203 as @xmath204 @xmath205 after these preparations , to solve the set of equations ( [ 4.1 ] )  ( [ 4.3 ] ) , we apply a method of fourier transform adopted by meyer - vernet and sicardy ( 1987 ) in their study of resonant disk - satellite interaction . here and hereafter , such subscripts as @xmath58 and @xmath136 to be attached to @xmath192 , @xmath206 , @xmath207 , and @xmath208 are sometimes omitted in order to avoid unnecessary complications , unless they are explicitly necessary .      in the present case ,
the fourier transform of equations ( [ 4.1  ] ) and ( [ 4.3 ] ) give , respectively , @xmath209 and @xmath210 where @xmath211 .
eliminating @xmath212 from equation ( [ 4.4 ] ) and ( [ 4.5 ] ) , we have an inhomogeneous differential equation of @xmath213 with respect to @xmath2 : @xmath214 where @xmath215 is defined by @xmath216 let us first consider the case of @xmath217 .
equation ( [ 4.6 ] ) is then solved to lead @xmath218 , \label{4.8}\ ] ] where @xmath219dk ' . \label{4.8a}\ ] ] in determining the integration range of equation ( [ 4.8a ] ) , we have used a boundary condition allowing @xmath220 for @xmath221 .
this boundary condition is required , since @xmath213 is the fourier transform of @xmath222 .    by using the inverse fourier
transform : @xmath223 we can express @xmath222 in an integration form as @xmath224dk , \label{xi - r - again}\ ] ] where @xmath225 is defined by @xmath226 since @xmath227 and @xmath228 are related to @xmath229 by [ see equations ( [ 4.2 ] ) and ( [ 4.3 ] ) ] @xmath230 the work integral , @xmath231 , in the case of the coupling through @xmath232 , is written as , using equation ( [ w+ ] ) , @xmath233         d_{(+)}(k ) { \rm exp}\biggr[-\frac{\epsilon}{j_1'}k+i\biggr(kx-\frac{1}{3}\alpha_{\rm p}^3k^3\biggr )         \biggr ] \nonumber \\         & & \times\biggr[{\breve a}_{r,+}^*(x)+i\frac{2\omega}{\omega-{\tilde m }               \omega}{\breve a}_{\varphi,+}^*(x )           + \frac{1}{j_2}ik\omega_\bot^2h{\breve a}_{z,+}^*(x)\biggr ] , \label{w+2}\end{aligned}\ ] ] where real parts on the right - hand side of equation ( [ w+2 ] ) have been omitted .    in the case of @xmath234 , the expression for @xmath235 is slightly changed from equation ( [ 4.8 ] ) as @xmath236 , \label{4.8'}\ ] ] where @xmath237 is given by @xmath238dk ' , \label{4.8a'}\ ] ] where @xmath239 is defined by equation ( [ 4.7 ] ) . the integration range in equation ( [ 4.8a ] )
is changed from equation ( [ 4.8a ] ) so that the boundary condition @xmath240 is satisfied at @xmath241 . by performing the same procedures as those in the case of @xmath217
, we can derive an expression for @xmath231 .
an expression for @xmath231 in this case is the same as equation ( [ w+2 ] ) , except that @xmath242 is now changed to @xmath243 , i.e. , @xmath233         d_{(-)}(k ) { \rm exp}\biggr[-\frac{\epsilon}{j_1'}k+i\biggr(kx-\frac{1}{3}\alpha_{\rm p}^3k^3\biggr )         \biggr ] \nonumber \\         & & \times\biggr[{\breve a}_{r,+}^*(x)+i\frac{2\omega}{\omega-{\tilde m }               \omega}{\breve a}_{\varphi,+}^*(x )           + \frac{1}{j_2}ik\omega_\bot^2h{\breve a}_{z,+}^*(x)\biggr ] , \label{w+2a}\end{aligned}\ ] ]    expressions of the work integral in the case of the coupling through @xmath244 can be derived by similar procedures as those in the case of the coupling through @xmath245 . in this case , the subscript + attached to @xmath246 s in equations ( [ w+2 ] ) and ( [ w+2a ] ) is changed to @xmath247 and @xmath248 .
first , we consider the case where the coupling terms , @xmath199 , @xmath200 , and @xmath201 , can be regarded to be spatially constant in the resonant region .
( this is really realized in a particular case , as mentioned in the next section . )
an another purpose of examining the case of constant coupling terms is to have a rough image on the width of resonant region in disks with pressure .
if @xmath249 , @xmath250 , and @xmath251 are spatially constant , their fourier transforms give @xmath252 in this case , @xmath253 given by equation ( [ 4.8a ] ) is reduced simply to @xmath254 , \label{d}\ ] ] where @xmath255 is the unit - step function , i.e. , @xmath256 for @xmath257 , while @xmath258 for @xmath259 . since @xmath253 does not involve @xmath195 , the integration with respect to @xmath195 in equation ( [ w+2 ] ) leads to @xmath260 .
hence , performing the integration with respect to @xmath2 in equation ( [ w+2 ] ) , we have @xmath261 in the case of @xmath217 , where @xmath40 and @xmath262 , depending on the type of couplings and the form of deformation .    in the case of @xmath234 , @xmath213
is given by equation ( [ 4.8 ] ) . by performing the same procedures as those in the case of @xmath217
, we can derive an expression for @xmath231 . the final expression for @xmath231 in this case
has the opposite sign from equation ( [ w+3 ] ) .
inequality @xmath217 implies that @xmath263 , since @xmath184 should be taken tentatively to be negative as mentioned before , while inequality @xmath264 means @xmath263 . considering this difference
, we can summarize @xmath231 in the above two cases of @xmath217 and @xmath234 as @xmath265 where @xmath40 .    in the case of coupling through @xmath266
, we have , after the similar processes as the above , @xmath267 where @xmath248 .    finally , from @xmath268 given above and the wave energy , @xmath66 , given by equation ( [ 2.12 ] ) , we find that the growth rate , equation ( [ 5.28 ] ) , is written as @xmath269 which is identical to equation ( [ 5.30 ] ) , and also equal to the growth rate derived in the previous paper ( kato 2008a ) , as will be mentioned at the end of this section .
the arguments in subsection 4.1.1 show that in the case where @xmath249 , @xmath250 , and @xmath251 are radially constant , the radial structure of resonant region is determined by [ see equation ( [ w+2 ] ) and notice that @xmath253 has the step - function ] @xmath270,\ ] ] where the term with @xmath271 can be practically neglected .
the @xmath195-dependence of this integral has been examined by meyer - vernet and sicardy ( 1987 ) [ see equation ( 41 ) and fig , 5 in their paper ] .
they show that the half width of the resonant region is @xmath272 , which is on the order of @xmath273 , i.e. , @xmath274 , in the present case .
in general , @xmath249 , @xmath250 , and @xmath251 radially oscillate with the wavelength determined by the product of @xmath275 and @xmath276 .
hence , in some cases the scale - length of the variation of @xmath277 s will be shorter than @xmath274 , although it depends on where the resonance occurs ( see the next section ) . to examine general cases ,
numerical calculations are needed .
however , in order to have a rough image of the effects of spatial variation of @xmath249 , @xmath250 , and @xmath251 on the work integral @xmath268 , we consider here a limiting case where @xmath249 , @xmath250 , and @xmath251 oscillate periodically in the radial direction .
that is , we assume @xmath278 where @xmath279 s and @xmath280 are the amplitudes and wavenumber of the spatial variations .
then , the fourier transform defined by equation ( [ 4.2b ] ) show that @xmath281    let us consider the case of @xmath217 . in this case ,
@xmath253 given by equation ( [ 4.8a ] ) is reduced to [ cf . , equation ( [ d ] ) ] @xmath282               { \rm exp}\biggr[\frac{\epsilon}{j_1'}k_0+i\frac{\alpha_{\rm p}^3}{3}k_0 ^ 3\biggr ] .
\label{newd}\ ] ]    since @xmath249 , @xmath250 , and @xmath251 are given by equation ( [ periodic ] ) , performance of integration with respect to @xmath195 in equation ( [ w+2 ] ) gives @xmath283 .
next , we perform the integration with respect to @xmath2 .
since @xmath253 is now given by equation ( [ newd ] ) , we have @xmath284    similar arguments can be made in the case of @xmath234 . in this case , @xmath231 has the sign opposite to equation ( [ w+4 ] ) .
expressions for @xmath231 in two cases of @xmath217 and @xmath234 can be combined into a single form as before .
similarly , we can consider the case where the coupling occurs through @xmath285 .
all of these results can be summarized as @xmath286 hence , the growth rate is written as @xmath287      in the case of the vertical resonance , from equations ( [ 4.1  ] ) and ( [ 4.3 ] ) , we have , respectively , @xmath288 and @xmath289 eliminating @xmath290 from the above two equations , we have an inhomogeneous differential equation with respect to @xmath291 as @xmath292 , \label{v4}\ ] ] where @xmath293 in deriving equation ( [ v4 ] ) , in order to avoid unnecessary complications , we have neglected @xmath294 compared with @xmath169 , since the magnitude of @xmath295 contributing to the present formulations is smaller than unity , i.e. , @xmath296 .    in the case of @xmath297 , equation ( [ v4 ] ) is solved as @xmath298 , \label{v6}\ ] ] where the coefficient @xmath299 is given by @xmath300 \nonumber \\         & & \times { \rm exp }         \biggr[\frac{\epsilon}{j_2'}k'+i{1\over 3}\alpha_{\rm p , v}^3k'^3\biggr]dk ' . \label{v7}\end{aligned}\ ] ] by using @xmath212 given above , we can write @xmath290 and @xmath301 as [ see equations ( [ 4.1  ] ) and ( [ 4.2 ] ) ] , @xmath302 and @xmath303 substitution of these relations into the work integral ( [ w+ ] ) gives @xmath304        d_{{\rm v},(+)}(k ) { \rm exp}\biggr[-\frac{\epsilon}{j_2'}k+         i\biggr(kx-\frac{1}{3}\alpha_{\rm p , v}^3k^3\biggr)\biggr ]   \nonumber \\       & & \times         \biggr[-i\frac{ik\omega_\bot^2 h}{j_1 }               \biggr({\bar a}_r^*+i\frac{2\omega}{\omega-{\tilde m}\omega}{\bar a}_\varphi^*\biggr )                  + \frac{1}{{\tilde n}}{\bar a}_z^*\biggr ] .        \label{v10}\end{aligned}\ ] ]    in the case of @xmath305 , @xmath306 has a similar form , but we should use , instead of @xmath299 , @xmath307 , which is given by @xmath308   \nonumber \\            & & \times { \rm exp }         \biggr[\frac{\epsilon}{j_2'}k'+i{1\over 3}\alpha_{\rm p , v}^3k'^3\biggr]dk ' . \label{v7'}\end{aligned}\ ] ]    as in the subsection 4.1.1 , we consider two cases where i ) @xmath251 is constant in the resonant region and ii ) @xmath251 is spatially periodic as @xmath309 . in the former case , summing the two cases of @xmath310 and @xmath311 , we have @xmath312 and in the latter case we have @xmath313    in the coupling through @xmath314 , the work integral , @xmath315 , has the same expression as equations ( [ v11a ] ) and ( [ v11 ] ) , except that the subscript @xmath124 is now changed to @xmath247 .
after the above considerations , we summarize the growth rate , @xmath316 , in the former case of @xmath201 being constant as @xmath317 in the case where @xmath251 changes as @xmath318 , we have @xmath319    finally , we note here that the growth rates obtained in the previous papers are the same as those obtained in subsection 4.1.1 , although the final expression in the previous papers are different from those in subsection 4.1.1 . in the case of horizontal resonance , for example ,
@xmath320 in kato ( 2008a ) is equal to @xmath321 in the present paper .
furthermore , if @xmath322 and @xmath323 are eliminated from equations ( 46 )  ( 48 ) in kato ( 2008a ) , we have @xmath324 where @xmath325 is @xmath326 .
this expression for @xmath229 shows that @xmath327 defined by equation ( 52 ) in kato ( 2008a ) is written as @xmath328 then , substituting @xmath329 and equation ( [ * * ] ) with @xmath330 into equation ( 70 ) in kato ( 2008a ) , we obtain an expression for growth rate in the case of horizontal resonances , which is identical with equation ( [ 4.19 ] ) in subsection 4.1.1 .
similar arguments easily show that the growth rate given by kato ( 2008a ) for vertical resonances , i.e. , equation ( 71 ) in kato ( 2008a ) , is identical with that in section 4.2 , i.e. , equation ( [ v11 ] ) .
in this paper we have derived a stability criterion on resonant excitation of disk oscillations in a deformed disk .
deformation of disks that can lead to excitation of disk oscillations are a warp and an eccentric deformation of disk plane , i.e. , @xmath331 , although this is not discussed in this paper ( see kato 2004 , 2008a ) . in previous papers ( kato 2004 , 2008a ) , the above resonant excitation problem was examined and a condition of excitation of disk oscillations was derived . from mathematical point of view , however , the treatment of disk oscillations around the resonant point was not suitable in the previous papers . in this paper , we corrected the inadequate points .
the results of the correction show that the stability condition and the growth rate derived in the previous papers are kept qualitatively unchanged . in this paper , however , the growth rate of oscillations is written in a form applicable more easily to numerical calculations of growth rate .    in disks with pressure ,
the resonant region is widened by pressure effects . based on arguments by meyer - vernet and sicardy ( 1987 )
, we showed in section 4 that the half - width of the resonant region is @xmath332 ( or @xmath333 ) , which is shorter than @xmath69 but longer than @xmath77 .
the coupling terms comes from nonlinear products of original oscillation , @xmath334 , and disk deformation , @xmath335 .
if the disk deformation is assumed to be global , the scale of spatial variation of the coupling terms is mainly governed by the spatial behavior of the original oscillation , @xmath334 .
hence , in general , the width of resonant region is comparable with the wavelength of the original oscillations and may not be narrow enough to regard the coupling terms @xmath125 s as constant in the region .
hence , in section 4 , we examined two limiting cases where i ) the coupling terms are constant in the resonant region and ii ) they vary sinusoidally in the region . analyses in section 4
show that in these two limiting cases there are no essential differences in growth rate , especially no difference in the criterion of stability . to examine what happens in intermediate cases of the above two limiting ones ,
however , numerical calculations will be needed .
it is noted that the case where the coupling terms , @xmath203 , are roughly constant in the resonant region is realized when resonance occurs near to the boundary between the propagation and evanescent regions of the the original oscillations .
this is because near to the boundary , the radial wavelength of the oscillations is long .
applications of the present disk oscillation model to high - frequency qpos by kato and fukue ( 2006 ) and kato ( 2008a ) are concerned in such cases .
the another limiting case where the coupling terms vary in the radial direction will be realized when resonance occurs in the propagation region far apart from the boundary between the propagation and evanescent regions .
there are two kinds of resonances , i.e. , horizontal and vertical resonances . in both resonances ,
the condition of excitation of oscillations is that the wave energy of the original oscillation , @xmath146 , and @xmath178 at the resonant radius have the opposite signs , i.e. , @xmath336 . for the resonance to occur efficiently
, the resonant radius must be within the region where both original and intermediate oscillations exists predominantly .
hence , it will be reasonable to suppose that the sign of @xmath177 roughly represents the sign of the energy of the intermediate oscillation .
[ notice that the sign of wave energy of the intermediate oscillation is determined by the sign of @xmath178 in the region where the wave exists predominantly , see equations ( [ 5.28 ] ) and ( [ 5.29 ] ) ] . in other words , the condition of @xmath337 is roughly equal to the condition that the original and intermediate oscillations have energies of the opposite signs .
this instability criterion thus can be understand as a result of energy exchange between the original and intermediate oscillations at the resonant point ( ferreira and ogilvie 2008 ) .
the resonant instability condition , i.e. , @xmath338 , however , is not exactly equal to the condition of both oscillaions ( original and intermediate one ) having the opposite signs .
this suggests that the energy exchange at the resonant point is not only between the above two oscillations , but the disk rotation also plays a part and gains or losses its energy . in the case of @xmath339 , the original oscillation is damped . roughly speaking ,
in this case , both oscillations ( original and intermediate ones ) have the same signs in their wave energy .
hence , in this case , the main energy flow will be between disk rotation and the oscillations so that the oscillations are damped .
it is noticed that in vertical resonance we always have @xmath339 in the case where a resonant radius can appear in disk plane and the oscillations are always damped ( see kato 2004 , 2008a ) .
finally , it is noted that the present resonant excitation mechanism has been applied to quasi - periodic oscillations observed in black - hole and neutron - star x - ray binaries , for example , by kato and fukue ( 2006 ) , kato ( 2008a ) and kato et al .