catastrophes as defined in the abstract are related to a class of phenomena sometimes called `` relaxation '' oscillations .
the observation of two widely separated time scales makes relaxation oscillations a good a priori subject of theoretical investigation , because one may suspect that their formation results from the existence of a ( more or less hidden ) small parameter .
this gives some hope of a `` general '' theory based upon the small size of the parameter .
to take an example , some major earthquakes lasting a few tens of seconds occur in the same general area about every hundred years .
they thus involve a ratio of ` typical ' times in the neighborhood of @xmath0 , which is a very small number . in the present study
a laboratory experiment was devised which displays similar relaxation oscillations , but acting time - wise on a scale that is convenient for investigation .
experimental observations are explained in light of a two - part _ dynamical bifurcation _ model of catastrophe . comprising the local form of a physical model shown to be valid when close to the catastrophe , the striking result is as follows . in response to an external source of noise the two models ( the local one and the physical one ) predict fluctuations with a correlation time that increases before the catastrophe , just as is observed experimentally .
the important point is that this well known critical slowing down phenomenon occurs significantly _ before _ the transition and could be used to forewarn it .
we study the creeping of a soft metal under constant stress .
accurate time records show that this creeping actually displays the following time dependent noisy component .
the wire typically lengthens slowly with background ( non - thermal ) noise , with sometimes a `` large '' sliding event , followed again by a noisy slow lengthening regime , etc .
plastic deformation of solids is a complex phenomenon , not yet fully understood @xcite .
it has long been observed to take place in a non - smooth manner .
most studies have focused on the portevin - le chatelier effect observed under constant strain rate conditions , whereas the present experiment is concerned with creep at constant stress . in both cases large steps ( cyclic slips )
are embedded in a noisy background .
there is a general agreement that this is due to the complex dynamics of networks of dislocations whose motion is a means for the stressed solid to flow . as reported in section [ sec : experiment ] below , careful observations of the creep in strained sn - pb quasi - eutectic material ( a solder wire ) show the following :    _ i ) _ on average a sample under constant stress lengthens at constant rate .    _
ii ) _ continuous monitoring shows time dependent fluctuations of this length superposed on its secular increase .
_ iii ) _ from time to time the length jumps by steps . afterwards a noisy and steady ( on average ) length increase is recovered until the next jump , etc .    explaining all this remains a challenge for the common models of creep .
nevertheless it is of great interest because it can be seen as a laboratory model of other far less accessible phenomena like earthquakes , where on average there is also a continuous slow sliding with random microseismic noise , interrupted by large fast sliding steps characterizing major earthquakes .
we recently introduced the idea @xcite that in such systems the slow to fast transition can be described by a saddle - node bifurcation in a dynamical system evolving slowly with time .
we pointed out the interest in such a modelization , which allows one to theoretically predict the response of this dynamical system to a noise source .
it was shown that the induced fluctuations drift toward low frequencies ( i.e. , toward large correlation times ) _ before _ the transition , and so could be used as a forewarning .    here
we show that the dynamical saddle - node bifurcation model is the reduced form of a set of equations previously used to describe the portevin - le chatelier effect in metals , or metallic alloys , i.e. , the ananthakrishna ( ak ) model @xcite .
the saddle - node model describes fast transitions resulting from the intrinsic dynamics of the original system which can be described locally ( close to the step ) as a slowly rocking potential system .
using the ak model as well as the saddle - node reduced model we show that , with an added external source of noise , the correlation time of the fluctuations increases before the transition , following the classical scenario of slowing down at bifurcation points . as was shown long ago by dorodnytsin @xcite
, the same local dynamics describe the transition from slow manifold to fast transients in relaxation oscillations of dynamical models , like the van der pol equation in the strongly nonlinear limit . by looking at the fluctuations of the length of our samples
, we found this behavior near the step - like transitions , with the characteristic drift to low frequencies before the transition .    unlike recent publications that have presented the idea @xcite that precursors of earthquakes could be found in the mechanical response to external perturbations ( such as the increase of the fluctuations and their slowing down near transitions ) our idea goes further , by giving an order of magnitude of the precursor time
the authors of @xcite do not introduce the effect of a given time dependence of the parameters and consider only systems with steady parameter values on both sides of the bifurcation . to quote reference @xcite
`` the suggested approach to analytical study of any kind of catastrophes is based essentially on the solution of a stationary problem of the possibility and conditions of the unstable equilibrium state in the system in question '' . without taking into account explicitly the time dependence of the parameters sweeping the bifurcation set ,
it is impossible to get the time scale for predictions . as we show
, this scale depends crucially on the rate of change of the parameters near the bifurcation , which may be estimated from the knowledge of the ratio of the two time scales ( fast and slow ones ) for saddle - node models .
an important challenge lies in the difficulty of stating a suitable model for a given catastrophic event .
actually there is more than one class of possible slow - to - fast transitions in dynamical systems .
the dynamical saddle - node can be seen as belonging to the class of systems with an equilibrium point losing stability as a parameter changes .
it is not a loss of stability , but rather a loss of existence of the equilibrium point , occurring at the folding point of the slow manifold .
however slow - to - fast transition may happen without any folding of the slow manifold .
as shown in @xcite exploring a very often used mathematical representation of stick - slip behavior , the dieterich - ruina equations , the slow - fast transition can originate from the finite time singularity of the slow dynamics itself . in that case , critical speed - up , or drift of fluctuations toward large frequencies , is found to replace the critical slowing - down effect .
this nonlinear phenomenon is obviously outside the class of phenomena explainable by a stability analysis of equilibria of dynamical systems .
we refer the interested reader to the paper on this subject @xcite . in the present paper
we do not consider this case , because it is clearly not the one observed in the creeping experiments .
note that at this time it is unknown if real earthquakes ( as well as other observed catastrophes ) belong to the saddle - node case with a slowing - down near the transition or to the finite time singularity case with speed - up expected .
the creeping experiment is described in section [ sec : experiment ] together with the striking spectral observations .
in the two sections following section [ sec : experiment ] , we present our theoretical approach . we introduce in section [ sec : dynsaddbif ] an abstract dynamical model showing fast jumps . in this model
, one assumes that , as a dynamical system , the jump follows a `` saddle - node '' bifurcation .
there a pair of fixed points , one locally stable , the other locally unstable , merge and disappear as a control parameter changes .
afterwards , the system moves quickly to a new equilibrium state that is at finite distance ( in phase space ) from previous equilibria , whence the jump . after reviewing the saddle - node bifurcation in this light , we assume that the parameter changes with time , namely that the parameter defining the bifurcation is itself a slow function of time .
when , by this change , the parameter crosses the bifurcation value , the dynamical system makes an abrupt transition and jumps `` generically '' from one equilibrium state to another . in section [ sec : ak ] we show that this `` abstract '' model is pertinent for describing the slow - fast transition in the relaxation regime of a set of equations derived by ananthakrishna for creeping in the relaxation regime .
this set of equations describes the dynamics of populations of dislocations in the creeping solid .
the model shows relaxation oscillations in which slow drift is interrupted by fast variations . near the slow - fast transition
, we show that this model reduces to the generic equation quoted above , for a certain range of parameters .
therefore it displays a typical critical slowing - down in its response to external noise , in agreement with the real data reported in section [ sec : experiment ] . of particular interest
is the fact that , from the experimental data , one can predict in advance a `` large slip '' event .
the event is preceded by a shift toward low frequencies in the random fluctuations of specimen length .
choose a uniform``wire '' that creeps under small stress at room temperature , fabricated from a much studied material , sn - pb solder , alloyed to be nearly eutectic .
to establish creep at virtually constant stress , let one end of the wire be fixed and at the other end establish a constant force of tension . by this
means we found the length of the wire to increase on average at constant rate .
by highly accurate monitoring of this length vs time we observed a small time dependent part with the following pattern : on a background of fluctuations , from time to time a large slip event was observed , after which the continuous lengthening with a small background noise was recovered .
the observed noise is non - thermal , since thermal noise has a far too small amplitude to be relevant .
it is assumed to result from rearrangement of defect structures in the poly - crystalline structure of the quasi - eutectic , involving dislocation dynamics .
we assume that this noise originates from a source that is to first - order independent of the overall lengthening of the wire .
this is equivalent to saying that it is due to ongoing micro - scale events triggered by the imposed stress , independent of the global creeping .
therefore we analyzed the response to this noise source according to the ak equations subject to an external noise ( see below ) .
the instrument used in these experiments , which is pictured in fig.[fig : solder ] , is an extensometer @xcite .
young s modulus can be accurately measured with a wire specimen , by placing different size masses on the weight pan .
the trace of the vertical wire holding up the boom / weight pan in fig.[fig : solder ] has been enforced ( colored in black ) to be visible in the image .
though not presently used , the black clamp was for purpose of holding a power resistor that was employed to measure the specimen s thermal coefficient of expansion . to measure temperature changes ,
a solid state thermometer was placed down into the sample space , through the hole seen near the top knurled clamp .
@xmath1   @xmath1    to calibrate the instrument a tungsten wire of diameter 0.1 mm was mounted in the extensometer , and signal output level changes were recorded as various gram - mass standards were placed on the weight pan . by using the known young s modulus for tungsten , the resulting measurements yielded a constant of 1.0 nm per analog to digital count , for the 24 bit adc employed , which is sold by symmetric research @xcite .
this constant is applicable to the data presently reported .
a different measurement technique yielded essentially the same calibration constant .
a he - ne laser was used with a mirror , operating as an optical lever , to measure boom position change as different masses were placed on the pan .    on the basis of two factors ,
ordinary solder was chosen for the present study .
first of all , an unusual property of tin is well known , when an ingot of the metal is strained by large amounts .
the sound which it then emits is well described by the german word  zinngeschrei " , which translates `` tin cries '' .
the initiative for this work was also influenced by the observation of unusual spectral features in the output from a novel seismograph @xcite .
the unusual low - frequency motions of the earth s crust that were then observed to precede an earthquake @xcite are readily seen by the volksmeter .
this is due to the instrument s use of a `` displacement '' sensor , rather than the `` velocity '' sensor used by conventional seismometers .
it was therefore natural to consider an alloy of tin , with the expectation that its defect properties should be more like those of the earth than is possible for a pure metal .
the tin alloy for our study was ordinary soft solder ( 60% sn/ 40% pb ) , used universally in the electronics industry .
although the pure ` eutectic ' alloy is actually 63% sn/37% pb , we will nevertheless use this word to describe our specimen in the discussions that follow .
the stress level due to the load placed on the pb - sn wire used in the present experiments was considerably smaller than the typically @xmath2 used in usual creep studies .
the present load was due solely to the weight of aluminum comprising the boom plus empty weight pan of the extensometer .
this stress value was estimated to be @xmath3 , based on considerations that include the distance of the wire s attachment point from the position of the ( fine ) roller bearing in the unpright housing , that supports the boom on its end opposite the pan .
@xmath1 ( a ) .
( b ) filtered signal @xmath4 at the exit of the single pole high - pass filter ( [ eq : filt ] ) .
, title="fig:",height=144 ] ( b ) .
( b ) filtered signal @xmath4 at the exit of the single pole high - pass filter ( [ eq : filt ] ) .
, title="fig:",height=144 ] @xmath1    during a typical avalanche , the elongation of a @xmath5 long , @xmath6 diameter specimen , would be about @xmath7 . before the avalanche
, typical lengthening velocity is @xmath8 and typical rms fluctuations observed in the wire length ( with secular term removed by high - pass filtering ) would be about @xmath9 .
an example of creep record ( length growth versus time ) is shown in fig.[fig : creep - exp](a ) .
the fluctuations of the raw signal are not visible because they are much smaller than the average length variations .
to eliminate the average growth , we use a standard technique of filtering . at the exit of a low - pass single pole filter , the filtered signal @xmath4 is given by the recursive formula    @xmath10    where @xmath11 stands for the raw creep signal .
in other words the filtered signal is the convolution product of the raw signal derivative by an exponential response function @xmath12 , where @xmath13 is the corner frequency of the filter . in equation ( [ eq : filt ] ) the parameter @xmath14 is given by @xmath15 where @xmath16 is the sampling time of the record , equal to @xmath17 in the experiment . using a corner frequency @xmath18
, the filtered signal takes the form shown in fig.[fig : creep - exp](b ) .
@xmath1 ( a ) before the second step ( the step starts at @xmath19 ) , for the two time intervals ( 1 - 2 ) indicated in figure ( [ fig : creep - exp])-a , at the abscissa @xmath20 and @xmath21 respectively , each for a time interval of @xmath22 counts .
( c - d ) cumulative spectra corresponding to figures ( a - b ) respectively ( e)width @xmath23 of the cumulative spectrum versus time ( at height equal to @xmath24 per cent of the maximum ) , the vertical arrows indicate the very beginning of the fast steps .
, title="fig:",height=115 ] ( b ) before the second step ( the step starts at @xmath19 ) , for the two time intervals ( 1 - 2 ) indicated in figure ( [ fig : creep - exp])-a , at the abscissa @xmath20 and @xmath21 respectively , each for a time interval of @xmath22 counts .
( c - d ) cumulative spectra corresponding to figures ( a - b ) respectively ( e)width @xmath23 of the cumulative spectrum versus time ( at height equal to @xmath24 per cent of the maximum ) , the vertical arrows indicate the very beginning of the fast steps .
, title="fig:",height=115 ] @xmath1    @xmath1 ( c ) before the second step ( the step starts at @xmath19 ) , for the two time intervals ( 1 - 2 ) indicated in figure ( [ fig : creep - exp])-a , at the abscissa @xmath20 and @xmath21 respectively , each for a time interval of @xmath22 counts .
( c - d ) cumulative spectra corresponding to figures ( a - b ) respectively ( e)width @xmath23 of the cumulative spectrum versus time ( at height equal to @xmath24 per cent of the maximum ) , the vertical arrows indicate the very beginning of the fast steps .
, title="fig:",height=115 ] ( d ) before the second step ( the step starts at @xmath19 ) , for the two time intervals ( 1 - 2 ) indicated in figure ( [ fig : creep - exp])-a , at the abscissa @xmath20 and @xmath21 respectively , each for a time interval of @xmath22 counts .
( c - d ) cumulative spectra corresponding to figures ( a - b ) respectively ( e)width @xmath23 of the cumulative spectrum versus time ( at height equal to @xmath24 per cent of the maximum ) , the vertical arrows indicate the very beginning of the fast steps . , title="fig:",height=115 ] @xmath1    @xmath1 ( e ) before the second step ( the step starts at @xmath19 ) , for the two time intervals ( 1 - 2 ) indicated in figure ( [ fig : creep - exp])-a , at the abscissa @xmath20 and @xmath21 respectively , each for a time interval of @xmath22 counts .
( c - d ) cumulative spectra corresponding to figures ( a - b ) respectively ( e)width @xmath23 of the cumulative spectrum versus time ( at height equal to @xmath24 per cent of the maximum ) , the vertical arrows indicate the very beginning of the fast steps .
, title="fig:",height=115 ] @xmath1    in order to analyze the spectral properties of the fluctuations , and see how they evolve with time , we must consider separate sequences of a given time duration .
one significant difficulty in this analysis is the proper choice for the magnitude of this time duration
. the window can not be too large , lest the spectrum becomes invariant with time .
conversely , if the window is too narrow , there will be excessive noise in the spectral ( or correlation ) signal .
we believe the time windows have been judiciously chosen ; so that there is a significant increase of the correlation time before a `` burst '' .
such increase of the correlation time ( actually the decrease of the frequency width of the signal ) was found to happen before every sliding event , giving some hope that this
critical slowing down " is pertinent for predicting the  catastrophe " before it occurs . for the present work a time duration of order @xmath25 of the interval between adjacent bursts was chosen .
the spectral density of a given sequence ( @xmath26 ) of the filtered signal    @xmath27    changes significantly during the creep process .
the striking effect is the shift of the spectral density toward low frequencies in the last stage of the slow regime , i.e. , just before the burst .
examples of this phenomenon are presented in the figures ( [ fig : creep - exp2])(a - b ) , which show the spectra corresponding to the time intervals ( @xmath28 indicated in figure ( [ fig : creep - exp])-a . the two spectra clearly differ .
the first spectral density ( a ) displays a large number of components ; whereas the second spectral density , corresponding to the time interval just before the burst ( figure b ) , is concentrated close to zero frequency . to quantify this effect
, we calculate the cumulative spectrum    @xmath29    which is a smooth curve whose asymptotic value ( for n = n ) gives the experimental variance of the fluctuations during the sequence considered    @xmath30    @xmath1 ( a ) .
the horizontal segments ending at each point indicates the time interval ( before this point ) over which the width is calculated.,title="fig:",height=115 ] ( b ) .
the horizontal segments ending at each point indicates the time interval ( before this point ) over which the width is calculated.,title="fig:",height=115 ] @xmath1    the shift of the density spectrum toward low frequencies is visually more clear in its associated cumulative spectrum , as seen in figures ( c - d ) .
we define the  low frequency extension " @xmath31 of a sequence by the width of the cumulative spectrum corresponding to @xmath24 per cent of its maximum value ,    @xmath32    shown in figure 3(e ) is the evolution of the spectral width for the whole experimental record .
a drop in @xmath31 is clearly seen to occur before each burst .
the important point is that the decrease of the width before the burst is experimentally foreseeable because it occurs during a  precursor " time which is noticeably larger than the duration of a sequence @xmath33 .
this is illustrated in fig.[fig : creep - exp3](b ) which displays a zoom of the spectral width evolution before the burst .
the two following sections are devoted to theoretical models , the first one shows the main features of what is measured in the experiment on creeping , namely the critical slowing down of the fluctuation spectrum occurring before the sliding event , the second one is a model of creeping found in the literature which consists in a set of nonlinear coupled ordinary differential equations , derived by ananthakrishna and collaborators @xcite .
this section explains why a saddle - node bifurcation with a slow sweeping of the bifurcation parameter exhibits a slowing down in its response to a source of noise in a window of time extending well before the bifurcation itself .
abstract model has no direct connection with the physical phenomenon of creeping . in the section
afterwards , however , we explain that a model of creeping shows the same slowing down in a range of parameters , linked to a local saddle - node bifurcation .
consider first the saddle - node bifurcation of a  gradient flow " , that is a damped dynamical system such that a coordinate @xmath34 is a solution of the equation of motion of the form    @xmath35    in this equation @xmath36 is a potential , and the dynamics tends to everywhere lower the value of @xmath36 .
the catastrophe theory of thom and arnold @xcite studies how steady equilibria of equations like ( [ eq : grad ] ) change under smooth deformations of the potential @xmath36 .
below we consider a different kind of question , namely what happens to the solution of equation ( [ eq : grad ] ) when the potential @xmath36 becomes itself a slowly varying function of time , and particularly when a pair of equilibrium points disappears by a saddle - node bifurcation . indeed
this question of the sweeping across bifurcations has been already widely studied @xcite with various applications in mind ( by `` sweeping '' we mean crossing of a transition point with a time dependent parameter in the equation(s ) of motion ) . however , to the best of our knowledge the occurrence of an intermediate time scale in the case of slow sweeping has been overlooked , although we believe it to be crucial for a strategy of foretelling catastrophes in the real world .
the equation ( [ eq : grad ] ) is too general to be very helpful .
however , it may describe a saddle - node bifurcation where a stable equilibrium disappears , assuming that @xmath37 depends slowly on time in a prescribed way , to become a function @xmath38 .
near the transition , one may use a mathematical picture which is correct for a short time around the transition if the potential @xmath38 is a smooth function ( see below for what happens afterwards ) .    assume first that @xmath36
does not depend explicitly on time and takes the form @xmath39    with @xmath40 real constant ( for the moment ) .
for @xmath40 negative @xmath36 has two real extrema ( _ i.e. _ the roots of @xmath41 ) , one @xmath42 is a stable equilibrium , the other , @xmath43 , is an unstable equilibrium . for @xmath44
the two equilibria merge and disappear for @xmath40 positive , see fig.[fig : pot](a ) .
this is the saddle - node bifurcation .
the shape of @xmath36 near @xmath45 and for @xmath40 small is universal : for a given smooth @xmath36 showing this saddle - node bifurcation , one can always rescale @xmath46 and the external parameter to obtain the " local problem in this form .
the extension to a time dependent control parameter @xmath40 goes as follows .
if @xmath40 is a smooth function of time , one can assume that @xmath47 crosses the critical value , _
i.e. _ zero in the present case , at time zero in such a way that @xmath48 with @xmath49 non zero constant , dots being for higher terms in the taylor expansion of @xmath47 . for @xmath50 and @xmath46 close to zero , after rescaling , one can represent the dynamical system ( [ eq : grad ] ) , close to the saddle - node bifurcation , by the `` universal '' parameterless equation    @xmath51    @xmath1 ( a ) .
( b ) quartic potential , @xmath52,title="fig:",height=144 ] ( b ) .
( b ) quartic potential , @xmath52,title="fig:",height=144 ] @xmath1    outside of the neighborhood of @xmath53 , the solution of ( [ eq : grad ] ) depends on other parameters defining @xmath54 for finite values of @xmath46 , as explained below .
although this is not obvious from its formulation , this model is also valid for the transition from slow to fast motion in the van der pol oscillator in the limit of large non linearities , as shown by dorodnitsyn @xcite .
let us sketch the proof of this ( interesting ) point .
actually we shall look at the formally more general situation of relaxation oscillations , namely at solutions of a set of coupled ode s ( ordinary differential equations ) with a large parameter in the form : @xmath55 and @xmath56 in this set of equations dots are for time derivatives and the functions @xmath57 and @xmath58 are smooth with values of order @xmath59 when their argument is also of order one .
moreover , @xmath60 is a large parameter .
the slow manifold is defined by the condition that , in the limit @xmath60 large , the function @xmath61 must be close to zero over at least part of the trajectory .
the cartesian equation @xmath62 defines a curve in the plane @xmath63 which allows to find @xmath64 as a function of @xmath46 , at least locally .
this defines the equation of motion ( along the slow manifold ) @xmath65 where @xmath66 is such that @xmath67 the slow trajectory so defined stops at `` folds '' where the function @xmath66 ceases to be well defined , namely for values of @xmath63 such that @xmath68 , @xmath62 , @xmath69 being not zero .
this defines a discrete set of points .
near those points , the equation of motion can be solved by taking @xmath70 and @xmath71 where @xmath72 are the cartesian coordinates of the point such that    @xmath73    let us look at the solution of the coupled equations ( [ eq : relax1 ] ) , ( [ eq : relax2 ] ) near @xmath72 .
the equation ( [ eq : relax2 ] ) is not singular at this point and so can be solved for small variation of @xmath64 and for @xmath50 small , like @xmath74 , @xmath75 being the ( arbitrary ) time where the trajectory is in @xmath76 and @xmath77 .
consider now the first equation ( [ eq : relax1 ] ) , and expand its right - hand side for @xmath78 and @xmath79 small . because @xmath80 at @xmath72 , the first nontrivial term in the taylor expansion of @xmath81 near @xmath82 with respect to @xmath46 is @xmath83 . on the other hand
the first term coming from the expansion with respect to @xmath79 is @xmath84 .
therefore , for @xmath50 and @xmath78 small , the equation for @xmath85 derived from ( [ eq : relax1 ] ) reads : @xmath86 one can check that all terms not written explicitly there are effectively negligible compared to the ones kept .
standard rescalings allows to transform the equation ( [ eq : relax3 ] ) into the `` universal '' equation ( [ eq : gradt ] ) , provided various constraints of sign are satisfied by the quantities @xmath87 , @xmath88 and @xmath89 all computed at @xmath72 .
this is a short version of the classical calculation by dorodnycsin , showing that the `` universal '' equation ( [ eq : gradt ] ) is also relevant for the slow to fast transitions in relaxation oscillations .
it is worth mentioning also that this derivation does not really assume that the overall motion is periodic , as it may be extended quite easily to a system of many more coupled ode s with only one fast variable .
this shows that such large jumps are also possible in non periodic dynamics , because the dynamics on the slow manifold may be chaotic in between the jumps if this slow manifold has a number of dimensions sufficiently large .
now we shall focus on the slowing down near the dynamical saddle - node bifurcation described by equation ( [ eq : gradt ] ) .
we shall first give its explicit solution .
we look for a solution of equation ( [ eq : gradt ] ) transiting from the `` stable '' fixed point at `` large '' negative times to the rolling down towards positive value of @xmath46 at positive times .
this solution behaves like @xmath90 at large negative times .
the equation ( [ eq : gradt ] ) is of the riccati type and can be integrated by introducing the function @xmath91 such that @xmath92 where @xmath93 and @xmath91 is a solution of airy s equation @xmath94    the solution of equation ( [ eq : gradt ] ) relevant with the given condition at @xmath95 is drawn on figure([fig : quartic]-a ) . in terms of the variable @xmath91
it is the airy function @xmath96 which writes    @xmath97    yet we have only solved the transient problem near the saddle - node bifurcation .
the transition ends - up when @xmath50 becomes equal to the first zero of the airy function @xmath96 , _ i.e. _ the smallest root of the equation    @xmath98 , a pure number , about @xmath99 .
it corresponds to a divergence of @xmath100 , which behaves as    @xmath101    just before this transition , as derived by expanding @xmath102 close to @xmath103 .
therefore the `` generic '' equation ( [ eq : gradt ] ) for the dynamical saddle - node bifurcation displays a finite time singularity .
let us precise the following mathematical subtlety .
this property of the _ local flow _ , which results from the folding of the slow manifold , differs qualitatively from the finite time singularity found in the dieterich - ruina equations @xcite , where it was a property of the flow _ reduced to the slow manifold _ which is everywhere convex , as discussed in the introduction .
in the case of the dynamical saddle - node bifurcation the singularity requires one to consider both the dynamics on and off the slow manifold , and this happens because the geometry forbids the continuation of an exact trajectory on this folded slow manifold .
actually the solution ( [ eq : singdn ] ) loses its physical meaning sometime before the singularity since the `` universal '' dynamical equation ( [ eq : gradt ] ) was derived under the assumption that @xmath46 remains close to zero .
this local theory can not deal with finite variations away from the critical values , therefore we shall need to add finite amplitude effects to limit the growth of the instability after the transition , see the subsection [ sec : quarticpot ] .
we shall study now two questions , first the response of this dynamical system to an external noise , then the dynamics of a system showing a saddle - node bifurcation of the type just studied and reaching a new stable fixed point after this bifurcation .
we explore first the response of our system to a small external noise , and look for qualitative changes in this response which could be a signal that occurs before the transition .
let us consider the equation ( [ eq : gradt ] ) with a small noise added , so that equation ( [ eq : gradt ] ) is replaced by    @xmath104    where @xmath105 is a random function of time , and @xmath106 a small coefficient .    in the limit @xmath106 small , one can solve equation ( [ eq : gradtn ] ) by expansion in powers of @xmath106 , @xmath107    where    @xmath108    the linear response to the noise is @xmath109 because @xmath110 tends rapidly to zero as @xmath111 tends to minus infinity
, one can take @xmath112 to get rid of the effect of the initial conditions .    to make the developments above more concrete ,
let us take a delta - correlated ( or white ) noise , such that @xmath113    the pair correlation of @xmath114 is given by @xmath115 where @xmath116 is the smallest of the two real numbers @xmath50 and @xmath117 .
the behavior of this pair correlation for large negative values of both @xmath50 and @xmath117 , is derived from the asymptotic expression of airy s function , @xmath118 .
setting @xmath119 , and @xmath120 , the variance of the fluctuations writes    @xmath121    in the limit @xmath122 the integral is concentrated near @xmath123 so that    @xmath124    which shows that the amplitude of the fluctuations increases some time before the transition itself .
as the transition approaches , the standard deviation of the fluctuations , @xmath125 increases close to the critical time @xmath103 , because @xmath126 .      as the zeroth order solution diverges at @xmath127
, it does not make sense to describe the dynamical behavior of the fluctuations due to the external noise very close to @xmath103 , as shown above .
as said before , this unbounded growth of the fluctuations is a consequence of the _ local _ cubic form of @xmath36 when expanded near @xmath45 , in obvious contradiction with the fact that @xmath34 tends to infinity . to suppress the divergence of @xmath34 after the saddle - node bifurcation
we add a stabilizing ( positive ) term to the potential @xmath36 which becomes quartic ,    @xmath128    as drawn in fig.[fig : pot](b ) . because of the growth of @xmath129 at infinity , like @xmath130
, the solution of the differential equation @xmath131 does not diverge at finite time .
notice that , formally this does not apply to cases where the local equation ( [ eq : gradt ] ) describes a slow - to - fast transition in a limit cycle , instead of a gradient flow dynamics .
however , unless one insists to look at what happens well after the saddle - node bifurcation , the details of the dynamics after the fast slide are not significant , because one can consider that the point of landing on the slow manifold after the fast drift is like a new equilibrium for a gradient flow , this neglecting the slow motion on this manifold .    the equation ( [ eq : dynquartic ] ) can be written in the given scaled form for any quartic potential provided the coefficient of @xmath130 is positive .
for such a potential one parameter only remains . in equation ( [ eq : quartic ] ) , we choosed to keep as explicit coefficient the coefficient @xmath40 of the linear term in equation ( [ eq : quartic ] ) . for @xmath132 the dynamical system ( [ eq : dynquartic ] )
is exactly at the saddle - node bifurcation , because at @xmath133 both the first and second derivative of @xmath129 vanish , but not the third derivative .
contrary to the case of the pure cubic potential , this system has always , that is for any value of @xmath40 , a stable fixed point beyond the pair of fixed points collapsing at the saddle - node bifurcation .
this makes it a fair candidate for describing the dynamical saddle - node bifurcation without blow - up .    as in the previous case
, we shall take now a time dependent @xmath40 , that will be taken as @xmath134 with @xmath49 positive constant . because of the rescaling of the cubic and quartic term , the parameter @xmath49 can not be eliminated ( another possibility would be to put a parameter in front of the cubic term ) . for the potential @xmath135
we shall analyze the solution of the dynamical equation @xmath136 tending at large negative and positive times to the equilibrium point @xmath137 , @xmath50 being considered as a parameter , see fig.[fig : quartic](a ) .
moreover we consider the limit @xmath49 small , @xmath49 being related to the ratio of small to large time scales , it can be estimated from experimental data .
we shall prove that in this limit @xmath49 small , there are three characteristic time intervals , depending how @xmath46 is close to zero .
the long time scale is the time lapse between successive major slips , typically of order @xmath138 in our experiments on creeping . in our model it is the time needed for the potential @xmath139 to change significantly , to move from a pair of fixed points to a saddle - node bifurcation . because time enters in @xmath139 through the combination @xmath140 ,
the a - dimensional time needed for a change of shape of @xmath141 is of order @xmath142    the short time @xmath143 is of order unity in units of our model equation ( [ eq : dynquarticscaling ] ) as shown in the next paragraph .
it is the duration of the abrupt change of slope of the function of time @xmath34 , the rising of @xmath34 at @xmath103 . in the creep experiment
reported above , it is about @xmath144 .
therefore the ratio of these two time scales @xmath145 is as small as @xmath146 in our experiments .
there is another time scale , @xmath147 , the time interval standing before the transition , and close to it , during which the potential is very flat , while the solution has not yet jumped . during this time , @xmath46 and @xmath148
are much smaller than unity , then the cubic term on the right - hand side of equation ( [ eq : dynquarticscaling ] ) is negligible . in this range one recovers the universal equation of the dynamical saddle - node bifurcation ( [ eq : gradt ] ) by taking @xmath149 and @xmath150 , with the boundary condition @xmath151 at @xmath152 tending to minus infinity .
this property concerns the rectangular domain drawn on fig .
[ fig : quartic]-(a ) , where @xmath46 is small , @xmath153 , and @xmath50 extends from @xmath154 to @xmath155 , located before the abrupt increase .
therefore the time extension of this domain introduces the intermediate time scale , @xmath156 long compared to the short time @xmath143 ( which is of order unity , see below ) and small compared to @xmath157 the average time between slips .
we now show that the short time is of order unity , by matching the solution @xmath158 of the universal equation to the solution of equation ( [ eq : dynquarticlate ] ) below , in the vicinity of the critical point @xmath159 .
because @xmath158 behaves like@xmath160 before it diverges , it follows that the solution @xmath34 behaves as @xmath161 for `` large '' values of @xmath162 before the critical time .
this becomes of order one when @xmath163 becomes also of order one .
when this happens , the term @xmath148 in equation ( [ eq : dynquarticscaling ] ) is negligible , therefore the solution of this equation , which can be matched with the solution near the bifurcation , is the solution of the integrable equation @xmath164 with the asymptotic behavior for very large negative times @xmath165 .
this equation shows that , in our model , the fast time scale is of order one , because it has no explicit dependence with respect to the small parameter @xmath49 .
this result is confirmed by the numerics , see fig.[fig : quartic](a ) .
more precisely defining the rising time of @xmath34 by @xmath166 of the width of the time derivative @xmath167 , we get @xmath168 independently of the value of @xmath49 , for @xmath49 small . in the experiment ,
the intermediate time scale is @xmath169 , which gives about @xmath170 for @xmath171 , @xmath172 .    in summary , by matching the two solutions in the range @xmath173
, we show that the catastrophe takes place during this time @xmath143 which is of order one , because the displacement is then of order one , compared to the small displacement of order @xmath174 taking place during time @xmath175 typical of the `` universal '' transition process .    from this understanding of the various scales in the deterministic part of the dynamical equations , we can now look at the response to noise of this system , particularly at the range of time where something like a `` critical slowing down '' could be observed , and which actually happens in our experiments .    with a noise source added ,
the dynamical equation ( [ eq : grad ] ) becomes , @xmath176    actually the effective noise amplitude is not equal to @xmath106 close to the saddle - node , as it depends on the value of the parameter @xmath49 . indeed for @xmath177 , the cubic term in equation ( [ eq : xxquartic ] ) is negligible , and the equation reduces to @xmath178 which may be written as @xmath179 , by setting @xmath180 , @xmath181 , and @xmath182 . therefore the effective noise is larger than @xmath106 , by a factor @xmath183 in the rectangular domain of fig.[fig : quartic ] .
@xmath1 ( a ) ) , with and without noise , @xmath184 ( smooth curve ) and large amplitude noise @xmath185 .
( b ) solution of equation ( [ eq : dynquarticscaling ] ) for @xmath171 .
the rectangle around the origin defines the region @xmath186 and @xmath187 , with @xmath188 .
the critical time is @xmath189 @xmath190 .
the two vertical lines inserted between the two arrows delimitate the large slope time duration , of order unity .
( c ) solution of equation ( [ eq : xxquartic ] ) with large amplitude noise @xmath185,title="fig:",height=120 ] ( b ) ) , with and without noise , @xmath184 ( smooth curve ) and large amplitude noise @xmath185 .
( b ) solution of equation ( [ eq : dynquarticscaling ] ) for @xmath171 .
the rectangle around the origin defines the region @xmath186 and @xmath187 , with @xmath188 .
the critical time is @xmath189 @xmath190 .
the two vertical lines inserted between the two arrows delimitate the large slope time duration , of order unity .
( c ) solution of equation ( [ eq : xxquartic ] ) with large amplitude noise @xmath185,title="fig:",height=144 ] ( c ) ) , with and without noise , @xmath184 ( smooth curve ) and large amplitude noise @xmath185 .
( b ) solution of equation ( [ eq : dynquarticscaling ] ) for @xmath171 .
the rectangle around the origin defines the region @xmath186 and @xmath187 , with @xmath188 . the critical time is @xmath189 @xmath190 .
the two vertical lines inserted between the two arrows delimitate the large slope time duration , of order unity .
( c ) solution of equation ( [ eq : xxquartic ] ) with large amplitude noise @xmath185,title="fig:",height=120]@xmath1    .
[ fig : quartic ]    let us consider now the fluctuations of the solution @xmath34 of equation ( [ eq : xxquartic ] ) . for a small noise source , the solution may be expanded in power of @xmath106 as above . at first order it gives @xmath191 x_1(t ) + \xi(t ) \mathrm { , } \label{eq : x1}\ ] ] whose solution is formally @xmath192 \mathrm{. } \label{eq : x1i}\ ] ] in general @xmath193 is the time integral of the second derivative of the potential @xmath194 , which yields with our choice of @xmath54 : @xmath195 \mathrm{.}\ ] ] the standard deviation @xmath196 has to be calculated numerically .
we expect it to display the same behavior as for the cubic case in the whole domain where @xmath197 , i.e. a little before the transition and close to it , because the potential is cubic in this range .
after the transition , we expect that the fluctuation decreases , because the solution without noise becomes quasi - steady .
this is confirmed by the numerics : the amplitude of the fluctuations strongly increases close to the critical time @xmath103 , its maximum occurring at time @xmath103 , then it decreases .
more precisely the standard deviation behaves exactly as @xmath198 , the red curve in fig.([fig : width](a ) , for small noise .
therefore the strong increase of the variance of the signal fluctuations can not be used as a precursor for predicting the transition because it occurs simultaneously with the signal itself close to the transition .
note that this observation seems to contradict the currently found statement that fluctuation enhancement precedes the transition and can be used as a precursor . in the case of the saddle - node bifurcation model ,
we have indeed observed a `` precursor '' growth of the fluctuations , but only in the case of `` large '' amplitude noise , see fig.[fig : quartic](b ) .
let us focus on the case of small noise .
in this case we show below that the correlation time of the fluctuations changes much earlier than the onset of amplitude growth .
this slowing - down can be used to foretell the event itself , in considerable advance of its occurrence .
it is seen to consistently happen in the creeping experiment described below .
@xmath1 ( a ) for the saddle - node model ( [ eq : xxquartic ] ) with @xmath171 .
the correlations functions are defined by equation ( [ eq : gamma ] ) and ( [ eq : gammap ] ) for curves ( a ) and ( b ) respectively ; ( c ) spectral width ( a. u. ) calculated with @xmath199 , for @xmath171 .
( d ) @xmath200 for @xmath171 ( blue),for @xmath201 ( red ) .
, title="fig:",height=153 ] ( b ) for the saddle - node model ( [ eq : xxquartic ] ) with @xmath171 .
the correlations functions are defined by equation ( [ eq : gamma ] ) and ( [ eq : gammap ] ) for curves ( a ) and ( b ) respectively ; ( c ) spectral width ( a. u. ) calculated with @xmath199 , for @xmath171 .
( d ) @xmath200 for @xmath171 ( blue),for @xmath201 ( red ) .
, title="fig:",height=153]@xmath1    @xmath1 ( c ) for the saddle - node model ( [ eq : xxquartic ] ) with @xmath171 .
the correlations functions are defined by equation ( [ eq : gamma ] ) and ( [ eq : gammap ] ) for curves ( a ) and ( b ) respectively ; ( c ) spectral width ( a. u. ) calculated with @xmath199 , for @xmath171 .
( d ) @xmath200 for @xmath171 ( blue),for @xmath201 ( red ) .
, title="fig:",height=168 ] ( d ) for the saddle - node model ( [ eq : xxquartic ] ) with @xmath171 .
the correlations functions are defined by equation ( [ eq : gamma ] ) and ( [ eq : gammap ] ) for curves ( a ) and ( b ) respectively ; ( c ) spectral width ( a. u. ) calculated with @xmath199 , for @xmath171 .
( d ) @xmath200 for @xmath171 ( blue),for @xmath201 ( red ) .
, title="fig:",height=144 ] @xmath1    consider the case of small effective noise , where the correlation function and the spectrum of the fluctuations @xmath202 are well described by the correlation function and spectrum of @xmath114 , respectively .
as noted in the previous section , the calculation of these functions requires some care ; because the system is not in a statistically steady state .
therefore the spectral density of the fluctuations depends on time and the correlation function @xmath203 depends on both @xmath50 and on the time difference @xmath204 .
more precisely , the correct definition of the correlation function is actually given by @xmath205 the latter definition of correlation is not readily accessible in experimental situations , because it requires knowledge of the time dependent variance , which is difficult to estimate from a single sample .
a more accessible tool is often used ; it is given by the expression @xmath206 which coincides with the correct expression if the variance is the same at time @xmath50 and @xmath207 only .
if the variance changes noticeably during the time interval of duration @xmath204 , the latter expression is biased , since @xmath208 .
we expect such a discrepancy to manifest itself close to the catastrophe , since the variance increases there by a large amount . to illustrate this point we show in fig.[fig : width](a - b ) the correlation functions of the fluctuation @xmath209 for the solution of the saddle - node model , as defined by equations ( [ eq : gamma ] ) and ( [ eq : gammap ] ) respectively . the strong increase of the correlation time before the catastrophe is noticeably truncated when using the biased expression ( [ eq : gammap ] ) . using the correct definition of the correlation function ( [ eq : gamma ] ) , the correlation time increases by a factor of ten over a time interval of order @xmath147 before the catastrophe , while the enhancement is only about @xmath210 when using the biased expression ( [ eq : gammap ] ) .
moreover the enhancement is followed by a drop in the latter case .
nevertheless both curves display well the critical slowing down effect , which manifests as a growth of the correlation time close to the transition .
the increase of the correlation length _ before _
the catastrophe is understood by looking at the formal expression ( [ eq : x1i ] ) .
the second derivative of the potential vanishes at @xmath211 , and that leads to the flatness of @xmath193 in the time domain @xmath212 , as shown in fig.[fig : width](d ) .
this time domain could therefore be identified as a `` precursor time '' , of order few @xmath147 .
now let us consider the spectrum of the fluctuations in order to compare with the experimental results . a time dependent spectrum can be defined formally by the ( real ) wigner transform @xmath213 that has to be modified for numerical applications , either by using a filtering procedure like the one used in the previous section or by introducing a slipping window . in this section we use the latter process .
choosing a gaussian window function of width @xmath214 , the numerical spectrum is given by @xmath215 where @xmath216 are the numerical integration time boundaries .
this expression does not take into account the variability of the variance ; therefore we expect the observed change in the spectral properties to be biased , as was the case for the correlation function @xmath217 .
the evolution of the spectral width @xmath218 ( half - height width ) is reported in fig.[fig : width](c ) in a range of a few times @xmath147 around the transition .
the solution @xmath219 is hightlighted by the solid red line .
the ( c ) part of the figure illustrates the same effect as the correlation function @xmath217 , but in the fourier space .
the spectral width decreases noticeably from negative time of order few @xmath147 , until the time @xmath220 , where it grows , because of the bias due to the variance increase close to the burst .
the decrease of @xmath218 corresponds to a shift of the spectrum towards low frequencies .
the important result is that this shift occurs well _ before the transition_. it occurs over a time interval of order few @xmath147 .
this result agrees with our experiment where @xmath147 was estimated to be about @xmath221 and the decrease of the spectral width keeps on for about @xmath222 ( see fig.[fig : creep - exp2 ] where @xmath223 corresponds to @xmath224 counts in the abscissa scale ) .
the growth of the fluctuations and their shift to lower frequencies can be understood as follows .
as the transition approaches , the potential @xmath38 becomes flatter and flatter , making weaker and weaker the restoring force toward equilibrium .
therefore , at constant noise source , the amplitude of the fluctuations driven by this noise source will grow because the damping is ever less efficient . moreover , the typical time scale for this damping will get increasingly larger because of the decreasing stiffness of the potential , thus favoring noise at lower and lower frequencies .
this section is devoted to an analysis of solutions of a set of equations devised for describing creeping in solids .
more precisely our purpose is to show that , in a range of parameters this equation exhibits a dynamical saddle node bifurcation . as in the  abstract " model of the previous section , this bifurcation
is also preceded by a slowing - down of the fluctuations triggered by an external source of noise .
creeping phenomena in real materials are complex and difficult to predict quantitatively , despite decades of efforts on theoretical models .
we have chosen to consider a model developed by ananthakrishna for creeping in strained solids . to make things simpler ,
we have only used its version without space dependence in the quantities involved . note that the introduction of space variables would lead to an aperiodic creep signal more realistic than the periodic signals of the present model
, however we conjecture that it should not affect the main result of our study ( the emphasis of a precursor signal over a given time interval ) .
the ak model considered here is a set of three coupled non linear ordinary equations with three dimensionless parameters ( @xmath49 , @xmath40 and @xmath225 , where the letters @xmath49 and @xmath40 have no connection with the same symbols used previously ) . the unknown time dependent quantities are three scaled variables corresponding to three density types of dislocations , @xmath34 , @xmath91 and @xmath226 , representing respectively mobile , immobile , and those with clouds of solute atoms that mimic cottrell s idea .
the model equations write ,    @xmath227    @xmath228    @xmath229    where the variable @xmath34 stands for the elongation rate .
as one can check , if the three variables are positive at time zero , they remain so at later times if @xmath49 , @xmath40 and @xmath225 are positive , as assumed .
the relative elongation ( or strain , or creep ) @xmath230 of the solder wire is the time integral of @xmath34 , @xmath231 .
solutions of equations ( [ eq : x1])-([eq : z1 ] ) have been extensively studied @xcite , their shape and duration versus the parameter values are given in @xcite .
recall that relaxation oscillations are depicted for small values of the parameter @xmath40 only .
for @xmath225 larger than a certain critical value @xmath232 , depending on @xmath233 and not written here , the solution is stable .
the fixed point coordinates are given by the expressions    @xmath234    at lowest order for the small parameter @xmath40 .
this fixed point becomes unstable by a poincar - andronov bifurcation for small values of @xmath40 , only under the condition that @xmath225 becomes smaller than @xmath232 . in a large range of parameter values
the limit cycle associated with the variable @xmath34 displays relaxation oscillations characterized by slow steps and fast bursts .    in this range of @xmath40 and @xmath225 small
, we focus on the case @xmath235 , for two reasons . first because it allows the set of equations to be reduced to the generic saddle - node equation introduced in the previous section . secondly because it leads to @xmath34 solutions looking approximately like our experimental data ( the strain rate @xmath34 has to be compared with the data files labeled @xmath4 in section [ sec : experiment ] ) .
an example of the numerical solution of equations ( [ eq : x1])-([eq : z1 ] ) is given in fig.[fig:2droite ] .
in this case the burst amplitude drawn in figure ( b ) ) is noticeably larger than the experimental one shown in fig.[fig : creep - exp]-(a ) , nevertheless we have chosen the parameter values for pedagogical reasons , in order to give a clear representation of the full cycle of the relaxation oscillations . consider the evolution of the elongation rate @xmath34 .
the limit cycle of duration @xmath236 displays four stages . in the first step , @xmath237
, the flow evolves slowly and @xmath34 takes values of order unity .
this stage is followed by a fast jump of @xmath34 at @xmath238 , then by a short slow stage @xmath239 with high values of @xmath34 , of order @xmath240 , finally followed by a fast decrease of @xmath34 at @xmath241 .
@xmath1 ( a))-([eq : z1 ] ) for @xmath242 , @xmath243 , and @xmath244 .
the three curves are @xmath245 ( blue ) , @xmath91 ( purple ) and @xmath246 ( yellow ) , ( b ) creep signal @xmath230 . , title="fig:",height=144 ] ( b))-([eq : z1 ] ) for @xmath242 , @xmath243 , and @xmath244 .
the three curves are @xmath245 ( blue ) , @xmath91 ( purple ) and @xmath246 ( yellow ) , ( b ) creep signal @xmath230 .
, title="fig:",height=144]@xmath1      we first consider the linear stability of the 3d flow , especially along the first stage described just above , which precedes the fast jump . in the next subsection we derive the local form of the model close to the burst
. we shall prove that both approaches provide the same information , namely an estimation of the precursor time value . in the slow regime , assuming a 3d flow of the form @xmath247 ( and similar expressions for @xmath248 ) close to the trajectory @xmath249 , the exponents @xmath250 are the eigenvalues of the jacobian matrix @xmath251    or solutions of the equation    @xmath252    with @xmath253 , @xmath254 , and @xmath255 .
@xmath1 ( a ) ( b)@xmath1    one of the eigenvalues is real and negative all along the trajectory , while the other two become complex conjugates in the slow regime , their real part crossing zero before the burst at time @xmath256 , see fig.[fig : lyap ] .
the time interval @xmath257 depends on the values of the parameters @xmath258 .
it is generally a small fraction of the limit cycle period .
below it is shown that the numerical value of @xmath257 is nearly equal to the intermediate time scale @xmath147 for the model .
we now consider the behavior of the 3d flow in the vicinity of the burst , and derive the normal form of the ak model close to @xmath259 .
we focus on the first part of the limit cycle , preceding the burst of @xmath34 , and try to understand how the trajectory leaves the slow manifold ( sm ) that we define now . because @xmath260 , the slow stages are described by canceling the right - hand side of equations ( [ eq : x1])-([eq : y1 ] ) , that reduces the 3d flow to the 1d flow    @xmath261    from the first equation the variables @xmath34 is an explicit function of @xmath91 , @xmath262 inserting this expression into the second equation ( [ eq : sm ] ) , allows to also express the variable @xmath226 in terms of @xmath91 ,    @xmath263    in the phase space @xmath264 expression ( [ eq : zs ] ) defines the slow manifold , which is illustrated in fig .
( [ fig : nullcline ] ) , red curve . on the positive slope parts of the slow manifold the 1d flow obeys the differential equation @xmath265 , which becomes singular at critical points defined ( on the sm ) by the relation @xmath266    the 3d flow ( closed blue curve ) follows the path @xmath267 . at the critical point @xmath259
, the trajectory leaves the slow manifold , jumps to point @xmath268 ( fast stage @xmath269 ) , then it follows the portion ( @xmath270 ) , and finally returns to the sm in @xmath271 .
we consider below the exit of the sm close to the critical point @xmath259 .
( red curve ) for @xmath242 , @xmath272 , and parametric plot of the 3d flow ( blue curve , numerical solution @xmath273 of equations ( [ eq : x1])-([eq : z1 ] ) for @xmath244 .
, height=144 ]      the coordinates of the critical point @xmath259 in the phase space @xmath274 are solutions of equations ( [ eq : xs ] ) and ( [ eq : df ] ) . from equation ( [ eq : xs ] ) we derive an expression for @xmath275 @xmath276    which should be identical to the expression of @xmath277 taken from the relation ( [ eq : df ] ) , @xmath278    these two expressions are identical if @xmath279 is a root of the cubic polynomial equation @xmath280 where @xmath281 .
the positive solution @xmath279 and the corresponding value of @xmath282 and @xmath283 satisfying equations ( [ eq : xs])-([eq : zs ] ) respectively are given at leading order with respect to the small parameter @xmath284 by    @xmath285    valid in the parameter range @xmath286 .
we next derive the normal form describing the dynamical behavior of the solution close to the critical point . assuming that the variable @xmath46 follows adiabatically the variable @xmath64 according to equation ( [ eq : xs ] )
, the original system reduces to @xmath287 defining the dynamics for the portion of the trajectory near @xmath259 , including the burst @xmath269 .
this is confirmed by the numerics : starting from the point b , with numerical initial conditions ( [ eq : turn ] ) , the behavior of the solutions of equation ( [ eq:2eqb ] ) agree well with the 3d flow , see fig . [
fig : closet ] .
@xmath1 ( a ) ) with relations ( [ eq : xs])-([eq : zs ] ) , ( b ) solution of the 3d flow equations ( [ eq : x1])-([eq : z1 ] ) . , title="fig:",height=120 ] ( b ) ) with relations ( [ eq : xs])-([eq : zs ] ) ,
( b ) solution of the 3d flow equations ( [ eq : x1])-([eq : z1 ] ) .
, title="fig:",height=144]@xmath1    since we are interested in the description of the solution before and at the burst , in the intermediate regime where @xmath46 remains close to @xmath279 , we shall pursue our analysis by canceling the terms @xmath288 in ( [ eq : xs ] ) and in @xmath289 , because @xmath290 .
this leads to the system    @xmath291    where @xmath292 .
close to @xmath259 and during the burst , the variable @xmath293 is essentially constant , due to the small value of the ratio @xmath294 , whereas the variable @xmath64 jumps toward smaller values , as shown in fig .
[ fig : nullcline ] .
therefore at leading order with respect to the small parameter @xmath294 , the solution of the last equation ( [ eq:2eqb ] ) is given by    @xmath295    in the vicinity of @xmath259 .
inserting this local solution ( [ eq : zbt ] ) for @xmath226 into the system ( [ eq:2eq ] ) , the exit from the sm is then described by a single equation for the local variation @xmath296 of @xmath64 of the variable @xmath64 , which is of the form    @xmath297    where @xmath298    and @xmath299 sets for the @xmath300 derivative of @xmath301 with respect to the variable @xmath64 , taken at point @xmath259 ( note that @xmath302 at the critical point @xmath259 ) . using the reduced system ( [ eq:2eq ] ) , for the first two derivatives of @xmath301 with respect to @xmath64 we write @xmath303    expanding all expressions close to @xmath259 at leading order with respect to the small parameter @xmath40 , one obtains @xmath304    at this stage we can limit the series expansion in equation ( [ eq : s1 ] ) .
a rough approximation for the series to converge is given by @xmath305 . using expressions ( [ eq : derf2 ] ) and the first equation ( [ eq:2eq ] ) gives the range of variation    @xmath306    in this range the local form of the ak equations close to @xmath259 becomes @xmath307 which is identical to equation ( [ eq : relax3 ] ) .
finally , setting @xmath308 , and @xmath309^{1/3 } \delta y$ ] , the relation ( [ eq : generic ] ) takes the generic form ( [ eq : gradt ] ) without any parameter recently proposed as a possible description of a signal before a catastrophe @xcite .
@xmath310    as shown in section [ sec : dynsaddbif ] the solution of equation ( [ eq : generic ] ) displays an intermediate time scale of order @xmath311 which could be used to predict the burst , because the critical slowing down effect occurs during this time interval . in the ak equations the intermediate time scale , or precursor time , is given by the relation    @xmath312    where the parameter @xmath313 is given by ( [ eq : gamma ] ) .
recall that @xmath314 was shown to be the time delay between the catastrophe and the time @xmath315 at which the instability builds up along the slow manifold .
for the parameter values of the above figures , we have @xmath316 , which agrees well with the lyapunov analysis , see fig .
[ fig : lyap ] , and also with the precursor time deduced from the spectral analysis presented in the next subsection .
in summary we have proven that the ak model , although formally different from the van der pol model , has a normal form close to the critical point that is consistent with the dynamical saddle - mode model equation studied in ( @xcite ) , with an intermediate time scale given by expression ( [ eq : prect ] ) .
consequently , as for the saddle - node model , the ak model should display a response to noise with a strong increase of the correlation time occurring with a few times @xmath147 before the burst .
the response to noise of the system ( [ eq : x1])-([eq : z1 ] ) is studied by setting ( @xmath317 , @xmath318 , @xmath319 ) , where @xmath249 is the solution of the noiseless ak equations , and the vector @xmath320 characterizes the fluctuations of the response to a noise source .
these fluctuations result from the introduction of noise terms ( either multiplicative or additive ) in the original system . in the case of multiplicative noise sources ,
the response of the ak equations is a solution of the system    @xmath321(1+\epsilon_x f_x(t))\\   \dot y ( t ) = [ b x(t)^2-x(t)y(t)-y(t)+az(t)](1 + \epsilon_y f_y(t))\\   \dot z ( t ) = ( c / b)[x(t)-z(t)](1 + \epsilon_z f_z(t ) )   \mathrm{. } \end{array } \right .
\label{eq : eqnoise}\ ] ]    we report below the result of the numerical study of the correlation function given by expressions ( [ eq : gamma ] ) and ( [ eq : gammap ] ) for small amplitude noise . in this case
the variance increases slowly along the sm , following approximately the evolution of the variable @xmath322 , as reported in fig .
[ fig : correl ] ( a ) . figures ( b - c ) displays the evolution of the correlation function before the burst , along the path @xmath323 .
the width of the correlation function @xmath324 clearly increases as the burst is approached , showing an evident critical slowing down effect .
the maximum growth of the correlation time occurs close to the time @xmath325 where the sm becomes linearly unstable , persisting until the burst .
this result is in agreement with the experimental result given in section [ sec : experiment ] showing a shift of the spectrum towards low frequencies before the burst .
@xmath1 ( a ) ( blue curve ) , in log scale .
( b ) correlation functions @xmath326 at times @xmath327 -green , @xmath328-orange and @xmath329-black curve versus @xmath204 .
( c ) half - height - half - with of the correlation function @xmath330 versus time @xmath50 .
the input data are those of figure ( [ fig:2droite ] ) , noise amplitudes are @xmath331 for the three independent noise sources @xmath332 . , title="fig:",height=120 ] ( b ) ( blue curve ) , in log scale .
( b ) correlation functions @xmath326 at times @xmath327 -green , @xmath328-orange and @xmath329-black curve versus @xmath204 .
( c ) half - height - half - with of the correlation function @xmath330 versus time @xmath50 .
the input data are those of figure ( [ fig:2droite ] ) , noise amplitudes are @xmath331 for the three independent noise sources @xmath332 . , title="fig:",height=120 ]
@xmath1    @xmath1 ( c ) ( blue curve ) , in log scale .
( b ) correlation functions @xmath326 at times @xmath327 -green , @xmath328-orange and @xmath329-black curve versus @xmath204 .
( c ) half - height - half - with of the correlation function @xmath330 versus time @xmath50 .
the input data are those of figure ( [ fig:2droite ] ) , noise amplitudes are @xmath331 for the three independent noise sources @xmath332 . , title="fig:",height=144 ]
@xmath1    for noise of additive type the width of the correlation function @xmath326 behaves similarly , as shown in fig .
[ fig : correladd ] . as pointed out in a previous section , using the expression ( [ eq : gammap ] ) to calculate the correlation time of the solution gives a biased result close to the burst where the variance is time dependent . in this case
the increase of the width before the burst is reduced as illustrated in fig .
[ fig : correladd](b , d ) .
@xmath1 ( a ) ( dashed , blue , red curves respectively ) , ( a ) @xmath326 , ( b ) @xmath333 versus @xmath204 .
( c - d ) evolution of the half - width of the correlation functions along the trajectory before the burst , ( c ) half - width of @xmath326 ( d ) half - width of @xmath334 both curves drawn versus @xmath50 .
same parameters and input data as in figures ( [ fig : correl ] ) .
noise amplitudes are @xmath335 for the three independent noise sources @xmath332 .
( e ) evolution of the variance along the path ( @xmath269 ) , title="fig:",height=120 ] ( b ) ( dashed , blue , red curves respectively ) , ( a ) @xmath326 , ( b ) @xmath333 versus @xmath204 .
( c - d ) evolution of the half - width of the correlation functions along the trajectory before the burst , ( c ) half - width of @xmath326 ( d ) half - width of @xmath334 both curves drawn versus @xmath50 .
same parameters and input data as in figures ( [ fig : correl ] ) .
noise amplitudes are @xmath335 for the three independent noise sources @xmath332 .
( e ) evolution of the variance along the path ( @xmath269 ) , title="fig:",height=120 ] @xmath1    @xmath1 ( c ) ( dashed , blue , red curves respectively ) , ( a ) @xmath326 , ( b ) @xmath333 versus @xmath204 .
( c - d ) evolution of the half - width of the correlation functions along the trajectory before the burst , ( c ) half - width of @xmath326 ( d ) half - width of @xmath334 both curves drawn versus @xmath50 .
same parameters and input data as in figures ( [ fig : correl ] ) .
noise amplitudes are @xmath335 for the three independent noise sources @xmath332 .
( e ) evolution of the variance along the path ( @xmath269 ) , title="fig:",height=120 ] ( d ) ( dashed , blue , red curves respectively ) , ( a ) @xmath326 , ( b ) @xmath333 versus @xmath204 .
( c - d ) evolution of the half - width of the correlation functions along the trajectory before the burst , ( c ) half - width of @xmath326 ( d ) half - width of @xmath334 both curves drawn versus @xmath50 .
same parameters and input data as in figures ( [ fig : correl ] ) .
noise amplitudes are @xmath335 for the three independent noise sources @xmath332 .
( e ) evolution of the variance along the path ( @xmath269 ) , title="fig:",height=120 ] ( e ) ( dashed , blue , red curves respectively ) , ( a ) @xmath326 , ( b ) @xmath333 versus @xmath204 .
( c - d ) evolution of the half - width of the correlation functions along the trajectory before the burst , ( c ) half - width of @xmath326 ( d ) half - width of @xmath334 both curves drawn versus @xmath50 .
same parameters and input data as in figures ( [ fig : correl ] ) .
noise amplitudes are @xmath335 for the three independent noise sources @xmath332 .
( e ) evolution of the variance along the path ( @xmath269 ) , title="fig : " ] @xmath1    in order to present the ak model results by using the same analytical tools as used for the experiment , we have calculated the cumulative spectrum of the response @xmath336 , @xmath337 where @xmath338 , is the spectrum of the fluctuations for the response signal sampled during the time interval @xmath339 , in analogy with equations ( [ eq : spec])-([eq : cumul ] ) . the characteristic spectral width @xmath340 defined by the relation ( [ eq : cumwidth ] )
evolves in time as illustrated in fig .
[ fig : specum ] for the ak model with additive noise .
the horizontal segments corresponding to abscissa ( @xmath341 ) , have ordinates @xmath340 , the numerical value of the spectral width calculated during this time interval .
the decrease of @xmath340 close to the burst shows well the expected shift of the spectrum towards low - frequencies , observed also in the experiments .
of the response to additive noise @xmath336 , versus time before the burst , same data as in other figures . , height=144 ]    in conclusion
the ak model displays a range of parameter values where one can see the critical slowing - down observed also in the experiments .
we have studied the ananthakrishna model for the case of small @xmath40 and @xmath225 parameter values , with @xmath342 . in this case
the burst occurs close to a critical point , as noticed in @xcite .
close to this critical point , the 3d flow ( namely a set of _ three _ coupled ode s ) can be reduced to the 2d system ( [ eq:2eq ] ) with a sort of langevin - like source term ( noise term with small amplitude noise)of the form ( [ eq : relax1])-([eq : relax2 ] ) discussed in section [ sec : dynsaddbif ] .
the present system differs from the van der pol equation ; however , close to the critical point both models take the form of the dynamical saddle - mode model equation studied in @xcite .
the ananthakrishna model displays a response to noise with a strong increase of the correlation time occurring close to the burst , or a shift of the spectrum towards low - frequency components . for this model
the  precursor time " is a time interval of order several times @xmath343 , also equal to several times the intermediate time scale which depends on the values of the parameters @xmath258 as given by equation ( [ eq : scgeneric ] ) .
a detailed quantitative comparison between the ak model and the creep experiment will require a more complete study , which is beyond the scope of the present work .
in particular it should be noted that including spatial variables in the ak model leads to chaotic solutions , which is in better agreement with our experiment where the limit cycle period may indeed vary by a factor of @xmath210 from one recorded data set to another .
we have shown that the dynamical model of saddle - node transitions recently proposed to foretell catastrophes is applicable to describe the physics of collective dislocations .
the experimental signal of the plastic deformation of the eutectic mixture of sn - pb subjected to a constant stress presented above clearly displays the well - known critical slowing down effect , with a precursor time of order @xmath344 of the relaxation oscillation period .
this observation is shown to agree with the response to noise of a physical model proposed by ananthakrishna for describing creeping in ductile materials .
this physical model has @xmath210 parameters , therefore it would be quasi - impossible to make a quantitative comparison between the experimental data and the theoretical model . here
we show the ak model has a range of parameter values for which the slow - fast transition is described by the dynamical saddle - node model mentioned above .
when adding a small noise to the ak model , we show that it indeed displays the critical slowing down scenario , and we derive the expression for the precursor time in terms of the @xmath210 parameter values .
this gives a relation between the parameter values of the ak model in which the precursor time ( i.e. the time during which the most intense spectral components clearly shift toward low frequencies ) is of order one tenth of the relaxation oscillation period , as observed in the experiment .
the important result is the following : for such systems which display relaxation oscillations with saddle - node transition , we show that there exists an intermediate time scale between the slow and fast regimes , which can be used to foretell catastrophes . unlike the majority of precursor tools that have focused on signals in the time domain , we have looked at fluctuations in the frequency domain .
these spectral precursors are prone to fewer false alarms .
g. ananthakrishna , physics reports * 440*,113 ( 2007 ) .
y. pomeau and m. le berre , http://arxiv.org/abs/1102.5637 .
1 ) ; y. pomeau , m. le berre , j .- l
. le moul , c. narteau et p. fromy , _ comptes - rendus de la 14e rencontre du non linaire _ , dits par c. josserand , m. lefranc et c. letellier , non - linaire publications , saint - etienne du rouvray ( 2011),p.135 .
isbn : 978 - 2 - 9538596 - 0 - 7 .
dubrovskiy and v.n .
sergeev , geos , moscow . * 1 * , 222 ( 2001 )
( in russian ) ; doklady phys . ,
* 49 * , 231 ( 2004 ) .
v.a . dubrovskiy and v.n .
sergeev , italian journal of engineering geology and environment , special issue * 1 * , 183 ( 2006 ) .
y.pomeau and m. le berre , submitted in situations similar to one of the van der pol equation in the limit of large nonlinearities , the slow to fast transition occurs at the point on the slow manifold where the large component of the velocity is tangent to the slow manifold .
this tangency never happens in the dr equations because , there , the large component of he velocity crosses everywhere the slow manifold at finite angle .
therefore the scenario of slow to fast transition is different in the dr equations than what it is in limit cycles of the linard - van der pol like equations .
this product came from tel - atomic , inc . , a company that specializes in  tools for teaching advanced physics " .
it is no longer available as a commercial product .
one reason this instrument is no longer manufactured was the loss , due to a factory fire , of readily available chips that constituted the key component ( ne5521n integrated circuit ) used in the support electronics for its displacement sensor@xcite . a same first example of fully differential capacitive sensors is also employed ( except employed as a pair , for mechanical common mode rejection of noise ) with different electronics , in another of the company s instruments , the computerized cavendish balance @xcite .    `` symmetric differential capacitance transducer employing cross coupled conductive plates to form equipotential pairs '' , u. s. patent no .
5461319 ( 1995 ) . a description of the sensor is provided at http://telatomic.com/mechanics/sensor.html    http://telatomic.com/mechanics/cavendish-balance.html .
the support electronics for this balance employs a novel capacitive to digital converter ( ad7745 ) , which was marketed by analog devices soon after the demise of the ne5521n .
sweeping bifurcation in dynamcal systems is a well - known topic of investigations . to the best of our knowledge
it has been considered in continuous transitions , either from the point of view of bifurcation theory in the time direction inverse of ours ( namely in the direction where a pair of equilibria appears in forward time imperfect bifurcation with a slowly - varying control parameter t.erneux and p. mandel , siam journal on applied mathematics , * 46*(1986 ) , pp . 1 - 15 ) or for continuous thermodynamical phase transitions ( sweeping of an instability : an alternative to self - organized criticality to get powerlaws without parameter tuning , d. sornette , j. phys .
l * 4*(1994 ) 209 - 221 didier sornette ) .
bouchaud and cont ( j - p .
bouchaud and r. cont , eur .
j. b * 6 * , 543 ( 1998 ) ) studied also the saddle - node transition with noise , but by considering the biurcation as constant at different values near the transition value , without explicit time dependence .
bouchaud and r. cont , eur .
j. b * 6 * , 543 ( 1998 ) .