Location-Aware Range-Error Correction for Improved UWB Localization.

In this paper, we present a novel localization scheme, location-aware ranging correction (LARC), to correct ranging estimates from ultra wideband (UWB) signals. Existing solutions to calculate ranging corrections rely solely on channel information features (e.g., signal energy, maximum amplitude, estimated range). We propose to incorporate a preliminary location estimate into a localization chain, such that location-based features can be calculated as inputs to a range-error prediction model. This way, we can add information to range-only measurements without relying on additional hardware such as an inertial measurement unit (IMU). This improves performance and reduces overfitting behavior. We demonstrate our LARC method using an open-access measurement dataset with distances up to 20 m, using a simple regression model that can run purely on the CPU in real-time. The inclusion of the proposed features for range-error mitigation decreases the ranging error 90th percentile (P90) by 58% to 15 cm (compared to the uncorrected range error), for an unseen trajectory. The 2D localization P90 error is improved by 21% to 18 cm. We show the robustness of our approach by comparing results to a changed environment, where metallic objects have been moved around the room. In this modified environment, we obtain a 56% better P90 ranging performance of 16 cm. The 2D localization P90 error improves as much as for the unchanged environment, by 17% to 18 cm, showing the robustness of our method. This method evolved from the first-ranking solution of the 2021 and 2022 International Conference on Indoor Position and Indoor Navigation (IPIN) Competition.

Robust IMU-Based Mitigation of Human Body Shadowing in UWB Indoor Positioning.

Ultra-wideband (UWB) indoor positioning systems have the potential to achieve sub-decimeter-level accuracy. However, the ranging performance degrades significantly under non-line-of-sight (NLoS) conditions. The detection and mitigation of NLoS conditions is a complex problem and has been the subject of many works over the past decades. When localizing pedestrians, human body shadowing (HBS) is a particular and specific cause of NLoS. In this paper, we present an HBS mitigation strategy based on the orientation of the body and tag relative to the UWB anchors. Our HBS mitigation strategy involves a robust range error model that interacts with a tracking algorithm. The model consists of a bank of Gaussian Mixture Models (GMMs), from which an appropriate GMM is selected based on the relative body-tag-anchor orientation. The relative orientation is estimated by means of an inertial measurement unit (IMU) attached to the tag and a candidate position provided by the tracking algorithm. The selected GMM is used as a likelihood function for the tracking algorithm to improve localization accuracy. Our proposed approach was realized for two tracking algorithms. We validated the implemented algorithms on dynamic UWB ranging measurements, which were performed in an industrial lab environment. The proposed algorithms outperform other state-of-the-art algorithms, achieving a 37% reduction of the p75 error.

Optimal Frequency and Wireless Power Budget for Miniature Receivers in Obese People.

This study investigates wireless power transfer for deep in-body receivers, determining the optimal frequency, power budget, and design for the transmitter and receiver. In particular, the focus is on small, in-body receivers at large depths up to 20 cm for obese patients. This enables long-term monitoring of the gastrointestinal tract for all body types. Numerical simulations are used to investigate power transfer and losses as a function of frequency and to find the optimal design at the selected frequency for an obese body model. From all ISM-frequencies in the investigated range (1 kHz-10 GHz), the value of 13.56 MHz yields the best performance. This optimum corresponds to the transition from dominant copper losses in conductors to dominant losses in conductive tissue. At this frequency, a transmitting and receiving coil are designed consisting of 12 and 23 windings, respectively. With a power transfer efficiency of 2.70×10-5, 18 µW can be received for an input power of 0.68 W while still satisfying exposure guidelines. The power transfer is validated by measurements. For the first time, efficiency values and the power budget are reported for WPT through 20 cm of tissue to mm sized receivers. Compared to WPT at higher frequencies, as commonly used for small receivers, the proposed system is more suitable for WPT to large depths in-body and comes with the advantage that no focusing is required, which can accommodate multiple receivers and uncertainty about receiver location more easily. The received power allows long-term sensing in the gastrointestinal tract by, e.g., temperature, pressure, and pH sensors, motility sensing, or even gastric stimulation.

