Compact model for light propagation simulation in fNIRS systems.

RTE-Sim is an open-source compact model designed for light propagation simulation in functional near-infrared spectroscopy (fNIRS) systems. The goals for this compact model are:

• The model offers a more efficient alternative to traditional numerical methods like the Finite Element Method (FEM) or Monte Carlo (MC) simulation for modeling light propagation in tissues.

• The model is compatible with mainstream circuit design tools, making it valuable for achieving comprehensive simulations of the fNIRS system's signal chain.

The model is constructed based on the analytical solution of the radiative transport equation (RTE) for uniform infinite media, where the RTE is expressed as:

$$ \frac{1}{c_\text{m}}\cdot\frac{\partial \Phi(\vec{r},t)}{\partial t}-\kappa\nabla^2\Phi(\vec{r},t)+\mu_a\Phi(\vec{r},t)=Q(\vec{r},t) $$

And the Green's function of the RTE $\hat{\rho}(\vec{r},\omega)$ is given by (in frequency domain):

$$ \left{\begin{array}{l} \displaystyle|\hat{\rho}(\vec{r},\omega)|=\frac{1}{4\pi\kappa r}\cdot\exp\left(-r\sqrt{\frac{\mu_a}{\kappa}}\right)\cdot\left(\frac{1}{\sigma r+\delta}\right)\\ \\ \displaystyle\angle\hat{\rho}(\vec{r},\omega)=\exp\left(-i\tau\omega\right),~\tau=\sqrt{\frac{\mu_a}{\kappa}}\cdot\frac{r}{2\mu_ac_\text{m}} \end{array}\right. $$

In this repository, we provides two implementations of this model: one in MATLAB and the other in Verilog-A.

To cite RTE-Sim in your research:

Chen W, Ma Z, Chen C, et al. An RTE-based Compact Model for Light Propagation Simulation in fNIRS Systems[C]//2023 IEEE Biomedical Circuits and Systems Conference (BioCAS). IEEE, 2023: 1-5.

The model is composed only of SIMULINK model file ( MATLAB\RTE_Sim.slx ), so the installation process for MATLAB users is as simple as downloading the code. Please ensure that you have installed SIMULINK before running the model.

Get the latest release of RTE-Sim (MATLAB)

The model constructed in Verilog-A is compatible with Cadence Virtuoso's simulation environment, so users of Verilog-A model can download the model files directly and import them into the library of Cadence Virtuoso.

Moreover, users interested in the Verilog-A code can just download the specific code file (you can find it in Verilog-A\<module_name>\veriloga\veriloga.va ) and independently build the model.

Get the latest release of RTE-Sim (Verilog-A)

The model can be used in the same way as a ordinary SIMULINK model.

The figure above illustrates the simulation of the light propagation process for a specific brain activity by configuring the parameters accordingly. Below is a list of some parameters that require configuration:

The model in Verilog-A comprises four sub-modules:

The model accepts a quantized voltage source as input, and the model's parameters are similar to the SIMULINK model which can be easily configured through Cadence Virtuoso's GUI:

@Wenyao Chen.

This project exists thanks to all the people who contribute.

MIT © Wenyao Chen