arxiv: 2605.13823 · v1 · submitted 2026-05-13 · 🧮 math.DS

Recognition: 2 theorem links

· Lean Theorem

Stability analysis of Richardson models with delay for confrontation between two countries

Anatoliy A.Martynyuk, Teresa Faria

Pith reviewed 2026-05-14 17:28 UTC · model grok-4.3

classification 🧮 math.DS

keywords delay differential equationsstability analysisRichardson modelHopf bifurcationglobal asymptotic stabilitynon-autonomous systemsinternational conflict models

0 comments

The pith

A delay differential equation model establishes global asymptotic stability for two-country confrontations with time-varying hostility.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper develops a non-autonomous mathematical model based on delay differential equations to capture the confrontation dynamics between two countries, incorporating a hostility factor. It provides comprehensive stability analysis for the autonomous case, including conditions for asymptotic stability and descriptions of Hopf bifurcations at critical delay values. For the non-autonomous version, it derives conditions that guarantee global asymptotic stability for both the linear approximation and the full nonlinear system. A sympathetic reader would care because these results offer a framework to predict when such confrontations might settle into stable equilibria despite changing conditions. The framework of special solutions for delay differential equations is used to support these findings.

Core claim

The central claim is that the proposed model with delay has an equilibrium state whose stability can be ensured by conditions involving the hostility factor, with global asymptotic stability holding for the non-autonomous system under appropriate bounds on the time-varying coefficients.

What carries the argument

The key machinery is the framework of special solutions for delay differential equations, which is applied to derive stability criteria for both autonomous and non-autonomous cases.

If this is right

If the stability conditions are met, the equilibrium remains globally attractive for all initial functions in the nonlinear model.
Hopf bifurcations occur when the delay parameter crosses critical values, leading to periodic solutions in the autonomous case.
Time-dependent hostility coefficients do not prevent global stability provided they satisfy the established inequalities.
The model extends classical Richardson models by including delay and non-autonomous terms while preserving stability properties.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

These stability results could inform models of other bilateral competitions, such as economic or technological rivalries between nations.
Future work might incorporate stochastic perturbations to test robustness against random events in international relations.
Empirical validation could involve fitting the model to historical arms expenditure data from past conflicts.

Load-bearing premise

The dynamics of country confrontation can be adequately represented by a deterministic delay differential equation system with coefficients that are either constant or explicitly time-dependent.

What would settle it

A counterexample would be a set of initial conditions and parameter values satisfying the derived stability conditions but where numerical simulations of the model show that solutions diverge or oscillate indefinitely away from the equilibrium.

read the original abstract

This article proposes a non-autonomous mathematical model with delay for confrontation between two countries, and examines the stability of its equilibrium state. Our criteria for stability take into account the influence of the factor of hostility between countries. For the autonomous case, the asymptotic stability is studied in a comprehensive way, and the Hopf bifurcations occurring as the delay crosses some critical values are described. For the non-autonomous model, conditions ensuring the global asymptotic stability for both the linear approximation and the nonlinear system are established. The framework of special solutions for delay differential equations is also applied.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper derives global stability conditions for a non-autonomous Richardson model with delays and time-varying hostility using standard DDE tools, extending the usual autonomous analysis without major new techniques.

read the letter

This paper takes the classic Richardson arms-race equations, adds constant and time-varying delays plus a hostility factor, and works out stability for both autonomous and non-autonomous versions. What it does cleanly is establish global asymptotic stability conditions for the non-autonomous linear and nonlinear systems via the special-solutions framework and Lyapunov-type arguments. The autonomous part also tracks Hopf bifurcations as delay crosses critical values. That gives explicit thresholds that could be checked numerically in similar models. The approach follows established DDE theory without circularity or hidden fitting, and the stress-test found no internal contradictions in the logical steps. The main soft spot is the modeling assumption that deterministic DDEs with explicitly given coefficients capture the confrontation dynamics well enough; real cases often involve stochastic elements or multi-agent feedback that are left out. Without seeing the full proof details it is hard to judge how tight the bounds end up being, but nothing in the stated results looks broken. This is mainly for people who build or analyze mathematical models of social conflict or who work on stability in non-autonomous delay systems. A reader wanting concrete conditions for comparable extensions would get some value. It is not something I would cite in my own work, but the claims are grounded enough that it deserves a serious referee rather than a desk reject.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a non-autonomous delay differential equation model for confrontation dynamics between two countries that incorporates a time-varying hostility factor. For the autonomous case it derives conditions for local asymptotic stability of the equilibrium and characterizes Hopf bifurcations at critical delay values. For the non-autonomous case it establishes global asymptotic stability of the equilibrium for both the linearized system and the full nonlinear system, employing the special-solutions framework together with Lyapunov-type or comparison arguments.

