19 August 2017 Miniaturized pulsed laser source for time-domain diffuse optics routes to wearable devices
Author Affiliations +
J. of Biomedical Optics, 22(8), 085004 (2017). doi:10.1117/1.JBO.22.8.085004
We validate a miniaturized pulsed laser source for use in time-domain (TD) diffuse optics, following rigorous and shared protocols for performance assessment of this class of devices. This compact source (12×6  mm2) has been previously developed for range finding applications and is able to provide short, high energy (∼100  ps, ∼0.5  nJ) optical pulses at up to 1 MHz repetition rate. Here, we start with a basic level laser characterization with an analysis of suitability of this laser for the diffuse optics application. Then, we present a TD optical system using this source and its performances in both recovering optical properties of tissue-mimicking homogeneous phantoms and in detecting localized absorption perturbations. Finally, as a proof of concept of in vivo application, we demonstrate that the system is able to detect hemodynamic changes occurring in the arm of healthy volunteers during a venous occlusion. Squeezing the laser source in a small footprint removes a key technological bottleneck that has hampered so far the realization of a miniaturized TD diffuse optics system, able to compete with already assessed continuous-wave devices in terms of size and cost, but with wider performance potentialities, as demonstrated by research over the last two decades.
Di Sieno, Nissinen, Hallman, Martinenghi, Contini, Pifferi, Kostamovaara, and Mora: Miniaturized pulsed laser source for time-domain diffuse optics routes to wearable devices



Time-domain (TD) diffuse optics was introduced 30 years ago as a promising tool for deep noninvasive probing of the human body and other highly scattering media.1,2 The measurement relies on the injection of a picosecond laser pulse in the location of interest (e.g., muscle, brain, or breast); re-emitted light is then collected at a distance on the order of a few centimeters by a time-resolved single-photon detector. The distribution of photon time of flights (DTOF), obtained using time-correlated single-photon counting (TCSPC) electronics,3 encodes information about tissue absorption and scattering, which are related respectively to tissue composition and microstructure. Hence, the technique allows the characterization of tumor lesions4 or the detection of brain activation through the change of oxy- and deoxyhemoglobin concentrations,5 just to mention a few common applications.

Diffuse optics systems have already reached a commercial level with compact and wearable devices, thus being proposed not only for medical applications but also for the consumer market such as the personal monitoring of muscle oxygenation for athletes. However, these devices are always based on continuous-wave (CW)6 or frequency-domain7 operation, due to the much higher scalability of components operating in such domains. TD systems are potentially quite attractive for the high information content,8 in particular for their unique feature of selecting different mean investigated depths using different arrival times of detected photons.9,10 During the last decade, relevant improvements in terms of cost and complexity have been achieved,11,12 squeezing systems from optical tables to rack cabinets,13 down to boxes (size: 200×160×50  mm3) fitting into a backpack.14 However, CW systems are still orders of magnitude smaller and cheaper, thus dominating the market despite their intrinsic limitations; indeed, they suffer from poor information content and high sensitivity to motion artifacts. Additionally, physics itself limits the development of improved CW systems since it has been demonstrated that an ultimate ideal CW system can achieve a significantly lower penetration depth in tissues with respect to an ideal TD one (3  cm instead of 6  cm for the same paradigmatic case) and orders of magnitude lower sensitivity to localized heterogeneities.15 These ultimate performances for TD are reachable, provided that large area detectors are available in contact with the tissue under investigation (to maximize the collected light) with a proper time-gating mechanism (to reject early-arriving photons),16 and that many pulsed sources can be arranged one close to another to maximize the injected light up to the safety limits.

While miniaturized probe-hosted detectors17 and also detectors with on-chip timing circuits18 have been already well demonstrated in the literature, thus almost reaching the level of compactness required to compete with CW detection chains, pulsed laser sources still lack the desired level of integration. In a previous work by some of us,15 we presented a probe directly hosting a time-gated single-photon detector and a pulsed laser die, thus emitting all the light generated directly into the tissue, without the need of optical fibers. However, while the probe hosted the laser die and its front-end circuitry needed to bias it, it still required an external benchtop ultrafast voltage pulse generator and a radio frequency 50  Ω cable to provide electrical pulses to the probe. The transition speed and shape of such pulses are crucial to properly operate the laser diode in gain-switching mode, thus requiring the bulky and expensive radio frequency pulse generator. Apart from the overall dimension and cost, this solution is not acceptable in particular for the design of multichannel systems since it requires an external pulse generator for each laser source of the system. Indeed, the fast electrical pulse cannot be split over different paths without affecting the pulse amplitude and shape, thus impairing the laser output power.

