Non-invasive measurement of hemodynamic change during 8 MHz transcranial focused ultrasound stimulation using near-infrared spectroscopy
BMC Neuroscience volume 20, Article number: 12 (2019)
Transcranial focused ultrasound (tFUS) attracts wide attention in neuroscience as an effective noninvasive approach to modulate brain circuits. In spite of this, the effects of tFUS on the brain is still unclear, and further investigation is needed. The present study proposes to use near-infrared spectroscopy (NIRS) to observe cerebral hemodynamic change caused by tFUS in a noninvasive manner.
The results show a transient increase of oxyhemoglobin and decrease of deoxyhemoglobin concentration in the mouse model induced by ultrasound stimulation of the somatosensory cortex with a frequency of 8 MHz but not in sham. In addition, the amplitude of hemodynamics change can be related to the peak intensity of the acoustic wave.
High frequency 8 MHz ultrasound was shown to induce hemodynamic changes measured using NIRS through the intact mouse head. The implementation of NIRS offers the possibility of investigating brain response noninvasively for different tFUS parameters through cerebral hemodynamic change.
Ultrasound is a widely used imaging tool for diagnosis and monitoring, which is gradually finding its way to therapeutic use. In particular, transcranial focused ultrasound (tFUS) has been demonstrated to be capable of modulating brain structures with considerable advantages over other neuromodulation techniques due to its combined noninvasive use, ability to penetrate deep brain structures, and precise targeting resolution [1, 2].
Ultrasound has already been implemented on animal models for neuroprotection after stroke , and in humans for neurological disorders , but the application of tFUS is still limited because of the incomplete understanding of its working mechanism and physiological effects on the brain . Moreover, the method of selection of stimulation parameters that would specify the ability to either excite or suppress neural activity is also not determined, and thus further investigation is in demand .
Simulations modelling the human head have shown that the optimal frequencies for US transmission through the skull is below 0.7 MHz . Higher frequencies may offer better spatial resolution but at the risk of increased attenuation and heating of the skull. For small animals however, the thinner skull reduces the effects of high frequency ultrasound absorption in the bone, potentially providing precise spatial stimulation specificity critical to animal studies. What remains is to determine whether the transmitted energy will be enough to induce brain activity in mice under safe sonication conditions.
Several studies measure physiological responses to assess the effects and to influence parameter selection for tFUS. The supporting techniques used in these studies have their strengths and weaknesses. The observation of a limb twitch along with electromyography is a common procedure to evaluate the effects of acoustic stimulation on brain circuits . Even though the approach is straightforward, it fails to measure the immediate effects of tFUS in the brain region that do not have a quantifiable behavior change. On the other hand, more sophisticated techniques like functional magnetic resonance imaging (fMRI)  and positron emission tomography (PET)  are already well-established as valuable tools to observe indirect brain activity induced by tFUS, although they can be technically demanding, expensive, and with low temporal resolution; in addition to exposure to ionizing radiation in PET. In this work, we propose to use near-infrared spectroscopy (NIRS) as an alternative noninvasive approach that addresses the weakness of previous techniques to investigate the cerebral hemodynamic changes induced by tFUS.
In a manner similar to fMRI, NIRS can be applied to record changes in cerebral blood oxygenation related to brain activity during the execution of a specific paradigm. Although NIRS may not provide deep brain information, it is an attractive and accessible approach to investigate brain activity changes with its high temporal resolution, robustness to electrical artifacts, and low cost . NIRS measures changes in intensity at different wavelengths, permitting estimation of cerebral hemoglobin concentration changes. The measure of change in oxygenated (HbO) and deoxygenated (RHb) hemoglobin could be used as an indirect way of recording brain activity due to the neurovascular coupling mechanism. NIRS has been used in various neuroscience studies including post-stroke rehabilitation , the anesthetic effect on the brain , functional brain connectivity , and investigation of other stimulation techniques . We believe NIRS can be a suitable noninvasive tool to observe cerebral hemodynamic changes induced by tFUS as well as to identify the optimal stimulation parameters for a desired application. In this study, we aim to present the feasibility of NIRS to detect the brain hemodynamic change induced by different acoustic intensities with a central frequency of 8 MHz.
The acoustic beam profile measured in water is shown on Fig. 1. The focal zone extends 1.76 mm in the lateral and 0.29 mm in the axial direction. Placing an excised mouse skull in between the transducer and hydrophone leads to an attenuation of 9 dB, maintaining 12% of the original acoustic power.
Cerebral hemodynamic changes were calculated from the NIRS signals obtained from 10 mice using the modified Beer–Lambert Law . The 5 s ultrasound stimulation with repetition rate of 2 Hz (Fig. 1c) induced a change in cerebral hemodynamics as shown from the averaged results from all animals in Fig. 2a–c. The dramatic increase of HbO and decrease of RHb was consistently observed from all the mice during the stimulation but not in sham. Two simulation conditions with different acoustic intensities (I-SPTA = 468 mW/cm2, 1077 mW/cm2) were delivered to the brain with the aim to observe the existence of any relationship between stimulation power and cerebral hemodynamic change. Increasing the acoustic intensity of stimulation also increases the observed hemodynamic changes in the mice. One-way repeated measures ANOVA show a statistically significant difference in maximum amplitude change of HbO (f(2, 18) = 19.9, p < 0.001) and RHb (f(2,18) = 15, p < 0.001), but not total hemoglobin THb (f(2,18) = 2.55, p = 0.1) between all three stimulation conditions. Tukey–Kramer post hoc analysis indicates that all three stimulation conditions produced hemodynamic profiles that are statistically distinct from each other (Fig. 2d). To assess the relationship between the hemodynamic profiles of the two stimulation conditions, a Pearson product-moment correlation was computed. Strong similarity for HbO and RHb was found with r2 > 0.9 for both.
The focus of this preliminary study was to determine if a reliable signal could be obtained with the simultaneous NIRS and tFUS. Although the fiber-based NIRS is difficult to implement with a larger number of channels due to the limited scalp space on a small animal model, it is still interesting due to the possibility of easy translation to human neuroimaging applications providing adequate spatial and temporal resolution .
A consistent pattern of cerebral hemodynamic change was observed across all animals tested during ultrasound stimulation. The results were in accordance with our previous report showing hemodynamic response initiated by tFUS using a minimally invasive neuroimaging technique . In the previous study, wide-field cerebral hemodynamics was monitored using optical intrinsic signal imaging through an intact skull cranial window. Although the parameters and ultrasound frequency used are different from the current study, the stimulation evoked a similar pattern of hemodynamics as observed non-invasively through the scalp and skull using NIRS. Both acoustic stimulation protocols induced a pronounced increase and decrease of HbO and RHb respectively, a pattern that could be referred to as the general hemodynamic response to neural activity. The results also exhibited a post-stimulus undershoot similarly observed in the common BOLD response. The post-stimulus undershoot is presumed to be related to a prolonged consumption of oxygen  and known to be affected by stimulation conditions .
Although typically, the sub-megahertz frequency is recommended to use for tFUS due to the acoustic transparency of the human skull in this range, higher frequencies up to 43 MHz have been demonstrated to modulate neural excitation . Higher frequencies offer better stimulation precision at the expense of attenuation through the skull. Therefore, finding a suitable frequency with sufficient transmission above 1 MHz may be more important for small animal studies. In a previous study, ultrasound with a frequency of 1.9 MHz has been shown to stimulate the motor cortex to produce repeatable contralateral limb movements in mice . In addition, another research team demonstrated that the anatomical specificity of ultrasound stimulation to induce a tail response was greatly improved by using 5 MHz compared to 1 MHz . However, since the motor response is limited to brain regions associated with inducing motion and can be affected by various parameters including the number of stimulated neurons, the absence of a motor response may not necessarily be related to the stimulation being insufficient to induce brain activity changes. In this work, high frequency 8 MHz tFUS applied to the somatosensory area with different acoustic intensities has been shown to be able to induce cerebral hemodynamic changes as probed by NIRS. Even with the observed changes in cerebral hemodynamics, no visible motor response was observed. Nonetheless, the mice still consistently showed motor response to fingers snapping due to the auditory startle reflex, in contrast with a recent study that elicited a startle response for both audible sound and low frequency ultrasound of 500 kHz . The high acoustic frequency also provides an improvement in stimulation targeting in terms of focal size and minimizing the interference of reflected waves inside the skull [25, 26]. Focused sonication based on 8 MHz is highly absorbed by the mouse skull but, with an f–number of 1, can provide a focal diameter of approximately 200 um (1,540,000 m*s−1/8,000,000 Hz), that is about the same stimulation area for optogenetics based on a 200 um fiber . Moreover, it affords the ability to obtain an adequate spatial resolution for imaging including functional ultrasound imaging, providing the possibility to implement a single ultrasound system for both neuromodulation and imaging.
In addition to central frequency, there are several important tFUS parameters that might affect the output of stimulation, among which is the peak acoustic intensity. The results using NIRS has shown the influence of acoustic intensity to the amplitude but not to the temporal profile of hemoglobin concentration change. The higher acoustic intensity induced a greater increase in cerebral blood oxygenation, which could be associated to an increase in neural activity and could play an important role in neuroprotection to minimize the effect of ischemic stroke .
Even though the obtained hemodynamic response agrees with the typical pattern induced by neural activity , the question of whether the obtained hemodynamic response originates from direct ultrasound neural stimulation or is only due to a cerebral vascular response is still unresolved. One possible solution would be the direct recording of neural activity through electrophysiology  or metabolic measurements using PET . NIRS has also been shown to measure changes in the oxidation state of cytochrome-c-oxidase (oxCCO) which would indicate cellular energy metabolism . However, the small concentration change of oxCCO compared to hemoglobin chromophores make the extraction of oxCCO a challenging task requiring additional computational complexity in the NIRS system .
The mechanical index of the applied ultrasound in this study (MI = 0.2) was much lower compared to the safety limit set for clinical ultrasound imaging applications (MI = 1.9), minimizing the likelihood of any cavitation-related brain tissue damage . After passing through the skull, the acoustic intensities of the stimulation (ISPTA = 129 and 56 mW/cm2) were also lower than the guidelines for ultrasound imaging (ISPTA = 720 mW/cm2). Furthermore, previous reports with higher acoustic energy (ISPTA = 6.4 W/cm2 in rabbit) did not show tissue damage . The maximum temperature change of brain tissue was estimated to be less than 0.01 °C using previously described equations valid for short exposure times . Obvious abnormalities in mouse behavior were not observed after stimulation.
The work shows conceptually the ability of NIRS to noninvasively measure brain changes related to tFUS. The method is sensitive enough to show that the high frequency (8 MHz) causes the hemodynamic change, where the amplitude of cerebral hemodynamics are positively correlated to acoustic peak power. The authors believe that NIRS can be an important supplementary tool to investigate the effect of tFUS on cerebral hemodynamic change, as well as to find optimal acoustic parameters for a desired application.
A total of ten female BALB/c mice weighing 17–20 g (9–10 weeks old) (Damul Science, Korea) were used for the experiment. The number of animals was determined using power analysis (power = 0.95, effect size = 1.3, η2 = 0.62 based on HbO data) . The mice were housed in cages with 3–4 occupants at a 12-h light/dark cycle, having access to feed and water ad libitum. One day before the measurement, the mice were anesthetized using an intraperitoneal injection of ketamine-xylazine cocktail (80:10 mg/kg, respectively) to depilate the scalp. The animals were allowed to recover at least 24 h before the tFUS experiment.
Ultrasound stimulation was achieved using a focused single element ultrasound transducer with a focal distance of 19 mm, and 8 MHz fundamental frequency. The input sinusoidal pulses were generated from a function generator (Agilent 33220A, Keysight, USA) amplified at 50 dB (E&I 240L, USA). A 3D printed acoustic guide was attached to the end of the ultrasound transducer to provide convenient access to the mouse head. The ultrasound intensity and beam profile were measured from the tip of the waveguide in water and through the mouse skull using an acoustic intensity measurement system (AIMS III, ONDA, USA). The beam profile measured in water showed that the focus extends 1.76 mm in the axial and 0.29 mm in the lateral directions (Fig. 1a, b). Passing through the skull caused a 9 dB attenuation of the acoustic wave, keeping only 12% of the original power.
The continuous wave near-infrared spectroscopy setup consisted of two 400 µm optical fibers, 4 mm apart, positioned over the mouse scalp. One fiber was connected to a halogen broadband light source (HL-2000, Ocean Optics) and the other to a spectrometer (USB 4000, Ocean Optics), providing a single channel system. The 3D printed fiber holder provided additional restraint to the mouse head and optical fiber tips located at 0 mm anteroposterior (AP), and + 1 mm mediolateral (ML) for the light source fiber; and − 4 mm AP, + 1 mm ML for the detector fiber. The approximate position of the bregma was found by naked eye under sufficient illumination. Light intensity from 700 to 900 nm in 1 nm increments was converted to hemodynamic concentration change (HbO and RHb) using the modified Beer–Lambert law , Eq. 1, where SD is the separation distance between the source and detector, ε is the absorption coefficient, and ΔA is the difference in absorption between a time point and at a designated baseline (time t = 0). A mean optical pathlength at 720 nm was estimated by fitting the second derivative of the water feature at 720 nm and HbO feature at 760 nm  assuming a mean water content of 80% in the rodent brain . The differential pathlength factor (DPF) at each wavelength was then derived by applying wavelength-dependent correction factors obtained from literature . The THb concentration change is the sum of the obtained HbO and RHb. The processed data were low pass filtered at 0.1 Hz to reduce physiological high frequency respiratory and cardiac noises.
Each animal was initially anesthetized with 3% isoflurane to restrain the animal in a stereotaxic frame. Once fixed, the animal was given at least 30 min to recover from anesthesia before the start of data acquisition. All animals received three different stimulation conditions, including sham, within a single acquisition experiment. The stimulation paradigms were selected in a way to keep the same pulse repetition frequency (PRF) of 1500 Hz and duty cycle of 9% but with different acoustic intensities (as measured in water, ISPTA) of 1077 and 468 mW/cm2. The acoustic parameters were chosen based on previous studies inducing a motor response for lower frequency ultrasound targeted to the motor cortex, and also considering sufficient acoustic power even after attenuation through the skull. The transmitted acoustic intensity (12%) is in a similar range with previously reported values for ultrasound brain stimulation in mice . A diagram illustrating the pulse modulation for one pulse train is shown in Fig. 1c.
The experiment consisted of 75 randomly shuffled trials, with 25 trials for each acoustic condition. Each acquisition trial lasts for a total of 17 s including 5 s of stimulation containing ten ultrasound pulse trains with a repetition rate of 2 Hz. The interval between the start of each trial was fixed to 1 min. The ultrasound probe coupled with the acoustic guide was filled with acoustic gel and placed on the top of the mouse head pointing towards the somatosensory cortex. The ultrasound wave was focused along with the optical fibers for NIRS in such a way that the stimulated and monitored areas overlap. Sham trials were applied using the same protocol and setup as stimulation trials but without any input to the transducer. After the experiment, the mice were routinely observed for 24 h but no apparent behavioral change were detected. The mice were then euthanized using CO2 inhalation followed by cervical dislocation.
differential pathlength factor
functional magnetic resonance imaging
oxidation state of cytochrome-c-oxidase
positron emission tomography
pulse repetition frequency
transcranial focused ultrasound
Tufail Y, Yoshihiro A, Pati S, Li MM, Tyler WJ. Ultrasonic neuromodulation by brain stimulation with transcranial ultrasound. Nat Protoc. 2011;6(9):1453.
Legon W, Ai L, Bansal P, Mueller JK. Neuromodulation with single-element transcranial focused ultrasound in human thalamus. Hum Brain Mapp. 2018;39(5):1995–2006.
Guo T, Li H, Lv Y, Lu H, Niu J, Sun J, et al. Pulsed transcranial ultrasound stimulation immediately after the ischemic brain injury is neuroprotective. IEEE Trans Bio-med Eng. 2015;62(10):2352–7.
Monti MM, Schnakers C, Korb AS, Bystritsky A, Vespa PM. Non-invasive ultrasonic thalamic stimulation in disorders of consciousness after severe brain injury: a first-in-man report. Brain Stimul. 2016;9(6):940–1.
Fini M, Tyler WJ. Transcranial focused ultrasound: a new tool for non-invasive neuromodulation. Int Rev Psychiatry. 2017;29(2):168–77.
Kubanek J. Neuromodulation with transcranial focused ultrasound. Neurosurg Focus. 2018;44(2):E14.
Tyler WJ. Noninvasive neuromodulation with ultrasound? A continuum mechanics hypothesis. Neuroscientist. 2011;17(1):25–36.
Ye PP, Brown JR, Pauly KB. Frequency dependence of ultrasound neurostimulation in the mouse brain. Ultrasound Med Biol. 2016;42(7):1512–30.
Lee W, Kim H-C, Jung Y, Chung YA, Song I-U, Lee J-H, et al. Transcranial focused ultrasound stimulation of human primary visual cortex. Sci Rep. 2016;6:34026.
Yoo S-S, Bystritsky A, Lee J-H, Zhang Y, Fischer K, Min B-K, et al. Focused ultrasound modulates region-specific brain activity. Neuroimage. 2011;56(3):1267–75.
Scholkmann F, Kleiser S, Metz AJ, Zimmermann R, Pavia JM, Wolf U, et al. A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology. Neuroimage. 2014;85:6–27.
Mihara M, Miyai I. Review of functional near-infrared spectroscopy in neurorehabilitation. Neurophotonics. 2016;3(3):031414.
Choi D-H, Shin TJ, Kim S, Bae J, Cho D, Ham J, et al. Monitoring of cerebral oxygenation and local field potential with a variation of isoflurane concentration in a rat model. Biomed Opt Express. 2016;7(10):4114–24.
Zhang Y-J, Lu C-M, Biswal BB, Zang Y-F, Peng D, Zhu C-Z. Detecting resting-state functional connectivity in the language system using functional near-infrared spectroscopy. J Biomed Opt. 2010;15(4):047003.
Yan J, Wei Y, Wang Y, Xu G, Li Z, Li X. Use of functional near-infrared spectroscopy to evaluate the effects of anodal transcranial direct current stimulation on brain connectivity in motor-related cortex. J Biomed Opt. 2015;20(4):046007.
Boas DA, Gaudette T, Strangman G, Cheng X, Marota JJ, Mandeville JB. The accuracy of near infrared spectroscopy and imaging during focal changes in cerebral hemodynamics. Neuroimage. 2001;13(1):76–90.
White BR, Culver JP. Quantitative evaluation of high-density diffuse optical tomography: in vivo resolution and mapping performance. J Biomed Opt. 2010;15(2):026006.
Kim E, Anguluan E, Kim JG. Monitoring cerebral hemodynamic change during transcranial ultrasound stimulation using optical intrinsic signal imaging. Sci Rep. 2017;7(1):13148.
Schroeter ML, Kupka T, Mildner T, Uludağ K, von Cramon DY. Investigating the post-stimulus undershoot of the BOLD signal—a simultaneous fMRI and fNIRS study. Neuroimage. 2006;30(2):349–58.
van Zijl PC, Hua J, Lu H. The BOLD post-stimulus undershoot, one of the most debated issues in fMRI. Neuroimage. 2012;62(2):1092–102.
Menz MD, Oralkan Ö, Khuri-Yakub PT, Baccus SA. Precise neural stimulation in the retina using focused ultrasound. J Neurosci. 2013;33(10):4550–60.
Kamimura HA, Wang S, Chen H, Wang Q, Aurup C, Acosta C, et al. Focused ultrasound neuromodulation of cortical and subcortical brain structures using 1.9 MHz. Med Phys. 2016;43(10):5730–5.
Li G-F, Zhao H-X, Zhou H, Yan F, Wang J-Y, Xu C-X, et al. Improved anatomical specificity of non-invasive neuro-stimulation by high frequency (5 MHz) ultrasound. Sci Rep. 2016;6:24738.
Sato T, Shapiro MG, Tsao DY. Ultrasonic neuromodulation causes widespread cortical activation via an indirect auditory mechanism. Neuron. 2018;98:1031–41.
Younan Y, Deffieux T, Larrat B, Fink M, Tanter M, Aubry JF. Influence of the pressure field distribution in transcranial ultrasonic neurostimulation. Med Phys. 2013;40(8):082902.
Constans C, Deffieux T, Pouget P, Tanter M, Aubry J-F. A 200–1380-kHz quadrifrequency focused ultrasound transducer for neurostimulation in rodents and primates: transcranial in vitro calibration and numerical study of the influence of skull cavity. IEEE Trans Ultrason Ferroelectr Freq Control. 2017;64(4):717–24.
Yizhar O, Fenno LE, Davidson TJ, Mogri M, Deisseroth K. Optogenetics in neural systems. Neuron. 2011;71(1):9–34.
Murry CE, Jennings RB, Reimer KA. Preconditioning with ischemia: a delay of lethal cell injury in ischemic myocardium. Circulation. 1986;74(5):1124–36.
Villringer A, Planck J, Hock C, Schleinkofer L, Dirnagl U. Near infrared spectroscopy (NIRS): a new tool to study hemodynamic changes during activation of brain function in human adults. Neurosci Lett. 1993;154(1–2):101–4.
Roche-Labarbe N, Wallois F, Ponchel E, Kongolo G, Grebe R. Coupled oxygenation oscillation measured by NIRS and intermittent cerebral activation on EEG in premature infants. Neuroimage. 2007;36(3):718–27.
Dierks T, Jelic V, Pascual-Marqui RD, Wahlund L-O, Julin P, Linden DE, et al. Spatial pattern of cerebral glucose metabolism (PET) correlates with localization of intracerebral EEG-generators in Alzheimer’s disease. Clin Neurophysiol. 2000;111(10):1817–24.
Bale G, Elwell CE, Tachtsidis I. From Jöbsis to the present day: a review of clinical near-infrared spectroscopy measurements of cerebral cytochrome-c-oxidase. J Biomed Opt. 2016;21(9):091307.
Cooper C, Cope M, Quaresima V, Ferrari M, Nemoto E, Springett R, et al. Measurement of cytochrome oxidase redox state by near infrared spectroscopy. In: Villringer A, et al, editors. Optical imaging of brain function and metabolism 2. Springer; 1997. p. 63–73.
Food and Drug Administration. Guidance for industry and FDA staff information for manufacturers seeking marketing clearance of diagnostic ultrasound systems and transducers. Rockville, MD: FDA; 2008.
O’Brien WD Jr. Ultrasound–biophysics mechanisms. Prog Biophys Mol Biol. 2007;93(1–3):212–55.
Festing MF, Altman DG. Guidelines for the design and statistical analysis of experiments using laboratory animals. ILAR J. 2002;43(4):244–58.
Matcher S, Cooper C. Absolute quantification of deoxyhaemoglobin concentration in tissue near infrared spectroscopy. Phys Med Biol. 1994;39(8):1295.
Wray S, Cope M, Delpy DT, Wyatt JS, Reynolds EOR. Characterization of the near infrared absorption spectra of cytochrome aa3 and haemoglobin for the non-invasive monitoring of cerebral oxygenation. Biochimica et Biophysica Acta (BBA) Bioenergetics. 1988;933(1):184–92.
Kohl M, Lindauer U, Dirnagl U, Villringer A. Separation of changes in light scattering and chromophore concentrations during cortical spreading depression in rats. Opt Lett. 1998;23(7):555–7.
Tufail Y, Matyushov A, Baldwin N, Tauchmann ML, Georges J, Yoshihiro A, et al. Transcranial pulsed ultrasound stimulates intact brain circuits. Neuron. 2010;66(5):681–94.
All authors made substantial contributions to the study concept, design, and/or the acquisition of data, analysis, and/or data interpretation. In detail, EK and EA prepared the setup and performed the animal experiments. SY, JK, JYH were in charge of the ultrasound transducer and performed the ultrasound beam characterization. EK and EA prepared the manuscript. JGK led the study. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The dataset used and analyzed during the current study are available from the corresponding author upon reasonable request.
Consent for publication
Ethics approval and consent to participate
All animals were cared for in accordance with the guidelines for the care and use of laboratory animals. The protocols used in the study has been approved by the Institutional Animal Care and Use Committee of the Gwangju Institute of Science and Technology (GIST-2017-018).
This work was supported by the National Research Foundation of Korea (2015R1D1A1A02062382, 2016M3C7A1905475, 2018R1A2B6006797) and the GIST Research Institute in 2019.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Kim, E., Anguluan, E., Youn, S. et al. Non-invasive measurement of hemodynamic change during 8 MHz transcranial focused ultrasound stimulation using near-infrared spectroscopy. BMC Neurosci 20, 12 (2019). https://doi.org/10.1186/s12868-019-0493-9
- Near-infrared spectroscopy
- Cerebral hemodynamics
- Ultrasound brain stimulation