_{1}

^{*}

*G*
_{0}
*/G*
_{1} “gaps” joint the
*S* phase and
* M* phase to form the cell cycle. The dynamics of enzyme reaction to drive the target protein production in
*M* phase is analyzed mathematically. Time delay is introduced since the signal transmission need time in
*G*
_{0}
*/G*
_{1} “gaps” phase. Hopf bifurcation of DDEs model is analyzed by applying geometrical analytical method. The instability oscillating periodic solutions arise as subcritical Hopf bifurcation occurs. The Hysteresis phenomena of the limit cycle are also observed underlying the saddle-node bifurcation of the limit cycle. Due to stability switching, interestingly, the bifurcating periodical solution dies out near the vicinity of Hopf lines. By Lyapunov-Schmidt reduction scheme, the normal form is computed on the center manifold. Finally, it is verified that the theory analytical results are in coincidence with the numerical simulation.

The important role relevance between modelling phase transition and complex dynamics of the cell cycle is found beyond the knowledge of biology background yet. The reason lies in that, on one respect, mathematical models of cell cycle can contribute to the understanding of its mechanism biologically; On another respect [_{0}/G_{1} phase [

In mathematical description of biological models, some key components and their mutant relationship with circumstance (inner-cell or inter-cells) are considered. People always translate their necessary biology knowledge into some differential equations to describe the growth stage of during cell cycle which include different phases [_{0}/G_{1} “gaps” joint the S phase and M phase to form the cell cycle. Hence, we introduce time delay to model the necessary time experienced to transpass “gaps”. Time delay is incorporated into the phosphate groups which while binding to target proteins to active protein phosphorylation process. In addition, to drive the downstream events to generate gene protein production, enzymes reaction formed by binding itself to Cdc2 monomer which is assumed no degradation in cyclin synthesis [

rate of loss of substrate = V max S K m + S

In cyclic synthesis, the rate of association with Cdk2 monomers depends on the probability of a collosion between two species. The enzyme activation/suppression leads to target protein production generated to form cyclin. Therefore, the protein phosphorylation and dephosphorylation activity dramatically altered the target protein which regulates kinase of cell growth.

Time delay describes the present state or state feedback dependent on the past history which becomes differential equations into DDEs of infinite dimensional [_{1} mitosis phase into a useful mathematical expression.

d y d t = − k 1 y + k 2 s 0 + y ( t − τ ) 2 θ 1 + y ( t − τ ) 2 (1.1)

Besides, one mother cell can be divided into two daughter cells during mitosis stage; the dynamics in G_{1}/M phase can also be expressed as

d x d t = − γ x + g ( x , y ) d y d t = δ − k 1 y + 2 e − r τ g ( x ( t − τ ) , y ( t − τ ) ) + k 3 ( y ( t − τ 2 ) − y ( t ) ) (1.2)

wherein x denotes enzyme concentration, and y represents gene protein concentration, τ is time delay with its valuable estimation between 2.8 h and 7.8 h. The associated function g ( x , y ) is written as

g ( x , y ) = k 2 s 0 + y 2 θ 1 + y 2 θ 2 θ 2 + x (1.3)

State feedback control is to introduce concentration difference into the model which reflects the time delay effects on dynamical evolution behavior. Herein τ 2 = τ + τ 1 , with τ 1 denotes time delay in “gap” junction as G_{0} phase.

System (1.2) is DDEs with multiple time delays which induce complex dynamical phenomena underlying Hopf bifurcation, such as the coexistence of four periodical solutions. It is observed that Neimake-Sacker bifurcation phenomena arise as varying time delay, and the quasi-periodic solution is also observed in coexistence regime of periodical solutions. By applying DDE-Biftool technique [

The saddle-node bifurcation of limit cycle occurs at some critical values of time delay τ . The stable periodical solutions are continued though the instability oscillation of bifurcating periodical solutions birth from both sub/super-critical Hopf point. The hysteresis phenomena of limit cycle are also observed beyond Hopf bifurcation. As is well-known, based on the strong continuous semigroup theory, the linear version of DDEs is described by the related form of ODE operator. Hopf bifurcation occurs as the infinitesimal generator has a pair of simple imaginary roots. A possible approach to investigate this bifurcation problem for DDEs involves the computation of normal form on center manifold. The preceeding work of normal form approach is famous with Faria and Magalhaes’s work, see also [

The whole paper is organized as follows: In Section 2, the bistable dynamics of steady states are analyzed. In Section 3, the imaginary roots of the related characteristic equation of the linear version of system (1.2) are analyzed to derive the critical values of Hopf bifurcation. With the introduction of a new time delay and define three different angle variables, the threshold of Hopf bifurcation lies in the intersection of two parameterized curves. The bifurcating periodical solution is continued as varying time delay continuously. In Section 4, the normal form near Hopf point is computed which further discovers that the bifurcating periodical solution is unstable, and this further verifies the numerical simulation results obtained in Section 3.

The bistability phenomena is observed as varying free parameter pairs ( k 1 , θ 1 ) to track fold bifurcation curves of the equilibrium solution. The cusp degenerate singularity is also tracked as the codimension-2 bifurcation point whilst fold curves on parameter plane immerge, as shown in

( Γ ) : { − k 1 x 3 − k 1 θ 1 x + x 2 + s 0 = 0 ; − 3 k 1 x 2 − k 1 θ 1 + 2 x = 0 ; (2.1)

As shown in

Hopf bifurcation triggers the phenomena of small amplitude periodic oscillation. Hopf bifurcation happens as one pair of characteristic roots cross the imaginary axis transversally as varying free parameter and time delay continuously. To study Hopf bifurcation, people used to linearize system near the positive equilibrium solution and further to analyze the roots with zero real part of the related characteristic equation.

In system (1.2), we assume the positive equilibrium solution is E * = ( x * , y * ) which satisfies

( − γ ( y * 2 + θ 1 ) x * + k 2 ( y * 2 + s 0 ) ) θ 2 − x * 2 γ ( y * 2 + θ 1 ) = 0 , 2 e − γ τ k 2 ( y * 2 + s 0 ) θ 2 − ( y * 2 + θ 1 ) ( k 1 y * − δ ) ( θ 2 + x * ) = 0 (3.1)

We plot x-nullcline and y-nullcline into x-y plane, the intersection of two nullclines is exactly the value level to describe the magnitude of the equilibrium solution, as shown in

d x d t = ( g x − γ ) x + g y y d y d t = 2 e − γ τ g x x ( t − τ ) − k 1 y + 2 e − γ τ g y y ( t − τ ) + k 3 ( y ( t − τ 2 ) − y ( t ) ) (3.2)

Therefore, the characteristic equation of the linearized version (3.2) is written as

Δ ( λ , k 1 , τ ) = | g x − γ − λ g y 2 e − ( λ + γ ) τ g x − k 1 − k 3 − λ + k 3 e − λ τ 2 + 2 e − ( λ + γ ) τ g y | = 0 (3.3)

Furthermore, one obtains that

Δ ( λ , k 1 , τ ) = P ( λ , c , τ ) e − ( λ + γ ) τ + S ( λ , c , τ ) e − λ τ 2 + Q ( λ , c , τ ) (3.4)

with

P ( λ , k 1 , τ ) = − 2 γ g y − 2 g y λ , S ( λ , c , τ ) = − k 3 λ , Q ( λ , k 1 , τ ) = γ k 1 + γ λ − g x k 1 − g x λ + k 1 λ + λ 2 − k 3 ( g x − γ − λ ) . (3.5)

where k 1 is the free parameter.

Suppose λ = i ω ( ω > 0 ) , substitute it into Equation (3.5) and separate the real part from the imaginary part to get

P R e − r τ cos ( ω τ ) + P I e − r τ sin ( ω τ ) + S R cos ( ω τ 2 ) + S I sin ( ω τ 2 ) + Q R = 0 , P I e − r τ cos ( ω τ ) − P R e − r τ sin ( ω τ ) − S R sin ( ω τ 2 ) + S I cos ( ω τ 2 ) + Q I = 0. (3.6)

with P R , P I , Q R , Q I represents the real part and the imaginary part of the polynomial P , Q respectively.

Solving cos ( ω τ ) , sin ( ω τ ) from Equations (3.6) to get

cos ( ω τ ) = − ( P I Q I + P R Q R ) + cos ( ω τ 2 ) ( P I S I + P R S R ) + sin ( ω τ 2 ) ( P R S I − P I S R ) e − r τ ( P I 2 + P R 2 ) , sin ( ω τ ) = − ( P I Q R − P R Q I ) + cos ( ω τ 2 ) ( P I S R − P R S I ) + sin ( ω τ 2 ) ( P I S I + P R S R ) e − r τ ( P I 2 + P R 2 ) (3.7)

Since in DDEs, some delays are specified which manifest very difficult dynamical properties as its corresponding characteristic equation is a transcendental equation. With the introduction of a new time delay, which is not really represented in DDEs, however, though undefined originally, a new time delay can be applied to solve unknown time delays. Accordingly, by introducing new time delay s and the associated angle variables θ = s ω , ϕ = ω τ , ψ = ω τ 2 , we define the new mapping

θ → θ 0 + 2 l π , ϕ = ϕ 0 + 2 k π , ψ = ψ 0 + 2 m π

with l , k , m ∈ N 0 + , then rewrite the polynomial P R , P I , Q R , Q I as

P R ( ω , k 1 , τ ) = P R ( θ s , k 1 , ϕ s θ ) , P I ( ω , k 1 , τ ) = P I ( θ s , k 1 , ϕ s θ ) , Q R ( ω , k 1 , τ ) = Q R ( θ s , k 1 , ϕ s θ ) , Q I ( ω , k 1 , τ ) = Q I ( θ s , k 1 , ϕ s θ ) , S R ( ω , k 1 , τ ) = S R ( θ s , k 1 , ϕ s θ ) , S I ( ω , c , τ ) = S I ( θ s , k 1 , ϕ s θ ) (3.9)

Define the new function

F ( θ , k 1 , s ) = ( P I Q I + P R Q R ) 2 + ( P I Q R − P R Q I ) 2 − e − 2 r ϕ s θ ( P I 2 + P R 2 ) 2 + 2 ( P I Q I + P R Q R ) ( P I S I + P R S R ) cos ( θ τ 2 s ) + 2 ( P I Q R − P R Q I ) ( P I S R − P R S I ) cos ( θ τ 2 s ) (3.10)

+ 2 ( P I Q I + P R Q R ) ( P R S I − P I S R ) sin ( θ τ 2 s ) + 2 ( P I Q R − P R Q I ) ( P I S I + P R S R ) sin ( θ τ 2 s ) + ( P I S I + P R S R ) 2 + ( P I S R − P R S I ) 2

By Equation (3.10), one defines curve L by

L : ϕ = − θ 2 r s ln H ( P R , P I , Q R , Q I , S R , S I , cos ( θ τ 2 s ) , sin ( θ τ 2 s ) ) ( P I 2 + P R 2 ) 2 (3.11)

with

H ( P R , P I , Q R , Q I , S R , S I , cos ( θ τ 2 s ) sin ( θ τ 2 s ) ) = ( P I Q I + P R Q R ) 2 + ( P I Q R − P R Q I ) 2 + 2 ( P I Q I + P R Q R ) ( P I S I + P R S R ) cos ( θ τ 2 s ) + 2 ( P I Q R − P R Q I ) ( P I S R − P R S I ) cos ( θ τ 2 s ) + 2 ( P I Q I + P R Q R ) ( P R S I − P I S R ) sin ( θ τ 2 s ) + 2 ( P I Q R − P R Q I ) ( P I S I + P R S R ) sin ( θ τ 2 s ) + ( P I S I + P R S R ) 2 + ( P I S R − P R S I ) 2

By Equation (3.11), we also get the description of curve S

S : ϕ = arctan U ( P R , P I , Q R , Q I , S R , S I , cos ( θ τ 2 s ) sin ( θ τ 2 s ) ) V ( P R , P I , Q R , Q I , S R , S I , cos ( θ τ 2 s ) sin ( θ τ 2 s ) ) − 2 k π (3.12)

with

U ( P R , P I , Q R , Q I , S R , S I , cos ( θ τ 2 s ) sin ( θ τ 2 s ) ) = ( P I Q R − P R Q I ) + cos ( θ τ 2 s ) ( P I S R − P R S I ) + sin ( θ τ 2 s ) ( P I S I + P R S R ) , V ( P R , P I , Q R , Q I , S R , S I , cos ( θ τ 2 s ) sin ( θ τ 2 s ) ) = ( P I Q I + P R Q R ) + cos ( θ τ 2 s ) ( P I S I + P R S R ) + sin ( θ τ 2 s ) ( P R S I − P I S R ) .

In geometry, Equation (3.11) and Equation (3.12) exhibit the infinite intersection or no intersection of curves L and S, and we list Hopf bifurcation condition as ϕ L = ϕ 0 + 2 l π and ϕ S + k π = ϕ L for some k and l. With the assumption that the intersection of two curves is expressed as ( θ * , ϕ * ) which satisfies ϕ S ( θ * ) + k π = ϕ L ( θ * ) , we compute the time delay τ * as the critical value of Hopf bifurcation, τ * = ϕ L s θ . Note that τ is represented by the introduction delay s. To calculate the Hopf bifurcation line, Equation (3.11) builts the inverse function c = c ( θ ) with θ 0 lying inside some subinterval of [ 0,2 π ] since we have defined the relationship between θ and θ 0 by mapping (3.8), and Hopf bifurcation occurs at

τ ( θ ) = s θ ( arctan U ( P R , P I , Q R , Q I , S R , S I , cos ( θ τ 2 s ) sin ( θ τ 2 s ) ) V ( P R , P I , Q R , Q I , S R , S I , cos ( θ τ 2 s ) sin ( θ τ 2 s ) ) − 2 k π ) (3.13)

Following we compute the transversal condition of Hopf bifurcation, differentiate the rightside of Equation (3.4) with respect to τ , one obtains

( ∂ P ∂ λ d λ d τ + ∂ P ∂ τ ) e − ( λ + r ) τ + P e − ( λ + r ) τ ( − d λ d τ τ − r − λ ) + ( ∂ S ∂ λ d λ d τ + ∂ S ∂ τ − τ 2 S d λ d τ ) e − λ τ 2 + ( ∂ Q ∂ λ d λ d τ + ∂ Q ∂ τ ) = 0 (3.14)

then by solving d λ d τ from Equation (3.14), we get

1 d λ d τ = ∂ P ∂ λ e − ( λ + r ) τ − P e − ( λ + r ) τ + ( ∂ S ∂ λ − τ 2 S ) e − λ τ 2 + ∂ Q ∂ λ ∂ P ∂ τ e − ( λ + r ) τ + P e − ( λ + r ) τ ( − r − λ ) + ∂ S ∂ τ e − λ τ 2 + ∂ Q ∂ τ (3.15)

differentiate the rightside of Equation (3.4) with respect to θ to further get

( ∂ P ∂ λ d λ d θ + ∂ P ∂ c d c d θ + ∂ P ∂ τ d τ d θ ) e − ( λ + r ) τ + P e − ( λ + r ) τ ( − d λ d θ τ − ( r + λ ) d τ d θ ) − τ 2 S d λ d θ e − λ τ 2 + ( ∂ S ∂ λ d λ d θ + ∂ S ∂ c d c d θ + ∂ S ∂ τ d τ d θ ) e − λ τ 2 + ( ∂ Q ∂ λ d λ d θ + ∂ Q ∂ c d c d θ + ∂ Q ∂ τ d τ d θ ) = 0 (3.16)

Noticed that λ = i ω and ω s = θ , we get d λ d θ = i 1 s , then substitute it into Equation (3.16) to get

( ∂ P ∂ λ 1 s i + ∂ P ∂ c d c d θ ) e − ( λ + r ) τ + P e − ( λ + r ) τ ( − i 1 s τ ) + S e − λ τ 2 ( − i 1 s τ 2 ) + ( ∂ S ∂ λ 1 s i + ∂ S ∂ c d c d θ ) e − λ τ 2 + ( ∂ Q ∂ λ 1 s i + ∂ Q ∂ c d c d θ ) = ( − ∂ P ∂ τ e − ( λ + r ) τ + P e − ( λ + r ) τ ( r + λ ) − ∂ S ∂ τ e − λ τ 2 − ∂ Q ∂ τ ) d τ d θ (3.17)

then by Equation (3.17) and Equation (3.15), one has

ℜ ( 1 d λ d τ | λ = i ω ) = − ℜ ( ( ∂ P ∂ λ e − ( λ + r ) τ − P τ e − ( λ + r ) τ − τ 2 S e − λ τ 2 + ∂ S ∂ λ e − λ τ 2 + ∂ Q ∂ λ ) ( a b c c . 1 s ( − i ) + s b c c . d c d θ ) ) (3.18)

herein a b c c . denotes the conjugate part of ∂ P ∂ λ e − ( λ + r ) τ − P τ e − ( λ + r ) τ − τ 2 S e − λ τ 2 + ∂ S ∂ λ e − λ τ 2 + ∂ Q ∂ λ , and s b c c . expresses the conjugate part of ∂ P ∂ c e − ( λ + r ) τ + ∂ S ∂ c e − λ τ 2 + ∂ Q ∂ c . Noticed the expression of the characteristic polynomial (3.4), one has

ℜ ( 1 d λ d τ | λ = i ω ) = − s i g n ℜ ( Δ λ ( λ , c , τ ) Δ ¯ c ( λ , c , τ ) ) s i g n ( d c d θ ) (3.19)

The stability property of the equilibrium solution changes into unstable state as ascending free parameter and stability switching phenomena occurs, as shown in

S satisfies Φ L = mod ( Φ S ,2 π ) , as shown in

Varying free parameter k 1 , Hopf bifurcation occurs at k 1 = 0.876 and k 1 = 0.157 , and periodical oscillation phenomena are observed which arise at Hopf points. By applying DDE-Biftool software, the bifurcating periodical solutions are calculated with high technique and which further simulate the periodical

solutions continuously as varying free parameter k 1 or time delay respectively. As shown in

Hopf bifurcation brings forth the stability switching phenomena of the equilibrium solution of system (1.2). On one respect, Hopf bifurcation partition the stable regime from the unstable regime hence the stable equilibrium solution

changes to be unstable when crossing over the margin of the stable region then switches back to be stable state again as further crossing over the margin of the unstable region. In addition, both the period-2 oscillating solution and the quasi-periodical solution manifests in system (1.2) as ascending time delay, as shown in Figures 7(a)-(d) and Figures 8(a)-(d). By Poincare section analysis, the switching phenomena between the period-2 solution and the quasi-periodical solution are also observed by numerical simulation technique, as shown in

It is verified that Hopf bifurcation occurs at parameter pairs ( k 3 , τ ) = ( k 3 * , τ c ) as shown in

With the introduction of 0 < ε ≪ 1 , set k 3 − k 3 c = ε k ε , τ − τ c = ε τ ε and scale x − x * = ε x , y − y * = ε y , then system (1.2) is transformed into its 3rd Taylor expansion truncation

d x d t = − γ x + g x x + g y y + ε 1 2 g x x x 2 + ε g x y x y + ε 1 2 g y y y 2 + ε 1 6 g x x x x 3 + ε 1 2 g x x y x 2 y + ε 1 2 g x y y x y 2 + ε 1 6 g y y y y 3

d y d t = 2 e − r τ c g x x ( t − τ c ) − 2 r ε τ ε e − r τ c g x x ( t − τ c ) − ( k 1 + k 3 c ) y

+ 2 e − r τ c g x ( x ( t − τ c − ε τ ε ) − x ( t − τ c ) ) + k 3 y ( t − τ 2 ) + ε k ε y ( t − τ 2 ) + 2 e − r τ c g y y ( t − τ c ) − 2 r ε τ ε g y y ( t − τ c ) + 2 e − r τ c g y ( y ( t − τ c − ε τ ε ) − y ( t − τ c ) )

+ ε e − r τ c g x x x ( t − τ c ) 2 + ε e − r τ c g x y x ( t − τ c ) y ( t − τ c ) + ε e − r τ c g y y y ( t − τ c ) 2 + ε 1 3 e − r τ c g x x x x ( t − τ c ) 3 + ε e − r τ c g x x y x ( t − τ c ) 2 y ( t − τ c ) + ε e − r τ c g x y y x ( t − τ c ) y ( t − τ c ) 2 + ε 1 3 e − r τ c g y y y y ( t − τ c ) 3 (4.1)

with initial value definition

( x ( t ) y ( t ) ) = ( ϕ 1 ϕ 2 ) , − τ m ≤ t ≤ 0

wherein τ m = max ( τ 2 , τ ) . According to the property of DDEs, solutions of Equation (4.1) are continuous on Banach space C = ( [ − τ m , 0 ] , R 2 ) . With few discontinuity jumps, the solution operator is differentiable on the extended phase space B C = ( [ − τ m , 0 ] , R 2 ) . With the assumption Z = ( x , y ) T , and denote Z ( t + s ) = Z t ( s ) for any s ∈ [ − τ m ,0 ] , then Equation (4.1) is written as

d Z ( t ) d t = L ( ε ) Z t + F ( Z t ( 0 ) , Z t ( − τ ) , z t ( − τ 2 ) , ε ) (4.2)

Based on the fundamental theory of DDEs, there is a bounded variational function to express the linear version of Equation (4.2) on the extended phase space. In addition, by Reize representation theorem, the corresponding 2 × 2 bounded matrix function is listed below to write the linear version into its integral operator form, that is, for any ϕ ∈ B C ,

L ( ε ) ϕ = ∫ − τ c 0 d η ( θ ) ϕ + ε ∫ − τ m 0 d η 1 ( θ ) ϕ + ∫ − τ m 0 d η 2 ( θ ) ϕ , (4.3)

with

d η ( θ ) = ( ( − γ + g x ) δ ( θ ) g y δ ( θ ) 2 e − r τ c g x δ ( θ + τ c ) − ( k 1 + k 3 c ) δ ( θ ) + k 3 δ ( θ + τ 2 ) + 2 e − r τ c g y δ ( θ + τ c ) )

d η 1 ( θ ) = ( 0 0 − 2 r τ ε e − r τ c g x δ ( θ + τ c ) k ε δ ( θ + τ 2 ) − 2 r τ ε g y δ ( θ + τ c ) )

d η 2 ( θ ) = ( 0 0 2 e − r τ c g x ( δ ( θ + τ c + ε τ ε ) − δ ( θ + τ c ) ) 2 e − r τ c g y ( δ ( θ + τ c + ε τ ε ) − δ ( θ + τ c ) ) )

and the nonlinear part F ( ⋅ ) is written as

F ( ϕ ( 0 ) , ϕ ( − τ ) ) = ( 1 2 g x x ϕ 1 ( 0 ) 2 + g x y ϕ 1 ( 0 ) ϕ 2 ( 0 ) + 1 2 g y y ϕ 2 ( 0 ) 2 + 1 6 g x x x ϕ 1 ( 0 ) 3 + 1 2 g x x y ϕ 1 ( 0 ) 2 ϕ 2 ( 0 ) + 1 2 g x y y ϕ 1 ( 0 ) ϕ 2 ( 0 ) 2 + 1 6 g y y y ϕ 2 ( 0 ) 3 , e − r τ c g x x ϕ 1 ( − τ c ) 2 + e − r τ c g x y ϕ 1 ( − τ c ) ϕ 2 ( − τ c ) + e − r τ c g y y ϕ 2 ( − τ c ) 2 + 1 3 e − r τ c g x x x ϕ 1 ( − τ c ) 3 + e − r τ c g x x y ϕ 1 ( − τ c ) 2 ϕ 2 ( − τ c ) + e − r τ c g x y y ϕ 1 ( − τ c ) ϕ 2 ( − τ c ) 2 + 1 3 e − r τ c g y y y ϕ 2 ( − τ c ) 3 ) (4.4)

Consider the linear operator (4.3), it’s an infinitesimal generator of the strong continuous semigroup in phase space BC, and for any ϕ ∈ B C , we define new linear operator A ( ε ) : B C → B C and its adjoint operator A * ( ε ) : B C * → B C * wherein B C * = C * ( [ 0, τ ] , R 2 ) , that is,

A ( ε ) ϕ = { d ϕ d θ , for − τ ≤ θ < 0 L ( ε ) ϕ , for θ = 0 (4.5)

and

A * ( ε ) ψ = { − d ψ d s , for 0 < s ≤ τ L * ( ε ) ψ , for s = 0 (4.6)

with

L * ( ε ) ψ = ∫ − τ 0 d η T ( s ) ψ ( − s ) (4.7)

For ϕ ∈ B C , ψ ∈ B C * , define the bilinear form

〈 ψ , ϕ 〉 = ψ ¯ T ( 0 ) ϕ ( 0 ) − ∫ − τ m 0 ∫ 0 θ ψ ¯ T ( ξ − θ ) d η ( θ ) ϕ ( ξ ) d ξ (4.8)

Based on the above analysis, Equation (4.2) can be written as its operator differential equations, that is

Z ′ t = A ( ε ) Z t + R Z t (4.8)

with

R ϕ = { 0 , for − τ ≤ θ < 0 F ( ϕ ( 0 ) , ϕ ( − τ ) , ε ) , for θ = 0 (4.9)

Considering the linear version of differential operator (4.8), Hopf bifurcation occurs at ε = 0 , and the associated characteristic roots set is denoted as Λ = { i ω , − i ω } as verrfied in Section 2. With the assumption of other eigenvalues with negative real parts, the phase space can be suspended on its restricted center manifold. Hence the eigenspace is decomposed into the eigenspace P corresponding to roots set Λ and its complementary subspace Q. We suppose eigenspace P is spanned by P = s p a n { q ( θ ) , q ¯ ( θ ) } , − τ ≤ θ ≤ 0 where q ( θ ) and its conjugation vector q ¯ ( θ ) are respectively the eigenvector of i ω and − i ω . We represent the eigenbasis Φ ( θ ) , Ψ ( s ) of Λ associated with linear operator A ( 0 ) and its adjoint operator A * ( 0 ) respectively, then

A ( 0 ) Φ ( θ ) = ( i ω 0 0 − i ω ) Φ ( θ ) , A * ( 0 ) Ψ T ( s ) = ( − i ω 0 0 i ω ) Ψ T ( s ) (4.10)

By the definition of Equation (4.8), we have 〈 Ψ , Φ 〉 = I .

With a possible discontinuity jump at θ = 0 , we also define the mapping on the phase space as the following

X 0 = { I for θ = 0 , 0 for − τ ≤ θ < 0 (4.11)

for any Z t ∈ B C , we can write it as Z t = ϕ + X 0 α , define the projection mapping Π : B C → P , then we have Π Z t = Φ ( θ ) [ 〈 Ψ , ϕ 〉 + Ψ ¯ T ( 0 ) w ] , therefore, we can write Z t as the expression Z t = Φ ( z z ¯ ) + Y t with Y t ∈ Q . Further computation express the differentiate relationship of axis variable z ( t ) given as

( z ′ ( t ) z ¯ ′ ( t ) ) = ( i ω 0 0 − i ω ) ( z z ¯ ) + Ψ ¯ T ( 0 ) ( A ( ε ) − A ( 0 ) ) Φ ( z z ¯ ) + Ψ ¯ T ( 0 ) R ( Φ ( z z ¯ ) + Y t ) ,

Y ′ t = A ( 0 ) Y t + ( I − Π ) X 0 R ( Φ ( z z ¯ ) + Y t ) = A ( 0 ) Y t + { − Φ ( θ ) Ψ ¯ T ( 0 ) F ( Φ ( 0 ) z + Y t ( 0 ) , Φ ( − τ ) z + Y t ( − τ ) ) , for − τ ≤ θ < 0 F ( Φ ( 0 ) z + Y t ( 0 ) , Φ ( − τ ) z + Y t ( − τ ) ) − Φ ( 0 ) Ψ ¯ T ( 0 ) F ( Φ ( 0 ) z + Y t ( 0 ) , Φ ( − τ ) z + Y t ( − τ ) ) , for θ = 0 (4.12)

Set B = ( i ω 0 0 − i ω ) , then Equation (4.12) is written into its Taylor expansion with 3rd truncation to be expressed as

( z ′ ( t ) z ¯ ′ ( t ) ) = B ( z z ¯ ) + f 2 ( 11 ) ( z , z ¯ , 0 , 0 , ε ) + f 2 ( 12 ) ( z , z ¯ , 0 , 0 , 0 ) + f 3 ( 1 ) ( z , z ¯ , y ( 0 ) , y ( − τ c ) , 0 ) , y ′ ( t ) = A ( 0 ) y + { − Φ ( θ ) f 2 ( 12 ) ( z , z ¯ , 0 , 0 , 0 ) , for − τ c ≤ θ < 0 f 2 ( 2 ) ( z , z ¯ , 0 , 0 , 0 ) − Φ ( 0 ) f 2 ( 12 ) ( z , z ¯ , 0 , 0 , 0 ) for θ = 0 (4.13)

and

f 2 ⋅ ( z , z ¯ , 0 , 0 , ε ) = ∑ i + j + k = 2 f i j 0 k ⋅ z i z ¯ j ε k , f 3 ⋅ ( z , z ¯ , 0 , 0 , 0 ) = ∑ i + j = 3 f i j 00 ⋅ z i z ¯ j

Define the operator M 2 11,12,2 on the homogeneous space H 2 ( z , z ¯ ) as

M 2 1 ⋅ ( p 1 ⋅ ( z , z ¯ ) p ¯ 1 ⋅ ( z , z ¯ ) ) = B ( p 1 ⋅ ( z , z ¯ ) p ¯ 1 ⋅ ( z , z ¯ ) ) − ( p z 1 ⋅ p z ¯ 1 ⋅ p ¯ z 1 ⋅ p ¯ z ¯ 1 ⋅ ) B ( z z ¯ ) , M 2 2 U ( z , z ¯ ) = A ( 0 ) U ( z , z ¯ ) − ( U z , U z ¯ ) B ( z z ¯ ) (4.14)

Then the normal form on the center manifold can be written as

( z ′ ( t ) z ¯ ′ ( t ) ) = B z + g 2 ( 1 ) ( z , z ¯ , 0 , 0 , ε ) + g 3 ( 1 ) ( z , z ¯ , 0 , 0 , 0 ) (4.15)

wherein

g 2 ( 1 ) ( z , z ¯ ,0,0, ε ) = P r o j I m M 2 11 c f 2 11 ( z , z ¯ ,0,0, ε ) g 3 ( 1 ) ( z , z ¯ ,0,0,0 ) = P r o j I m M 3 1 c f ˜ 3 ( 1 ) ( z , z ¯ ,0,0,0 ) (4.16)

with

f ˜ 3 ( 1 ) ( z , z ¯ , 0 , 0 , 0 ) = f 3 ( 1 ) ( z , z ¯ , 0 , 0 , 0 ) − ( p z 12 p z ¯ 12 p ¯ z 12 p ¯ z ¯ 12 ) f 2 ( 1 ) ( z , z ¯ , 0 , 0 , 0 ) + f 2 y 0 ( 12 ) ( z , z ¯ , 0 , 0 , 0 ) U ( z , z ¯ ) ( 0 ) + f 2 y τ ( 12 ) ( z , z ¯ , 0 , 0 , 0 ) U ( z , z ¯ ) ( − τ c ) (4.17)

Hence by (4.15), one obtains the formal norm near Hopf point

z ′ ( t ) = i ω z + ε k 10 z + k 12 z 2 z ¯ (4.18)

The bases of projection I m M 2 11 , I m M 2 12 are easily to compute as

( 2 i ω z ¯ ε 0 ) , ( 0 − 2 i ω z ε ) , ( − i ω z 2 0 ) , ( i ω z z ¯ 0 ) , ( 3 i ω z ¯ 2 0 ) , ( 0 − 3 i ω z 2 ) , ( 0 − i ω z z ¯ ) , ( 0 i ω z ¯ 2 )

The bases of I m ( M 1 3 ) are expressed as

( − 2 i ω z 3 0 ) , ( 2 i ω z z ¯ 2 0 ) , ( 4 i ω z ¯ 3 0 ) , ( 0 − 4 i ω z 3 ) , ( 0 − 2 i ω z 2 z ¯ ) , ( 0 2 i ω z ¯ 3 )

Hence, we get

k 10 = f 10001 , 1 ( 11 ) ( z , z ¯ , 0 , 0 ) (4.19)

By choosing

( p ( 12 ) ( z , z ¯ ) p ¯ ( 12 ) ( z , z ¯ ) ) = f 20000 ( 12 ) − i ω z 2 + f 11000 ( 12 ) i ω z z ¯ + f 02000 ( 12 ) 3 i ω z ¯ 2 , (4.20)

we eliminate the term f 2 ( 12 ) ( z , z ¯ ,0,0,0 ) . Set

U ( z , z ¯ ) = H 2000 z 2 + H 1100 z z ¯ + H 0200 z ¯ 2 (4.21)

then by Equation (4.13) and Equation (4.14), for − τ ≤ θ < 0 , we have

A ( 0 ) H 2000 ( θ ) = 2 i ω H 2000 ( θ ) − Φ ( θ ) f 20000 ( 12 ) , A ( 0 ) H 1100 ( θ ) = − Φ ( θ ) f 1100 ( 12 ) , A ( 0 ) H 0200 ( θ ) = − 2 i ω H 0200 ( θ ) − ϕ ( θ ) f 02000 ( 12 ) (4.22)

The initial condition of Equation (4.22) is also computed as

L H 2000 = 2 i ω H 2000 ( 0 ) + f 20000 ( 2 ) − Φ ( 0 ) f 20000 ( 12 ) , L H 1100 = f 11000 ( 2 ) − Φ ( 0 ) f 11000 ( 12 ) , L H 0200 = − 2 i ω H 0200 ( 0 ) + f 02000 ( 2 ) − Φ ( 0 ) f 02000 ( 12 ) (4.23)

Based on the initial condition (4.23), one can derive the coefficients H 2000 ( θ ) , H 1100 ( θ ) , H 0200 ( θ ) by the integral Equation (4.22) with respect to θ . Therefore, we get the coefficient k 12 in formal norm (4.18) to be listed as

k 12 = f 21000 , 1 ( 1 ) + f 20000 , 1 ( 12 ) f 11000 , 1 ( 12 ) i ω − | f 11000 , 1 ( 12 ) | ( 12 ) i ω + 2 | f 02000 , 1 ( 12 ) | 2 3 i ω + f 10100 , y 0 ( 12 ) H 1100 ( 0 ) + f 10010 ( 12 ) H 1100 ( − τ c ) + f 01100 , y 0 ( 12 ) H 2000 ( 0 ) + f 01010 , y τ ( 12 ) H 2000 ( − τ c ) (4.24)

Cell grows in size by double in mass then divide to form into two daughter cells. In such a way, cells segregate to form more descendant cells to make lifetime perpetual both in plants and animals. The mutation relationship between different reaction components during cells growth and signal transaction between phase transitions is ubiquitous in regulating cells growth and division. The introduction of time delay is a crucial factor in cells cycle growth model since signal transaction is experienced in G_{1}/G_{2} phase. The cells growth model was set forth by the enzyme reaction formed by binding itself to Cdk2 monomer to drive the downstream events to generate gene protein production. In addition, the protein phosphorylation and dephosphorylation activity dramatically altered the target protein which regulates kinase of cell growth. We focused on the stability analysis and bifurcation phenomena in cells growth system. Hopf bifurcation occurs as varying free parameter and time delay continuously. The co-existence periodical solutions were observed as varying parameter continuously near Hopf bifurcation point. The manifest stability switching phenomena also brings forth the quasi-periodical solution appearing in system. The complex dynamics as switching routes between period-2 solution and quasi-periodical solution were explored. We applied center manifold analytical scheme combined with the dimensional reduction method to analyze the normal form near Hopf bifurcation, which further verifies that Hopf bifurcation is sub-critical and the bifurcating periodical solution is unstable.

The author declares no conflicts of interest regarding the publication of this paper.

Ma, S.Q. (2021) Bifurcation Analysis of the Regulatory Modules of the Mammalian M/G_{1} Phase with Time Delay. International Journal of Modern Nonlinear Theory and Application, 10, 29-48. https://doi.org/10.4236/ijmnta.2021.102003