Modal Analysis of Cerebrovascular Effects for Digital Health Integration of Neurostimulation Therapies—A Review of Technology Concepts

Transcranial electrical stimulation (tES) is increasingly recognized for its potential to modulate cerebral blood flow (CBF) and evoke cerebrovascular reactivity (CVR), which are crucial in conditions like mild cognitive impairment (MCI) and dementia. This study explores the impact of tES on the neurovascular unit (NVU), employing a physiological modeling approach to simulate the vascular response to electric fields generated by tES. Utilizing the FitzHugh–Nagumo model for neuroelectrical activity, we demonstrate how tES can initiate vascular responses such as vasoconstriction followed by delayed vasodilation in cerebral arterioles, potentially modulated by a combination of local metabolic demands and autonomic regulation (pivotal locus coeruleus). Here, four distinct pathways within the NVU were modeled to reflect the complex interplay between synaptic activity, astrocytic influences, perivascular potassium dynamics, and smooth muscle cell responses. Modal analysis revealed characteristic dynamics of these pathways, suggesting that oscillatory tES may finely tune the vascular tone by modulating the stiffness and elasticity of blood vessel walls, possibly by also impacting endothelial glycocalyx function. The findings underscore the therapeutic potential vis-à-vis blood-brain barrier safety of tES in modulating neurovascular coupling and cognitive function needing the precise modulation of NVU dynamics. This technology review supports the human-in-the-loop integration of tES leveraging digital health technologies for the personalized management of cerebral blood flow, offering new avenues for treating vascular cognitive disorders. Future studies should aim to optimize tES parameters using computational modeling and validate these models in clinical settings, enhancing the understanding of tES in neurovascular health.


