Floquet–Magnus Expansion in Periodic Quantum Systems

• Floquet–Magnus expansion is a framework that derives a time-independent effective Hamiltonian for periodic quantum drives using nested commutators.
• It quantifies noncommutativity and micromotion effects with rigorous error control and convergence criteria in high-frequency regimes.
• Applications span quantum control, solid-state NMR, and inverse Floquet engineering, enabling effective Hamiltonian design in various systems.

The Floquet–Magnus expansion is a perturbative construction of a time-independent Floquet Hamiltonian for a periodically driven system with $H(t)=H(t+T)$, defined through the one-period propagator

$$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$

In periodic quantum dynamics it provides an effective stroboscopic generator, and in formulations that keep micromotion explicit it yields a factorization $U(t)=P(t)e^{-iH_F t}$ with $P(t+T)=P(t)$. The modern literature treats the Floquet–Magnus expansion not only as a formal high-frequency series, but also as an asymptotic tool with rigorous finite-order error control, and as an object whose divergence can coexist with localized or metastable dynamics rather than implying immediate heating or ergodicity (Arnal et al., 2024, Kuwahara et al., 2015, Haga, 2019).

1. Definition and canonical formulas

For a $T$-periodic Hamiltonian, the Floquet Hamiltonian $H_F$ is defined by

$$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$

with the short-period requirement $H_F\to H_0$ as $T\to0$ in the convention used for periodically driven quantum systems (Haga, 2019). A common form of the Floquet–Magnus expansion writes

$H_F=\sum_{n=0}^\infty \Omega_n T^n.$

The first coefficients are

$$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$0

$$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$1

$$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$2

(Haga, 2019).

A closely related convention expands the Magnus generator $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$3 in

$$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$4

and then identifies $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$5. In that notation the first term is $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$6, and higher orders are ordered time integrals of nested commutators (Krondorfer et al., 3 Jul 2025, Wei et al., 14 Mar 2026). The two conventions differ only by indexing and normalization; both encode the same nested-commutator structure.

The leading contribution is the time average of the drive over one period. Higher orders quantify noncommutativity between Hamiltonians at different times. In this sense the Floquet–Magnus expansion is more than a static averaging prescription: it resolves the extent to which time ordering fails to collapse to a simple period average.

2. Floquet decomposition, micromotion, and coordinate transformations

Floquet theory permits a decomposition of the exact propagator as

$$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$7

where $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$8 is unitary and $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$9 is a time-independent Hermitian operator (Arnal et al., 2024). One constructive route introduces a time-dependent unitary transformation

$U(t)=P(t)e^{-iH_F t}$0

and chooses $U(t)=P(t)e^{-iH_F t}$1 to be periodic. In this formulation the expansion determines both the effective Hamiltonian $U(t)=P(t)e^{-iH_F t}$2 and the micromotion generator $U(t)=P(t)e^{-iH_F t}$3 order by order (Arnal et al., 2024).

This representation matters because a finite truncation preserves unitarity by construction: $U(t)=P(t)e^{-iH_F t}$4 so the approximation remains a product of unitary exponentials (Arnal et al., 2024). In practice, this is one of the main distinctions between Floquet–Magnus constructions and perturbative schemes that do not preserve the Lie-group structure term by term.

The same logic extends beyond linear finite-dimensional quantum systems. For general nonlinear periodic ODEs $U(t)=P(t)e^{-iH_F t}$5, one may seek a near-identity periodic change of variables $U(t)=P(t)e^{-iH_F t}$6 that transforms the original system into an autonomous equation $U(t)=P(t)e^{-iH_F t}$7. In the continuous-changes-of-variables approach, the recursive step is

$U(t)=P(t)e^{-iH_F t}$8

with $U(t)=P(t)e^{-iH_F t}$9 fixed by period averaging and Lie brackets replacing operator commutators (Casas et al., 2019). In the linear limit this reduces to the standard Floquet–Magnus structure.

In solid-state NMR, an allied formulation writes

$P(t+T)=P(t)$0

with periodic frame-functions $P(t+T)=P(t)$1. That formulation emphasizes an operational point: unlike average Hamiltonian theory restricted to integer multiples of $P(t+T)=P(t)$2, the Floquet–Magnus representation can track dynamics between stroboscopic observation points (Mananga, 2024).

3. Convergence, asymptoticity, and finite-order control

The classical sufficient condition for convergence of the Magnus series on one period is

$P(t+T)=P(t)$3

or, in equivalent notation, $P(t+T)=P(t)$4 for the generator $P(t+T)=P(t)$5 (Krondorfer et al., 3 Jul 2025, Arnal et al., 2024). This criterion is sufficient rather than necessary, and it is often too restrictive for many-body or strong-drive settings.

For generic few-body periodically driven many-body Hamiltonians with local energy scale $P(t+T)=P(t)$6, Kuwahara, Mori, and Saito proved

$P(t+T)=P(t)$7

