🔬 mesoscale physics

Semiclassical Analysis of Wave Dynamics via Wigner Transforms: Applications to the Quantum Harmonic, Anharmonic, and Dirac Oscillator

This Master's thesis provides a pedagogical semiclassical analysis of wave dynamics using Wigner transforms across single-band and multi-band systems, demonstrating how first-order corrections capture wave packet deformations and comparing the method's accuracy and limitations against the Sundaram-Niu formalism and exact quantum evolution for harmonic, anharmonic, and Dirac oscillators.

Original authors: August Burg

Published 2026-10-19

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: August Burg

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 microscopic world of atoms and electrons, the rules of physics look very different from the world we see every day. Particles do not follow neat, predictable paths like planets orbiting a sun; instead, they behave like fuzzy clouds of probability, existing in many places at once until they are measured. To understand how these quantum clouds move and change over time, scientists often try to bridge the gap between the strange quantum realm and the familiar classical world we experience. They look for a "semiclassical" approach, a method that treats these particles as if they have a definite position and speed, but adds small corrections to account for their quantum nature. This is a delicate balancing act. If the corrections are too simple, the prediction fails; if they are too complex, the math becomes impossible to solve. The challenge is to find a way to describe the motion of these quantum wave packets that is accurate enough to be useful, yet simple enough to provide a clear picture of what is actually happening.

A recent study by August Burg, a researcher working between institutions in Germany and France, tackles this challenge by examining how wave packets evolve in three specific types of systems: a standard vibrating system, a more complex vibrating system, and a system that mimics the behavior of relativistic particles. The researcher used a mathematical tool called the Wigner transform, which acts like a map of the quantum state, showing where the particle is likely to be and how fast it is moving simultaneously. By tracking how this map changes over time, the study reveals how well semiclassical approximations can predict the future of these quantum systems. The work focuses on understanding when these approximations work, what causes them to fail, and how the shape of the wave packet itself plays a crucial role in the dynamics.

The first system studied was the quantum harmonic oscillator, a model for a particle bouncing back and forth in a smooth, bowl-shaped potential. In this specific case, the semiclassical map turns out to be perfect; it matches the exact quantum behavior exactly. This served as a reliable starting point, confirming that the mathematical tools were working correctly. The researcher then moved to the anharmonic oscillator, a system where the potential is not a perfect bowl but has a slightly different shape. Here, the behavior becomes much more interesting. When a wave packet is set in motion, it does not simply oscillate back and forth forever. Instead, the different parts of the wave packet move at slightly different speeds, causing the packet to stretch, twist, and eventually lose its initial shape. This process is known as decoherence. The study found that the semiclassical map could successfully predict the initial stretching and the decay of the oscillation amplitude. However, it could not predict a phenomenon called "revival," where the wave packet suddenly snaps back to its original shape after a long time. The semiclassical method, by its nature, smooths out the fine details required for this reassembly, meaning it works well for the initial chaos but fails to capture the eventual return to order.

The third system, the Dirac oscillator, is more complex because it involves two different "bands" or energy states that can mix, similar to how an electron can have different spin orientations. In this system, the wave packet is not just a single blob but a combination of two components. The study introduced a new type of starting state called a "single-branch almost coherent state." This is a carefully constructed wave packet that stays within just one of the two energy bands, avoiding the messy mixing that usually happens. By using this specific state, the researcher could apply the semiclassical map to a two-band system for the first time in this context. The results showed that the semiclassical approach could accurately predict the oscillation period and the time it takes for the wave packet to spread out and lose its coherence. However, just like in the anharmonic case, the method could not predict the partial relocalization of the wave packet that occurs later in the quantum evolution. The map showed the packet spreading out and staying that way, while the exact quantum calculation showed it eventually gathering itself up again, a nuance the simplified method missed.

A significant portion of the research was dedicated to testing the limits of a common method used in physics to describe these wave packets. This traditional method, known as the Lagrangian approach, treats the wave packet as a single point moving through space, ignoring its width and shape. The study demonstrated that this point-particle view is insufficient for capturing the full dynamics. By comparing the point-particle prediction with the full map-based prediction, the researcher showed that ignoring the finite width of the wave packet leads to errors in predicting how fast the oscillations decay and how the packet deforms. The study argues that to get an accurate picture of wave dynamics, one must track the entire shape of the probability map, not just its center. This distinction is vital for understanding how quantum systems behave over time, especially when they are not perfectly simple.

The work also identified the specific conditions under which these semiclassical approximations hold true. For the anharmonic oscillator, the accuracy depends on the size of the initial wave packet; larger packets follow the semiclassical rules more closely for longer. For the Dirac oscillator, the accuracy depends on the energy gap between the two bands and the specific shape of the wave packet. The study confirmed that when the wave packet stays within a single energy band, the semiclassical map is a powerful tool. However, the moment the system tries to jump between bands or the wave packet becomes too distorted, the approximation begins to break down. The research concludes that while these semiclassical methods cannot replace exact quantum calculations in every situation, they provide a remarkably clear and intuitive window into the complex dance of quantum wave packets, revealing the mechanisms of spreading and decay that are otherwise hidden in the math.

Technical Summary: Semiclassical Analysis of Wave Dynamics via Wigner Transforms