In this work, we take advantage of a technology introduced by some of us in the field of range finding:19 a compact application-specific integrated circuit (ASIC) CMOS driver coupled to a double heterostructure laser diode operated in enhanced gain-switching mode, initially developed for time-of-flight range finding. Despite its low maximum repetition rate (1 MHz) that limits the TCSPC acquisition, diffuse optics can benefit from the high energy (0.5  nJ) and narrow (100  ps) pulses provided. Hence, the goal of this work is to demonstrate the suitability of this technology for TD measurements into highly scattering media, thus removing the last technological barrier preventing the design of wearable TD systems. Also, substantial increase in the laser repetition rate is definitely feasible, provided that the laser driver is redesigned taking into account this requirement that was not addressed in conceiving it for range finding applications.

The paper is organized as follows: Sec. 2 introduces a basic level characterization of the laser (pulse shape, average power, wavelength, and stability over time) with focus on TD diffuse optics requirements. Section 3 deals with validation on tissue-mimicking phantoms both in recovering the optical properties of homogeneous media and in detecting a localized absorbing perturbation in a scattering sample. Section 4 reports a proof of concept of in vivo use in monitoring hemodynamic changes occurring in the arm of three volunteers during a venous cuff occlusion. Finally, Sec. 5 presents conclusions and future perspectives.


Basic Characterization of the Laser


Laser Description

In Refs. 2021.22.23, it has been shown that the energy of short (100  ps) optical pulse emission in gain-switching mode can be markedly increased using a GaAs/GaAlAs quantum well laser diode with a large ratio of dact/Γ>2  μm, where dact and Γ are the active layer thickness and confinement factor, respectively. Here a compact CMOS laser driver fabricated in a 0.35-high-voltage (HV) CMOS technology was used to drive a specially grown quantum well semiconductor laser diode utilizing the enhanced gain-switching mode. The laser diode has a stripe width and length of 7  μm and 1.5 mm, respectively, and it generates optical pulses with a peak power of 0.5  W. The amplitude of the driving current pulse generated by the CMOS laser driver is adjusted with the supply voltage of the driver (VHV). The pulsing rate achieved with the CMOS driver is 1 MHz. For most of the tests presented in the following, VHV was set to 18 V, which leads to a current pulse amplitude of 1.3  A for about 1 ns. The laser driver (CMOS driver in Fig. 1) is mounted together with the semiconductor laser diode (LD in Fig. 1) and a few other passive discrete components on a small printed circuit board (12×6  mm2) as shown in Fig. 1(a). In addition, the metal can of the laser transmitter and all needed connections except ground, which is connected at the bottom of the board by using vias, are shown in Figs. 1(b) and 1(a), respectively. More details can be found in Ref. 19. To trigger the light emission at the desired rate, we generated a square wave synchronism signal with 5 V pulses lasting a few tens of nanoseconds at 1-MHz repetition rate.

Fig. 1

(a) A photograph of the laser transmitter and (b) a packaged laser transmitter with a coin for scale.



Laser Characterization

In TD diffuse optics, performances in retrieving optical properties of diffusive media and in localizing deep perturbation can be affected by the instrument response function shape, which is given by the convolution of the laser pulse shape, the detector response, and other jitter contributions.16 Its shape is usually evaluated not only in terms of full-width at half maximum (FWHM, i.e., width at 0.5 of the peak value) but also down to some orders of magnitude below the peak value (i.e., 0.1 or 0.01 of the peak value), where the effect of reflections, bumps, or decay tails can enlarge the response, thus affecting performances. A first step of this work is therefore the evaluation of the laser pulse shape at different supply voltages VHV of the laser, ranging from 17 to 35 V. The setup used was composed of the laser described earlier, whose optical pulses were collected by a 1-mm silica fiber (NA=0.39, Thorlabs GmbH, Germany) that was connected to the input of a U-bracket. Another 1-mm silica fiber exiting from the U-bracket focused the collected photons onto the active area of a photomultiplier tube (MCP-PMT with S1 photocathode, Hamamatsu Photonics). Due to a fast temporal response (<35  ps FWHM), free from slow tails, this detector guarantees the reconstruction of the generated optical pulse with negligible distortions over a wide dynamic range (4 orders of magnitude). To adjust the power injected in the collection fiber, we used a variable optical attenuator (VOA) (NDC-100C-4, Thorlabs GmbH, Germany) that was inserted in the clear path of the U-bracket. The “start” signal to the TCSPC board (SPC-130, Becker&Hickl GmbH, Germany) was fed by the PMT whereas the “stop” signal was provided by the main synchronism. For each voltage, we recorded the DTOF for 60 repetitions of 1 s each.

In Fig. 2(a), some of the recorded output pulse shapes for VHV=18, 22, 26, 30, and 34 V are reported, corresponding with current pulse amplitudes of 1.27, 1.53, 1.76, 2.00, and 2.22 A, respectively.

Fig. 2

Characterization of the laser: (a) optical pulse shapes at different VHV values, (b) temporal characteristics, (c) emission spectrum, (d) emission spectrum characteristics, (e) emitted average power at different VHV, and (f) optical pulse stability over time.


