Nonlinear economic growth dynamics in the context of a military arms race
Daniel Metz and Adrian Viorel
Abstract. In the present contribution, we propose and analyze a dynamical eco- nomic growth model for two rival countries that engage an arms race. Under natural assumptions, we prove that global solutions exist and discuss their as- ymptotic long-time behavior. The results of our stability analysis support the recurring hypothesis in Cold War political science that engaging in an arms race with a technologically superior and hence faster growing adversary has damaging economic consequences. Numerical findings illustrate our claims.
Mathematics Subject Classification (2010):91B62, 91F10.
Keywords:Solow-Swan model, arms race, asymptotic behavior.
The Cold War has ended three decades ago and regional antagonisms have re- placed the previous colossal struggle. Nevertheless, our understanding of the United States vs. USSR Arms Race and its wider economic consequences remains far from complete (see ).
A plausible, and often repeated explanation attributes the Eastern Block’s col- lapse to an economic crisis triggered by unsustainable military ambitions. Given the ever growing and ever more visible gap in technological and economic capabilities, matching American military development was possible only at the expense of eco- nomic growth and stability.
The aim of the present contribution is to examine this hypothesis from an an- alytic point of view by developing a model that, at least qualitatively, reproduces economic stagnation caused by a prolonged military rivalry with a faster developing adversary.
Arms races have a long history that goes back far beyond the Cold War Era.
The ancient Greeks and Romans built fleets to match their Persian and Carthaginian rivals, but the naval race that followed the 1889 Naval Defense Act calling for the Royal
Navy to be as strong as the world’s next two largest navies combined, is probably the most intensively studied predecessor as it led to World War I (WWI).
For decades, Arms races have been a topical subject in Political Science, such that the scarcity of treatments from a mathematical perspective comes as a real surprise - all the more so given L. F. Richardson’s pioneering contributions to the field synthesized in Arms and Insecurity: A Mathematical Study of the Causes and Origins of War  andStatistics of Deadly Quarrels .
The classical Richardson model, which relies on a system of two coupled linear differential equations, has dominated theoretical debates for more than half a century.
If x(t) and y(t) denote the levels of arms for two rival states, with rates of change driven by the sum of a positive reaction to the other country’s arms, a negative
’fatigue’ reaction to own military level and a constant ’grievance’ term, then the time evolution is described by
dt(t) =−β11x(t) +β12y(t) +γ1, dy
dt(t) = β21x(t)−β22y(t) +γ2.
The unique equilibrium point of the system, which exists provided that the two straight lines defined by the right hand side of (1.1) are not parallel, may be un- stable and Richardson related exponentially diverging solutions with the outbreak of war. Nonlinear extensions of the classical Richardson model have been considered by Hill .
We take a similar approach but augment the model by adding an economic dimension described in terms of Solow-Swan dynamics discussed below based on .
A different line of thought, that we don’t pursue here, deals with arms races or, more generally, strategic interactions from a game theoretical perspective. Two or more actors play a (repeated) game in which the strategies that they can choose from are to arm or not to arm (see, for example ). For a recent contribution that is somewhat pertaining to the present work, we refer to . The direction contrary to our study, that is disarmament models has also been pursued (cf. ), while a strongly misleading use of the term arms race in a biological context has been rendered popular by Dawkins and Krebs in .
The Solow-Swan model, originating from the independent works  and , explains long-run economic growth in a neoclassical framework by relating capital, labor and technology. The model relies on three fundamental assumptions. The first assumption is anexponential population (or labor) growth
L(t) =L0ent (equivalently dL dt =nL).
The second assumption concerns aCobb-Douglas production functionconnecting the economic outputY to the laborL, capitalKand the level of technologyA
Y =AKαL1−α. Here, α∈[0,1] is the returns to scale constant.
The third assumption of the model asserts thatchange in capitalK(t) is due to the positive capital output saving (with saving rate σ ∈(0,1)) and to the negative capital depreciation (at a rateδ)
dt =σY −δK.
By combining these three assumptions and expressing them in terms of the capital intensity
k(t) = K(t) L(t) one obtains
dt L+k nL=dK
dt =σA(kL)αL1−α−δ(kL) and deduces the fundamental equation of the Solow-Swan model
dt =σAkα−(n+δ)k. (1.2)
Observe that here we have reached an explicitly solvable Bernoulli equation which defines a dynamical system with two equilibria,k∗= 0 being unstable in contrast to
which is asymptotically stable (attractor).
2. An economic growth model with arms race military expenses
The simplicity and lack of specificity proved to be both a strength and weak- ness of the classical Richardson model which has become a cornerstone of strategic thinking despite the somewhat imprecise concept of arms not allowing rigorous fit- ting to measurable data. It turns out that replacing weapon quantities by an abstract
’security’ concept that can be linked to economic factors is more lucrative. Loosely following discrete models in both Krabs  and Larrosa  we consider an augmented arms race model in underlying economic growth context
dt (t) =−k2(t)s1(t) +k1(t)s2(t), ds2
dt (t) = k2(t)s1(t)−k1(t)s2(t), dk1
dt (t) =a1k1(t)α−bk1(t)−cs2(t)k1(t), dk2
dt (t) =a2k2(t)α−bk2(t)−cs1(t)k2(t).
Here, si(t) describes the level of security of the state i at time t, ki(t) being the country’s capital intensity. Security levels obey Richardson-type equations but with
time-varying coefficients. The competitive nature of the model is reflected in the fact that an increase in one actor’s security is its adversary’s security loss as
dt (t) =−ds2
dt (t). (2.2)
In other words the total security is constant s1(t) +s2(t) = const and imposing s1(0) +s2(0) = 1 will assure, as we will see in the next section, 0≤s1(t), s2(t)≤1 meaning that security levels range from 0 (totally insecure) to 1 (totally secure).
On the other hand, both economies grow according to a Solow-Swan model with an additional term explicitly accounting for military expenses. These exchange terms including the adversary’s security might look surprising at first glance, but in view of (2.2) one country’s security is the other’s insecuritysj= 1−siand military expenses are proportional precisely to the insecurity 1−si. The coefficientc∈(0,1) represents a budget constraint and expresses the percentual limit which military spendings cannot exceed in a functional peacetime economy.
Returning to the security equations, one can now see that the right hand side terms are actually proportional to military costs, insecurity rising based on rival spending and decreasing based on own spendings.
The parametersα, bandairetain their original Solow model meaning and only a1, a2 differ from country to country. In view of (1.2), this difference is essential to our model and accounts for the technological gap separating the two economies.
3. Analysis of the model
We start our analysis by discussing an uncoupled Solow-Swan model with vari- able military expenditures. Quite naturally,the best and worst case scenarios, namely zero or maximal military spending, provide upper and lower bounds for the dynamics.
Lemma 3.1 (upper and lower bounds). Let us consider the initial value problem dk
dt =akα−bk−cs(t)k, k(0) =k0 (3.1) with coefficientsa >0, α∈(0,1), b >0, c∈[0,1] ands: [0,∞)→R a given smooth function withs(t)∈[0,1]for any t≥0. Ifk0>0 then the solution of (3.1)exists, is positive and satisfies for all times
a) k(t)≤k(t), wherekis the solution of dk
dt =akα−bk, k(0) =k0; (3.2)
a) k(t)≤k(t), wherekis the solution of dk
dt =akα−bk−ck, k(0) =k0. (3.3) Proof. In (3.1) we are dealing with a Bernoulli equation which is exactly solvable, hence the global existence using the usual change substitutionz(t) =k(t)1−α. Using
the variation of constants formula one has the desired positivity from k(t)1−α=z(t) =e−R0tb+cs(ρ)1−α dρ
k1−α0 + a 1−α
Z t 0
eR0τ b+cs(ρ)1−α dρdτ
To obtain both the upper and lower bounds, one can rely on standard sub and super- solution arguments. Since 0≤s(t)≤1 andk(t)≥0
and the conclusion follows.
Remark 3.2. From a dynamical systems point of view, both autonmous equations in Lemma 3.1 are Solow-Swan equations and have the same stability behavior albeit with different nonzero asymptotically stable equilibria namely
and k∗= a
After this helpful preliminaries we are in position to prove the existence of global solutions to (2.1).
Theorem 3.3 (global existence). Let us consider the growth under arms race rivalry model (2.1) with a1, a2 > 0, α ∈ (0,1), b > 0 and c ∈ [0,1]. Then for any initial conditions k1(0), k2(0) > 0 and s1(0), s2(0) >0 with s1(0) +s2(0) = 1 there exists a unique classical solution of the initial value problem associated to the system (2.1) which remains bounded for allt≥0.
Proof. We divide the proof in several steps.
Step 1. Local existence. As the right hand side of the system has good regularity (only local Lipschitz continuity is actually required), a standard Banach fixed point argument guarantees the existence of local in time solutions, defined on a maximal intervalt∈[0, T),T =T(s1(0), s2(0), k1(0), k2(0).
Step 2. Positivity ofk1andk2.Based on the representation formula (3.4) which holds on their maximal interval of existencet∈[0, T), one can see that for positive initial statesk1(0), k2(0)>0, bothk1(t) andk2(t) must be positive for t∈[0, T).
Step 3. Positivity ofs1ands2.Using the fact that
s1(t) +s2(t) = 1 for all t∈[0, T), (3.5) one can rewrite the evolution equations fors1 ands2as
dt =k1(t)−(k1(t) +k2(t))s1, ds2
dt =k2(t)−(k1(t) +k2(t))s2,
such that aplying the variation of constants formula one again s1(t) =e−R0t(k1(ρ)+k2(ρ))dρ
Z t 0
which is positive for positivek1. Similarly, one can show thats2has the same property.
In view of (3.5) and the positivity ofs1, s2we have
0≤s1(t), s2(t)≤1 for all t∈[0, T). (3.8) Step 4. Boundedness of k1, k2 and global solutions. The estimates in (3.8) not only assure the boundedness fors1, s2but also allow us to apply Lemma 3.1, more precisely the upper bound in a), and hence deduce the boundedness ofk1 and k2. A classical result (see Barbu ) now assures that the local in time solutions can be extend to
arbitrary positive times.
4. Asymptotic behavior of the model
We start by determining the equilibrium points of the system, that is the solu- tions of
0 =−k2∗s∗1+k∗1s∗2, 0 = k2∗s∗1−k∗1s∗2,
0 =a1(k1∗)α−bk1∗−cs∗2k1∗, 0 =a2(k2∗)α−bk2∗−cs∗1k2∗.
One can reduce this to a 3 by 3 nonlinear system by assuming that (3.5) holds. The resulting equilibrium equations are
0 =a1(k∗1)α−bk1∗−c(1−)s∗1)k1∗, 0 =a2(k∗2)α−bk2∗−cs∗1k∗2.
Trivial equilibria, that is with k∗i = 0, exist but are not interesting from a modeling perspective as they would indicate the disappearance of an economy. Nevertheless, we note without going into details, that all such equilibria are unstable, as a natural consequence of the lower bound b) in Lemma 3.1 means that both k1 and k2 are pushed away from zero even when starting arbitrarily close.
However, there exists also a nontrivial equilibrium point.
4.1. The unique nontrivial equilibrium In terms of the convenient notation
R∗=s∗1 s∗2 = k1∗
k2∗, from (3.5) we have
1 +R∗ and s∗2= 1
1 +R∗. (4.3)
Such that inserting this in the 3rd and 4thequation of (4.1) leads to k∗1
·b+ (b+c)R∗ (b+c) +bR∗,
F(R) =b+ (b+c)R (b+c) +bR.
Hence, finding the nontrivial equilibrium reduces to solving the coincidence problem (4.4) (or equivalently the fixed point problemR∗=aa1
2F(R∗)1/1−α). The existence of a unique nontrivial coincidence pointR∗follows from geometric considerations.
Lemma 4.1. The coincidence problem(4.4)witha1≥a2has a unique solutionR∗≥1 provided that α∈
and 2b+cc ≤1−α. If a1> a2 thenR∗> aa1
Proof.We start by noting some geometric properties ofF. One can easily check that F(0) = b
b+c, F(1) = 1 and F(∞) = b+c b while
F0(R) = (2b+c)c
((b+c) +bR)2 >0 and F00(R)<0.
In other words, F is monotonically increasing, convex and bounded from above by b+c/b. Consequently, the range ofR7→ aa1
2F(R) is [aa1
b+c b ].
Now, observe that for aa1
2 = 1,R= 1 is a solution of R1−α=F(R).
For the moment, we assume that there are no other solutions in (0,1) and later give a sufficient condition for this to hold true.
IfR1−α =F(R) has no solutions is (0,1), that is,F(R)> R1−α forR ∈(0,1) due to F(0) =b/(b+c) >0, then aa1
2R1−α =F(R) has no solutions in (0, a1/a2).
Indeed, on one hand
R1−α< F(R)≤ a1
a2F(R) for 0< R <1 while on the other hand
R1−α< R < a1 a2
< a1 a2
F(R) for 1< R < a1 a2
AsR7→R1−αis increasing and unbounded while F is increasing and bounded, the two curves will cross at a unique pointR∗>aa1
The necessary and sufficient condition forR1−α=F(R) to have no solutions in (0,1) is again geometric in nature. Actually, the slope ofF atR= 1 must not exceed that ofG(R) =R1−α, that is precisely
Remark 4.2. Both essential conditions for the uniqueness of the coincidence point α < 1
2 and c
2b+c ≤1−α (4.5)
are natural and in accordance with econometric data. The returns to scale constant is generally considered to beα≈1/3 while the depreciation constant isb≈0.05 (see Acemoglu ). On the other hand, even at the hight of the Cold War, according to the World Bank1 military expenses have not exceeded 10% of GDP, so roughly c ≈2b, which satisfies the coincidence condition.
0 2 4 6 8 10
0 1 2 3 4
R1 − α F(R) 1.5F(R) 1.75F(R) 2F(R)
Figure 1. Qualitative behavior of the coincidence problem under the uniqueness assumptions (4.5).
In the sequel, we analyze the stability of the equilibrium point corresponding to this uniqueR∗, that is of (4.3) together with
and k2∗= a2
To this end, we compute the Jacobi matrix of the (3 by 3) system which is
J(s∗1, k∗1, k2∗) =
−(k∗1+k∗2) 1−s∗1 −s∗1
ck∗1 T1 0
−ck∗2 0 T2
1 +R∗ and T2= (α−1)b+ (b+c)R∗ 1 +R∗ .
We discuss the eigenvalues of this matrix in two different, parameter-dependent cases.
4.2. The catch-up scenarioa1=a2=a
This is the simpler yet less realistic situation in which there exists no difference between the parameters describing the two countries, this especially means that both economies have the same technology level, and only their initial states may differ.
Returning to (4.4), one can see that it reduces to the simpler (R∗)1−α=F(R∗)
which has the unique coincidence pointR∗= 1. As a consequence s∗1=s∗2= 1
2 and k∗1=k∗2=k∗= a
and straightforward but rather tedious computation show that all three eigenvalues of the Jacobian J(12, k∗, k∗) have negative real part, so the equilibrium is locally asymptotically stable.
From a modeling perspective, this describes a catch-up evolution in which the country with the initially weaker economy will recover the deficit in the long-run and stabilize at the same level as its rival, as depicted in Figure 2.
0 20 40 60 80 100
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 k1(t)
0 20 40 60 80 100
0.0 0.2 0.4 0.6 0.8 1.0
Figure 2. The catch-up scenarioa1 =a2. The two countries expe- rience a convergent economic growth with the initialy weaker econ- omy catching up to the stronger (Left panel). Security levels also converge towards a balanced stationary state (Right panel). Simula- tions correspond to a1 =a2 = 0.15, α = 1/3, b = 0.06, c = 0.1 and s1(0) =s2(0) = 1/2, k1(0) = 0.3, k2(0) =k1(0)/4.
4.3. The increasing gap scenarioa1> a2
From our point of view, the more interesting and realistic situation is that of unequal coefficients a1 > a2. This describes a technological gap between the two contenders, and we will show that the quotient a1/a2 plays a decisive role in the long-term dynamics as its affects the equilibrium quotient R∗ of the two economies (see Figure 3).
0 20 40 60 80 100 0.0
0.5 1.0 1.5 2.0 2.5
0 20 40 60 80 100
0.0 0.2 0.4 0.6 0.8 1.0
Figure 3. The increasing gap scenarioa1> a2. Despite growth for both countries, the gap separating them widens as excessive military spending harms the slower developing one (Left panel). Security levels also separate in the stronger economy’s favour (Right panel). Simu- lations correspond toa1= 0.18, a2= 0.15, α= 1/3, b= 0.06, c= 0.1 ands1(0) =s2(0) = 1/2, k1(0) = 0.3, k2(0) =k1(0)/4.
Indeed, from the proof of Lemma 4.1, we know that the coincidence pointR∗ must lie above the aa1
2 threshold. This means that the equlibirum quotient exceeds the quotient of coefficients
Again, the eigenvalues of the Jacobian at the equilibrium point corresponding toR∗> aa1
2 have negative real part and hence the equilibrium point is asymptotically stable. We omit the details of this technical computation, but in order to strike a bal- ance between the abstract and the concrete level, we mention that given the realistic values
a1= 0.18, a2= 0.15, α= 1/3, b= 0.06, c= 0.1 for the parameters, one obtains
R∗= 2.208, k1∗= 2.773 and k2∗= 1.256 such that the eigenvalues of the Jacobi matrix are all negative
λ1,2,3=−4.081, −0.076 and −0.019.
In order to describe the economic implications of a prolonged military rivalry, we have constructed a nonlinear dynamical model that merges the classical Richardson arms race evolution with economic growth in the sense of Solow’s pioneering work.
The ensuing nonstandard model turns out to be well-posed and in accordance with both political and economic intuitions. More precisely, when considering different levels of technology for the two competing powers, the model predicts that in the long
run, due to nonlinear effects, the sizes of the two economies will be separated by a gap that exceeds the technology gap.
The reality of Cold War dynamics has been far more complex than the relatively simple model that we propose can describe. Many extensions are possible and, actually desirable. The most natural extension would be to consider the augmented human capital version of the Solow-Swan model due to Mankiw, Romer and Weil  not the Solow-Swan economic growth model itself.
Furthermore, the Cold War arms race is just a prototype for more general eco- nomic rivalry phenomena. Trade or economic wars provide very interesting challenges form a modeling perspective.
 Acemoglu, D., Introduction to Modern Economic Growth, Princeton University Press, 2011.
 Barbu, V., Differential Equations, Springer Undergraduate Mathematics Series, Springer, Cham, 2016. Translated from the 1985 Romanian original by Liviu Nicolaescu.
 Barro, R.J., Sala-i Martin, X.,Economic Growth, MIT Press, 2004.
 Bueno de Mesquita, B., Lalman, D.,Arms races and the opportunity for peace, Formal Analysis in International Relations,76(1988), 263-283.
 Dawkins, R., Krebs, J.R., Arms races between and within species, Proceedings of the Royal Society of London, Series B, Biological Sciences,205(1161)(1979), 489-511.
 Hill, W.W., Several sequential augmentations of Richardson’s arms race model, Math.
 Intriljgator, M.D., Brito, D.L.,Arms races, Defence and Peace Economics,11(2000), no.
 Krabs, W.,Dynamische Systeme: Steuerbarkeit und Chaotisches Verhalten, B. G. Teub- ner, Stuttgart, 1998.
 Larrosa, J.M.C.,Arms build-up and arms race in optimal economic growth, Int. J. Econ.
Theory,12(2016), no. 2, 167-182.
 Majeski, S.J., Arms races as iterated prisoner’s dilemma games, Math. Social Sci., 7(3)(1984), 253-266.
 Mankiw, N.G., Romer, D., Weil, D.N., A contribution to the empirics of economic growth, The Quarterly Journal of Economics,107(1992), no. 2, 407-437.
 Richardson, L.F.,Arms and Insecurity: A Mathematical Study of the Causes and Ori- gins of War, Edited by Nicolas Rashevsky and Ernesto Trucco, The Boxwood Press, Pittsburgh, Pa.; Quadrangle Books, Chicago, Ill., 1960.
 Richardson, L.F.,Statistics of Deadly Quarrels, Boxwood Press, 1960.
 Solow, R.M.,A contribution to the theory of economic growth, The Quarterly Journal of Economics,70(1956), no. 1, 65-94.
 Swan, T.W., Economic growth and capital accumulation, Economic Record,32(1956), no. 2, 334-361.
NTT DATA Romania
19-21, Constant¸a Street, 400158 Cluj-Napoca, Romania e-mail:[email protected]
Faculty of Mathematics and Computer Sciences 1, Kog˘alniceanu Street, 400084 Cluj-Napoca, Romania e-mail:[email protected]