Mapping of the electrical activity of human atria. Multiscale modelling and simulations

La fibrilación auricular es una de las arritmias cardíacas más comunes observadas en la práctica clínica. Por lo tanto, es de vital importancia desarrollar nuevas tecnologías destinadas a diagnosticar y acabar con este tipo de arritmia, para mejorar la calidad de vida de los pacientes y reducir los costes de los sistemas nacionales de salud. En los últimos años ha aumentado el uso de las nuevas técnicas de mapeo auricular, basadas en sistemas multi-electrodo para mapear la actividad eléctrica en humanos. Dichas técnicas permiten localizar y ablacionar los impulsores de la fibrilación auricular, como son las fuentes focales o los rotores. Sin embargo, todavía existe incertidumbre sobre su precisión y los procedimientos experimentales para su análisis están limitados debido a su carácter invasivo. Por lo tanto, las simulaciones computacionales son una herramienta muy útil para superar estas limitaciones, al permitir reproducir con fidelidad las observaciones experimentales, dividir el problema bajo estudio en sub-estudios más simples, y realizar investigaciones preliminares imposibles de llevar a cabo en el práctica clínica. Esta tesis doctoral se centra en el análisis de la precisión de los sistemas de mapeo multi-electrodo a través de modelos y simulaciones computacionales. Para ello, desarrollamos modelos realistas multi-escala con el objetivo de simular actividad eléctrica auricular reentrante, en primer lugar en una lámina de tejido auricular, y en segundo lugar en las aurículas completas. Posteriormente, analizamos los efectos de las configuraciones geométricas multi-electrodo en la precisión de la localización de los rotores, mediante el uso de agrupaciones multi-electrodo con distancias inter-electrodo equidistantes, así como a través de catéteres de tipo basket con distancias inter-electrodo no equidistantes. Después de calcular los electrogramas unipolares intracavitarios, realizamos mapas de fase, detecciones de singularidad de fase para rastrear los rotores, y mapas de frecuencia dominantes. Finalmente, descubrimos que la precisión de los sistemas de mapeo multi-electrodo depende de su posición dentro de la cavidad auricular, de la distancia entre los electrodos y el tejido, de la distancia inter-electrodo, y de la contribución de las fuentes de campo lejano. Además, como consecuencia de estos factores que pueden afectar a la precisión de los sistemas de mapeo multi-electrodo, observamos la aparición de rotores falsos que podrían contribuir al fracaso de los procesos de ablación de la fibrilación auricular.Atrial fibrillation is one of the most common cardiac arrhythmias seen in clinical practice. Therefore, it is of vital importance to develop new technologies aimed at diagnosing and terminating this kind of arrhythmia, to improve the quality of life of patients and to reduce costs to national health systems. In the last years, new atrial mapping techniques based on multi-electrode systems are increasingly being used to map the atrial electrical activity in humans and localise and target atrial fibrillation drivers in the form of focal sources or rotors. However, significant concerns remain about their accuracy and experimental approaches to analyse them are limited due to their invasive character. Therefore, computer simulations are a helpful tool to overcome these limitations since they can reproduce with fidelity experimental observations, permit to split the problem to treat into more simple substudies, and allow the possibility of performing preliminary investigations impossible to carry out in the clinical practice. This PhD thesis is focused on the analysis for accuracy of the multielectrode mapping systems through computational models and simulations. For this purpose, we developed realistic multiscale models in order to simulate atrial electrical reentrant activity, first in a sheet of atrial tissue and, then, in the whole atria. Then, we analysed the effects of the multi-electrode geometrical configurations on the accuracy of localizing rotors, by using multi-electrode arrays with equidistant inter-electrode distances, as well as multi-electrode basket catheters with non-equidistant inter-electrode distances. After computing the intracavitary unipolar electrograms, we performed phase maps, phase singularity detections to track rotors, and dominant frequency maps. We finally found out that the accuracy of multi-electrode mapping systems depends on their position inside the atrial cavity, the electrode-to-tissue distance, the inter-electrode distance, and the contribution of far field sources. Furthermore, as a consequence of these factors, false rotors might appear and could contribute to failure of atrial fibrillation ablation procedures.La fibril·lació auricular és una de les arítmies cardíaques més comuns observades en la pràctica clínica. Per tant, és de vital importància desenvolupar noves tecnologies destinades a diagnosticar i acabar amb aquest tipus d'arítmia, per tal de millorar la qualitat de vida dels pacients i reduir els costos dels sistemes nacionals de salut. En els últims anys, ha augmentat l'ús de les noves tècniques de mapeig auricular, basades en sistemes multielèctrode per a mapejar l'activitat elèctrica auricular en humans. Aquestes tècniques permeten localitzar i ablacionar els impulsors de la fibril·lació auricular, com són les fonts focals o els rotors. No obstant això, encara hi ha incertesa sobre la seua precisió i els procediments experimentals per al seu anàlisi estan limitats a causa del seu caràcter invasiu. Per tant, les simulacions computacionals són una eina molt útil per a superar aquestes limitacions, en permetre reproduir amb fidelitat les observacions experimentals, dividir el problema sota estudi en subestudis més simples, i realitzar investigacions preliminars impossibles de dur a terme en el pràctica clínica. Aquesta tesi doctoral es centra en l'anàlisi de la precisió del sistemes de mapeig multielèctrode mitjançant els models i les simulacions computacionals. Per a això, desenvolupàrem models realistes multiescala per tal de simular activitat elèctrica auricular reentrant, en primer lloc en una làmina de teixit auricular, i en segon lloc a les aurícules completes. Posteriorment, analitzàrem els efectes de les configuracions geomètriques multielèctrode en la precisió de la localització dels rotors, mitjançant l'ús d'agrupacions multielèctrode amb distàncies interelèctrode equidistants, així com catèters de tipus basket amb distàncies interelèctrode no equidistants. Després de calcular els electrogrames unipolars intracavitaris, vam realitzar mapes de fase, deteccions de singularitat de fase per a rastrejar els rotors, i mapes de freqüència dominants. Finalment, vam descobrir que la precisió dels sistemes de mapeig multielèctrode depèn de la seua posició dins de la cavitat auricular, de la distància entre els elèctrodes i el teixit, de la distància interelèctrode, i de la contribució de les fonts de camp llunyà. A més, com a conseqüència d'aquests factors, es va observar l'aparició de rotors falsos que podrien contribuir al fracàs de l'ablació de la fibril·lació auricular.Martínez Mateu, L. (2018). Mapping of the electrical activity of human atria. Multiscale modelling and simulations [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/104604TESI

