Capture into parametric autoresonance in the presence of noise
Abstract. System of differential equations describing the initial stage of the capture of oscillatory systems into the parametric autoresonance is considered. Of special interest are solutions whose amplitude increases without bound with time. The possibility of capture the system into the autoresonance is related with the stability of such solutions. We study the stability of autoresonant solutions with respect to persistent perturbations of white noise type, and we show that under certain conditions on the intensity of the noise, the capture into parametric autoresonance is preserved with probability tending to one.
Keywords: nonlinear oscillations, autoresonance, perturbation, white noise, stability
Mathematics Subject Classification: 34C15, 70K20, 37B25
Autoresonance is the phenomenon of continuous phase locking of nonlinear oscillator with slowly varying parametric pumping that leads to significant growth of the energy of the oscillator. This phenomenon was first suggested in the problem of acceleration of relativistic particles [1, 2]. Later, it was observed that autoresonance occurs widely in nature and plays the important role in many problems of nonlinear physics [3, 4, 5]. A wide range of applications requires the study of the effect of perturbations on different mathematical models of autoresonance. The influence of additive noise on the capture into the autoresonance (non-parametric) was analysed in , where the effect of perturbations was considered only at the initial time. The problem of capture into the parametric autoresonance for a quantum anharmonic oscillator with initial disturbances was discussed in . The effect of persistent perturbations with random jumps on the stability of autoresonance models was investigated in . In this paper we consider the deterministic model of parametric autoresonance  and we study the effect of persistent perturbations of white noise type on the captured solutions.
The paper is organized as follows. In section 1, we give the mathematical formulation of the problem. In section 2 we discuss the stability of autoresonant solutions with respect to perturbations of initial data. Section 3 deals with a more general problem of stochastic stability of a class of locally stable dynamical systems. In section 4 these results are applied in the study of the capture into parametric autoresonance in the presence of stochastic perturbations.
1. Problem statement
We consider the system of primary parametric autoresonance equations
where and are parameters. This system appears after the averaging of equations, describing the behaviour of nonlinear oscillators in the presence of small slowly changing parametric pumping. The real-valued functions and represent the amplitude and phase shift of harmonic oscillations. Solutions with an infinitely growing amplitude and bounded phase shift as correspond to the capture into the parametric autoresonance. The existence and asymptotic behaviour of captured solutions of system (1) were discussed in [9, 10, 11, 12]. In this paper we study the effect of white noise on the stability of such solutions.
The asymptotic solution of system (1) with growing amplitude at infinity can be constructed in the following form
where and are constant coefficients. Substituting these series in system (1) and grouping the expressions of same power of give the following recurrence relations for determining the coefficients , : , , , , etc. Note that there are no free parameters in the asymptotic series. The existence of exact particular solutions of system (1) with constructed asymptotics follows from . The stability of these isolated solutions determines the presence of capture into autoresonance. We show that in the case of stability, solutions with asymptotics (2) attract many other autoresonant solutions with more complicated asymptotic expansions (see Fig. 1). Note that system (1) has also non-resonant solutions with the slipping phase and the bounded amplitude. The existence of such solutions excludes the global stability of autoresonant solutions for all initial data. We also note that the structure of the capture region (the set of initial points such that the corresponding solutions possess the unboundedly growing amplitude) for system (1) remains unknown, and we do not discuss this problem here.
Consider the perturbed system in the form
where the stochastic processes and defined on a probability space play the role of perturbations. It is assumed that , for all and , where is the Kronecker delta, is the Dirac delta function, and the deterministic functions together with the small parameter are used to control the intensity of the perturbations. Let , where and are independent Wiener processes. Then we can consider the perturbed system (3) in the form of Itô stochastic differential equations. Our goal is to find constraints on the functions , , such that the capture into parametric autoresonance is preserved in the perturbed system with probability tending to one. Since the persistent perturbation of white noise type leads to the loss of stability of solutions for all (see [14, 15, 16] and Fig. 2), we consider a weaker problem. Specifically, our goal is to find the largest possible time interval on which the stability of autoresonant solutions is preserved.
Note that system (1) is of universal character in the description of parametric autoresonance in nonlinear systems. It describes long-term evolution of different nonlinear oscillations under small parametric driving. As but one example let us consider the following equation
where , , . The functions , play the role of perturbations. Solutions of equation (4) with whose amplitude increase with time from small values to quantities of order one are associated with the capture into parametric autoresonance. For the asymptotic description of such solutions at the initial stage of the capture we use the method of two scales. We introduce a slow time and a fast variable . Then the asymptotic substitution
2. Perturbations of initial data for autoresonant solutions
Note that the unperturbed system (1) has two different asymptotic solutions in the form (2) distinguished by the choice of a root to the equation . It can easily be checked that the solution with is linearly unstable. However, linear stability analysis fails for the captured solution , with asymptotics (2), . To study the stability of this solution, we need to take into account high-order terms of the equations. In our analysis we use only the first terms of the asymptotic expansion for the solution,
In system (1) we make the change of variables , , and for new functions , we study the stability of the trivial solution , to the following system close to Hamiltonian system
By taking into account the asymptotics of the captured solution , one can readily write out the asymptotics of the function as and :
The asymptotic estimates are uniform with respect to in the domain with positive constants and . It is clear that the function is positive definite function in the neighbourhood of the equilibrium . A Lyapunov function candidate for system (7) is constructed of the form
From the properties of the function it follows that there exist and such that
for all . Integrating the last expression with respect to , we obtain the following estimates
From the inequality for the total derivative of the Lyapunov function along the trajectories of system (7) it follows that
as , where . The change-of-variables formula implies the following asymptotic estimates for solutions to system (1) with initial data from the -neighbourhood of the isolated autoresonant solution , as . ∎
3. Stochastic perturbations of locally stable systems
In the study of stochastic perturbations of system (7) and other similar equations with locally stable solutions it is convenient to consider the system of differential equations
where , and . Suppose that the vector-valued function is continuous and for all satisfies a Lipschitz condition: for all , with positive constant . Assume that there exists a local Lyapunov function for system (9) satisfying the inequalities:
in the domain with parameters , . The existence of such Lyapunov function guaranties that the trivial solution is locally stable with respect to variables (if , the trivial solution is stable with respect to all variables). Note that if , then the solution is stable with respect to variables in some weighted norm. Let us remark that the Lyapunov function, constructed in the previous section for the system of primary parametric autoresonance, possesses the similar estimates (cp. (8) with (10)). Note also that such Lyapunov functions are constructed in the stability analysis of nonlinear non-autonomous systems of differential equations (see, for example, [17, 18, 19]).
Together with system (9) we consider the perturbed system in the form of Itô stochastic differential equations
where is -dimensional Wiener process defined on a probability space , is a continuous matrix which is independent of and for all satisfies the following conditions , for all , with positive constants . We assume that does not depend on . These constraints on the coefficients of system (11) guarantee the existence and uniqueness of solution for all and for all (see, for instance, [20, §5.2], [21, §3.3]). We assume that the perturbed system does not preserve the trivial solution, . Define the class of perturbations as a set of matrices such that for all and , where .
We study the stability of the solution of system (9) with respect to stochastic perturbations on a finite time interval. One variant of this approach is to find the largest possible time interval on which solutions to the perturbed system (11) are close to the equilibrium of the deterministic system (9) (see, for instance, [14, Chap. 9] and [22, Chap. 7]). We have
Let us fix the parameters , , and . Let be a solution of system (11) with and initial data , , and let be the first exit time of the solution from the domain
We define the function , then is the process stopped at first exit time from the domain . Positive parameters , will be specified later. Let us first consider the case . The Lyapunov function for the stochastic system (11) is constructed in following the form
In the study of stability of solutions to stochastic differential equations the following operator plays the role of the total derivative along the trajectories [21, §3.6]: It easy to see that and
for all . These estimates guarantee that is a nonnegative supermartingale [21, §5.2]. Using the properties of the function and Doob’s inequality for supermartingales, we get the following estimates
Define and the parameters and . Then for all , . If , this estimate holds for all . Taking into account (13), we obtain the estimate
The stability for is proved by using the Lyapunov function in the following form 
where . It is easy to check that the following inequalities hold
for all and any natural number . Since for all , we see that the function is a nonnegative supermartingale and the following estimates hold
Now we define ; then as and with a positive constant . Therefore, for all there exist such that for all and (if , we can choose and ). Taking into account (15), we obtain (14). Thus, for any natural and for all the trivial solution to system (9) is stable with respect to variables under stochastic perturbations on the interval uniformly for . Note that if , then from similar arguments it follows that
Note that if with and , then it can be proved that the stochastic stability of the trivial solution holds on the time interval . In this case, the Lyapunov function has the form with a positive constant .
4. Stochastic perturbations of stable autoresonant solution
In this section we study the stability of the autoresonant solution , with asymptotics (5) under stochastic perturbations. In system (3) we make the change of variables , ; then for the functions , we have the following system of stochastic differential equations
Note that system (7) has the Lyapunov function satisfying (10) with , and . The restrictions of the coefficients imply that . If we combine this with Theorem 2, we get the stochastic stability of the trivial solution to system (7) and the estimates (12) with and . By means of change of variables we derive the inequalities (17). ∎
Thus, the stability of the isolated autoresonant solution , is preserved in the perturbed system on asymptotically long time interval , . Therefore, the perturbation of white noise type with moderate intensity cannot destroy the stability of the capture into parametric autoresonance.
-  V. I. Veksler, A new method of acceleration of relativistic particles, J. Phys. USSR, 9, 153–158 (1945)
-  E. M. McMillan, The synchrotron - a proposed high energy particle accelerator, Phys. Rev., 68, 143–144 (1945)
-  J. Fajans, E. Gilson, and L. Friedland, Second harmonic autoresonant control of the diocotron mode in pure-electron plasmas, Phys. Rev. E, 62, 4131 (2000)
-  J. Fajans and L. Friedland, Autoresonant (nonstationary) excitation of pendulums, Plutinos, plasmas, and other nonlinear oscillators, Am. J. Phys., 69, 1096–1102 (2001)
-  L. Friedland, Scholarpedia 4, 5473 (2009).
-  I. Barth, L. Friedland, E. Sarid, and A. G. Shagalov, Autoresonant transition in the presence of noise and self-fields, Phys. Rev. Lett., 103, 155001 (2009)
-  I. Barth and L. Friedland, Quantum phenomena in a chirped parametric anharmonic oscillator, Phys. Rev. Lett., 113, 040403 (2014)
-  O. A. Sultanov, Stability of autoresonance models subject to random perturbations for systems of nonlinear oscillation equations, Comput. Math. and Math. Phys., 54, 59–73 (2014)
-  E. Khain and B. Meerson, Parametric autoresonance, Phys. Rev. E, 64, 036619 (2001)
-  M. Assaf and B. Meerson, Parametric autoresonance in Faraday waves, Phys. Rev. E, 72, 016310 (2005)
-  O.M. Kiselev and S.G. Glebov, The capture into parametric autoresonance, Nonlinear Dynam., 48, 217–230 (2007)
-  L. A. Kalyakin, Asymptotic analysis of autoresonance models, Russian Math. Surveys, 63, 791–857 (2008)
-  A. N. Kuznetsov, Existence of solutions entering at a singular point of an autonomous system having a formal solution, Funct. Anal. Appl., 23, 308–317 (1989)
-  M. I. Freidlin and A. D. Wentzell, Random Perturbations of Dynamical Systems. New York, Heidelberg, Berlin : Springer-Verlag, 1998.
-  R. Khasminskii, Stochastic Stability of Differential Equations. Berlin, Heidelberg: Springer-Verlag, 2012.
-  A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization. A Universal Concept in Nonlinear Sciences. Cambridge: Cambridge University Press, 2001.
-  V. I. Vorotnikov, Partial Stability and Control. Boston, Basel, Berlin : Birkhäuser, 1998.
-  H. K. Khalil, Nonlinear systems. Prentice Hall, Upper Saddle River, NJ (2002)
-  L. A. Kalyakin, Synchronization in a nonisochronous nonautonomous system, Theoret. and Math. Phys., 181:2, 1339–1348 (2014)
-  B. Øksendal, Stochastic Differential Equations. An Introduction with Applications. New York, Heidelberg, Berlin : Springer-Verlag, 1998.
-  R. Khasminskii, Stochastic Stability of Differential Equations. Berlin, Heidelberg : Springer-Verlag, 2012.
-  M. M. Hapaev, Averaging in Stability Theory: A Study of Resonance Multi-Frequency Systems. Dordrecht, Boston : Kluwer Academic Publishers, 1993.
-  O. Sultanov, White noise perturbation of locally stable dynamical systems, Stochastics and Dynamics, 17:1, 1750002 (2017)