Numerical Periodic Normalization at Codim 1 Bifurcations of Limit Cycles in DDEs
M. M. Bosschaert, B. Lentjes, L. Spek, Yu. A. Kuznetsov
TL;DR
This work derives explicit, computable formulas for the critical normal form coefficients of codimension-one bifurcations of limit cycles in delay differential equations by combining periodic center manifolds with sun-star calculus and a novel characteristic operator Δ(z). The approach yields compact expressions for fold, period-doubling, and Neimark-Sacker bifurcation coefficients, enabling robust numerical evaluation via orthogonal collocation on periodic BVPs and bordered inverses, implemented in the Julia package PeriodicNormalizationDDEs. It avoids Poincaré maps, works directly with the center manifold reductions, and provides a unifying framework that mirrors finite-dimensional ODE results while addressing periodic spectral challenges through the adjoint characterstic operator. Applications to neural-field and active-control models demonstrate correctness, efficiency gains over pseudospectral methods, and practical utility for bifurcation analysis of limit cycles in DDEs. The methodology paves the way for codimension-two extensions and broader applicability to periodic evolution equations beyond classical DDEs.
Abstract
Recent work in [53, 54] by the authors on periodic center manifolds and normal forms for bifurcations of limit cycles in delay differential equations (DDEs) motivates the derivation of explicit computational formulas for the critical normal form coefficients of all codimension one bifurcations of limit cycles. In this paper, we derive such formulas via an application of the periodic normalization method in combination with the functional analytic perturbation framework for dual semigroups (sun-star calculus). The explicit formulas allow us to distinguish between nondegenerate, sub- and supercritical bifurcations. To efficiently apply these formulas, we introduce the characteristic operator as this enables us to use robust numerical boundary-value algorithms based on orthogonal collocation. Although our theoretical results are proven in a more general setting, the software implementation and examples focus on discrete DDEs. The actual implementation is described in detail and its effectiveness is demonstrated on various models.
