[Paper Review] Evidence for an unfolded border-collision bifurcation in paced cardiac tissue

Carolyn Berger, Zhao, Xiaopeng+4 more|ArXiv.org|Jul 13, 2006

Nonlinear Dynamics and Pattern Formation4 references3 citations

TL;DR

This paper proposes an 'unfolded border-collision bifurcation' model to explain hybrid dynamics in paced cardiac tissue, where alternans transitions exhibit smooth-like behavior near the bifurcation point and border-collision-like behavior farther away. Using alternate pacing experiments and a modified skew-tent map with a small smoothing parameter ε, the study demonstrates that the system’s sensitivity (Γ) decreases with δ (perturbation size) near the bifurcation, consistent with smooth dynamics, but increases with δ farther from it, consistent with border-collision behavior—offering a unified explanation for conflicting prior models.

ABSTRACT

We investigate, both experimentally and theoretically, the bifurcation to alternans in heart tissue. Previously, this phenomenon has been modeled either as a smooth or as border-collision period-doubling bifurcation. Using a new experimental technique, we find a hybrid behavior: very close to the bifurcation point the dynamics are smooth-like, whereas further away they are border-collision-like. This behavior is captured by a new type of model, called an unfolded border-collision bifurcation.

Motivation & Objective

• To resolve the long-standing debate on whether alternans in paced cardiac tissue arises from a smooth or border-collision period-doubling bifurcation.
• To develop a new theoretical framework that reconciles experimental observations showing both smooth-like and border-collision-like behavior in the same system.
• To investigate the sensitivity of action potential duration (APD) to periodic perturbations using alternate pacing, a method that reveals bifurcation type through gain Γ.
• To test whether the clinical application of alternate pacing for predicting arrhythmic risk is viable, based on bifurcation dynamics.

Proposed method

• Employed alternate pacing with periodic perturbations of the basic cycle length (B₀) by ±δ to measure the system’s response via the gain Γ = (APD_long - APD_short)/(2δ).
• Applied a modified skew-tent map: x_{n+1} = μ - αx_n - β√(x_n² + ε²), which unfolds the singular border-collision map (3) into a smooth, differentiable system for ε > 0.
• Used numerical simulations of the unfolded map (4) to compute Γ as a function of bifurcation parameter μ and perturbation size δ, comparing trends to experimental data.
• Defined the bifurcation point at μ_bif = -0.0091 for the unfolded map (4), with ε = 0.01, and analyzed Γ behavior across varying distances from μ_bif.
• Mapped the transition in Γ behavior: decreasing with δ near μ < 0 (smooth-like), increasing with δ at μ > 0 (border-collision-like), indicating a scale-dependent bifurcation type.
• Validated experimental Γ trends against theoretical predictions, showing qualitative agreement between simulation and data across multiple pacing intervals.

Experimental results

Research questions

• RQ1 Can experimental data distinguish between a smooth and a border-collision period-doubling bifurcation in paced cardiac tissue?
• RQ2 Why do previous models of alternans bifurcation yield conflicting interpretations despite similar bifurcation diagrams?
• RQ3 Does the system’s sensitivity (Γ) to perturbations change qualitatively depending on proximity to the bifurcation point?
• RQ4 Can a single model explain both smooth-like and border-collision-like behavior observed in the same experimental system?
• RQ5 Is alternate pacing a viable clinical tool for predicting alternans and arrhythmic risk, given the underlying bifurcation dynamics?

Key findings

• Experimental Γ values decrease with increasing δ when close to the bifurcation point (μ < 0), indicating smooth-like dynamics, consistent with a smooth period-doubling bifurcation.
• Experimental Γ values increase with increasing δ when farther from the bifurcation point (μ > 0), indicating border-collision-like dynamics, inconsistent with smooth models.
• Theoretical simulations of the unfolded map (4) with ε = 0.01 reproduce the hybrid Γ behavior: decreasing Γ with δ near μ = -0.0073 and increasing Γ with δ at μ = 0.01.
• The bifurcation point for the unfolded map (4) occurs at μ_bif = -0.0091, which is shifted from the original border-collision point at μ = 0, due to the smoothing parameter ε.
• The region of smooth-like behavior (Γ decreasing with δ) is narrow—on the order of ε in width—making it difficult to detect experimentally, which may explain prior misclassification.
• The study challenges the clinical use of alternate pacing to predict arrhythmic risk, as Γ remains small until pacing rates are very close to the bifurcation, limiting its predictive utility.