It is shown that the Husimi representations of chaotic eigenstates are strongly correlated along classical trajectories. These correlations extend across the whole system size and, unlike the corresponding eigenfunction correlations in configuration space, they persist in the semiclassical limit. A quantitative theory is developed on the basis of Gaussian wavepacket dynamics and random-matrix arguments. The role of symmetries is discussed for the example of time-reversal invariance.PACS numbers: 03.65. Sq, 05.45.Mt Chaotic eigenfunctions and in particular their localization and correlation statistics are a topic of continuing interest [1,2,3,4,5,6,7]. Applications include classical, mesoscopic and pure quantum systems such as optical, mechanical and microwave resonators [8,9,10], electron transmission and interaction in chaotic quantum dots [11,12], and decay and fluctuations of heavy nuclei [13]. One foundation of eigenfunction statistics is the randomwave model of Berry [1] which is essentially equivalent to random-matrix theory (RMT) [2]. Within RMT the eigenfunction components in an arbitrary basis are uncorrelated Gaussians. Current research is frequently aiming at deviations from RMT due to the specific dynamics. Prominent examples are scarring by periodic orbits [3] or long-range correlations [6]. As classical dynamics takes place in phase space, representations of eigenstates via Husimi or Wigner functions seem appropriate, and recently some of their statistical properties have attracted a lot of attention [7]. Surprisingly this does not apply to dynamically induced correlations although numerous studies of the corresponding spatial correlations demonstrate their relevance [4,5,6], and although there are systems where a direct relation between phase-space correlations and measurable quantities must be expected. For example, in optical resonators [8] the power emitted at a certain point of the boundary depends strongly on the angle of incidence of the wave (total internal reflection). Also in quantum dots effects of eigenfunction directionality have been measured using tilted leads [12].In this paper we analyze for the first time dynamical correlations between points in phase space. Our results are surprising in view of the fact that spatial two-point correlators of eigenfunction amplitudes or densities vanish in the semiclassical limith → 0 for any x = x ′ [14]. In contrast we find strong and semiclassically persistent correlations between phase-space points ξ = ξ ′ (ξ = (x, p)) where the distance between x and x ′ can be of the order of the system size. This is no contradiction since ψ(x)ψ(x ′ ) = 0 does not imply that ψ(x) and ψ(x ′ ) are statistically independent; it just means that there are no linear correlations in configuration space. In other words, although the existence of dynamically induced correla- tions cannot be a matter of the chosen basis they turn out to be most relevant in phase space.In our approach we make use of methods which proved successful in studies of eigenfunction scarring...