so that $P(t+T)=P(t)$8. The factorial growth implies that the Floquet–Magnus series is generically asymptotic: partial sums decrease until $P(t+T)=P(t)$9 and then diverge (Kuwahara et al., 2015). In the high-frequency regime $T$0, the optimal truncation

$T$1

satisfies

$T$2

and the corresponding stroboscopic approximation remains accurate up to times $T$3 (Kuwahara et al., 2015). This is the rigorous basis of Floquet prethermalization: for $T$4 with $T$5, dynamics are governed by a quasi-conserved truncated Floquet Hamiltonian, whereas at much later times a generic nonintegrable system heats toward the infinite-temperature ensemble (Kuwahara et al., 2015).

A complementary nonperturbative route derives effective Hamiltonians by iterated integration by parts. In that framework, if $T$6 denotes the Floquet–Magnus truncation through order $T$7, then at stroboscopic times

$T$8

without assuming convergence of the full infinite series (Dey et al., 29 Apr 2025). Applied to the semiclassical Rabi model, this analysis shows that the rotating-wave approximation can outperform the Bloch–Siegert Hamiltonian in most regimes, while the third-order approximation ultimately outperforms both over very long timescales (Dey et al., 29 Apr 2025).

4. Divergence, resonances, and unbounded spectra

For infinite-dimensional systems with unbounded Hamiltonians, the operator norms entering the standard convergence criteria may be ill-defined. In a driven anharmonic oscillator,

$T$9

a basis-truncation procedure gives a practical definition of the convergence radius. One chooses an orthonormal basis, truncates to dimension $H_F$0, computes $H_F$1 and its Magnus coefficients $H_F$2, forms the infinite-$H_F$3 limits of matrix elements, and defines

$H_F$4

(Haga, 2019).

For the purely harmonic case $H_F$5, the ratio $H_F$6 approaches $H_F$7, so $H_F$8. For any $H_F$9, however, the asymptotic behavior is

$$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$0

hence $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$1 and $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$2: the Floquet–Magnus expansion diverges for all driving frequencies once any nonlinearity is present, even if the anharmonicity is arbitrarily small (Haga, 2019). The same work shows that this divergence coexists with exponentially localized Floquet eigenstates in energy space, bounded long-time average energy, and Poissonian rather than Wigner–Dyson level statistics. The divergence therefore does not require quantum ergodicity or heating to infinite temperature (Haga, 2019).

The interpretation proposed for this coexistence invokes quantum resonances: extremely narrow resonant peaks in $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$3 versus $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$4 that proliferate arbitrarily close to $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$5, analogously to the dense resonances of the kicked rotor at rational $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$6. This suggests that the analytic continuation of $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$7 can fail at every nonzero $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$8 even when those resonances are unobservable at finite resolution (Haga, 2019).

Related conclusions appear in the periodically driven Friedrichs models. In the discrete model, where the spectrum is bounded, the Floquet–Magnus expansion converges for $$U(T)=\mathcal T e^{-i\int_0^T H(t)\,dt}=e^{-iH_FT},$$9, and the zeroth-order effective Hamiltonian correctly predicts an exact Floquet bound state. In the continuous model, where the spectrum is unbounded, no true bound state exists at finite $H_F\to H_0$0; instead the Floquet–Magnus prediction survives as a metastable resonant state whose lifetime scales as $H_F\to H_0$1 and diverges as $H_F\to H_0$2 (Mori, 2014). In the low-frequency regime, the same model exhibits Floquet resonances with exponentially small imaginary parts of quasi-energy,

$H_F\to H_0$3

which are interpreted as tunneling in energy space (Mori, 2014).

A recent operator-theoretic extension places these observations on a broader footing. For time-periodic unbounded Hamiltonians, the expansion can be reformulated on the invariant core $H_F\to H_0$4 using relative bounds with respect to a reference self-adjoint operator $H_F\to H_0$5, together with spectral clustering and conjugation symmetry. Under these hypotheses, the finite-order effective Hamiltonian

$H_F\to H_0$6

has coefficients that coincide order by order with the bounded Floquet–Magnus series, and the stroboscopic error obeys

$H_F\to H_0$7

again without requiring convergence of the full series (Burgarth et al., 22 May 2026).

5. Classical, dissipative, and nonlinear generalizations

The Floquet–Magnus expansion is not confined to closed quantum systems. For classical driven dynamics, including stochastic equations, one can pass from the equations of motion to a periodic Fokker–Planck operator $H_F\to H_0$8 and define an effective static generator through

$H_F\to H_0$9

In Fourier components $T\to0$0, the first three terms are

$T\to0$1

$T\to0$2

$T\to0$3

(Higashikawa et al., 2018).

This master-equation formulation applies to the Langevin equation, the Gross-Pitaevskii equation, and the time-dependent Ginzburg-Landau equation (Higashikawa et al., 2018). In many-body classical systems it is again asymptotic rather than convergent in the strict infinite-order sense: $T\to0$4 with effective truncation order $T\to0$5, and the truncated generator approximately commutes with one-period evolution up to exponentially small error (Higashikawa et al., 2018).

