🔢 mathematics

On removing orders from amplitude equations

This paper introduces a modified renormalization group method that utilizes new homogeneous functions to systematically remove removable terms from amplitude equations, revealing a non-removable "core" that prevents linear growth while enabling higher numerical accuracy without increasing complexity.

Original authors: David Juhasz, Per Kristen Jakobsen

Published 2026-08-11

📖 6 min read🧠 Deep dive

CC0 1.0

Original authors: David Juhasz, Per Kristen Jakobsen

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.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 you are trying to predict the future of a complex system, like a swarm of bees, a beating heart, or a planet orbiting a star. These systems are messy; they wiggle, wobble, and interact in ways that are hard to pin down. In the world of mathematics, scientists use a tool called "perturbation theory" to make sense of this chaos. Think of it like trying to describe a bumpy road. You start by describing the flat, smooth highway (the easy part), and then you add little bumps one by one to account for the roughness. The problem is, as you add more and more bumps to get a perfect picture, the math gets so incredibly complicated that it becomes impossible to solve. It's like trying to write a novel where every sentence requires a footnote explaining the grammar of the previous sentence.

For decades, mathematicians have used a clever trick called the "Renormalization Group" (RG) method to tidy up these messy equations. It's like a magic eraser that wipes away the confusing, infinite parts of the math, leaving behind a clean, manageable "amplitude equation" that tells you how the system behaves over the long haul. But even with this magic eraser, getting a super-precise prediction usually meant the final equation had to become huge and unwieldy. You had to choose: either have a simple equation that was less precise, or a super-accurate equation that was too complex to use.

Now, two mathematicians, David Juhasz and Per Kristen Jakobsen, have come up with a new twist on this old trick. They asked a simple question: "What if we could keep the equation simple and make it super accurate at the same time?" Their answer is a modified version of the RG method that acts like a surgical scalpel. Instead of just erasing the messy parts, they introduce a new, flexible ingredient at every step of the calculation. This ingredient allows them to surgically remove the complicated terms from the final equation without losing any precision. However, they discovered a fascinating limit to this power: you can't remove everything. There is a tiny, stubborn "core" of the equation that must stay, or else the math breaks and starts growing uncontrollably. It's as if the system has a skeleton that holds it together; you can strip away the fat and muscle to make it lighter, but you can't remove the bones.

The Story of the "Core" and the Magic Eraser

In their paper, Juhasz and Jakobsen tested this new method on several mathematical puzzles, ranging from a simple bouncing ball (the Duffing equation) to a model of predators and prey (the Lotka-Volterra system).

Here is how their method works, using a playful analogy. Imagine you are trying to describe the motion of a swing.

The Old Way (Classical RG): You start with a basic description. To make it more accurate, you add more and more complex terms to your description. Eventually, your description becomes a giant, tangled knot of words that is hard to read.
The New Way (Modified RG): The authors realized that when they were solving the math, they were ignoring a "ghost" solution—a part of the math that was technically there but usually thrown away. They decided to keep this ghost. By introducing this ghost as a flexible variable, they could use it to cancel out the messy, complicated terms in their final equation.

The result? They managed to strip the amplitude equation down to its bare essentials. For the Duffing equation (a cubic oscillator), they reduced a complex equation with many terms down to a single, simple line: $A' = -3/2 i \epsilon A^2 A^*$ . This is the "core."

But here is the catch, and it's the most important part of their discovery: You cannot remove the core.

As they tried to remove even more terms, they found that the math started to behave badly. The solution began to grow linearly, meaning it would shoot off to infinity, which doesn't make sense for a swinging pendulum or a stable planet. They realized that the "core" is the mathematical skeleton of the system. It's the specific nonlinear term that stops the system from just being a simple, boring oscillation. If you try to remove this core, you break the physics.

Testing the Theory: Does it Actually Work?

The authors didn't just stop at the theory; they put their method to the test with numbers.

The Duffing Equation: They compared their new, simplified solution against a high-precision computer simulation. They found that their method was incredibly accurate, matching the computer's results almost perfectly up to $t = 10,000$ (which is $1/\epsilon^4$ for their test parameter $\epsilon = 0.1$ ). Even better, their error was smaller than the old, more complex method. It seems that by simplifying the equation, they actually made it more robust.
The Van der Pol Oscillator: This one was trickier. When they first applied their method, the error was actually worse than the old method—about 100 times worse! Why? Because their new method introduced "free parameters" (like knobs you can turn) that were initially set to zero. But once they used a computer optimization technique called "gradient descent" to tune these knobs, the error dropped by a factor of 100. Suddenly, their simplified equation was just as good as the complex one.
Predator-Prey (Lotka-Volterra): In this system, their method didn't just simplify the equation; it made the solution itself simpler and reduced the error by nearly 10 times compared to the classical method.
The Selkov Model: This was a test involving a "Hopf bifurcation," a point where a system changes from stable to unstable. Again, they found that by tuning their free parameters, they could match the accuracy of the complex method while keeping the equation much cleaner.

The Bottom Line

The main finding of this paper is that you can simplify the "amplitude equations" used to predict complex systems without sacrificing accuracy. By introducing a new type of flexible function at each step, the authors can strip away the mathematical clutter.

However, they explicitly rule out the idea that you can simplify everything. There is a hard limit. If you try to remove the "core" term (the first nonlinear term that prevents simple oscillation), the math fails. This core is the irreducible heart of the system.