Vehicle Localization Using Doppler Shift and Time of Arrival Measurements in a Tunnel Environment.

Most applications and services of Cooperative Intelligent Transport Systems (C-ITS) rely on accurate and continuous vehicle location information. The traditional localization method based on the Global Navigation Satellite System (GNSS) is the most commonly used. However, it does not provide reliable, continuous, and accurate positioning in all scenarios, such as tunnels. Therefore, in this work, we present an algorithm that exploits the existing Vehicle-to-Infrastructure (V2I) communication channel that operates within the LTE-V frequency band to acquire in-tunnel vehicle location information. We propose a novel solution for vehicle localization based on Doppler shift and Time of Arrival measurements. Measurements performed in the Beveren tunnel in Antwerp, Belgium, are used to obtain results. A comparison between estimated positions using Extended Kalman Filter (EKF) on Doppler shift measurements and individual Kalman Filter (KF) on Doppler shift and Time of Arrival measurements is carried out to analyze the filtering methods performance. Findings show that the EKF performs better than KF, reducing the average estimation error by 10 m, while the algorithm accuracy depends on the relevant RF channel propagation conditions and other in-tunnel-related environment knowledge included in the estimation. The proposed solution can be used for monitoring the position and speed of vehicles driving in tunnel environments.

Numerically simulated exposure of children and adults to pulsed gradient fields in MRI.

PURPOSE: To determine exposure to gradient switching fields of adults and children in a magnetic resonance imaging (MRI) scanner by evaluating internal electric fields within realistic models of adult male, adult female, and child inside transverse and longitudinal gradient coils, and to compare these results with compliance guidelines. MATERIALS AND METHODS: Patients inside x-, y-, and z-gradient coils were simulated using anatomically realistic models of adult male, adult female, and child. The induced electric fields were computed for 1 kHz sinusoidal current with a magnitude of 1 A in the gradient coils. Rheobase electric fields were then calculated and compared to the International Commission on Non-Ionizing Radiation Protection (ICNIRP) 2004 and International Electrotechnical Commission (IEC) 2010 guidelines. The effect of the human body, coil type, and skin conductivity on the induced electric field was also investigated. RESULTS: The internal electric fields are within the first level controlled operating mode of the guidelines and range from 2.7V m-1 to 4.5V m-1 , except for the adult male inside the y-gradient coil (induced field reaches 5.4V m-1 ).The induced electric field is sensitive to the coil type (electric field in the skin of adult male: 4V m-1 , 4.6V m-1 , and 3.8V m-1 for x-, y-, and z-gradient coils, respectively), the human body model (electric field in the skin inside y-gradient coil: 4.6V m-1 , 4.2V m-1 , and 3V m-1 for adult male, adult female, and child, respectively), and the skin conductivity (electric field 2.35-4.29% higher for 0.1S m-1 skin conductivity compared to 0.2S m-1 ). CONCLUSION: The y-gradient coil induced the largest fields in the patients. The highest levels of internal electric fields occurred for the adult male model. J. Magn. Reson. Imaging 2016;44:1360-1367.

Analysis of eddy currents induced by transverse and longitudinal gradient coils in different tungsten collimators geometries for SPECT/MRI integration.

PURPOSE: We investigated the temporal variation of the induced magnetic field due to the transverse and the longitudinal gradient coils in tungsten collimators arranged in hexagonal and pentagonal geometries with and without gaps between the collimators. METHODS: We modeled x-, y-, and z-gradient coils and different arrangements of single-photon emission computed tomography (SPECT) collimators using FEKO, a three-dimensional electromagnetic simulation tool. A time analysis approach was used to generate the pulsed magnetic field gradient. The approach was validated with measurements using a 7T MRI scanner. RESULTS: Simulations showed an induced magnetic field representing 4.66% and 0.87% of the applied gradient field (gradient strength = 500 mT/m) for longitudinal and transverse gradient coils, respectively. These values can be reduced by 75% by adding gaps between the collimators for the pentagonal arrangement, bringing the maximum induced magnetic field to less than 2% of the applied gradient for all of the gradient coils. CONCLUSION: Characterization of the maximum induced magnetic field shows that by adding gaps between the collimators for an integrated SPECT/MRI system, eddy currents can be corrected by the MRI system to avoid artifact. The numerical model was validated and was proposed as a tool for studying the effect of a SPECT collimator within the MRI gradient coils.

