skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Integral equation method for the 1D steady-state Poisson-Nernst-Planck equations
Abstract An integral equation method is presented for the 1D steady-state Poisson-Nernst-Planck equations modeling ion transport through membrane channels. The differential equations are recast as integral equations using Green’s 3rd identity yielding a fixed-point problem for the electric potential gradient and ion concentrations. The integrals are discretized by a combination of midpoint and trapezoid rules, and the resulting algebraic equations are solved by Gummel iteration. Numerical tests for electroneutral and non-electroneutral systems demonstrate the method’s 2nd order accuracy and ability to resolve sharp boundary layers. The method is applied to a 1D model of the K$$^+$$ +  ion channel with a fixed charge density that ensures cation selectivity. In these tests, the proposed integral equation method yields potential and concentration profiles in good agreement with published results.  more » « less
Award ID(s):
1819193 2110869 2110767 1819094
PAR ID:
10526822
Author(s) / Creator(s):
; ;
Publisher / Repository:
Springer
Date Published:
Journal Name:
Journal of Computational Electronics
Volume:
22
Issue:
5
ISSN:
1569-8025
Page Range / eLocation ID:
1396 to 1408
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We introduce a family of Finsler metrics, called the$$L^p$$ L p -Fisher–Rao metrics$$F_p$$ F p , for$$p\in (1,\infty )$$ p ( 1 , ) , which generalizes the classical Fisher–Rao metric$$F_2$$ F 2 , both on the space of densities$${\text {Dens}}_+(M)$$ Dens + ( M ) and probability densities$${\text {Prob}}(M)$$ Prob ( M ) . We then study their relations to the Amari–C̆encov$$\alpha $$ α -connections$$\nabla ^{(\alpha )}$$ ( α ) from information geometry: on$${\text {Dens}}_+(M)$$ Dens + ( M ) , the geodesic equations of$$F_p$$ F p and$$\nabla ^{(\alpha )}$$ ( α ) coincide, for$$p = 2/(1-\alpha )$$ p = 2 / ( 1 - α ) . Both are pullbacks of canonical constructions on$$L^p(M)$$ L p ( M ) , in which geodesics are simply straight lines. In particular, this gives a new variational interpretation of$$\alpha $$ α -geodesics as being energy minimizing curves. On$${\text {Prob}}(M)$$ Prob ( M ) , the$$F_p$$ F p and$$\nabla ^{(\alpha )}$$ ( α ) geodesics can still be thought as pullbacks of natural operations on the unit sphere in$$L^p(M)$$ L p ( M ) , but in this case they no longer coincide unless$$p=2$$ p = 2 . Using this transformation, we solve the geodesic equation of the$$\alpha $$ α -connection by showing that the geodesic are pullbacks of projections of straight lines onto the unit sphere, and they always cease to exists after finite time when they leave the positive part of the sphere. This unveils the geometric structure of solutions to the generalized Proudman–Johnson equations, and generalizes them to higher dimensions. In addition, we calculate the associate tensors of$$F_p$$ F p , and study their relation to$$\nabla ^{(\alpha )}$$ ( α )
    more » « less
  2. Abstract The Markoff graphs modulopwere proven by Chen (Ann Math 199(1), 2024) to be connected for all but finitely many primes, and Baragar (The Markoff equation and equations of Hurwitz. Brown University, 1991) conjectured that they are connected for all primes, equivalently that every solution to the Markoff equation moduloplifts to a solution over$$\mathbb {Z}$$ Z . In this paper, we provide an algorithmic realization of the process introduced by Bourgain et al. [arXiv:1607.01530] to test whether the Markoff graph modulopis connected for arbitrary primes. Our algorithm runs in$$o(p^{1 + \epsilon })$$ o ( p 1 + ϵ ) time for every$$\epsilon > 0$$ ϵ > 0 . We demonstrate this algorithm by confirming that the Markoff graph modulopis connected for all primes less than one million. 
    more » « less
  3. Abstract We consider the Cauchy problem for the logarithmically singular surface quasi-geostrophic (SQG) equation, introduced by Ohkitani,$$\begin{aligned} \begin{aligned} \partial _t \theta - \nabla ^\perp \log (10+(-\Delta )^{\frac{1}{2}})\theta \cdot \nabla \theta = 0, \end{aligned} \end{aligned}$$ t θ - log ( 10 + ( - Δ ) 1 2 ) θ · θ = 0 , and establish local existence and uniqueness of smooth solutions in the scale of Sobolev spaces with exponent decreasing with time. Such a decrease of the Sobolev exponent is necessary, as we have shown in the companion paper (Chae et al. in Illposedness via degenerate dispersion for generalized surface quasi-geostrophic equations with singular velocities,arXiv:2308.02120) that the problem is strongly ill-posed in any fixed Sobolev spaces. The time dependence of the Sobolev exponent can be removed when there is a dissipation term strictly stronger than log. These results improve wellposedness statements by Chae et al. (Comm Pure Appl Math 65(8):1037–1066, 2012). This well-posedness result can be applied to describe the long-time dynamics of the$$\delta $$ δ -SQG equations, defined by$$\begin{aligned} \begin{aligned} \partial _t \theta + \nabla ^\perp (10+(-\Delta )^{\frac{1}{2}})^{-\delta }\theta \cdot \nabla \theta = 0, \end{aligned} \end{aligned}$$ t θ + ( 10 + ( - Δ ) 1 2 ) - δ θ · θ = 0 , for all sufficiently small$$\delta >0$$ δ > 0 depending on the size of the initial data. For the same range of$$\delta $$ δ , we establish global well-posedness of smooth solutions to the dissipative SQG equations. 
    more » « less
  4. Abstract We construct an example of a group$$G = \mathbb {Z}^2 \times G_0$$ G = Z 2 × G 0 for a finite abelian group $$G_0$$ G 0 , a subsetEof $$G_0$$ G 0 , and two finite subsets$$F_1,F_2$$ F 1 , F 2 of G, such that it is undecidable in ZFC whether$$\mathbb {Z}^2\times E$$ Z 2 × E can be tiled by translations of$$F_1,F_2$$ F 1 , F 2 . In particular, this implies that this tiling problem isaperiodic, in the sense that (in the standard universe of ZFC) there exist translational tilings ofEby the tiles$$F_1,F_2$$ F 1 , F 2 , but no periodic tilings. Previously, such aperiodic or undecidable translational tilings were only constructed for sets of eleven or more tiles (mostly in $$\mathbb {Z}^2$$ Z 2 ). A similar construction also applies for$$G=\mathbb {Z}^d$$ G = Z d for sufficiently large d. If one allows the group$$G_0$$ G 0 to be non-abelian, a variant of the construction produces an undecidable translational tiling with only one tile F. The argument proceeds by first observing that a single tiling equation is able to encode an arbitrary system of tiling equations, which in turn can encode an arbitrary system of certain functional equations once one has two or more tiles. In particular, one can use two tiles to encode tiling problems for an arbitrary number of tiles. 
    more » « less
  5. Abstract The skew mean curvature flow is an evolution equation forddimensional manifolds embedded in$${\mathbb {R}}^{d+2}$$ R d + 2 (or more generally, in a Riemannian manifold). It can be viewed as a Schrödinger analogue of the mean curvature flow, or alternatively as a quasilinear version of the Schrödinger Map equation. In an earlier paper, the authors introduced a harmonic/Coulomb gauge formulation of the problem, and used it to prove small data local well-posedness in dimensions$$d \geqq 4$$ d 4 . In this article, we prove small data local well-posedness in low-regularity Sobolev spaces for the skew mean curvature flow in dimension$$d\geqq 2$$ d 2 . This is achieved by introducing a new, heat gauge formulation of the equations, which turns out to be more robust in low dimensions. 
    more » « less