Open main menu
Home
Random
Recent changes
Special pages
Community portal
Preferences
About Wikipedia
Disclaimers
Incubator escapee wiki
Search
User menu
Talk
Dark mode
Contributions
Create account
Log in
Editing
Quantum chaos
(section)
Warning:
You are not logged in. Your IP address will be publicly visible if you make any edits. If you
log in
or
create an account
, your edits will be attributed to your username, along with other benefits.
Anti-spam check. Do
not
fill this in!
== Semiclassical methods == === Periodic orbit theory === [[File:Bhevendk1.jpg|right|thumb|Even parity recurrence spectrum ([[Fourier transform]] of the [[density of states]]) of diamagnetic hydrogen showing peaks corresponding to periodic orbits of the classical system. Spectrum is at a scaled energy of −0.6. Peaks labeled R and V are repetitions of the closed orbit perpendicular and parallel to the field, respectively. Peaks labeled O correspond to the near circular periodic orbit that goes around the nucleus.]] [[File:Bhevendk2.jpg|350px|right|thumb|Relative recurrence amplitudes of even and odd recurrences of the near circular orbit. Diamonds and plus signs are for odd and even quarter periods, respectively. Solid line is A/cosh(nX/8). Dashed line is A/sinh(nX/8) where A = 14.75 and X = 1.18.]] Periodic-orbit theory gives a recipe for computing spectra from the periodic orbits of a system. In contrast to the [[Einstein–Brillouin–Keller method]] of action quantization, which applies only to integrable or near-integrable systems and computes individual eigenvalues from each trajectory, periodic-orbit theory is applicable to both integrable and non-integrable systems and asserts that each periodic orbit produces a sinusoidal fluctuation in the density of states. The principal result of this development is an expression for the density of states which is the trace of the semiclassical Green's function and is given by the Gutzwiller trace formula: : <math>g_c(E) = \sum_k T_k \sum_{n=1}^\infty \frac{1}{2\sinh{(\chi_{nk}/2)}}\, e^{i(nS_k - \alpha_{nk} \pi/2)}.</math> Recently there was a generalization of this formula for arbitrary matrix Hamiltonians that involves a [[Berry phase]]-like term stemming from spin or other internal degrees of freedom.<ref>{{Cite journal|last1=Vogl|first1=M.|last2=Pankratov|first2=O.|last3=Shallcross|first3=S.|date=2017-07-27|title=Semiclassics for matrix Hamiltonians: The Gutzwiller trace formula with applications to graphene-type systems|url=https://link.aps.org/doi/10.1103/PhysRevB.96.035442|journal=Physical Review B|volume=96|issue=3|pages=035442|doi=10.1103/PhysRevB.96.035442|arxiv=1611.08879|bibcode=2017PhRvB..96c5442V|s2cid=119028983 }}</ref> The index <math>k</math> distinguishes the primitive [[periodic orbit]]s: the shortest period orbits of a given set of initial conditions. <math>T_k</math> is the period of the primitive periodic orbit and <math>S_k</math> is its classical action. Each primitive orbit retraces itself, leading to a new orbit with action <math>nS_k</math> and a period which is an integral multiple <math>n</math> of the primitive period. Hence, every repetition of a periodic orbit is another periodic orbit. These repetitions are separately classified by the intermediate sum over the indices <math>n</math>. <math>\alpha_{nk}</math> is the orbit's [[Maslov index]]. The amplitude factor, <math>1/\sinh{(\chi_{nk}/2)}</math>, represents the square root of the density of neighboring orbits. Neighboring trajectories of an unstable periodic orbit diverge exponentially in time from the periodic orbit. The quantity <math>\chi_{nk}</math> characterizes the instability of the orbit. A stable orbit moves on a [[torus]] in phase space, and neighboring trajectories wind around it. For stable orbits, <math>\sinh{(\chi_{nk}/2)}</math> becomes <math>\sin{(\chi_{nk}/2)}</math>, where <math>\chi_{nk}</math> is the winding number of the periodic orbit. <math>\chi_{nk} = 2\pi m</math>, where <math>m</math> is the number of times that neighboring orbits intersect the periodic orbit in one period. This presents a difficulty because <math>\sin{(\chi_{nk}/2)} = 0</math> at a classical [[Bifurcation theory|bifurcation]]. This causes that orbit's contribution to the energy density to diverge. This also occurs in the context of photo-[[absorption spectrum]]. Using the trace formula to compute a spectrum requires summing over all of the periodic orbits of a system. This presents several difficulties for chaotic systems: 1) The number of periodic orbits proliferates exponentially as a function of action. 2) There are an infinite number of periodic orbits, and the convergence properties of periodic-orbit theory are unknown. This difficulty is also present when applying periodic-orbit theory to regular systems. 3) Long-period orbits are difficult to compute because most trajectories are unstable and sensitive to roundoff errors and details of the numerical integration. Gutzwiller applied the trace formula to approach the [[anisotropic]] [[Kepler]] problem (a single particle in a <math>1/r</math> potential with an anisotropic mass [[tensor]]) semiclassically. He found agreement with quantum computations for low lying (up to <math>n = 6</math>) states for small anisotropies by using only a small set of easily computed periodic orbits, but the agreement was poor for large anisotropies. The figures above use an inverted approach to testing periodic-orbit theory. The trace formula asserts that each periodic orbit contributes a sinusoidal term to the spectrum. Rather than dealing with the computational difficulties surrounding long-period orbits to try to find the density of states (energy levels), one can use standard quantum mechanical perturbation theory to compute eigenvalues (energy levels) and use the Fourier transform to look for the periodic modulations of the spectrum which are the signature of periodic orbits. Interpreting the spectrum then amounts to finding the orbits which correspond to peaks in the Fourier transform. ==== Rough sketch on how to arrive at the Gutzwiller trace formula ==== # Start with the semiclassical approximation of the time-dependent Green's function (the Van Vleck propagator). # Realize that for caustics the description diverges and use the insight by Maslov (approximately Fourier transforming to momentum space (stationary phase approximation with h a small parameter) to avoid such points and afterwards transforming back to position space can cure such a divergence, however gives a phase factor). # Transform the Greens function to energy space to get the energy dependent Greens function (again approximate Fourier transform using the stationary phase approximation). New divergences might pop up that need to be cured using the same method as step 3 # Use <math>d(E)=-\frac{1}{\pi}\Im(\operatorname{Tr}(G(x,x^\prime,E))</math> (tracing over positions) and calculate it again in stationary phase approximation to get an approximation for the density of states <math>d(E)</math>. Note: Taking the trace tells you that only closed orbits contribute, the stationary phase approximation gives you restrictive conditions each time you make it. In step 4 it restricts you to orbits where initial and final momentum are the same i.e. periodic orbits. Often it is nice to choose a coordinate system parallel to the direction of movement, as it is done in many books. === Closed orbit theory === [[File:Cotcomp.jpg|right|thumb|Experimental recurrence spectrum (circles) is compared with the results of the closed orbit theory of John Delos and Jing Gao for lithium [[Rydberg atom]]s in an electric field. The peaks labeled 1–5 are repetitions of the electron orbit parallel to the field going from the nucleus to the classical turning point in the uphill direction.]] Closed-orbit theory was developed by J.B. Delos, M.L. Du, J. Gao, and J. Shaw. It is similar to periodic-orbit theory, except that closed-orbit theory is applicable only to atomic and molecular spectra and yields the oscillator strength density (observable photo-absorption spectrum) from a specified initial state whereas periodic-orbit theory yields the density of states. Only orbits that begin and end at the nucleus are important in closed-orbit theory. Physically, these are associated with the outgoing waves that are generated when a tightly bound electron is excited to a high-lying state. For [[Rydberg atoms]] and molecules, every orbit which is closed at the nucleus is also a periodic orbit whose period is equal to either the closure time or twice the closure time. According to closed-orbit theory, the average oscillator strength density at constant <math>\epsilon</math> is given by a smooth background plus an oscillatory sum of the form : <math> f(w) = \sum_k \sum_{n=1}^{\infty} D^{i}_{\it nk} \sin(2\pi nw\tilde{S_k} - \phi_{\it nk}). </math> <math>\phi_{\it nk}</math> is a phase that depends on the Maslov index and other details of the orbits. <math>D^i_{\it nk}</math> is the recurrence amplitude of a closed orbit for a given initial state (labeled <math>i</math>). It contains information about the stability of the orbit, its initial and final directions, and the matrix element of the dipole operator between the initial state and a zero-energy Coulomb wave. For scaling systems such as [[Rydberg atoms]] in strong fields, the [[Fourier transform]] of an oscillator strength spectrum computed at fixed <math>\epsilon</math> as a function of <math>w</math> is called a recurrence spectrum, because it gives peaks which correspond to the scaled action of closed orbits and whose heights correspond to <math>D^i_{\it nk}</math>. Closed-orbit theory has found broad agreement with a number of chaotic systems, including diamagnetic hydrogen, hydrogen in parallel electric and magnetic fields, diamagnetic lithium, lithium in an electric field, the <math>H^{-}</math> ion in crossed and parallel electric and magnetic fields, barium in an electric field, and helium in an electric field. === One-dimensional systems and potential === For the case of one-dimensional system with the boundary condition <math> y(0)=0 </math> the density of states obtained from the Gutzwiller formula is related to the inverse of the potential of the classical system by <math> \frac{d^{1/2}}{dx^{1/2}} V^{-1}(x)=2 \sqrt \pi \frac{dN(x)}{dx} </math> here <math> \frac{dN(x)}{dx} </math> is the density of states and V(x) is the classical potential of the particle, the [[half derivative]] of the inverse of the potential is related to the density of states as in the [[Wu–Sprung potential]].
Edit summary
(Briefly describe your changes)
By publishing changes, you agree to the
Terms of Use
, and you irrevocably agree to release your contribution under the
CC BY-SA 4.0 License
and the
GFDL
. You agree that a hyperlink or URL is sufficient attribution under the Creative Commons license.
Cancel
Editing help
(opens in new window)