Bohmian mechanics is an alternative formulation of quantum mechanics that incorporates the familiar and intuitive picture of particles moving along trajectories and yet predicts the same results as the more widely accepted Copenhagen interpretation. In recent years there has been renewed interest in this Bohmian view, in part for the novel approach that it suggests to certain problems, such as decay processes, both from a theoretical and computational stand point. In this thesis we focus on using the concepts introduced by the Bohmian framework as a practical computational tool. I evaluate a number of implementations of the Bohmian method, get a sense of their strengths and weaknesses and attempt to overcome some stability issues that arise. For problems in one-dimension (lD), accurate solutions of the time-dependent Schrodinger equation produce a wave function from which Bohmian trajectories can be computed by integrating along flux lines. For direct integration of the quantum Hamilton-Jacobi equations, the main problems that arise are related to evauating the quantum potential (QP), especially in regions of low probability density. Sufficient accuracy is required to avoid unphysical trajectory crossings. A number of interpolation schemes were investigated, and smoothed splines with special treatment of edge effects gave the best results. For problems in 2D the alternating direction implicit (ADI) method was employed to produce the wave function. Ways of dealing with unphysical reflections from the boundaries of a finite size domain were studied. The use of cellular automata, especially the lattice-Boltzmann method (LBM) were also considered. Here Bohm trajectories would be propagated by following a small set of rules. The main problem identified is that, unless a scheme can be found in which the quantum potential is self-generating from an equation of continuity, the overhead of computing the QP at each time step, is prohibitive.
Bohmian mechanics is an alternative formulation of quantum mechanics that incorporates the familiar and intuitive picture of particles moving along trajectories and yet predicts the same results as the more widely accepted Copenhagen interpretation. In recent years there has been renewed interest in this Bohmian view, in part for the novel approach that it suggests to certain problems, such as decay processes, both from a theoretical and computational stand point. In this thesis we focus on using the concepts introduced by the Bohmian framework as a practical computational tool. I evaluate a number of implementations of the Bohmian method, get a sense of their strengths and weaknesses and attempt to overcome some stability issues that arise. For problems in one-dimension (lD), accurate solutions of the time-dependent Schrodinger equation produce a wave function from which Bohmian trajectories can be computed by integrating along flux lines. For direct integration of the quantum Hamilton-Jacobi equations, the main problems that arise are related to evauating the quantum potential (QP), especially in regions of low probability density. Sufficient accuracy is required to avoid unphysical trajectory crossings. A number of interpolation schemes were investigated, and smoothed splines with special treatment of edge effects gave the best results. For problems in 2D the alternating direction implicit (ADI) method was employed to produce the wave function. Ways of dealing with unphysical reflections from the boundaries of a finite size domain were studied. The use of cellular automata, especially the lattice-Boltzmann method (LBM) were also considered. Here Bohm trajectories would be propagated by following a small set of rules. The main problem identified is that, unless a scheme can be found in which the quantum potential is self-generating from an equation of continuity, the overhead of computing the QP at each time step, is prohibitive.
Bohmian mechanics, developed by Bohm in 1952, provides an alternative interpretation to nonrelativistic quantum mechanics [1-3]. In the hydrodynamic formulation of quantum mechanics, the continuity equation and the quantum Hamilton-Jacobi equation (QHJE) are obtained by substituting the wave function expressed in terms of the real amplitude and the real action function into the time-dependent Schrödinger equation. Bohm’s analytical approach has been used to compute and interpret realvalued quantum trajectories from a precomputed wave function for a diverse range of physical processes [3-8]. In the synthetic approach, the quantum trajectory method has been developed as a computational tool to generate the wave function by evolving ensembles of real-valued quantum trajectories through the integration of the hydrodynamic equations on the fly [9]. Remarkable progress has been made in the use of T&F Cat # K11224, Chapter 18, Page 284, 15-9-2010 real-valued quantum trajectories for solving a wide range of quantum mechanical problems [10].
Several theoretical methods for the computation of quantum dynamical quantities are formulated, implemented, and applied, in the overall context of overcoming the exponential scaling property of quantum mechanics. The convergence of the approximate forward-backward semiclassical dynamics (FBSD) method has been accelerated by working with only the imaginary part of the real-time position and momentum time correlation functions. Simple manipulations allow the computationally favorable limit of a single factorization (“bead”) of the Boltzmann operator to yield greatly improved results. Numerical examples for a one-dimensional anharmonic oscillator and Lennard-Jones fluids are presented. Bohm's trajectory formulation of quantum mechanics is an exact mapping of the time-dependent Schrödinger equation onto classical fluid dynamics, plus one extra term responsible for all non-classical effects, the dreaded quantum potential. A computational procedure has been introduced to stably solve Bohm's hydrodynamic equations of motion in the forward-backward quantum dynamics (FBQD) formulation of time correlation functions, and applied to one-dimensional bound systems at zero and low temperature. A key step in the formulation of this procedure was the use of Hamilton's law of varying action, a generalization of Hamilton's principle of stationary action applicable to initial-value problems. This has resulted in the first reported synthetic computation of non-stationary quantum trajectories for a fully anharmonic oscillator. Success in one dimension has laid the groundwork for future extension to a rigorous quantum-classical scheme, which is briefly discussed. Other unpublished methods for obtaining time correlation functions and propagators in quantum dynamics are also discussed.
Bohmian mechanics as formulated originally in 1952, has been useful in the implementation of numerical methods applied to quantum mechanics. The scientific community though has had ever since a critical thought about it. Therefore, there are still points to be clarified and rectified. The two main problems are basically: Bohmian mechanics gives a privilege role to the position representation. Secondly, the current interpretation of Bohmian trajectories has been recently proven wrong. In this context, in Chapter 2, new complex Bohmian quantities are defined; so that they allow the capacity to formulate Bohmian mechanics in any arbitrary continuous representation, for instance, the momentum representation. This Chapter is fully based on two articles, regarding the proposed complex Bohmian formulation and its extension into momentum space. Chapter 3 deals with a redefinition and reinterpretation of the Bohmian trajectories from the handling of the continuity equation, this is done without any need of additional postulates or interpretations. Also, it is proved that Bohmian mechanics is actually more than a projective aspect of the Wigner function. As a third point, Chapter 4 presents a sytematic treatment of the hydrodynamic scheme of Bohmian mechanics. Then, a brief summary of the transport equations in Bohmian mechanics is done. Next, a unified hydrodynamic treatment is found for the Bohmian mechanics. This treatment is useful to sketch, a Bohmian treatment to efficiently find the steady value of the transmission integral. In Chapter 5 conclusions of this thesis are drawn.
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.
This chapter provides a comprehensive overview of the Bohmian formulation of quantum mechanics. It starts with a historical review of the difficulties found by Louis de Broglie, David Bohm, and John S. Bell to convince the scientific community about the validity and utility of Bohmian mechanics. Then, a formal explanation of Bohmian mechanics for nonrelativistic, single-particle quantum systems is presented. The generalization to many-particle systems, where the exchange interaction and the spin play an important role, is also presented. After that, the measurement process in Bohmian mechanics is discussed. It is emphasized that Bohmian mechanics exactly reproduces the mean value and temporal and spatial correlations obtained from the standard, that is the Copenhagen or orthodox, formulation. The ontological characteristics of Bohmian mechanics provide a description of measurements as another type of interaction without the need for introducing the wave function collapse. Several solved problems are presented at the end of the chapter, giving additional mathematical support to some particular issues. A detailed description of computational algorithms to obtain Bohmian trajectories from the numerical solution of the Schrodinger or the Hamilton-Jacobi equations are presented in an appendix. The motivation of this chapter is twofold: first, as a didactic introduction to Bohmian formalism, which is used in the subsequent chapters, and second, as a self-contained summary for any newcomer interested in using Bohmian mechanics in his or her daily research activity.
In the last few years the hydrodynamic formulation of quantum mechanics, equivalent to the Bohmian equations of motion, has been used to obtain numerical solutions of the Schrodinger equation. Problems, however, have been experienced near wave function nodes (or low probability regions). Here we attempt to compute wave functions and Bohmian trajectories for the interference of one particle or of two identical particles. It turns out that the large number of nodes (i.e. interference minima) makes the hydrodynamic equations impractical, whereas a more straightforward solution of the Schrodinger equation gives very good results.