⚛️ quantum physics

Floquet Theory and Average Hamiltonian Theory Revisited: Equivalence, Convergence and Applications to NMR

This paper establishes the mathematical equivalence between Floquet theory and average Hamiltonian theory for periodically driven quantum systems by utilizing the Floquet-Magnus expansion to derive a more robust, integration-free calculation scheme that, when applied to solid-state NMR experiments, yields significantly higher accuracy through the consistent separation of secular and non-secular contributions.

Original authors: Antonia J. Bock, Matthias Ernst, Götz S. Uhrig

Published 2026-09-16

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Antonia J. Bock, Matthias Ernst, Götz S. Uhrig

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.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

In the quiet, controlled world of a laboratory, scientists often study the tiny magnetic spins of atoms to understand the structure of matter. To do this, they place samples in strong magnetic fields and bombard them with radio waves, listening to the faint signals the atoms send back. However, in solid materials, these signals are often blurred and hard to read because the atoms interact with each other in complex ways. To sharpen the picture, researchers spin the sample rapidly, like a top, while hitting it with precise pulses of radio energy. This spinning and pulsing creates a system that repeats itself over and over, a rhythm that scientists must decode to see the true nature of the material. For decades, two different mathematical toolkits have been used to predict how these spinning atoms will behave. One approach, known as average Hamiltonian theory, simplifies the problem by calculating a single, steady average of the forces at play. The other, Floquet theory, treats the repeating motion as a series of distinct steps, accounting for the rapid, rhythmic fluctuations that happen between the spins. While both methods aim to describe the same physical reality, they often produce slightly different results, leaving researchers to wonder which tool is truly reliable and why they sometimes disagree.

A team of physicists and chemists recently revisited this long-standing puzzle to determine exactly how these two theories relate and which one offers the most accurate path forward. By tracing both methods back to a common mathematical root, the researchers demonstrated that they are fundamentally equivalent in their deepest structure, yet they diverge in how they handle the details when the calculations are stopped at a certain level of complexity. The team developed a new, streamlined way to perform these calculations that avoids the messy, error-prone algebra that has plagued previous attempts. Using this cleaner method, they worked out the first four levels of detail for both theories, a feat that had been considered too cumbersome to achieve before. Their analysis revealed that while the two theories are mathematically twins in the limit of perfect precision, one approach consistently outperforms the other when applied to real-world scenarios. Specifically, the method that keeps track of the rapid, rhythmic "kicks" the system receives during each cycle of motion provides results that are roughly three times more accurate than the method that ignores these rapid fluctuations.

The researchers tested their findings on two classic problems in nuclear magnetic resonance. First, they looked at a single atom interacting with a strong radio wave, a situation known to cause a subtle shift in the atom's signal called the Bloch–Siegert shift. They found that to predict the exact location of this shift correctly, especially when the radio waves are strong, it is essential to include the rapid, rhythmic adjustments in the calculation. Without them, the predicted signal flips its sign entirely, leading to a complete misunderstanding of the experiment. Next, they examined a more complex system of atoms spinning inside a solid sample, mimicking the behavior of protons in a chemical chain. Here, they compared the predictions of both theories against a computer simulation that solved the problem exactly. The results showed that the method including the rhythmic adjustments not only matched the exact solution more closely but also converged to the correct answer much faster as the calculations became more detailed. Even when the researchers looked at the fine details of the spectrum, such as the faint sidebands that appear due to the spinning motion, the method that ignored the rapid kicks failed to reproduce them, while the complete method captured them perfectly.

The study concludes that for scientists working with solid materials, the most robust approach is to use the theory that separates the slow, steady evolution of the system from the fast, rhythmic jitters, keeping both parts in the final calculation. This separation ensures that the predictions remain stable and accurate, even when the experimental conditions are challenging. The researchers emphasize that while the two theories are formally the same in an ideal, infinite calculation, the practical necessity of stopping at a finite number of steps makes the inclusion of the rapid rhythmic terms a critical advantage. By providing a clear, error-free way to calculate these terms up to the fourth level of detail, the team has offered a practical guide for improving the accuracy of experiments that reveal the hidden structures of the solid world. Their work does not just settle a theoretical debate; it provides a concrete, more reliable tool for interpreting the signals that tell us how matter is built.

Technical Summary: Floquet Theory and Average Hamiltonian Theory Revisited

Problem Statement
Periodically driven quantum systems are central to fields ranging from quantum control to solid-state Nuclear Magnetic Resonance (NMR). The dynamics of these systems are governed by the time-dependent Schrödinger equation with a Hamiltonian $H(t) = H(t+T)$ . While exact solutions exist formally via Dyson's formula, analytical integration is often impossible, necessitating approximate methods. Two dominant frameworks are Average Hamiltonian Theory (AHT) and Floquet-based methods (specifically the Floquet–Van Vleck or FVV approach). Despite their shared goal of approximating system dynamics, they often yield seemingly different or conflicting results in practice. The NMR community frequently utilizes AHT as a standard tool, while Floquet methods offer an alternative formulation suited for periodic Hamiltonians and sideband structures. However, the precise mathematical equivalence, the origin of discrepancies, and the conditions under which one approach outperforms the other remain unclear. Furthermore, the algebraic complexity of higher-order expansions in both theories often leads to calculation errors and obscures their formal relationship.

