Quantum states can be described equivalently by density matrices, Wigner functions, or quantum tomograms. We analyze the accuracy and performance of three related semiclassical approaches to quantum dynamics, in particular with respect to their numerical implementation. As test cases, we consider the time evolution of Gaussian wave packets in different one-dimensional geometries, whereby tunneling, resonance, and anharmonicity effects are taken into account. The results and methods are benchmarked against an exact quantum mechanical treatment of the system, which is based on a highly efficient Chebyshev expansion technique of the time evolution operator.
Quantum states can be described equivalently by density matrices, Wigner functions, or quantum tomograms. We analyze the accuracy and performance of three related semiclassical approaches to quantum dynamics, in particular with respect to their numerical implementation. As test cases, we consider the time evolution of Gaussian wave packets in different one-dimensional geometries, whereby tunneling, resonance, and anharmonicity effects are taken into account. The results and methods are benchmarked against an exact quantum mechanical treatment of the system, which is based on a highly efficient Chebyshev expansion technique of the time evolution operator.
We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.
In this Master's thesis, we give a pedagogical introduction to the semiclassical analysis of wave dynamics in single-band and multi-band systems. We show how the semiclassical evolutions of Wigner functions can approximate the dynamics of wave packet centers and wave packet deformations. Three systems of increasing complexity are considered: the quantum harmonic oscillator and the quantum anharmonic oscillator as single-band systems, and the Dirac oscillator as a two-band system. We explain the origins of the first-order semiclassical corrections (Moyal, Littlejohn-Flynn, and Berry) in the equations of motion and examine their effects on predictive accuracy.
Another focus is set on the conditions for validity of the semiclassical analysis. For both the quantum anharmonic and the Dirac oscillators, a dimensionless parameter controlling the accuracy of the approximation is identified. In multi-band systems, the formalism applies only to states within a single adiabatic subspace, which, for the Dirac oscillator, corresponds to an energy branch. This motivates the construction of single-branch almost coherent states: approximately Gaussian wave packets formed by superpositions of states from a single energy branch. Their centers undergo oscillations with an aperiodic envelope. These states enable comparison between semiclassical Wigner function evolution and the commonly used Sundaram-Niu formalism. We show that semiclassical Wigner function evolution captures deformation-related dynamics, including the initial decay of oscillations, which the Sundaram-Niu formalism cannot describe by construction.
Additionally, exact quantum evolution reveals partial relocalization of these single-branch almost coherent states, a phenomenon that cannot be captured by the semiclassical evolution.
In this thesis, we study semiclassical phase-space methods for quantum evolution in Hermitian and non-Hermitian systems. We first present the dynamics of Gaussian wave packets under non-Hermitian Hamiltonians and interpret them as a classical dynamics for complex Hamiltonians. We use these to derive exact dynamics for wave packets in the quadratic Swanson oscillator. We show that in the case of unbroken PT-symmetry there can be periodic divergences in this system and relate this to the fact that any operator mapping the system to a Hermitian counterpart is unbounded. We apply the semiclassical wave-packet dynamics to two further anharmonic example systems: a PT-symmetric wave guide, a version of which we propose as a filtering device for optical beams, and a non-Hermitian single-band tight-binding model, for which we use classical equations of motion to model both narrowly and widely distributed initial states. We further develop an exact quantum propagator for non-Hermitian dynamics using lattices of wave packets, whose evolution is governed by the semiclassical equations of motion. We demonstrate that this accurately reproduces quantum dynamics compared to the split operator method. Finally we study a Hermitian two-mode many-particle model for bosonic atom-molecule conversion, for which the classical phase-space structure is an orbifold. We show that standard semiclassical tools may be applied to recover features such as the dynamics, spectrum and density of states of the many-particle system.
The time-dependent WKB approximation for a coherent state is expanded to third order around a guiding real trajectory, allowing for the novel treatment of nonlinearity in its semiclassical dynamics, which is generally present in all physical systems far from the classical regime. The result is a closed-form solution consisting of a linear combination of Airy functions and their derivatives multiplied by an exponential. The expression's ability to capture nonlinearity is demonstrated by examining the autocorrelation of initial coherent states in anharmonic systems with few to many contributing periodic orbits. Its accuracy is compared to the quadratic expansion and found to be superior in regimes of ℏ where the curvature begins to be significant, as expected. Moreover, the expression is shown to be a real-trajectory uniformization over two coalescing saddle points that are emblematic of significant curvature. This extends real-trajectory time-dependent wave-packet semiclassical methods to highly anharmonic systems for the first time and establishes their regime of validity.
We investigate the suitability of the Wigner method as a tool for semiclassical dynamics. In spite of appearances, the dynamical time evolution of Wigner phase space densities is found not to reduce to classical dynamics in most circumstances, even as h→0. In certain applications involving highly ’coherent’ density matrices, this precludes direct h-expansion treatment of quantum corrections. However, by selective resummation of terms in the Wigner–Moyal series for the quantum phase space propagation it is possible to arrive at a revised or renormalized classicallike dynamics which solves the difficulties of the direct approach. In this paper, we review the Wigner method, qualitatively introduce the difficulties encountered in certain semiclassical applications, and derive quantitative means of surmounting these difficulties. Possible practical applications are discussed.
In this work we will analyse the capabilities of several numerical techniques for the description of different physical systems. Thereby, the considered systems range from quantum over semiclassical to classical and from few- to many-particle systems. For each case we address an interesting, partly unsolved question. Despite the different topics we address in the individual chapters, the problems under study are somehow related because we focus on the time evolution of the system. In chapter 1 we investigate the behaviour of a single quantum particle in the presence of an external disordered background (static potentials). Starting from the quantum percolation problem, we address the fundamental question of a disorder induced (Anderson-) transition from extended to localised single-particle eigenstates. Distinguishing isolating from conducting states by applying a local distribution approach for the local density of states (LDOS), we detect the quantum percolation threshold in two- and three-dimensions. Extending the quantum percolation model to a quantum random resistor model, we comment on the possible relevance of our results to the influence of disorder on the conductivity in graphene sheets. Furthermore, we confirm the localisation properties of the 2D percolation model by calculating the full quantum time evolution of a given initial state. For the calculation of the LDOS as well as for the Chebyshev expansion of the time evolution operator, the kernel polynomial method (KPM) is the key numerical technique. In chapter 2 we examine how a single quantum particle is influenced by retarded bosonic fields that are inherent to the system. Within the Holstein model, these bosonic degrees of freedom (phonons) give rise to an infinite dimensional Hilbert space, posing a true many-particle problem. Constituting a minimal model for polaron formation, the Holstein model allows us to study the optical absorption and activated transport in polaronic systems. Using a two-dimensional variant of the KPM, we calculate for the first time quasi-exactly the optical absorption and dc-conductivity as a function of temperature. Concerning the numerical technique, the close relation to the time evolution in the other chapters get clear if we identify temperature with an imaginary time. In chapter 3 we come back to the time evolution of a quantum particle in an external, static potential and investigate the capability of semiclassical approximations to it. Considering various one-dimensional geometries, we address basic quantum effects as tunneling, interference and anharmonicity. The question is, to which extend and at which numerical costs, several semiclassical methods can reproduce the exact result for the quantum dynamics, calculated by Chebyshev expansion. To this end we consider the linearised semiclassical propagator method, the Wigner-Moyal approach and the recently proposed quantum tomography. A conceptually very interesting aspect of the compared semiclassical methods is their relation to different representations of quantum mechanics (wave function/density matrix, Wigner function, quantum tomogram). Finally, in chapter 4 we calculate the dynamics of a classical many-particle system under the influence of external fields. Considering a low-temperature rf-plasma, we investigate the interplay of the plasma dynamics and the motion of dust particles, immersed into the plasma for diagnostic reasons. In addition to the huge number of involved particles, the numerical description of this systems faces the challenge of a large range of involved time and length scales. Exploiting the mass differences of plasma constituents and dust particles allows for separating the PIC description of the plasma from the MD simulation of the dust particles in the effective surrounding plasma.