It is possible to notice that the higher the voltage (i.e., the larger the current pulse amplitude), the steeper the rising edge of the emitted laser pulse, and the wider the temporal shape of the laser.

The trend is more marked at voltages lower than 25 V. Indeed, looking at Fig. 2(b) (right axis), the center of gravity (CoG, i.e., the first moment) of the recorded curves, computed at different thresholds (0.5, 0.1, and 0.01 of the peak amplitude), after a VHV of about 25 V reaches an almost constant value (maximum deviation <2%), if data at 27 and 28 V are not taken into account. Indeed, for those VHV, an unexpected behavior, not in line with the overall trend is noticed and in such points a maximum deviation from average value (computed in all other points of the stable region, VHV>2.5  V) is 4.5%. The abnormal discontinuity might be caused by the relaxation oscillation of the semiconductor laser diode, which is a normal behavior when working in the gain-switching mode. At these specific supply voltages, the second peak of the relaxation oscillation has exceeded a certain level, thus increasing the pulse width, and then suddenly dropped because of an electrical disturbance of the driving current of the laser diode. This behavior can be avoided with better transmitter circuit board design.

The variation of the CoG position upon increasing the bias voltage highlights whether the waveform dynamic range is dominated by a smooth rise of the power with time or by decay tails or bumps after the peak power. In particular, from Fig. 2(b), we can see that at VHV<20  V the position of the CoG is almost independent of the considered threshold, in agreement with the fact that the pulse shape is quite symmetric in the first two-orders of magnitude from the peak. On the contrary, at VHV>20  V, the CoG computed at the threshold of 0.5 is 100  ps lower than those computed at lower thresholds. This is due to the fact that the pulse shape has a steeper rise on the left and a higher tail on the right, which affects the CoG dependence with the threshold.

Upon an increase in the applied voltage, an increase of the width at level 0.5 (i.e., the FWHM) of the temporal curve can be noticed (left axis). Indeed, the FWHM is 122 ps at VHV=17  V while it reaches 253 ps at VHV=35  V. In all cases, the FWHM is lower than 260 ps. These performances are compatible with TD diffuse optics measurements, where lasers with FWHM of about 100 to 200 ps are typically used.11 Additionally, laser pulses of a few hundred picoseconds width at 0.01 of the peak value have been already demonstrated suitable for the diffuse optics application in a previous work.16 The VHV-dependent shift in the CoG of the laser shape is not a problem for diffuse optics applications, where the temporal reference is given by the acquisition of the instrument response function (see Ref. 24).

We also analyzed the emission spectrum of the laser as a function of VHV by means of a spectrometer (OceanOptics, NIRQuest series), as shown in Fig. 2(c). The mean central wavelength is about 797 nm. In Fig. 2(d), both the width of the emission spectrum at different levels from the peak wavelength and the centers of gravity as a function of VHV are shown. The spectral width increases upon higher VHV, while for the center of gravity there is no clear trend. Indeed, in the range of VHV between 17 and 20 V, the CoG moves toward shorter wavelengths for increasing VHV, while for VHV>20  V it shifts toward longer wavelengths. The spectral width (FWHM) of the laser [Fig. 2(d), left axis) is always lower than 8 nm. This value is very close to the one of the state-of-the-art TD spectroscopy setup25 that features a spectral bandwidth of 7 nm at 1350 nm. For this reason, the proposed laser can be considered suitable also for spectroscopy measurements.

Using a power meter (Thorlabs GmbH, calibration: 800 nm), we measured the average power emitted from the laser at the different VHV values. Results are shown in Fig. 2(e). The measured power has a clear linear behavior, as underlined by the linear fit (black dotted line) that approximates with high accuracy the experimental data (R2>0.999). Due to heat dissipation issues of the driver chip, it was not possible to explore the laser power behavior for VHV>34  V. Even at this low repetition rate (1 MHz), the average power is suitable for diffuse optical measurements, being of the same order of magnitude of the output power of state-of-the-art fiber-based TD clinical systems.26 The measured power consumption at 1 MHz and at VHV=18  V is 50 mW (160 mW @ VHV=34  V), thus making this source also compatible with battery operation.

Finally, we tested the output power stability of the laser system following the standardized procedures defined by the basic instrumental performances protocol.24 The working point was set at VHV=18  V, where the laser pulse features an FWHM of 100 ps and an average output power of 62  μW. Even if the emitted power is only 12% of the maximum value obtained at VHV=34  V, the narrow pulse (99 ps FWHM) should ensure better performance in recovering optical properties of highly scattering media. In this case, we decided to use an SiPM module as a detector (1  mm2 active area, see Ref. 27 for details), which is insensitive to thermal perturbations, thus ensuring the possibility to properly monitor possible fluctuations of the laser optical power. Hence, photons exiting from the laser were coupled using a collimator into a silica fiber (core diameter: 1 mm); a VOA hosted in a U-bracket was used to properly attenuate the power to fit the single-photon statistics;3 attenuated light was collected by another 1-mm-core silica fiber (whose propagation modes were completely filled using a thin layer of Teflon) coupled to the SiPM module due to a doublets of lens performing 1:1 imaging (f=30  mm). We recorded the DTOF for about 100 min with an integration period of 10 s.