Introduction
Transcranial electrical stimulation (tES) offers a promising approach to evoke CBF [1] to study cerebrovascular reactivity (CVR), which can be partly effected by the somatoautonomic responses to transcutaneous stimulation [2].Autonomic reflexes play a vital role in maintaining homeostasis as they regulate blood pressure, breathing, and tissue oxygenation.After an autonomic reflex occurs, the autonomic nervous system (ANS) transmits signals from internal organs (such as the vasculature and heart) to the central nervous system (CNS).The vagus nerve is primarily responsible for conveying these signals to regions like the medulla, pons, and hypothalamus, including during the orienting reflex to a change in the environment.The noradrenergic (NA) nucleus locus coeruleus (LC) is a brain region that becomes activated in parallel with the autonomic nervous system [3].The LC is a brainstem nucleus involved in arousal, attention, and responses to novel stimuli.LC-NA responds to biological imperatives, which can be either spontaneous reactions to unexpected salient or threatening stimuli, or conditioned responses to anticipated behaviorally relevant stimuli.So, the activation of the LC can be associated with the general orienting reflex [4]; when we encounter something new or unexpected, the LC becomes active, helping us focus our attention and respond appropriately.However, the role of this autonomic reflex including LC-NA in regulating CBF is not fully understood, and the contribution of the baroreflex in CBF regulation is debated [5].Recent studies in animal models and the mathematical analysis of cardiovascular signals in humans provided insights into the interplay between the arterial baroreceptor reflex (baroreflex) and arousal [6].Mild baroreceptor stimulation, especially under anesthesia, may inhibit cortical arousal, while significant increases or decreases in baroreflex activation induce arousal in both animals and humans under normal physiological conditions.Additionally, cardiovascular changes during autonomic arousals and transitions between wakefulness and sleep involve adjustments in the baroreflex set point and its equilibrium with central autonomic commands.While acute increases in systemic blood pressure trigger peripheral vasodilation via the baroreflex, cerebral vasculature may need to constrict to protect the blood-brain barrier [5].Here, regional differences in the autonomic outflow between systemic and cerebral blood vessels are possible.Additionally, the impact of direct (autonomic regulation) and indirect (systemic blood pressure regulation) baroreflex influences on CBF may vary [7].In hypertensive patients, sympathetic nerve activation leads to cerebral vasoconstriction, unlike in normotensive individuals [8].The complex interplay of the arterial baroreflex with other factors like the cardiac output, respiratory chemoreflex, and physiological conditions in comorbid cardiovascular disease complicates understanding the CBF regulation and CVR to tES.Despite these challenges, the autonomic regulation including the LC-NA system response and cardiac baroreflex's role in maintaining adequate CBF should be acknowledged and investigated vis-à-vis CVR dysfunction, for example, in mild cognitive impairments (MCI) [9] that increase the risk of progression to dementia [10].Here, LC degenerates early during Alzheimer's disease (AD), and the loss of norepinephrine may contribute to AD pathogenesis [11].Aside from neurons, LC terminals densely innervate brain intraparenchymal arterioles/capillaries, and norepinephrine modulates astrocyte functions.The neurovascular unit (NVU), which involves interactions between neurons, glial cells, and brain vessels, plays a fundamental role in coupling the energy demand with the regional cerebral blood flow.It includes the blood-brain barrier (BBB) and participates in neuroinflammation and the glymphatic system.NVU alteration is implicated in AD pathophysiology due to relative oligemia in activated brain regions and impaired BBB integrity, contributing to the intracerebral accumulation of insoluble amyloid.The interaction between LC and NVU in AD pathogenesis remains an area of study with therapeutic tES applications [2].
We presented a physiological modeling approach [12] based on the physiology of neurovascular tissues for assessing the vascular response to electric fields generated by tES through various pathways in NVU.The studies [2] presented a physiological model that incorporated the NVU components, including a vascular smooth muscle, perivascular space, synaptic space, and astrocyte cell.The model aimed to capture the effects of transcranial electrical stimulation (tES), specifically transcranial alternating current stimulation (tACS), on direct and indirect vascular responses.Four nested NVU compartmental pathways were proposed, allowing the simulation of tES-induced vessel volume responses.In computational analysis, the tES current density, acting as an input pulse, perturbed state variables in each NVU compartment.This study considered four simulated pathways for vessel response modulation: synaptic potassium, astrocytic membrane potential, perivascular potassium currents, and voltage-gated ion channels on smooth muscle cells.These pathways were designed to simulate vessel oscillations within the frequency range <0.2 Hz.Modal analysis [13,14], a technique commonly used in structural and fluid mechanics [14], was applied to derive the characteristic dynamics of the NVU model.Modal analysis involves determining the system's natural frequencies and damping factors, allowing the development of a mathematical model describing the system's behavior.While modal analysis is traditionally used in engineering fields, the study applied this approach to analyze the NVU model, specifically focusing on evaluating neurovascular coupling modes induced by tACS.Then, the design of the controls is imperative to modify the natural behavior of interconnected synchronous generators in NVU systems.Despite the inherent nonlinearity of NVU systems, accurately predicting oscillations around an operating point is possible through linearized system models [12], which justifies the application of linear control theory [15] for the design of tES controllers.
The examination of tES effects, both immediate and prolonged, and the design of controllers is an area of interest of neuroscientific research [12,16,17].We examined an fMRI-tES dataset [18] with a TR of 3.36 s, revealing a similar finite impulse response hemodynamic response function "FIR HRF" model (with RH, TTP, and FWHM using the rsHRF toolbox [19]) in anodal and sham tDCS conditions at the electrode locations FC5 and PZ ROIs, but differing in the tDCS electrode location FP2 ROI.This discrepancy may be linked to local cortical inhibitory circuits, perivascular nerves, and astrocyte stimulation [18].Here, prior computational analyses proposed direct perivascular nerve and astrocyte stimulation during tDCS onset leading to an "initial dip" [18].Prior findings also indicate perivascular space changes (https://www.ismrm.org/workshops/2022/Neuromodulation/program.php;accessed on 12 May 2024), which indicated vascular effects but may also raise safety concerns especially for higher intensity 4 mA stimulation.Then, tES is postulated to also impact blood-brain barrier permeability in NVU, influencing the neuronal function [20], which also raises safety concerns.Here, tDCS and tACS differ in their current profiles, with different therapeutic and safety implications.Indeed, short-duration tDCS can have physiological effects, impacting the autonomic and hemodynamic response [2].Then, tES sympathoexcitation, likely from somato-autonomic reflexes, reflected in pupil dilation [2], may impact glucose regulation via noradrenaline's impact on cellular energy processes [21] and neurometabolic coupling.Here, fNIRS's total hemoglobin (blood volume) changes due to neurovascular and neurometabolic coupling needs further investigation vis-à-vis their correlation with pupil dilation [2].In animal studies, the pupil diameter has been shown to be inversely correlated with cortical hemodynamics during rest (and non-REM sleep), while when the mouse is alert, moving, or stimulated (as well as during REM sleep), positive correlations were found between the pupil diameter and blood volume, possibly linked to the bistability [22] of norepinephrine-linked astrocytic activity [23].Therefore, the correlation analysis of pupil dilation vis-à-vis fNIRS's total hemoglobin (blood volume) changes may delineate tES effects via LC-NA [24].Also, a mechanistic understanding of the glucoseneurovascular tissue interaction during tES is crucial [2].Here, the tES sympathoexcitation (indirect effect) and direct modulation of the membrane potential of smooth muscle cells [12], particularly with specific oscillatory frequencies, is examined via frequency domain modal analysis [2].
It is postulated that orthosympathetic activation can affect pial and perforant arteries (>0.05 Hz), while LC-NA affects intraparenchymal arterioles and capillaries (<0.05 Hz) [25].Understanding these frequency domain modulatory effects of tES on blood vessels necessitates exploring multiple pathways within the NVU.Unravelling these signaling pathways is crucial for comprehending the tES effects on neurons and blood vessels to develop therapeutics, as discussed by Arora et al. [12].The stimulation-evoked CBF changes depend on multiple pathways including the sympathetic vascular tone (stimulus-related norepinephrine released from sympathetic efferent nerves modulating the vascular tone [26]), against which the neurovascular coupling need to act to dilate the blood vessels [2].Indeed, it has been shown that the long-term administration of the noradrenaline (norepinephrine) reuptake inhibitor reboxetine (RBX) extended the effects of anodal tDCS on long-term potentiation-like plasticity for over 24 h [27].RBX also transformed cathodal tDCS-induced long-term depression-like plasticity into facilitation for 120 min.Here, tES effects on autonomic innervation [2] and the activation of noradrenergic receptors can have stimulatory effects on both energy-requiring and energy-yielding processes [21].Specifically, tES electric field distribution [28] and its activating function [2] can affect the three nerve arrangements in the human cerebral arteries [29]: (a) paravascular nerve bundles outside the tunica adventitia; (b) a meshwork-like perivascular plexus in the outer or middle ad-ventitial zone; and (c) a deep intrinsic perivascular plexus at the adventitial-medial border, oriented transversely.Then, different subtypes of adrenoceptors often activate distinct processes, and the stimulation can allow for the simultaneous enhancement of the oxidative metabolism and/or glycogenolysis, e.g., beneficial in type 2 diabetes [2].Unlike classical mechanisms, tES effects may enable the stimulation of the energy metabolism without preceding decreases in ATP (for activating ATP-sensitive potassium channels) thereby leading to beneficial effects via the relationship between the metabolism, excitability, and Aβ pathology [30].The stimulation of glycogenolysis is noteworthy, as it is considered an integral part of glucose breakdown via a significant 'glucose-glycogen shunt.'While an increase in mitochondrial Ca 2+ has been observed in astrocytes, the direct stimulation of the oxidative metabolism by elevated intracellular calcium has been extensively studied, mainly in muscle and the liver.This includes the direct stimulation of mitochondrial dehydrogenases and oxidative phosphorylation, contributing to the understanding of noradrenaline's impact on cellular energy processes [21].
Indeed, neurovascular coupling itself may be tES-modulated [21,31] since the noradrenaline release from LC axons induces vessel tone in arteriolar smooth muscle and contractile capillary pericytes [32], in addition to small tES field-induced smooth muscle cell membrane potential hyperpolarization [33].Vessel tone is crucial since it enables neuronal activity to trigger vasodilation via neurovascular coupling [32], enhancing local cerebral blood flow regulation-see Figure 1a.In the brain, a significant portion of vascular resistance is in the capillaries, and LC axons release noradrenaline closer to pericytes than to arterioles.Then, the cerebral adrenoreceptor distribution exhibits heterogeneity, indicating the region-specific autonomic regulation of CBF [34].Cerebral circulation features unique compensatory responses, involving chemo-and autoregulatory mechanisms, and interactions with parasympathetic nerve activity.This interplay between sympathetic and parasympathetic reflexes ensures the optimal perfusion of CBF in response to changing perfusion pressures, aiming to optimize oxygen and nutrient delivery to the brain while maintaining the blood volume and intracranial pressure.Here, tES effects on the autonomic innervation of the cranial circulation from the external noradrenergic innervation can include a predominant sympathetic component from the superior cervical ganglion and a cranial parasympathetic component that passes through the pterygopalatine (sphenopalatine) and otic ganglion [11].Claassen and colleagues [35] discuss the mechanisms contributing to myogenic responses, focusing on the regulation of vascular tone, intravascular pressure distribution in the brain, and the autoregulation of CBF.Also, factors include the sensitivity of the vascular diameter to tES changes in the cellular membrane potential, influenced by voltage-dependent calcium channels (e.g., CaV2.1) and large conductance potassium channels (BK Ca ) [36].Then, mechanotransducers are proposed as sensors for pressure changes leading to the depolarization of vascular muscle and increased intracellular Ca 2+ .This, in turn, activates contractile proteins, resulting in vasoconstriction.Also, the local release of Ca 2+ can activate BK Ca , inducing hyperpolarization and limiting vasoconstriction [36].Therefore, the autonomic modulation of the vascular tone needs to be well-coordinated with the neurovascular coupling-related effects for an adequate CBF response to quick task-related metabolic needs (this may be compromised in type 2 diabetes [37]) that may be facilitated with tES, e.g., burst tES [38].In individuals over 60 years of age, tDCS has been shown to induce cardiovascular and autonomic improvements, enhancing ventricular repolarization dispersion timing, decreasing sympathetic activity and peripheral resistance, and increasing vagal sinus activity and baroreflex sensitivity [39].Although the authors [39] hypothesized that the anodal tDCS effects on the neuronal networks in the temporal cortex and insular cortex are linked to autonomic nervous system regulation and the perception of emotional sensations within the body, they did not compare the effects of anodal tDCS at an unrelated brain area and so the mechanism of action (e.g., autonomic reflex arc) is unclear.Another In individuals over 60 years of age, tDCS has been shown to induce cardiovascular and autonomic improvements, enhancing ventricular repolarization dispersion timing, decreasing sympathetic activity and peripheral resistance, and increasing vagal sinus activity and baroreflex sensitivity [39].Although the authors [39] hypothesized that the anodal tDCS effects on the neuronal networks in the temporal cortex and insular cortex are linked to autonomic nervous system regulation and the perception of emotional sensations within the body, they did not compare the effects of anodal tDCS at an unrelated brain area and so the mechanism of action (e.g., autonomic reflex arc) is unclear.Another study [40] investigated the impact of repetitive tDCS on cerebral perfusion-anodal left prefrontal tDCS was administered over three consecutive days, leading to widespread increases in perfusion, suggesting a heightened metabolism.In contrast, a matched group receiving sham tDCS showed general perfusion decreases.The active stimulation group exhibited significantly greater perfusion increases in various brain areas, particularly in the LC and neocortex regions associated with object recognition and attentional modulation.The observed changes in the neocortex may result directly from stimulation or indirectly through altered noradrenergic system activity from LC modulation, which can also modulate neurovascular coupling [31].Another study [24] suggests that the effects of electrical stimulation, specifically occipital nerve tDCS (ON-tDCS), involve ascending fibers of the occipital nerve synapsing with neurons in the nucleus tractus solitarius.These neurons then project to the LC, promoting NA release and enhancing functional connectivity with the hippocampus.Here, we suggest a short-duration ON-OFF tDCS time series (see Figure 1 in [41]) acting as slow transcranial oscillating-current stimulation (tOCS) [2] leveraging tissue low-pass filtering which may act through superficial nerves, noradrenergic axons, and efferent innervation to evoke a beneficial mural cell response for neurovascular coupling [42].Unlike pharmaceutical approaches, ON-tDCS offers the potential for persistent, stimulus-specific changes to neural circuits with minimal side effects.Here, tES under portable brain imaging [17] could be an individualized therapy option for conditions involving insufficient LC-NA function [11], including late-life depression leading to mild cognitive impairment [43,44], where the effects may be mediated by the individual's cerebrovascular and NVU status.

Materials and Methods
The major factors involved in the design of tES dosage are: current amplitude, waveform, polarity, duration, montage, and electrode specifications [45], as depicted in Figure 1b.These factors are crucial in neuromodulating specific characteristics.For instance, tOCS has been shown to facilitate corticospinal excitability phase independently, both on-and off-line, similar to tDCS [46].Meanwhile, tACS was more likely to entrain neuronal activity while blocking sensory input [47].To comprehend the mechanistic aspects of tES techniques on hemodynamics, we used mathematical model [12] based on the physiology of neurovascular tissue for evaluating the vascular response through various paths that are susceptible to the electric fields generated by tES, as shown in Figure 1a.Our simulation model was constructed with four compartments, drawing from the existing literature: synaptic space, astrocyte space, perivascular space, and arteriole smooth muscle cell space.To simulate the vessel volume response within the physiological model, we designed four nested neurovascular unit (NVU) compartmental pathways.Each state variable in these pathways could be influenced by the transcranial electrical stimulation (tES) current density, acting as the input pulse.The simulations of the model considered different tES-induced perturbations: synaptic potassium release from active neurons for Pathway 1, astrocytic transmembrane current for Pathway 2, perivascular potassium concentration for Pathway 3, and voltage-gated ion channel current on the smooth muscle cells (SMC) for Pathway 4. Detailed information regarding the model implementation and computational analysis can be found in publication [12].
Physiologically detailed models [12] were executed using the 'ode23tb' solver in Simulink Release 2019b (MathWorks, Inc., Natick, MA, USA).These models simulated oscillations ranging from 0 to 0.2 Hz in response to tDCS perturbations.Subsequently, we subjected the four nested NVU compartmental pathways [12] to perturbations from transcranial oscillating current stimulation (tOCS), tDCS, and transcranial alternating current stimulation (tACS) at varying frequencies (0.1 Hz to 10 Hz), conducting a sensitivity analysis for blood vessel diameter changes.This comprehensive approach considered the vascular effects of tES, incorporating both neuronal and non-neuronal mechanisms with distinct sensitivity levels.Notably, within the frequency range of 0.1 Hz to 10 Hz, we observed that vessel oscillations exhibited greater sensitivity to tOCS than to tACS, and entrainment effects were more pronounced at lower frequencies.Subsequently, modal analysis was performed on the physiological model where we applied ten random tES perturbations, utilizing bandpass-filtered (0.01-1 Hz) white noise inputs, to the four implemented physiologically detailed models [12].The input and output time series were recorded using the time-domain data object ('iddata' in MATLAB Release 2019b, MathWorks, Inc., USA).For modal analysis, we focused on the oscillatory component of the vessel response, excluding the initial 50 s of time series data.Modal analysis functions, including 'modalfrf,' 'modalfit,' and 'modalsd,' were applied to the data object to generate frequency-response functions, natural frequencies, and stabilization diagrams, respectively.The natural frequencies of the four system modes, determined from the measured frequency-response functions (frf) at frequencies (f) and a sample rate of 10 samples per second, were calculated using the peak-picking method.The peak-picking method posits that each notable peak in the frf represents a single natural mode; however, this does not ensure the orthogonality of the modes in a multi-degree of freedom (MDOF) system of nxn matrices such that it can be more simply described as a series of n decoupled single-degree of freedom (SDOF) system.The natural modes of a discrete linear system exhibit orthogonality in relation to the mass and stiffness matrices in case of structural mechanics, but this orthogonality is understood in a generalized form.Generally, natural modes are not orthogonal in the usual sense, meaning the dot product of two modal vectors typically is not zero.However, certain simple conditions can be introduced to achieve ordinary orthogonality of these modes in fluid mechanics [48] to address the loss of orthogonality issue.Then, 'modalsd' estimates the natural frequencies and damping ratios and generates the stabilization diagram using the least-squares complex exponential method that calculates the impulse response for each frequency-response function and fits it with a series of complex damped sinusoids using Prony's method.
Figure 1a outlines proposed mechanism through which tES influences the perivascular space, as theorized by Arora et al. [12].The hypothesis posits that tES modulates the vasculature via the perivascular pathway (including autonomic innervation and activation of noradrenergic receptors-see Figure 1a), leading to vasoconstriction increasing and vasodilation decreasing the perivascular space volume.In a computational model [12], an immediate vascular response was captured through the perivascular pathway, involving the interaction between perivascular potassium and calcium concentrations, resulting in steady-state vessel oscillations below 0.1 Hz.These oscillations may potentially phase synchronize with neuronal oscillations, motivating the exploration of neurovascular coupling through joint imaging with fNIRS-EEG [17].Acute tES within a short duration (<150 s), as indicated by Arora et al. [12], can impact the vasculature for the immediate control of blood vessel response using model predictive control (MPC) [2].Then, the cerebrospinal fluid (CSF) circulation involves two distinguishable components: bulk flow (circulation) and pulsatile flow (back and forth motion) that can be measured using near-infrared transillumination-backscattering sounding methods [49].The pulsatile flow in CSF is influenced by both cardiac and respiratory cycles reflected in the fNIRS oscillations [50] and there is growing interest in the arterial pulsations driving CSF flow [51,52] that needs model-based analysis [2] of the role of neurovascular coupling [53,54] via synaptic-like transmission between neural axons and arteriolar smooth muscle cells [55].MPC employs an internal model that considers cortical activity, local metabolic factors, and vascular response to optimize tES control actions over a predefined prediction horizon, operating in a receding horizon fashion for online operation.Then, rapid vascular response observed at the onset of tES [12] could modulate bulk flow of cerebrospinal fluid [56] when applied as slow oscillations according to the resonance frequencies creating an opportunity for human-in-the-loop optimization based on blood volume (total hemoglobin) feedback from fNIRS [2].Previous research [2] utilizing modal analysis and a case study in a healthy human suggested that an "optimal" oscillatory frequency warrants further investigation on the effects of bulk flow of cerebrospinal fluid [56] vis-à-vis the state of the astrocytes, extracellular glucose [57], and interstitial potassium (tES modulation of neurovascular coupling [42]).
Our recent findings from Minager et al. (https://www.ismrm.org/workshops/2022/Neuromodulation/program.php,accessed on 12 May 2024) provide support for the immediate alteration of perivascular space morphology by tES that may be linked to the vascular response [2,15].In a recent magnetic resonance imaging (MRI) study [52], the effects of low-frequency hemodynamics (0.01-0.1 Hz) on cerebrospinal fluid (CSF) movement were explored, with a focus on the correlation between vascular low-frequency oscillations (LFOs) in the neck and CSF directed toward the spinal canal.This study established a temporal relationship between the internal carotid artery and caudally directed CSF movement.The origin of vascular LFOs remained unclear, with potential sources including blood CO2 variations and vasomotion.Additionally, the study investigated the influence of respiration on CSF dynamics, revealing a sequential and continual impact on the CSF-venous system.Cross-frequency coupling suggested that respiration affects CSF dynamics through both pressure-related venous return changes and blood gas-related vessel volume changes.Although low-frequency hemodynamics and respiration contributed to CSF flow, power spectrum analysis indicated that the LFOs play a more substantial role in regulating CSF dynamics than respiration.Cardiac pulsations were found to have a comparatively weaker influence on CSF movement.With the convenience of fNIRS neuroimaging [58] and nearinfrared transillumination-backscattering sounding methods [49] in a point-of-care setting compared to MRI [52], the portable neuroimaging-based Model Predictive Control (MPC) of tES [2] is proposed to be applied during sleep to facilitate glymphatic clearance [59] in MCI to protect from AD.Then, decreased metabolite clearance [60] due to sleep deprivation can lead to decreased neurovascular coupling responses [61]; however, alterations in functional connectivity due to sleep deprivation are likely attributable to impaired transitions within brain networks during task performance, possibly related to neurovascular coupling status of different brain regions.To gain a comprehensive understanding of the mechanisms behind sleep deprivation-induced lack in neurovascular homeostasis [62] and their impact on cognitive performance in real-life situations, further longitudinal research is essential [61].Here, our innovative approach enables the optimization of tES patterns at the point of care to elicit the optimal blood volume response and CSF dynamics, potentially reducing risks associated with heightened metabolic demand with neuronal activation in pathological tissues, such as ischemia in vascular dementia.
Yashika et al. [12] employed mathematical modeling and hypothesis testing to investigate systems dynamics within four compartments of the neurovascular unit (NVU): Pathway 1-synaptic space, Pathway 2-intracellular astrocyte space, Pathway 3-perivascular space, and Pathway 4-intracellular space of the arteriolar smooth muscle cells (aSMCs).A system identification approach using a physiologically constrained linear model was used to analyze NIRS-based CVR measured during anodal high-definition tDCS in healthy individuals.Yashika et al. [12] found that the perivascular Pathway 3 exhibited the best fit with lowest mean square error and Akaike information criterion within the NVU as the primary mechanism for transient CVR at the onset of tDCS.Here, the "initial dip" or vasoconstriction at the start of the anodal tDCS was postulated to be related to afferent stimulation of nerves from the autonomic and sensory ganglia [12].Researchers debated whether perivascular nerves directly signal aSMCs in Pathway 4 via neurotransmitters and the consensus leaned towards vasodilation resulting from molecules produced by neural activity, often by astrocytes, rather than direct neurotransmission as modeled computationally [12,42].Recently, Zhang et al. [55] investigated direct neuron-vessel interactions in mice and found that the gaps in astrocyte end feet revealed dendrites and axons forming neuromuscular junctions with aSMCs.Specifically, a single glutamatergic axon dilated the arterioles they innervate through synaptic-like transmission at neural-aSMC junctions.The presynaptic bouton connected with both postsynaptic dendrites and aSMCs, and these aSMCs expressed various neuromediator receptors, including a low level of the glutamate NMDA receptor subunit 1 [55].Here, following "initial dip," the vasodilation due to tDCS [12] can be driven by various neuromediator receptors including aSMCs NMDA receptor for tDCS after effects [63].Therefore, we augmented Pathway 4 model [12] to include direct neural-aSMC interactions for disentangling the 'onset response' of tES [18].We hypothesize that following "initial dip" or vasoconstriction from pial arterial vasculature downstream to penetrating arterioles, an upstream vasodilation from aSMCs of penetrating arterioles and first-order vessels due to tES needs investigation.Here, under the assumption of expression of glutamate receptors in human aSMCs, Pathway 4 [12] can capture perivascular glutamatergic signaling at the neurovascular junctions.Then, the activation of NMDA receptors in neurovascular junctions is postulated to trigger the influx of calcium ions that bind to and activate BK channels, leading to an increased potassium ion efflux.This process induces membrane hyperpolarization, ultimately causing the relaxation of aSMCs and upstream penetrating arterioles.So, tES electric field perturbs BK current ( K 4 s/τ+1 I tdcs , K 4 constant, τ time constant) that was added to other currents including I L , I K , I Ca , and I KIR that represent leak, K + , Ca 2+ , and KIR channel currents, respectively, in the aSMC compartment from prior work [12].Then, aSMC membrane potential, V SMC , is given: Temporal interference stimulation (tIS) as well as deep transcranial magnetic stimulation (TMS) can target deep brain regions [64].Neuronal effects of tIS and TMS can be captured with FitzHugh-Nagumo model [65] that was used to model perivascular nerves of the penetrating arterioles-the bottleneck of neocortex perfusion [66].For tES, we employ the FitzHugh-Nagumo model to simulate neuronal activity using a sub-threshold stimulus, meaning the stimulus is not strong enough to trigger action potentials (spikes).The microvascular bed consists of specific territories, each supplied by a single penetrating arteriole [66], so facilitating blood flow in the penetrating arterioles with tES/TMS can significantly improve perfusion.Based on prior work [67], tES/TMS-driven perfusion was modeled with an open source fluid-structure interaction software based on Peskin's Immersed Boundary Method [68].Since complete blockade of voltage-dependent sodium channels eliminated the excitability changes [63,69], we hypothesize tES modulation of the fast dynamics (sodium current) of the FitzHugh-Nagumo model.In the default model from [68], diffusion coefficient = 10, threshold potential = 0.3, resetting rate = 1, blocking strength = 0.001, and activation strength = 0.05.Then, a Force-Length-Velocity smooth muscle model [68] was used with maximum isometric force produced at the optimum length of the muscle fibers = 100,000, length of the muscle fibers = 1, length at which the muscle fibers exert their maximum tension = 1, and constant specific for each muscle = 0.3.Vascular stiffness was represented by a damped spring mode [68] where a frictional damping force is assumed to be proportional to the velocity of the oscillation that we neglected.Here, we changed the modulus of elasticity from 0.5 × 2 × 10 8 to 2 × 2 × 10 8 Pa to study the effects of bending among three consecutive Lagrangian points, such as X1, X2, and X3 (see Figure 2a), employing a noninvariant beam that links these three successive nodes [68].Figure 2a also shows tES voltage at the top (V in (X 1 ), V in (X 2 ), V in (X 3 )) and the bottom (V out (X 1 ), V out (X 2 ), V out (X 3 )) adventitia layer (with nerves) of the axial cross-section of the blood vessel.The J T and J N denote the tangential current at those locations and all R denote the lumped electrical resistances of the ohmic model.    in the pial arterial vasculature, prompting blood flow downstream towards the penetrating arteriolesee panel A. Subsequently, tES-induced metabolic demands in the capillary bed activate the upstream penetrating arteriole within the brain parenchyma, resulting in delayed vasodilation-see panel B. We investigated the voltage perturbation required to activate the perivascular nerve using the FitzHugh-Nagumo model along the adventitia layer of the penetrating arteriole (i.e., the bottleneck [66]).The distribution of pressure and velocity in a 300-micron length of the penetrating arteriole in response to perivascular nerve stimulation over a duration of 0.01 s is shown starting from (c) to ending in (d) due to peristaltic vessel wall motion (data cursor 'Level' shows the values).The stimulation begins at the interface to pial arterial vasculature (near cortical surface) at X = 100 microns and extends to X = 400 microns into the brain parenchyma over a duration of 0.01 s, corresponding to a flow speed of 30 mm/s.Here, the arteriole's modulus of elasticity is 200 million Pascals.(e) Distribution of lumen pressure and velocity resulting from peristalsis induced by perivascular nerve stimulation in the vessel wall, illustrated under two different arteriole's modulus of elasticity.In (e), the first scenario (panels A,B) spans from time t = 0 to t = 0.01 s with arteriole's modulus of elasticity at 0.5 × 2 × 10 8 Pa.The second scenario (panels C,D) also covers from time t = 0 to t = 0.01 s, but with arteriole's modulus of elasticity at 2 × 2 × 10 8 Pa.(data cursor 'Level' shows the values).

Results
Figure 2b shows how pial tES can facilitate downstream blood flow.Figure 2b shows the hypothesis that tES "initial dip" or vasoconstriction from the pial arterial to the penetrating arteriole can force blood downstream, while the tES metabolic need-based upstream activation [70] of the penetrating arteriole leads to delayed vasodilation.Figure 2b shows the tES-evoked electropotential (mV) of the perivascular nerve of a simulated 500-micron-long mesoscopic penetrating arteriole [66] with a 100-micron diameter.Then, a mesoscopic 300-micron-long penetrating arteriole segment with a 100-micron diameter was simulated for fluid-structure-nerve interactions using the FitzHugh-Nagumo neuron model with tES-evoked V in (X 1 = 100 µm) = 1.08 mV near the cortical surface and V in (X 2 = 250 µm) = V in (X 3 = 400 µm) = 0 downstream.Figure 2c shows a FitzHugh-Nagumo model of the penetrating arteriole segment with V in (X 1 = 100 µm) = 1.08 mV at the pial-penetrating arteriole neck region, then V in (X 2 = 250 µm) = V in (X 3 = 400 µm) = 0 downstream.Also, vasodilation takes place downstream from the penetrating arteriole (X 2 = 250 µm) and roughly to the first-order vessel depth (X 3 = 400 µm), which are all represented by a single mesoscopic penetrating arteriole shown in pink in Figure 2c,d.The pressure and velocity magnitude in the lumen of the mesoscopic penetrating arteriole due to vasoconstriction ("initial dip") are shown in Figure 2c where the initially flow velocity is high near X = 100 µm (Figure 2c) and moves downstream to X = 400 µm (Figure 2d) in 0.01 s due to the peristaltic vessel wall motion which is important for 'coordination' with the upstream vasodilation-see the hypothesis in Figure 2b.The pressure and velocity magnitude in the penetrating arteriole with peristaltic vessel wall motion for two different moduli of elasticity, 0.5 × 2 × 10 8 Pa (Figure 2e(A,B)) to 2 × 2 × 10 8 Pa, are shown in Figure 2e(C,D).Initially, the flow velocity is high, near X = 100 µm (Figure 2e(A,C)), which moves downstream to X = 400 µm (Figure 2e(B,D)) in 0.01 s due to the peristaltic motion (hypothesized coordination with the upstream vasodilation-see Figure 2b).There is an increase in the flow velocity magnitude with an increase in the moduli of elasticity, which can lead to downstream turbulence affecting the endothelial and glycocalyx function [67].
In our conceptual proposition on tES for digital health integration at the point-ofcare settings, we augmented the model for Pathway 4 [12] that was motivated by Zhang et al. [55].We explored the impact of tES on perivascular axons through synaptic-like transmission at neural-aSMC junctions, leading to changes in the vessel circumference.The Pathway 4 model [12] was augmented with tES effects via the FitzHugh-Nagumo model of the perivascular axon at neural-aSMC junctions.This is relevant in vascular cognitive impairments where tES effects on the bottleneck penetrating arterioles may be therapeutic, which can be mediated by vessel stiffness in SVD.Increased vascular stiffness can also be due to elevated glucose levels in diabetes that lead to structural and functional changes in proteins and the implications of protein aggregation and advanced glycation end products [37] in the context of neurovascular coupling in dementia [25].Indeed, vessel stiffness emerges as a crucial factor in determining hemodynamic oscillations that can be probed using modal analysis [37].Then, the human-in-the-loop optimization of tES/TMS under NIRS imaging may individualize the dose-response by supporting higher oxygen availability in the capillary bed and paravascular waste clearance [2] via the optimal peristaltic motion of the penetrating arteriole [25].The current study augmented the postulated mechanism [12] for tES modulation by involving direct effects via neuron-vessel interactions [55].Computationally immersed boundary method modeling for different vessel stiffness captured an immediate vascular peristaltic response, that may lead to a turbulent flow downstream in the stiffer vessels.Oscillatory tES perturbed vessel oscillations [71] and the effects on neurovascular coupling effects can be explored using joint imaging with NIRS-EEG [72].
Here, oscillatory tES hypothesized to acutely coordinate bottleneck penetrating arteriole peristalsis to satisfy metabolic needs, which is highlighted in Figure 2b.We posit in this paper that digital health integration with neurostimulation therapies for the therapeutic use of non-invasive brain stimulation might transform the complex neuroenergetic processes and the bulk flow of cerebrospinal fluid, enhancing favorable cognitive outcomes [42] in mild cognitive impairments to early dementia, such as diminished cognitive fatigue [73], attributed to improved extracellular clearance [74].
For oscillatory modal analysis, the linear model of the four physiologically detailed tES perturbation pathways was established using the Model Linearizer tool in the Simulink (MathWorks, Inc., USA) linear analysis package [12].The damping ratio, natural frequency, and time constant of the poles were derived using the 'damp' function on the linear model system model of the four physiologically detailed tES perturbation pathways following the application of the Model Linearizer tool.Pathway 4 (see Figure 1a), from the tES perturbation effecting the voltage gated ion channels of smooth muscle cells in parts of the brain [75] to the change in the blood vessel circumference [12], i.e., the transfer function (TF4) in the Laplace domain, is TF4 = (s + 2.962)/(s 6  Note that for Pathway 4 for the tES perturbation effects, the natural frequency below 0.2 Hz is 0.04 Hz.Then, if we have the nested Pathway 3 (see Figure 1a) from the tES perturbation effecting the perivascular potassium concentration to the change in the blood vessel circumference [12], the transfer function (TF3) in the Laplace domain is TF3 = (s 2 + 2.371 × 10 7 s + 7.023 × 10 7 )/(s 8 + 9.624 × 10 6 s 7 + 2.857 × 10 11 s 6 + 1 × 10 Note that for Pathway 3 for the tES perturbation effects, the natural frequencies below 0.2 Hz are 0.04 Hz and 0.16 Hz.Then, if we have the nested Pathway 2 (see Figure 1a) from the tES perturbation effecting the astrocytic membrane potential to the change in the blood vessel circumference [12], the transfer function (TF2) in the Laplace domain is TF2 = (s 3 + 2.371 × 10 7 s 2 + 1.172 × 10 9 s + 3.262 × 10 9 )/(s 10  Note that for the Pathway 2 for the tES perturbation effects, the natural frequencies below 0.2 Hz are 0.04 Hz and 0.16 Hz.Then, if we have the nested Pathway 1 (see Figure 1a) from the tES perturbation effecting the synaptic potassium to the change in the blood vessel circumference [12], the transfer function (TF1) in the Laplace domain is TF1 = (s 3 + 2.371 × 10 7 s 2 + 1.172 × 10 9 s + 3.262 × 10 9 )/(s 11  Note that for the Pathway 1 for the tES perturbation effects, the natural frequencies below 0.2 Hz are 0.04 Hz, 0.06 Hz, and 0.16 Hz.In order to test the modal analysis functions, including 'modalfrf', 'modalfit', and 'modalsd', using simulated data from the non-linear models of four physiologically detailed tES perturbation pathways [12], the input of transcranial electrical stimulation (tES) and the corresponding vessel responses for the proposed pathways were investigated.We subjected the four compartmental pathways to tOCS (combined tDCS and tACS) perturbations with varying frequencies (ranging from 0.1 Hz to 10 Hz) and direct current (DC) offsets (ranging from 0 to 2 mA).Subsequently, we conducted a sensitivity analysis to assess changes in the blood vessel circumference.
We considered three model input waveforms as follows.
tES inputs : In Equation ( 1), we considered a sinusoidal amplitude (a) of 1 mA, a sinusoidal frequency (f ) ranging from 0.1 to 10 Hz, and a DC offset (c) ranging from 0 to 2 times the amplitude for conducting a sensitivity analysis on the physiologically detailed mathematical model of the neurovascular unit.This analysis, facilitated by the Sensitivity Analyzer tool in MATLAB Simulink (MathWorks, Inc., USA), enabled the exploration of the tES design space, identifying the most influential model parameters.Figure 3 illustrates the impact of frequency and DC offset on the vessel response with tOCS @ 0.1 Hz performing the best across all the pathways.Then, Pathway 4 was most responsive to tACS and tOCS @ 1 Hz.Pathways 1 and 4 were entrained with a positive correlation while Pathways 2 and 3 were entrained with a negative correlation for the tACS part of the tOCS-see Figure 3.
Brain Sci.2024, 14, x FOR PEER REVIEW 15 of 24 proposed pathways were investigated.We subjected the four compartmental pathways to tOCS (combined tDCS and tACS) perturbations with varying frequencies (ranging from 0.1 Hz to 10 Hz) and direct current (DC) offsets (ranging from 0 to 2 mA).Subsequently, we conducted a sensitivity analysis to assess changes in the blood vessel circumference.We considered three model input waveforms as follows. : In Equation ( 1), we considered a sinusoidal amplitude (a) of 1 mA, a sinusoidal frequency (f) ranging from 0.1 to 10 Hz, and a DC offset (c) ranging from 0 to 2 times the amplitude for conducting a sensitivity analysis on the physiologically detailed mathematical model of the neurovascular unit.This analysis, facilitated by the Sensitivity Analyzer tool in MATLAB Simulink (MathWorks, Inc., USA), enabled the exploration of the tES design space, identifying the most influential model parameters.Figure 3 illustrates the impact of frequency and DC offset on the vessel response with tOCS @ 0.1 Hz performing the best across all the pathways.Then, Pathway 4 was most responsive to tACS and tOCS @ 1 Hz.Pathways 1 and 4 were entrained with a positive correlation while Pathways 2 and 3 were entrained with a negative correlation for the tACS part of the tOCS-see Figure 3. Table 1 presents the natural frequencies derived through the modal analysis of the non-linear models of four physiologically detailed tES perturbation pathways under a Bandpass-filtered White Noise Input and Figure 4 illustrates a boxplot representation of these frequencies.Additionally, Figure 5 displays the stabilization diagram obtained for the first case, where the model input was bandpass-filtered white noise with the default seed value for all the four pathways.Table 1 lists the system parameters associated with a linearized model for the four pathways-Pathway 1 involves tES perturbing the vessel response through the synaptic potassium pathway, Pathway 2 through the astrocytic pathway, Pathway 3 through the perivascular potassium pathway, and Pathway 4 through the smooth muscle cell pathway.Figure 5 showed the natural frequencies of the non-linear models of four physiologically detailed tES perturbation pathways where we Table 1 presents the natural frequencies derived through the modal analysis of the non-linear models of four physiologically detailed tES perturbation pathways under a Bandpass-filtered White Noise Input and Figure 4 illustrates a boxplot representation of these frequencies.Additionally, Figure 5 displays the stabilization diagram obtained for the first case, where the model input was bandpass-filtered white noise with the default seed value for all the four pathways.Table 1 lists the system parameters associated with a linearized model for the four pathways-Pathway 1 involves tES perturbing the vessel response through the synaptic potassium pathway, Pathway 2 through the astrocytic pathway, Pathway 3 through the perivascular potassium pathway, and Pathway 4 through the smooth muscle cell pathway.Figure 5 showed the natural frequencies of the non-linear models of four physiologically detailed tES perturbation pathways where we can find stable modes ~0.05 Hz for all the pathways that were compatible with the linearized model analysis (0.04 Hz, 0.06 Hz).
Table 1.Natural frequencies (Hz) obtained from the modal analysis (see also Figure 4) of the four physiologically detailed tES perturbation pathways for ten different seed values for the Bandpassfiltered White Noise Input (1* presents the results from the default seed value).

Discussion
Our results demonstrate the potential of tES in manipulating perivascular nerve activity, thereby modulating blood flow dynamics in cerebral arterioles.Specifically, our findings from Figure 2b support the hypothesis that tES can induce initial vasoconstriction ("initial dip") in pial arteries, which facilitates downstream blood flow enhancement.This effect is likely complemented by a delayed vasodilation response due to metabolic demands in upstream penetrating arterioles, a phenomenon that may be particularly beneficial in managing vascular cognitive impairments.Moreover, our study incorporates the FitzHugh-Nagumo model to simulate the electrophysiological behavior of perivascular nerves along penetrating arterioles, highlighting how tES can influence the vascular diameter through synaptic-like interactions at neural-aSMC junctions.These interactions may be critical in conditions characterized by increased vascular stiffness, such as small vessel disease (SVD) seen in diabetes.Elevated glucose levels in diabetes contribute to structural changes in vessel proteins, which can affect overall neurovascular coupling-a key factor in dementia progression.Our computational modeling also suggests that variations in the modulus of the elasticity of vessel walls can significantly affect hemodynamic responses, such as the flow velocity and turbulence, potentially impacting endothelial and glycocalyx functions.Such insights are crucial for understanding the mechanistic underpinnings of tES in therapeutic contexts, particularly its role in enhancing neurovascular coupling and cognitive function.Indeed, a hallmark of dementia syndromes is an energy deficit in neurons, leading to synaptic loss and resulting in cognitive decline and behavioral changes [76]; therefore, an energy boost is necessary [77], e.g., via facilitating neurovascular [78] and neurometabolic [79] coupling rather than solely excitatory tDCS that may increase the energy deficit leading to synaptic loss.Also, Gamma entrainment may positively impact excitation/inhibition imbalances by preventing neurochemical changes linked to Aβ, such as hyperexcitability, and triggering neuroprotective mechanisms [80].However, we should be aware of the energy cost to avoid a deficit (and cell death) due to Gamma entrainment since fast-spiking interneurons, with their membrane potentials around 10 to 15 mV, which is more depolarized than those of pyramidal cells, may incur higher energy costs due to the need for chloride extrusion [81].Additionally, the dynamic balance between excitatory and inhibitory synaptic inputs can challenge ion homeostasis.For example, an increase in inhibition typically coincides with an increase in excitation, maintaining excitability within a functional range.However, in high membrane conductance states, this balance leads to heightened ion fluxes for transmission, which can significantly increase the energy demands on both inhibitory interneurons and excitatory pyramidal cells during periods of intense synaptic activity due to Gamma entrainment.
There are various neurostimulation approaches to facilitate CBF: Animal studies have revealed the presence of spontaneous "hypoxic pockets" in the brains of awake, active mice, linked to interruptions in the local capillary flow [82].Moreover, it was found that exercise reduced these hypoxic areas by 52% compared to when at rest.This research offers new understanding of the oxygen distribution in the brain and introduces interventions to target oxygen levels in physiological functions and cerebrovascular disorders.Here, peripheral neuromuscular electrical stimulation (NMES) can have beneficial cerebral hemodynamic effects [83] that can be adjuvant to exercise interventions where peripheral NMES effects can be partly due to the somato-autonomic reflex.We postulate, for multi-level non-invasive electrical stimulation interventions, that the orthosympathetic and LC-NA system effects need to be delineated since they can have opposite effects on CBF.Also, different forms of tES have varying current profiles for tES-transcranial direct current stimulation (tDCS) uses a monophasic, non-oscillating baseline, while transcranial alternating current stimulation (tACS) involves rhythmically reversing the electrical current.Other methods include transcranial oscillating current stimulation (tOCS), using a direct component to guide oscillations, and transcranial random noise stimulation (tRNS), injecting an alternating current with bounded [84,85].Neurovascular modulation occurs in various stimulation protocols, with mechanisms not fully understood that are postulated to be partly driven by somato-autonomic [86] and the LC-NA [12] effects.As tES affects blood vessels through neuronal or non-neuronal cells, a deeper understanding of signaling pathways is crucial [12].Here, tACS, unique in manipulating cell-class specific [87] intrinsic oscillations possibly through ephaptic coupling [88], may hold promise for addressing vascular dementia [89,90] and related microvascular dysfunction.In the systematic review by Machado and colleagues [91] and meta-analysis of 22 studies involving 393 participants, the effects of transcranial Direct Current Stimulation (tDCS) on exercise performance were examined.Weak evidence suggested a significant positive effect of anodal tDCS (a-tDCS) over the motor cortex (M1) on time to exhaustion (TTE) in cycling, but the results were influenced by a single study.No significant effects were found for cathodal tDCS (c-tDCS) on TTE.For isometric muscle strength, no significant effects were observed for a-tDCS applied before or during exercise.Mixed results were reported for isokinetic muscle strength.A quantitative synthesis indicated a significant improvement in cycling performance with a-tDCS over M1, but caution was advised due to the influence of a single study.Commercial tDCS devices for exercise performance were not addressed, raising safety and efficacy concerns.Methodological aspects, including individual variability and optimal tDCS parameters, need further exploration in future research.Here, tES has a role in improving exercise performance when individually customized for priming the NVU for optimal hemodynamic responses in a closed-loop manner [2] while addressing detrimental cerebral effects of the metaboreflex [92,93] using peripheral NMES.Human in the loop approach [2] via ephaptic coupling [88] may leverage the synaptic transmission between the neural axons and the smooth muscle cells in arterioles [55] thereby modulating the neurovascular pathway [12].
Limitations: In our mathematical model, we did not present a spatio-temporal tES approach to facilitate long-wavelength traveling waves of vasomotion, which is a major limitation and needs to be addressed in the future.Brain arterioles, oscillating at approximately 0.1 Hz, are active multicellular complexes [94].In awake mice, these vaso-oscillations in penetrating arterioles significantly impact neocortex perfusion, with resting-state activity modulating blood flow more than stimulus-induced activity.The weak relationship between arteriole diameter changes and perfusion suggests that the capillary bed predominantly controls brain vasculature resistance.The phase of these oscillations evolves slowly along arterioles, creating cortical areas with uniform phases which may be measured using whole-head multi-distance fNIRS with spatio-temporal oscillatory modal analysis [95] and spatiotemporal dynamic mode decomposition [96]-see the Supplementary Materials for the theoretical development [48,97] using coupled complex Ginzburg-Landau equations [98][99][100].Here, phase-gradient supports bidirectional traveling waves, indicating that penetrating arteriole waves can mix, but not directional transport of interstitial fluids, which needs further investigation in future works.

Conclusions
This study provides important insights into the vascular mechanisms activated by tES and their potential therapeutic implications for neurovascular and cognitive dysfunctions.The ability of tES to induce coordinated peristaltic motion in penetrating arterioles could be a significant factor in enhancing cerebral blood flow and, consequently, cognitive functions in patients with vascular impairments.Then, the integration of tES with digital health technologies at the point-of-care could revolutionize the management of cognitive impairments in early dementia stages, potentially reducing cognitive fatigue and enhancing extracellular clearance through improved neurovascular coupling.Future research should focus on refining these neurostimulation strategies, optimizing their parameters, and validating their effectiveness in clinical settings to fully harness their therapeutic potential.These findings open new avenues for the use of neuromodulation in treating vascular and neurodegenerative diseases, emphasizing the need for further investigation into the optimal integration of these technologies into clinical practice.

Figure 1 .
Figure 1.(a) Effect of tES dosage on neurovascular tissue: physiological modeling depiction.(b) tES time series perturbations for model evaluation.tDCS has a monophasic, non-oscillating baseline, while tACS rhythmically reverses electron flow.Additional methods include tOCS, guiding oscillations with a direct component, and tRNS, injecting bounded stochastic alternating current.

Figure 1 .
Figure 1.(a) Effect of tES dosage on neurovascular tissue: physiological modeling depiction.(b) tES time series perturbations for model evaluation.tDCS has a monophasic, non-oscillating baseline, while tACS rhythmically reverses electron flow.Additional methods include tOCS, guiding oscillations with a direct component, and tRNS, injecting bounded stochastic alternating current.

Figure 2 .
Figure 2. (a) The picture displays the transcranial electrical stimulation (tES) voltage at t (Vin(X1), Vin(X2), Vin(X3)) and bottom (Vout(X1), Vout(X2), Vout(X3)) adventitia layer of a vessel's axial cross-section, along with JT and JN denoting tangential current and R repres lumped electrical resistances in the ohmic model.Investigation of the fluid-structure inte due to smooth muscle cell's tES activation with altered modulus of elasticity of bending among three consecutive Lagrangian points (X1, X2, and X3) along the penetrating arteriole u

Figure 2 .
Figure 2. (a) The picture displays the transcranial electrical stimulation (tES) voltage at the top (Vin(X1), Vin(X2), Vin(X3)) and bottom (Vout(X1), Vout(X2), Vout(X3)) adventitia layer of a blood vessel's axial cross-section, along with JT and JN denoting tangential current and R representing lumped electrical resistances in the ohmic model.Investigation of the fluid-structure interaction due to smooth muscle cell's tES activation with altered modulus of elasticity of bending effects among three consecutive Lagrangian points (X1, X2, and X3) along the penetrating arteriole using a noninvariant beam connecting these nodes.(b) We propose that tES initially triggers vasoconstriction

Figure 3 .
Figure 3.Effect of frequency (tDCS has frequency 0) and DC offset of tOCS on vessel circumference based sensitivity analysis for the proposed pathways.

Figure 3 .
Figure 3.Effect of frequency (tDCS has frequency 0) and DC offset of tOCS on vessel circumference based sensitivity analysis for the proposed pathways.

Figure 4 .Figure 4 .Figure 5 .
Figure 4. Boxplot illustrates the distribution of natural frequencies resulting from modal analysis (with different seed values) for the four transcranial electrical stimulation (tES) perturbation model pathways that are <0.1 Hz.Each box represents the interquartile range, with the central mark denoting the median.The bottom and top edges of the box indicate the 25th and 75th percentiles, respectively.Whiskers extend to the most extreme data points that are not considered outliers, while outliers are individually marked with a red '+' symbol.

Table 1 .
find stable modes ~0.05 Hz for all the pathways that were compatible with the linearized model analysis (0.04 Hz, 0.06 Hz).Natural frequencies (Hz) obtained from the modal analysis (see also Figure4) of the four physiologically detailed tES perturbation pathways for ten different seed values for the Bandpassfiltered White Noise Input (1* presents the results from the default seed value). can