Significance. If the derivations are correct, the work supplies explicit, verifiable stability criteria for a deterministic DDE model of bilateral conflict that includes both constant and explicitly time-dependent coefficients. The global-stability results for the non-autonomous system constitute the strongest contribution, as they avoid stochastic or multi-agent extensions while still handling time-varying hostility.

major comments (2)

[§3] §3 (autonomous case): the characteristic equation and the transversality condition for the Hopf bifurcation are stated but the algebraic verification that the real part of the eigenvalue crosses the imaginary axis with nonzero speed is only sketched; an explicit computation of d(Re λ)/dτ at the critical value is required to confirm the bifurcation is supercritical or subcritical as claimed.
[§4] §4 (non-autonomous case): the global asymptotic stability proof for the nonlinear system relies on a comparison principle applied to a special solution; the argument assumes the hostility factor remains bounded, yet the manuscript does not state an explicit uniform bound or show that the comparison function remains positive for all initial data, which is load-bearing for the global claim.

minor comments (2)

[§2] The model equations (1)–(2) should be displayed with all parameters labeled (including the hostility coefficient) immediately after the introduction for easy reference.
[§5] Figure 1 (phase portraits) lacks axis labels and a caption indicating the parameter values used; this reduces readability of the numerical illustrations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate the revisions we will make to strengthen the proofs.

read point-by-point responses

Referee: [§3] §3 (autonomous case): the characteristic equation and the transversality condition for the Hopf bifurcation are stated but the algebraic verification that the real part of the eigenvalue crosses the imaginary axis with nonzero speed is only sketched; an explicit computation of d(Re λ)/dτ at the critical value is required to confirm the bifurcation is supercritical or subcritical as claimed.

Authors: We agree that the transversality condition requires an explicit algebraic verification rather than a sketch. In the revised manuscript we will compute d(Re λ)/dτ at the critical delay τ₀ by differentiating the characteristic equation with respect to τ, substituting the purely imaginary root, and simplifying to obtain the sign of the real-part derivative. This will rigorously confirm the crossing direction and allow us to state whether the Hopf bifurcation is supercritical or subcritical. revision: yes
Referee: [§4] §4 (non-autonomous case): the global asymptotic stability proof for the nonlinear system relies on a comparison principle applied to a special solution; the argument assumes the hostility factor remains bounded, yet the manuscript does not state an explicit uniform bound or show that the comparison function remains positive for all initial data, which is load-bearing for the global claim.

Authors: We acknowledge the need for explicit statements. The hostility factor is continuous and positive by model construction (Section 2, equation (1)); we will add an explicit uniform bound derived from the parameter assumptions and prove that the comparison function remains strictly positive for all admissible initial data. These clarifications will be inserted into the revised Section 4 to make the global-stability argument fully rigorous. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper starts from an explicitly stated system of non-autonomous delay differential equations modeling bilateral confrontation, then applies standard Lyapunov-type arguments, comparison principles, and the special-solutions framework to derive explicit stability criteria. No parameters are fitted to data, no quantity is renamed as a prediction after being used as an input, and no load-bearing step reduces to a self-citation whose validity is presupposed by the present work. The autonomous Hopf-bifurcation analysis follows directly from the characteristic equation of the linearized system, while the global-stability results for the non-autonomous case are obtained by bounding the time-varying hostility coefficient within the same model equations. The derivation chain is therefore self-contained and independent of its own outputs.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claims rest on the assumption that the confrontation process is faithfully represented by a deterministic delay differential system whose interaction terms incorporate a scalar hostility factor; no independent empirical validation of this modeling choice is provided.

free parameters (2)

hostility factor
Scalar parameter scaling the strength of mutual antagonism between the two countries; its specific functional form and range are not detailed in the abstract.
delay tau
Fixed time lag in each country's response to the other's actions; treated as a bifurcation parameter.

axioms (1)

domain assumption The dynamics of bilateral confrontation admit a deterministic representation by a system of delay differential equations with linear and possibly nonlinear interaction terms.
Invoked by the proposal of the model itself.

pith-pipeline@v0.9.0 · 5385 in / 1288 out tokens · 50067 ms · 2026-05-14T17:28:32.835352+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
Theorem 3.2 … global asymptotic stability via comparison functions G,H and lim sup Bi(t)<1

Reference graph

Works this paper leans on

27 extracted references · 1 canonical work pages

[1]

Arino O., Gy¨ ori I., Pituk M., Asymptotically diagonal delay differential systems,J. Math. Anal. Appl.204 (1996), 701–728

1996
[2]

Berezansky L., Dibl´ ık J., Svoboda Z., Smarda Z., Exponential stability of linear delayed differ- ential systems,Appl. Math. Comput.320 (2018), 474–484

2018
[3]

Bohner M., Martynyuk A.A., Equilibrium stability under nuclear confrontation,Differ. Equ. Dyn. Syst.34 (2024), 209–223

2024
[4]

Burton T.A.,Stability and Periodic Solutions of Ordinary and Functional Differential Equations, 2nd Ed., Dover Publications Inc., Mineola, NY, 2005

