The accurate estimation of skin and skull optical properties over a wide wavelength range of laser radiation has great importance in optogenetics and other related applications. In the present work, using the Kubelka-Munk model, finite-element solution of the diffusion equation, inverse adding-doubling (IAD), and Monte-Carlo simulation, we estimated the refractive index, absorption and scattering coefficients, penetration depth, and the optical fluence distribution in rabbit head tissues ex vivo, after dividing the heads into three types of tissues with an average thickness of skin of 1.1 mm, skull of 1 mm, and brain of 3 mm. The total diffuse reflectance and transmittance were measured using a single integrating sphere optical setup for laser radiation of 532, 660, 785, and 980 nm. The calculated optical properties were then applied to the diffusion equation to compute the optical fluence rate distribution at the boundary of the samples using the finite element method. Monte-Carlo simulation was implemented for estimating the optical fluence distribution through a model containing the three tissue layers. The scattering coefficient decreased at longer wavelengths, leading to an increase in optical fluence inside the tissue samples, indicating a higher penetration depth, especially at 980 nm. In general, the obtained results show good agreement with relevant literature.
Keywords: Inverse adding-doubling (IAD); Kubelka–Munk model; Monte-Carlo simulation; brain; diffuse reflectance; finite-element solution; optical fluence rate; optical transmittance; tissue optical properties; transcranial laser irradiation.