Quantitative Comparison of Analytical Solution and Finite Element Method for Investigation of Near-infrared Light Propagation in Brain Tissue Model.

Introduction: Functional near-infrared spectroscopy (fNIRS) is an imaging method in which a light source and detector are installed on the head; consequently, the re-emission of light from human skin contains information about cerebral hemodynamic alteration. The spatial probability distribution profile of photons penetrating tissue at a source spot, scattering into the tissue, and being released at an appropriate detector position, represents the spatial sensitivity.

Methods: Modeling light propagation in a human head is essential for quantitative near-infrared spectroscopy and optical imaging. The specific form of the distribution of light is obtained using the theory of perturbation. An analytical solution of the perturbative diffusion equation (DE) and finite element method (FEM) in a Slab media (similar to the human head) makes it possible to study light propagation due to absorption and scattering of brain tissue.

Results: The simulation result indicates that sensitivity is slowly decreasing in the deep area, and the sensitivity below the source and detector is the highest. The depth sensitivity and computation time of both analytical and FEM methods are compared. The simulation time of the analytical approach is four times larger than the FEM.

Conclusion: In this paper, an analytical solution and the performance of FEM methods when applied to the diffusion equation for heterogeneous media with a single spherical defect are compared. The depth sensitivity along with the computation time of simulation has been investigated for both methods. For simple and Slab modes of the human brain, the analytical solution is the right candidate. Whenever the brain model is sophisticated, it is possible to use FEM methods, but it costs a higher computation time.

Highlights: Analytical and finite element method (FEM) depth sensitivity are almost the same.FEM requires more computation time, but can handle complicated head models.The analytical method is proposed for the first step and simple head models.

Plain Language Summary: The functional near-infrared spectroscopy (fNIRS) is a type of neuromonitoring that uses near-infrared light to measure brain activity indirectly and is similar to electroencephalography (EEG). A single-channel fNIRS system contains a near-infrared light source, which emits near-infrared light (NIR), and a detector is placed near the source. A light intensity change received by detectors indicates brain activity when NIR light penetrates into the gray matter. It is necessary to have a prior understanding of light penetration depth in order to measure brain activity more accurately. fNIRS can be better understood, optimized, and investigated through modeling light propagation in brain tissue. In order to study light in tissues, analytical and numerical methods can be used. In this work, we compared these two approaches quantitatively in a simple slab medium. We concluded that the numerical method takes too much time to calculate the results, but it can be applied to complicated head models. The results of these studies provide researchers with new insights into the modeling and simulation of fNIRS and diffuse optical tomography.