Restitution Curve

See code in GitLab.
Author: Cristian Barrios Espinosa cristian.espinosa@kit.edu

Finding Conduction Velocity Restitution Parameters

This is a tutorial about obtaining conduction velocity (CV) restitution curve paramters to be used in Diffusion Reaction Eikonal Alternans Model (DREAM). This tutorial demonstrates how to run monodomain simulations to obtain a conduction velocity (CV) restitution curve for a selected ionic model, a specific parameter set, and a predefined cellular pre-pacing interval. The resulting conduction velocity (CV) restitution data will be used to derive parameters required as input for DREAM simulations.

To run the experiments described in this tutorial, navigate to the following directory:

cd ${TUTORIALS}/03_Eikonal/03C_CV_Restitution_Fitting

Simulation Setup

To find the best fit for the CV restitution curve parameters, a monodomain simulation is first run using the desired ionic model and its parameters in a 2D mesh. The ionic model, the ionic model parameters and the desired prepace cellular interval should be provided

python3 run.py --help
...
Tutorial specific options:
--ionic_model IONIC_MODEL   Set ionic model (default: Courtemanche)
--model_par MODEL_PAR       Set ionic model parameters (default: )
--prepace_bcl PREPACE_BCL   Basic cycle length for cellular prepacing (ms) (default: 1000)

Use the runtime option --visualize to plot the CV restitution curve.

Experiment

Find the cv restitution parameters for the Courtemanche model with parameters for atrial fibrillation as described in [Loewe2016].

This experiment finds the CV restitution parameters for the Courtemanche model with atrial fibrillation parameters and a cellular prepacing at an interval of 250ms [Courtemanche1998] [Loewe2016].

python3 run.py --model_par 'GCaL*0.45,GK1*2,GNa*1,factorGKur*0.5,Gto*0.35,GKs*2,maxIpCa*1.5,maxINaCa*1.6' --prepace_bcl 250

The output will provide the fitted parameters for the CV restitution formula used in DREAM that corresponds to the ionic model and set of parameters used:

Fitted parameters:
    rho = 2108654.8268
    kappa = 105.1744
    theta = 62.8912

Recent questions tagged experiments, examples, tissue, eikonal, dream, restitution