Non-Invasive Reconstruction of Cardiac Activation Dynamics Using Physics-Informed Neural Networks

TL;DR

A physics-informed neural network framework for recovering cardiac activation patterns, active tension propagation, deformation fields, and hydrostatic pressure from measurable deformation data in simplified left ventricular geometries is presented.

Abstract

Cardiac arrhythmogenesis is governed by complex electromechanical interactions that are not directly observable in vivo, motivating the development of non-invasive computational approaches for reconstructing three-dimensional activation dynamics. We present a physics-informed neural network framework for recovering cardiac activation patterns, active tension propagation, deformation fields, and hydrostatic pressure from measurable deformation data in simplified left ventricular geometries. Our approach integrates nonlinear anisotropic constitutive modeling, heterogeneous fiber orientation, weak formulations of the governing mechanics, and finite-element-based loss functions to embed physical constraints directly into training. We demonstrate that the proposed framework accurately reconstructs spatiotemporal activation dynamics under varying levels of measurement noise and reduced spatial resolution, while preserving global propagation patterns and activation timing. By coupling mechanistic modeling with data-driven inference, this method establishes a pathway toward patient-specific, non-invasive reconstruction of cardiac activation, with potential applications in digital phenotyping and computational support for arrhythmia assessment.

Figures (6)

• Figure 1: Schematic overview of the PINN framework. The input consists of a 3D volumetric mesh together with available deformation data defined on the vertices. Laplacian eigenfunctions $f_{n}(\mathbf{r})$ are then calculated as a preprocessing step and serve as the input coordinates for the PINN. The spatial-temporal fields of the deformation vector $\mathbf{U}(\mathbf{r},t)$, active tension $T_a(\mathbf{r},t)$ and hydrostatic pressure $p(\mathbf{r},t)$ are represented by three fully connected NNs. The optimization is done by minimizing two terms that make up the total loss function: the data loss which is evaluated on $N_d$ data points and the physics loss enforced on $N_c$ collocation points.
• Figure 2: PINN results for the estimated deformation field $U$ and hydrostatic pressure $p$, optimized on all available data points. The fields are visualized by slicing the 3D volumetric mesh with the $x=0$ mm plane (longitudinal, left) and the $z=-5$ mm plane (radial, right). Both planes slice through the source point, see Fig. \ref{['fig:noise']} for the 3D orientation. Panels (a) and (b) show the y-component of the deformation, while panels (c) and (d) consist of the hydrostatic pressure $p$, used in the forward simulation to satisfy the incompressibility constraint.
• Figure 3: PINN results for the estimated active tension $T_a$, optimized on all available data points. The fields are visualized by slicing the 3D volumetric mesh with the $x=0$ mm plane (longitudinal, panel (a)) and the $z=-5$ mm plane (radial, panel (b)). Both planes slice through the source point, see Fig. \ref{['fig:noise']} for the 3D orientation. In addition to the $T_a$-field, the last column shows the corresponding LAT map for the activation wave.
• Figure 4: PINN results for the estimated active tension $T_a$, optimized on all available data points. Panel (a) shows polar maps of the endo- and epicardial surfaces, calculated from the longitudinal and circumferential coordinates. The continuous field is visualized by linearly interpolating the vertex values in the new space. Panel (b) presents a cross-section of the volume by slicing the 3D volumetric mesh with the $x=-5$ mm plane which, unlike other figures, does not go through the source point, revealing the transmural differences in wave propagation.
• Figure 5: PINN results for the estimated active tension $T_a$ and corresponding LAT map, optimized on all available data points with Gaussian noise levels of $5\;\%$ and $10\;\%$ of the maximum values. The datasets are illustrated in the first column by the y-component of the deformation at one timestep. The plane in the first row indicates the 3D orientation of the cross-sections, panel (a) for the longitudinal view (($x=0$ mm) and panel (b) for the radial ($z=-5$ mm). Both planes go through the source point.