# Quantum chaos

Quantum chaos is a branch of physics focused on how chaotic classical dynamical systems can be described in terms of quantum theory. The primary question that quantum chaos seeks to answer is: "What is the relationship between quantum mechanics and classical chaos?" The correspondence principle states that classical mechanics is the classical limit of quantum mechanics, specifically in the limit as the ratio of the Planck constant to the action of the system tends to zero. If this is true, then there must be quantum mechanisms underlying classical chaos (although this may not be a fruitful way of examining classical chaos). If quantum mechanics does not demonstrate an exponential sensitivity to initial conditions, how can exponential sensitivity to initial conditions arise in classical chaos, which must be the correspondence principle limit of quantum mechanics?
In seeking to address the basic question of quantum chaos, several approaches have been employed:
* Development of methods for solving quantum problems where the perturbation cannot be considered small in perturbation theory and where quantum numbers are large.
* Correlating statistical descriptions of eigenvalues (energy levels) with the classical behavior of the same Hamiltonian (system).
* Study of probability distribution of individual eigenstates (see scars and quantum ergodicity).
* Semiclassical methods such as periodic-orbit theory connecting the classical trajectories of the dynamical system with quantum features.
* Direct application of the correspondence principle.

## History

During the first half of the twentieth century, chaotic behavior in mechanics was recognized (as in the three-body problem in celestial mechanics), but not well understood. The foundations of modern quantum mechanics were laid in that period, essentially leaving aside the issue of the quantum-classical correspondence in systems whose classical limit exhibit chaos.

## Approaches

Questions related to the correspondence principle arise in many different branches of physics, ranging from nuclear to atomic, molecular and solid-state physics, and even to acoustics, microwaves and optics. However, classical-quantum correspondence in chaos theory is not always possible. Thus, some versions of the classical butterfly effect do not have counterparts in quantum mechanics.
Important observations often associated with classically chaotic quantum systems are spectral level repulsion, dynamical localization in time evolution (e.g. ionization rates of atoms), and enhanced stationary wave intensities in regions of space where classical dynamics exhibits only unstable trajectories (as in scattering). In the semiclassical approach of quantum chaos, phenomena are identified in spectroscopy by analyzing the statistical distribution of spectral lines and by connecting spectral periodicities with classical orbits. Other phenomena show up in the time evolution of a quantum system, or in its response to various types of external forces. In some contexts, such as acoustics or microwaves, wave patterns are directly observable and exhibit irregular amplitude distributions.
Quantum chaos typically deals with systems whose properties need to be calculated using either numerical techniques or approximation schemes (see e.g. Dyson series). Simple and exact solutions are precluded by the fact that the system's constituents either influence each other in a complex way, or depend on temporally varying external forces.

## Quantum mechanics in non-perturbative regimes

For conservative systems, the goal of quantum mechanics in non-perturbative regimes is to find the eigenvalues and eigenvectors of a Hamiltonian of the form
    $H=H_{s}+\varepsilon H_{ns},\,$
where $H_{s}$ is separable in some coordinate system, $H_{ns}$ is non-separable in the coordinate system in which $H_{s}$ is separated, and $\epsilon$ is a parameter which cannot be considered small. Physicists have historically approached problems of this nature by trying to find the coordinate system in which the non-separable Hamiltonian is smallest and then treating the non-separable Hamiltonian as a perturbation.
Finding constants of motion so that this separation can be performed can be a difficult (sometimes impossible) analytical task. Solving the classical problem can give valuable insight into solving the quantum problem. If there are regular classical solutions of the same Hamiltonian, then there are (at least) approximate constants of motion, and by solving the classical problem, we gain clues how to find them.
Other approaches have been developed in recent years. One is to express the Hamiltonian in different coordinate systems in different regions of space, minimizing the non-separable part of the Hamiltonian in each region. Wavefunctions are obtained in these regions, and eigenvalues are obtained by matching boundary conditions.
Another approach is numerical matrix diagonalization. If the Hamiltonian matrix is computed in any complete basis, eigenvalues and eigenvectors are obtained by diagonalizing the matrix. However, all complete basis sets are infinite, and we need to truncate the basis and still obtain accurate results. These techniques boil down to choosing a truncated basis from which accurate wavefunctions can be constructed. The computational time required to diagonalize a matrix scales as $N^{3}$, where $N$ is the dimension of the matrix, so it is important to choose the smallest basis possible from which the relevant wavefunctions can be constructed. It is also convenient to choose a basis in which the matrix is sparse and/or the matrix elements are given by simple algebraic expressions because computing matrix elements can also be a computational burden.
A given Hamiltonian shares the same constants of motion for both classical and quantum dynamics. Quantum systems can also have additional quantum numbers corresponding to discrete symmetries (such as parity conservation from reflection symmetry). However, if we merely find quantum solutions of a Hamiltonian which is not approachable by perturbation theory, we may learn a great deal about quantum solutions, but we have learned little about quantum chaos. Nevertheless, learning how to solve such quantum problems is an important part of answering the question of quantum chaos.

