2026-09-07 | | Total: 2
An alternative, streamlined methodology is presented for deriving the isotropic Cartesian tensor bases under the special orthogonal group $\text{SO}(3)$. By shifting the traditional interpretation of the isotropy condition to analyze it as an algebraic system of equations for the rotation matrices themselves rather than the tensor components, lower-rank bases can be explicitly established in just a few lines without requiring finite coordinate rotations or complicated contractions of infinitesimal generators. Furthermore, a transparent combinatorial interpretation is provided for the number of tensors in the general higher-rank spanning sets. Finally, the Gram-matrix method is advocated as an efficient computational sieve to resolve the subsequent linear dependence.
We present a Hessian-based approach to numerically continue periodic orbits in dynamical systems. A loop (periodic orbit candidate) is parametrized as a Fourier series; a loss function is defined based on the deviation of the loop from the physical differential equations. Unlike previous work relying on hand-derived Jacobians, our method automates the process by leveraging automatic differentiation, a common machine learning technique. The continuation direction can be determined by the flat directions of the loss landscapes (directions with zero eigenvalues), making the search of periodic orbits efficient and guided. Our method is integrator-free, precisely initializes oscillations around unstable fixed points, and efficiently detects orbit family intersections and subharmonic bifurcations. As a demonstration, we present full continuations of periodic double pendulum oscillations from fixed points, showing bifurcations along orbit families and categorizing branches of periodic orbits. In particular, we find periodic orbits where both pendulum masses are never simultaneously at rest, which to our knowledge has been missing in the literature.