arxiv: 2604.03682 · v1 · submitted 2026-04-04 · 🧮 math.AP

Recognition: 2 theorem links

· Lean Theorem

A boundary integral approach to the eigenvalue problem for the anisotropic bidomain operator with perfect contact conditions

Raul Felipe-Sosa , Yofre H. Garc\'ia-G\'omez

Authors on Pith no claims yet

Pith reviewed 2026-05-13 17:30 UTC · model grok-4.3

classification 🧮 math.AP

keywords bidomain operatoranisotropic conductivityboundary integral equationseigenvalue problemlayer potentialsHelmholtz operatorBessel functionstransmission conditions

0 comments

The pith

Solutions to the anisotropic bidomain eigenvalue problem are reduced to a system of Fredholm boundary integral equations via single- and double-layer potentials, with explicit Bessel-function kernels for the anisotropic Helmholtz operator.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper reformulates the eigenvalue problem for the anisotropic bidomain operator across three subdomains (left ventricle, septum, right ventricle) by representing the solution with single- and double-layer potentials. This converts the original transmission problem with discontinuous coefficients into a boundary integral system on the interfaces. Explicit fundamental solutions and kernels are derived in terms of Bessel functions. A discretization scheme for the resulting integral equations is proposed to approximate the eigenvalues numerically.

Core claim

By expressing the solution in terms of single- and double-layer potentials, the original boundary value problem is reduced to a system of Fredholm-type boundary integral equations. Explicit expressions are derived for the fundamental solution of the associated anisotropic Helmholtz operator and for the corresponding kernels, given in terms of Bessel functions. A numerical scheme is then proposed based on discretization of the resulting integral system.

What carries the argument

Single- and double-layer potentials for the anisotropic Helmholtz operator, which reduce the transmission problem with perfect contact conditions to boundary integral equations via standard jump relations.

Load-bearing premise

The fiber orientations are piecewise constant within each subdomain and the interfaces satisfy perfect contact conditions, allowing standard jump relations to handle the transmission conditions.

What would settle it

Compute the first few eigenvalues via the discretized integral system and check whether they match independent finite-element or finite-difference solutions of the original PDE on the same three-domain geometry.

Figures

Figures reproduced from arXiv: 2604.03682 by Raul Felipe-Sosa, Yofre H. Garc\'ia-G\'omez.