In-to-out body antenna-independent path loss model for multilayered tissues and heterogeneous medium.

In this paper, we investigate multilayered lossy and heterogeneous media for wireless body area networks (WBAN) to develop a simple, fast and efficient analytical in-to-out body path loss (PL) model at 2.45 GHz and, thus, avoid time-consuming simulations. The PL model is an antenna-independent model and is validated with simulations in layered medium, as well as in a 3D human model using electromagnetic solvers.

Simulation study on high spatio-temporal resolution acousto-electrophysiological neuroimaging.

Objective.Acousto-electrophysiological neuroimaging (AENI) is a technique hypothesized to record electrophysiological activity of the brain with millimeter spatial and sub-millisecond temporal resolution. This improvement is obtained by tagging areas with focused ultrasound (fUS). Due to mechanical vibration with respect to the measuring electrodes, the electrical activity of the marked region will be modulated onto the ultrasonic frequency. The region's electrical activity can subsequently be retrieved via demodulation of the measured signal. In this study, the feasibility of this hypothesized technique is tested.Approach.This is done by calculating the forward electroencephalography response under quasi-static assumptions. The head is simplified as a set of concentric spheres. Two sizes are evaluated representing human and mouse brains. Moreover, feasibility is assessed for wet and dry transcranial, and for cortically placed electrodes. The activity sources are modeled by dipoles, with their current intensity profile drawn from a power-law power spectral density.Results.It is shown that mechanical vibration modulates the endogenous activity onto the ultrasonic frequency. The signal strength depends non-linearly on the alignment between dipole orientation, vibration direction and recording point. The strongest signal is measured when these three dependencies are perfectly aligned. The signal strengths are in the pV-range for a dipole moment of 5 nAm and ultrasonic pressures within Food and Drug Administration (FDA)-limits. The endogenous activity can then be accurately reconstructed via demodulation. Two interference types are investigated: vibrational and static. Depending on the vibrational interference, it is shown that millimeter resolution signal detection is possible also for deep brain regions. Subsequently, successful demodulation depends on the static interference, that at MHz-range has to be sub-picovolt.Significance.Our results show that mechanical vibration is a possible underlying mechanism of acousto-electrophyisological neuroimaging. This paper is a first step towards improved understanding of the conditions under which AENI is feasible.

Validation of experimental whole-body SAR assessment method in a complex indoor environment.

Experimentally assessing the whole-body specific absorption rate (SAR(wb) ) in a complex indoor environment is very challenging. An experimental method based on room electromagnetics theory (accounting only the line-of-sight as specular path) is validated using numerical simulations with the finite-difference time-domain method. Furthermore, the method accounts for diffuse multipath components (DMC) in the total absorption rate by considering the reverberation time of the investigated room, which describes all the losses in a complex indoor environment. The advantage of the proposed method is that it allows discarding the computational burden because it does not use any discretizations. Results show good agreement between measurement and computation at 2.8 GHz, as long as the plane wave assumption is valid, that is, at large distances from the transmitter. Relative deviations of 0.71% and 4% have been obtained for far-field scenarios, and 77.5% for the near field-scenario. The contribution of the DMC in the total absorption rate is also quantified here, which has never been investigated before. It is found that the DMC may represent an important part of the total absorption rate; its contribution may reach up to 90% for certain scenarios in an indoor environment.

Personal distributed exposimeter for radio frequency exposure assessment in real environments.