For each DTOF, the background noise was subtracted and then both the overall number of counts (i.e., area of the curve) and the moments (at 0.5, 0.1, and 0.01 levels) were computed. Results are shown in Fig. 2(f). After an initial warm-up of 20 min, the number of counts and the center of gravity are almost stable (within ±1% and ±10  ps, respectively), apart from fast changes possibly related to instabilities that could be improved by employing proper heat dissipation strategies.


Validation on Phantom

With the aim of using a technology that already allows the fabrication of wearable devices, for the phantom measurements we decided to use the probe-hosted SiPM detector reported in Ref. 17, thus also maximizing the collection of diffused light from the phantom due to the large numerical aperture of the SiPM and also avoiding the use of both optical fibers or lenses. While the detector was placed directly in contact with the phantom, the laser was inserted into a metallic case (with a hole drilled for laser output) to avoid electromagnetic interferences between the pulsed laser current and the electrical output of the detector, which can introduce distortions in the DTOF. A VOA was inserted between the laser and the sample surface to adjust the injected power in order not to exceed the maximum count rate allowed by the TCSPC technique to avoid pile-up distortions (i.e., a count rate of a few percent of the laser pulse rate of 1 MHz). VHV was set to 18 V.


MEDPHOT Protocol

To assess the suitability of the proposed compact setup in retrieving the optical properties—absorption (μa) and reduced scattering coefficient (μs)—of homogeneous phantoms, we adopted the linearity test defined by the MEDPHOT protocol.28 The linearity test consists of the measurement of four series of solid phantoms featuring different reduced scattering coefficients (labeled A, B, C, and D with μs=2.8, 5.6, 8.5, and 11.27  cm1 at 800 nm, respectively). Each series is composed of eight phantoms with constant μs and covering a wide range of absorption coefficients (labeled from 1 to 8 and spanning μa from 0 to 0.32  cm1 at step of 0.05  cm1 at 800 nm).28 Despite the low repetition rate, the high laser power permitted measuring all phantoms in transmittance geometry (thickness of the phantom around 45 mm) and we collected the DTOF with 60 repetitions of 1 s each.

To obtain the values of μa and μs, the measured DTOFs were fitted using an analytical solution of the radiative transport equation under the diffusion approximation,29 convoluted with the instrument response function (previously measured, data not shown).

Results are shown in Fig. 3 where measured μa and μs are shown against their conventionally true values for the kit of phantoms employed.28 The first graph [Fig. 3(a)] shows the linearity in recovering the absorption coefficient of the phantoms against their conventionally true values. The behavior is linear for all phantoms series (minimum R2=0.9977 for phantoms with μs=2.8  cm1). For all of the series, the slope of the fit is only a few percent higher than expected (from worst case: C-series with a recovered slope of 1.16 instead of 1, thus leading to an error of 16%).

Fig. 3

Results from MEDPHOT protocol. (a) Displays absorption linearity, (b) the coupling of the absorption coefficient to the scattering one, (c) the scattering linearity, (d) the coupling of the scattering coefficient to the absorption one. For all measurements the laser driver was operating at VHV of 18 V.


Figure 3(b) shows the dependence of the recovered μa on the conventionally true μs, thus highlighting a possible scattering-to-absorption coupling, while Fig. 3(c) enlightens the absorption-to-scattering coupling. For both cases, ideally, a horizontal line should be displayed, meaning that there is no cross talk between μa and μs. As expected, the coupling is higher for higher absorption due to deviation from the diffusion approximation.

The last graph in Fig. 3(d) reports the linearity for the recovered reduced scattering coefficient displaying overall a good linearity (with R2>0.9895). Apart from the “null” absorption case with a slope of 1.32 instead of 1 due to some criticalities in dealing with very low absorption, for the other seven series the discrepancy in slope is always <15%.


nEUROPt Protocol

With the same setup described in Sec. 3.1, we tested the capability to detect a perturbation buried inside a homogeneous background. For this purpose, we used the “contrast” (C) and “contrast-to-noise ratio” (CNR) tests defined in the “nEUROPt protocol.”30 For its implementation, we used a switchable solid phantom (described in Ref. 31) that allows one to move an absorbing perturbation along one direction. The optical properties of the background phantom were μa=0.1  cm1 and μs=8  cm1 while the perturbation features a Δμa of 0.17  cm1 in a volume of 1  cm3.