**Figure 1.** Figure 1: Geometric representation of the Ω region As we will see later, it is convenient to introduce the following notation: • We define ∂ΩI = ΓI ∪ ΓIS, where ΓI = (x1, b) : x1 ∈ [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗

**Figure 2.** Figure 2: Geometric representation of the longitudinal fiber direction [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗

**Figure 3.** Figure 3: Numerical representation of the solutions [PITH_FULL_IMAGE:figures/full_fig_p025_3.png] view at source ↗

read the original abstract

In this work, we study the eigenvalue problem associated with the bidomain operator in an anisotropic heterogeneous domain composed of three subregions representing the left ventricle, the septum, and the right ventricle. The anisotropic conductivity, together with the different orientations of the fiber directions in each subdomain, leads to an elliptic boundary value problem with discontinuous coefficients and transmission conditions across the interfaces. Our main contribution consists in reformulating this problem using potential theory. By expressing the solution in terms of single- and double-layer potentials, we reduce the original boundary value problem to a system of Fredholm-type boundary integral equations. We derive explicit expressions for the fundamental solution of the associated anisotropic Helmholtz operator, as well as for the corresponding kernels, which are given in terms of Bessel functions. Finally, we propose a numerical scheme for approximating the eigenvalues of the bidomain operator based on the discretization of the resulting integral system. This approach provides an efficient framework for the analysis of anisotropic boundary value problems with interface conditions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives explicit Bessel kernels for the anisotropic Helmholtz operator in a three-subdomain bidomain setup and reduces the eigenvalue problem to a Fredholm boundary system, which is a clean but narrow advance.

read the letter

The main point is that the authors derive explicit fundamental solutions and layer-potential kernels for the piecewise-anisotropic Helmholtz operator, then turn the transmission eigenvalue problem into a system of boundary integral equations on the interfaces between the left ventricle, septum, and right ventricle. Because the conductivity tensors are constant inside each piece, an affine change of variables reduces each local problem to the isotropic case, and the kernels come out in terms of Bessel functions. That explicit step is the part that is not completely standard and could be useful for people who want to avoid volume meshes in cardiac simulations. The perfect-contact conditions are handled with the usual jump relations, so the reduction itself follows classical potential theory without new inventions. The stress-test note is right that nothing in the construction looks internally inconsistent. The soft spot is that the abstract only sketches the discretization scheme and supplies no numerical examples, convergence checks, or comparisons with existing finite-element codes. Without those, it is hard to tell how practical the method becomes on realistic heart geometries. The theoretical claim that the resulting operator is Fredholm of index zero is asserted but not shown in the provided summary, so a referee would want to see the details of that argument. This paper is for a small group of researchers working on boundary-integral methods for anisotropic elliptic problems in physiology. It does not reorganize general PDE theory, but it gives a concrete computational route for this specific setting. I would send it to peer review; the core derivation is worth a careful check even if the numerics still need work.

Referee Report

1 major / 2 minor

Summary. The paper studies the eigenvalue problem for the anisotropic bidomain operator on a heterogeneous domain partitioned into left ventricle, septum, and right ventricle subregions with piecewise-constant fiber orientations. It reformulates the transmission problem via single- and double-layer potentials for the associated anisotropic Helmholtz operator, reduces it to a system of Fredholm boundary integral equations whose kernels are expressed explicitly in terms of Bessel functions after local affine transformations, and outlines a numerical discretization scheme for the resulting integral operator.

Significance. If the derivations hold, the explicit kernel expressions and potential-theoretic reduction supply a concrete, parameter-free route to eigenvalue computation for piecewise-anisotropic elliptic problems with perfect-contact interfaces. This is a useful technical contribution for cardiac bidomain modeling, where such eigenvalue problems arise in stability analysis; the approach avoids volume discretization and exploits the standard jump relations for layer potentials once the fundamental solution is obtained.

major comments (1)

§3 (reduction step): the manuscript asserts that the transmission conditions are satisfied by the usual jump relations for single- and double-layer potentials, but does not explicitly verify that the resulting boundary integral operator is Fredholm of index zero on the chosen trace spaces; a short reference to the standard theory for elliptic transmission problems or a direct estimate of the principal symbol would strengthen the claim.

minor comments (2)

Notation: the conductivity tensors are denoted differently in the differential formulation and in the integral kernels; a single consistent symbol table would improve readability.
Numerical section: the discretization scheme is sketched but no convergence rates or quadrature error analysis for the Bessel kernels is supplied; adding a brief remark on the handling of the logarithmic singularity would be helpful.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive comment on our manuscript. We address the point raised below and will incorporate the suggested clarification in the revised version.

read point-by-point responses

Referee: §3 (reduction step): the manuscript asserts that the transmission conditions are satisfied by the usual jump relations for single- and double-layer potentials, but does not explicitly verify that the resulting boundary integral operator is Fredholm of index zero on the chosen trace spaces; a short reference to the standard theory for elliptic transmission problems or a direct estimate of the principal symbol would strengthen the claim.

Authors: We agree that an explicit reference would strengthen the presentation. In the revised manuscript we will add a short paragraph in §3 that recalls the standard theory for elliptic transmission problems with piecewise-constant coefficients. This theory guarantees that the boundary integral operator obtained from the single- and double-layer potentials is Fredholm of index zero on the natural trace spaces (H^{1/2} and H^{-1/2} on the interfaces). We will cite the relevant results from the literature on boundary integral methods for transmission problems (e.g., the framework developed in works on anisotropic elliptic operators with jump conditions) and note that the principal symbol remains elliptic after the local affine transformations used to express the kernels in terms of Bessel functions. revision: yes

Circularity Check

0 steps flagged

No significant circularity in the derivation

full rationale

The paper starts from the piecewise-anisotropic elliptic eigenvalue problem with perfect-contact transmission conditions and reduces it to an equivalent system of Fredholm boundary integral equations by representing the solution via single- and double-layer potentials. The fundamental solution of the anisotropic Helmholtz operator is obtained from the isotropic kernel by an affine change of variables (preserving the required singularity) and expressed using Bessel functions; the kernels follow directly from this construction. Jump relations for the layer potentials then enforce the transmission conditions without additional assumptions or fitted quantities. No step equates a derived quantity to its own input by definition, no parameters are calibrated to data, and no load-bearing self-citation chain is required. The derivation is therefore a direct, self-contained mathematical equivalence.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard results from potential theory for elliptic operators with discontinuous coefficients. No free parameters, ad-hoc constants, or new postulated entities are introduced in the abstract.

axioms (2)

standard math The anisotropic Helmholtz operator admits a fundamental solution expressible in closed form via Bessel functions when the conductivity tensor is constant in each subdomain.
Invoked to obtain the explicit kernels used in the layer potentials.
domain assumption Jump relations for single- and double-layer potentials hold across interfaces with perfect contact transmission conditions.
Required to convert the volume transmission problem into boundary integral equations.

pith-pipeline@v0.9.0 · 5478 in / 1449 out tokens · 50633 ms · 2026-05-13T17:30:53.846689+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

By expressing the solution in terms of single- and double-layer potentials, we reduce the original boundary value problem to a system of Fredholm-type boundary integral equations. We derive explicit expressions for the fundamental solution of the associated anisotropic Helmholtz operator, as well as for the corresponding kernels, which are given in terms of Bessel functions.
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

the fundamental solution to the Helmholtz problem (12) is given by u(∗,I)(x,y;λ)=1/(4√(σl σt)) N0(λ r(I)(x,y))

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

[1]

A bidomain model for describing ischemic myocardial D-C potentials

L. Tung. “A bidomain model for describing ischemic myocardial D-C potentials.” In: Ph.D Ph.D. dissertation, Massachusetts Inst. Technol,(1978)

work page 1978
[2]

Sundnes et al.Computing the Electrical Activity in the Heart

J. Sundnes et al.Computing the Electrical Activity in the Heart. Monographs in Computational Science and Engineering. Springer, Berlin Heidelberg, 2006. 25

work page 2006
[3]

Evolution Equations, Semigroups and Functional Anal- ysis: In Memory of Brunello Terreni

P . Colli-Franzone and G. Savaré. “Evolution Equations, Semigroups and Functional Anal- ysis: In Memory of Brunello Terreni”. In: Evolution Equations, Semigroups and Func- tional Analysis, 2002. Chap. Degenerate Evolution Systems Modeling the Cardiac Electric Field at Micro- and Macroscopic Level

work page 2002
[4]

Reaction-diffusion systems for the macroscopic bidomain model of the cardiac electric field

M. Veneroni. “Reaction-diffusion systems for the macroscopic bidomain model of the cardiac electric field”. In:Nonlinear Anal., Real World Appl.(2009)

work page 2009
[5]

Existence and uniqueness of the solution for the bidomain model used in cardiac electrophysiology

Y. Bourgault, Y. Coudière, and C. Pierre. “Existence and uniqueness of the solution for the bidomain model used in cardiac electrophysiology.” In:Nonlinear Anal., Real World Appl.(2009)

work page 2009
[6]

Existence of a T-Periodic solution for the monodomain model corre- sponding to an isolated ventricle due to ionic-diffusive relations

A. Fraguela et al. “Existence of a T-Periodic solution for the monodomain model corre- sponding to an isolated ventricle due to ionic-diffusive relations”. In:Acta Applicandae Mathematicae(2022)

work page 2022
[7]

On the strong convergence of the Faedo-Galerkin approximations to a strong t-periodic solution of the torso-coupled bidomain model

R. Felipe-Sosa, A. Fraguela-Collar, and Y. H. García-Gómez. “On the strong convergence of the Faedo-Galerkin approximations to a strong t-periodic solution of the torso-coupled bidomain model”. In:Math. Model. Nat. Phenom.(2023)

work page 2023
[8]

Mathematical modeling of electriocardiograms: A numerical study

M. Boulakia et al. “Mathematical modeling of electriocardiograms: A numerical study”. In:Annals of Biomedical Engineering(2010)

work page 2010
[9]

Simulating patterns of excitation, repo- larization and action potential duration with cardiac Bidomain and Monodomain mod- els

P . Colli Franzone, L.F. Pavarino, and B. Taccardi. “Simulating patterns of excitation, repo- larization and action potential duration with cardiac Bidomain and Monodomain mod- els”. In:Mathematical Biosciences(2006)

work page 2006
[10]

Isogeometric approximation of cardiac electrophysiology models on surfaces: An accuracy study with application to the human left atrium

Alessandro S. Patelli et al. “Isogeometric approximation of cardiac electrophysiology models on surfaces: An accuracy study with application to the human left atrium”. In: Computer Methods in Applied Mechanics and Engineering(2017)

work page 2017
[11]

Multipatch Isogeometric Analysis for electrophysiology: Simula- tion in a humanheart

Michele Bucelli et al. “Multipatch Isogeometric Analysis for electrophysiology: Simula- tion in a humanheart”. In:Computer Methods in Applied Mechanics and Engineering(2021)

work page 2021
[12]

Isogeometric overlapping Schwarz preconditioners for the Bidomain reaction–diffusion system

Lara Antonella Charawi. “Isogeometric overlapping Schwarz preconditioners for the Bidomain reaction–diffusion system”. In:Computer Methods in Applied Mechanics and En- gineering(2017)

work page 2017
[13]

An experimental contribution to the concept of “Jumping wave

G. A Medrano, A. de Micheli, and R. Chávez Domínguez A. Aranda P . Iturralde. “An experimental contribution to the concept of “Jumping wave” phenomenon in the inter- ventricular septum”. In:Archives of Cardiology of Mexico(2002)

work page 2002
[14]

The activation of the interventricular septum in the dog’s heart under normal conditions and in bundle-branch block

G. Medrano et al. “The activation of the interventricular septum in the dog’s heart under normal conditions and in bundle-branch block”. In:Ann N Y Acad Sci(1957)

work page 1957
[15]

Springer, 2002

A.M.Linkov.Boundary Integral Equations in Elasticity Theory. Springer, 2002

work page 2002
[16]

Anisotropic Cardiac Conduction

Kotadia I et al. “Anisotropic Cardiac Conduction”. In:Arrhythm Electrophysiol Rev.(2020)

work page 2020
[17]

Theory and Practice

Norbert Ortner and Peter Wagner.Fundamental Solutions of Linear Partial Differential Op- erators. Theory and Practice. London: Springer, 2015

work page 2015
[18]

Bell.Special functions for and engineers

William W. Bell.Special functions for and engineers. London: Van Nostrand, 1968

work page 1968
[19]

Korenev.Bessel functions and their applications

Boris G. Korenev.Bessel functions and their applications. Boca Raton: Chapman & Hall/CRC, 2002

work page 2002
[20]

Philadelphia: SIAM, 2013

David Colton and Rainer Kress.Integral equation methods in scattering theory. Philadelphia: SIAM, 2013

work page 2013
[21]

C. A. Brebbia, J. C. F. Telles, and L. C. Wrobel.Boundary Element Techniques. Berlin: Springer, 1984. 26

work page 1984