For the first time, a personal distributed exposimeter (PDE) for radio frequency (RF) measurements is presented. This PDE is designed based on numerical simulations and is experimentally evaluated using textile antennas and wearable electronics. A prototype of the PDE is calibrated in an anechoic chamber. Compared to conventional exposimeters, which only measure in one position on the body, an excellent isotropy of 0.5 dB (a factor of 1.1) and a 95% confidence interval of 7 dB (a factor of 5) on power densities are measured.

Quantitative analysis of the optogenetic excitability of CA1 neurons.

Introduction: Optogenetics has emerged as a promising technique for modulating neuronal activity and holds potential for the treatment of neurological disorders such as temporal lobe epilepsy (TLE). However, clinical translation still faces many challenges. This in-silico study aims to enhance the understanding of optogenetic excitability in CA1 cells and to identify strategies for improving stimulation protocols. Methods: Employing state-of-the-art computational models coupled with Monte Carlo simulated light propagation, the optogenetic excitability of four CA1 cells, two pyramidal and two interneurons, expressing ChR2(H134R) is investigated. Results and discussion: The results demonstrate that confining the opsin to specific neuronal membrane compartments significantly improves excitability. An improvement is also achieved by focusing the light beam on the most excitable cell region. Moreover, the perpendicular orientation of the optical fiber relative to the somato-dendritic axis yields superior results. Inter-cell variability is observed, highlighting the importance of considering neuron degeneracy when designing optogenetic tools. Opsin confinement to the basal dendrites of the pyramidal cells renders the neuron the most excitable. A global sensitivity analysis identified opsin location and expression level as having the greatest impact on simulation outcomes. The error reduction of simulation outcome due to coupling of neuron modeling with light propagation is shown. The results promote spatial confinement and increased opsin expression levels as important improvement strategies. On the other hand, uncertainties in these parameters limit precise determination of the irradiance thresholds. This study provides valuable insights on optogenetic excitability of CA1 cells useful for the development of improved optogenetic stimulation protocols for, for instance, TLE treatment.

Influence of Temporal Interference Stimulation Parameters on Point Neuron Excitability.

Temporal interference (TI) stimulation is a technique in which two high frequency sinusoidal electric fields, oscillating at a slightly different frequency are sent into the brain. The goal is to achieve stimulation at the place where both fields interfere. This study uses a simplified version of the Hodgkin - Huxley model to analyse the different parameters of the TI-waveform and how the neuron reacts to this waveform. In this manner, the underlying mechanism of the reaction of the neuron to a TI -signal is investigated. Clinical relevance- This study shows the importance of the parameter choice of the temporal interference waveform and provides insights into the underlying mechanism of the neuronal response to a beating sine for the application of temporal interference stimulation.

Membrane Charge Oscillations During Ultrasonic Neuromodulation by Intramembrane Cavitation.

OBJECTIVE: To investigate the importance of membrane charge oscillations and redistribution in multi-compartmental ultrasonic neuromodulation (UNMOD) intramembrane cavitation models. METHODS: The Neuronal Intramembrane Cavitation Excitation (NICE) model and multiScale Optimized model of Neuronal Intramembrane Cavitation (SONIC) of UNMOD are compared for a nanoscale multi-compartmental and point neuron approximation of the bilayer sonophore and surrounding proteins. The temporal dynamics of charge oscillations and their effect on the resulting voltage oscillations are investigated by fourier series analysis. RESULTS: Comparison of excitation thresholds and neuronal response between nanoscale multi-compartmental and point models, implemented in the SONIC and NICE framework, demonstrates that the explicit modeling of fast spatial charge redistribution is critical for an accurate multi-compartmental UNMOD-model. Furthermore, the importance of modeling partial protein coverage is quantified by the excitability thresholds. Subsequently, we establish by fourier analysis that these charge oscillations are slowly changing in time. CONCLUSION: Fast charge redistribution significantly alters neuronal excitability in a multi-compartmental nanoscale UNMOD-model. Also the mutual exclusivity between protein and sonophore coverage should be taken into account, when simulating the dependency of neuronal excitability on coverage fractions. Charge oscillations are periodic and their fourier components change on a slow timescale. Furthermore, the resulting voltage oscillations decrease in energy with overtone number, implying that an extension of the existing multiscale model (SONIC) to multi-compartmental neurons is possible by taking into account a limited number of fourier components. SIGNIFICANCE: First steps are taken towards a morphologically realistic and computationally efficient UNMOD-model, improving our understanding of the underlying ultrasonic neuromodulation mechanisms.

