[Paper Review] Geometric analysis of Oscillations in the Frzilator model

Hadi Taghvafard, Hildeberto Jardón-Kojakhmetov+2 more|arXiv (Cornell University)|Dec 2, 2019

Nonlinear Dynamics and Pattern Formation29 references5 citations

TL;DR

This paper proves the existence of a strongly attracting limit cycle in the Frzilator model—a biochemical oscillator governing myxobacteria development—using geometric singular perturbation theory and blow-up methods. The limit cycle arises as a relaxation oscillation due to small Michaelis-Menten constants, with a well-defined structure across multiple timescales, confirmed numerically and analytically for a range of parameters including γ ∈ (0.0561, 0.1177).

ABSTRACT

A biochemical oscillator model, describing developmental stage of myxobacteria, is analyzed mathematically. Observations from numerical simulations show that in a certain range of parameters, the corresponding system of ordinary differential equations displays stable and robust oscillations. In this work, we use geometric singular perturbation theory and blow-up method to prove the existence of a strongly attracting limit cycle. This cycle corresponds to a relaxation oscillation of an auxiliary system, whose singular perturbation nature originates from the small Michaelis-Menten constants of the biochemical model. In addition, we give a detailed description of the structure of the limit cycle, and the timescales along it.

Motivation & Objective

• To rigorously establish the existence of a strongly attracting limit cycle in the Frzilator model, a biochemical oscillator governing myxobacteria development.
• To analyze the oscillatory dynamics of the Frzilator model through a two-timescale (slow-fast) system derived from rescaling the original ODEs.
• To characterize the geometric structure and timescale decomposition of the limit cycle, particularly focusing on transitions across non-hyperbolic manifolds.
• To identify parameter ranges—specifically γ ∈ (0.0561, 0.1177)—where the qualitative dynamics of the fast subsystem remain consistent, ensuring robust oscillations.
• To extend analytical techniques from the Goldbeter minimal model to a more complex, biologically relevant system with non-hyperbolic structures.

Proposed method

• Rescale the original Frzilator ODE system to reveal a singularly perturbed slow-fast structure, with ε representing the ratio of timescales.
• Apply geometric singular perturbation theory (GSPT) to analyze the layer and reduced problems, identifying critical manifolds and their stability.
• Use the blow-up method to resolve non-hyperbolic dynamics at fold points and transition layers, particularly where the fast dynamics cross critical surfaces.
• Construct the full limit cycle by stitching together segments of slow flow on attracting critical manifolds and fast jumps between them.
• Numerically solve transcendental equations arising from boundary conditions on the critical manifolds to determine parameter ranges where the cycle structure is preserved.
• Verify the robustness of the limit cycle structure through numerical simulations and parameter continuation, confirming convergence of almost all trajectories to the cycle.

Experimental results

Research questions

• RQ1 Does the Frzilator model exhibit a unique, strongly attracting limit cycle for biologically relevant parameter values?
• RQ2 How does the slow-fast structure of the Frzilator model give rise to relaxation oscillations with distinct timescales?
• RQ3 What is the geometric and dynamical structure of the limit cycle, particularly in regions with non-hyperbolic behavior?
• RQ4 For which parameter ranges is the qualitative behavior of the fast dynamics preserved, ensuring the existence of a stable limit cycle?
• RQ5 Can the blow-up method effectively resolve the dynamics near non-hyperbolic lines where standard GSPT fails?

Key findings

• A strongly attracting limit cycle exists in the Frzilator model for a range of parameters, confirmed via geometric singular perturbation theory and blow-up analysis.
• The limit cycle is a relaxation oscillation, with distinct fast and slow phases, arising from the singular perturbation structure induced by small Michaelis-Menten constants.
• For γ ∈ (0.0561, 0.1177), the fast dynamics consistently transition from initial points in 𝒮₀,₃ᵃ to 𝒮₀,₆ᵃ and from 𝒮₀,₆ᵃ to 𝒮₀,₁ᵃ, preserving the cycle structure.
• The blow-up method successfully resolves non-hyperbolic dynamics at fold points, enabling a complete description of the limit cycle’s geometry.
• Numerical simulations confirm that almost all trajectories converge to the same limit cycle, indicating robust oscillatory behavior.
• The analysis provides a detailed decomposition of the cycle into slow segments on critical manifolds and fast jumps across transition layers, with explicit timescale separation.