Regional electric field induced by electroconvulsive therapy in a realistic finite element head model: influence of white matter anisotropic conductivity

We present the first computational study investigating the electric field (E-field) strength generated by various electroconvulsive therapy (ECT) electrode configurations in specific brain regions of interest (ROIs) that have putative roles in the therapeutic action and/or adverse side effects of ECT. This study also characterizes the impact of the white matter (WM) conductivity anisotropy on the E-field distribution. A finite element head model incorporating tissue heterogeneity and WM anisotropic conductivity was constructed based on structural magnetic resonance imaging (MRI) and diffusion tensor MRI data. We computed the spatial E-field distributions generated by three standard ECT electrode placements including bilateral (BL), bifrontal (BF), and right unilateral (RUL) and an investigational electrode configuration for focal electrically administered seizure therapy (FEAST). The key results are that (1) the median E-field strength over the whole brain is 3.9, 1.5, 2.3, and 2.6 V/cm for the BL, BF, RUL, and FEAST electrode configurations, respectively, which coupled with the broad spread of the BL E-field suggests a biophysical basis for observations of superior efficacy of BL ECT compared to BF and RUL ECT; (2) in the hippocampi, BL ECT produces a median E-field of 4.8 V/cm that is 1.5-2.8 times stronger than that for the other electrode configurations, consistent with the more pronounced amnestic effects of BL ECT; and (3) neglecting the WM conductivity anisotropy results in E-field strength error up to 18% overall and up to 39% in specific ROIs, motivating the inclusion of the WM conductivity anisotropy in accurate head models. This computational study demonstrates how the realistic finite element head model incorporating tissue conductivity anisotropy provides quantitative insight into the biophysics of ECT, which may shed light on the differential clinical outcomes seen with various forms of ECT, and may guide the development of novel stimulation paradigms with improved risk/benefit ratio.     

EEG source analysis of epileptiform activity using a 1 mm anisotropic hexahedra finite element head model

The major goal of the evaluation in presurgical epilepsy diagnosis for medically intractable patients is the precise reconstruction of the epileptogenic foci, preferably with non-invasive methods. This paper evaluates whether surface electroencephalography (EEG) source analysis based on a 1 mm anisotropic finite element (FE) head model can provide additional guidance for presurgical epilepsy diagnosis and whether it is practically feasible in daily routine. A 1 mm hexahedra FE volume conductor model of the patient's head with special focus on accurately modeling the compartments skull, cerebrospinal fluid (CSF) and the anisotropic conducting brain tissues was constructed using non-linearly co-registered T1-, T2- and diffusion-tensor-magnetic resonance imaging data. The electrodes of intra-cranial EEG (iEEG) measurements were extracted from a co-registered computed tomography image. Goal function scan (GFS), minimum norm least squares (MNLS), standardized low resolution electromagnetic tomography (sLORETA) and spatio-temporal current dipole modeling inverse methods were then applied to the peak of the averaged ictal discharges EEG data. MNLS and sLORETA pointed to a single center of activity. Moving and rotating single dipole fits resulted in an explained variance of more than 97%. The non-invasive EEG source analysis methods localized at the border of the lesion and at the border of the iEEG electrodes which mainly received ictal discharges. Source orientation was towards the epileptogenic tissue. For the reconstructed superficial source, brain conductivity anisotropy and the lesion conductivity had only a minor influence, whereas a correct modeling of the highly conducting CSF compartment and the anisotropic skull was found to be important. The proposed FE forward modeling approach strongly simplifies meshing and reduces run-time (37 ms for one forward computation in the model with 3.1 million unknowns), corroborating the practical feasibility of the approach.     

An extended Kalman filtering approach for the estimation of human head tissue conductivities by using EEG data: a simulation study

In this study, we propose an extended Kalman filter approach for the estimation of the human head tissue conductivities in vivo by using electroencephalogram (EEG) data. Since the relationship between the surface potentials and conductivity distribution is nonlinear, the proposed algorithm first linearizes the system and applies extended Kalman filtering. By using a three-compartment realistic head model obtained from the magnetic resonance images of a real subject, a known dipole assumption and 32 electrode positions, the performance of the proposed method is tested in simulation studies and it is shown that the proposed algorithm estimates the tissue conductivities with less than 1% error in noiseless measurements and less than 5% error when the signal-to-noise ratio is 40 dB or higher. We conclude that the proposed extended Kalman filter approach successfully estimates the tissue conductivities in vivo.     

Detection of bridging veins rupture and subdural haematoma onset using a finite element head model