2005
[5]

Studies Quarterly11 (1967), 63–88

Caspary W.R., Richardson’s model of arms rates: description, critique, and an alternative model,Int. Studies Quarterly11 (1967), 63–88

1967
[6]

Chow S.-N., Hale J.K.,Methods of Bifurcation Theory, Springer-Verlag, New York, 1982

1982
[7]

629, Springer, Berlin, 1978

Coppel W.A.,Dichotomies in Stability Theory, Lecture Notes in Mathematics, Vol. 629, Springer, Berlin, 1978

1978
[8]

Driver R.D., On Ryabov’s asymptotic characterization of the solutions of quasi-linear differential equations with small delays,SIAM Review10 (1968), 329–341

1968
[9]

Differential Equations21 (1976), 148–166

Driver R.D., Linear differential systems with small delays,J. Differential Equations21 (1976), 148–166

1976
[10]

Faria T., Stability for non-autonomous linear differential systems with infinite delay,J. Dyn. Differ. Equ.34 (2022), 747–773

2022
[11]

Differential Equations122 (1995), 181-200

Faria T., Magalh˜ aes L.T., Normal forms for retarded functional differential equations with parameters and applications to Hopf singularity,J. Differential Equations122 (1995), 181-200

1995
[12]

Faria T., Oliveira J.J., Global stability for functional differential equations with infinite delay and applications to Hopfield-type neural networks,J. Dyn. Differ. Equ.(2026). https://doi.org/10.1007/s10884-025-10468-w

work page doi:10.1007/s10884-025-10468-w 2026
[13]

Kluwer, Dordrechit, 1986

Fiedler M.,Special Matrices and Their Applications in Numerical Mathematics, Martinus Nijhoff Publ. Kluwer, Dordrechit, 1986

1986
[14]

Thesis, Department of Mathematics, Jimma University, 2017

Guye E.,Mathematical Modeling and Stability Analysis of Conflict among Three Nations, MSc. Thesis, Department of Mathematics, Jimma University, 2017

2017
[15]

Gy¨ ori I., Pituk M., Stability criteria for linear delay differential equations,Differential Integral Equations10 (1997), 841–852

1997
[16]

Hale J.K.,Theory of Functional Differential Equations, 2nd Ed., Applied Mathematical Sciences, Springer-Verlag, New York-Heidelberg, 1977. 25

1977
[17]

Hale J.K., Lunel S.M.V.,Introduction to FunctionalDifferential Equations, Applied Mathemat- ical Sciences, Springer-Verlag, New York, 1993

1993
[18]

Hill W.W., A time lagged arms race model,J.Peace Sci.3 (1978), 55–62

1978
[19]

Kuang Y.,Delay Differential Equations with Applications in Population Dynamics, Mathematics in Science and Engineering, Academic Press Inc., Boston, 1993

1993
[20]

Lyapunov A.M.,The general problem of the stability of motion(English title), original: Kharkov, 1892, reprinted by ONTI, 1935; Translator: A. T. Fuller. Translation edition: Int. J. Control, Vol. 55, 1992

1935
[21]

A.,Integral Methods in Nonlinear Dynamics of Systems, Ser

Martynyuk A. A.,Integral Methods in Nonlinear Dynamics of Systems, Ser. Adv. Math. Appl. Sci. 96, World Scientific Publishing Co., 2026

2026
[22]

Martynyuk A.A, Stamov G., Stamova I., Novel estimates of the impact of hostility on the equilibrium stability in an arms race model,Electron. J. Qual. Theory Differ. Equ. 2025, Paper No. 58, 1–19

2025
[23]

Rashevsky and E

Richardson L.F.,Arms and Insecurity: A mathematical study of the causes and origins of war, N. Rashevsky and E. Trucco Eds., The Boxwood Press, Chicago, 1960

1960
[24]

Ryabov Y.A., Certain asymptotic properties of linear systems with small time-lag,Trudy Sem. Teor. Differencial. Uravnenii s Otklon. Argumentom Univ. Druzby Narodov Patrisa Lumumby 3 (1965), 153–164 (in Russian)

1965
[25]

In: Gillespie J.V., Zinnes D.A

Siljak D.D., On the stability of arms race. In: Gillespie J.V., Zinnes D.A. (eds.) Mathematical Systems in International Relations Research, New York: Praeger Publishers, 1977

1977
[26]

In: Gleditsch, N.P

Smith R.P., The Influence of the Richardson Arms Race Model. In: Gleditsch, N.P. (eds.) Lewis Fry Richardson: His Intellectual Legacy and Influence in the Social Sciences. Pioneers in Arts, Humanities, Science, Engineering, Practice, vol 27. Springer, Cham. (2020)

2020
[27]

Wallace M.D., Arms races and the balance of power: a mathematical model,Appl. Math. Modelling, 1 (1976), 83–92. 26

1976