Abstract
The construction of a statistical model for eigenfunctions of the Ising model in transverse and longitudinal fields is discussed in detail for the chaotic case. When the number of spins is large, each wave function coefficient has the Gaussian distribution with zero mean and variance calculated from the first two moments of the Hamiltonian. The main part of the paper is devoted to the discussion of various corrections to the asymptotic result. One type of correction is related to higher order moments of the Hamiltonian, and can be taken into account by Gibbs-like formulae. Other corrections are due to symmetry contributions, which manifest as different numbers of non-zero real and complex coefficients. The statistical model with these corrections included agrees well with numerical calculations of wave function moments.