## Correlating statistical descriptions of quantum mechanics with classical behaviour

Statistical measures of quantum chaos were born out of a desire to quantify spectral features of complex systems. Random matrix theory was developed in an attempt to characterize spectra of complex nuclei. The remarkable result is that the statistical properties of many systems with unknown Hamiltonians can be predicted using random matrices of the proper symmetry class. Furthermore, random matrix theory also correctly predicts statistical properties of the eigenvalues of many chaotic systems with known Hamiltonians. This makes it useful as a tool for characterizing spectra which require large numerical efforts to compute.
A number of statistical measures are available for quantifying spectral features in a simple way. It is of great interest whether or not there are universal statistical behaviors of classically chaotic systems. The statistical tests mentioned here are universal, at least to systems with few degrees of freedom (Berry and Tabor have put forward strong arguments for a Poisson distribution in the case of regular motion and Heusler et al. present a semiclassical explanation of the so-called Bohigas–Giannoni–Schmit conjecture which asserts universality of spectral fluctuations in chaotic dynamics). The nearest-neighbor distribution (NND) of energy levels is relatively simple to interpret and it has been widely used to describe quantum chaos.
Qualitative observations of level repulsions can be quantified and related to the classical dynamics using the NND, which is believed to be an important signature of classical dynamics in quantum systems. It is thought that regular classical dynamics is manifested by a Poisson distribution of energy levels:
    $P(s)=e^{-s}.\$
In addition, systems which display chaotic classical motion are expected to be characterized by the statistics of random matrix eigenvalue ensembles. For systems invariant under time reversal, the energy-level statistics of a number of chaotic systems have been shown to be in good agreement with the predictions of the Gaussian orthogonal ensemble (GOE) of random matrices, and it has been suggested that this phenomenon is generic for all chaotic systems with this symmetry. If the normalized spacing between two energy levels is $s$, the normalized distribution of spacings is well approximated by
    $P(s)={\frac {\pi }{2}}se^{-\pi s^{2}/4}.$
Many Hamiltonian systems which are classically integrable (non-chaotic) have been found to have quantum solutions that yield nearest neighbor distributions which follow the Poisson distributions. Similarly, many systems which exhibit classical chaos have been found with quantum solutions yielding a Wigner-Dyson distribution, thus supporting the ideas above. One notable exception is diamagnetic lithium which, though exhibiting classical chaos, demonstrates Wigner (chaotic) statistics for the even-parity energy levels and nearly Poisson (regular) statistics for the odd-parity energy level distribution.