Measurements were done in reflectance geometry (source–detector distance: 3 cm) with the inclusion in the center between the laser and the detector. The perturbation was moved in depth, from 5 to 35 mm (steps of 5 mm) with respect to the surface of the phantom. We computed contrast and CNR in time windows (i.e., subdividing the whole DTOF curve into portions) as defined in the nEUROPt protocol.30 In particular, the contrast was calculated summing up all the counts acquired within selected time windows (tw) of the DTOFs for each position of the inclusion (z) using


in which N0 is the number of photons collected in the unperturbed case (inclusion far—i.e., deeper than 7 cm—from source and detector, thus its effect is completely negligible), while N is the number of photons collected with the perturbation (i.e., inhomogeneous medium). Each time window was 500 ps wide and their temporal positions (i.e., opening of the time window) were referred to as the peak of the IRF. This calculation was performed for all depths z of the inclusion.

For CNR, we made use of the definition given in Ref. 30


CNR=  ΔN  (tw)σ[N0(tw)]  =N0(tw)  N(tw)σ[N0(tw)],
where ΔN was taken over the 60 repetitions for the time window under consideration. The denominator σ[N0(tw)] was computed as the standard deviation (over the 60 repetitions) of the counts summed in the homogeneous acquisitions. Figures 4(a) and 4(b) show the graphs of contrast and CNR as a function of depth, respectively, computed at different delays of the time window.

Fig. 4

(a) Contrast and (b) CNR computed for different temporal windows with a perturbation Δμa of 0.17  cm1 in a volume of 1  cm3.


As expected: (i) the contrast decreases upon increasing depth of the inclusion due to the broadening of the sensitivity profile and (ii) deeper inclusions (z>5  mm) are more visible at longer delays of the time window (>500  ps) due to the longer path-length of late-arrival photons.

Considering a CNR=1 as a threshold for detectability of a perturbation, the perturbation can be detected at depths up to 25 mm (with time window starting at 2000 ps: C=0.03 and CNR=1.77 in Fig. 3). Notwithstanding the compactness of both the laser system and the detector, these values are comparable to other bulkier state-of-the-art systems.13


Validation In Vivo

As a last proof of concept of the proposed laser device, we used the same setup described for the protocols on phantom to perform in vivo measurements of venous occlusions on healthy volunteers. The protocol was approved by the ethical committee of Politecnico di Milano and was conducted in compliance with the declaration of Helsinki. All volunteers gave their written informed consent to participate.

We varied the separation between source and detector to have a level of signal on the detector down to the single-photon level. In this way, the count rate was about 120, 35, and 150 kcps at a source–detector distance of 5.5, 4.5, and 5 cm, respectively, for subjects 1, 2, and 3. The minimum separation between injection and detection points was 4.5 cm because the dimension of the cases hosting the source and detector, which were mounted on two different printed circuit boards.

The occlusion was induced by a sphygmomanometer cuff placed on the subject’s forearm. The measurement protocol was made of three subsequent occlusions, more precisely we had: 120 s of baseline (i.e., no obstruction of blood flow in the arm) followed by a block repeated three times consisting of 60 s of occlusion followed by 60 s of recovery, and then a final recovery phase of 120 s. To induce a venous occlusion (i.e., complete block of venous blood circulation), the pressure applied to the cuff had a value between the own diastolic and systolic ones of each subject measured before starting the experiment.

We made use of the definition given in Eq. (1) to compute contrast, where N0 is the mean number of counts acquired during the first 120 s of rest state. The time window used for analysis is 500 ps wide and started 2.6 ns after the peak of the IRF. For this experiment, we decided to adopt the time-gated analysis (see Ref. 32) to increase the signal-to-noise ratio with respect to the fitting procedure, which would have instead lead to the recovery of absorption and reduced scattering coefficients.

Figure 5 shows the contrasts obtained for the three subjects. In almost all cases, there is a clear task-related trend, with an increase of contrast (i.e., a decrease of counts) during the phases in which the occlusion was performed (gray rectangles in the figure) meaning an increase in absorption. Then, when the occlusion was removed, the contrast went back close to baseline values, as expected.

Fig. 5

Results for three different subjects for the venous occlusion measurements.


Considering the central wavelength of the laser of 796 nm, which is close to the isosbestic point of oxy- and deoxyhemoglobin (located at 797 nm), the system is sensible to the overall concentration of hemoglobin, thus the increase in contrast is in agreement with the expected increase in the blood volume due to the venous occlusion. Indeed, during the occlusion, the deoxyhemoglobin cannot be washed out, while the arterial blood flow is not altered. For subject 1, results are less clear, possibly due to a suboptimally working occlusion.



