A new approach to the determination of good action-angle variables for coupled oscillator systems

Lynn C. Geiger, George C Schatz

Research output: Contribution to journalArticle

1 Citation (Scopus)

Abstract

A theory of action-angle variables for coupled oscillator systems is developed which involves solving the Schrödinger equation using a basis of WKB eigenfunctions, then using the logarithm of the resulting wavefunction to define the generator for the canonical transformation which determines the action-angle variables. This theory is based on the marriage between Miller's method for solving the Hamilton-Jacobi equation using the logarithm of the generating function, and the Ratner-Buch-Gerber method for solving the Schrödinger equation using WKB basis functions. A perturbation-theory analysis of this theory indicates that the semiclassical eigenvalues and canonical transformations obtained from it should become identical to their exact classical counterparts in the limit of large actions for each vibrational mode. Two methods for systematically improving the theory for the lower eigenstates are also proposed. Numerical applications of the theory are presented for two systems, the Morse oscillator and the Henon-Heiles two-mode hamiltonian. The resulting semiclassical eigenvalues are in excellent agreement with their exact quantum counterparts, with the magnitude of the error roughly independent of the energy of the eigenstate. Analogous good agreement is found in comparing the approximate and exact classical canonical transformations. In particular, for the Morse oscillator, good results are obtained for certain higher energy states where second-order classical perturbation theory makes serious errors. Other information examined includes surfaces of section for the Henon-Heiles system (comparing the analytical functions obtained from the present theory with results based on exact trajectory calculations) and vibrational distributions chosen to simulate trajectory calculations (using the present theory to determine bin boundaries for a histogram calculation). Again, the comparison in each case with accurate results is excellent, with maximum errors in action calculations of 0.02 h, and in angle calculations of 0.01 rad.

Original languageEnglish
Pages (from-to)431-447
Number of pages17
JournalChemical Physics
Volume79
Issue number3
DOIs
Publication statusPublished - Sep 15 1983

Fingerprint

oscillators
eigenvectors
logarithms
eigenvalues
perturbation theory
Trajectories
trajectories
Hamilton-Jacobi equation
Hamiltonians
Bins
Wave functions
histograms
Eigenvalues and eigenfunctions
Electron energy levels
vibration mode
generators
energy

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Spectroscopy
  • Atomic and Molecular Physics, and Optics

Cite this

A new approach to the determination of good action-angle variables for coupled oscillator systems. / Geiger, Lynn C.; Schatz, George C.

In: Chemical Physics, Vol. 79, No. 3, 15.09.1983, p. 431-447.

Research output: Contribution to journalArticle

@article{fa0eaaa5c37042bc82892f4b2ed06354,
title = "A new approach to the determination of good action-angle variables for coupled oscillator systems",
abstract = "A theory of action-angle variables for coupled oscillator systems is developed which involves solving the Schr{\"o}dinger equation using a basis of WKB eigenfunctions, then using the logarithm of the resulting wavefunction to define the generator for the canonical transformation which determines the action-angle variables. This theory is based on the marriage between Miller's method for solving the Hamilton-Jacobi equation using the logarithm of the generating function, and the Ratner-Buch-Gerber method for solving the Schr{\"o}dinger equation using WKB basis functions. A perturbation-theory analysis of this theory indicates that the semiclassical eigenvalues and canonical transformations obtained from it should become identical to their exact classical counterparts in the limit of large actions for each vibrational mode. Two methods for systematically improving the theory for the lower eigenstates are also proposed. Numerical applications of the theory are presented for two systems, the Morse oscillator and the Henon-Heiles two-mode hamiltonian. The resulting semiclassical eigenvalues are in excellent agreement with their exact quantum counterparts, with the magnitude of the error roughly independent of the energy of the eigenstate. Analogous good agreement is found in comparing the approximate and exact classical canonical transformations. In particular, for the Morse oscillator, good results are obtained for certain higher energy states where second-order classical perturbation theory makes serious errors. Other information examined includes surfaces of section for the Henon-Heiles system (comparing the analytical functions obtained from the present theory with results based on exact trajectory calculations) and vibrational distributions chosen to simulate trajectory calculations (using the present theory to determine bin boundaries for a histogram calculation). Again, the comparison in each case with accurate results is excellent, with maximum errors in action calculations of 0.02 h, and in angle calculations of 0.01 rad.",
author = "Geiger, {Lynn C.} and Schatz, {George C}",
year = "1983",
month = "9",
day = "15",
doi = "10.1016/0301-0104(83)85266-5",
language = "English",
volume = "79",
pages = "431--447",
journal = "Chemical Physics",
issn = "0301-0104",
publisher = "Elsevier",
number = "3",

}

