Nonlinear dynamics in photon capture and uptake at the photosystem level may have a strong effect on photosynthetic yield. However, the magnitude of such effects is difficult to estimate theoretically because nonlinear systems often cannot be represented accurately using equations. A nonanalytical simulation was developed that used a simple decision tree and Monte Carlo methods, instead of equations, to model how a population of photosystems absorbs and utilizes photons from an ambient light field. This simulation replicated realistic kinetics in the closure and variable fluorescence yield of PSII on the single-turnover timescale, as well as the saturating behavior in light-driven electron flow that is observed in nature with increasing irradiance. This simulation indicated that the transfer of absorbed photon energy among PSII units can introduce strong nonlinear enhancement in light-driven electron flow. However, this effect was seen only in populations with particular photosynthetic states as determined by physiological properties of PSII. Other populations with the same degree of energy transfer but with different photosynthetic states exhibited little enhancement in electron flow and, in some cases, a reduction. This nonanalytical approach provides a simple means to quantify theoretically how nonlinear dynamics in photosynthesis arise at the photosystem level and how these dynamics may act to enhance or constrain photosynthetic rates and yields. Such simulations can provide quantitative insight into different physiological bases of nonlinear light-harvesting dynamics and identify those that would have the strongest theoretical influence and thus warrant closer experimental examination.