Photoplethysmography (PPG) is a low-cost, non-invasive optical technology that takes physiological measurements on the surface of the skin. One of its most widespread applications is the wearable heart rate sensor included in commercially available smartwatches and sports bracelets, that provides comfortable and continuous pulse monitoring during everyday tasks. This blog post demonstrates how to model the human skin in OpticStudio for physiological measurements and illustrates a time-dependent simulation of PPG-based heart rate sensors using ZOS-API.

PPG devices consist of infrared or visible range light-emitting diodes (LEDs) and photodetectors. They provide a simple optical technique to detect blood volume changes in tissues. Light is more strongly absorbed and scattered by blood than by the surrounding tissue. Therefore, the pulsation of blood causes a variation of the opposite phase in the signal of the detector. In this blog, we will explore how to implement a human skin tissue model in OpticStudio and how to simulate the measured signal of a PPG device over time using a ZOS-API application.

Basic design

PPG sensors can be used either in reflection or in transmission mode. As the penetration depth of the light depends on its wavelength, green and yellow LEDs are most suitable to take measurements in the superficial blood flow and are typically used in reflection mode. On the other hand, infrared and near-infrared wavelengths are better suited for measurements of deep-tissue blood flow and can be used in transmission mode as well. In this example, we present a reflection PPG device.

Our goal was to develop a realistic skin model based on data published in the relevant literature. Thus, we intended to apply a wavelength at which the optical parameters of both the skin and the blood are widely available in the literature, and which is also close to the most frequently used wavelengths in commercial devices. As a result, we chose the modelling wavelength of 575 nm, and we used a QSMF-C160 LED (Avago Technologies) as the source. The model of this LED can be directly downloaded from the Radiant Source Models library, and a Source File can be created by generating rays from the Radiant Source Model file. This feature is available in the Premium edition of OpticStudio, and the process is discussed in detail in the knowledgebase article How to generate a ray set from an RSMX Source Model.

Human skin model

To simulate light transport in tissue media, we implemented a layered skin model, which accounts for the epidermis, the dermis, and the subcutaneous fat. As the primary goal for this use case was to simulate a PPG-based heart rate sensor, where the key point is to measure the changes caused by the pulsation of the blood, we focused on accurately modelling the layers where this pulsation can be observed. Therefore, we separately modelled the sub-layers of the dermis with different blood content values, namely the papillary dermis, subpapillary dermis, upper blood net dermis, reticular dermis, and deep blood net dermis. On the other hand, as there is no blood content in the epidermis - to keep the model simple - we decided to use only one thick epidermis layer, which accounts for all the stratum corneum, stratum granulosum, stratum spinosum, and stratum basale. Finally, similarly to most of the published skin models, we represented the subcutaneous fat with one single layer too.

All the above-mentioned layers are modelled as Rectangular Volumes in OpticStudio. The thickness values of the layers were based on literature data, and the cross-sectional size was determined in such a way that there is no light leakage at the sides. The subsequent layers were placed by using the previous layer as Reference Object and applying Pickup solve on the Z Position value from the Z Length column of the previous layer. This solution ensured that the layers were located closely one after the other without any gap between them.

Customizing the tissue layers

As this case study relies solely on data published in the literature, we did not carry out any new measurements throughout the research process. Although the model parameters were based on published data, it is important to note that the optical parameters of the human skin may vary significantly across the population. Thus, different parameters might be required for specific subject groups. Therefore, if more accurate data are available for your specific applications, then do not hesitate to customize the model accordingly! 

Representing all the individual blood vessels in the skin would require the addition of several hundreds of objects with the complex spatial arrangement, and would reduce the generality of the model, thus this type of modelling is not widespread in the literature. Therefore, we did not apply this technique either. Instead, we took the blood content of the different skin layers into account by calculating the weighted average of the optical parameters of the blood and the surrounding tissue.

Accordingly, we modelled the material of the skin layers using the Model material solve type in OpticStudio, based on the following raw data:

Skin layer	n @ 575 nm	Blood content [%]
living epidermis	1.45	0
papillary dermis	1.40	2
upper blood net dermis	1.40	5
reticular dermis	1.40	1
deep blood net dermis	1.40	5
subcutaneous fat	1.44	5
blood	1.35	n/a

Bulk scattering in tissue

Light scattering by small particles in turbid media, such as in biological tissues, can be accurately described by the Henyey-Greenstein distribution function. The Henyey-Greenstein model has only one free parameter, the anisotropy factor g. The domain of this parameter is the [-1, 1] interval, where g=-1 corresponds to backscattering, g=0 describes isotropic scattering, and g=1 represents forward scattering. The angular distribution of the scattered light is defined as:

In the Non-Sequential Mode of OpticStudio, the Henyey-Greenstein bulk scattering model is available in a DLL (Henyey-Greenstein-bulk.DLL), included in the OpticStudio installation files. Further explanations of the model together with comprehensive analyses of the DLL are presented in the knowledgebase article Using the Henyey-Greenstein distribution to model bulk scattering. Besides, a detailed discussion on the surface and bulk scattering models in OpticStudio can be found in this knowledgebase article: What scattering models are available in OpticStudio?

In the attached multilayer skin model, the scattering parameters of each layer were set based on realistic values presented in the literature. While the input parameters of the Henyey-Greenstein scattering DLL are the Mean Path, the Transmission fraction, and the g anisotropy parameter, in the literature, typically the scattering and absorption coefficients, µ_s and µ_a respectively, are published with the anisotropy factor. Therefore, we used the following equations to compute the input parameters of the model:

Similarly, as for the refractive index, the scattering and absorption coefficients, as well as the anisotropy factor of the different skin layers were calculated as the weighted average of the corresponding values for the blood and the rest of the tissue, based on the following raw data corresponding to the applied wavelength of 575 nm:

Skin layer	µa [1/mm]	µs [1/mm]	g [-]
living epidermis	1	20	0.79
papillary dermis	0.28	21.5	0.79
upper blood net dermis	0.28	21.5	0.79
reticular dermis	0.28	21.5	0.79
deep blood net dermis	0.28	21.5	0.79
subcutaneous fat	0.081	1.396	0.75
blood	32.4	50	0.98

The multi-layered skin model and light transport in the tissue are shown in the 3D layout plots below. To illustrate the individual scattering events, rays are coloured by segments on the plots.

To provide numerical results in addition to the illustration, we added three Detector Rectangle objects to the design. They are separated by a thin air gap from the surface of the skin. Two detectors have the same cross-sectional size as the skin layers, one facing towards the source and the other towards the skin model, to measure all the incoming and back-scattered light respectively for references. The third detector is a small one (2 mm by 2 mm), again facing towards the skin, which represents a typical photodetector in PPG devices.

The above-presented design can be used as an off-the-shelf skin model when the time-dependence of the measurements/simulations is not relevant.

To access the simulation files and learn about modeling time-dependent effects using ZOS-API, such as how to mimic pulsating blood flow in the tissue for heart rate monitoring, read the full article here.

To try out Zemax capabilities please download a free trial here.

Author:

Csilla Timar-Fulep
Senior Application Engineer
Zemax