Methodology
The authors revisit AHT and Floquet theory through a unified perspective based on the Floquet–Magnus Expansion (FME). The FME serves as the general framework from which both AHT and FVV are derived as specific representations.

Unification: The paper demonstrates that AHT and FVV arise from the same underlying FME structure but differ in their choice of gauge (initial conditions for the kick operator $\Lambda(t)$ ). AHT corresponds to setting $\Lambda(t_0) = 0$ , while FVV corresponds to setting the zeroth Fourier coefficient of the kick operator to zero ( $\int_0^T \Lambda(t) dt = 0$ ).
Integration-Free Derivation: To address the high algebraic complexity and susceptibility to errors in traditional integration-based derivations, the authors introduce a calculation scheme that avoids explicit time integration. By exploiting the periodicity of the kick operator and applying a Fourier transformation to the defining differential equations, they derive the first four orders of both the effective Hamiltonian ( $F$ ) and the kick operator ( $\Lambda$ ) directly in the frequency domain.
Numerical Validation: The derived formulas are numerically validated against exact time-propagators for random and specific time-periodic Hamiltonians (including dipolar-coupled spin systems) across various dimensions (2, 4, 8, 16). The error is measured in terms of eigenvalue and eigenvector deviations.
Application to NMR: The theories are applied to three specific NMR scenarios:
1. The Bloch–Siegert shift (a single spin-1/2 system under a resonant radio-frequency field).
2. Magic-Angle Spinning (MAS) spectra of homonuclear dipolar-coupled two- and three-spin systems (both single crystals and powder averages).
3. Non-stroboscopic measurements in MAS, where intra-stroboscopic sampling is used to resolve sidebands.

Key Contributions

Proof of Equivalence and Divergence: The work establishes that AHT and FVV are perturbatively equivalent up to the order of truncation in $1/\omega$ . However, it clarifies that their truncated forms yield different results due to the handling of non-secular terms and the "kick operator" (the micromotion operator $\Lambda$ ).
Higher-Order Formulas: The authors provide the first four orders of the effective Hamiltonian and kick operator for both AHT and FVV. These are presented in a consistent notation derived via the integration-free method, correcting a known typo in the third-order AHT expression found in existing literature.
Role of the Kick Operator: A critical finding is that the kick operator is essential for mathematical rigor when calculating the full stroboscopic evolution operator and, crucially, for modeling intra-stroboscopic dynamics (sidebands). While the kick operator has a negligible effect on the eigenvalues (resonance positions) in some stroboscopic NMR spectra, it significantly affects the eigenvectors (transition intensities) and is strictly required to reproduce sidebands in non-stroboscopic measurements.
Superiority of FVV: The analysis reveals that the FVV approach is consistently more accurate and numerically robust than AHT. Specifically, the deviation of FVV from the exact solution is approximately three times smaller than that of AHT for the same truncation order.

Results

Accuracy: Numerical validation shows that FVV provides a more accurate approximation of the stroboscopic evolution operator $U(t_0+T, t_0)$ than AHT. The ratio of deviations ( $\Delta_{AHT} / \Delta_{FVV}$ ) averages around 3.1 for eigenvalues and 2.6 for eigenvectors.
Bloch–Siegert Shift: In the Bloch–Siegert shift example, neglecting the kick operator in FVV leads to a sign flip in the spectrum for certain parameter regimes. Third-order FVV is required to accurately predict the resonance shift, matching the numerically exact benchmark.
MAS Spectra: For dipolar-coupled spin systems under MAS, FVV converges faster to the exact line shape than AHT. While the difference between FVV with and without the kick operator is small for stroboscopic powder spectra (affecting only fine details of line shape), the inclusion of the kick operator is necessary for the correct calculation of sidebands in non-stroboscopic measurements.
Sidebands: The study confirms that the kick operator is indispensable for modeling sidebands arising from intra-stroboscopic sampling. Neglecting it erases the micromotion contributions at higher frequencies.

Significance and Claims
The paper claims to provide a unified theoretical background that resolves apparent contradictions between AHT and Floquet theory. By identifying the Floquet–Magnus expansion as the unifying framework, the authors clarify that the differences between AHT and FVV are artifacts of truncation and gauge choice rather than fundamental theoretical incompatibilities.

The authors recommend the Floquet–Van Vleck (FVV) approach, including both the effective Hamiltonian and the kick operator, as the preferred method for analyzing solid-state NMR experiments. They argue that the consistent separation of secular and non-secular contributions in FVV offers superior numerical robustness and accuracy (approximately three times better than AHT) despite the formal equivalence of the theories in the limit of infinite order. The work also highlights that while the kick operator may sometimes be negligible for resonance positions in specific stroboscopic NMR contexts, it is crucial for the correct description of transition intensities and for any analysis involving non-stroboscopic data or sidebands.

The authors modestly note that while their empirical results strongly support the superiority of FVV, a rigorous mathematical proof bounding the remaining deviations for general convergence mechanisms remains an open question for future research.

Technical Summary: Floquet Theory and Average Hamiltonian Theory Revisited

More like this