ABSTRACT
SIGNIFICANCE: Diffuse correlation spectroscopy (DCS) is an optical technique that measures blood flow non-invasively and continuously. The time-domain (TD) variant of DCS, namely, TD-DCS has demonstrated a potential to improve brain depth sensitivity and to distinguish superficial from deeper blood flow by utilizing pulsed laser sources and a gating strategy to select photons with different pathlengths within the scattering tissue using a single source-detector separation. A quantitative tool to predict the performance of TD-DCS that can be compared with traditional continuous wave DCS (CW-DCS) currently does not exist but is crucial to provide guidance for the continued development and application of these DCS systems. AIMS: We aim to establish a model to simulate TD-DCS measurements from first principles, which enables analysis of the impact of measurement noise that can be utilized to quantify the performance for any particular TD-DCS system and measurement geometry. APPROACH: We have integrated the Monte Carlo simulation describing photon scattering in biological tissue with the wave model that calculates the speckle intensity fluctuations due to tissue dynamics to simulate TD-DCS measurements from first principles. RESULTS: Our model is capable of simulating photon counts received at the detector as a function of time for both CW-DCS and TD-DCS measurements. The effects of the laser coherence, instrument response function, detector gate delay, gate width, intrinsic noise arising from speckle statistics, and shot noise are incorporated in the model. We have demonstrated the ability of our model to simulate TD-DCS measurements under different conditions, and the use of our model to compare the performance of TD-DCS and CW-DCS under a few typical measurement conditions. CONCLUSION: We have established a Monte Carlo-Wave model that is capable of simulating CW-DCS and TD-DCS measurements from first principles. In our exploration of the parameter space, we could not find realistic measurement conditions under which TD-DCS outperformed CW-DCS. However, the parameter space for the optimization of the contrast to noise ratio of TD-DCS is large and complex, so our results do not imply that TD-DCS cannot indeed outperform CW-DCS under different conditions. We made our code available publicly for others in the field to find use cases favorable to TD-DCS. TD-DCS also provides a promising way to measure deep brain tissue dynamics using a short source-detector separation, which will benefit the development of technologies including high density DCS systems and image reconstruction using a limited number of source-detector pairs.
Subject(s)
Hemodynamics , Photons , Computer Simulation , Monte Carlo Method , Spectrum AnalysisABSTRACT
Significance: Time domain diffuse correlation spectroscopy (TD-DCS) can offer increased sensitivity to cerebral hemodynamics and reduced contamination from extracerebral layers by differentiating photons based on their travel time in tissue. We have developed rigorous simulation and evaluation procedures to determine the optimal time gate parameters for monitoring cerebral perfusion considering instrumentation characteristics and realistic measurement noise. Aim: We simulate TD-DCS cerebral perfusion monitoring performance for different instrument response functions (IRFs) in the presence of realistic experimental noise and evaluate metrics of sensitivity to brain blood flow, signal-to-noise ratio (SNR), and ability to reject the influence of extracerebral blood flow across a variety of time gates to determine optimal operating parameters. Approach: Light propagation was modeled on an MRI-derived human head geometry using Monte Carlo simulations for 765- and 1064-nm excitation wavelengths. We use a virtual probe with a source-detector separation of 1 cm placed in the pre-frontal region. Performance metrics described above were evaluated to determine optimal time gate(s) for different IRFs. Validation of simulation noise estimates was done with experiments conducted on an intralipid-based liquid phantom. Results: We find that TD-DCS performance strongly depends on the system IRF. Among Gaussian pulse shapes, â¼ 300 ps pulse length appears to offer the best performance, at wide gates (500 ps and larger) with start times 400 and 600 ps after the peak of the TPSF at 765 and 1064 nm, respectively, for a 1-s integration time at photon detection rates seen experimentally (600 kcps at 765 nm and 4 Mcps at 1064 nm). Conclusions: Our work shows that optimal time gates satisfy competing requirements for sufficient sensitivity and sufficient SNR. The achievable performance is further impacted by system IRF with â¼ 300 ps quasi-Gaussian pulse obtained using electro-optic laser shaping providing the best results.