TY - JOUR

T1 - A new approach to the determination of good action-angle variables for coupled oscillator systems

AU - Geiger, Lynn C.

AU - Schatz, George C

PY - 1983/9/15

Y1 - 1983/9/15

N2 - A theory of action-angle variables for coupled oscillator systems is developed which involves solving the Schrödinger equation using a basis of WKB eigenfunctions, then using the logarithm of the resulting wavefunction to define the generator for the canonical transformation which determines the action-angle variables. This theory is based on the marriage between Miller's method for solving the Hamilton-Jacobi equation using the logarithm of the generating function, and the Ratner-Buch-Gerber method for solving the Schrödinger equation using WKB basis functions. A perturbation-theory analysis of this theory indicates that the semiclassical eigenvalues and canonical transformations obtained from it should become identical to their exact classical counterparts in the limit of large actions for each vibrational mode. Two methods for systematically improving the theory for the lower eigenstates are also proposed. Numerical applications of the theory are presented for two systems, the Morse oscillator and the Henon-Heiles two-mode hamiltonian. The resulting semiclassical eigenvalues are in excellent agreement with their exact quantum counterparts, with the magnitude of the error roughly independent of the energy of the eigenstate. Analogous good agreement is found in comparing the approximate and exact classical canonical transformations. In particular, for the Morse oscillator, good results are obtained for certain higher energy states where second-order classical perturbation theory makes serious errors. Other information examined includes surfaces of section for the Henon-Heiles system (comparing the analytical functions obtained from the present theory with results based on exact trajectory calculations) and vibrational distributions chosen to simulate trajectory calculations (using the present theory to determine bin boundaries for a histogram calculation). Again, the comparison in each case with accurate results is excellent, with maximum errors in action calculations of 0.02 h, and in angle calculations of 0.01 rad.

AB - A theory of action-angle variables for coupled oscillator systems is developed which involves solving the Schrödinger equation using a basis of WKB eigenfunctions, then using the logarithm of the resulting wavefunction to define the generator for the canonical transformation which determines the action-angle variables. This theory is based on the marriage between Miller's method for solving the Hamilton-Jacobi equation using the logarithm of the generating function, and the Ratner-Buch-Gerber method for solving the Schrödinger equation using WKB basis functions. A perturbation-theory analysis of this theory indicates that the semiclassical eigenvalues and canonical transformations obtained from it should become identical to their exact classical counterparts in the limit of large actions for each vibrational mode. Two methods for systematically improving the theory for the lower eigenstates are also proposed. Numerical applications of the theory are presented for two systems, the Morse oscillator and the Henon-Heiles two-mode hamiltonian. The resulting semiclassical eigenvalues are in excellent agreement with their exact quantum counterparts, with the magnitude of the error roughly independent of the energy of the eigenstate. Analogous good agreement is found in comparing the approximate and exact classical canonical transformations. In particular, for the Morse oscillator, good results are obtained for certain higher energy states where second-order classical perturbation theory makes serious errors. Other information examined includes surfaces of section for the Henon-Heiles system (comparing the analytical functions obtained from the present theory with results based on exact trajectory calculations) and vibrational distributions chosen to simulate trajectory calculations (using the present theory to determine bin boundaries for a histogram calculation). Again, the comparison in each case with accurate results is excellent, with maximum errors in action calculations of 0.02 h, and in angle calculations of 0.01 rad.

UR - http://www.scopus.com/inward/record.url?scp=48749143614&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=48749143614&partnerID=8YFLogxK

U2 - 10.1016/0301-0104(83)85266-5

DO - 10.1016/0301-0104(83)85266-5

M3 - Article

AN - SCOPUS:48749143614

VL - 79

SP - 431

EP - 447

JO - Chemical Physics

JF - Chemical Physics

SN - 0301-0104

IS - 3

ER -