instabilities in plastic flow has been an object of attention for a long time in metallurgical literature .
conceptually simplest form of this [ lubahn & felgar , 1961 ; hall , 1970 ] manifests when the material is subjected to a creep test wherein a force is applied and the response in the form of elongation of the material is measured , which under normal conditions , is a smooth strain - time curve
. however , under certain metallurgical conditions , one sees steps on creep curve suggesting a form of instability [ lubahn & felgar , 1961 ; hall , 1970 ; da silveira & monteiro , 1979 ] .
a better known form of the instability having the same physical origin , but conceptually more difficult , appears when the material is deformed under tensile deformation [ bodner & rosen , 1967 ; brindley & worthington , 1970 ; penning , 1972 ] . here
, the material is subjected to a predetermined response namely a constant rate of deformation , and the force or the stress developed in the sample is sought to be measured . even in this case , one finds a smooth stress - strain curve under normal conditions .
however , when the system is in the regime of instability , for some values of material parameters the stress - strain curve exhibits multiple load drops .
each of the load drops is related to the formation and propagation of dislocation bands [ chihab _ et al _ , 1987 ] .
it is in the latter type of testing where plastic instability manifests much more easily than in the former .
the phenomenon is referred to as the portevin - le chatelier ( plc ) effect or the jerky flow and is seen in several metals such as commercial aluminium , brass , alloys of aluminium and magnesium [ brindley & worthington , 1970 ] .
it is observed only in a window of strain rates and temperature .
in contrast , the phenomenon of steps on creep curve , which is the subject of the present discussion , is seen in few instances [ da silveira & monteiro , 1979 ; zagarukuyko _ et al _ , 1977 ; stejskalova _ et al _ , 1981 ] .
the reason attributed to this is that it is difficult to obtain proper control on metallurgical parameters wherein this form of instability can be observed .
there is one more form of instability which manifests under constant stress test where similar strain jumps are seen as a function of time .
the origin of instability in all these three modes of testing is known to be the same .
( for a current status of both experiments and theory see kubin _
et al _ , 1993 ) .
eventhough steps on creep curve are seen in a limited number of experiments , it is often straight forward to translate many of the experimental results of constant strain rate test into that of creep test .
it is generally agreed that the microscopic origin of the instabilities arises due to the interaction of dislocations with mobile point defects and is referred to as dynamic strain ageing .
this leads to negative strain rate characteristic .
the basic idea was formulated by cottrell [ 1953 ] few decades ago .
the early phenomenological models are all static since they do not deal with time development .
in contrast , methodology of dynamical systems addresses precisely this aspect . until about a decade and a half ago
, there were no models which looked at the problem from the point of view of bifurcation theory .
an attempt to understand the problem in the above perspective was first made by our group several years ago [ ananthakrishna & sahoo , 1981b ; ananthakrishna & valsakumar , 1982 , 1983 ; valsakumar & ananthakrishna , 1983 ] . in a series of papers starting from an extended fokker - planck equation for the distribution function of the velocity of dislocation segments , and then splitting this into a mobile and an immobile component , we arrived at a model which consisted of three types of dislocations and some transformations between them [ ananthakrishna & sahoo , 1981a , 1981b ; sahoo & ananthakrishna , 1982 ; valsakumar & ananthakrishna , 1983 ]
. the basic idea could be summarized by stating that the phenomena is due to a hopf bifurcation resulting from nonlinear interactions between three different types of dislocations , suggesting a new mathematical mechanism for the instability .
eventhough the spatial inhomogeneous structure was ignored and only the temporal oscillatory state was sought to be described , the model and its extensions to the case of constant strain rate test , proved to be very successful in that it could explain most of the experimentally observed features ( for critical reviews see kubin & lepinoux , 1988 and kubin _ et al _ , 1993 ) .
_ for instance , the existence of negative strain rate sensitivity in a regime of strain rates , comes out naturally as a consequence of the hopf bifurcation _ [ ananthakrishna & valsakumar , 1982 ] .
it must be emphasized , that this feature has been measured in most experiments on the plc effect , and is _ assumed _ in all other theoretical models [ cottrell , 1953 ; bodner & rosen , 1967 ; penning , 1972 ; kubin and estrin 1990 ; see also kubin & martin , 1988 ; martin & kubin , 1992 ; kubin _ et al _ 1993 ] . other experimentally observed features such as the existence of bounds on strain rate for the plc effect to occur , the existence of critical strain and its dependence on applied strain rate , the dependence of the amplitude on the strain etc .
, also follow .
_ one other important prediction which is direct consequence of the dynamical basis of the model is the existence of chaotic stress drops in a window of strain rates _ [ ananthakrishna & valsakumar , 1983 ; ananthakrishna & john , 1990 ] .
_ recently this prediction has been verified which in turn implies that only a few degrees of freedom are required for a dynamical description of the phenomenon .
this offers justification for the use of only a few degrees of freedom for the description of the temporal aspect , eventhough the system is spatially extended . _
[ ananthakrishna _ et al_. , 1995 ; ananthakrishna & noronha , 1995 ; quaouire & fressengeas , 1995 ; venkadesan _ et al _ 1995 ; noronha _ et al _ 1996 & 1997 ] .
( note that a spatially extended system implies infinite degrees of freedom . )
description of the phenomenon which includes the initiation and propagation of the bands during the plc effect has also been recently attempted [ ananthakrishna , 1993 ] .    since the introduction of bifurcation theory into this field several years ago by our group
, there has been a resurgence of interest in plastic instabilities in the light of introduction of new methodology borrowed from the theory of dynamical systems .
this has further helped to obtain new insights hitherto not possible [ kubin & martin , 1988 ; estrin & kubin , 1989 ; kubin & estrin , 1990 ; see also kubin & martin , 1988 ; martin & kubin , 1992 ; kubin , 1993 ; ananthakrishna _ et al_. , 1995a & b and references therein ] .
one of the aims of such theories is to be able to relate the microscopic dislocation mechanisms to the measurable macroscopic quantities . however , in the process , we feel that finer aspects of dynamical systems have been glossed over in this field .
for instance , one often finds that casual remarks are made about fast and slow modes without actually going through the procedure of demonstrating the existence of such modes and eliminating the fast modes in favour of the slow ones [ aifantis , 1988 ; hahner & kubin , 1992 ; hahner , 1993 ] .
in addition , under the adiabatic elimination , the slow modes are complicated functions of the original modes . yet ,
hand waiving arguments have been used in building models which we believe are technically suspect .
the purpose of this article is at least three fold .
first , we wish to illustrate the power and utility of dynamical methods to the study of dynamical aspects of the model .
we demonstrate the technical aspects of adiabatic elimination ( largely addressed to metallurgists working in the area of plastic instabilities ) and then derive the equation for the slow modes having the form of time - dependent ginzburg - landau ( tdgl ) equation .
second , we wish to compare the results of the dependence of the amplitude and the period of the strain jumps on stress and temperature with experiments on creep curve as also with the results translated from the constant strain rate case to the creep case .
this will allow us to relate the theoretically introduced parameters to the macroscopically measured quantities .
third , we wish to compare these results and the detailed numerical results via power series expansion with the approximate solutions for the limit cycles obtained earlier [ ananthakrishna & sahoo , 1981b ; valsakumar & ananthakrishna , 1983 ] .
the principal results of this work are as follows .
we show that the derived quintic tdgl equation for the slow modes is valid in the mid - range of the instability domain . within this range
, there is a narrow sub - range of the parameter where the system exhibits a supercritical bifurcation , beyond which subcritical bifurcation is exhibited .
the results of this calculations on the amplitude and period of the strain jumps are compared with experiments .
in addition , analytical work based on the tdgl equation and a detailed comparative numerical study via the power series solution shows that the model exhibits unusual properties in a certain domain of the parameter .
the study further shows that very high order nonlinearities actually govern the nature of the subcritical bifurcation beyond the domain handled by the quintic tdgl equation .
the plan of the paper is as follows .
in what follows ( section 2 ) we present a brief summary of the model . in section 3 , the technique of extracting the tdgl equation using reductive perturbative approach is applied to the model . in section 4 ,
approximate limit cycle solution obtained through the tdgl equation is compared with experiments and with the numerical solution of the model as well as earlier studies on the model .
section 5 contains summary and discussion .
we start with a brief summary of the model .
the basic idea of the model is that instabilities in plastic flow are a consequence of nonlinear interactions between different types of dislocation populations .
spatial dependence is ignored and only temporal aspects are described , the idea being that the few degrees of freedom used for describing the phenomenon ( dislocation populations ) correspond to the collective modes of dislocations in the spatially extended system . as mentioned in the introduction , this finds support from the recent experimental verification of chaos in the plc effect .
the details of the model can be found in the original references [ ananthakrishna & sahoo , 1981b ; valsakumar & ananthakrishna , 1983 ] .
the model consists of mobile dislocations @xmath0 and immobile dislocations @xmath1 and another type which mimics the cottrell s type @xmath2 , which are dislocations with clouds of solute atoms .
let the corresponding densities be @xmath3 , @xmath4 and @xmath5 , respectively .
the basic dislocation mechanisms included are the following : ( a ) production of dislocations by cross glide ( @xmath6 , @xmath7 is the cross glide coefficient , @xmath8 is the velocity of mobile dislocations ; @xmath9 [ alexander , 1986 ] where @xmath10 is the applied stress ) , ( b ) immobilization of two mobile dislocations ( @xmath11 , @xmath12 is the rate at which two mobile dislocations get immobilized ) and annihilation of two mobile dislocations ( @xmath13 ) , and annihilation of a mobile with an immobile dislocation ( @xmath14 ) , ( c ) once a mobile dislocation is immobilized , it can be reactivated by athermal or thermal means ( @xmath15 ) , ( d ) lastly , a mobile dislocation can acquire a cloud of solute atoms and move with them .
we consider such dislocations as distinct from both the mobile as well as the immobile and denote them @xmath2 .
this process is represented by @xmath16 , @xmath17 being a function of the concentration of solute atoms . however , as the solute atoms gather progressively around dislocations , these dislocations eventually will be immobilized ( @xmath18 ) .
( note that in principle , @xmath12 , @xmath19 and @xmath20 are different . ) these mechanisms lead to the rate equations for the densities of dislocations : @xmath21    ( for the sake of simplicity , we take @xmath22 and @xmath23 . ) one can give a more transparent interpretation for eq .
( 3 ) which represents whole process of acquiring a cloud of solute atoms by a mobile dislocation and eventual immobilization .
this can be seen by defining @xmath24 , where the kernel @xmath25 which represents the fact that solute atoms are arriving at a certain rate has been chosen to be @xmath26 .
it can be easily checked that eq .
( 3 ) is the differential form of the definition of @xmath5 .
one last remark about some of the terms appearing in these equations .
the first two terms plus the last term in eq.(1 ) , and equivalent terms in eqs.(2 ) , were derived starting from a velocity distribution function of the dislocation segments [ ananthakrishna & sahoo 1981a ; sahoo & ananthakrishna , 1982 ] .
thus , this model is an improved model to explain the instability in plastic flow .
equations ( 1 - 3 ) can be cast into a dimensionless form by using scaled variables : @xmath27 to get @xmath28 where @xmath29 .
the dot represents differentiation with respect to @xmath30 .
equations ( 5 - 7 ) are coupled set of nonlinear equations which support limit cycle solutions for a range of parameters @xmath31 and @xmath32 , that are physically relevant .
@xmath33 refers to the concentration of the solute atoms , @xmath34 refers to the reactivation of immobile dislocations and @xmath35 to the time scales over which the slowing down occurs .
we have demonstrated the existence of limit cycle solutions and also obtained approximate closed form solutions for the limit cycles [ ananthakrishna & sahoo , 1981b ; valsakumar & ananthakrishna , 1983 ] .
in addition , the model has been studied numerically . using the orowan equation which relates the rate of change of strain(@xmath36 ) to mobile dislocation density and the mean velocity : @xmath37 , with @xmath38 as the burger s vector , steps on the creep curve follow automatically since the densities of dislocations are oscillatory .
several experimental results are reproduced [ ananthakrishna & sahoo , 1981b ; valsakumar & ananthakrishna , 1983 ] .
for the sake of completeness , we outline the reductive perturbative approach . near the point of hopf bifurcation of the system ( eqs.5 - 7 ) , corresponding to a critical value of the drive parameter , a pair of complex conjugate eigenvalues and another real negative eigenvalue exist for the linearized system of equations around the steady state . as we approach the critical value from the stable side , the real parts of the pair of complex conjugate eigenvalues approach zero and
hence the corresponding eigenvectors get slower and slower .
in contrast , the effect of the change in the drive parameter on the real negative eigenvalue is negligible . thus , _ the eigenvector corresponding to the real negative eigenvalue is the fast mode_. for this reason , the slow modes are the ones which determine the formation of new states of order . _
the reductive perturbative method is a method where the asymptotic equation is extracted in a systematic way _ [ taniuti & wei , 1968 ; newell & whitehead , 1969 ; kuramoto & tsuzuki , 1974 ; mashiyama _ et al_. , 1975 ; richter _ et al_. , 1981 ] .
the method involves in first finding the eigenvectors corresponding to the slow modes and looking for a solution in the subspace spanned by these vectors .
the effect of the nonlinearity is handled progressively as a perturbation of the linear solution in power series of the deviation from the critical value of the drive parameter .
it may be worth emphasizing that this method is essentially the same as reduction to center manifold .
indeed , the equivalence of the center manifold theory [ carr , 1981 ; guckenheimer & holmes , 1983 ; troger & steindl , 1991 ] with the reductive perturbation has been established [ chen _ et al _ , 1996 ] .
other techniques of extracting amplitude equations have been devised whose end results are basically the same .
for instance , perturbative renormalization group method [ goldenfeld _ et al _ 1989 ; chen _ et al _ , 1996 ] and its recent extension on the basis of envelope theory [ kunihiro , 1995 & 1996 ] has also been developed as a tool for global asymptotic analysis which can be used to extract the amplitude equations .    in the present model ,
there is one simplification arising from the fact that @xmath39 can be identified as the fast ( stable ) mode , and hence this will be eliminated adiabatically .
then we derive the equation for the slow dynamics through a power series expansion .
although this adiabatic elimination is approximate , the result is found to be in a reasonably good agreement with the accurate method ( which will be reported elsewhere ) except for finer details .
consider the eqs.(5 - 7 ) .
it can be shown that there is a domain of instability for the parameters @xmath33 , @xmath34 and @xmath35 .
the range of values of @xmath34 is @xmath40 , that of @xmath35 is @xmath41 while @xmath42 .
thus , @xmath43 and @xmath44 @xmath45 .
thus , there are three time scales that are in principle different from each other .
physically , they correspond to the mobile ( fastest @xmath33 ) , the cottrell type ( @xmath35 ) and the immobile ( slowest @xmath34 ) . under these conditions
we can eliminate one full equation adiabatically .
this aspect becomes transparent if the above equations are written in terms of variables which are deviations from the steady state .
there is only one steady state defined by : @xmath46^{1/2}}\over{4b_0 } } , { \rm and},y_{a}=1/2.\ ] ]    defining new variables which are deviations from the steady state @xmath47 eqns .
( 5 - 7 ) take the form @xmath48 where @xmath49 now , we can rescale the time - like variable by @xmath50 and get @xmath51 it is clear that in the limit of small @xmath34 , @xmath52 , unless the left hand side vanishes identically .
thus , we can eliminate @xmath53 in favour of the other two variables and obtain : @xmath54 , \\
\mbox{and}\hspace{.5 cm }   { { d{z}}\over{d\tau } } & = & c\left[x -z\right],\end{aligned}\ ] ] where @xmath55^{\frac{1}{2 } } \right).\ ] ] ( the other root for @xmath53 is unphysical since it corresponds to negative dislocation density . )
it must be emphasized that this adiabatic elimination becomes more exact as the value of the parameter @xmath34 gets smaller .
equations ( 17 - 18 ) will be solved reductive perturbatively .
writing these equations as a matrix equation where the nonlinear part appears separately from the linear part , we get @xmath56 where @xmath57 @xmath58 with @xmath59 and the nonlinear part , @xmath60 , is given by @xmath61 the coefficients @xmath62 and @xmath63 appearing in the last expression for @xmath60 are functions of @xmath64 , @xmath65 and @xmath35 whose expressions for the required first few are given in eqs.(a8-a15 ) of appendix a.    consider stability of the fixed point as a function of the parameter @xmath35 . the eigenvalues , @xmath66 , of the matrix @xmath67 are @xmath68 the fixed point becomes unstable when @xmath35 is less than @xmath69 ( @xmath35 is non - negative ) and the discriminant of eq.(24 ) is negative giving a pair of complex conjugate eigenvalues .
this holds when the inequality @xmath70 is satisfied .
the instability region in the @xmath71 versus @xmath33 plane for a fixed value of @xmath72 is shown in fig .
it extends approximately between @xmath73 and @xmath74 .    to get approximate analytical solution of eq.(20 ) , we follow the reductive perturbative approach similar to that used by mashiyama _
et al_. [ 1975 ] and richter _ et al_. [ 1981 ] .
we choose @xmath75 with @xmath76 and write the matrix @xmath67 as a sum of two matrices , @xmath77 , where @xmath78 is the matrix @xmath67 evaluated for @xmath79 and @xmath80 the eigenvalues of @xmath78 are @xmath81 , where @xmath82 . considering that the solution for @xmath83 grows continuously out of the critical eigenmodes
, we can express it as @xmath84 where @xmath85 and @xmath86 are time - dependent complex amplitudes while @xmath87 are right eigenvectors defined by @xmath88 .
similarly , we introduce left eigenvectors , @xmath89 defined by @xmath90 .
( note that @xmath91 .
t refers to transpose operation of the matrix concerned ) .
substituting this expression for @xmath83 in eq.(@xmath92 ) and multiplying it by @xmath93 from the left - hand side gives the equation governing @xmath85 : @xmath94 @xmath95 and @xmath96 are complex coefficients which are functions of @xmath33 , @xmath34 and @xmath35 .
their expressions are , respectively , given in eqs .
( a1 ) and ( a4 ) in appendix a. we express @xmath85 in power series of @xmath97 : @xmath98 and introduce slowly varying multiple time scales @xmath99 , @xmath100 , ... replacing @xmath101 by @xmath102 . with this substitution in eq.(@xmath103 ) , it can be solved successively by equating terms of equal power in @xmath104 on both sides of the equation . at first , @xmath105(@xmath97 ) terms give rise to the equation @xmath106 which means that , on the time scale of @xmath30 , @xmath107 is constant . the next higher order , @xmath108 ,
terms give the equation @xmath109 which can be integrated to get @xmath110 .
@xmath105(@xmath111 ) terms give the equation @xmath112 + p_3 ( \psi_1 e^{i\omega\tau } + c.c . ) ^3 \right ] e^{-i\omega\tau},\ ] ] from which we extract the slow dynamics : @xmath113 since to @xmath114 @xmath115 , eq.(@xmath116 ) gives the dynamics for the complex amplitude @xmath85 : @xmath117 the expression for the complex coefficient @xmath118 is given in eq .
( a2 ) in appendix a. equation ( @xmath119 ) is a _
cubic _ tdgl equation where @xmath85 is the complex amplitude for the slow mode which is also referred to as complex order parameter ( of the new state of temporal order ) .
its steady state solution gives the amplitude squared as @xmath120 and the associated frequency , @xmath121 , ( with @xmath122 ) as @xmath123 note that @xmath85 describes the limit cycle oscillation .
this solution exists provided @xmath124 is negative since @xmath125 is positive ( see expression of @xmath126 in eq.(a1 ) in appendix a ) .
@xmath127 is found to be negative in a narrow region ( @xmath128 ) as can be seen in fig .
2 . in this case , since the amplitude of the slow modes grows continuously in proportion to @xmath97 ( see eq.(@xmath129 ) ) the transition is supercritical bifurcation ( continuous or ` second order ' ) for other values of @xmath33 ( within the instability range ) , @xmath124 is positive implying that the transition is subcritical bifurcation ( discontinuous or ` first order ' ) .
one has to then go to quintic or even higher terms in the tdgl equation for obtaining an expression for the slow modes ( order parameter ) .
we have carried out the reductive perturbative method further to derive the _
quintic _ tdgl equation : @xmath130 the steady state solution of this equation gives the amplitude squared as @xmath131 and the frequency as @xmath132 we found @xmath133 to be negative in the range @xmath134 ( see fig .
_ this means a narrow region : @xmath135 within the subcritical bifurcation exists over which the quintic tdgl equation has a solution _ ( see fig .
to start with , consider the supercritical regime where expressions for the amplitude and period of the limit cycle are simple .
using the steady state solution of the cubic tdgl equation , eqs.(@xmath129 ) & ( @xmath136 ) , the dependence of @xmath137 and period , @xmath138 , of the limit cycle on @xmath33 are , to the leading order , found to be @xmath139 and @xmath140 , respectively .
the parameter @xmath141 is a function of the applied stress ( @xmath10 ) and temperature ( t ) . as remarked earlier , @xmath17 is proportional to the concentration of the solute atoms , and therefore @xmath142 . using the standard expression [ alexander 1986 ] @xmath143 ( with @xmath144 ) , we get @xmath145 from this we see that the amplitude of the limit cycle has a strong increasing dependence on stress and an increasing dependence on temperature if we assume that @xmath146 . in contrast , the period of the limit cycle has a weak decreasing dependence on stress ( as the leading contribution is constant ) and a decreasing dependence on temperature .
since stress and temperature are measurable quantities , our predictions can be compared with experimental results .
the amplitude and period of the limit cycle are related , respectively , to the amount of strain jumps and the period of the jumps on the creep curve through the orowan equation .
there are very few experiments in this mode of testing as mentioned in the introduction .
the only experiment where this dependence on stress and temperature has been measured is that by zagoruyko _
et al _ [ 1977 ] .
however , it is possible to translate the results from experiments in constant strain rate case to the creep case and compare them with the results of this calculation . according to zagoruyko _
et al _ [ 1977 ] , the amplitude of the strain jumps increases with stress while its period has a weak but decreasing dependence on stress .
experiments from constant strain rate case also exhibit the same trend when the results are translated in terms of constant stress experiments .
it is well known that the amplitude of the stress drops decreases with applied strain rate .
in fact , even the present model predicts this behaviour for the constant strain rate mode [ ananthakrishna & valsakumar 1982 ] .
this implies that the dependence of the amplitude of strain jumps on stress should increase [ kubin _ et al _ 1993 ] .
[ this relation can be seen as follows . in constant strain rate case ,
the deformation rate is fixed and the stress developed in the sample is measured . when the contribution to the plastic strain rate increases due to increased dislocation motion ( for whatever reasons ) , the stress has to fall in order to keep the applied strain rate constant .
thus , the relation between strain rate and stress is opposite . ] clearly , the general trend is consistent with the experimental results .
3 shows the actual dependence of the amplitude and the period on @xmath33 in the domain of supercritical bifurcation ( not the leading one as given by the above expressions ) .
zagarukuyko @xmath147 [ 1977 ] also report that the amplitude of the strain jumps increases while its period decreases with temperature , which is consistent with our result ( provided @xmath148 ) .    in the case of subcritical bifurcation , the dependence of the amplitude of the strain jumps and the period on @xmath33 have similar behaviour as in the case of the supercritical bifurcation .
eventhough , in this case , leading order dependence on @xmath33 is difficult to obtain we have their actual dependence plotted in fig .
it may be noted that they have qualitatively the same behaviour as in fig .
3 with the degree of dependence more pronounced in this case .    experiments in constant strain rate case show that the stress drops are seen to arise both abruptly as well as continuously [ kubin _ et al _ , 1988 ] . translating this result to the constant stress case , it implies that the strain jumps can arise both abruptly and continuously .
this feature again follows from our calculations .
having derived the tdgl equation , we would first like to compare its result with the numerical solutions obtained via eqs .
( 17 - 18 ) and via eq.(@xmath92 ) . secondly , due to the fact that the quintic tdgl equation ( eq .
@xmath149 ) is valid in a limited domain , one suspects that higher order nonlinearities are controlling the subcritical bifurcation , which is quite unusual .
therefore , we would like to investigate the numerical solution obtained by solving eq.(@xmath92 ) keeping successive leading powers in @xmath150 .
this analysis confirms the above suspicion that higher order nonlinearities are important in this model outside the domain of validity of the quintic tdgl equation .    using the steady state solution of the quintic tdgl equation , eq.(@xmath149 ) , in eq.(@xmath151 )
, we get an approximate expression for the limit cycle which is usually called the secular motion . the equations governing the secular motion are given by eqs .
( b5 ) and ( b6 ) ( hereafter called secular equations ) in appendix b. this can be compared with the numerical solutions obtained using eqs .
( 17 - 18 ) ( hereafter called the reduced model ) and with that obtained using eq.(@xmath92 ) _ where _ the expansion in @xmath150 of eq . ( 23 ) is limited up to a certain power ( hereafter called the @xmath152-truncated model , n signifying the highest power of @xmath150 after truncation ) .
_ at the outset we state that the reduced model has bounded solution over the entire instability domain_.    in the case of the numerical solutions obtained from eq.(@xmath92 ) , we have found that while truncating the @xmath150-expansion in eq .
( 23 ) at @xmath153 is sufficient to give a bounded solution for the range of supercritical bifurcation ( @xmath154 ) , it fails beyond this . on the other hand , truncating the @xmath150-expansion at @xmath155 does extend the domain of bounded solution ( @xmath156 ) which matches with the domain where the quintic tdgl equation works .
figure 5 shows the plots of the solutions to the secular equations ( eqs .
b5 and b6 ) , the reduced model ( eqs .
17 - 18 ) and the @xmath153-truncated model ( eq .
@xmath92 ) for @xmath157 and @xmath158 .
all of them are numerically almost indistinguishable .
we have verified that the solutions obtained by these three methods show that the amplitude of the limit cycle scales with @xmath97 .
figure 6 shows the solutions obtained by the three ways for @xmath159 and @xmath158 where , this time , the truncation is at @xmath155 .
the agreement between the three solutions is reasonably good .
part of the discrepancy between the secular motion and the two numerical solutions can be attributed to the range of values of the various terms ( @xmath160 s , @xmath161 s , @xmath0 s and @xmath162 s of appendix c ) contributing to @xmath163 and @xmath164 .
( they range from @xmath165 to @xmath166 ) .
we have also verified that the amplitude of the limit cycle has a finite jump for this case .
for @xmath167 , where the quintic tdgl equation is inadequate , @xmath155-truncated model is also inadequate .
however , @xmath168-truncated model has a bounded solution for the region @xmath169 .
in particular , fig .
7 shows the plots of the solutions obtained by the @xmath168-truncated model and that by the reduced model for @xmath170 and @xmath158 .
_ this implies that seventh order tdgl equation need to be considered for this range of @xmath33_. we have verified that @xmath171 term extends the domain of bounded solutions upto @xmath172 . from this trend
, it appears that very high powers in @xmath150 need to be retained to cover the entire domain upto @xmath173 .
the region @xmath174 is even more interesting .
while the quintic tdgl equation is inadequate to give a bounded solution , we found that retaining even upto the ninth power in @xmath150 does not give a bounded solution for the entire range of values @xmath175 ( even for those close to @xmath176 ) .
this suggests that this region is one where the nonlinearity is strong . on the other hand ,
as pointed out earlier the reduced model has limit cycle solution in this range as well .
in particular , fig .
8 shows such a plot for @xmath177 and @xmath158
.    thus , the results of the tdgl equation , the @xmath152-truncated model and the reduced model are consistent with each other .
further , the numerical solutions via the power series in @xmath150 throws light on the nature of nonlinearities governing the solution in different regions of instability .
we have carried out the reductive perturbative approach to the problem of steps on creep curve and shown that the dynamics of the system is described by a tdgl equation for the amplitude of the slow modes ( complex order parameter ) @xmath85 in the neighbourhood of the bifurcation point .
since the above derivation is valid only in the neighbourhood of the critical value , the expression for @xmath85 is valid only for small @xmath104 .
this has been exemplified by the reasonable agreement between the secular motion and the numerical solution of the reduced model within the domain of the validity of the quintic tdgl equation .
we have shown that both subcritical as well as supercritical bifurcations are seen in the instability range of the parameter @xmath33 .
_ while the subcritical bifurcation ( abrupt or ` first order ' transition ) seen over most of the values of @xmath33 is consistent with our earlier calculation _ [ ananthakrishna & sahoo , 1981b ; valsakumar & ananthakrishna , 1983 ] , _ we see a supercritical bifurcation ( continuous or ` second order ' transition ) over a narrow mid - range of _ @xmath33 . the existence of the supercritical bifurcation was missed by us earlier . _ the existence of both supercritical and subcritical bifurcations is consistent with experiments . in addition , the dependence of the amplitude and the periodicity of the strain jumps on stress and temperature are consistent with experiments_. from this point of view , even though the derived tdgl equation is valid in a limited domain , the present calculation allows a direct mapping of these quantities .
this calculation also helps us to demonstrate the complicated dependence of the slow modes on the original modes .
this will serve as a warning to those using hand waiving arguments for declaring certain modes as fast modes and others as slow modes in modelling of such of these problems .    looked at from the point of view of the properties of the model , there are some interesting features . for the region outside the validity of the quintic tdgl equation , on the basis of a careful numerical analysis , we found that high order nonlinearities govern the nature of the subcritical bifurcation .
it must be emphasized that this feature is unusual and this insight would not have been possible but for the numerical solutions obtained by keeping successive higher powers in @xmath150 , which in itself is the basis for the derivation of the tdgl equation .
this unusual feature of the model interpreted in the language of phase transitions give better insight and could also find applicability . for conventional models ,
the ` free energy ' is described by an expansion in power series of the order parameter . even in the case of dynamic transition such as the present one , it is possible to associate a ` free energy ' like function such that @xmath178 @xmath179 , where @xmath180 $ ] is the ` free energy ' like function .
while the ` free energy ' for ` second order ' phase transitions is described by retaining up to fourth power in the order parameter , up to sixth power is conventionally required for ` first order ' phase transitions ( with the appropriate signs for the coefficients in the expansion ) . in the present case , however , we need to go to as high as twelfth power ( or more ) of the order parameter to cover the interval @xmath181 ( or equivalently twelfth power or higher in @xmath150 ) . in this case , including successively higher powers in @xmath150 increases the domain of description . on the other hand , in the interval of @xmath175 , even retaining upto tenth power in @xmath150 ( which is the highest power we have checked ) does not work .
this must be contrasted with the existence of bounded solution for the reduced model over the entire interval . _
this feature is rather unusual and , to the best of our knowledge , we are not aware of any other model which exhibits this property_.      one of us ( m.b . ) would like to thank the international program in physical sciences , uppsala university ( sweden ) for offering a fellowship to study at indian institute of science .
part of this work has been supported by ifcpar grant no . 1108 - 1 .
* references *    aifantis , e. c. [ 1988 ] `` on the problem of dislocation patterning and persistent slip bands '' in _ non linear phenomena in materials science i _ , eds .
kubin , l. p. & martin , g. ( trans tech , switzerland ) , pp .
397 - 406 .
alexander , h. [ 1986 ] `` dislocations in covalent crystals '' in _ dislocations in solids _ ed .
nabarro , f. r. n. ( north - holland ) , 151 .
ananthakrishna , g. [ 1993 ] `` formation , propagation of bands and chaos in jerky flow '' , _ scripta .
_ * 29 * , 1183 - 1188 .
ananthakrishna , g. & john , t. m. [ 1990 ] `` instabilities and chaos in plastic flow '' , in _ directions in chaos _ vol .
3 , ed . hao bai - lin ( world scientific , singapore ) , pp .
133 - 148 .
ananthakrishna , g. , fressengeas c. , grosbras , m. , vergnol , j. , engelke , c. , plessing , j. , neuhaeuser , h. , bouchaud , e. , planes , j. & kubin , l. p. [ 1995a ] `` on the existence of chaos in jerky flow , '' _ scripta metall . _ * 32 * , 1731 - 1737 .
ananthakrishna , g. , kubin , l. p. & martin , g. ( eds . ) [ 1995b ] _ solid state phenomena _ vol .
42&43 ( scitec , switzerland ) .
ananthakrishna , g. & noronha , s. j. [ 1995 ] `` chaos in jerky flow : theory and experiment '' in _ solid state phenomena _ vol .
42&43 , eds .
ananthakrishna , g. , kubin , l. p. & martin , g. ( scitec , switzerland ) , pp .
277 - 286 .
ananthakrishna g. & sahoo , d. [ 1981a ] `` a statistical dislocation dynamics with application to lif '' , _
d _ * 14 * ) , 699 - 713 .
ananthakrishna , g. & sahoo , d. [ 1981b ] `` a model based on nonlinear oscillations to explain jumps on creep curves '' , _
d _ * 14 * , 2081 - 2090 .
ananthakrishna , g. & valsakumar , m. c. [ 1982 ] `` repeated yield drop phenomena : a temporal dissipative structure '' , _ j. phys .
d _ * 15 * , l171-l175 .
ananthakrishna , g. & valsakumar , m. c. [ 1983 ] `` chaotic flow in a model for repeated yielding '' , _ phys .
a _ * 95 * , 69 - 71 .
ananthakrishna , g. , fressengeas c. , grosbras , m. , vergnol , j. , engelke , c. , plessing , j. , neuhaeuser , h. , bouchaud , e. , planes , j. & kubin , l. p. [ 1995
] `` on the existence of chaos in jerky flow , '' _ scripta metall .
_ * 32 * , 1731 - 1737 .
bodner , s. r. & rosen , a. [ 1967 ] `` discontinuous yielding of commercially - pure aluminium '' , _ j. mech . phys . solids _ * 15 * , 63 - 77 .
brindley , b. j. & worthington , p. j. [ 1970 ] `` yield - point phenomena in substitutional alloys '' , _ metall .
reviews _ * 145 * , 101 - 114 .
carr , j. [ 1981 ] _ applications of center manifold theory _ ( springer - verlag , berlin ) .
chen , l. y. , goldenfeld , n. & oono , y. , [ 1996 ] `` the renormalization group and singular perturbations : multiple- scales , boundary layers and reductive perturbation theory '' , _ phys .
e _ * 54 * , 376 - 394 .
chihab , k. , estrin , y. , kubin , l. p. & vergnol , j. [ 1987 ] `` the kinetics of the portevin - le chatelier effect in an al-5 at% mg alloy '' , _ scripta .
_ * 21 * , 203 - 208 .
cottrell , a. h. [ 1953 ] `` a note on the portevin - le chatelier effect '' , _ phil . mag .
_ * 44 * , 829 - 832 .    da silveira , t. l. & s. n. monteiro , s. n. [ 1979 ] `` jumps in the creep curve of austentic stainless steels '' , _ met .
a _ * 10 * , 1795 - 1796 .
estrin , y. & kubin , l. p. [ 1989
] `` collective dislocation behavior in dilute alloys and the portevin - le chatelier effect '' , _ j. mech .
behaviour of materials _ * 2 * , 255 - 292 .
goldenfeld , n. , martin , o. & oono , y. [ 1989 ] `` intermediate asymptotics and renormalization group theory '' , _
comp . _ * 4 * , 355 .
guckenheimer , j. & holmes , p. [ 1983
] _ nonlinear oscillations , dynamical systems , and bifurcations of vector fields _ vol .
42 of _ applied math .
( springer - verlag , berlin ) , 123 - 145 .
hahner , p. & kubin , l. p. [ 1992
] `` coherent propagative structures in plastic deformation : a theory of ludlers bands in polycrystals '' in _ solid state phenomena _ vol .
3&4 , eds .
kubin , l. p. & martin , g. ( trans tech , switzerland ) , pp .
385 - 402 .
hahner , p. [ 1993
] `` modelling the spatiotemporal aspects of the portevin - le chatelier effect '' , _ mater .
sci . and eng .
_ * 164 * , 23 - 34
.    hall , e. o. [ 1970 ] _ yield point phenomena in metals and alloys _ ( macmillan , london ) .
kubin , l. p. , chihab , k. & estrin , y. [ 1988 ] `` the rate dependence of the portevin - le chatelier effect '' , _ acta metall . _ * 36 * , 2707 - 2718 .
kubin , l. p. & estrin , y. [ 1990 ] `` evolution of dislocation densities and the critical conditions for the portevin - le chatelier effect '' , _ acta metall . _ * 38 * , 697 - 708 .
kubin , l. p. & estrin , y. ( organizers ) [ 1993 ] `` viewpoint set no . 21 : propagative instabilities '' , scripta metall .
* 29 * , 1147 - 1188 .
kubin , l. p. & lepinoux , j. [ 1988 ] `` the dynamic organization of dislocation structures '' , _ strength of metals and alloys _
( icsma 8) , vol .
kettunen , p. o. , lepisto , t. k. , & lehtonen , m. e. , ( pergamon press , oxford ) , pp .
35 - 59 .
kubin , l. p. & martin , g. ( eds . ) [ 1988 ] _ solid state phenomena _ vol .
3&4 ( trans tech , switzerland ) .
kunihiro , t. [ 1995 ] `` a geometrical formulation of the renormalization group method for global analysis '' , _ prog .
phys . _ * 94 * 503 - 514 .
kunihiro , t. [ 1996 ] `` the renormalization group method applied to asymptotic analysis of vector fields '' , in hep - th/9609045 .
kuramoto , y. & tsuzuki , t. [ 1974 ] `` reductive perturbative approach to chemical instabilities '' , _ prog .
_ * 52 * , 1399 - 1401 .
lubahn , j. d. & felgar , r. p. [ 1961
] _ plasticity and creep of metals _ ( john wiley , new york ) .
martin , g. & kubin , l. p. ( eds . ) [ 1992 ] _ solid state phenomena _ vol.23&24 ( trans tech , switzerland ) .
mashiyama , h. , ito , a. & ohta , t. `` fluctuations and phase transitions far from equilibrium '' , _ prog .
* 54 * , 1050 - 1066 .
newell , a. c. & whitehead , j. a. [ 1968 ] `` finite bandwidth , finite amplitude convection '' , _ j. fluid mech . _ * 38 * , 279 - 303 .
noronha , s. j. , ananthakrishna , g. , quaouire , l. & fressengeas , c. [ 1997 ] `` chaos in jerky flow- experimental verification of a theoretical prediction '' , _ pramana _ , in press .
noronha , s. j. , ananthakrishna , g. , quaouire , l. , fressengeas , c. & kubin , l. p. [ 1996
] `` chaos in the portevin - le chatelier effect '' , submitted to int .
j. bif . & chaos .    penning , p. [ 1972
] `` mathematics of the portevin - le chatelier effect '' , _ acta metall . _ * 20 * , 1169 - 1175 .
quaouire , l. & fressengeas , c. [ 1995 ] `` dynamical analysis of the portevin - le chatelier effect '' in _ solid state phenomena _ vol .
3&4 ( trans tech , switzerland ) , eds .
ananthakrishna , g. , kubin , l. p. & martin , g. ( scitec , switzerland ) , pp .
293 - 302 .
richter , p. h. , procaccia , i. & ross , j. [ 1981 ] `` chemical instabilities '' in _ advances in chemical physics _
prigogine , i , & rice , s. ( wiley - interscience , new york ) , pp .
217 - 268 .
sahoo , d. & ananthakrishna , g. [ 1982 ] `` a phenomenological dislocation transformation model for the mobile fraction in simple systems '' , _ j. phys .
d _ * 15 * , 1439 - 1449 .
stejskalova , v. , hammersky , m. , luckac , p. , vostry , p. & spursil , b. [ 1981 ] _ czech .
* b31 * , 195 .
taniuti , t. & wei , c. c. [ 1968 ] `` reductive perturbative method in nonlinear wave propagation i '' , _ j. phys .
japan _ * 24 * , 941 - 946 .
troger , h. & steindl , a. [ 1991 ] _ nonlinear stability and bifurcation theory _
( springer - verlag , wien , new york ) , pp .
52 - 68 .
valsakumar , m. c. & ananthakrishna , g. [ 1983 ] `` a model based on nonlinear oscillations to explain jumps on creep curves : ii .
approximate solutions '' , _ j. phys .
d _ * 16 * , 1055 - 1068 .
venkadesan , s. , murthy , k. p. n. , & valsakumar , m. c. [ 1995 ] `` chaotic flow in al - mg alloy '' in _ solid state phenomena _ vol .
42&43 , eds .
ananthakrishna , g. , kubin , l. p. & martin , g. ( scitec , switzerland ) , pp .
287 - 292 .
zagoruyko , l. n. , osetskiy a. i. & soldatov , v. p. [ 1977
] `` discontinuous deformation of zinc single crystals under creep conditions '' , _ phys . met .
metallogr . _ * 43*(5 ) , 156 - 164 .
in this appendix we give the expressions for the coefficients appearing in the tdgl equation up to the quintic term , i.e. the expressions for @xmath126 , @xmath163 and @xmath164 in the equation @xmath182 .
they are given by    @xmath183    @xmath184    and @xmath185 where @xmath186,\ ] ] and @xmath187 @xmath188 @xmath189 the expressions for @xmath190 s and @xmath191 s used in eq.(a5 ) are given by @xmath192,\ ] ] @xmath193 @xmath194,\ ] ] @xmath195,\ ] ] and @xmath196 @xmath197 @xmath198 @xmath199 where @xmath200    expressions for @xmath160 s , @xmath161 s , @xmath0 s and @xmath162 s are given in appendix c.
in this appendix the equations describing the secular motion are derived for the case when the order parameter , @xmath85 , is found from the quintic tdgl equation ( eq .
44 ) .    from eq .
( 39 ) ,    @xmath201    and @xmath202 where the complex amplitude is @xmath203 with @xmath204 given in terms of functions of @xmath107 and phase factors ( see appendix c ) .
on the other hand , to @xmath114 , @xmath205 which , upon inverting , gives @xmath107 in terms of @xmath85 : @xmath206 substituting the @xmath207-expressions from appendix c in eqs .
( b1 ) and ( b2 ) and , then , using eq .
( b4 ) leads to obtaining the secular equations given by : @xmath208 and @xmath209 where @xmath210    with    @xmath211    where @xmath212 stands for either @xmath213 , @xmath214 , @xmath215 or @xmath216 , while @xmath217 stands for either @xmath160 , @xmath161 , @xmath0 , or @xmath162 correspondingly .
expressions for @xmath213 , @xmath214 , @xmath215 and for @xmath216 can be obtained from @xmath160 s , @xmath161 s , @xmath0 s and @xmath162 s which are given in appendix c.
expressions for @xmath110 , @xmath218 , @xmath219 , @xmath220 ( with a factor @xmath221 ) are    @xmath222    where @xmath223,\\ \ell_{3 } & \equiv &         \frac{i}{2\omega}\,\left[2p_{2}(k_{1}^{\star}+k_{2}^{\star}+k_{2 } +        k_{3})+3p_{3}\right],\\ \ell_{4 } & \equiv &         \frac{i}{4\omega}\,\left[2p_{2}(k_{1}^{\star}+k_{3})+p_{3}\right],\\ m_{1 } & \equiv & -\frac{i}{\omega}\left[-\mu\ , k_{1 } - \mu^{\star}\ ,        k_{3}^{\star } + 2 p_{2 } \ , \ell_{1}^{\star } \right],\\ m_{2 } & \equiv & \frac{i}{\omega}\left[-\mu^{\star}\,(k_{2}+ k_{2}^{\star})+        2 p_{2 } ( \ell_{1}+\ell_{1}^{\star } ) \right],\\ m_{3 } & \equiv & \frac{i}{3\omega}\left[\mu\ , k_{3 } -        \mu^{\star}(k_{1}^\star+ 2k_{3 } ) + 2p_2 \ell_{1}\right],\\ m_{4 } & \equiv & -        \frac{i}{3\omega } \left\{p_{2}\left[(k_{1}+k_{3}^{\star})^{2}+ 2        ( \ell_{2}+\ell_{4}^{\star})\right]+ 3p_{3 } ( k_{1}+k_{3}^{\star } ) + p_{4 }        \right\},\\ m_{5 } & \equiv & -        \frac{i}{\omega}\left\{2p_{2}\left[(k_{1}+k_{3}^{\star})(k_{2}+k_{2}^{\star } )        + \ell_{2}+\ell_{3}^{\star } + \ell_{4}^{\star}\right]+ 3p_{3 }        \left[2(k_{1}+k_{3}^{\star } ) + k_{2 } + k_{2}^{\star}\right]\right .
\nonumber \\   &   &      + 4p_{4}- 2k_{1}\,\eta   \left .
\right\},\\ m_{6 } & \equiv &   \frac{i}{\omega}\left\{p_{2}\left[2(k_{1}+k_{3}^{\star } )                       ( k_{1}^\star+k_{3})+(k_{2}+k_{2}^{\star})^{2 }                        + 2(\ell_{3 } + \ell_{3}^{\star})\right ]   \right .
\nonumber\\ &   &       + 3p_{3}\left[k_{1}+k_{1}^{\star } + 2(k_{2 }                  + k_{2}^{\star})+k_{3}+k_{3}^{\star}\right ]                  + 6p_{4}-2k_{2}\eta_1   \left .
\right\},\\ m_{7 } & \equiv &          \frac{i}{3\omega}\left\{2p_{2}\left[(k_{1}^{\star}+k_{3})(k_{2}^{\star}+k_{2 } )          + \ell_{2}^{\star}+\ell_{3 } + \ell_{4}\right]+ 3p_{3 }         \left[2(k_{1}^{\star}+k_{3 } ) + k_{2 }         + k_{2}^{\star}\right]\right .
\nonumber \\   & &      + 4p_{4}-2k_{3 } \eta^{\star}\left . \right\},\\ m_{8 } & \equiv &          \frac{i}{5\omega}\left\{p_{2}\left[(k_{1}^{\star}+k_{3})^{2 }         + 2(\ell_{2}^{\star } + \ell_{4})\right ]         + 3p_{3}(k_{1}^{\star } + k_{3})+p_{4 } \right\},\\ n_{1 } & \equiv & \frac{i}{2 \omega}\left(\mu - \mu^{\star}\right )         \ell_{1},\\ n_{2 } & \equiv & -\frac{i}{2 \omega } \left\ { 2 p_{2}\left[\ell_{1}^{\star } ( k_{1}+ k_{3}^{\star } ) +         m_{1 } + m_{3}^{\star}\right ]   + 3p_{3}\ , \ell_{1}^{\star } -2\mu\ell_{2 } - \mu^{\star }         \ell^{\star}_{4 } \right\},\\ n_{3 } & \equiv & \frac{i}{2\omega } \left\ { 2 p_{2}\left[\ell_{1}(k_{2}+ k_{2}^{\star } ) +         \ell^{\star}_{1}(k_{1}^{\star } + k_{3 } ) + m_{1}^{\star } + m_{2}^{\star }         + m_{2}+ m_{3}\right]\right .
\nonumber\\ &   &    + 3p_{3}(\ell_{1}^{\star}+2\ell_{1})\left.-\eta^{\star } \ , \ell_{1 } - 2         \mu^{\star}\,\ell_{3 } \right\},\\ n_{4 } & \equiv & \frac{i}{4\omega } \left\ { 2 p_{2}\left[\ell_{1}(k_{1}^{\star}+ k_{3 } ) +        m_{1}^{\star } + m_{3}\right ] + 3p_{3}\ , \ell_{1}-\mu^{\star } \ell_{2}^{\star } +         ( \mu - 3\mu^{\star})\ell_{4 } \right\},\\ n_{5 } & \equiv & -         \frac{i}{4\omega}\left\{2p_{2}\left[m_{4}+m_{8}^{\star } +         ( k_{1}+k_{3}^{\star } )         ( \ell_{2}+\ell_{4}^{\star})\right ] + 3p_{3 }         \left[(k_{1}+k_{3}^{\star})^2 + \ell_{2 } + \ell_{4}^{\star}\right ] \right.\nonumber \\ &       &   + 4p_{4}(k_{1}+k_{3}^{\star } ) + p_{5}\left . \right\},\\ n_{6 } & \equiv & -\frac{i}{2\omega}\left\ { 2p_{2}\left[m_{4 } +        m_{5 } + m_{7}^{\star } +         m_{8}^{\star } + \ell_{3}^{\star}(k_{1}+k_{3}^{\star } ) + ( \ell_{2 }        + \ell_{4}^{\star } ) ( k_{2 } + k_{2}^{\star } ) \right ] \right .
\nonumber \\   &    & + 3p_{3 }         \left[\ell_{3}^{\star } + 2(\ell_{2 } + \ell_{4}^{\star } ) + ( k_{1 } +          k_{3}^{\star})^{2 } + 2(k_{1}+k_{3}^{\star})(k_{2}+k_{2}^{\star } ) \right]\nonumber \\   &   &   + 4p_{4 } \left [ k_{2 } + k_{2}^{\star } + 3(k_{1}+k_{3}^{\star } ) \right ]         + 5 p_{5}\left.-3\eta \ell_{2 } \right\},\\ n_{7 } & \equiv & \frac{i}{2\omega}\left\ { 2p_{2}\left[m_{5}^{\star }   + m_{6}^{\star}+         m_{6 } + m_{7 } + ( \ell_{2}^{\star}+\ell_{4})(k_{1}+k_{3}^{\star } ) + \ell_{3 }         ( k_{2 } + k_{2}^{\star})+\ell_{3}^{\star}(k_{1}^{\star}+k_{3 } ) \right ] \right .
\nonumber\\ &   & + 3p_{3 } \left [ ( k_{2 } + k_{2}^{\star})^{2 } +         2(k_{1}^{\star}+k_{3})(k_{1}+k_{2}+k_{2}^{\star}+k_{3}^{\star})+\ell_{2}^{\star }         + \ell_{3}^{\star } + 2\ell_{3}+\ell_{4}\right ] \nonumber\\ &    & + 4p_{4 } \left [ k_{1 } + k_{3}^{\star } + 3(k_{1}^{\star}+k_{2}^{\star}+k_{2 }         + k_{3})\right ] + 10 p_{5}\left.-(\eta          + 2\eta^{\star})\ell_{3 } \right\},\\ n_{8 } & \equiv & \frac{i}{4\omega}\left\{2p_{2}\left[m_{4}^{\star } + m_{5}^{\star } + m_{7 } +         m_{8 } + \ell_{3}(k_{1}^{\star}+k_{3 } ) + ( \ell_{2}^{\star }         + \ell_{4 } ) ( k_{2 } + k_{2}^{\star } ) \right ] \right .
\nonumber \\ &   &   + 3p_{3}\left [   ( k_{1}^{\star } +        k_{3})^{2 } + 2(k_{2}+k_{2}^{\star})(k_{1}^\star+k_{3})+2 ( \ell_{2}^{\star } +        \ell_{4 } ) + \ell_{3}\right ] \nonumber \\ &   & + 4p_{4 } \left [ k_{2 } + k_{2}^{\star } + 3(k_{1}^{\star}+k_{3 } ) \right ]        + 5 p_{5}\left.-3\eta^{\star}\ell_{4 }          \right\},\\ n_{9 } & \equiv & \frac{i}{6\omega}\left\{2p_{2}\left[m_{4}^{\star } +        m_{8 } +   ( \ell_{2}^{\star }        + \ell_{4 } ) ( k_{1}^{\star } + k_{3}\right ] + 3p_{3 }        \left [   ( k_{1}^{\star } +         k_{3})^{2 } + \ell_{4}^{\star } +         \ell_{4}\right ] \right.\nonumber \\ &     &   + 4p_{4 } ( k_{1}^{\star}+k_{3 } ) +   p_{5}\left . \right\}.\end{aligned}\ ] ]    note that expressions for @xmath224, ... ,@xmath225 are given in eq .
( a5 ) in appendix a.
fig . 1 : :    plot of @xmath226 versus @xmath33 with    @xmath227 .
the hatched portion shows the instability    region which extends between @xmath73 and    @xmath74 approximately .
fig . 2 : :    plots of @xmath127 and @xmath228 versus    @xmath33 ( @xmath227 ) .
note that    @xmath127 is negative ( but finite ) in the interval    @xmath229 .
similarly , @xmath228 is    negative ( but finite ) in the interval    @xmath134 .
fig . 3 : :    plots of @xmath230 and @xmath231 versus    @xmath33 in the region of supercritical bifurcation . fig . 4 : :    plots of @xmath230 and @xmath231 versus    @xmath33 in a subregion of subcritical bifurcation , i. e. the    region where quintic tdgl equation holds .
fig . 5 : :    plots of the ( numerical ) limit cycle solutions obtained from the    secular equations ( dotted line ) , the reduced model ( dashed line ) and    the @xmath153-truncated model ( crossed line ) for    @xmath157 and @xmath158 . fig . 6
: :    plots of the ( numerical ) limit cycle solutions obtained from the    secular equations ( dotted line ) , the reduced model ( dashed line ) and    the @xmath155-truncated model ( crossed line ) for    @xmath159 and @xmath158 . fig . 7 : :    plots of the ( numerical ) limit cycle solutions obtained from the    reduced model ( dashed line ) and the @xmath168-truncated model    ( crossed line ) for @xmath170 and    @xmath158 .
: :    plot of the ( numerical ) limit cycle solution obtained from the reduced    model for @xmath177 and @xmath158 .