#### Abstract

We introduce a time-to-build technology in a Solow model with nonconstant population. Our analysis shows that the population dynamics may be a source of stability switches and Hopf bifurcations. The analytical results are obtained using the recent technique introduced by Beretta and Kuang (2002) in the studying of delayed differential equations with delay-dependent coefficients in characteristic equation. Numerical simulations are performed in order to illustrate the main dynamic features of the model.

#### 1. Introduction

The relationship between economic activity and population dynamics is a multifaceted issue. On the one hand, demographic changes affect many economic variables such as labour force [1] or public spending for pension system [2–4]; on the other hand, macroeconomic environment may influence fertility choices [5–7], life time expectation [8, 9], and emigration and immigration decisions [10].

Taking into account these aspects, many authors have built models in which nonlinear mechanisms are sources of oscillating trajectories of the main demoeconomic variables. Day at al. [11] show the possibility of chaotic dynamics in a discrete overlapping generations model in which agents are assumed to determine the number of children in a trade-off between the cost of childbearing and private consumption; Feichtinger and Sorger [12] introduce, according to Easterlin theory, some feedback mechanisms into population dynamics able to generate oscillating birth trajectories; Benhabib and Nishimura [13] and Feichtinger and Sorger [12] show the existence of nonlinear dynamics in neoclassical models of fertility; Manfredi and Fanti [14, 15] study the occurrence of limit cycles or complex behaviours in Goodwin-type models in which the dynamics of the population of workers is introduced; Fanti and Gori [16] study the possibility of cyclical instability in an overlapping generations model with pay-as-you-go pension scheme.

In this paper we analyze another source of endogenous fluctuations in demoeconomic variables. In particular we focus on the relation between neoclassical production, population dynamics, and stability of equilibrium path using a Solow-type framework in which population evolves at a constant (positive or negative) rate and a time-to-build technology is introduced. From a theoretical point of view, it is interesting to note that in standard versions of the Solow and Ramsey models without delays but with a positive rate of population growth, because of the assumption of constant returns to scale, (i) the size of population does not influence the long run growth rate of the economy (this property is not shared by many endogenous growth models à la Romer that suffer from an unrealistic scale effect (see the criticism and the alternative approach proposed by Jones [26] who studies a model with results really similar to the Solow model)); (ii) the rate of growth of population determines the long run rate of capital accumulation to which the economy converges; (iii) there exists a unique globally attracting equilibrium, and trajectories converge monotonically to it (see [17, 18]). Instead, if population is decreasing and the interior equilibrium acts as a watershed for the trajectories of the model and the ratio of capital to labour and according to initial conditions may converge to zero or be divergent. Nevertheless, in all cited works, it is worthy to note that the time series seems to be too regular in order to describe the real evolutions of economic variables.

Instead, if temporal lags in the production of capital goods are introduced and the rate of growth of population is different from zero, the dynamical system that describes the time evolution of per capita physical capital and the associated characteristic equation involve delay-dependent coefficients. In this case, the results are really different with respect to models (e.g., [19–22]) with fixed population size (and then with delay independent parameters). In particular stability switches as well as Hopf bifurcations may arise for realistic values of time delay.

From an analytical point of view, in order to characterize the local dynamics of the system, we will use the recent approach introduced by Beretta and Kuang [23] and based on the existence of real zeros of particular functions . We will show that if then the system is able to produce stability switches and Hopf bifurcations when time delay varies. The paper is organised as follows. Section 2 presents the model. Section 3 studies the dynamic phenomena; Section 4 concludes the paper.

#### 2. The Model

We consider a growth model à la Solow [24] in which there exists a delay of periods before capital becomes productive after the investment. This implies that, at time , the productive capital stock is . By assuming a Cobb-Douglas technology, we have that the total income at time is , where represents labor force at time and is the capital’s share. Departing from Zak [22] we assume that population grows at a constant rate . Let be the constant saving rate and the constant proportion of the capital stock that depreciates during production. From the identities and , where , , and are consumption, saving, and investment, respectively, all evaluated at time , it follows that the equation of capital accumulation is In order to simplify the notation we omit the indication of time dependence for variables and derivatives referred at time , and we use to indicate the state of the variable at time . Setting , the evolution of per capita physical capital over time is given by Normalizing the number of people at time zero to one yields . Then, substituting this in the above equation, we derive that our model is described by the following delay differential equation with delay-dependent coefficients

The equilibria of (3) are obtained by setting and for all . Doing this, one finds that the nontrivial equilibria of (3) must satisfy

Lemma 1. *For every , there exist a unique positive equilibrium if and a unique positive equilibrium if and .*

Using Taylor expansion on the right-hand side of (3) around , we get the linearized equation From this, we obtain the characteristic equation Equation (6) is a transcendental equation which has an infinite number of complex roots. For , the characteristic equation (6) becomes . Therefore, the equilibrium point is locally asymptotically stable. We have discussed (3) in the absence of delay. In the next section, we will investigate our model when delay is present and try to analyze the effect of the gestation time lag.

#### 3. Stability and Existence of Hopf Bifurcation

Now we assume and regard it as a bifurcation parameter to obtain results on the stability of . It is well known that the change of stability of will occur for . First, is not a root of the characteristic equation (6). Otherwise, we would have , yielding the absurd . Second, we note that (6) takes the form of an exponential polynomial with a coefficient depending on . Thus, we use the method introduced by Beretta and Kuang [23], which gives the existence of purely imaginary roots of a characteristic equation with delay-dependent coefficients. To find out the position of change of stability we substitute in (6) and then separating the real and imaginary parts we get By taking square of both sides of system (7) and then adding them up, one obtains the following equation:

Lemma 2. *Let . Then (8) has a positive root given by
**
only if the following condition holds:
*

*Proof. *Solving (8) for leads to (9), which is well defined if . When , this means ; that is, and . The first of these two inequalities yields . Let . Then the above condition reduces to . The conclusion follows.

Corollary 3. *(1) Let . Then condition (10) is equivalent to .**(2) Let . Then condition (10) is equivalent to
*

*Proof. *Let . Then (10) is . Hence, we get , which is always satisfied being in this case. Let . Then (10) becomes . A direct calculation yields the statement.

Define with Note that since the logarithm argument is greater than one. According to Beretta and Kuang’s [23] procedure, we now substitute (9) into system (7) and define the angle for as solution of A necessary and sufficient condition for , solution of (8), to be characteristic roots of (6) is that the arguments in system (7) and in system (14) are in the relationship Then, we can define the functions whose explicit expressions can be obtained from (14) as follows: Therefore, , , is a purely imaginary root of (6) if and only if is a zero of the functions , where

Proposition 4. *The characteristic equation (6) has a pair of simple pure imaginary roots , , at , provided for some . Moreover, this pair of simple conjugate purely imaginary roots crosses the imaginary axis from left to right if and crosses the imaginary axis from right to left if , where
*

*Proof. *First, we prove that is a simple root of (6). If is a repeated root, then holds true, leading to , which is a contradiction. Differentiating (6) with respect to , we find
Using (6) and (8), one has
Noticing that (8) yields
the above rewrites as
Hence,
Differentiating (18) with respect to and using (17), it follows
Since gives , comparing (25) evaluated at with (24), the conclusion holds.

We can easily observe the following properties of .(1), for and . Therefore, has no zeros in for any if has no zeros in .(2)From (9), we have . Hence, .(3)As , , we derive as .(4)As , , we get as .

Without loss of generality, we may suppose that if has positive zeros, then at these points; that is, the zeros are simple. Applying Hopf bifurcation theorem for functional differential equations (see Hale’s book [25]), we can conclude the existence of a Hopf bifurcation.

Summing up all the above analysis, we can state the main results in the next two Theorems.

Theorem 5. *Let .*(i)*A stability switch always exists in since the function has at least one positive simple zero in .*(ii)*The function has a finite number of positive simple zeros in , for any and the number of these roots is odd; say , respectively.(ii _{1})Equation (3) undergoes a Hopf bifurcation at the equilibrium when , .(ii_{2}) is locally asymptotically stable for and unstable if .*

Theorem 6. *Let .*(i)*If holds or if and the function has no positive simple zeros in , then the equilibrium is locally asymptotically stable for all .*(ii)*If holds, then the function has a finite number of positive simple zeros in , for any fixed and the number of these roots is even; say , respectively.(ii _{1})Equation (3) undergoes a Hopf bifurcation at the equilibrium when , .(ii_{2}) is locally asymptotically stable for and and unstable if .*

In order to clarify the interesting dynamic properties of the system stated in Theorems 5 and 6 we now resort to some numerical experiments by taking the following configurations of parameters:

In the whole section may be interpreted as the number of months before capital can be used for production. We first consider (26). Figure 1 shows the existence of a unique zero for . According to Theorem 5 the equilibrium is locally asymptotically stable for and undergoes a Hopf bifurcation for . Just after this value a closed invariant curve exists. For large enough value of time delay no feasible trajectories exist.

**(a)**

**(b)**

**(c)**

**(d)**

We now consider (27). If is large enough () we are in case (i) of Theorem 6 and no zero for exists and the equilibrium is locally asymptotically stable. Numerically we verify that no zero exists until , where is the value of alpha such that is tangent to horizontal axis (see Figure 2(a)).

**(a)**

**(b)**

For has two positive zeros. We set (see Figure 2(b)). The zeros of are and .

According to case (ii) in Theorem 6, we have that, for all , the equilibrium is stable; for the equilibrium undergoes a supercritical Hopf bifurcation and for a closed curve captures trajectories starting close enough to the unstable equilibrium; finally, for equilibrium undergoes a second Hopf bifurcation and the equilibrium becomes stable again (see the set of Figure 3). No other bifurcation occurs for larger values of .

**(a)**

**(b)**

**(c)**

**(d)**

#### 4. Conclusions

In this paper we have presented a theoretical explanation of economic cycles due to demographic reasons. Our model makes use of the Solow model with exogenous saving rate, exogenous population growth, and time-to-build technology. The economic and demographic variables are related through production function. The model proposed in this paper has several features that distinguish it from other approaches to nonlinear evolution of demoeconomic variables. Different from some contributions on the topic (see, e.g., [1]), we do not introduce ad hoc feedback mechanisms in order to describe non linear phenomena. Furthermore, we do not model other sources of nonlinearities such as the age structure (see [12]) of population or behavioural assumptions on fertility rate or optimizing decisions on the number of children (see [13]). Nonetheless, comparable results in terms of dynamical properties are obtained. Thus, we can conclude that the existence of persistent cycles seems to be a quite general result because it relays on the coexistence of capital and nonconstant labour in production of goods. From a mathematical point of view, we have shown that the introduction of population dynamics drastically changes the results of the model with respect to the case (see [22]). If , the system is able to produce stability switches and Hopf bifurcations when time delay varies not observable for . In particular we have seen that when , time delay may play a destabilising role for intermediate values of the parameter and eventually it may stabilise.