Kapitza’s pendulum provides a concrete bridge between classical and quantum uses of the expansion. For a vertically driven pendulum with $T\to0$6, the third-order Floquet–Magnus calculation yields an effective Hamiltonian

$T\to0$7

or equivalently the effective potential

$T\to0$8

The inverted position becomes stable when

$T\to0$9

(Krondorfer et al., 3 Jul 2025). In the dissipative classical treatment with friction, the effective equation of motion becomes

$H_F=\sum_{n=0}^\infty \Omega_n T^n.$0

reproducing the stabilization criterion $H_F=\sum_{n=0}^\infty \Omega_n T^n.$1 and matching exact stroboscopic dynamics up to long times (Higashikawa et al., 2018).

For genuinely nonlinear deterministic systems, the continuous-transformation formalism extends Floquet–Magnus methods to the Van der Pol oscillator and the nonlinear Schrödinger equation on the torus. In that setting the universal commutator integrals are replaced by Lie-bracket expressions for vector fields, while the periodic near-identity map eliminates explicit time dependence order by order (Casas et al., 2019).

6. Engineering and domain-specific applications

One contemporary use of the Floquet–Magnus expansion is inverse Floquet engineering. For a periodically driven oscillator in the rotating frame, Guo and Xu construct a real-space drive

$H_F=\sum_{n=0}^\infty \Omega_n T^n.$2

so that $H_F=\sum_{n=0}^\infty \Omega_n T^n.$3 produces the desired target Hamiltonian at zeroth order, $H_F=\sum_{n=0}^\infty \Omega_n T^n.$4 cancels the induced first-order correction, $H_F=\sum_{n=0}^\infty \Omega_n T^n.$5 cancels the net second-order correction, and so on (Xu et al., 2024). Their six-step iteration computes the error Hamiltonian at order $H_F=\sum_{n=0}^\infty \Omega_n T^n.$6, converts it through a non-commutative Fourier transform, and defines a correction drive whose rotating-wave contribution precisely cancels that error (Xu et al., 2024).

A key technical ingredient is the NcFT “bracket” operation, which translates commutator hierarchies into low-dimensional c-number integrals without explicit operator reordering. This allows the systematic engineering of target Hamiltonians with discrete rotational and chiral symmetries in phase space (Xu et al., 2024). In the example of

$H_F=\sum_{n=0}^\infty \Omega_n T^n.$7

the zeroth-order rotating-wave Hamiltonian has $H_F=\sum_{n=0}^\infty \Omega_n T^n.$8-fold rotation and chiral symmetry, while $H_F=\sum_{n=0}^\infty \Omega_n T^n.$9 and higher terms break those symmetries. Adding the analytically derived $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$00 restores the pairwise quasi-energy degeneracies; numerically, the Floquet eigenstate fidelities increase from $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$01 to $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$02, and the remaining micromotion and $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$03-dependence are suppressed to $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$04 (Xu et al., 2024).

Two-level systems remain a major testing ground because the $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$05 algebra collapses nested commutators to vector cross products. For

$$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$06

the semiclassical Rabi model has $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$07, vanishing second-order correction in the specific convention used there, and third-order correction

$$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$08

(Wei et al., 14 Mar 2026). A parameter scan over $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$09 and $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$10 shows absolute quasienergy errors decreasing from $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$11 to $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$12 to $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$13 at $$U(T,0)=\mathcal T\exp\!\Bigl[-\,i\int_0^T H(t)\,dt\Bigr]=e^{-iH_FT}.$$14 as one goes from first to third to fifth order, with analogous improvements at lower frequency and stronger drive (Wei et al., 14 Mar 2026). Suitable interaction-picture and adiabatic-basis transformations enlarge the regime of satisfactory agreement with exact results, and symmetry under half-period shifts can force all even orders to vanish (Wei et al., 14 Mar 2026).

In solid-state NMR, the Floquet–Magnus expansion has been revisited for the Triple Oscillating Field Technique. There it is emphasized that the method operates directly in the finite-dimensional system Hilbert space rather than in an infinite-dimensional Sambe space, and that its frame-function formulation supplies propagators and effective Hamiltonians that describe the evolution between stroboscopic points (Mananga, 2024). This places Floquet–Magnus methods at the intersection of effective-Hamiltonian theory, spectroscopy, and time-dependent quantum control.

Taken together, these applications show that the Floquet–Magnus expansion serves two technically distinct roles. First, it is a high-frequency or asymptotic reduction tool that exposes prethermal generators, suppressed heating, and metastable resonant structures. Second, it is a constructive design formalism for building target effective Hamiltonians and symmetry-protected Floquet phases. The current literature indicates that these roles remain meaningful even when the full infinite series diverges, provided that truncation, micromotion, and error control are handled with the appropriate operator-theoretic or model-specific framework (Dey et al., 29 Apr 2025, Burgarth et al., 22 May 2026, Haga, 2019).