A New Approach to Solving Periodic Differential Systems
© The Physical Society of Japan
This article is on
A remark on renormalization group theoretical perturbation in a class of ordinary differential equations
Prog. Theor. Exp. Phys. 2021, 013A02 (2021).
Mathematicians and physicists are well acquainted with secondorder ordinary differential equations (ODE), the most prominent of them being the class of equations that govern oscillatory motion.
Mathematicians and physicists are well acquainted with secondorder ordinary differential equations (ODE), the most prominent of them being the class of equations that govern oscillatory motion. Often, these oscillatorytype equations involve an inhomogeneous “potential” term coupled with a perturbation term and are usually solved as a power series in the perturbation with coefficients that are polynomials in the independent variable. However, when describing timedependent behavior of a system, i.e., the independent variable is time (t), the series solution becomes problematic due to the coefficients diverging in time, which limits the valid description of the system to time scales of the order of (1/perturbation).
In a new study, I report on an elegant (yet unnoticed) functional identity in renormalization group (RG) theoretical perturbation theory, which is useful to get rid of the divergent timedependence and to replace it instead with a description based on “renormalized amplitudes” and their time dynamics.
I begin with a secondorder ODE in time with a potential term depending on the perturbation, the independent variable and its derivative, and the “resonant harmonics,” e^{±it}, and its naïve perturbation series solution. Next, I write the dependent variable, say y, as a perturbation series expansion and substitute it in the ODE to get an equation for each power of the perturbation. I define Y as the resulting formal solution of the power series expansion of y and write it as a Fourier expansion with coefficients that are functions of what I call “bare amplitudes,” A, B. I will formally refer to them as “secular coefficients.”
I go on to prove an identity and a relation between the secular coefficients that are key to my results. I next define the “renormalized amplitudes” A_{r}(t) and B_{r}(t) as a function of the “resonant” secular coefficient and use the proved relation to express the bare amplitudes as a function of A_{r}(t) and B_{r}(t). Using the identity proved before, I further establish the renormalized series expansion in terms of A_{r }and B_{r}, which is free from the secular time dependence.
Further, the dynamics of A_{r}(t) and B_{r}(t) is solely determined by the resonant secular coefficients to all orders of perturbation.
This work has important implications in studies involving the classic cases of Van der Pol, Mathieu, Duffing, and Rayleigh equations and may trigger new theoretical development of a RGbased approach to solving differential equations.
A remark on renormalization group theoretical perturbation in a class of ordinary differential equations
Prog. Theor. Exp. Phys. 2021, 013A02 (2021).
Share this topic
Fields
Related Articles

Random Numbers Can Help Solve Difficult Problems in Manybody Physics
Mathematical methods, classical and quantum physics, relativity, gravitation, numerical simulation, computational modeling
Electronic transport in condensed matter
Magnetic properties in condensed matter
2021329
Theorists review a random state vectorbased description of quantum manybody systems which helps greatly reduce the computational burden involved in their numerical simulations, opening doors to applications in quantum computing.

Hybrid Quantum–Classical Algorithms: At the Verge of Useful Quantum Computing
Mathematical methods, classical and quantum physics, relativity, gravitation, numerical simulation, computational modeling
Crossdisciplinary physics and related areas of science and technology
2021322
Scientists discuss the recent progress in algorithms that have enabled hybrid quantum–classical computers, which has brought the quest to realize useful quantum computing much closer to its finish line.

Solving Quantum Equations with Gauge Fields: How Explicit Integrators Based on a Bipartite Lattice and Affine Transformations Can Help
Mathematical methods, classical and quantum physics, relativity, gravitation, numerical simulation, computational modeling
Superconductivity
2021315
We proposed an explicit numerical integrator consisting of affine transformation pairs resulting from the checkerboard lattice for spatial discretization. It can efficiently solve time evolution equations that describe dynamical quantum phenomena under gauge fields, e.g., generation, motion, interaction of quantum vortices in superconductors or superfluids.