In this paper, we evaluated the performance of a compact laser device based on the combination of a CMOS ASIC driver coupled to a heterostructure laser diode operated in enhanced gain-switching mode in the TD diffuse optics application. Originally, this laser transmitter was developed for time-of-flight range finding applications. Due to its miniaturized size, this technology is potentially able to allow the fabrication of pulsed sources with the same compactness of continuous-wave sources employed in currently commercially available diffuse optics systems, but allowing much better performances due to inherent advantages of measurements performed in the TD, such as (i) possibility to disentangle absorption from scattering when measuring optical properties of highly scattering media, (ii) reduced sensitivity to motion artifacts, (iii) higher penetration depth in the tissue, and (iv) possibility to disentangle optical properties of different layers of the measured tissue due to the unique feature of the TD to encode the information of different depths to different arrival delays of photons with respect to the injection time. Additionally, the low power consumption of 50 to 160 mW (depending on the operative conditions) makes this component compatible with battery operation.

To offer a wide overview of the performances, we characterized the laser system at different levels, starting from its basic performances (optical pulse shape, emission spectrum, and average optical power). Additionally, we combined the source with one of the most recent developments of a compact time-resolved detector based on SiPM technology, which can be easily hosted directly in the probe in contact with the tissue under investigation due its small footprint. The resulting TD diffuse optics system was characterized on phantoms following rigorous protocols for performance assessment of this kind of a system both in recovering the optical properties of homogeneous diffusive media and in detecting localized absorption perturbations. The obtained performances are comparable with state-of-the-art systems that make use of extremely more bulky and expensive lasers and detectors. As a final validation, we tested the system in vivo during a venous cuff occlusion applied on the forearm of three healthy volunteers. Again, the system performance in detecting signal changes due to variations of the total blood concentration during the occlusion is in line with state-of-the-art diffuse optics systems based on previous generation not-scalable technologies.26

The results reported here could represent the overcoming of the last bottleneck preventing the design of compact, wearable, and inexpensive TD diffuse optical devices. Together with other advances in photonics components, such as fully integrated TCSPC detection chains,18 and probe-hosted solid-state detectors, it is expected within the next few years the development of miniaturized fully operative high-performance TD systems, opening new market scenarios for medical devices but also for consumer appliances.


No conflicts of interest, financial or otherwise, are declared by the authors.


This project has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement no. 654148 Laserlab-Europe, the Academy of Finland (Centre of Excellence in Laser Scanning Research, contract nos. 272196, 255359, 283075, and 251571) and the Finnish Funding Agency for Innovation (TEKES).


1. D. T. Delpy et al., “Estimation of optical pathlength through tissue from direct time of flight measurement,” Phys. Med. Biol. 33(12), 1433–1442 (1988).PHMBA70031-9155 http://dx.doi.org/10.1088/0031-9155/33/12/008 Google Scholar

2. B. Chance et al., “Comparison of time-resolved and -unresolved measurements of deoxyhemoglobin in brain,” Proc. Natl. Acad. Sci. U. S. A. 85(14), 4971–4975 (1988). http://dx.doi.org/10.1073/pnas.85.14.4971 Google Scholar

3. W. Becker, Advanced Time-Correlated Single-Photon Counting, Springer, Berlin (2005). Google Scholar

4. P. Taroni et al., “Noninvasive assessment of breast cancer risk using time-resolved diffuse optical spectroscopy,” J. Biomed. Opt. 15(6), 060501 (2010).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.3506043 Google Scholar

5. L. Di Sieno et al., “Functional near-infrared spectroscopy at small source-detector distance by means of high dynamic-range fast-gated SPAD acquisitions: first in-vivo measurements,” Proc SPIE 8804, 880402 (2013). http://dx.doi.org/10.1117/12.2032492 Google Scholar

6. D. T. Delpy and M. Cope, “Quantification in tissue near-infrared spectroscopy,” Philos. Trans. R. Soc. B: Biol. Sci. 352(1354), 649–659 (1997). http://dx.doi.org/10.1098/rstb.1997.0046 Google Scholar

7. T. Durduran et al., “Bulk optical properties of healthy female breast tissue,” Phys. Med. Biol. 47(16), 2847–2861 (2002).PHMBA70031-9155 http://dx.doi.org/10.1088/0031-9155/47/16/302 Google Scholar

8. T. Durduran et al., “Diffuse optics for tissue monitoring and tomography,” Rep. Prog. Phys. 73(7), 076701 (2010).RPPHAG0034-4885 http://dx.doi.org/10.1088/0034-4885/73/7/076701 Google Scholar

9. J. Steinbrink et al., “Determining changes in NIR absorption using a layered model of the human head,” Phys. Med. Biol. 46(3), 879–896 (2001).PHMBA70031-9155 http://dx.doi.org/10.1088/0031-9155/46/3/320 Google Scholar

10. F. Martelli et al., “There’s plenty of light at the bottom: statistics of photon penetration depth in random media,” Sci. Rep. 6, 27057 (2016).SRCEC32045-2322 http://dx.doi.org/10.1038/srep27057 Google Scholar