Improved alpha-beta power reduction via combined electrical and ultrasonic stimulation in a parkinsonian cortex-basal ganglia-thalamus computational model.

Objective. To investigate computationally the interaction of combined electrical and ultrasonic modulation of isolated neurons and of the parkinsonian cortex-basal ganglia-thalamus loop.Approach. Continuous-wave or pulsed electrical and ultrasonic neuromodulation is applied to isolated Otsuka plateau-potential generating subthalamic nucleus (STN) and Pospischil regular, fast and low-threshold spiking cortical cells in a temporally alternating or simultaneous manner. Similar combinations of electrical/ultrasonic waveforms are applied to a parkinsonian biophysical cortex-basal ganglia-thalamus neuronal network. Ultrasound-neuron interaction is modelled respectively for isolated neurons and the neuronal network with the NICE and SONIC implementations of the bilayer sonophore underlying mechanism. Reduction inα-ßspectral energy is used as a proxy to express improvement in Parkinson's disease by insonication and electrostimulation.Main results. Simultaneous electro-acoustic stimulation achieves a given level of neuronal activity at lower intensities compared to the separate stimulation modalities. Conversely, temporally alternating stimulation with50 Hzelectrical and ultrasound pulses is capable of eliciting100 HzSTN firing rates. Furthermore, combination of ultrasound with hyperpolarizing currents can alter cortical cell relative spiking regimes. In the parkinsonian neuronal network, continuous-wave and pulsed ultrasound reduce pathological oscillations by different mechanisms. High-frequency pulsed separated electrical and ultrasonic deep brain stimulation (DBS) reduce pathologicalα-ßpower by entraining STN-neurons. In contrast, continuous-wave ultrasound reduces pathological oscillations by silencing the STN. Compared to the separated stimulation modalities, temporally simultaneous or alternating electro-acoustic stimulation can achieve higher reductions inα-ßpower for the same safety contraints on electrical/ultrasonic intensity.Significance. Focused ultrasound has the potential of becoming a non-invasive alternative of conventional DBS for the treatment of Parkinson's disease. Here, we elaborate on proposed benefits of combined electro-acoustic stimulation in terms of improved dynamic range, efficiency, spatial resolution, and neuronal selectivity.

Double Two-State Opsin Model With Autonomous Parameter Inference.

Optogenetics has a lot of potential to become an effective neuromodulative therapy for clinical applications. Selecting the correct opsin is crucial to have an optimal optogenetic tool. With computational modeling, the neuronal response to the current dynamics of an opsin can be extensively and systematically tested. Unlike electrical stimulation where the effect is directly defined by the applied field, the stimulation in optogenetics is indirect, depending on the selected opsin's non-linear kinetics. With the continuous expansion of opsin possibilities, computational studies are difficult due to the need for an accurate model of the selected opsin first. To this end, we propose a double two-state opsin model as alternative to the conventional three and four state Markov models used for opsin modeling. Furthermore, we provide a fitting procedure, which allows for autonomous model fitting starting from a vast parameter space. With this procedure, we successfully fitted two distinctive opsins (ChR2(H134R) and MerMAID). Both models are able to represent the experimental data with great accuracy and were obtained within an acceptable time frame. This is due to the absence of differential equations in the fitting procedure, with an enormous reduction in computational cost as result. The performance of the proposed model with a fit to ChR2(H134R) was tested, by comparing the neural response in a regular spiking neuron to the response obtained with the non-instantaneous, four state Markov model (4SB), derived by Williams et al. (2013). Finally, a computational speed gain was observed with the proposed model in a regular spiking and sparse Pyramidal-Interneuron-Network-Gamma (sPING) network simulation with respect to the 4SB-model, due to the former having two differential equations less. Consequently, the proposed model allows for computationally efficient optogenetic neurostimulation and with the proposed fitting procedure will be valuable for further research in the field of optogenetics.