Problem Statement
Exact analytical solutions for the time evolution of quantum states are rarely available, and numerical solutions often obscure physical intuition. While semiclassical approximations—derived from Hamilton's equations with correction terms—offer a bridge between quantum and classical dynamics, standard derivations often rely on a Lagrangian ansatz that assumes wave packets move as rigid bodies with well-defined centers ( $x_c, p_c$ ) but never change form. This thesis investigates the limitations of such rigid approximations and proposes a more robust semiclassical framework based on the evolution of the full Wigner function. The work addresses three specific questions: (1) What dimensionless parameter controls the validity of the semiclassical approximation when $\hbar$ is a constant? (2) Which states correspond to the adiabatic subspaces required for multiband systems? (3) What are the specific roles and origins of first-order correction terms (Moyal, Littlejohn–Flynn, and Berry curvature) in predictive accuracy?

Methodology
The thesis employs the Wigner-Weyl phase space representation, where quantum operators are mapped to scalar or matrix-valued symbols via the Wigner transform. The dynamics are analyzed by expanding the Moyal bracket (the phase-space commutator) in powers of a small parameter ( $\lambda$ or $\epsilon$ ) representing the scale of non-commutativity relative to system quantities.

The study proceeds through three systems of increasing complexity:

Quantum Harmonic Oscillator (Single-band): Used as a pedagogical baseline where the semiclassical flow is exact.
Quantum Anharmonic Oscillator (Single-band): A system with a non-trivial Moyal term and higher-order corrections, but no internal degrees of freedom (spin/bands).
Dirac Oscillator (Two-band): A relativistic generalization involving a $2 \times 2$ matrix Hamiltonian acting on pseudo-spin, requiring the identification of adiabatic subspaces and the inclusion of Berry curvature and the Littlejohn–Flynn term.

For the multiband Dirac Oscillator, the author constructs "single-branch almost coherent states" ( $|\psi_{\alpha\pm}\rangle$ ). These are superpositions of eigenstates from a single energy branch designed to approximate minimal uncertainty wave packets, ensuring the state remains within a specific adiabatic subspace. The evolution of these states is compared against exact quantum solutions and two semiclassical approaches:

Wigner Function Evolution: Tracking the full phase-space distribution $w(x,p)$ under a semiclassical flow $\phi_t$ .
Lagrangian Ansatz: Tracking only the wave packet center $(x_c, p_c)$ using equations of motion derived from a variational principle (Sundaram–Niu formalism).

Key Contributions and Results

Control Parameters and Validity:
- For the Anharmonic Oscillator, the validity of the semiclassical approximation for initially coherent states is controlled by $\epsilon = 1/(2|\alpha|^2)$ . The approximation holds well for the initial decay of oscillations but fails to predict the "revival" of the wave packet, as the assumption that only phase space regions near the initial radius contribute breaks down during relocalization.
- For the Dirac Oscillator, the control parameter is the rescaled commutator $\epsilon = \hbar\beta c / p_y^2$ . The validity is further enhanced by the effective band gap, which increases with the wave packet's radius in phase space.
Role of Correction Terms:
- Moyal Term: In the Anharmonic Oscillator, neglecting the first-order Moyal term (arising from the non-commutativity of $\hat{x}$ and $\hat{p}$ in the symbol Hamiltonian) significantly degrades the prediction of the oscillation period and decay time.
- Berry Curvature and Littlejohn–Flynn Term: In the Dirac Oscillator, these terms arise from the projection onto adiabatic subspaces. The thesis finds that while including these terms improves the prediction of the evolution of individual spin polarizations ( $\uparrow$ and $\downarrow$ ), their impact on the total expectation value $\langle \hat{x} \rangle$ is less distinct, sometimes even appearing to worsen the prediction at short times due to a "compromise" between the differing dynamical phases of the spin components.
Wave Packet Deformation vs. Rigid Motion:
- The Wigner function evolution successfully captures the deformation of wave packets, including the initial "banana-shaped" stretching in phase space, the decay of oscillation amplitudes, and the squeezing of uncertainties ( $\Delta x, \Delta p$ ).
- The Lagrangian ansatz, which assumes a rigid packet, fails to predict the decay of oscillation amplitudes and the evolution of wave packet width. It only tracks the trajectory of the Wigner function's maximum, missing the physics of the packet's finite width.
- Revivals and Relocalization: Neither semiclassical approach can fully capture the fractional revivals in the Anharmonic Oscillator or the partial relocalization in the Dirac Oscillator, as these phenomena rely on the interference of phase-space regions that the semiclassical flow averages out or neglects.
Single-Branch Almost Coherent States:
- The thesis demonstrates that standard coherent states in the Dirac Oscillator (constructed from all branches) cannot be described by a single-band semiclassical flow. By constructing states restricted to a single adiabatic branch, the author enables a valid comparison between exact quantum dynamics and the semiclassical Wigner evolution, revealing aperiodic oscillations with decaying envelopes.

Significance and Claims
The paper claims that the semiclassical evolution of the Wigner function provides a superior and more intuitive framework for analyzing wave packet dynamics compared to the widely used Lagrangian ansatz for wave packet centers. Specifically:

It correctly identifies the dimensionless parameters that govern the semiclassical limit, moving beyond the formal $\hbar \to 0$ limit.
It explicitly delineates the origins of first-order corrections (Moyal vs. geometric terms) and quantifies their impact on predictive accuracy in both single-band and multi-band systems.
It highlights that the "rigid packet" assumption inherent in the Lagrangian approach is a significant source of error, particularly for predicting the decay of oscillations and the evolution of uncertainty.
It proposes a modified semiclassical ansatz for wave packets that incorporates the evolution of the full Wigner function, thereby capturing deformation-related dynamics that the Sundaram–Niu formalism cannot describe by construction.

The work concludes that while semiclassical methods cannot capture non-perturbative effects like full quantum revivals, they offer a powerful tool for characterizing the initial dynamics, oscillation periods, and decay times of wave packets in complex systems, provided the appropriate adiabatic subspaces and control parameters are identified.

More like this