Parabolic equation numerical simulations of waveguide acoustical beam propagation in an ocean of Garrett–Munk internal waves are used to examine the range evolution of beam properties such as beamwidth (both spectral and spatial), Shannon entropy, and scintillation index, as a function of beam angle. Simulations are carried out at 250- and 125-Hz acoustic frequencies. The ray trajectories associated with these beams are predominantly chaotic or exponentially sensitive to initial conditions and/or medium perturbations. At long range near saturation, the finite-frequency beams show a constant rate of change of Shannon entropy with range, independent of acoustic frequency. This full-wave rate of entropy is of the same order of magnitude as the average rate of entropy for the ray trajectories associated with this beam. Finite-range Lyapunov exponents provide the estimates of ray entropy rate or Kolmogorov–Siani entropy. The correspondence between full-wave and ray entropies suggests a full-wave manifestation of ray chaos, but only once statistical saturation is obtained. In spite of this correspondence, the simulated acoustical beams expand diffusively not exponentially (or explosively).