Research and Education in Computational Science and Engineering

Over the past two decades the field of computational science and engineering (CSE) has penetrated both basic and applied research in academia, industry, and laboratories to advance discovery, optimize systems, support decision-makers, and educate the scientific and engineering workforce. Informed by centuries of theory and experiment, CSE performs computational experiments to answer questions that neither theory nor experiment alone is equipped to answer. CSE provides scientists and engineers of all persuasions with algorithmic inventions and software systems that transcend disciplines and scales. Carried on a wave of digital technology, CSE brings the power of parallelism to bear on troves of data. Mathematics-based advanced computing has become a prevalent means of discovery and innovation in essentially all areas of science, engineering, technology, and society; and the CSE community is at the core of this transformation. However, a combination of disruptive developments---including the architectural complexity of extreme-scale computing, the data revolution that engulfs the planet, and the specialization required to follow the applications to new frontiers---is redefining the scope and reach of the CSE endeavor. This report describes the rapid expansion of CSE and the challenges to sustaining its bold advances. The report also presents strategies and directions for CSE research and education for the next decade.Comment: Major revision, to appear in SIAM Revie

Studies on the dynamics of chaotic multi-wavelet reentrant propagation using a hybrid cellular automaton model of excitable tissue

There is a compelling body of evidence implicating continuous propagation (reentry) sustained by multiple meandering wavelets in the pathology of advanced human atrial fibrillation (AF). This forms the basis for many current therapies such as the Cox MAZE procedure and its derivatives, which aim to create non-conducting lesions in order to "transect" these circuits before they form. Nevertheless, our ability to successfully treat persistent and permanent AF using catheter ablation remains inadequate due to current limitations of clinical mapping technology as well as an incomplete understanding of how to place lesions in order to maximize circuit transection and, more importantly, minimize AF burden. Here, we used a hybrid cellular automaton model to study the dynamics of chaotic, multi-wavelet reentry (MWR) in excitable tissue. First, we used reentry as an exemplar to investigate a hysteretic disease mechanism in a multistable nonlinear system. We found that certain interactions with the environment can cause persistent changes to system behavior without altering its structure or properties, thus leading to a disconnect between clinical symptoms and the underlying state of disease. Second, we developed a novel analytical method to characterize the spatiotemporal dynamics of MWR. We identified a heterogeneous spatial distribution of reentrant pathways that correlated with the spatial distribution of cell activation frequencies. Third, we investigated the impact of topological and geometrical substrate alterations on the dynamics of MWR. We demonstrated a multi-phasic relationship between obstacle size and the fate of individual episodes. Notably, for a narrow range of sizes, obstacles appeared to play an active role in rapidly converting MWR to stable structural reentry. Our studies indicate that reentrant-pathway distributions are non-uniform in heterogeneous media (such as the atrial myocardium) and suggest a clinically measurable correlate for identifying regions of high circuit density, supporting the feasibility of patient-specific targeted ablation. Moreover, we have elucidated the key mechanisms of interaction between focal obstacles and MWR, which has implications for the use of spot ablation to treat AF as some recent studies have suggested

Computer assisted optimization of cardiac resynchronization therapy

The efficacy of cardiac resynchronization therapy (CRT) through biventricular pacing (BVP) has been demonstrated by numerous studies in patients suffering from congestive heart failure. In order to achieve a guideline for optimal treatment with BVP devices, an automated non-invasive strategy based on an electrophysiological computer model of the heart is presented. The presented research investigates an off-line optimization algorithm based on different electrode positioning and timing delays

Doctor of Philosophy

dissertationComputational simulation has become an indispensable tool in the study of both basic mechanisms and pathophysiology of all forms of cardiac electrical activity. Because the heart is comprised of approximately 4 billion electrically active cells, it is not possible to geometrically model or computationally simulate each individual cell. As a result computational models of the heart are, of necessity, abstractions that approximate electrical behavior at the cell, tissue, and whole body level. The goal of this PhD dissertation was to evaluate several aspects of these abstractions by exploring a set of modeling approaches in the field of cardiac electrophysiology and to develop means to evaluate both the amplitude of these errors from a purely technical perspective as well as the impacts of those errors in terms of physiological parameters. The first project used subject specific models and experiments with acute myocardial ischemia to show that one common simplification used to model myocardial ischemia-the simplest form of the border zone between healthy and ischemic tissue-was not supported by the experimental results. We propose a alternative approximation of the border zone that better simulates the experimental results. The second study examined the impact of simplifications in geometric models on simulations of cardiac electrophysiology. Such models consist of a connected mesh of polygonal elements and must often capture complex external and internal boundaries. A conforming mesh contains elements that follow closely the shapes of boundaries; nonconforming meshes fit the boundaries only approximately and are easier to construct but their impact on simulation accuracy has, to our knowledge, remained unknown. We evaluated the impact of this simplification on a set of three different forms of bioelectric field simulations. The third project evaluated the impact of an additional geometric modeling error; positional uncertainty of the heart in simulations of the ECG. We applied a relatively novel and highly efficient statistical approach, the generalized Polynomial Chaos-Stochastic Collocation method (gPC-SC), to a boundary element formulation of the electrocardiographic forward problem to carry out the necessary comprehensive sensitivity analysis. We found variations large enough to mask or to mimic signs of ischemia in the ECG