BackgroundOne of the most severe traumatic brain injuries, the subdural haematoma, is related to damage and rupture of the bridging veins, generating an abnormal collection of blood between the dura mater and arachnoid mater. Current numerical models of these vessels rely on very simple geometries and material laws, limiting its accuracy and bio-fidelity.MethodsIn this work, departing from an existing human head numerical model, a realistic geometry for the bridging veins was developed, devoting special attention to the finite elements type employed. A novel and adequate constitutive model including damage behavior was also successfully implemented.FindingsResults attest that vessel tearing onset was correctly captured, after comparison against experiments on cadavers.InterpretationDoing so, the model allow to precisely predict the individual influence of kinematic parameters such as the pulse duration, linear and rotational accelerations in promoting vessel tearing.     

Multimodal integration of high-resolution EEG and functional magnetic resonance imaging data: a simulation study

Previous simulation studies have stressed the importance of the use of fMRI priors in the estimation of cortical current density. However, no systematic variations of signal-to-noise ratio (SNR) and number of electrodes were explicitly taken into account in the estimation process. In this simulation study we considered the utility of including information as estimated from fMRI. This was done by using as the dependent variable both the correlation coefficient and the relative error between the imposed and the estimated waveforms at the level of cortical region of interests (ROI). A realistic head and cortical surface model was used. Factors used in the simulations were the different values of SNR of the scalp-generated data, the different inverse operators used to estimated the cortical source activity, the strengths of the fMRI priors in the fMRI-based inverse operators, and the number of scalp electrodes used in the analysis. Analysis of variance results suggested that all the considered factors significantly afflict the correlation and the relative error between the estimated and the simulated cortical activity. For the ROIs analyzed with simulated fMRI hot spots, it was observed that the best estimation of cortical source currents was performed with the inverse operators that used fMRI information. When the ROIs analyzed do not present fMRI hot spots, both standard (i.e., minimum norm) and fMRI-based inverse operators returned statistically equivalent correlation and relative error values.     

Vibration induced refrigeration and energy harvesting using piezoelectric materials: a finite element study

In this study, the bi-functional performance of a small-scale piezoelectric cantilever, which coupled piezoelectric and elastocaloric phenomena in a single device to produce energy harvesting as well as refrigeration effects due to vibration, has been investigated. Finite element modeling has been used to examine the performance of the device. The basic structure of the device is a cantilever that vibrates between two thermal bodies (hot and cold). The properties of BaTiO₃ (single crystal) were used to examine the bi-functional performance of piezoelectric cantilevers. In this study, different cases have been investigated, which are based on a number of cantilevers between hot and cold thermal bodies. When the number of cantilevers is one, the net cooling is 0.3 K and the power is 0.03 μW, while for four cantilevers, the net cooling is 1.2 K and 0.13 μW of power is produced. The results show that as we increase the number of cantilevers, a greater refrigeration effect is produced and higher power across the electrical load is achieved.

A device utilizing ambient vibrations to produce the required cooling effect as well as energy harvesting through elastocaloric and piezoelectric effects is demonstrated.     

Transcranial direct current stimulation in patients after decompressive craniectomy: a finite element model to investigate factors affecting the cortical electric field

ObjectiveTo simulate the process of transcranial direct current stimulation (tDCS) on patients after decompressive craniectomy (DC), and to model cortical electric field distributions under different electrode montages, we constructed a finite element model that represented the human head at high resolution.MethodsUsing computed tomography images, we constructed a human head model with high geometrical similarity. The removed bone flap was simplified to be circular with a diameter of 12 cm. We then constructed finite element models according to bioelectrical parameters. Finally, we simulated tDCS on the finite element models under different electrode montages.ResultsInward current had a linear relationship with peak electric field value, but almost no effect on electric field distribution. If the anode was not over the skull hole (configuration 2), there was almost no difference in electric field magnitude and focality between the circular and square electrodes. However, if the anode was right over the hole (configuration 1), the circular electrodes led to higher peak electric field values and worse focality. In addition, configuration 1 significantly decreased focality compared with configuration 2.ConclusionOur results might serve as guidelines for selecting current and electrode montage settings when performing tDCS on patients after DC.     

A three-dimensional,anatomically detailed foot model: a foundation for a finite element simulation and means of quantifying foot-bone position

We generated an anatomically detailed, three-dimensional (3-D) reconstruction of a human foot from 286 computerized topographic (CT) images. For each bone, 2-D cross-sectional data were obtained and aligned to form a stacked image model. We calculated the inertial matrix of each bone from the stacked image model and used it to determine the principal axes. Relative angles between the principal axes of the bones were employed to describe the shape of the foot, i.e., the relationships between the bones of the foot. A 3-D surface model was generated from the stacked image models and a detailed 3-D mesh for each bone was created. Additionally, the representative geometry of the plantar soft tissue was obtained from the CT scans, while the geometries of the cartilage between bones were obtained from the 3-D surface bone models. This model served dual purposes: it formed the anatomical foundation for a future finite element model of the human foot and we used it to objectively quantify foot shape using the relationships between the principal axes of the foot bones.