Asymptotic Behavior of a Delay Differential Neoclassical Growth Model

A neoclassical growth model is examined with a special mound-shaped production function. Continuous time scales are assumed and a complete steady state and stability analysis is presented. Fixed delay is then assumed and it is shown how the asymptotic stability of the steady state is lost if the delay reaches a certain threshold, where Hopf bifurcation occurs. In the case of continuously distriubuted delays, we show that with small average delays stability is preserved, then lost at a threshold, then it is regained if the average delay becomes sufficiently large. The occurence of Hopf bifurcation is shown at both critical values.


Introduction
The examination of economic growth models is one of the most frequently discussed issues in mathematical economics.Day [1,2] has investigated a neoclassical growth model, and a productivity and population growth model and showed the emergence of complex behavior even under simple economic structure.His models were based on discrete time scales and a mound-shaped production function that represented the negative effect of pollution resulting from increasing capital.It was demonstrated by numerical computations that these models could generate cyclic and even chaotic OPEN ACCESS behavior.Following Day's pioneering works, a lot of effort has been given to the understanding of complex economic dynamics.Day [3], Puu [4] and Bischi et al. [5] present the earlier contributions of this field.A large number of studies assumed discrete time scales.Li and Yorke [6] have introduced the "period-three condition" to detect chaos, which has many applications in first-order nonlinear difference equations.The papers collected by Rosser [7] offer many applications.Only a few studies are devoted to the case of continuous time scales, since there is no general criterion to detect chaos and the system must have at least three dimensions.
In this paper, we will examine an extension of the neoclassical growth model, which can be traced back to the early works of Solow [8] and Swan [9].The neoclassical growth model is constructed with the two (usually implicitly mentioned) assumptions; one is the full-employment of labor and capital and the other is instantaneous adjustment in the output market.Thus, it is suitable for describing the long-run behavior of the economy.Due to the well-behaved production function, the steady state of the model is usually asymptotically stable.However, it is often observed in reality that growth path exhibits persistent fluctuations.The neoclassical model could be the good point of departure to show how such persistent behavior can emerge when nonlinearities and a production lag are present.Matsumoto and Szidarovszky [10] attempt to fill the gap and have introduced a neoclassical model with a mound-shaped production function that was assumed to be a Cobb-Douglas type function of the form  with  being the capital per unit labor.Although they show emergence of erratic fluctuations in the capital accumulation process, the production function is restrictive in the sense that it is defined only in the unit interval.This paper modifies this drawback, considers another type of mound-shaped production function and will examine the stability of the steady state with and without time delays in the continuous-time framework.Two kinds of delays will be discussed, fixed and continuously distributed (continuously hereafter) delays.We keep the relatively simple model of Matsumoto and Szidarovszky [10] in order to be able to compare the results and to illustrate that complex dynamics can be generated under simple economic assumptions with both function types.This paper develops as follows.First the mathematical model is formulated without time delays, and complete steady state and stability analysis is presented.Then, models with fixed delays and then with continuous delays are introduced and complete stability analysis is given.The last section concludes the paper.

The Mathematical Model
Matsumoto and Szidarovszky [10] have introduced a special growth model of the form where x is the capital per labor, s and  are positive parameters where The drawbacks of the neoclassical production function that ignores natural resource or energy are partly remedied in our function. of the exponential term reflecting a strength of a "negative effect" caused by increasing concentration of capital.The value of  is exogenously given, however, can be thought to be determined by a damaging degree of natural environment or energy resources.With the modified function, the mathematical model becomes where , ,    and   s = are positive parameters.The number of steady states and their locations depend on the specific values of the model parameters.We will consider three different cases:


. We can give the following interpretation for the value of the parameter γ of the production function; γ can be thought as a proxy for measuring returns to scale of the production function.Indeed, when x is small (i.e., the exponential term is close to unity), output increases more than unity, exactly unity and less than unity if γ > 1, γ = 1 and γ < 1, respectively.Let now ) (x f denote the right hand side of equation (1).If 0, = (0) x then the identically zero function is a solution which case is not interesting from the economic point of view and is eliminated from further considerations

Case I.
Assume first that 1. <


The steady states are the solutions of converges to  as x tends to zero with positive values.Hence so the positive steady state is the unique solution of equation The left hand side is zero at 0 = x and strictly increasing, furthermore, converges to  as x tends to infinity.The right hand side is so zero is a steady state and there is a unique root of the second factor, then the value of is decreasing and converges to zero with arbitrary then x decreases and converges to , x and if . Hence x is globally asymptotically stable.

Case III.
Consider finally the case of


The steady state equation has now the form so zero is a steady state again, and any other steady state is the solution of equation Notice that increases for x x < and decreases for .

> x x
Now we have three sub-cases.
then there is no positive steady state and with arbitrary 0, > (0) x ) (t x decreases and converges to zero.
(ii) if 0, < ) (x g then x x = is the unique positive steady state and   < 0 for all decreases and converges to 0, and if  0 >  then x decreases again and now converges to then equation ( 7) has two positive solutions, decreases and converges to .x That is, 1 x is locally unstable and 2 x is locally asymptotically stable.
remains at that steady state level for all 0. > t

Model with Fixed Delay
The fixed delay 0 > T is assumed in the second term of the right hand side of equation ( 1), so we have the following equation: where The local asymptotic behavior of the trajectory can be examined by linearization.Let x be a positive steady state.Then the linearized equation has the form is the deviation of Then x is locally asymptotically stable.
and at the steady state so the characteristic equation ( 11) can be rewritten as We also mention that the condition of Lemma 1 can be rewritten as In the special case of 1, =


this condition has the form In order to give a complete stability analysis, we have to find the possible stability switches.Substituting any stability switch, Separating the real and imaginary parts gives two equations, '  (16) and Adding the squares of these equations gives In order to have solution we have to assume now that that is, ( 14) is violated with strict inequalities.Concerning this assumption, we can give the following interpretation.Let F(x) be y.Then So this means that the absolute value of the elasticity of output with respect to capital is larger than unity.From (16) we have that if 0, > ) (x h '  then and by ( 13) and ( 19), cannot be zero.
By selecting T as the bifurcation parameter and implicitly differentiating the characteristic equation with respect to , T we have Therefore if a steady state is unstable with 0, = T then it remains unstable for all 0, > T and if a steady state is asymptotically stable at 0, = T then this stability is lost at and cannot be regained later.
That is, if the steady state is unstable without delay, then it remains unstable with any delay of positive length.If the steady state is asymptotically stable without delay, then it remains asymptotically stable until the delay reaches a certain threshold, and then becomes unstable and the stability cannot be regained later.Taking,  and 1 =  , we give an illustrative numerical example in Figure 1.The critical value and then the corresponding critical value of the delay is In Figure 1(A), the steady state is locally asymptotically stable in the dark-gray region with 1 >  z due to Lemma 1.It is also locally asymptotically stable in the light-gray region, which is under the critical curve and it is unstable in the white region above the curve.Setting c z z = and increasing T along the vertical dotted line in Figure 1(A), we can see that the steady state loses stability at c T T = .
Further increasing T , as observed in Figure 1(B), generates complex dynamics through a quasi period- doubling bifurcation in which T increases from 0.05 

Model with Continuously Distributed Delay
Assuming continuously distributed delays in the second term of equation ( 1) gives the following Volterra-type integro-differential equation: where 0 > T is a positive parameter, the average delay and 0  m is an integer.The kernel function has the form A similar model is investigated by Fanti and Manfredi [11] where 2  m is selected and the stability of the system with a cubic characteristic polynomial is examined based on the Routh-Hurwitz criterion.Stamova and Stamov [12] consider a generalized Solow model with endegenous labor growth and impulsive perturbations.Their stability analysis is based on the Lyapunov-Razumikhin sufficient stability conditions, which is a different approach than ours.
Linearizing equation ( 22), we have is the deviation of from the steady state level .
x We are looking for the solution in the usual exponential form This equation can be rewritten as Then similarly to the case of fixed delay we can prove the following result: Then x is locally asymptotically stable.
It is well known that the Routh-Hurwitz stability theorem provides necessary and sufficient conditions for all roots of a polynomial equation with real coefficients to have negative real parts.It is also known that it is difficult to locate the eigenvalues with analytic methods in general.However in some special cases, analytic results are still possible to obtain, as it will be next demonstrated.
then all coefficients are with a positive steady state, is locally asymptotically stable.Assume next that then the constant term is zero indicating that one eigenvalue is zero and the other is negative.So x is marginally stable in the linearized model, so no conclusion can be drawn about its asymptotical behavior in the nonlinear model.If  Then we have a cubic characteristic polynomial: then zero is an eigenvalue and the other two eigenvalues have negative real parts implying that x in the linearized system is marginally stable.Therefore no conclusion can be drawn about the stability of x in the nonlinear system.If , 1)/ ( < x then all coefficient of (25) are positive.In this case and when 1 


the Routh-Hurwitz criterion implies that the real parts of the eigenvalues are negative if and only if which can be reduced to a quadratic inequality in T : For the sake of simplicity, we re-introduce the notation which are positive and .< 2 and (26) holds if and only if and there are equal roots , so (26) holds and x is asymptotically stable.The instability region is shown in Figure 2(A) where z is the horizontal axis and T is the vertical axis.If we start with a very small value of T with any given 8, <  z then x is asymptotically stable.If we gradually increase , T then x remains asymptotically stable until it reaches the critical value ,

T T
If the steady state is unstable without delay, then it remains unstable with continuous delay with any T and m.If it is asymptotically stable without delay, then either it remains asymptotically stable with all T and m, or loses stability at a certain value of the average delay T and stability is regained with an even larger value of T and the steady state remains asymptotically stable afterwards.In such cases small and large average delays lead to asymptotically stable steady states.
We will next show that at the critical values  1 T and , 2  T Hopf bifurcation occurs giving the possibility of the birth of limit cycles.We select T as the bifurcation parameter.At the critical values (26) is satisfied with equality, so and the characteristic equation ( 25) can be rewritten as showing that there is a negative eigenvalue 1 and a pair of pure complex eigenvalues Consider  as a function of the bifurcation parameter T and differentiate implicitly equation (25) to have By simple calculation we can see that at at both values Hopf bifurcation occurs giving the possibility of the birth of limit cycles.We perform numerical simulations to illustrate the results obtained above.In Figure 2(A), the steady state is locally asymptotically stable in the dark-gray region with 1 >  z due to Lemma 2. It is also locally asymptotically stable in the light gray region and unstable in the white region when 1.

<  z
The appearance and disappearance of a limit cycle can be observed in Figure 2

T T
Under these specifications, the Volterra-type integro-differential equation ( 22) can be written as a 3D system of differential equations, with an increment of 0.01 , the steady state loses stability at point A and regains stability at point .

Conclusions
In this paper, a special neoclassical growth model was introduced and examined.A mound-shaped production function for capital growth was assumed in the dynamic equation.Zero is always a steady state, and depending on model parameters there is either no positive steady state, or one, or two positive steady states.A complete steady state analysis was followed by the derivation of stability conditions.By introducing fixed delay we demonstrated that stability can be lost at a certain value of the delay and the equilibrium remains unstable afterwards.In the case of continuously distributed delays it has been shown how stability can be lost at a certain value of the average delay and by further increasing the average delay it can be regained.At the critical values, Hopf bifurcation occurs giving the possibility of the birth of limit cycles.In our further study, more complex kernel functions will be considered and their effect on the asymptotic behavior of the steady state will be examined.
It is now well-known that any discrete-time dynamic model can generate complex dynamics involving chaos.However, in an aggregate model like the neoclassical growth model, it is natural to treat economic phenomena as continuous since the variables to be examined are the outcomes of a great number of decisions taken by different agents at different points of time.This is the main reason why this paper considers the neoclassical growth model in the continuous-time framework.Implication of the analysis is the following.Coexistence of nonlinearity and delay in production can be a source of persistent fluctuations in the capital-labor ration.

2
an increment of 0.01 and the local maximum and minimum of the corresponding trajectory are plotted against each value of T .
the most current value has the highest weight, which is then decreasing exponentially.If , 1  m then the most current value has zero weight which is then increasing until , = T s t and is decreasing exponentially afterwards.As m increases, the weighting function becomes more peaked around converges to the Dirac delta function.
reduces equation (22) with delays to equation (1) without delays.The asymptotic properties of this equation were already discussed earlier.kernel function becomes exponentially declining.Then characteristic equation (23) becomes quadratic,

x
then x is unstable and if the shape of the kernel function takes a bell-shaped form.
holds implying the asymptotical stability of the steady state.So we can assume that the left hand side of (26) has two roots
value of