Approach to ergodicity in quantum wave functions

According to theorems of Shnirelman and followers, in the semiclassical limit the quantum wavefunctions of classically ergodic systems tend to the microcanonical density on the energy shell. We here develop a semiclassical theory that relates the rate of approach to the decay of certain classical fluctuations. For uniformly hyperbolic systems we find that the variance of the quantum matrix elements is proportional to the variance of the integral of the associated classical operator over trajectory segments of length $T_H$, and inversely proportional to $T_H^2$, where $T_H=h\bar\rho$ is the Heisenberg time, $\bar\rho$ being the mean density of states. Since for these systems the classical variance increases linearly with $T_H$, the variance of the matrix elements decays like $1/T_H$. For non-hyperbolic systems, like Hamiltonians with a mixed phase space and the stadium billiard, our results predict a slower decay due to sticking in marginally unstable regions. Numerical computations supporting these conclusions are presented for the bakers map and the hydrogen atom in a magnetic field.