SECONIC: Towards multi-compartmental models for ultrasonic brain stimulation by intramembrane cavitation.

OBJECTIVE: To design a computationally efficient model for ultrasonic neuromodulation (UNMOD) of morphologically realistic multi-compartmental neurons based on intramembrane cavitation. APPROACH: A Spatially Extended Neuronal Intramembrane Cavitation model that accurately predicts observed fast Charge Oscillations (SECONIC) is designed. A regular spiking cortical Hodgkin-Huxley type nanoscale neuron model of the bilayer sonophore and surrounding proteins is used. The accuracy and computational efficiency of SECONIC is compared with the Neuronal Intramembrane Cavitation Excitation (NICE) and multiScale Optimized model of Neuronal Intramembrane Cavitation (SONIC). MAIN RESULTS: Membrane charge redistribution between different compartments should be taken into account via fourier series analysis in an accurate multi-compartmental UNMOD-model. Approximating charge and voltage traces with the harmonic term and first two overtones results in reasonable goodness-of-fit, except for high ultrasonic pressure (adjusted R-squared ≥0.61). Taking into account the first eight overtones results in a very good fourier series fit (adjusted R-squared ≥0.96) up to 600 kPa. Next, the dependency of effective voltage and rate parameters on charge oscillations is investigated. The two-tone SECONIC-model is one to two orders of magnitude faster than the NICE-model and demonstrates accurate results for ultrasonic pressure up to 100 kPa. SIGNIFICANCE: Up to now, the underlying mechanism of UNMOD is not well understood. Here, the extension of the bilayer sonophore model to spatially extended neurons via the design of a multi-compartmental UNMOD-model, will result in more detailed predictions that can be used to validate or falsify this tentative mechanism. Furthermore, a multi-compartmental model for UNMOD is required for neural engineering studies that couple finite difference time domain simulations with neuronal models. Here, we propose the SECONIC-model, extending the SONIC-model by taking into account charge redistribution between compartments.

FIELD ENHANCEMENT AND SIZE OF RADIO-FREQUENCY HOTSPOTS INDUCED BY MAXIMUM RATIO FIELD COMBINING IN FIFTH GENERATION NETWORK.

The goal of this paper is to experimentally assess the field enhancement and hotspot size of radio frequency electromagnetic fields created by the Maximum Ratio Combining (MRC) precoding scheme using lab measurements at 3.5, 5.5 and 11 GHz. MRC is an adaptive precoding scheme used by Massive Multiple Input Multiple Output systems, one of the enabling techniques of the fifth generation of telecommunications (5G). A virtual antenna array was used to compare MRC with two passive precoding schemes: the Random Phase Model (RPM) and the Centerline Beam Model (CBM). The field enhancement going from CBM to MRC was largest in obstructed line of sight (OLOS), ranging from 1.9 to 7.4 dB. The field enhancement going from RPM to MRC was about 9.5 dB across frequency bands in both line of sight (LOS) and OLOS. The hotspot size, quantified by the full width at half maximum (FWHM), ranged from 0.5 wavelengths to one wavelength.

Sensitivity Study of Neuronal Excitation and Cathodal Blocking Thresholds of Myelinated Axons for Percutaneous Auricular Vagus Nerve Stimulation.

