News

Jul 08, 2023

Portable, high speed blood flow measurements enabled by long wavelength, interferometric diffuse correlation spectroscopy (LW

Scientific Reports volume 13, Article number: 8803 (2023) Cite this article Diffuse correlation spectroscopy (DCS) is an optical technique that can be used to characterize blood flow in tissue. The

Scientific Reports volume 13, Article number: 8803 (2023) Cite this article

Diffuse correlation spectroscopy (DCS) is an optical technique that can be used to characterize blood flow in tissue. The measurement of cerebral hemodynamics has arisen as a promising use case for DCS, though traditional implementations of DCS exhibit suboptimal signal-to-noise ratio (SNR) and cerebral sensitivity to make robust measurements of cerebral blood flow in adults. In this work, we present long wavelength, interferometric DCS (LW-iDCS), which combines the use of a longer illumination wavelength (1064 nm), multi-speckle, and interferometric detection, to improve both cerebral sensitivity and SNR. Through direct comparison with long wavelength DCS based on superconducting nanowire single photon detectors, we demonstrate an approximate 5× improvement in SNR over a single channel of LW-DCS in the measured blood flow signals in human subjects. We show equivalence of extracted blood flow between LW-DCS and LW-iDCS, and demonstrate the feasibility of LW-iDCS measured at 100 Hz at a source-detector separation of 3.5 cm. This improvement in performance has the potential to enable robust measurement of cerebral hemodynamics and unlock novel use cases for diffuse correlation spectroscopy.

Diffuse correlation spectroscopy (DCS) is an established optical technique that allows for the non-invasive measurement of tissue blood flow1. Through the measurement of diffusely backscattered light, DCS relates the temporal fluctuations of the collected signals to the motion of blood cells through the vasculature. Bedside clinical blood flow monitoring2, especially cerebral blood flow monitoring3, has exploded as a use case for DCS, with DCS having been used to estimate metrics of cerebral perfusion during surgical procedures4,5,6,7,8, cerebral autoregulation9,10, cerebrovascular reactivity11, intracranial pressure12,13,14, and critical closing pressure15,16. While a number of studies including DCS monitoring have been demonstrated in adult populations, due to limitations in cerebral sensitivity and signal-to-noise ratio17, the standard DCS technique is better suited to measure blood flow in neonates and children, where the extracerebral tissue (scalp and skull) is significantly thinner than in adults18,19. To improve the performance of DCS in adult populations, many groups have developed modifications on DCS which provide improvements to cerebral sensitivity, signal-to-noise ratio, or both. These methods include interferometric detection20,21,22,23,24,25, parallelized speckle detection26,27,28, acousto-optic modulation29,30,31, pathlength resolved methods32,33,34,35,36,37, speckle contrast methods38,39,40, and long wavelength approaches41,42. Recent work in our group has shown the utility of the use of long wavelength DCS applied at 1064 nm, though in practice for clinical measurements, currently available commercial detectors either do not have reasonable noise performance for measurements sensitive to deep flow (InGaAs/InP single-photon avalanche diodes (SPADs))43 or are too bulky to be applied clinically (superconducting nanowire single-photon detectors (SNSPD)). To address this gap in detector technology, we have developed long wavelength, interferometric DCS (LW-iDCS), which takes advantage of all the benefits of working at 1064 nm and sidesteps the negative aspects of detector technologies sensitive to light at 1064 nm using interferometric detection together with a highly parallel line scan camera sensor (inspired by the work done at shorter wavelengths by Zhou et al.21,44). In this work, we directly compare the performance of LW-DCS and LW-iDCS in a pilot human subject study to verify the equivalence of blood flow estimation by the new LW-iDCS technique and compare the quality of the measured signals.

Diffuse correlation spectroscopy estimates the flow in tissue through the analysis of the normalized intensity autocorrelation function, \({g}_{2}\left(\tau \right)\). The autocorrelation of the detected signal is connected to the dynamics of the tissue by the Siegert relation45, expressed in Eq. (1),

where \({g}_{1}\left(\tau \right)\) is the normalized electric field, temporal autocorrelation function and β is the coherence parameter46, which is related to the coherence length of the source, the geometry of the measurement, number of modes detected, and the degree of environmental light contamination. The Siegert relation connects the measured signals to the underlying fluctuations of the electric field due dynamic scattering events. The electric field autocorrelation function in the DCS measurement can be described as an integral of individual pathlength-specific correlation functions over the pathlength distribution detected. This form, given in Eq. (2)47, allows for the connection between the measured intensity autocorrelation function and the dynamics in the tissue.

where P(s) is the distribution of pathlengths, s, taken by photons in the tissue, k0 is the wavenumber of the detected light in a vacuum, n is the index of refraction of the sample, \(\langle \Delta {r}^{2}\left(\tau \right)\rangle\) is the mean squared displacement of the scattering particles, and l* is the reduced, mean free path of photons in the tissue which is described as the inverse of the tissue’s reduced scattering coefficient \(\left({l}^{*}=\frac{1}{{\mu }_{s}^{^{\prime}}}\right)\). For DCS measurements in tissue, the mean squared displacement term is assumed to reflect diffusive motion48 \(\left(\langle \Delta {r}^{2}\left(\tau \right)\rangle =6\alpha {D}_{b}\tau =6B{F}_{i}\tau \right)\), where the blood flow index (BFi) describes the effective diffusion coefficient, which reflects the true diffusion coefficient (Db) multiplied by the probability of scattering events occurring at moving scatterers (α). While this description of flow in vessels as a diffusive process is slightly puzzling, multiple theoretical and simulation studies have examined the appropriateness of the model to describe the detected signals, and found the diffusive process as a good description under standard DCS measurement conditions48,49,50,51, though some conflicting theories have been proposed52. When fitting correlation curves in this study, the model selected for \({g}_{1}\left(\tau \right)\) is that reflecting a semi-infinite sample measured in the reflectance geometry, given in Eq. (3),

where \(K\left(\tau \right)=\sqrt{3{\mu }_{a}{\mu }_{s}^{\mathrm{^{\prime}}}+6{k}_{0}^{2}{n}^{2}{\mu }_{s}^{\mathrm{^{\prime}}2}B{F}_{i}\tau }\), μa is the optical absorption coefficient, \({r}_{1}=\sqrt{{\rho }^{2}+{{l}^{*}}^{2}}\), ρ is the distance between the source and detector, \({r}_{b}=\sqrt{{\rho }^{2}+{\left({l}^{*}+2{z}_{b}\right)}^{2}}\), \({z}_{b}=\frac{2}{3{\mu }_{s}^{\mathrm{^{\prime}}}}\frac{\left(1+{R}_{eff}\right)}{\left(1-{R}_{eff}\right)}\), and \({R}_{eff}\left(n\right)=-1.440{n}^{-2}+0.71{n}^{-1}+0.668+0.0636n\).

To improve the signal-to-noise ratio of blood flow estimates made by DCS, our group and others have implemented DCS system which utilize interferometric detection20,21,22,23,24. These approaches are attractive as they allow for an intrinsic improvement to the signal-to-noise ratio of DCS derived blood flow signals53 as well as allow for the use of less sensitive, nosier detectors. In this study, we utilize a Mach–Zehnder interferometer, which combines the diffusely scattered light from the sample with a reference signal split from the laser. This results in a form of \({g}_{2}\left(\tau \right)\) which is different than the homodyne \({g}_{2}\left(\tau \right)\), and is given in Eq. (4),

where \(\langle {I}_{S}\rangle\) is the average light intensity collected from the sample, \(\langle {I}_{R}\rangle\) is the average light intensity of the reference arm of the interferometer, and \(\langle {I}_{T}\rangle\) is the average total light intensity falling on the detector.

A graphical representation of the optical instrumentation is given in Fig. 1. To make direct, collocated comparisons between homodyne LW-DCS and heterodyne LW-iDCS, a custom fiber optic probe was constructed, similar to the probe reported previously42, to deliver light from the laser and return light to the detectors. The fiber optic probe contained two adjacent source fibers (> 3.5 mm center-to-center distance), 1 single mode fiber for short-separation DCS (5 mm) and several co-localized long-separation detection fibers: 4 single mode fibers (LW-DCS), and 7 multimode detection fibers (LW-iDCS). A high coherence (lc > 10 km), fiber (MFD 6.6 µm) laser source emitting ~ 125 mW at 1064 nm (RFLM-125-0-1064, NP Photonics) was fusion spliced (S185HS Fusion Splicer, Fitel) to a 90:10, polarization maintaining fused fiber coupler (MFD 6.6 µm, PN1064R2A1, Thorlabs). The 10% arm of the coupler was used as the input for a fiber amplifier (MAKO-AMP1064, Cybel), and was connected via an FC/APC connector. The amplifier output fiber (MFD 10 µm) was fusion spliced to the input of a 50:50, 105 µm, multimode fused fiber coupler (TW1064R5A1B, Thorlabs). The two outputs of the fiber coupler were spliced to two 105 µm multimode source fibers connected to the probe. The light was amplified to allow for two MPE limited spots54 (1 W/cm2 at 1064 nm, 3.6 mm spot size diameter, 102 mW each spot) to increase the achievable signal-to-noise ratio. The 90% output arm of the polarization maintaining coupler was connected to the reference arm input of the LW-iDCS interferometer. All spliced connections were confirmed by the fusion splicer to have losses less than 0.03 dB.

Optical instrumentation used in this work. A long coherence, 1064 nm laser was coupled to a 90%/10% fused fiber coupler to split the light into a reference arm for the interferometer (90%) and the seed source for the fiber amplifier (10%). The amplified source light was split by a 50%/50% fused fiber coupler to supply two MPE limited sources. Single mode fibers were placed at 5 mm (1) and at 35 mm (4) to bring the light to the SPAD detector and the SNSPDs, respectively. The single photon detection events were time tagged at 6.67 ns resolution and transferred to the computer via USB 3.0. Multimode fibers were also placed at 35 mm (7), which brought light to the sample arm of the interferometer. Light from both the reference and sample arms was shaped to match the size of the linescan camera array (12.5 µm × 25.6 mm), and the intensity signals from the camera were digitized at 300 kHz.

For the homodyne LW-DCS measurements, a custom LW-DCS detection system based on single-photon detectors [one silicon SPAD (SPCM-NIR-14, Excelitas) and four SNSPDs tuned to 1064 nm (Opus One, Quantum Opus)] and a custom FPGA-based software correlator board with a timing resolution of 6.67 ns was used55. Single mode detector fibers were used to bring the backscattered light from the optical probe to the detectors. The silicon SPAD was used to detect light collected at a source-detector separation of 5 mm (i.e. short separation), and the four SNSPDs were used to detect light collected at a source-detector separation of 35 mm (i.e. long separation). For the short separation signal, autocorrelation functions were calculated at a rate of 10 Hz due to the limited photon counts available. For long separation signals, autocorrelation functions for each channel were calculated at a rate of 100 Hz and subsequently averaged.

For the heterodyne LW-iDCS measurements, a free space Mach Zehnder interferometer was constructed to project the light from both the sample and reference optical fibers onto a fast InGaAs linescan camera (Manx 2048 SQ CXP 390, Xenics). The 90% arm of the polarization maintaining coupler was connected to a variable optical attenuator (VOA1064-APC, Thorlabs), which was connected to a fiber collimator (F220APC-1064, Thorlabs) placed in a kinematic mount providing XY translation (CXY1A, Thorlabs) as well as tip and tilt (KAD11F, Thorlabs). A Powell lens (LGL130, Thorlabs) was used to homogenize the intensity of the reference arm along the dimension of pixel array of the linescan camera. The diverging output of the Powell lens was collimated by a cylindrical lens (LJ1765L1-C, Thorlabs). For the sample arm, the seven 200 µm, multi-mode fibers in the optical probe placed at a source-detector separation of 35 mm were aligned linearly and bundled into an SMA connector (FG200LEA, BFL200LS02, Thorlabs) and collimated using a fiber collimator (F220SMA-1064, Thorlabs). To correct the remaining angular divergence in the sample arm, a spherical lens (LA1131-B-ML, Thorlabs) was used. Both the sample and reference beams were sent to a 50:50 non-polarizing beam splitter (CCM1-BS014, Thorlabs) and the combined output was focused onto the linear array of the camera using a cylindrical lens (LJ1328L2-B, Thorlabs). The ratio of sample arm intensity to reference arm intensity during the measurements was estimated to be 1 to 5 × 107. This was estimated by scaling the estimated average power per fiber mode at 3.5 cm at 1064 nm (~ 1.2 × 10−12 mW per source fiber (~ 6 kcps)42) by the number of modes contained in the 7 detector fibers (~ 59,000) and dividing by the total reference arm power used (~ 6 mW). The interferometer was built on an optical breadboard (MBH1224, Thorlabs) with vibration isolating feet (AV4, Thorlabs) and placed on a rolling cart (61 × 46 × 122 cm3) for mobility. The collected signals were digitized at a line rate of 300 kHz. Raw data were captured by a frame grabber (Coaxlink Quad G3, Euresys) and saved directly to disk at a rate of 1.2 GB/s and post processed to estimate BFi. Several signal processing steps were completed to maximize measurement SNR and remove camera induced distortions prior to converting the recorded pixel intensity data into autocorrelation functions. These steps included processing to address hardware induced distortions caused by the integrate-while-read mode, quadratically detrending the signals over the analysis interval, averaging adjacent pixel signals, and removal of common noise signals across the camera. These steps are detailed both graphically and in text in the supplement. Correlation functions from the LW-iDCS instrument were computed at a rate of 100 Hz. For multilayer tissue, in the process of assessing the simulated performance of the different DCS implementations, it was discovered that if the semi-infinite model is used to fit the data and if the same percentage of the decay is used for fitting (i.e. correlation function decays to 5% of the value of the plateau), the sensitivity of the measurement to the cerebral signal differs between DCS, based on fitting a function that is proportional to \({\left|{g}_{1}\left(\tau \right)\right|}^{2}\), and iDCS, based on fitting a function that is proportional to \(\left|{g}_{1}\left(\tau \right)\right|\). To address this, LW-iDCS correlation functions, proportional to \({g}_{1}\left(\tau \right)\), are fit with a weighted fitting approach to match the cerebral sensitivity of the LW-DCS measurement, which reflect a correlation function based on \({\left|{g}_{1}\left(\tau \right)\right|}^{2}\). The applied weighting approach was optimized using a Monte Carlo simulation of light transport and momentum transfer in a multi-layer, slab geometry, detailed in Table S1, that is meant to represent the typical tissue geometry of a measurement made over the forehead.

For this study we enrolled five healthy subjects (3 female, 2 male, age 38 ± 19, included participants of Middle Eastern (1), European (2), Southeast Asian (1), and East Asian (1) descent, 2 with dark skin pigmentation) to compare and validate the performance of the LW-iDCS system against the standard 1064 nm DCS system across several physiological manipulations. This study was reviewed and approved by the Mass General Brigham Institutional Review Board (#2019P003074). All participants gave written informed consent prior to the measurements. All methods were performed in accordance with the relevant guidelines and regulations. The measurement protocol included three tasks: breath holding, hyperventilation, and tourniquet pressure modulation.

Due to the lack of head coverage in this study, systemic physiological manipulations were selected to induce large, repeatable changes in the measured blood flow. Breath holding and hyperventilation were performed to perturb heart rate56, blood pressure57,58, and vasoactive state59,60 to cause changes in both scalp and brain blood flow. The tourniquet pressure modulation technique is used to selectively reduce the blood flow in the scalp to increase the specificity of the blood flow signal to cerebral blood flow61. By comparing the reduction in blood flow at the short separation channel and the long separation channel, this method also allows the for the assessment of the long channel’s sensitivity to the cerebral signal. A list of the timing of the activity and recovery intervals can be seen below in Table 1. In addition to the LW-DCS optical instrumentation, systemic physiological monitoring was performed, and included electrocardiography (ECG), pulse oximetry (SpO2), continuous, non-invasive blood pressure monitoring (Nova, Finapres), and respiratory monitoring. These signals were digitized at 1 kHz by a Powerlab ADC (ADInstruments). Data collection was synchronized between the measurement devices using an external trigger box which was connected to an auxiliary analog input channel of the LW-DCS FPGA correlator and a channel on the Powerlab ADC. Rough synchronization for the LW-iDCS instrument was accomplished by providing a trigger to the other instruments when the LW-iDCS acquisition was started. Fine tuning of the synchronization was accomplished by determining the time offset between the LW-DCS and LW-iDCS around the trigger which maximized the cross-correlation between the BFi signals.

To compare the blood flow changes during the physiological manipulations, we first removed the cardiac pulsation from the blood flow signals. To remove the influence of the cardiac pulsatility, using the previously identified RR intervals from the ECG, averaged beat-to-beat blood flow index values were calculated62. The data were then resampled back to the original sampling rate (i.e. 100 Hz for long separation signals, 10 Hz for short separation signals), and the periods of physiological manipulation were identified and separated. To compare between trials and between subjects, individual trial BFi values were scaled to a relative blood flow index (rBFi) by dividing by the average BFi in the 20 s preceding the activity. The individual trials for each experimental paradigm were then averaged across subjects.

To investigate the difference in measured BFi between correlation functions proportional to \(\left|{g}_{1}\left(\tau \right)\right|\) and \({\left|{g}_{1}\left(\tau \right)\right|}^{2}\), a baseline and activated condition (+50% brain BFi) are simulated63. The observed difference in the baseline BFi and sensitivity to cerebral blood flow changes between iDCS and DCS simulations of the three-layer geometry can be seen in Fig. 2A and B. One considered method to make the two measurements equivalent was to square the collected iDCS correlation function to have a signal proportional to \({\left|{g}_{1}\left(\tau \right)\right|}^{2}\). While on noise-less correlation functions, this would be the preferred method, with noisy correlation functions, taking the square of the curves could result in distortions to the noise properties of the curves, which could further affect the fitting. The alternative approach that was taken was to increase the weighting given to the earlier part of the correlation function using a data driven weighting scheme. The optimized objective function is given in Eq. (5), where the weighting coefficients are taken from the overall average correlation function from the entire measurement.

where the \({\langle \rangle }_{T}\) is the average over the measurement interval, and x is a factor that was optimized to reach equivalent BFi from the simulations. A range of factors, x, were investigated using the simulated data from the comparison of DCS and iDCS performance. Based on the tissue geometry simulated, the optimal value was found to be 2.5, and this value was fixed and used to fit the data in the manuscript. The corrected baseline blood flow and resolved changes can be seen in Fig. 2A and B.

Comparison of the BFi fit from simulated multilayer DCS measurements. (A) The BFi fit from the baseline simulations with different fitting functions can be seen. The discrepancy between the fits of \({g}_{1}\left(\tau \right)\) and \({g}_{1}{\left(\tau \right)}^{2}\) would be observed as discrepancies between the fits of iDCS and DCS, respectively. (B) The change in BFi measured in response to a 50% increase in the brain BFi is shown. Without the weighted fitting, the iDCS measurement resolves ~ 50% of the changes that DCS does, reducing the sensitivity to the cerebral signal. With weighted fitting, the fit based on \({g}_{1}\left(\tau \right)\) is seen to be equivalent to the fit based on \({g}_{1}{\left(\tau \right)}^{2}\).

To assess the equivalency of the measurements taken by the LW-iDCS instrument, blood flow traces are compared between the LW-DCS and LW-iDCS instruments. To assess both noise performance as well as equivalence of the BFi measured for the pulsatile signal, cardiac gated averaging was performed using data from the baseline periods where the variability in the shape of the cardiac pulsation and modification in heart rate can be assumed to be minimal. RR intervals were identified in the time aligned ECG traces, and LW-DCS and LW-iDCS BFi values were averaged in sync with the cardiac cycle, as was done previously15. In Fig. 3A, the time aligned blood flow trace measured by each instrument is shown for three cardiac cycles (n = 20 averaged three-cycle traces) in an example subject with the standard deviation of the estimates of blood flow shown as well. The equivalence of measured blood flow index is also confirmed across all experimental conditions and subjects. Average BFi values were calculated for each cardiac pulse identified in all subject BFi time traces across all tasks, and the BFi values from each instrument are plotted against each other in Fig. 3C, showing excellent agreement over the range of measured BFi. In Fig. 3D, a Bland–Altman plot comparing the agreement of the measured BFi between the two instruments shows a bias (4.27 × 10−10 cm2/s) for the LW-iDCS instrument to measure a faster blood flow with a standard deviation of the difference equal to 8.39 × 10−10 cm2/s. The bias and spread of the difference are relatively small compared to the range of values typically measured with DCS (Fig. 3C), and this result indicates good agreement between the two instruments. To assess the noise performance across subjects, the coefficient of variation of the measured pulsatile BFi at each point in the cardiac cycle is calculated for each instrument. These values are compared as violin plots in Fig. 3B and show an average reduction in the coefficient of variation given by LW-iDCS of ~ 2.25×. This matches quite well with the estimated improvement in contrast-to-noise ratio (CNR) estimated from Monte Carlo simulation, shown in the supplement.

Comparison of the characteristics of BFi time traces measured at 100 Hz from the LW-iDCS and LW-DCS instruments. (A) An example of a single subject, pulsatile cardiac signal is shown for both instruments, demonstrating equivalency in the measured blood flow index as well as the reduced noise of the blood flow trace measured by the LW-iDCS instrument. (B) The coefficient of variation \(\left({\sigma }_{B{F}_{i}}/{\mu }_{B{F}_{i}}\right)\) was computed for each point in the cardiac cycle, and the results for each subject for each measurement modality are shown in violin plots. On average, the reduction in coefficient of variation provided by the LW-iDCS instrument is ~ 2.25× when compared to the 4 channel LW-DCS instrument. Equivalency of the measured BFi values beyond the pulsatile signals between the two instruments is also shown across subjects and tasks using the cardiac filtered, BFi signals. (C) The measured BFi values are plotted against each other, and cluster nicely around the line of unity. (D) The Bland–Altman plot shows a narrow distribution of the differences in the measured BFi, characterized by a mean difference of 4.27 × 10−10 cm2/s and a standard deviation of 8.39 × 10−10 cm2/s, demonstrating good agreement between the two blood flow measurements.

In response to an end expiratory breath hold, typical physiological responses include an increase in blood pressure57,58 and a hypercapnic state64,65. In Fig. 4, we show the subject averaged response to a 30 s breath holding trial. The mean arterial pressure (MAP) is seen to increase 25% ± 9% by the end of the breath hold, while the heart rate remains relatively constant. The BFi increases measured at the long separation are matched well between LW-DCS and LW-iDCS (32% ± 17% at end breath hold) and are seen to be distinct from the increase measured at the short separation (51% ± 17% at end breath hold). The observed increases in blood flow are consistent with the expectation given an increase in blood pressure and the mild hypercapnic state.

Subject averaged responses to breath holding. (A) Comparison of the measured blood flow responses to the 30 s breath hold. The relative change in flow in the long separation measurements is seen to be slightly lower than the change observed at the short separation, which has been previously observed42. (B) Comparison of relative changes in blood pressure and heart rate, respectively in response to the 30 s breath hold. A progressive increase in blood pressure was observed throughout the breath holding period, while the heart rate remains relatively constant.

For the hyperventilation task, subjects performed one minute of paced breathing at seventy breaths per minute. The expected hypocapnia due to the over breathing is expected to cause vasoconstriction and a decrease in blood flow. Because cerebral and scalp tissue metabolism is maintained during this period, the reduced blood flow results in a decrease in tissue hemoglobin saturation (SO2), which induces a vasoactive reaction to return the blood flow to the baseline level59. This biphasic response was observed in both the long and short channel measurements, as seen in Fig. 5A. The observed systemic physiologic response to the bout of hyperventilation was a significant increase in the heart rate (38% ± 15%) and a significant decrease in the mean arterial pressure (−18% ± 10%), as seen in Fig. 5B. The latency of the return to baseline of these physiologic parameters was longer as compared to the latency of the return to baseline of the measured blood flow. The long separation blood flow signal responses to hyperventilation were seen to be consistent, and additionally match well to the responses observed from the long separation blood flow measurements previously reported in our group42.

Subject averaged response to hyperventilation maneuver. (A) Measured hemodynamic response to 60 s of hyperventilation. As in the breath holding task, the short separation measurement shows a more exaggerated response to the physiologic manipulation, exhibiting a decrease of 30% in BFi after the onset of hyperventilation. The matching long separation responses show a lesser degree of BFi reduction, and all blood flow can be seen to return to the baseline before the end of the hyperventilation trial. (B) For this maneuver, the heart rate increased significantly following the start of the trial, while the blood pressure was seen to reduce.

The tourniquet tightening causes selective reduction of blood flow in the scalp, leaving the brain blood flow unmodified61. By comparing the differential response between a short separation channel, sensitive almost exclusively to the blood flow in the scalp, and a long separation channel, exhibiting sensitivity to both the scalp and cerebral hemodynamics, estimates of the long separation channel’s sensitivity to each compartment can be estimated. The group averaged blood flow traces from the pressure modulation trials can be seen in Fig. 6A. In this group of subjects, the average reduction in BFi measured was 85.3% and 39.2% for the short and long channels, respectively. This pair of measurements corresponds to a long separation channel with a 46% sensitivity to the superficial blood flow signal, which implies a cerebral sensitivity of over 50%, based on a recent simulation study66. As was expected for this maneuver, the systemic physiology did not respond to the change in torniquet pressure, seen in Fig. 6B.

Subject averaged response to the pressure modulation maneuver. (A) Measured hemodynamic response to 30 s of tourniquet tightening. Using the ratio of the relative decrease between the long channel (39.2%) and the short channel (85.3%), the sensitivity of the long channel to the scalp blood flow can be estimated to be 46%. For long separation DCS measurements, brain sensitivity has been shown to be inversely proportional to scalp sensitivity66, and we can estimate that the 35 mm separation measurement has a brain sensitivity > 50%. (B) For this maneuver, as expected, the systemic physiology was not significantly affected by the tightening of the tourniquet on the forehead.

In this work we have demonstrated the development of long wavelength, interferometric diffuse correlation spectroscopy. Using a fiber optic probe with collocated detection fibers, we were able to directly compare measured blood flow and signal-to-noise ratio between the LW-DCS and LW-iDCS instruments. For both single photon detection and interferometric detection, the availability of cost-effective, high-power fiber optic amplifiers, an engineering benefit of making measurements at 1064 nm, allows for an increased SNR beyond the intrinsic benefits of using light at 1064 nm. By sacrificing spatial resolution, multiple sources spaced > 3.5 mm apart could be used, allowing for an even higher SNR for high quality pulsatile blood flow measurements. The SNR of the LW-iDCS measurement seen in the high-speed pulsatile measurements was 4.5× the SNR of the SNSPD LW-DCS measurement when making single channel comparisons, representing an enabling improvement to the quality of blood flow measured. In the context of the DCS systems currently used for translational research, this improvement is especially significant considering that even the single illumination SNSPD LW-DCS has an SNR gain of 16× over conventional DCS42, and that measurements at 3.5 cm are not feasible with conventional NIR DCS. The use of a camera which is sensitive to light at 1064 nm takes advantage of both the higher number of photons per mode as compared to traditional NIR wavelengths as well as the slower decay of the autocorrelation function. For cerebral blood flow measurements made at long source-detector separations, the autocorrelation decay for traditional NIR DCS can happen in 1–10 s of microseconds, and a significant portion of the decay could be missed if not sampled quickly enough. The use of both heterodyne detection, measuring the slower decaying \({g}_{1}\left(\tau \right)\) as opposed to \({g}_{2}\left(\tau \right)\), and 1064 nm relaxes the sampling rate needed to effectively sample the correlation function. The longer source-detector separation achievable with these advanced DCS systems enables measurements with reduced sensitivity to the upper tissue layers relative to the sensitivity of currently applied DCS systems in the traditional NIR wavelength range (explored in the supplement). The decreased sensitivity to extracerebral signals is greatly beneficial to DCS measurements, especially in clinical applications where systemic physiological fluctuations are more likely to occur and the timing of relevant cerebral hemodynamic changes is not as well defined. We also see good agreement with the estimated noise performance given by Monte Carlo simulation (Figure S3). Additionally, the cost of the system is greatly reduced compared to LW-DCS based on SNSPDs. For this implementation of the LW-iDCS system, the detector used is ~ 7× less expensive (~ $25 k total, camera + frame grabber: ~ $20 k, assorted lenses, opto-mechanics, and fibers: ~ $5 k) as compared to the SNSPDs (~ $180 k total, cryostat: ~ $100 k, individual nanowire detectors: ~ $20 k each). The LW-iDCS cart-based system is also more mobile than the SNSPD based LW-DCS system. These improvements in cost, SNR, and mobility are promising for the clinical usability of LW-iDCS measurements of CBF in adults. The signal processing approach used to extract the correlation function from the raw data stream points to potential pitfalls in the development of iDCS instruments using multimode fiber and free space interferometers though. The motion of fibers and vibrations in the environment have the potential to corrupt the iDCS signals, however, these challenges are manageable, and the use of the custom data analysis pipeline, described in supplementary information, was successful in removing artifacts from the data. The use of a weighted fitting approach allowed for equivalent blood flow indices to be fit from both the LW-DCS and LW-iDCS correlation functions, evidenced by the results shown in Fig. 3C and D. While the results presented matched well, investigation of the generalizability of the weighting factor selected in this study is warranted given the influence that tissue layer thicknesses, optical properties, and ratios of scalp and brain blood flow are known to have on fitting autocorrelation functions67,68. Another challenge posed by the implementation of massively parallel multi-speckle detection is the raw data rate of the instruments. Recent publications on massively parallelized detection have quoted raw data rates between 0.24 GB/s (0.864 TB/hr) and 9.0 GB/s (32.4 TB/hr)22,25,26,27,28,44,69. For clinical blood flow measurements, these data rates could result in untenably large data files, though real time processing utilizing GPUs or FPGAs have been explored as a solution to address this challenge28,69. The increased SNR provided by the LW-iDCS instrument presented here enabled high sensitivity to the cerebral blood flow signal as well as a high rate of BFi calculation. These factors will be highly enabling for the clinical translation of DCS as a noninvasive cerebral blood flow monitor.

The data that support the findings of this study are available from the corresponding authors upon reasonable request.

Boas, D. A. & Yodh, A. G. Spatially varying dynamical properties of turbid media probed with diffusing temporal light correlation. J. Opt. Soc. Am. A 14, 192 (1997).

Article ADS Google Scholar

Shang, Y., Li, T. & Yu, G. Clinical applications of near-infrared diffuse correlation spectroscopy and tomography for tissue blood flow monitoring and imaging. Physiol. Meas. 38, R1–R1 (2017).

Article ADS PubMed PubMed Central Google Scholar

Buckley, E. M., Parthasarathy, A. B., Grant, P. E., Yodh, A. G. & Franceschini, M. A. Diffuse correlation spectroscopy for measurement of cerebral blood flow: Future prospects. Neurophotonics 1, 11009 (2014).

Article Google Scholar

Kaya, K. et al. Intraoperative cerebral hemodynamic monitoring during carotid endarterectomy via diffuse correlation spectroscopy and near-infrared spectroscopy. Brain Sci. 12, 1025 (2022).

Article PubMed PubMed Central Google Scholar

Zavriyev, A. I. et al. The role of diffuse correlation spectroscopy and frequency-domain near-infrared spectroscopy in monitoring cerebral hemodynamics during hypothermic circulatory arrests. JTCVS Tech. 7, 161–177 (2021).

Article PubMed PubMed Central Google Scholar

Shang, Y. et al. Cerebral monitoring during carotid endarterectomy using near-infrared diffuse optical spectroscopies and electroencephalogram. Phys. Med. Biol 56, 3015–3032 (2011).

Article PubMed Google Scholar

Rajaram, A. et al. Cerebral perfusion and metabolic neuromonitoring during cardiopulmonary bypass. In Optical Tomography and Spectroscopy of Tissue (Vol. 11639, pp. 39). https://doi.org/10.1117/12.2578986 (2021).

Busch, D. R. et al. Continuous cerebral hemodynamic measurement during deep hypothermic circulatory arrest. Biomed. Opt. Express 7, 3461 (2016).

Article PubMed PubMed Central Google Scholar

Selb, J. et al. Prolonged monitoring of cerebral blood flow and autoregulation with diffuse correlation spectroscopy in neurocritical care patients. Neurophotonics 5, 1 (2018).

Article Google Scholar

Busch, D. R. et al. Detection of brain hypoxia based on noninvasive optical monitoring of cerebral blood flow with diffuse correlation spectroscopy. Neurocrit. Care 30, 72–80 (2019).

Article CAS PubMed PubMed Central Google Scholar

Milej, D. et al. Characterizing dynamic cerebral vascular reactivity using a hybrid system combining time-resolved near-infrared and diffuse correlation spectroscopy. Biomed. Opt. Express 11, 4571–4585 (2020).

Article CAS PubMed PubMed Central Google Scholar

Ruesch, A. et al. Estimating intracranial pressure using pulsatile cerebral blood flow measured with diffuse correlation spectroscopy. Biomed. Opt. Express 11, 1462 (2020).

Article PubMed PubMed Central Google Scholar

Flanders, T. M. et al. Optical detection of intracranial pressure and perfusion changes in neonates with hydrocephalus. J. Pediatr. 236, 54-61.e1 (2021).

Article PubMed PubMed Central Google Scholar

Tabassum, S. M. et al. Clinical translation of intracranial pressure sensing with diffuse correlation spectroscopy. J. Neurosurg. 1, 20 (2022).

Google Scholar

Wu, K.-C. et al. Validation of diffuse correlation spectroscopy measures of critical closing pressure against transcranial Doppler ultrasound in stroke patients. J. Biomed. Opt. 26, 36008–36009 (2021).

Article Google Scholar

Baker, W. B. et al. Noninvasive optical monitoring of critical closing pressure and arteriole compliance in human subjects. J. Cereb. Blood Flow Metab. 37, 2691–2705 (2017).

Article PubMed PubMed Central Google Scholar

Selb, J. et al. Sensitivity of near-infrared spectroscopy and diffuse correlation spectroscopy to brain hemodynamics: Simulations and experimental findings during hypercapnia. Neurophotonics 1, 15005 (2014).

Article Google Scholar

Beauchamp, M. S. et al. The developmental trajectory of brain-scalp distance from birth through childhood: Implications for functional neuroimaging. PLoS ONE 6, e24981–e24981 (2011).

Article ADS CAS PubMed PubMed Central Google Scholar

Davis, N. J. Variance in cortical depth across the brain surface: Implications for transcranial stimulation of the brain. Eur. J. Neurosci. https://doi.org/10.1111/ejn.14957 (2020).

Article PubMed Google Scholar

Xu, J., Jahromi, A. K., Brake, J., Robinson, J. E. & Yang, C. Interferometric speckle visibility spectroscopy (ISVS) for human cerebral blood flow monitoring. APL Photonics 5, 126102 (2020).

Article ADS CAS Google Scholar

Zhou, W. et al. Functional interferometric diffusing wave spectroscopy of the human brain. Sci. Adv. 7, eabe0150–eabe0150 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Zhou, W. et al. Multi-exposure interferometric diffusing wave spectroscopy. Opt. Lett. 46, 4498–4501 (2021).

Article ADS PubMed PubMed Central Google Scholar

Robinson, M. B., Boas, D. A., Sakadzic, S., Franceschini, M. A. & Carp, S. A. Interferometric diffuse correlation spectroscopy improves measurements at long source–detector separation and low photon count rate. J. Biomed. Opt. 25, 97004 (2020).

Article CAS Google Scholar

James, E., Powell, S. & Munro, P. Performance optimisation of a holographic Fourier domain diffuse correlation spectroscopy instrument. Biomed. Opt. Express 13, 3836 (2022).

Article PubMed PubMed Central Google Scholar

Samaei, S., Nowacka, K., Gerega, A., Pastuszak, Ż & Borycki, D. Continuous-wave parallel interferometric near-infrared spectroscopy (CW πNIRS) with a fast two-dimensional camera. Biomed. Opt. Express 13, 5753 (2022).

Article PubMed PubMed Central Google Scholar

Liu, W. et al. Fast and sensitive diffuse correlation spectroscopy with highly parallelized single photon detection. APL Photonics 6, 26106 (2021).

Article Google Scholar

Sie, E. J. et al. High-sensitivity multispeckle diffuse correlation spectroscopy. Neurophotonics 7, 35010 (2020).

Article CAS Google Scholar

Wayne, M. A. et al. Massively parallel, real-time multispeckle diffuse correlation spectroscopy using a 500 × 500 SPAD camera. Biomed. Opt. Express 14, 703 (2023).

Article PubMed PubMed Central Google Scholar

Ling, H., Gui, Z., Hao, H. & Shang, Y. Enhancement of diffuse correlation spectroscopy tissue blood flow measurement by acoustic radiation force. Biomed. Opt. Express 11, 301 (2020).

Article CAS PubMed Google Scholar

Robinson, M. B. et al. Characterization of continuous wave ultrasound for acousto-optic modulated diffuse correlation spectroscopy (AOM-DCS). Biomed. Opt. Express 11, 3071 (2020).

Article PubMed PubMed Central Google Scholar

Tsalach, A. et al. Depth selective acousto-optic flow measurement. Biomed. Opt. Express 6, 4871–4886 (2015).

Article PubMed PubMed Central Google Scholar

Sutin, J. et al. Time-domain diffuse correlation spectroscopy. Optica 3, 1006 (2016).

Article ADS CAS PubMed PubMed Central Google Scholar

Ozana, N. et al. Functional time domain diffuse correlation spectroscopy. Front. Neurosci. 16, 1123 (2022).

Article Google Scholar

Zhao, M., Zhou, W., Aparanji, S., Mazumder, D. & Srinivasan, V. Interferometric diffusing wave spectroscopy imaging with an electronically variable time-of-flight filter. Optica 10, 42–52 (2022).

Article ADS CAS Google Scholar

Borycki, D., Kholiqov, O. & Srinivasan, V. J. Interferometric near-infrared spectroscopy directly quantifies optical field dynamics in turbid media. Optica 3, 1471 (2016).

Article ADS PubMed PubMed Central Google Scholar

Pagliazzi, M. et al. Time resolved speckle contrast optical spectroscopy at quasi-null source-detector separation for non-invasive measurement of microvascular blood flow. Biomed. Opt. Express 12, 1499 (2021).

Article PubMed PubMed Central Google Scholar

Poon, C.-S. et al. First-in-clinical application of a time-gated diffuse correlation spectroscopy system at 1064 nm using superconducting nanowire single photon detectors in a neuro intensive care unit. Biomed. Opt. Express 13, 1344 (2022).

Article CAS PubMed PubMed Central Google Scholar

Zilpelwar, S. et al. A model of dynamic speckle evolution for evaluating laser speckle contrast measurements of tissue dynamics. Biomed. Opt. Express 13, 6533–6549 (2022).

Article PubMed PubMed Central Google Scholar

Dragojević, T. et al. Compact, multi-exposure speckle contrast optical spectroscopy (SCOS) device for measuring deep tissue blood flow. Biomed. Opt. Express 9, 322 (2018).

Article PubMed Google Scholar

Valdes, C. P. et al. Speckle contrast optical spectroscopy, a non-invasive, diffuse optical method for measuring microvascular blood flow in tissue. Biomed. Opt. Express 5, 2769 (2014).

Article PubMed PubMed Central Google Scholar

Carp, S. A. et al. Diffuse correlation spectroscopy measurements of blood flow using 1064 nm light. J. Biomed. Opt. 25, 97003–97004 (2020).

Article CAS Google Scholar

Ozana, N. et al. Superconducting nanowire single-photon sensing of cerebral blood flow. Neurophotonics 8, 35006 (2021).

Article CAS Google Scholar

Robinson, M. B. et al. Diffuse correlation spectroscopy beyond the water peak enabled by cross-correlation of the signals from InGaAs/InP single photon detectors. IEEE Trans. Biomed. Eng. 69, 1943–1953 (2022).

Article PubMed PubMed Central Google Scholar

Zhou, W., Kholiqov, O., Chong, S. P. & Srinivasan, V. J. Highly parallel, interferometric diffusing wave spectroscopy for monitoring cerebral blood flow dynamics. Optica 5, 518 (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Siegert, A. J. F. On the fluctuations in signals returned by many independently moving scatterers. (Radiation Laboratory, Massachusetts Institute of Technology, 1943).

Bellini, T., Glaser, M. A., Clark, N. A. & Degiorgio, V. Effects of finite laser coherence in quasielastic multiple scattering. Phys. Rev. A (Coll Park) 44, 5215 (1991).

Article ADS CAS Google Scholar

Pine, D. J., Weitz, D. A., Chaikin, P. M. & Herbolzheimer, E. Diffusing wave spectroscopy. Phys. Rev. Lett. 60, 1134–1137 (1988).

Article ADS CAS PubMed Google Scholar

Boas, D. A. et al. Establishing the diffuse correlation spectroscopy signal relationship with blood flow. Neurophotonics 3, 31412 (2016).

Article Google Scholar

Verdecchia, K., Diop, M., Morrison, L. B., Lee, T.-Y. & Lawrence, K. S. Assessment of the best flow model to characterize diffuse correlation spectroscopy data acquired directly on the brain. Biomed. Opt. Express 6, 4288 (2015).

Article PubMed PubMed Central Google Scholar

Sakadžic, S. et al. Theoretical model of blood flow measurement by diffuse correlation spectroscopy. J. Biomed. Opt. 22, 27006 (2017).

Article PubMed Google Scholar

Carp, S. A. et al. Due to intravascular multiple sequential scattering, diffuse correlation spectroscopy of tissue primarily measures relative red blood cell motion within vessels. Biomed. Opt. Express 2, 2047 (2011).

Article PubMed PubMed Central Google Scholar

Du Le, V. N. & Srinivasan, V. J. Beyond diffuse correlations: Deciphering random flow in time-of-flight resolved light dynamics. Opt. Express 28, 11191 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Koppel, D. Statistical accuracy in FCS. Phys. Rev. A (Coll Park) 10, 1938–1945 (1974).

Article ADS Google Scholar

American National Standard for Safe Use of Lasers. ANSI Z136.1-2007. (2007).

Farzam, P. et al. Fast diffuse correlation spectroscopy (DCS) for non-invasive measurement of intracranial pressure (ICP)(Conference Presentation). In Clinical and Translational Neurophotonics (eds. Madsen, S. J. & Yang, V. X. D.) vol. 10050 100500U-100500U (International Society for Optics and Photonics, 2017).

Fischer, K., Guensch, D. P. & Friedrich, M. G. Response of myocardial oxygenation to breathing manoeuvres and adenosine infusion. Eur. Heart J. Cardiovasc. Imaging 16, 395–401 (2015).

Article PubMed Google Scholar

Parkes, M. J., Green, S., Stevens, A. M. & Clutton-Brock, T. H. Assessing and ensuring patient safety during breath-holding for radiotherapy. Br. J. Radiol. 87, 20140454 (2014).

Article CAS PubMed PubMed Central Google Scholar

Perini, R. et al. Heart rate and blood pressure time courses during prolonged dry apnoea in breath-hold divers. Eur. J. Appl. Physiol. 104, 1–7 (2008).

Article ADS PubMed Google Scholar

Wilson, D. F. et al. Effect of hyperventilation on oxygenation of the brain cortex of newborn piglets. J. Appl. Physiol. 70, 2691–2696 (1991).

Article CAS PubMed Google Scholar

Meyer, J. S., Gotoh, F., Takagi, Y. & Kakimi, R. Cerebral hemodynamics, blood gases, and electrolytes during breath-holding and the Valsalva maneuver. Circulation 33, II–35 (1966).

Article CAS Google Scholar

Baker, W. B. et al. Pressure modulation algorithm to separate cerebral hemodynamic signals from extracerebral artifacts. Neurophotonics 2, 35004 (2015).

Article Google Scholar

Skytioti, M., Søvik, S. & Elstad, M. Respiration-related cerebral blood flow variability increases during control-mode non-invasive ventilation in normovolemia and hypovolemia. Eur. J. Appl. Physiol. 117, 2237–2249 (2017).

Article PubMed Google Scholar

Fang, Q. & Boas, D. A. Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units. Opt. Express 17, 20178–20190 (2009).

Article ADS CAS PubMed Google Scholar

Larsen, J. et al. Breath holding for 20 s following extended expiration is a practical, effective and robust standard when measuring cerebrovascular reactivity in healthy adults using BOLD fMRI at 3 T. Neuroimage Rep. 1, 100021 (2021).

Article Google Scholar

Zerweck, L., Hauser, T.-K., Roder, C. & Klose, U. Investigation of the BOLD-based MRI signal time course during short breath-hold periods for estimation of the cerebrovascular reactivity. SN Compr. Clin. Med. 2, 1551–1562 (2020).

Article CAS Google Scholar

Wu, M. M. et al. Complete head cerebral sensitivity mapping for diffuse correlation spectroscopy using subject-specific magnetic resonance imaging models. Biomed. Opt. Express 13, 1131 (2022).

Article PubMed PubMed Central Google Scholar

Wu, M. M. et al. Improved accuracy of cerebral blood flow quantification in the presence of systemic physiology cross-talk using multi-layer Monte Carlo modeling. Neurophotonics 8, 15001 (2021).

Article Google Scholar

Zhao, H., Sathialingam, E. & Buckley, E. M. Accuracy of diffuse correlation spectroscopy measurements of cerebral blood flow when using a three-layer analytical model. Biomed. Opt. Express 12, 7149–7161 (2021).

Article PubMed PubMed Central Google Scholar

James, E. & Powell, S. Fourier domain diffuse correlation spectroscopy with heterodyne holographic detection. Biomed. Opt. Express 11, 6755 (2020).

Article PubMed PubMed Central Google Scholar

Download references

Funding was provided by National Institute of Biomedical Imaging and Bioengineering (Grant Nos. T32EB001680, R01EB033202, U01EB028660, R21EB028626 (NIBIB), F31NS118753 (NINDS)).

Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, USA

Mitchell B. Robinson, Marco Renna, Nisan Ozana, Alyssa N. Martin, Nikola Otic, Stefan A. Carp & Maria Angela Franceschini

Bar-Ilan University, Tel Aviv District, Ramat Gan, Israel

Nisan Ozana

Department of Biomedical Engineering, Boston University, Boston, MA, USA

Nikola Otic

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

M.B.R, S.A.C., and M.A.F conceptualized the study, M.B.R, M.R., N.N.O., and M.A.F planned the experiments, M.B.R., M.R., A.N.M., and N.O. carried out the experiments, M.B.R analyzed the data, all authors interpreted and discussed the data, M.B.R., M.R., S.A.C. and M.A.F prepared the manuscript, and all authors reviewed the manuscript.

Correspondence to Mitchell B. Robinson.

M.A.F had a financial interest in 149 Medical, Inc., a company developing DCS technology for assessing and monitoring cerebral blood flow in newborn infants. M.A.F.’s interests are managed by Mass General Hospital, and Mass General Brigham in accordance with their conflict-of-interest policies. M.B.R, M.R., N.N.O., A.N.M, N.O., and S.A.C. have nothing to report.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Robinson, M.B., Renna, M., Ozana, N. et al. Portable, high speed blood flow measurements enabled by long wavelength, interferometric diffuse correlation spectroscopy (LW-iDCS). Sci Rep 13, 8803 (2023). https://doi.org/10.1038/s41598-023-36074-8

Download citation

Received: 26 February 2023

Accepted: 29 May 2023

Published: 31 May 2023

DOI: https://doi.org/10.1038/s41598-023-36074-8

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.