🔢 mathematics

An intrinsic characterization of the Bogdanov-Takens normal-form coefficients and a mixed-volume obstruction to non-isolated degeneracies

This paper provides a coordinate-free characterization of Bogdanov-Takens normal-form coefficients as directional derivatives of Jacobian invariants and establishes a combinatorial obstruction based on mixed volume that limits the types of nilpotent singularities achievable in Kolmogorov systems.

Original authors: Víctor Castellanos, Ramón Eduardo Chan-López

Published 2026-08-17

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Víctor Castellanos, Ramón Eduardo Chan-López

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, invisible dance floor where everything from planets to bacteria is constantly moving. In the world of mathematics, scientists use "vector fields" to map out the steps of this dance. These maps show where things are going and how fast. Sometimes, dancers get stuck in one spot, called an "equilibrium." Usually, if you nudge a dancer slightly, they wobble and find a new rhythm. But sometimes, the dance floor has a weird glitch where a nudge doesn't just change the rhythm; it completely reshapes the dance. This is where things get tricky. Mathematicians study these "bifurcations"—moments where the rules of the dance suddenly change. One of the most famous glitches is called the Bogdanov–Takens bifurcation. It's like a double-jointed move where the system is so balanced on a knife-edge that two different ways of breaking the balance can happen at once. Understanding exactly how and when this happens is crucial for predicting everything from how animal populations boom and crash to how chemical reactions explode or fizzle out.

Now, imagine you are trying to describe this glitchy dance move to a friend. Usually, you'd have to use a complicated, coordinate-heavy language, like giving directions using a specific grid map that only works if you stand in a certain spot. If you move the map, the instructions get messy. This new paper by V. Castellanos and E. Chan-López is like discovering a secret, universal language for this dance. They found a way to describe the "glitch" using only the intrinsic properties of the dance floor itself, without needing any external maps or complicated extra steps. They proved that the two key numbers that tell us if the glitch is "normal" or "broken" are actually just simple measurements of how the dance floor's slope changes in a specific direction. It's like realizing that instead of measuring the wind with a complex weather station, you can just feel the breeze on your cheek to know exactly how strong it is.

The authors also discovered a fascinating "combinatorial" rule, which is like a game of Tetris for math. They looked at a specific type of system used to model nature (like predators eating prey) and found that the shape of the mathematical "blocks" (called Newton polytopes) that make up the system acts as a strict bouncer. If the blocks are too small or simple, the dance floor simply cannot support certain types of weird glitches. Specifically, they proved that if the "mixed volume" (a way of measuring how much space these blocks take up together) is small, the system is forced to have an isolated equilibrium. It's as if the rules of the game say, "If your building blocks are this small, you cannot build a twisted, multi-level maze where the dancers get stuck in a whole line; you are forced to build a single, isolated point."

