Propagation of Parametric Uncertainty in Aliev-Panfilov Model of Cardiac Excitation

Unsupervised stochastic learning and reduced order modeling for global sensitivity analysis in cardiac electrophysiology models

BACKGROUND AND OBJECTIVE

Numerical simulations in electrocardiology are often affected by various uncertainties inherited from the lack of precise knowledge regarding input values including those related to the cardiac cell model, domain geometry, and boundary or initial conditions used in the mathematical modeling. Conventional techniques for uncertainty quantification in modeling electrical activities of the heart encounter significant challenges, primarily due to the high computational costs associated with fine temporal and spatial scales. Additionally, the need for numerous model evaluations to quantify ubiquitous uncertainties increases the computational challenges even further.

METHODS

In the present study, we propose a non-intrusive surrogate model to perform uncertainty quantification and global sensitivity analysis in cardiac electrophysiology models. The proposed method combines an unsupervised machine learning technique with the polynomial chaos expansion to reconstruct a surrogate model for the propagation and quantification of uncertainties in the electrical activity of the heart. The proposed methodology not only accurately quantifies uncertainties at a very low computational cost but more importantly, it captures the targeted quantity of interest as either the whole spatial field or the whole temporal period. In order to perform sensitivity analysis, aggregated Sobol indices are estimated directly from the spectral mode of the polynomial chaos expansion.

RESULTS

We conduct Uncertainty Quantification (UQ) and global Sensitivity Analysis (SA) considering both spatial and temporal variations, rather than limiting the analysis to specific Quantities of Interest (QoIs). To assess the comprehensive performance of our methodology in simulating cardiac electrical activity, we utilize the monodomain model. Additionally, sensitivity analysis is performed on the parameters of the Mitchell-Schaeffer cell model.

CONCLUSIONS

Unlike conventional techniques for uncertainty quantification in modeling electrical activities, the proposed methodology performs at a low computational cost the sensitivity analysis on the cardiac electrical activity parameters. The results are fully reproducible and easily accessible, while the proposed reduced-order model represents a significant contribution to enhancing global sensitivity analysis in cardiac electrophysiology.

Effect of Shock Vector Orientation in Modulating and Terminating Rotors - a Simulation Study

The main treatment option for Ventricular Fibrillation (VF), especially in out-of-hospital cardiac arrests (OHCA) is defibrillation. Typically, the survival-to-discharge rates are very poor for OHCA. Existing studies have shown that rotors may be the sources of arrhythmia and ablating them could modulate or terminate VF. However, tracking rotors and ablating them is not a feasible solution in a OHCA scenario. Hence, if the sources (or rotors) can be regionally localized non-invasively and this information can be used to direct the orientation of the shock vectors, it may aid the termination of rotors and defibrillation success. In this work, using computational modeling, we present our initial results on testing the effect of shock vector orientation on modulating (or) terminating rotors. A combination of Sovilj's and Aliev Panfilov's monodomain cardiac models were used in inducing rotors and testing the effect of shock vector magnitude and direction. Based on our simulation results on an average with four experimental trials, a shock vector directed in the perpendicular direction along the axis of the rotor terminated the rotor with 16% lesser magnitude than parallel direction and 38% lesser magnitude than in oblique direction.Clinical Relevance- A rotor localization dependent defibrillation strategy may aid the defibrillation protocol procedures to improve the survival rates. Based on the four experimental trials, the results indicate shock vectors oriented perpendicular to the axis of the rotors were efficient in modulating or terminating rotors with lower magnitude than other directions.