The authors are confident in their results because they have tested them numerically on multiple different systems. They show that their method is not just a theoretical curiosity but a practical tool that produces smaller errors and simpler equations. They suggest that this approach could be a game-changer for solving even harder problems, like the propagation of light through certain materials, where current methods struggle with complexity.

In short, Juhasz and Jakobsen have found a way to declutter the math of the universe. They showed us that while we can't remove the bones of the system, we can definitely trim the fat, making our predictions sharper and our equations easier to read.

Technical Summary: On Removing Orders from Amplitude Equations

Problem Statement
The paper addresses the complexity inherent in perturbation methods used to analyze the long-term behavior and stability of complex dynamical systems. Specifically, it targets the Renormalization Group (RG) method, a technique used to derive amplitude equations that describe the slow evolution of system solutions. In the classical RG approach, achieving higher precision in the asymptotic solution typically requires adding more terms to the amplitude equation, thereby increasing its algebraic complexity. The authors identify a limitation in standard perturbation hierarchies where homogeneous solutions are often ignored at higher orders, potentially missing degrees of freedom that could simplify the resulting amplitude equations without sacrificing accuracy.

Methodology
The authors propose a modified RG method that fundamentally alters the standard perturbation hierarchy. The core innovation is the systematic introduction of arbitrary homogeneous functions at every order of the perturbation expansion.

Modified Perturbation Hierarchy: Unlike the classical method, which typically discards homogeneous solutions after the zeroth order, this approach retains them. New arbitrary constants (e.g., $\alpha_0, \beta_0, \gamma_0$ ) are introduced at each order $\epsilon^n$ .
Renormalization of Homogeneous Constants: These new constants are treated as slowly varying functions of time, similar to the primary amplitudes. They are renormalized using multiplicative constants to absorb secular (divergent) terms arising from the time variable $\xi = \mu - t_0$ .
Derivation of Amplitude Equations: The method derives amplitude equations by ensuring the solution is independent of the arbitrary time scale $\mu$ . This process yields a system of equations where the derivatives of the new homogeneous functions appear.
Strategic Removal of Terms: The authors solve quasi-linear partial differential equations (PDEs) for these homogeneous functions. By choosing specific forms for these functions, they actively cancel out higher-order nonlinear terms in the amplitude equation.
Identification of the "Core": The authors demonstrate that there is a theoretical limit to this simplification. Removing all nonlinear terms would reintroduce linear growth (secular terms proportional to $\log(A)$ ), rendering the expansion invalid. Consequently, a "core" remains: the lowest-order nonlinear term that cannot be removed without breaking the asymptotic ordering. This core corresponds to the system's hypernormal form.

Key Contributions

Simplification of Amplitude Equations: The method allows for the derivation of amplitude equations with significantly fewer terms (often reduced to the "core" term) while maintaining the same order of precision in the final solution.
Enhanced Accuracy: In several cases, the modified method yields solutions with lower numerical error compared to the classical RG method, even when the amplitude equation is simpler.
Optimization via Free Parameters: The homogeneous functions introduce free complex constants. The paper demonstrates that these constants can be tuned (e.g., via gradient descent) to minimize the error of the asymptotic solution, offering a mechanism to improve accuracy beyond standard perturbation limits.
Constructive Approach: The authors frame their method as an algebraic, constructive rediscovery of hypernormal form theory, achieved without relying on the heavy machinery of Lie algebras or differential geometry typically associated with normal form theory.

Results
The method was tested on scalar Ordinary Differential Equations (ODEs) and systems of ODEs:

Duffing Equation: The modified RG method reduced the amplitude equation to a single nonlinear term ( $A' \propto A^2A^*$ ), whereas the classical method retained higher-order terms. Numerical tests showed the modified solution had smaller errors than the classical solution over long time scales ( $t \sim 1/\epsilon^4$ ).
Van der Pol Oscillator: The amplitude equation was simplified to contain only linear and one nonlinear core term. While the initial modified solution (with arbitrary constants set to zero) had higher error than the classical method, optimizing the free constants via gradient descent reduced the error to match or exceed the classical performance.
Lotka-Volterra System: The modified method simplified the amplitude equations and reduced the error of the solution by nearly an order of magnitude compared to the classical RG solution.
Second-Order System (Two Amplitudes): The method successfully simplified the amplitude equations for a coupled system, removing non-core terms and reducing numerical error.
Selkov Model (Hopf Bifurcation): The method was applied to a system near a Hopf bifurcation. Similar to the Van der Pol case, the free constants in the core term were optimized to minimize error, resulting in a solution that performed comparably to the classical method over specific time intervals.

Significance and Claims
The paper claims that the modified RG method offers a more efficient pathway to high-precision asymptotic solutions by decoupling the complexity of the amplitude equation from the precision of the solution. The authors assert that their method provides a "core" representation of the system that is equivalent to the hypernormal form but derived through a more accessible, algebraic RG framework.

The significance lies in the ability to:

Maintain high accuracy while keeping the governing amplitude equations minimal.
Utilize the freedom of homogeneous solutions to actively minimize numerical errors.
Provide a constructive alternative to classical normal form theory that requires less specialized mathematical background.

The authors modestly propose that this approach could be generalized to nonlinear Partial Differential Equations (PDEs), specifically citing the potential application to light propagation equations with cubic nonlinearity (e.g., the nonlinear Schrödinger equation), where analytical solutions are difficult to obtain with higher-order terms. However, the paper focuses primarily on ODEs and does not present results for PDEs.

The Story of the "Core" and the Magic Eraser

Testing the Theory: Does it Actually Work?

The Bottom Line

More like this