Masters Theses

Date of Award

8-2008

Degree Type

Thesis

Degree Name

Master of Science

Major

Chemistry

Major Professor

Robert J. Hinde

Committee Members

Jeffrey Kovac, Charles S. Feigerle

Abstract

Because the Schrödinger equation cannot be solved analytically for systems larger than hydrogen, approximate solution methods must be applied to study systems of chemical interest. Quantum Monte Carlo is one of these methods. Quantum Monte Carlo simulations involve the diffusion of walkers through configuration space and necessitate the use of one of two imaginary time propagators, 2nd or 4th order. The choice of the propagator implemented depends on a balance between accuracy and efficiency. Obtaining accurate solutions is the primary objective; once this is accomplished, finding ways to make the simulations efficient is the next intention. Using a model system of a multidimensional harmonic oscillator, this study intends to compare and contrast both 2nd and 4th order propagators' diminished accuracy as the dimensionality of the system is increased. Two phases of the propagation are sampled to determine accuracy: a single propagation step and the asymptotic limit. To test the accuracy after one step, we start a simulation by placing walkers at the origin, moving them forward one time step in imaginary time, creating normalized histograms of their positions and comparing these with the known result for the imaginary time propagation of the delta function centered at x=0. The convergence of the simulation is determined by observing the change in distribution from the known asymptotic limit distribution. This is done by distributing the walkers in the known asymptotic distribution and observing the change in average weight upon further propagation. A converged distribution’s average weight will not change.

Download
Files over 3MB may be slow to open. For best results, right-click and select "save as..."

Included in

Chemistry Commons

Share

COinS