In this paper, we consider the inverse problem of space dependent multiple ionic parameters identification in cardiac electrophysiology modelling from a set of observations. We use the monodomain system known as a state-of-the-art model in cardiac electrophysiology and we consider a general Hodgkin-Huxley formalism to describe the ionic exchanges at the microscopic level. This formalism covers many physiological transmembrane potential models including those in cardiac electrophysiology. Our main result is the proof of the uniqueness and a Lipschitz stability estimate of ion channels conductance parameters based on some observations on an arbitrary subdomain. The key idea is a Carleman estimate for a parabolic operator with multiple coefficients and an ordinary differential equation system.

In this paper, we consider an inverse problem of determining multiple ionic parameters of a 2 × 2 strongly coupled parabolic–elliptic reaction–diffusion system arising in cardiac electrophysiology modeling. We use the bidomain model coupled to an ordinary differential equation (ODE) system and we consider a general formalism of physiologically detailed cellular membrane models to describe the ionic exchanges at the microscopic level. Our main result is the uniqueness and a Lipschitz stability estimate of the ion channels conductance parameters of the model using subboundary observations over an interval of time. The key ingredients are a global Carleman-type estimate with a suitable observations acting on a part of the boundary.

In this work, we present an optimal control formulation for the bidomain model in order to estimate maximal conductances parameters in the physiological ionic model. We consider a general Hodgkin-Huxley formalism to describe the ionic exchanges at the microcopic level. We consider the parameters as control variables to minimize the mismatch between the measured and the computed potentials under the constraint of the bidomain system. The solution of the optimization problem is based on a gradient descent method, where the gradient is obtained by solving an adjoint problem. We show through some numerical examples the capability of this approach to estimate the values of sodium, calcium and potassium ion channels conductances in the Luo Rudy phase I model.