# American Institute of Mathematical Sciences

January  2021, 8(1): 59-97. doi: 10.3934/jcd.2021004

## A general framework for validated continuation of periodic orbits in systems of polynomial ODEs

 1 VU Amsterdam, Department of Mathematics, De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 2 Rutgers University, Department of Mathematics, 110 Frelinghusen Rd, Piscataway, NJ 08854-8019, USA

Received  February 2019 Revised  November 2019 Published  October 2020

Fund Project: The authors are partially supported by NWO-VICI grant 639033109

In this paper a parametrized Newton-Kantorovich approach is applied to continuation of periodic orbits in arbitrary polynomial vector fields. This allows us to rigorously validate numerically computed branches of periodic solutions. We derive the estimates in full generality and present sample continuation proofs obtained using an implementation in Matlab. The presented approach is applicable to any polynomial vector field of any order and dimension. A variety of examples is presented to illustrate the efficacy of the method.

Citation: Jan Bouwe van den Berg, Elena Queirolo. A general framework for validated continuation of periodic orbits in systems of polynomial ODEs. Journal of Computational Dynamics, 2021, 8 (1) : 59-97. doi: 10.3934/jcd.2021004
##### References:

show all references

##### References:
Representation of the matrix introduced in Equation (25)
The shape of the operator $Q$, as defined in Lemma 5.1. The vertical columns represent the vectors $y^{[k]}$, where $k$ denotes the column index
The boundaries of the validation interval $(r_{\min},r_{\max})$ versus the number of modes $K$ for various values of $\nu$
On the right: the minimal number of modes $K$ necessary for the validation of a periodic orbit of (52) depending on $\mu$. The stepsize in $\mu$ is not chosen uniformly, because the problem is more sensitive to changes in $\mu$ for $\mu$ bigger, therefore requesting a steeper change in the number of modes. One the left: several of the validated periodic orbits
The computation time versus the number of modes $K$. The reference line has a slope corresponding to $K^3$
The computation time versus the dimension $N$ of the system (53). The reference lines have slopes corresponding to quadratic and cubic dependence on $N$
The computation time versus the order $D$ of the system (54). The reference lines have slopes corresponding to quadratic and cubic dependence on $D$
Validation of a periodic solution for the Lorenz system (55) with the classical parameter values $\sigma = 10$, $\beta = 8/3$ and $\rho = 28$. The solution has a period close to 25. This solution was validated with $\nu = 1+10^{-6}$ and $K = 800$
Step size versus $\mu$. Some oscillation is noticeable in the step size. This is due to the fact that the stepsize is decreased when the number of modes is changed, see Section 9. During this validation, the number of modes has been increased automatically from 4 to about 400
Validated continuation for the Lorenz system (55) from $\rho = 28$ to $\rho = 14.8$. In the left graph, the lower orbit still has two distinct swirls in the left lobe, but they are too close to each other to be seen in this graph
Validated continuation for $10$ coupled Lorenz systems (77) from $\epsilon = 0.14205$ to $\epsilon = 0.14252$. In the left graph we have depicted in blue the orbit of the first three component, in green the orbit of the last three component, with the biggest validated $\epsilon$, in the right graph we depicted the values of $\epsilon$, where we can notice the adaptation of the stepsize
Continuation through the fold of the Rychkov system (78). The validation started at the top, where a small initial stepsize was chosen. One can see how the stepsize was automatically adjusted along the validated curve. As can be seen from the values along the horizontal axis, we only depict a validated continuation for parameter values very close to saddle-node
 [1] Alain Jacquemard, Weber Flávio Pereira. On periodic orbits of polynomial relay systems. Discrete & Continuous Dynamical Systems, 2007, 17 (2) : 331-347. doi: 10.3934/dcds.2007.17.331 [2] J. P. Lessard, J. D. Mireles James. A functional analytic approach to validated numerics for eigenvalues of delay equations. Journal of Computational Dynamics, 2020, 7 (1) : 123-158. doi: 10.3934/jcd.2020005 [3] Jan Bouwe van den Berg, Ray Sheombarsing. Rigorous numerics for ODEs using Chebyshev series and domain decomposition. Journal of Computational Dynamics, 2021, 8 (3) : 353-401. doi: 10.3934/jcd.2021015 [4] Paolo Perfetti. An infinite-dimensional extension of a Poincaré's result concerning the continuation of periodic orbits. Discrete & Continuous Dynamical Systems, 1997, 3 (3) : 401-418. doi: 10.3934/dcds.1997.3.401 [5] Tiago de Carvalho, Rodrigo Donizete Euzébio, Jaume Llibre, Durval José Tonon. Detecting periodic orbits in some 3D chaotic quadratic polynomial differential systems. Discrete & Continuous Dynamical Systems - B, 2016, 21 (1) : 1-11. doi: 10.3934/dcdsb.2016.21.1 [6] Maxime Breden, Jean-Philippe Lessard. Polynomial interpolation and a priori bootstrap for computer-assisted proofs in nonlinear ODEs. Discrete & Continuous Dynamical Systems - B, 2018, 23 (7) : 2825-2858. doi: 10.3934/dcdsb.2018164 [7] Ana Cristina Mereu, Marco Antonio Teixeira. Reversibility and branching of periodic orbits. Discrete & Continuous Dynamical Systems, 2013, 33 (3) : 1177-1199. doi: 10.3934/dcds.2013.33.1177 [8] Ilie Ugarcovici. On hyperbolic measures and periodic orbits. Discrete & Continuous Dynamical Systems, 2006, 16 (2) : 505-512. doi: 10.3934/dcds.2006.16.505 [9] Katrin Gelfert, Christian Wolf. On the distribution of periodic orbits. Discrete & Continuous Dynamical Systems, 2010, 26 (3) : 949-966. doi: 10.3934/dcds.2010.26.949 [10] Jacky Cresson, Christophe Guillet. Periodic orbits and Arnold diffusion. Discrete & Continuous Dynamical Systems, 2003, 9 (2) : 451-470. doi: 10.3934/dcds.2003.9.451 [11] Annalisa Iuorio, Christian Kuehn, Peter Szmolyan. Geometry and numerical continuation of multiscale orbits in a nonconvex variational problem. Discrete & Continuous Dynamical Systems - S, 2020, 13 (4) : 1269-1290. doi: 10.3934/dcdss.2020073 [12] Peter Albers, Jean Gutt, Doris Hein. Periodic Reeb orbits on prequantization bundles. Journal of Modern Dynamics, 2018, 12: 123-150. doi: 10.3934/jmd.2018005 [13] Peter Giesl, Holger Wendland. Approximating the basin of attraction of time-periodic ODEs by meshless collocation. Discrete & Continuous Dynamical Systems, 2009, 25 (4) : 1249-1274. doi: 10.3934/dcds.2009.25.1249 [14] Chiara Zanini, Fabio Zanolin. Periodic solutions for a class of second order ODEs with a Nagumo cubic type nonlinearity. Discrete & Continuous Dynamical Systems, 2012, 32 (11) : 4045-4067. doi: 10.3934/dcds.2012.32.4045 [15] László Mérai, Igor E. Shparlinski. Unlikely intersections over finite fields: Polynomial orbits in small subgroups. Discrete & Continuous Dynamical Systems, 2020, 40 (2) : 1065-1073. doi: 10.3934/dcds.2020070 [16] Qiaoyi Hu, Zhijun Qiao. Analyticity, Gevrey regularity and unique continuation for an integrable multi-component peakon system with an arbitrary polynomial function. Discrete & Continuous Dynamical Systems, 2016, 36 (12) : 6975-7000. doi: 10.3934/dcds.2016103 [17] Amadeu Delshams, Pere Gutiérrez. Exponentially small splitting for whiskered tori in Hamiltonian systems: continuation of transverse homoclinic orbits. Discrete & Continuous Dynamical Systems, 2004, 11 (4) : 757-783. doi: 10.3934/dcds.2004.11.757 [18] Jan Sieber, Matthias Wolfrum, Mark Lichtner, Serhiy Yanchuk. On the stability of periodic orbits in delay equations with large delay. Discrete & Continuous Dynamical Systems, 2013, 33 (7) : 3109-3134. doi: 10.3934/dcds.2013.33.3109 [19] Nicola Guglielmi, Christian Lubich. Numerical periodic orbits of neutral delay differential equations. Discrete & Continuous Dynamical Systems, 2005, 13 (4) : 1057-1067. doi: 10.3934/dcds.2005.13.1057 [20] Jorge Rebaza. Bifurcations and periodic orbits in variable population interactions. Communications on Pure & Applied Analysis, 2013, 12 (6) : 2997-3012. doi: 10.3934/cpaa.2013.12.2997

Impact Factor: