Modern optical technologies encompass classical light phenomena and non-linear effects, crucial for biomedical imaging and therapies. Despite substantial interest and many experimental studies, non-linear optical effects in inhomogeneous scattering media require extensive computer modeling to understand light-material interactions. This report focuses on Monte Carlo simulations of Raman scattering in turbid media, such as biological tissues, essential for clinical applications. We developed a novel, open-source, power-efficient Monte Carlo algorithm optimized for Apple's M-series chips. By incorporating single-photon tracking, the algorithm accurately models time-resolved light-tissue interactions. Validated against experiments, it matches traditional methods' accuracy while offering superior efficiency, advancing Raman functional imaging.