In the first part of their work, the authors tackle the problem of how to calculate the "coefficients" (the numbers $a$ and $b$ ) that define the Bogdanov–Takens glitch. Traditionally, finding these numbers is like trying to solve a puzzle by first building a giant, unnecessary scaffolding around it. You have to find special "generalized eigenvectors" and do complex multi-step calculations. The authors show that this scaffolding is unnecessary. They prove that these two numbers are simply the "directional derivatives" of two basic properties of the system: the trace and the determinant of the Jacobian matrix (which are like the system's total energy and its tendency to expand or contract). Imagine you are standing on a hill (the equilibrium point). The "trace" and "determinant" are like the height and the slope of the hill. The authors show that the numbers $a$ and $b$ are just how fast the height and slope change if you take a single step in the direction of the "kernel" (the specific direction where the system is stuck). This is a huge simplification. It means you don't need to build the whole scaffolding; you just need to look at the hill and take one step. They also show that this method is "invariant," meaning it works no matter how you rotate or stretch your coordinate system, making it a much more robust tool for scientists.

The second part of the paper is a "mixed-volume obstruction." This sounds scary, but it's actually a very elegant rule about what is impossible in certain systems. The authors focus on "Kolmogorov systems," which are mathematical models often used in ecology to describe how species interact (like lions and zebras). These systems are built from polynomials, which are just sums of terms like $x$ , $y$ , $xy$, $x^2$ , etc. Each polynomial has a "Newton polytope," which is a geometric shape formed by connecting the dots representing the powers of the variables. The authors prove that the "order" of the glitch (how complex the dance move is) is limited by the "mixed volume" of these shapes.

Here is the punchline: If the mixed volume of the shapes is 2 or less, the system must have an isolated equilibrium where the glitch coefficient $a$ is not zero. In plain English, if the mathematical building blocks are simple enough, the system is forced to have a "cusp" type of glitch (a sharp, pointy turn) rather than a "saddle" or "focus" type (which are more complex and unstable). The paper explicitly rules out the possibility of finding complex, non-isolated glitches in systems with small mixed volumes. It's a hard limit, like a law of physics. If you try to build a system with these simple blocks and expect a complex, non-isolated glitch (where the dancers get stuck in a line), the math says, "Nope, that's impossible." The only way for the coefficient $a$ to vanish (which would allow for a non-isolated line of equilibria) is if the building blocks are complex enough to support it.

The authors tested this theory on two famous ecological models: one involving "predator satiation" (where predators get full and stop eating) and another involving "mutual interference" (where predators get in each other's way). In both cases, they found that the mixed volume was exactly 2. As predicted by their new rule, the coefficient $a$ never vanished (meaning the equilibrium was always isolated), and the only time the system became "degenerate" (meaning the glitch got weird) was when the second coefficient, $b$ , became zero. This confirmed that the systems were stuck in the "cusp" category and could not reach the more complex "saddle" or "focus" categories.

The paper is rigorous and mathematically proven, not just a suggestion or a simulation. The authors provide exact formulas and logical proofs that hold true for these classes of systems. They don't just say "it looks like this"; they show that the geometry of the equations forces this outcome. The only thing they don't determine is the exact "codimension" (a measure of how rare the glitch is) beyond the fact that it's at least three. They leave that door slightly open, noting that while the mixed volume controls the first part of the glitch, it doesn't tell the whole story about the second part. But for the main question of whether these complex glitches can exist in simple systems, the answer is a definitive "no."

In summary, this paper gives us two powerful new tools. First, a simple, coordinate-free way to measure the "glitchiness" of a system without doing messy calculations. Second, a combinatorial rule that acts as a gatekeeper, telling us that certain types of complex behaviors are simply impossible in systems with simple mathematical structures. It's a beautiful example of how deep mathematical structures can reveal hidden limits in the natural world, showing us that sometimes, the shape of the equation itself decides the fate of the dance.

Technical Summary: Intrinsic Characterization of Bogdanov–Takens Coefficients and Mixed-Volume Obstructions

Problem Statement
The Bogdanov–Takens (BT) bifurcation is the generic codimension-two local bifurcation of a planar equilibrium with a nilpotent linear part of rank one. Its normal form on the center manifold is given by $\dot{u} = v, \dot{v} = a u^2 + b u v + o(|(u, v)|^2)$ . The nondegeneracy of the bifurcation requires $ab \neq 0$ . The paper addresses two specific challenges:

Intrinsic Computation: The standard calculation of the coefficients $a$ and $b$ typically relies on generalized eigenvectors, biorthonormalization, and second-order multilinear forms, which are coordinate-dependent and computationally cumbersome.
Structural Obstructions: For a specific class of Kolmogorov systems (predator–prey models), it is unclear whether certain degeneracies (specifically $a=0$ ) can occur at isolated equilibria, or if they are structurally forbidden by the algebraic form of the vector field.

Methodology
The paper employs a combination of differential geometry, singularity theory, and algebraic combinatorics (specifically Newton polytopes and mixed volumes).

Intrinsic Derivatives: The authors derive an identity expressing the BT coefficients $a$ and $b$ directly as directional derivatives of the Jacobian's invariants (determinant and trace) along the kernel of the linearization. This avoids the construction of normal forms or generalized eigenvectors.
Geometric Interpretation: The equilibrium manifold of a parameterized family is analyzed. The authors relate the vanishing of $a$ and $b$ to the tangency of the kernel line field with the "fold set" (where the determinant vanishes) and the "neutrality set" (where the trace vanishes).
Combinatorial Bounds: For Kolmogorov systems of the form $\dot{x} = x A/g_1, \dot{y} = y B/g_2$ , the authors utilize the Bernstein–Kushnirenko theorem. They relate the local multiplicity of an isolated equilibrium to the mixed volume ($MV$) of the Newton polytopes of the polynomials $A$ and $B$ . The multiplicity of a nilpotent zero is shown to equal the order of the function $f$ in the Takens normal form ( $\text{ord } f$ ).

Key Contributions and Results

1. Intrinsic Characterization of Coefficients (Theorem A)
The paper establishes that for a planar vector field $X$ with a nilpotent rank-one Jacobian $J = DX(p)$ and kernel spanned by $q_0$ , the BT coefficients are:
$a = -\frac{1}{2} \langle \nabla \det DX(p), q_0 \rangle, \quad b = \langle \nabla \text{tr } DX(p), q_0 \rangle$

Significance: This formula is coordinate-free and invariant under phase-space coordinate changes. It requires only the $C^2$ regularity of the field and the 2-jet of the vector field, eliminating the need for generalized eigenvectors or solving linear systems for biorthonormal bases.
Consequence: If the vector field factors as $X = hZ$ where $h(p)=0$ , then $a=0$ automatically. This implies that non-isolated equilibria (curves of equilibria) always satisfy $a=0$ .

2. Geometric Interpretation of Degeneracies (Theorem B)
The coefficients $a$ and $b$ measure the failure of the kernel line field $\ell$ to be transverse to specific hypersurfaces in the equilibrium manifold $E$ :

$a = 0$ if and only if $\ell$ is tangent to the fold set $\Sigma = \{\det DX = 0\}$ .
$b = 0$ if and only if $\ell$ is tangent to the neutrality set $\tilde{H} = \{\text{tr } DX = 0\}$ .
The paper clarifies that $b=0$ does not imply a loss of transversality between $\Sigma$ and $\tilde{H}$ themselves, but rather a specific alignment of the kernel direction relative to $\tilde{H}$ .

3. Mixed-Volume Obstruction (Theorem C)
For Kolmogorov systems with Laurent polynomial numerators $A$ and $B$ , the paper proves a combinatorial obstruction:

If $A$ and $B$ share a common nonconstant factor vanishing at the equilibrium, the equilibrium is non-isolated and $a=0$ .
If the equilibrium is isolated, the order of the Takens normal form term $f$ satisfies $\text{ord } f \leq MV(\text{Newt } A, \text{Newt } B)$ .
Critical Bound: If $MV(\text{Newt } A, \text{Newt } B) \leq 2$ , then $a \neq 0$ for any isolated equilibrium. Consequently, nilpotent singularities of saddle, focus, or elliptic type (which require $\text{ord } f \geq 3$ ) are unreachable at isolated equilibria in such systems. Only cusp-type singularities ( $\text{ord } f = 2$ ) can occur.
The bound constrains the first coefficient ( $a$ ) but imposes no restriction on the second coefficient ( $b$ ), meaning $b=0$ (degenerate cusp) remains possible even when $MV \leq 2$ .

4. Application to Specific Classes

Cross-Product Cubic Systems: The paper analyzes systems where cubic terms are restricted to cross-products ( $x^2y, xy^2$ ). Here, $MV=2$. The authors explicitly identify the parameter values where $a=0$ ; these correspond exactly to the loss of isolation (common factors in $P$ and $Q$ ).
Ecological Models: Two classical models are examined:
1. Predator Satiation (Holling Type II): $MV \leq 2$ . The analysis confirms $a \neq 0$ for isolated equilibria. A unique degenerate point exists where $b=0$ .
2. Mutual Interference (Beddington–DeAngelis): $MV \leq 2$ . Similarly, $a \neq 0$ is guaranteed for isolated points, and a unique degenerate point with $b=0$ is found.
  In both cases, the degeneracy $b=0$ yields a cusp-type singularity of codimension at least three.

Significance and Claims
The paper claims to provide a unified framework for understanding Bogdanov–Takens nondegeneracy through three complementary lenses:

Intrinsic: The coefficients are directional derivatives of Jacobian invariants, simplifying computation and removing reliance on specific coordinate normalizations.
Geometric: Nondegeneracy is determined by the relative position of the kernel line field with respect to the fold and neutrality sets.
Combinatorial: For Kolmogorov systems, the vanishing of $a$ is not merely an algebraic coincidence but a consequence of the mixed volume of the Newton polytopes.

The authors emphasize a genuine asymmetry revealed by their results: while the first coefficient ( $a$ ) is controlled by the combinatorial structure of the system's supports (Newton polytopes), the second coefficient ( $b$ ) remains beyond the reach of this combinatorial bound. The paper concludes that for the class of systems with $MV \leq 2$ , the only possible isolated nilpotent singularities are topological cusps; the more complex saddle, focus, and elliptic types are structurally forbidden at isolated equilibria.

More like this