11. A. Pifferi et al., “New frontiers in time-domain diffuse optics, a review,” J. Biomed. Opt. 21(9), 091310 (2016).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.21.9.091310 Google Scholar

12. L. Di Sieno et al., “Diffuse optics using a dual window fast-gated counter,” Appl. Opt. 53(31), 7394–7401 (2014). http://dx.doi.org/10.1364/AO.53.007394 Google Scholar

13. A. Torricelli et al., “Time domain functional NIRS imaging for human brain mapping,” NeuroImage 85(1), 28–50 (2014).NEIMEF1053-8119 http://dx.doi.org/10.1016/j.neuroimage.2013.05.106 Google Scholar

14. M. Buttafava et al., “A compact two-wavelength time-domain NIRS system based on SiPM and pulsed diode lasers,” IEEE Photonics J. 9(1), 7800114 (2017). http://dx.doi.org/10.1109/JPHOT.2016.2632061 Google Scholar

15. A. Dalla Mora et al., “Towards next-generation time-domain diffuse optics for extreme depth penetration and sensitivity,” Biomed. Opt. Express 6(5), 1749–1760 (2015).BOEICL2156-7085 http://dx.doi.org/10.1364/BOE.6.001749 Google Scholar

16. D. Contini et al, “Effects of time-gated detection in diffuse optical imaging at short source-detector separation,” J. Phys. D: Appl. Phys. 48(4), 045401 (2015).JPAPBE0022-3727 http://dx.doi.org/10.1088/0022-3727/48/4/045401 Google Scholar

17. R. Re et al., “Probe-hosted silicon photomultipliers for time-domain functional near-infrared spectroscopy: phantom and in vivo tests,” Neurophotonics 3(4), 045004 (2016). http://dx.doi.org/10.1117/1.NPh.3.4.045004 Google Scholar

18. C. Veerappan et al., “A 160×128 single-photon image sensor with on-pixel 55 ps 10b time-to-digital converter,” in IEEE Int. Solid-State Circuits Conf. Digest of Technical Papers (ISSCC 2011), pp. 312–314 (2011). http://dx.doi.org/10.1109/ISSCC.2011.5746333 Google Scholar

19. J. Nissinen and J. Kostamovaara, “A high repetition rate CMOS driver for high-energy sub-NS laser pulse generation in SPAD-based time-of-flight range finding,” IEEE Sens. J. 16(6), 1628–1633 (2016).ISJEAZ1530-437X http://dx.doi.org/10.1109/JSEN.2015.2503774 Google Scholar

20. B. Ryvkin et al., “Asymmetric-waveguide laser diode for high-power optical pulse generation by gain switching,” J. Lightwave Technol. 27(12), 2125–2131 (2009).JLTEDG0733-8724 http://dx.doi.org/10.1109/JLT.2008.2009075 Google Scholar

21. B. S. Ryvkin et al., “Quantum well laser with an extremely large active layer width to optical confinement factor ratio for high-energy single picosecond pulse generation by gain switching,” Semicond. Sci. Technol. 26(4), 045010 (2011).SSTEET0268-1242 http://dx.doi.org/10.1088/0268-1242/26/4/045010 Google Scholar

22. B. S. Ryvkin et al., “Vertical cavity surface emitting lasers with the active layer position detuned from standing wave antinode for picosecond pulse generation by gain switching,” J. Appl. Phys. 110(12), 123101 (2011).JAPIAU0021-8979 http://dx.doi.org/10.1063/1.3668321 Google Scholar

23. B. Lanz et al., “Performance improvement by a saturable absorber in gain-switched asymmetric-waveguide laser diodes,” Opt. Express 21(24), 29780–29791 (2013).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.21.029780 Google Scholar

24. H. Wabnitz et al., “Performance assessment of time-domain optical brain imagers, part 1: basic instrumental performance protocol,” J. Biomed. Opt. 19(8), 086010 (2014).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.19.8.086010 Google Scholar

25. S. Konugolu Venkata Sekar et al., “Broadband (600–1350 nm) time-resolved diffuse optical spectrometer for clinical use,” IEEE J. Sel. Top. Quantum Electron. 22(3), 7100609 (2016). http://dx.doi.org/10.1109/JSTQE.2015.2506613 Google Scholar

26. D. Contini et al., “Multi-channel time-resolved system for functional near infrared spectroscopy,” Opt. Express 14(12), 5418–5432 (2006).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.14.005418 Google Scholar

27. E. Martinenghi et al., “Time-resolved single-photon detection module based on silicon photomultiplier: a novel building block for time-correlated measurement systems,” Rev. Sci. Instrum. 87(7), 073101 (2016).RSINAK0034-6748 http://dx.doi.org/10.1063/1.4954968 Google Scholar

