“…[19][20][21] If one wishes to treat accurately the GHz and THz regions, the complete phase-space distribution W͑ , p , t͒ ͑p = m˙͒ must be used, giving rise, [15][16][17][18]30 on expansion of the momentum part of the distribution in orthogonal Hermite polynomials, to a hierarchy of partial differential recurrence relations in configuration space. The actual configuration-space distribution P͑ , t͒ must then be extracted from the hierarchy, usually by continued fractions, with the much simpler QSE naturally emerging 15,17 from the hierarchy in the high-damping limit …”