OBJECTIVE: Excitation of myelinated nerve fibers is investigated by means of numerical simulations, for the application of percutaneous auricular vagus nerve stimulation (pVNS). High sensitivity to axon diameter is of interest regarding the goal of targeting thicker fibers. METHODS: Excitation and blocking thresholds for different pulse types, phase durations, axon depths, axon-electrode distances, temperatures and axon diameters are investigated. The used model consists of a 50 mm long axon and a centrally located needle electrode in a layered medium representing the auricle. Neuronal excitation is simulated using the Frankenhaeuser-Huxley equations for all combinations of parameter values. RESULTS AND CONCLUSION: Multiple modes and locations of excitation along the axon were observed, depending on the pulse type and amplitude. When increasing the axon-electrode distance from 1 mm to 2 mm, sensitivity of thresholds to axon depth decreased with ca. 50%, while sensitivity to axon-electrode distance, axon diameter and phase duration each increased with ca. 15% to 20%, except from monophasic anodal pulses, showing a 45% decrease for axon-electrode distance. These trends for axon diameter and axon-electrode distance allow for more selective stimulation of thicker target fibers using monophasic anodal pulses at higher axon-electrode distances. Cathodal monophasic pulses did not perform well due to blocking of the thicker fibers, which was only rarely seen for other pulse types. SIGNIFICANCE: Sensitivities of stimulation thresholds to these parameters by numerical simulation reveal how the stimulation parameters can be changed in order to increase therapeutic effect and comfort during pVNS by enabling more selective stimulation.

Stimulation Pattern Efficiency in Percutaneous Auricular Vagus Nerve Stimulation: Experimental Versus Numerical Data.

OBJECTIVE: Percutaneous electrical stimulation of the auricular vagus nerve (pVNS) is an electroceutical technology. The selection of stimulation patterns is empirical, which may lead to under-stimulation or over-stimulation. The objective is to assess the efficiency of different stimulation patterns with respect to individual perception and to compare it with numerical data based on in-silico ear models. METHODS: Monophasic (MS), biphasic (BS) and triphasic stimulation (TS) patterns were tested in volunteers. Different clinically-relevant perception levels were assessed. In-silico models of the human ear were created with embedded fibers and vessels to assess different excitation levels. RESULTS: TS indicates experimental superiority over BS which is superior to MS while reaching different perception levels. TS requires about 57% and 35% of BS and MS magnitude, respectively, to reach the comfortable perception. Experimental thresholds decrease from non-bursted to bursted stimulation. Numerical results indicate a slight superiority of BS and TS over MS while reaching different excitation levels, whereas the burst length has no influence. TS yields the highest number of asynchronous action impulses per stimulation symbol for the used tripolar electrode set-up. CONCLUSION: The comparison of experimental and numerical data favors the novel TS pattern. The analysis separates excitatory pVNS effects in the auricular periphery, as accounted by in-silico data, from the combination of peripheral and central pVNS effects in the brain, as accounted by experimental data. SIGNIFICANCE: The proposed approach moves from an empirical selection of stimulation patterns towards efficient and optimized pVNS settings.

Computational Modeling of Ultrasonic Subthalamic Nucleus Stimulation.

OBJECTIVE: To explore the potential of ultrasonic modulation of plateau-potential generating subthalamic nucleus neurons (STN), by modeling their interaction with continuous and pulsed ultrasonic waves. METHODS: A computational model for ultrasonic stimulation of the STN is created by combining the Otsuka-model with the bilayer sonophore model. The neuronal response to continuous and pulsed ultrasonic waves is computed in parallel for a range of frequencies, duty cycles, pulse repetition frequencies, and intensities. RESULTS: Ultrasonic intensity in continuous-wave stimulation determines the firing pattern of the STN. Three observed spiking modes in order of increasing intensity are low frequency spiking, high frequency spiking with significant spike-frequency and spike-amplitude adaptation, and a silenced mode. Continuous-wave stimulation has little capability to manipulate the saturated spiking rate in the high frequency spiking mode. In contrast, STN firing rates induced by pulsed ultrasound insonication will saturate to the pulse repetition frequency with short latencies, for sufficiently large intensity and repetition frequency. CONCLUSION: Computational results show that the activity of plateau-potential generating STN can be modulated by selection of the stimulus parameters. Low intensities result in repetitive firing, while higher intensities silence the STN. Pulsed ultrasonic stimulation results in a shorter saturation latency and is able to modulate spiking rates. SIGNIFICANCE: Stimulation or suppresion of the STN is important in the treatment of Parkinson's disease, e.g., in deep brain stimulation. This explorative study on ultrasonic modulation of the STN, could be a step in the direction of minimally invasive alternatives to conventional deep brain stimulation.