Reconstructing in vivo spatially offset Raman spectroscopy of human skin tissue using a GPU-accelerated Monte Carlo platform

As one type of spatially offset Raman spectroscopy (SORS), inverse SORS is particularly suited to in vivo biomedical measurements due to its ring-shaped illumination scheme. To explain inhomogeneous Raman scattering during in vivo inverse SORS measurements, the light–tissue interactions when excitation and regenerated Raman photons propagate in skin tissue were studied using Monte Carlo simulation. An eight-layered skin model was first built based on the latest transmission parameters. Then, an open-source platform, Monte Carlo eXtreme (MCX), was adapted to study the distribution of 785 nm excitation photons inside the model with an inverse spatially shifted annular beam. The excitation photons were converted to emission photons by an inverse distribution method based on excitation flux with spatial offsets Δ s of 1 mm, 2 mm, 3 mm and 5 mm. The intrinsic Raman spectra from separated skin layers were measured by continuous linear scanning to improve the simulation accuracy. The obtained results explain why the spectral detection depth gradually increases with increasing spatial offset, and address how the intrinsic Raman spectrum from deep skin layers is distorted by the reabsorption and scattering of the superficial tissue constituents. Meanwhile, it is demonstrated that the spectral contribution from subcutaneous fat will be improved when the offset increases to 5 mm, and the highest detection efficiency for dermal layer spectral detection could be achieved when Δ s = 2 mm. Reasonably good matching between the calculated spectrum and the measured in vivo inverse SORS was achieved, thus demonstrating great utility of our modeling method and an approach to help understand the clinical measurements.