Determination of Quantum Yield in Scattering Media Using Monte Carlo Photoluminescence Cascade Simulation and Integrating Sphere Measurements

Accurate determination of the quantum yield (Φf) in scattering media is essential for numerous scientific and industrial applications, but it remains challenging due to re-absorption and scattering-induced biases. In this study, we present a GPU-accelerated Monte Carlo simulation framework that solves the full fluorescence radiative transfer equation (FRTE), incorporating spectrally dependent absorption, scattering, and fluorescence cascade processes. The model accounts for re-emission shifts, energy scaling due to the Stokes shift and implements a digital optical twin of the experimental setup, including the precise description of the applied integrating sphere. Using Rhodamine 6G in both ethanol and PDMS matrices, we demonstrate the accuracy of the method by comparing simulated reflectance and transmission spectra with independent experimental measurements. Φf and emission distributions are optimized using a Levenberg-Marquardt algorithm. The obtained quantum yields agree well with literature values for Rhodamine 6G. This approach eliminates the need for empirical correction factors, enabling the reliable determination of actual, undistorted emission spectra and the Φf in complex scattering media.