28. A. Pifferi et al., “Performance assessment of photon migration instruments: the MEDPHOT protocol,” Appl. Opt. 44(11), 2104–2114 (2005).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.44.002104 Google Scholar

29. D. Contini et al., “Photon migration through a turbid slab described by a model based on diffusion approximation. I. Theory,” Appl. Opt. 36(19), 4587–4599 (1997).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.36.004587 Google Scholar

30. H. Wabnitz et al., “Performance assessment of time-domain optical brain imagers, part 2: nEUROPt protocol,” J. Biomed. Opt. 19(8), 086012 (2014).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.19.8.086012 Google Scholar

31. A. Pifferi et al., “Mechanically switchable solid inhomogeneous phantom for performance tests in diffuse imaging and spectroscopy,” J. Biomed. Opt. 20(12), 121304 (2015).JBOPFO1083-3668 http://dx.doi.org/10.1117/1.JBO.20.12.121304 Google Scholar

32. D. Contini et al., “Novel method for depth-resolved brain functional imaging by time-domain NIRS,” Proc. SPIE 6629, 662908 (2007).PSISDG0277-786X http://dx.doi.org/10.1117/12.728104 Google Scholar


Laura Di Sieno received her master’s degree in electronics engineering at the Politecnico di Milano, Italy, in 2011, and her PhD in physics at the same university in 2015. Since January 2015, she has been a postdoctoral fellow. Her activity is mainly focused on the study and application of instrumentation for time-resolved optical spectroscopy of highly scattering media.

Jan Nissinen received his MSc Engineering, Licentiate of Tech. and Dr. Tech. degrees in electrical engineering from the University of Oulu, Finland, in 2002, 2007, and 2011, respectively. He has been a research scientist from 2002 to 2011 and postdoctoral researcher since 2011 in Electronics Laboratory at the University of Oulu. His research interests include the design of mixed-signal integrated circuits for pulsed time-of-flight laser rangefinders and for pulsed Raman spectroscopy.

Lauri Hallman received his MSc and Dr. Tech. degrees in electrical engineering from the University of Oulu, Oulu, Finland, in 2006 and 2015, respectively. He has done research in laser transmitter development and pulsed time-of-flight range finding. Since 2016, he has been a postdoctoral researcher at the Circuits and Systems Research Unit, University of Oulu.

Edoardo Martinenghi received his master’s degree in electronics engineering from Politecnico di Milano in 2013 and his PhD in physics in 2016. Currently, he is working as a postdoc researcher. His research interests include design of photodetectors for time-domain diffuse optics and its application to biological media.

Davide Contini received his master’s degree in electronics engineering and his PhD in physics from the Politecnico di Milano, Milan, Italy, in 2003 and 2007, respectively. Since 2014, he has been an associate professor in physics at the Department of Physics, Politecnico di Milano. His research activity is mainly focused on the study of time-resolved spectroscopy of highly diffusive media for applications in biology and medicine.

Antonio Pifferi is a full professor in the Department of Physics, Politecnico di Milano, Italy. He received his MS degree in nuclear engineering from Politecnico di Milano in 1991 and his PhD in physics from Politecnico di Torino in 1995. He was a visiting scientist at Lund Institute of Technology, Sweden, from 1995 to 1996. His research focuses on the development of new laser techniques and instrumentation for diagnosis, and on the study of light propagation in diffusive media, with applications to optical biopsy, optical mammography, and functional brain imaging.

Juha Kostamovaara received his degrees of Dipl. Eng. and Dr. Tech in electrical engineering from the University of Oulu, Finland, in 1980 and 1987, respectively. In 1995, he was invited to become full professor of electronics at the University of Oulu, where he is currently also the head of the Circuits and Systems Research Unit. His main interests are in the development of high-speed electronic circuits and systems and their applications in electronic and optoelectronic measurements.

Alberto Dalla Mora is an assistant professor of physics at Politecnico di Milano, Italy, since 2011. From Politecnico di Milano, he received cum laude his master’s degree in electronics engineering and his PhD in information and communication technology in 2006 and 2010, respectively. He authored more than 40 papers in international peer-reviewed journals and more than 40 conference proceedings. His research interests include time-resolved diffuse optics techniques and instrumentation for biomedical applications.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Laura Di Sieno, Jan Nissinen, Lauri Hallman, Edoardo Martinenghi, Davide Contini, Antonio Pifferi, Juha Kostamovaara, Alberto D. Mora, "Miniaturized pulsed laser source for time-domain diffuse optics routes to wearable devices," Journal of Biomedical Optics 22(8), 085004 (19 August 2017). http://dx.doi.org/10.1117/1.JBO.22.8.085004 Submission: Received 4 April 2017; Accepted 27 July 2017
Submission: Received 4 April 2017; Accepted 27 July 2017

Diffuse optical imaging

Pulsed laser operation


Laser sources

Picosecond phenomena


Semiconductor lasers

Back to Top