LETTER TO THE EDITOR: Exact energy distribution function in a time-dependent harmonic oscillator

Following a recent work by Robnik and Romanovski (2006 J. Phys. A: Math. Gen. 39 L35, 2006 Open Syst. Inf. Dyn. 13 197-222), we derive an explicit formula for the universal distribution function of the final energies in a time-dependent 1D harmonic oscillator, whose functional form does not depend on the details of the frequency ω(t) and is closely related to the conservation of the adiabatic invariant. The normalized distribution function is P(x) = \pi^{-1} (2\mu^2 - x^2)^{-\frac{1}{2}} , where x=E_1- \skew3\bar{E}_1 ; E₁ is the final energy, \skew3\bar{E}_1 is its average value and µ² is the variance of E₁. \skew3\bar{E}_1 and µ² can be calculated exactly using the WKB approach to all orders.