SciELO - Scientific Electronic Library Online

vol.64 número6Cosmological analysis of a Dvali-Gabadadze-Porrati stable model with H (z) observationsElectronic excitation of atoms by positron impact using the scaling born approach índice de autoresíndice de assuntospesquisa de artigos
Home Pagelista alfabética de periódicos  

Serviços Personalizados




Links relacionados

  • Não possue artigos similaresSimilares em SciELO


Revista mexicana de física

versão impressa ISSN 0035-001X

Rev. mex. fis. vol.64 no.6 México Nov./Dez. 2018  Epub 05-Nov-2019 


Mathematical modeling of DNA vibrational dynamics and its solitary wave solutions

Reza Abazaria  * 

Shabnam Jamshidzadehb 

Gangwei Wangc  d 

a Departamento de Matemática, Facultad de Ciencias, Universidad de Mohaghegh Ardabili, P. O. Box: 56199-179, Ardabil, Irán,

b Young Researchers and Elite Club, Ardabil Branch, Islamic Azad University, Ardabil, Iran.

c School of Mathematics and Statistics, Hebei University of Economics and Business, Shijiazhuang 050061, PR China.

d School of Mathematics and Statistics, Beijing Institute of Technology, Beijing 100081, PR China.


In this work, the traveling wave solutions of a mathematical modeling of DNA vibration dynamics proposed by Peyrard-Bishop, that takes into consideration the inclusion of nonlinear interaction between adjacent displacements along the Hydrogen bonds, is investigated by both (G'/G)-expansion and F-expansion methods. Using these methods, some new explicit forms of traveling wave solutions of present nonlinear equation are given. The methods come in to be easier and faster by means of a symbolic computation and yield powerful mathematical tools for solving nonlinear evolution equations in many branches of sciences, especially Physics, Biology, etc.

Keywords: (G'=G)-expansion method; F-expansion method; generalized DNA vibration dynamics; hyperbolic function solutions; trigonometric function solutions; solitary wave solutions

PACS: 02.30.Jr; 03.65.Ge; 04.20.Jb; 87.14.Gg

1. Introduction

In recent decades, a new attitude with regard to the exploration of nonlinear evolution equations (NLEEs) has been actively progressing in various branches of Sciences. Nonlinear evolution equations have been the important subject of investigation in various branches of mathematical and physical sciences such as physics, fluid mechanics, chemistry, biology and etc. Obtaining the analytical solutions of NLEEs is an important topic in various areas of Science, since many of mathematical and physical models are explained by NLEEs. Among the possible solutions to NLEEs, specific form of solutions may contingent only on a single combination of variables such as solitons. In mathematics and physics, a soliton is a self-reinforcing solitary wave that sustains its form through it travels with constant speed. Solitons are the special solutions of certain nonlinear partial differential equations, with interesting properties. Because of a balance between nonlinear and linear effects, the shape of soliton wave pulses does not change during propagation in a medium. The soliton phenomenon firstly described by John Scott Russell (1808-1882). The main idea of this phenomenon occurred to him when he observed a solitary wave in the Union Canal in Scotland. Russell reproduced same phenomenon in a wave tank and named it the “Wave of Translation" (also known as solitary wave or soliton) [1].

The mathematical modeling of DNA vibration dynamics is proposed by Peyrard-Bishop, that takes into consideration the inclusion of nonlinear interaction between adjacent displacements along the Hydrogen bonds [8]. To study the main important properties of DNA structure, we have investigated the oscillator-chain of Peyrard-Bishop model (PB) [8] which has successfully predicated the appearance of solitonic structures. As it is well known, the balance between weak nonlinearity and dispersion in DNA dynamic model yields the typical derivation of model equations mathematically. This dispersion in many cases enters only in the linear level. Unfortunately, this approach eliminates other possibilities related to the presence of nonlinear dispersion [9,10]. Some mathematical structure of DNA has been described in various interesting lines of research, to see more structure and physical properties of DNA, the authors can Ref. to [2-6] and their references. The mathematical and physical modelling of DNA dynamics has reduced to an important nonlinear structures. The nonlinearity of the DNA dynamic model causes it to form localized waves. The localized waves have some interesting properties, such as the ability of transporting energy without dissipation [2-7].

Taking into consideration the inharmonic potential, Aguero et al. [6] and Najera et al. [11] have studied the following modified PB model:

utt-l21+215l2ux2uxx-2αDe-αu(e-αu-21)=0, (1)

where l21,l2,D and α are constants. Krumhansl et al., [12], suggested a theory of soliton excitations as an explanation of the open states of DNA modeling system. In [12], they firstly developed the possibility that nonlinear effects might concentrate vibrational energy in DNA into localized soliton like objects. Some other similar worksare studied about solitons of DNA dynamics model, such as Yomosa in [13], proposed a soliton theory using a plane base-rotor model. This model was further refined by Takeno and Homma [14], who allowed discreteness effects to be taken into account, and by Zhang [15], who improved the model for base coupling. Zayed and Arnous also applied Homogeneous Balance Method in [35] and generalized Riccati equation mapping method in [36,37] to find traveling wave solutions of PB model (1.1).

In the recent decade, many powerful and direct methods have been proposed to find special solutions of nonlinear evolution equations (NLEEs), such as the Bácklund transformation [16] and Hirota bilinear method [17]. With the help of the computer implementations, some other algebraic method proposed, such as tanh/coth method [18], homogeneous balance method [19], the Miura transformation [20], sine/cosine method [21] and Exp-function method [22] and some other, see [23]. But most of the methods may sometimes fail or can only lead to a kind of special solution and the solution procedures become very complex as the degree of nonlinearity increases.

Since 2008, the (G'/G)-expansion method, firstly introduced by Wang et al. [24], has become widely used to search for various exact solutions of NLEEs [24-31]. The worth of the (G'/G)-expansion method is that it reduced the nonlinear PDEs to ODEs by some essentially linear methods, i.e. the method is based on the explicit linearization of NLEEs for traveling waves which leads to a second-order differential equation with constant coefficients. The F-expansion method is also an effective and direct algebraic method for constructing the exact solutions of nonlinear evolution equations [32]. To do this, it reduces the nonlinear evolution equations to a simple algebraic equation. Many of NLEEs have been solved by F-expansion method, (for a deeper discussion we refer the reader to [33,34] and the references given there).

In this paper, we first describe briefly the mathematical concept of both (G'/G)-expansion method and F-expansion method, then we use them to investigate the traveling wave solutions of a mathematical modeling of PB model (1.1).

2. Description of the (G'/G)-expansion method

In this section, we introduce briefly the (G'/G)-expansion method for solving certain nonlinear partial differential equations (PDEs). For a deeper discussion of (G'/G)-expansion method we refer the reader to [24-27]. Suppose we have a nonlinear PDE for u(x,t), in the form

P(u,ut,ux,uxt,uxx,utt,)=0, (2)

where P is a polynomial in its arguments, which includes nonlinear terms and the highest order derivatives. Using the transformation u(x,t)=U(ξ),ξ=x-ωt, Eq. (2) reduces to the ordinary differential equation (ODE)

P(U,-ωU',U',-ωU,U,ω2U,)=0, (3)

where U=U(ξ), and prime denotes derivative with respect to ξ. We assume that the solution of Eq. (3) can be expressed by a polynomial in (G'/G) as follows:

U(ξ)=n=21mαnG'Gn+α0,αm0 (4)

where αn,n=0,21,2,,m, are constants to be determined later and G(ξ) satisfies a second order linear ordinary differential equation (LODE):

d2G(ξ)dξ2+λdG(ξ)dξ+μG(ξ)=0 (5)

where λ and μ are arbitrary constants. Using the general solutions of Eq. (5), we have

G'(ξ)G(ξ)=λ2-21μ2C21sinhλ2-21μ2ξ+C2coshλ2-21μ2ξC21coshλ2-21μ2ξ+C2sinhλ2-21μ2ξ-λ2,λ2-21μ>0,21μ-λ22-C21sin21μ-λ22ξ+C2cos21μ-λ22ξC21cos21μ-λ22ξ+C2sin21μ-λ22ξ-λ2,λ2-21μ<0, (6)

and it follows, from (4) and (5), that


U=n=1mnαn((n+1)(G'G)n+2+(2n+1)λ(G'G)n+1+n(λ2+2μ)(G'G)n+(2n-1)λμ(G'G)n-1+(n-1)μ2(G'G)n-2), (7)

and so on, here the prime denotes the derivative with respective to ξ. To determine u explicitly, we take the following four steps:

  • Step 1. Determine the integer m by substituting Eq. (4) along with Eq. (5) into Eq. (3), and balancing the highest order nonlinear term(s) and the highest order partial derivative.

  • Step 2. Substitute Eq. (4) give the value of m determined in Step 1, along with Eq. (5) into Eq. (3) and collect all terms with the same order of (G'/G) together, the left-hand side of Eq. (3) is converted into a polynomial in (G'/G). Then set each coefficient of this polynomial to zero to derive a set of algebraic equations for ω,λ,μ,αn for n=0,1,2,,…,,m.

  • Step 3. Solve the system of algebraic equations obtained in Step 2, for ω,λ,μ,α0,…,,αm by use of computer programs such Matlab, Maple and Mathematica.

  • Step 4. Use the results obtained in above steps to derive a series of fundamental solutions u(ξ) of Eq. (3) depending on (G'/G), since the solutions of Eq. (5) have been well known for us, then we can obtain exact solutions of Eq. (2).

2.1. Application

Let us consider the Peyrard-Bishop DNA dynamic model equation

utt-(l21+215l2ux2)uxx-2αDe-αu(e-αu-21)=0, (8)

where l21,l2,α and D are constants.

We would like to use our method to obtain new and more general, travelling wave solutions of Eq. (8) by assuming the solution in the following frame:

u=U(ξ),ξ=x-ωt, (9)

where ω is arbitrary constants generally termed the wave velocity, and prime denotes derivative with respect to ξ. Using the wave variable ξ in (8) we find

ω2U-(l1+3l2(U')2)U-2αDe-αU(e-αU-1)=0, (10)

Multiplying Eq. (10) by U' and integrating once with respect to ξ, we obtain

12(ω2-l1)(U')2-34l2(U')4+De-αU(e-αU-2)+R=0, (11)

whereℛ is an integration constant. Introducing the new variable

ϕ(ξ)=e-αU(ξ), (12)

reduces Eq. (11) into the following form:

12α2(ω2-l1)ϕ2(ϕ')2-34α4l2(ϕ')4+Dϕ5(ϕ-2)+Rϕ4=0. (13)

According to Step 1, considering the homogeneous balance between (ϕ')21 and ϕ6 we get 4m+4=6m, hence m=2 We then suppose that Eq. (13) has the following formal solutions:

ϕ=α2(G'/G)2+α21(G'/G)+α0, (15)

where α20,α1, and α0, are constants which are unknown, to be determined later.

Substituting Eq. (14) along with Eq. (5) into Eq. (13) and collecting all terms with the same order of (G'/G), together, the left-hand sides of Eq. (13) are converted into a polynomial in (G'/G). Equating each coefficient of this polynomial to zero yields a set of simultaneous algebraic equations for λ,μ,ω,α0,α1 and α2. Solving the system of algebraic equations with the aid of Maple 12, we obtain the following general results

μ=±3α26l2l2D±D+D2-DR,λ=0,ω=±23l2DD2-DR+l1,α0=±D+D2-DRD,α1=0,α2=±23α2l2D (15)

Therefore, substitute the obtained results (15) in (14), we get

ϕ(ξ)=±23α2l2DG'G±2±D+D2-DRD (16)

Substituting the general solutions (6) into Eq. (14), we obtain following two types of traveling wave solutions of Eq. (8):

2.1.1. Hyperbolic type function solutions

When μ<0, using the relationship u=(-21/α)lnϕ, we obtain hyperbolic type function solution uH, of Peyrard-Bishop DNA dynamic model (8) as follows:

u±H(x,t)=-1αln(±23α2l2D(G'G)±2±D+D2-DRD), (17)





are arbitrary constants and ℛ is integration constant. It is easy to see that the hyperbolic type solution (15) can be rewritten at C212>C22, as follows

u±H1(x,t)=-1αln(D+D2-DRDcosh2(-66-3α2l2l2D(D+D2-DR)ξ+ηH)), (18a)

while at C212<C22, one can obtain

u±H2(x,t)=-1αln(-D-D2-DRDcosh2(-66-3α2l2l2D(D+D2-DR)ξ+ηH)-1). (18b)


ξ=x23l2DD2-DR+l1) t,μ=±3α26l2l2D(±D+D2-DR),ηH=tanh-1(C1C2),

are arbitrary constants and ℛ is integration constant. By setting appropriate parameters, the solution ([18a]) can be reduce to solution (31) in .

2.1.2. Trigonometric type function solutions

Now when μ>0, using the relationship (-1/α)lnϕ, the trigonometric type function solution uT, of Peyrard-Bishop DNA dynamic model (8) can be obtain as follows:

u±T(x,t)=-1αln(±23α2l2D(G'G)±2±D+D2-DRD), (19)





are arbitrary constants and ℛ is integration constant. It is easy to see that the trigonometric type solution (17) can be rewritten at C12>C22, as follows

u±T1(x,t)=-1αln(D+D2-DRD(tan2(663α2l2l2D(D+D2-DR)ξ+ηT)+1)), (20a)

and for C12<C22, one can obtain

u±T2(x,t)=-1αln(D+D2-DRD(cot2(663α2l2l2D(D+D2-DR)ξ+ηT)+1)). (20b)


ξ=x(23l2DD2-DR+l1) t,μ=±3α26l2l2D(±D+D2-DR),ηT=tan-1(C1C2),

are arbitrary constants and ℛ is integration constant.

3. F-expansion method

In order to get more solutions of (1), we use the F-expansion method [38] to deal with (13). Supposing that Eq. (13) has the following formal solutions:

ϕ=α2(φ)2+α217(φ)+α0, (21)

where α20,α1, and α0, are constants to be determined further. And 𝜑 satisfy

φ'=C(φ)2+B(φ)+A (22)

Substituting (21) and (22) into (13) and collecting all terms with the same order of φ together equating each coefficient of this polynomial to zero, one can get

Case 1:

R=-4AC(-α2(D)l2D3+3ACl2)α-4,α=α,ω=ω,α0=2l2DAC3α-2,l1=(12ACl2-2α2(D)l2D3+α2ω)α-2,l2=l2,A=A,C=C,D=D,B=0,α1=0,α2=23C2α2l2D. (23)

Case 2:

R=-4AC(α2(D)l2D3+3ACl2)α-4,α=α,ω=ω,α0=-2l2DAC3α-2,l1=(12ACl2+2α2(D)l2D3+α2ω)α-2,l2=l2,A=A,C=C,D=D,B=0,α1=0,α2=-23C2α2l2D. (24)

Thus, substitute the obtained results (23) and (24) in (22), we get

ϕ(ξ)=±23C2α2l2Dφ2±2l2DAC3α-2 (25)

Substituting the general solutions (22) into Eq. (21), we obtain the following two types of traveling wave solutions of Eq. (8):

In view of (22) has a lot of fundamental solutions (twenty seven solutions) [33], one can find a number of exact travelling wave solutions for (1), which are listed some special solutions as follows.

Family 1: -4AC>0 and AC0,

u=-1αln(23C2α2l2D((-12C[-4ACtanh(-4AC2ξ)]))2+2l2DAC3α-2). (26)

u=-1αln(23C2α2l2D((-12C[-4ACcoth(-4AC2ξ)]))2+2l2DAC3α-2). (27)

u=-1αln(23C2α2l2D(-12C[-4AC(tanh(-4ACξ)±isech(-4ACξ))])2+2l2DAC3α-2). (28)

u=-1αln(23C2α2l2D(-12C-4AC(coth(-4ACξ)±icsch(-4ACξ)))2+2l2DAC3α-2). (29)

u=-1αln(23C2α2l2D(-14C-4AC(tanh(-4AC4ξ)+coth(-4AC4ξ)))2+2l2DAC3α-2). (30)

u=-1αln(23C2α2l2D(12C[(E2+F2)(-4AC)-E-4ACcosh(-4ACξ)Esinh(-4ACξ)+F])2+2l2DAC3α-2). (31)

u=-1αln(23C2α2l2D(12C[-(F2-E2)(-4AC)+E-4ACsinh(-4ACξ)Ecosh(-4ACξ)+F])2+2l2DAC3α-2). (32)

where E and F are two non-zero real constants and satisfies F2-E2>0.

u=-1αln(23C2α2l2D([2Acosh(-4ACξ)-4ACsinh(-4ACξ)±i-4ACξ])2+2l2DAC3α-2). (33)

u=-1αln(23C2α2l2D([2Asinh(-4ACξ)-4ACcosh(-4ACξ)±-4ACξ])2+2l2DAC3α-2). (34)

u=-1αln(23C2α2l2D([4Asinh(-4AC4ξ)cosh(-4AC4ξ)2-4ACcosh2(-4AC4ξ)--4AC])2+2l2DAC3α-2). (35)

Family 2 : When -4AC<0 and AC0,

u=-1αln(23C2α2l2D(12C[4ACtan(4AC2ξ)])2+2l2DAC3α-2). (36)

u=-1αln(23C2α2l2D(12C[4ACcot(4AC2ξ)])2+2l2DAC3α-2). (37)

u=-1αln(23C2α2l2D(12C[4AC(tan(4ACξ)±sec(4ACξ))])2+2l2DAC3α-2). (38)

u=-1αln(23C2α2l2D(-12C[4AC(cot(4ACξ)±csc(4ACξ))])2+2l2DAC3α-2). (39)

u=-1αln(23C2α2l2D(14C[4AC(tan(4AC4ξ)-cot(4AC4ξ))])2+2l2DAC3α-2). (40)

u=-1αln(23C2α2l2D(12C[±(F2-E2)(4AC)-E4ACcos(4ACξ)Esin(4ACξ)+F])2+2l2DAC3α-2). (41)

u=-1αln(23C2α2l2D(12C[±(F2-E2)(4AC)+E4ACsin(4ACξ)Ecos(4ACξ)+F])2+2l2DAC3α-2). (42)

where E and F are two non-zero real constants and satisfies F2-E2>0.

u=-1αln(23C2α2l2D([-2Acos(4ACξ)4ACsin(4ACξ)±i4ACξ])2+2l2DAC3α-2). (43)

u=-1αln(23C2α2l2D([2Asin(4ACξ)4ACcos(4ACξ)±4ACξ])2+2l2DAC3α-2). (44)

u=-1αln(23C2α2l2D([4Asin(4AC4ξ)cos(4AC4ξ)24ACcos2(4AC4ξ)-4AC])2+2l2DAC3α-2). (45)

Remark 1. Using Case. 2, one can get other exact solutions of Eq. (1). The details are omitted here.

4. Conclusions

This study shows that both (G'/G)-expansion method and F-expansion method are quite efficient and practically well suited for use in finding exact solutions for a mathematical modeling of DNA vibration dynamics. The reliability of the methods and the reduction in the size of computational do- main give these methods a wider applicability. With the aid of Maple 12, we have assured the correctness of the obtained solutions by putting them back into the original equation. We hope that they will be useful for further studies in applied sciences and engineering.


The authors wish to thank the referee for his interesting suggestions and comments.


1 J. Scott Russell, Report on waves, Fourteenth meeting of the British Association for the Advancement of Science, (1844). [ Links ]

2 L. Yakushevich, Nonlinear Physics of DNA, (Wiley and Sons, 1998). [ Links ]

3 Mika Gustafsson, Coherent waves in DNA within the Peyrard-Bishop model, Master thesis, (Linköpings Universitet, 2003). [ Links ]

4 E. Villagran, Estructuras solitonicas y su influencia en la dinamica vibraional del ADN, PhD Thesis, (Universidad Autonoma del Estado de Mexico, Mexico, 2007). [ Links ]

5 Zidong Wang et al., Comput Methods Programs Biomed 111 (2013) 189-198. [ Links ]

6 M. Aguero, M. Najera and M. Carrillo, Int. J. Mod. Phys. B 22 (2008) 2571-2582. [ Links ]

7 G. Gaeta, J. Biol. Phys. 24 (1999) 81-96. [ Links ]

8 M. Peyrard, A.R. Bishop, Phys. Rev. Lett. 62 (1989) 2755-2758. [ Links ]

9 S. Dusuel, P. Michaux and M. Remoissenet, Physical Review E 57 (1998) 2320-2326. [ Links ]

10 A. Alvarez, S.R. Romero, J.F.R. Archilla, J. Cuevas, and P.V. Larsen, Eur. Phys. J. B 51 (2006) 119-130. [ Links ]

11 M. Najera , M. Carrillo , andM. Aguero , Adv. Studies Theor. Phys. 4 (2010) 495-510. [ Links ]

12 S.W. Englander, N.R. Kallenbach, A.J. Heeger, J.A. Krumhansl, and S. Litwin, Proc. Natl. Acad. Sci (U.S.A.) 77 (1980) 7222-7226. [ Links ]

13 S. Yomosa, J. Phys. Soc. Japan 51 (1982) 3318-3324. [ Links ]

14 S. Takeno, S. Homma, Prog. Theor. Phys. 70 (1983) 308-311. [ Links ]

15 C-T. Zhang, Phys. Rev. A 35 (1987) 886-891. [ Links ]

16 M. Wadati, Prog. Theor. Phys. Suppl. 59 (1976) 36-63. [ Links ]

17 H.W. Tam, X.B. Hu, Appl. Math. Lett. 15 (2002) 987-993. [ Links ]

18 A.M. Wazwaz, Comm. Nonlinear. Sci. Numer. Simulat. 13 (2008) 584-592. [ Links ]

19 Fan En-gui, Zhang Hong-qing, Phys. Lett. A 246 (1998) 403-406. [ Links ]

20 H. Hirota, A. Ramani, Phys. Lett. A 76 (1980) 95-96. [ Links ]

21 C.T. Yan, Phys. Lett. A 224 (1996) 77-84. [ Links ]

22 J.H. He, X.H. Wu, Chaos Solitons Fractals 30 (2006) 700-708. [ Links ]

23 Sheng Fu Deng, Bo Ling Guo, Ting Chun Wang, Sci. China Math 54 (2011) 555-572. [ Links ]

24 M. Wang, X. Li, J. Zhang, Phys. Lett. A 372 (2008) 417-423. [ Links ]

25 Shabnam Jamshidzadeh, Reza Abazari, Nonlinear Dyn 88 (2017) 2797-2805. [ Links ]

26 E. M. E. Zayed, J. Appl. Math. Comput 30 (2009) 89-103. [ Links ]

27 Reza Abazari , Comput. Math. Math. Phys. 53 (2013) 1371-1376. [ Links ]

28 Reza Abazari , Japan J. Indust. Appl. Math 31 (2014) 125-136. [ Links ]

29 E. M. E. Zayed , and Khaled A. Gepreel , J. Math. Phys. 50 013502 (2009); [ Links ]

30 E. M. E. Zayed , and Khaled A. Gepreel , Appl. Math. Comput. 212 (2009) 1-13. [ Links ]

31 E. M. E. Zayed , Journal of Information and Computing Science 8 (2013) 003-012. [ Links ]

32 W.-W. Li, Y. Tian, and Z. Zhang, Appl. Math. Comput. 219 (2012) 1135-1143. [ Links ]

33 JB Liu, KQ Yang, Chaos, Solitons and Fractals 22 (2004) 111-121. [ Links ]

34 E Yomba, Phys Lett A. 340 (2005) 149-160. [ Links ]

35 E. M. E. Zayed , A. H. Arnous, CHIN. PHYS. LETT. 29 (2012) 080203. [ Links ]

36 E. M. E. Zayed , A. H. Arnous , Scientific Research and Essays 8 (2013) 340-346. [ Links ]

37 Zayed EME, Arnous A. H., J. Part. Differ. Equ. 26 (2013) 373-384. [ Links ]

38 Fuding Xiea, Ying Zhanga, Zhuosheng Lüb, Chaos, Solitons Fractals 24 (2005) 257-263. [ Links ]

Received: July 02, 2018; Accepted: August 08, 2018

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License