Adaptive Parameter Selection for Deep Brain Stimulation in Parkinson’s Disease

Each year, around 60,000 people are diagnosed with Parkinson’s disease (PD) and the economic burden of PD is at least $14.4 billion a year in the United States. Pharmaceutical costs for a Parkinson’s patient can be reduced from$12,000 to $6,000 per year with the addition of neuromodulation therapies such as Deep Brain Stimulation (DBS), transcranial Direct Current Stimulation (tDCS), Transcranial Magnetic Stimulation (TMS), etc. In neurodegenerative disorders such as PD, deep brain stimulation (DBS) is a desirable approach when the medication is less effective for treating the symptoms. DBS incorporates transferring electrical pulses to a specific tissue of the central nervous system and obtaining therapeutic results by modulating the neuronal activity of that region. The hyperkinetic symptoms of PD are associated with the ensembles of interacting oscillators that cause excess or abnormal synchronous behavior within the Basal Ganglia (BG) circuitry. Delayed feedback stimulation is a closed loop technique shown to suppress this synchronous oscillatory activity. Deep Brain Stimulation via delayed feedback is known to destabilize the complex intermittent synchronous states. Computational models of the BG network are often introduced to investigate the effect of delayed feedback high frequency stimulation on partially synchronized dynamics. In this work, we developed several computational models of four interacting nuclei of the BG as well as considering the Thalamo-Cortical local effects on the oscillatory dynamics. These models are able to capture the emergence of 34 Hz beta band oscillations seen in the Local Field Potential (LFP) recordings of the PD state. Traditional High Frequency Stimulations (HFS) has shown deficiencies such as strengthening the synchronization in case of highly fluctuating neuronal activities, increasing the energy consumed as well as the incapability of activating all neurons in a large-scale network. To overcome these drawbacks, we investigated the effects of the stimulation waveform and interphase delays on the overall efficiency and efficacy of DBS. We also propose a new feedback control variable based on the filtered and linearly delayed LFP recordings. The proposed control variable is then used to modulate the frequency of the stimulation signal rather than its amplitude. In strongly coupled networks, oscillations reoccur as soon as the amplitude of the stimulus signal declines. Therefore, we show that maintaining a fixed amplitude and modulating the frequency might ameliorate the desynchronization process, increase the battery lifespan and activate substantial regions of the administered DBS electrode. The charge balanced stimulus pulse itself is embedded with a delay period between its charges to grant robust desynchronization with lower amplitude needed. The efficiency and efficacy of the proposed Frequency Adjustment Stimulation (FAS) protocol in a delayed feedback method might contribute to further investigation of DBS modulations aspired to address a wide range of abnormal oscillatory behaviors observed in neurological disorders. Adaptive stimulation can open doors towards simultaneous stimulation with MRI recordings. We additionally propose a new pipeline to investigate the effect of Transcranial Magnetic Stimulation (TMS) on patient specific models. The pipeline allows us to generate a full head segmentation based on each individual MRI data. In the next step, the neurosurgeon can adaptively choose the proper location of stimulation and transmit accurate magnetic field with this pipeline

Brain and Human Body Modelling 2021

This open access book describes modern applications of computational human modelling to advance neurology, cancer treatment, and radio-frequency studies including regulatory, safety, and wireless communication fields. Readers working on any application that may expose human subjects to electromagnetic radiation will benefit from this book’s coverage of the latest models and techniques available to assess a given technology’s safety and efficacy in a timely and efficient manner. This is an Open Access book

Brain and Human Body Modelling 2021

This open access book describes modern applications of computational human modelling to advance neurology, cancer treatment, and radio-frequency studies including regulatory, safety, and wireless communication fields. Readers working on any application that may expose human subjects to electromagnetic radiation will benefit from this book’s coverage of the latest models and techniques available to assess a given technology’s safety and efficacy in a timely and efficient manner. This is an Open Access book

A high-performance boundary element method and its applications in engineering

As a semi-numerical and semi-analytical method, owing to the inherent advantage, of boundary-only discretisation, the boundary element method (BEM) has been widely applied to problems with complicated geometries, stress concentration problems, infinite domain problems, and many others. However, domain integrals and non-symmetrical and dense matrix systems are two obstacles for BEM which have hindered the its further development and application. This thesis is aimed at proposing a high-performance BEM to tackle the above two drawbacks and broaden the application scope of BEM. In this thesis, a detailed introduction to the traditional BEM is given and several popular algorithms are introduced or proposed to enhance the performance of BEM. Numerical examples in heat conduction analysis, thermoelastic analysis and thermoelastic fracture problems are performed to assess the efficiency and correction of the algorithms. In addition, necessary theoretical derivations are embraced for establishing novel boundary integral equations (BIEs) for specific engineering problems. The following three parts are the main content of this thesis. (1) The first part (Part II consisting of two chapters) is aimed at heat conduction analysis by BEM. The coefficient matrix of equations formed by BEM in solving problems is fully-populated which occupy large computer memory. To deal with that, the fast multipole method (FMM) is introduced to energize the line integration boundary element method (LIBEM) to performs better in efficiency. In addition, to compute domain integrals with known or unknown integrand functions which are caused by heat sources or heterogeneity, a novel BEM, the adaptive orthogonal interpolation moving least squares (AOIMLS) method enhanced LIBEM, which also inherits the advantage of boundary-only discretisation, is proposed. Unlike LIBEM, which is an accurate and stable method for computing domain integrals, but only works when the mathematical expression of integral function in domain integrals is known, the AOIMLS enhanced LIBEM can compute domain integrals with known or unknown integral functions, which ensures all the nonlinear and nonhomogeneous problems can be solved without domain discretisation. In addition, the AOIMLS can adaptively avoid singular or ill-conditioned moment matrices, thus ensuring the stability of the calculation results. (2) In the second part (Part III consisting of four chapters), the thermoelastic problems and fracture problems are the main objectives. Due to considering thermal loads, domain integrals appear in the BIEs of the thermoelastic problems, and the expression of integrand functions is known or not depending on the temperature distribution given or not, the AOIMLS enhanced LIBEM is introduced to conduct thermoelasticity analysis thereby. Besides, a series of novel unified boundary integral equations based on BEM and DDM are derived for solving fracture problems and thermoelastic fracture problems in finite and infinite domains. Two sets of unified BIEs are derived for fracture problems in finite and infinite domains based on the direct BEM and DDM respectively, which can provide accurate and stable results. Another two sets of BIEs are addressed by employing indirect BEM and DDM, which cannot ensure a stable result, thereby a modified indirect BEM is proposed which performs much more stable. Moreover, a set of novel BIEs based on the direct BEM and DDM for cracked domains under thermal stress is proposed. (3) In the third part (Part IV consisting of one chapter), a high-efficiency combined BEM and discrete element method (DEM) is proposed to compute the inner stress distribution and particle breakage of particle assemblies based on the solution mapping scheme. For the stress field computation of particles with similar geometry, a template particle is used as the representative particle, so that only the related coefficient matrices of one template particle in the local coordinate system are needed to be calculated, while the coefficient matrices of the other particles, can be obtained by mapping between the local and global coordinate systems. Thus, the combined BEM and DEM is much more effective when modelling a large-scale particle system with a small number of distinct possible particle shapes. Furthermore, with the help of the Hoek-Brown criterion, the possible cracks or breakage paths of a particle can be obtained