## Semiclassical methods
### Periodic orbit theory

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:
    $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)}.$
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. The index $k$ distinguishes the primitive periodic orbits: the shortest period orbits of a given set of initial conditions. $T_{k}$ is the period of the primitive periodic orbit and $S_{k}$ is its classical action. Each primitive orbit retraces itself, leading to a new orbit with action $nS_{k}$ and a period which is an integral multiple $n$ 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 $n$. $\alpha _{nk}$ is the orbit's Maslov index. The amplitude factor, $1/\sinh {(\chi _{nk}/2)}$, 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 $\chi _{nk}$ 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, $\sinh {(\chi _{nk}/2)}$ becomes $\sin {(\chi _{nk}/2)}$, where $\chi _{nk}$ is the winding number of the periodic orbit. $\chi _{nk}=2\pi m$, where $m$ is the number of times that neighboring orbits intersect the periodic orbit in one period. This presents a difficulty because $\sin {(\chi _{nk}/2)}=0$ at a classical 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 $1/r$ potential with an anisotropic mass tensor) semiclassically. He found agreement with quantum computations for low lying (up to $n=6$) 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 $d(E)=-{\frac {1}{\pi }}\Im (\operatorname {Tr} (G(x,x^{\prime },E))$ (tracing over positions) and calculate it again in stationary phase approximation to get an approximation for the density of states $d(E)$.
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

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 $\epsilon$ is given by a smooth background plus an oscillatory sum of the form
    $f(w)=\sum _{k}\sum _{n=1}^{\infty }D_{\it {nk}}^{i}\sin(2\pi nw{\tilde {S_{k}}}-\phi _{\it {nk}}).$
$\phi _{\it {nk}}$ is a phase that depends on the Maslov index and other details of the orbits. $D_{\it {nk}}^{i}$ is the recurrence amplitude of a closed orbit for a given initial state (labeled $i$). 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 $\epsilon$ as a function of $w$ is called a recurrence spectrum, because it gives peaks which correspond to the scaled action of closed orbits and whose heights correspond to $D_{\it {nk}}^{i}$.
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 $H^{-}$ 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 $y(0)=0$ the density of states obtained from the Gutzwiller formula is related to the inverse of the potential of the classical system by ${\frac {d^{1/2}}{dx^{1/2}}}V^{-1}(x)=2{\sqrt {\pi }}{\frac {dN(x)}{dx}}$ here ${\frac {dN(x)}{dx}}$ 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.

## Recent directions

One open question remains understanding quantum chaos in systems that have finite-dimensional local Hilbert spaces for which standard semiclassical limits do not apply. Recent works allowed for studying analytically such quantum many-body systems.
The traditional topics in quantum chaos concerns spectral statistics (universal and non-universal features), and the study of eigenfunctions of various chaotic Hamiltonian. For example, before the existence of scars was reported, eigenstates of a classically chaotic system were conjectured to fill the available phase space evenly, up to random fluctuations and energy conservation (Quantum ergodicity). However, a quantum eigenstate of a classically chaotic system can be scarred: the probability density of the eigenstate is enhanced in the neighborhood of a periodic orbit, above the classical, statistically expected density along the orbit (scars). In particular, scars are both a striking visual example of classical-quantum correspondence away from the usual classical limit, and a useful example of a quantum suppression of chaos. For example, this is evident in the perturbation-induced quantum scarring: More specifically, in quantum dots perturbed by local potential bumps (impurities), some of the eigenstates are strongly scarred along periodic orbits of unperturbed classical counterpart.
Further studies concern the parametric ($R$) dependence of the Hamiltonian, as reflected in e.g. the statistics of avoided crossings, and the associated mixing as reflected in the (parametric) local density of states (LDOS). There is vast literature on wavepacket dynamics, including the study of fluctuations, recurrences, quantum irreversibility issues etc. Special place is reserved to the study of the dynamics of quantized maps: the standard map and the kicked rotator are considered to be prototype problems.
Works are also focused in the study of driven chaotic systems, where the Hamiltonian $H(x,p;R(t))$ is time dependent, in particular in the adiabatic and in the linear response regimes. There is also significant effort focused on formulating ideas of quantum chaos for strongly-interacting many-body quantum systems far from semi-classical regimes as well as a large effort in quantum chaotic scattering.

## Berry–Tabor conjecture

In 1977, Berry and Tabor made a still open "generic" mathematical conjecture which, stated roughly, is: In the "generic" case for the quantum dynamics of a geodesic flow on a compact Riemann surface, the quantum energy eigenvalues behave like a sequence of independent random variables provided that the underlying classical dynamics is completely integrable.