In this work we consider a family of cubic, with respect to the first derivative, nonlinear oscillators. We obtain the equivalence criterion for this family of equations and a non-canonical form of Ince Ⅶ equation, where as equivalence transformations we use generalized nonlocal transformations. As a result, we construct two integrable subfamilies of the considered family of equations. We also demonstrate that each member of these two subfamilies possesses an autonomous parametric first integral. Furthermore, we show that generalized nonlocal transformations preserve autonomous invariant curves for the equations from the studied family. As a consequence, we demonstrate that each member of these integrable subfamilies has two autonomous invariant curves, that correspond to irreducible polynomial invariant curves of the considered non-canonical form of Ince Ⅶ equation. We illustrate our results by two examples: An integrable cubic oscillator and a particular case of the Liénard (4, 9) equation.
Citation: Dmitry Sinelshchikov. On an integrability criterion for a family of cubic oscillators[J]. AIMS Mathematics, 2021, 6(11): 12902-12910. doi: 10.3934/math.2021745
[1] | P. Pirmohabbati, A. H. Refahi Sheikhani, H. Saberi Najafi, A. Abdolahzadeh Ziabari . Numerical solution of full fractional Duffing equations with Cubic-Quintic-Heptic nonlinearities. AIMS Mathematics, 2020, 5(2): 1621-1641. doi: 10.3934/math.2020110 |
[2] | SAIRA, Wenxiu Ma, Suliman Khan . An efficient numerical method for highly oscillatory logarithmic-algebraic singular integrals. AIMS Mathematics, 2025, 10(3): 4899-4914. doi: 10.3934/math.2025224 |
[3] | Kai Wang, Guicang Zhang . Curve construction based on quartic Bernstein-like basis. AIMS Mathematics, 2020, 5(5): 5344-5363. doi: 10.3934/math.2020343 |
[4] | Taher S. Hassan, Amir Abdel Menaem, Hasan Nihal Zaidi, Khalid Alenzi, Bassant M. El-Matary . Improved Kneser-type oscillation criterion for half-linear dynamic equations on time scales. AIMS Mathematics, 2024, 9(10): 29425-29438. doi: 10.3934/math.20241426 |
[5] | Dexin Meng . Wronskian-type determinant solutions of the nonlocal derivative nonlinear Schrödinger equation. AIMS Mathematics, 2025, 10(2): 2652-2667. doi: 10.3934/math.2025124 |
[6] | Samia BiBi, Md Yushalify Misro, Muhammad Abbas . Smooth path planning via cubic GHT-Bézier spiral curves based on shortest distance, bending energy and curvature variation energy. AIMS Mathematics, 2021, 6(8): 8625-8641. doi: 10.3934/math.2021501 |
[7] | Chunli Li, Wenchang Chu . Remarkable series concerning (3nn) and harmonic numbers in numerators. AIMS Mathematics, 2024, 9(7): 17234-17258. doi: 10.3934/math.2024837 |
[8] | Beatriz Campos, Alicia Cordero, Juan R. Torregrosa, Pura Vindel . Dynamical analysis of an iterative method with memory on a family of third-degree polynomials. AIMS Mathematics, 2022, 7(4): 6445-6466. doi: 10.3934/math.2022359 |
[9] | A. Palanisamy, J. Alzabut, V. Muthulakshmi, S. S. Santra, K. Nonlaopon . Oscillation results for a fractional partial differential system with damping and forcing terms. AIMS Mathematics, 2023, 8(2): 4261-4279. doi: 10.3934/math.2023212 |
[10] | Tongzhu Li, Ruiyang Lin . Classification of Möbius homogeneous curves in R4. AIMS Mathematics, 2024, 9(8): 23027-23046. doi: 10.3934/math.20241119 |
In this work we consider a family of cubic, with respect to the first derivative, nonlinear oscillators. We obtain the equivalence criterion for this family of equations and a non-canonical form of Ince Ⅶ equation, where as equivalence transformations we use generalized nonlocal transformations. As a result, we construct two integrable subfamilies of the considered family of equations. We also demonstrate that each member of these two subfamilies possesses an autonomous parametric first integral. Furthermore, we show that generalized nonlocal transformations preserve autonomous invariant curves for the equations from the studied family. As a consequence, we demonstrate that each member of these integrable subfamilies has two autonomous invariant curves, that correspond to irreducible polynomial invariant curves of the considered non-canonical form of Ince Ⅶ equation. We illustrate our results by two examples: An integrable cubic oscillator and a particular case of the Liénard (4, 9) equation.
We consider the following family of nonlinear oscillators
yzz+k(y)y3z+h(y)y2z+f(y)yz+g(y)=0, | (1.1) |
where k, h, f≠0 and g≠0 are arbitrary sufficiently smooth functions. Particular members of (1.1) are used for the description of various processes in physics, mechanics and so on and they also appear as invariant reductions of nonlinear partial differential equations [1,2,3].
Integrability of (1.1) was studied in a number of works [4,5,6,7,8,9,10,11,12,13,14,15,16]. In particular, in [15] linearization of (1.1) via the following generalized nonlocal transformations
w=F(y),dζ=(G1(y)yz+G2(y))dz. | (1.2) |
was considered. However, equivalence problems with respect to transformations (1.2) for (1.1) and its integrable nonlinear subcases have not been studied previously. Therefore, in this work we deal with the equivalence problem for (1.1) and its integrable subcase from the Painlevé-Gambier classification. Namely, we construct an equivalence criterion for (1.1) and a non-canonical form of Ince Ⅶ equation [17,18]. As a result, we obtain two new integrable subfamilies of (1.1). What is more, we demonstrate that for any equation from (1.1) that satisfy one of these equivalence criteria one can construct an autonomous first integral in the parametric form. Notice that we use Ince Ⅶ equation because it is one of the simplest integrable members of (1.1) with known general solution and known classification of invariant curves.
Moreover, we show that transformations (1.2) preserve autonomous invariant curves for equations from (1.1). Since the considered non-canonical form of Ince Ⅶ equation admits two irreducible polynomial invariant curves, we obtain that any equation from (1.1), which is equivalent to it, also admits two invariant curves. These invariant curves can be used for constructing an integrating factor for equations from (1.1) that are equivalent to Ince Ⅶ equation. If this integrating factor is Darboux one, then the corresponding equation is Liouvillian integrable [19]. This demonstrates the connection between nonlocal equivalence approach and Darboux integrability theory and its generalizations, which has been recently discussed for a less general class of nonlocal transformations in [20,21,22].
The rest of this work is organized as follows. In the next Section we present an equivalence criterion for (1.1) and a non-canonical form of the Ince Ⅶ equation. In addition, we show how to construct an autonomous first integral for an equation from (1.1) satisfying this equivalence criterion. We also demonstrate that transformations (1.2) preserve autonomous invariant curves for (1.1). In Section 3 we provide two examples of integrable equations from (1.1) and construct their parametric first integrals, invariant curves and integrating factors. In the last Section we briefly discuss and summarize our results.
We begin with the equivalence criterion between (1.1) and a non-canonical form of the Ince Ⅶ equation, that is [17,18]
wζζ+3wζ+ϵw3+2w=0. | (2.1) |
Here ϵ≠0 is an arbitrary parameter, which can be set, without loss of generality, to be equal to ±1.
The general solution of (1.1) is
w=e−(ζ−ζ0)cn{√ϵ(e−(ζ−ζ0)−C1),1√2}. | (2.2) |
Here ζ0 and C1 are arbitrary constants and cn is the Jacobian elliptic cosine. Expression (2.2) will be used below for constructing autonomous parametric first integrals for members of (1.1).
The equivalence criterion between (1.1) and (2.1) can be formulated as follows:
Theorem 2.1. Equation (1.1) is equivalent to (2.1) if and only if either
(I)25515lgp2qy+2352980l10+(3430q−6667920p3)l5−14580qp3−10q2−76545lgqppy=0, | (2.3) |
or
(II)343l5−972p3=0, | (2.4) |
holds. Here
l=9(fgy−gfy+fgh−3kg2)−2f3,p=gly−3lgy+l(f2−3gh),q=25515gylp2−5103lgppy+686l5−8505p2(f2−3gh)l+6561p3. | (2.5) |
The expression for G2 in each case is either
(I)G2=126l2qp2470596l10−(1333584p3+1372q)l5+q2, | (2.6) |
or
(II)G22=−49l3G2+9p2189pl. | (2.7) |
In all cases the functions F and G1 are given by
F2=l81ϵG32,G1=G2(f−3G2)3g. | (2.8) |
Proof. We begin with the necessary conditions. Substituting (1.2) into (2.1) we get
yzz+k(y)y3z+h(y)y2z+f(y)yz+g(y)=0, | (2.9) |
where
k=FG31(ϵF2+2)+3G21Fy+G1Fyy−FyG1,yG2Fy,h=G2Fyy+(6G1G2−G2,y)Fy+3FG2G21(ϵF2+2)G2Fy,f=3G2(Fy+FG1(ϵF2+2))Fy,g=FG22(ϵF2+2)Fy. | (2.10) |
As a consequence, we obtain that (1.1) can be transformed into (2.1) if it is of the form (2.9) (or (1.1)).
Conversely, if the functions F, G1 and G2 satisfy (2.10) for some values of k, h, f and g, then (1.1) can be mapped into (2.1) via (1.2). Thus, we see that the compatibility conditions for (2.10) as an overdertmined system of equations for F, G1 and G2 result in the necessary and sufficient conditions for (1.1) to be equivalent to (2.1) via (1.2).
To obtain the compatibility conditions, we simplify system (2.10) as follows. Using the last two equations from (2.10) we find the expression for G1 given in (2.8). Then, with the help of this relation, from (2.10) we find that
81ϵF2G32−l=0, | (2.11) |
and
567lG32+(243lgh−81lf2−81gly+243lgy)G2−7l2=0,243lgG2,y+324lG32−81glyG2+2l2=0, | (2.12) |
Here l is given by (2.5).
As a result, we need to find compatibility conditions only for (2.12). In order to find the generic case of this compatibility conditions, we differentiate the first equation twice and find the expression for G22 and condition (2.3). Differentiating the first equation from (2.12) for the third time, we obtain (2.6). Further differentiation does not lead to any new compatibility conditions. Particular case (2.4) can be treated in the similar way.
Finally, we remark that the cases of l=0, p=0 and q=0 result in the degeneration of transformations (1.2). This completes the proof.
As an immediate corollary of Theorem 2.1 we get
Corollary 2.1. If coefficients of an equation from (1.1) satisfy either (2.3) or (2.4), then an autonomous first integral of this equation can be presented in the parametric form as follows:
y=F−1(w),yz=G2wζFy−G1wζ. | (2.13) |
Here w is the general solution of (2.1) given by (2.2). Notice also that, formally, (2.13) contains two arbitrary constants, namely ζ0 and C1. However, without loss of generality, one of them can be set equal to zero.
Now we demonstrate that transformations (1.2) preserve autonomous invariant curves for equations from (1.1).
First, we need to introduce the definition of an invariant curve for (1.1). We recall that Eq (1.1) can be transformed into an equivalent dynamical system
yz=P,uz=Q,P=u,Q=−ku3−hu2−fu−g. | (2.14) |
A smooth function H(y,u) is called an invariant curve of (2.14) (or, equivalently, of (1.1)), if it is a nontrivial solution of [19]
PHy+QHu=λH, | (2.15) |
for some value of the function λ, which is called the cofactor of H.
Second, we need to introduce the equation that is equivalent to (1.1) via (1.2). Substituting (1.2) into (1.1) we get
wζζ+˜kw3ζ+˜hw2ζ+˜fwζ+˜g=0, | (2.16) |
where
˜k=kG32−gG31+(G1,y−hG1)G22+(fG1−G2,y)G1G2F2yG22,˜h=(hFy−Fyy)G22−(2fG1−G2,y)G2Fy+3gG21FyF2yG22,˜f=fG2−3gG1G22,˜g=gFyG22. | (2.17) |
An invariant curve for (2.16) can be defined in the same way as that for (1.1). Notice that, further, we will denote wζ as v.
Theorem 2.2. Suppose that either (1.1) possess an invariant curve H(y,u) with the cofactor λ(y,u) or (2.16) possess an invariant curve ˜H(w,v) with the cofactor ˜λ(w,v). Then, the other equation also has an invariant curve and the corresponding invariant curves and cofactors are connected via
H(y,u)=˜H(F,FyuG1u+G2),λ(y,u)=(G1u+G2)˜λ(F,FyuG1u+G2). | (2.18) |
Proof. Suppose that ˜H(w,v) is an invariant curve for (2.16) with the cofactor ˜λ(w,v). Then it satisfies
v˜Hw+(−˜kv3−˜hv2−˜fv−˜g)˜Hv=˜λ˜H. | (2.19) |
Substituting (1.2) into (2.19) we get
uHy+(−ku3−hu2−fu−g)H=(G1u+G2)˜λ(F,FyuG1u+G2)H. | (2.20) |
This completes the proof.
As an immediate consequence of Theorem 2.2 we have that transformations (1.2) preserve autonomous first integrals admitted by members of (1.1), since they are invariant curves with zero cofactors.
Another corollary of Theorem 2.2 is that any equation from (1.1) that is connected to (2.1) admits two invariant curves that correspond to irreducible polynomial invariant curves of (2.1). This invariant curves of (2.1) and the corresponding cofactors are the following (see, [23] formulas (3.18) and (3.19) taking into account scaling transformations)
˜H=±i√−2ϵ(v+w)+w2,˜λ=±√−2ϵw−2. | (2.21) |
Therefore, we have that the following statement holds:
Corollary 2.2. If coefficients of an equation from (1.1) satisfy either (2.3) or (2.4), then is admits the following invariant curves with the corresponding cofactors
H=±i√−2ϵ(FyuG1u+G2+F)+F2,λ=(G1u+G2)(±√−2ϵF−2). | (2.22) |
Let us remark that connections between (2.1) and non-autonomous variants of (1.1) can be considered via a non-autonomous generalization of transformations (1.2). However, one of two nonlocally related equations should be autonomous since otherwise nonlocal transformations do not map a differential equation into a differential equation [5].
In this Section we have obtained the equivalence criterion between (1.1) and (2.1), that defines two new completely integrable subfamilies of (1.1). We have also demonstrated that members of these subfamilies posses an autonomous parametric first integral and two autonomous invariant curves.
In this Section we provide two examples of integrable equations from (1.1) satisfying integrability conditions from Theorem 2.1.
Example 1. One can show that the coefficients of the following cubic oscillator
yzz−12ϵμy(ϵμ2y4+2)2y3z−6μyyz+2μ2y3(ϵμ2y4+2)=0, | (3.1) |
satisfy condition (2.3) from Theorem 2.1. Consequently, Eq (3.1) is completely integrable and its general solution can be obtained from (2.2) by inverting transformations (1.2). However, it is more convenient to use Corollary 2.1 and present the autonomous first integral of (3.1) in the parametric form as follows:
y=±√wμ,yz=w(ϵw2+2)wζ2wζ+w(ϵw2+2), | (3.2) |
where w is given by (2.2), ζ is considered as a parameter and ζ0, without loss of generality, can be set equal to zero. As a result, we see that (3.1) is integrable since it has an autonomous first integral.
Moreover, using Corollary 2.2 one can find invariant curves admitted by (3.1)
H1,2=y4[(√2±√−ϵμy2)2(√2∓√−ϵμy2)+2(ϵμy2∓√−2ϵ)u]2μ2y2(ϵμ2y4+2)−4u,λ1,2=±2(μy2(ϵμ2y4+2)−2u)(√−2ϵμy2∓2)y(ϵμ2y4+2) | (3.3) |
With the help of the standard technique of the Darboux integrability theory [19], it is easy to find the corresponding Darboux integrating factor of (3.1)
M=(ϵμ2y4+2)94(2ϵu2+(ϵμ2y4+2)2)34(μy2(ϵμ2y4+2)−2u)32. | (3.4) |
Consequently, equation is (3.1) Liouvillian integrable.
Example 2. Consider the Liénard (1, 9) equation
yzz+(biyi)yz+ajyj=0,i=0,…4,j=0,…,9. | (3.5) |
Here summation over repeated indices is assumed. One can show that this equation is equivalent to (2.1) if it is of the form
yzz−9(y+μ)(y+3μ)3yz+2y(2y+3μ)(y+3μ)7=0, | (3.6) |
where μ is an arbitrary constant.
With the help of Corollary 2.1 one can present the first integral of (3.6) in the parametric form as follows:
y=3√−2ϵμw2−√−2ϵw,yz=7776√2ϵμ5wwζ(√−2ϵw−2)5(2√−ϵwζ+(√2ϵw+2√−ϵ)w), | (3.7) |
where w is given by (2.2). Thus, one can see that (3.5) is completely integrable due to the existence of this parametric autonomous first integral.
Using Corollary 2.2 we find two invariant curves of (3.6):
H1=y2[(2y+3μ)(y+3μ)4−2u)](y+3μ)2[(y+3μ)4y−u],λ1=6μ(u−y(y+3μ)4)y(y+3μ), | (3.8) |
and
H2=y2(y+3μ)2y(y+3μ)4−u,λ2=2(2y+3μ)(u−2y(y+3μ)4)y(y+3μ). | (3.9) |
The corresponding Darboux integrating factor is
M=[y(y+3μ)4−u]−32[(2y+3μ)(y+3μ)4−2u]−34. | (3.10) |
As a consequence, we see that Eq (3.6) is Liouvillian integrable.
Therefore, we see that equations considered in Examples 1 and 2 are completely integrable from two points of view. First, they possess autonomous parametric first integrals. Second, they have Darboux integrating factors.
In this work we have considered the equivalence problem between family of Eqs (1.1) and its integrable member (2.1), with equivalence transformations given by generalized nonlocal transformations (1.2). We construct the corresponding equivalence criterion in the explicit form, which leads to two new integrable subfamilies of (1.1). We have demonstrated that one can explicitly construct a parametric autonomous first integral for each equation that is equivalent to (2.1) via (1.2). We have also shown that transformations (1.2) preserve autonomous invariant curves for (1.1). As a consequence, we have obtained that equations from the obtained integrable subfamilies posses two autonomous invariant curves, which corresponds to the irreducible polynomial invariant curves of (2.1). This fact demonstrate a connection between nonlocal equivalence approach and Darboux and Liouvillian integrability approach. We have illustrate our results by two examples of integrable equations from (1.1).
The author was partially supported by Russian Science Foundation grant 19-71-10003.
The author declares no conflict of interest in this paper.
[1] | A. A. Andronov, A. A. Vitt, S. E. Khaikin, Theory of oscillators, New York: Dover Publications, 2011. |
[2] |
S. Ghosh, D. S. Ray, Chemical oscillator as a generalized Rayleigh oscillator, J. Chem. Phys., 139 (2013), 164112. doi: 10.1063/1.4826169
![]() |
[3] | A. L. Kazakov, P. A. Kuznetsov, A. A. Lempert, On a heat wave for the nonlinear heat equation: An existence theorem and exact solution, In: G. Demidenko, E. Romenski, E. Toro, M. Dumbser, Continuum mechanics, applied mathematics and scientific computing: Godunov's legacy, Cham: Springer, 2020,223–228. |
[4] | C. Muriel, J. L. Romero, Second-order ordinary differential equations and first integrals of the form A(t,x)˙x+B(t,x), J. Nonlinear Math. Phys., 16 (2009), 209–222. |
[5] | W. Nakpim, S. V. Meleshko, Linearization of second-order ordinary differential equations by generalized sundman transformations, Symmetry, Integr. Geom. Methods Appl., 6 (2010), 051. |
[6] | M. C. Nucci, K. M. Tamizhmani, Lagrangians for dissipative nonlinear oscillators: The method of Jacobi last multiplier, J. Nonlinear Math. Phys., 17 (2010), 167–178. |
[7] | C. Muriel, J. L. Romero, Second-order ordinary differential equations with first integrals of the form C(t)+1/(A(t,x)˙x+B(t,x)), J. Nonlinear Math. Phys., 18 (2011), 237–250. |
[8] |
Y. Y. Bagderina, Invariants of a family of scalar second-order ordinary differential equations for Lie symmetries and first integrals, J. Phys. A Math. Theor., 49 (2016), 155202. doi: 10.1088/1751-8113/49/15/155202
![]() |
[9] | A. Ruiz, C. Muriel, On the integrability of Liénard I-type equations via λ-symmetries and solvable structures, Appl. Math. Comput., 339 (2018), 888–898. |
[10] | D. I. Sinelshchikov, N. A. Kudryashov, On the Jacobi last multipliers and Lagrangians for a family of Liénard-type equations, Appl. Math. Comput., 307 (2017), 257–264. |
[11] |
M. V. Demina, Novel algebraic aspects of Liouvillian integrability for two-dimensional polynomial dynamical systems, Phys. Lett. A, 382 (2018), 1353–1360. doi: 10.1016/j.physleta.2018.03.037
![]() |
[12] | J. Giné, C. Valls, Liouvillian integrability of a general Rayleigh-Duffing oscillator, J. Nonlinear Math. Phys., 26 (2019), 169–187. |
[13] |
J. Giné, C. Valls, On the dynamics of the Rayleigh-Duffing oscillator, Nonlinear Anal. Real World Appl., 45 (2019), 309–319. doi: 10.1016/j.nonrwa.2018.07.007
![]() |
[14] |
M. R. Cândido, J. Llibre, C. Valls, Non-existence, existence, and uniqueness of limit cycles for a generalization of the Van der Pol-Duffing and the Rayleigh-Duffing oscillators, Phys. D Nonlinear Phenom., 407 (2020), 132458. doi: 10.1016/j.physd.2020.132458
![]() |
[15] |
D. I. Sinelshchikov, Linearizability conditions for the Rayleigh-like oscillators, Phys. Lett. A, 384 (2020), 126655. doi: 10.1016/j.physleta.2020.126655
![]() |
[16] |
D. I. Sinelshchikov, On linearizability via nonlocal transformations and first integrals for second-order ordinary differential equations, Chaos, Solitons Fractals, 141 (2020), 110318. doi: 10.1016/j.chaos.2020.110318
![]() |
[17] | E. L. Ince, Ordinary differential equations, Dover, New York, 1956. |
[18] |
F. J. Bureau, Differential equations with fixed critical points, Ann. Mat. Pur. Appl., 64 (1964), 229–364. doi: 10.1007/BF02410054
![]() |
[19] | X. Zhang, Integrability of dynamical systems: Algebra and analysis, Vol. 47, Singapore: Springer, 2017. |
[20] |
M. Demina, D. Sinelshchikov, Integrability properties of cubic Liénard oscillators with linear damping, Symmetry, 11 (2019), 1378. doi: 10.3390/sym11111378
![]() |
[21] |
M. V. Demina, D. I. Sinelshchikov, On the integrability of some forced nonlinear oscillators, Int. J. Nonlin. Mech., 121 (2020), 103439. doi: 10.1016/j.ijnonlinmec.2020.103439
![]() |
[22] |
M. V. Demina, D. I. Sinelshchikov, Darboux first integrals and linearizability of quadratic-quintic Duffing-van der Pol oscillators, J. Geom. Phys., 165 (2021), 104215. doi: 10.1016/j.geomphys.2021.104215
![]() |
[23] |
M. V. Demina, Liouvillian integrability of the generalized Duffing oscillators, Anal. Math. Phys., 11 (2021), 1–18. doi: 10.1007/s13324-020-00437-5
![]() |
1. | Dmitry I. Sinelshchikov, Linearizabiliy and Lax representations for cubic autonomous and non-autonomous nonlinear oscillators, 2023, 01672789, 133721, 10.1016/j.physd.2023.133721 | |
2. | Jaume Giné, Xavier Santallusia, Integrability via algebraic changes of variables, 2024, 184, 09600779, 115026, 10.1016/j.chaos.2024.115026 | |
3. | Meryem Belattar, Rachid Cheurfa, Ahmed Bendjeddou, Paulo Santana, A class of nonlinear oscillators with non-autonomous first integrals and algebraic limit cycles, 2023, 14173875, 1, 10.14232/ejqtde.2023.1.50 | |
4. | Jaume Giné, Dmitry Sinelshchikov, Integrability of Oscillators and Transcendental Invariant Curves, 2025, 24, 1575-5460, 10.1007/s12346-024-01182-x |