Abstract-Several electrically large Phased Array Feed (PAF) reflector systems are modeled to examine the mechanism of multiple reflections between parabolic reflectors and low-and highscattering feeds giving rise to frequency-dependent patterns and impedance ripples. The PAF current is expanded in physics-based macro domain basis functions (CBFs), while the reflector employs the Physical Optics (PO) equivalent current. The reflector-feed coupling is systematically accounted for through a multiscattering Jacobi approach. An FFT expands the reflector radiated field in only a few plane waves, and the reflector PO current is computed rapidly through a near-field interpolation technique. The FEKO software is used for several cross validations, and the convergence properties of the hybrid method are studied for several representative examples showing excellent numerical performance. The measured and simulated results for a 121-element Vivaldi PAF, which is installed on the Westerbork Synthesis Radio Telescope, are in very good agreement.Index Terms-phased array feeds, radio astronomy, method of moments, characteristic basis function method, physical optics.