Publications
2007
In recent years, the use of joint time-frequency representations to characterize and interpret shaped femtosecond laser pulses has proven to be very useful. However, the number of points in a joint time-frequency representation is daunting as compared with those in either the frequency or time representation. In this article we introduce the use of the von Neumann representation, in which a femtosecond pulse is represented on a discrete lattice of evenly spaced time-frequency points using a non-orthogonal Gaussian basis. We show that the information content in the von Neumann representation using a lattice of √N points in time and √N points in frequency is exactly the same as in a frequency (or time) array of N points. Explicit formulas are given for the forward and reverse transformation between an iV-point frequency signal and the von Neumann representation. We provide numerical examples of the forward and reverse transformation between the two representations for a variety of different pulse shapes; in all cases the original pulse is reconstructed with excellent precision. The von Neumann representation has the interpretational advantages of the Husimi representation but requires a bare minimum number of points and is stably and conveniently inverted; moreover, it avoids the periodic boundary conditions of the Fourier representation.
We present a unified derivation of Bohmian methods that serves as a common starting point for the derivative propagation method (DPM), Bohmian mechanics with complex action (BOMCA), and the zero-velocity complex action method (ZEVCA). The unified derivation begins with the ansatz Ψ = eiSlh where the action (S) is taken to be complex, and the quantum force is obtained by writing a hierarchy of equations of motion for the phase partial derivatives. We demonstrate how different choices of the trajectory velocity field yield different formulations such as DPM, BOMCA, and ZEVCA. The new derivation is used for two purposes. First, it serves as a common basis for comparing the role of the quantum force in the DPM and BOMCA formulations. Second, we use the new derivation to show that superposing the contributions of real, crossing trajectories yields a nodal pattern essentially identical to that of the exact quantum wavefunction. The latter result suggests a promising new approach to deal with the challenging problem of nodes in Bohmian mechanics.
We present a new semiclassical method that yields an approximation to the quantum mechanical wavefunction at a fixed, predetermined position. In the approach, a hierarchy of ODEs are solved along a trajectory with zero velocity. The new approximation is local, both literally and from a quantum mechanical point of view, in the sense that neighboring trajectories do not communicate with each other. The approach is readily extended to imaginary time propagation and is particularly useful for the calculation of quantities where only local information is required. We present two applications: the calculation of tunneling probabilities and the calculation of low energy eigenvalues. In both applications we obtain excellent agrement with the exact quantum mechanics, with a single trajectory propagation.
In recent years, intensive effort has gone into developing numerical tools for exact quantum mechanical calculations that are based on Bohmian mechanics. As part of this effort we have recently developed as alternative formulation of Bohmian mechanics in which the quantum action S is taken to be complex [Y. Goldfarb et al., J. Chem. Phys. 125, 231103 (2006)]. In the alternative formulation there is a significant reduction in the magnitude of the quantum force as compared with the conventional Bohmian formulation, at the price of propagating complex trajectories. In this paper we show that Bohmian mechanics with complex action is able to overcome the main computational limitation of conventional Bohmian methods - the propagation of wave functions once nodes set in. In the vicinity of nodes, the quantum force in conventional Bohmian formulations exhibits rapid oscillations that present a severe numerical challenge. We show that within complex Bohmian mechanics, multiple complex initial conditions can lead to the same real final position, allowing for the accurate description of nodes as a sum of the contribution from two or more crossing trajectories. The idea is illustrated on the reflection amplitude from a one-dimensional Eckart barrier. We believe that trajectory crossing, although in contradiction to the conventional Bohmian trajectory interpretation, provides an important new tool for dealing with the nodal problem in Bohmian methods.