Generalized Time-Coarse Graining via an Operator Cumulant Expansion

Abstract

We introduce a general framework for deriving effective dynamics from arbitrary time-dependent generators, based on a systematic operator cumulant expansion. Unlike traditional approaches, which typically assume periodic or adiabatic driving, our method applies to systems with general time dependencies and is compatible with any dynamics generated by a linear operator -- Hamiltonian or not, quantum or classical, open or closed. This enables modeling of systems exhibiting strong modulation, dissipation, or non-adiabatic effects. Our approach unifies Hamiltonian techniques such as Lie-transform Perturbation Theory (LPT) with averaging-based methods like Time-Coarse Graining (TCG), revealing their structural equivalence through the lens of generalized cumulants. It also clarifies how non-Hamiltonian terms naturally emerge from averaging procedures, even in closed systems. We illustrate the power and flexibility of the method by analyzing a damped, parametrically driven Kapitza pendulum, a system beyond the reach of standard tools, demonstrating how accurate effective equations can be derived across a wide range of regimes.

Figures (10)

• Figure 1: Illustration of Hamiltonian vs. averaging methods. The $t$ and $\tau$ axes denote slow and fast time-scales. Hamiltonian methods split phase space into slow and fast parts, each governed by its own Hamiltonian. Averaging methods smooth over the fast scale, yielding non-Hamiltonian dynamics.
• Figure 2: $T_w$ should act like a filtering operator, retaining only the dynamics of interest.
• Figure 3: Shaded regions indicate stable inverted pendulum regimes. Blue lines show standard analytical boundaries; orange lines are computed from our seventh-order effective generator.
• Figure 4: (top) Phase-space dynamics obtained from the exact (numeric) time-coarse grained dynamics versus the effective time-coarse grained dynamics obtained using our method. (bottom) The phase-space area obtained using the method. The time averaging causes the phase-space area to shrink, a manifestly non-Hamiltonian effect. Our numerics assume the following model parameters: $\Gamma^2 = 0.02, \lambda = 0.3, \nu = 20, \beta = 0.05, A = 0.2, T = 5/2\pi$, with a time-coarse graining scale of $\tau=0.4$.
• Figure 5: As the drive strength $\lambda$ is ramped up, limit cycles first emerge due to the parametric resonance at side-band frequencies $\nu/2$ and $3\nu/2$. The limit cycles are not symmetrical in general with nonzero $\alpha_{1}$ and $\alpha_{2}$ which comes from the dissipative dynamics.