pith. sign in

arxiv: 2606.28223 · v1 · pith:4KCMP6FUnew · submitted 2026-06-26 · 🧮 math.AP

Analysis and Numerics of a Stationary Drift-Diffusion Model for Electrical Discharge in MEMS

Runan He , Wenjia Xie This is my paper

Pith reviewed 2026-06-29 03:06 UTC · model grok-4.3

classification 🧮 math.AP
keywords drift-diffusion modelelectrical dischargeMEMSweak solutionsimpact ionizationfinite element methodPoisson equationTownsend term
0
0 comments X

The pith

The stationary drift-diffusion model for electrical discharge in MEMS admits weak solutions with uniformly bounded carrier densities.

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

This paper proves existence of weak solutions for a coupled Poisson-continuity system that includes an exponential Townsend-type impact ionization term. The proof proceeds by regularization and approximation of the nonlinearities, followed by application of monotone operator theory and Stampacchia truncation to obtain the required bounds. The result supplies a rigorous foundation for the model under physically relevant assumptions. A finite-element solver is constructed and applied to two-dimensional domains and a three-dimensional axisymmetric geometry to compute the potential and carrier profiles.

Core claim

The authors prove the existence of weak solutions to the stationary drift-diffusion model under physically relevant assumptions and establish uniform bounds on the carrier densities. The proof relies on a regularization-approximation scheme with truncated nonlinearities, monotone operator theory (Browder-Minty), iterative energy estimates, and Stampacchia-type truncation arguments. They also develop a robust finite element solver to simulate the carrier density and electrostatic potential profiles for several geometries, including two-dimensional domains and a three-dimensional axisymmetric geometry.

What carries the argument

Regularization-approximation scheme combined with monotone operator theory and Stampacchia truncation arguments applied to the exponential Townsend-type impact ionization source term.

If this is right

  • Weak solutions exist for the stationary model equations.
  • Carrier densities remain uniformly bounded under the stated assumptions.
  • The finite-element discretization produces computable profiles for the tested two- and three-dimensional geometries.
  • The analysis supplies a foundation for further numerical exploration of discharge behavior in MEMS devices.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The uniform bounds may furnish a priori estimates usable in a time-dependent version of the same system.
  • The numerical solver could be extended to treat electrode-shape optimization problems in MEMS design.
  • Techniques developed here may transfer to other models whose source terms share the same exponential dependence on field strength.

Load-bearing premise

The existence result depends on the specific exponential form of the Townsend-type impact ionization source term together with the physically relevant assumptions that enable the regularization scheme and truncation arguments to close.

What would settle it

A concrete parameter set satisfying the physically relevant assumptions for which either the approximating sequence fails to converge or the carrier densities become unbounded would falsify the existence and bound claims.

Figures

Figures reproduced from arXiv: 2606.28223 by Runan He, Wenjia Xie.

Figure 1
Figure 1. Figure 1: Sketch of a MEMS capacitor undergoing sparking. The domain Ω is depicted by the shaded region. Its boundary decomposes into a Dirichlet portion, formed by electrodes A and B, where the electric potential ϕ together with the concentrations p and n are prescribed, and a Neumann portion C = ∂ΩN , where the artificial boundary is placed far enough away that homogeneous Neumann conditions can reasonably be impo… view at source ↗
Figure 2
Figure 2. Figure 2: Contour plots of the electric potential ϕ, positive ion density p and electron density n on the L-shaped domain Annular Domain (Corona Geometry) The annulus r ∈ [1, 5], with inner radius rin = 1 and outer radius rout = 5, reproduces the coaxial electrode arrangement analysed by Budd [6] for positive corona discharges. This geometry is fundamental to many industrial applications: high-voltage power transmis… view at source ↗
Figure 3
Figure 3. Figure 3: Solutions on the annular domain r ∈ [1, 5]. Top row: contour plots of ϕ, p and n. Bottom row: radial profiles showing ϕ(r) compared with the analytic Laplace solution ϕ(r) ∝ log(rout/r)/ log(rout/rin) 3.2 Axisymmetric Corona Test Many discharge configurations of practical relevance, such as needle-to-plane gaps, coaxial cylin￾ders, and curved-electrode coronas, possessing rotational symmetry. Exploiting ax… view at source ↗
Figure 4
Figure 4. Figure 4: shows the converged solution at V0 = 1.0 in the meridian section. The potential decays smoothly from the dome toward the floor and side wall. The positive ion density is sharply concentrated near the dome, where ions are injected, and decreases monotonically into the bulk. The electron density develops a weak interior maximum in a region of elevated electric field, where impact ionisation is most active [… view at source ↗
Figure 5
Figure 5. Figure 5: |∇ϕ| in the meridian section. Maximum is at the rim corner (r = 5, y = 10.55) due to the Dirichlet jump between A and C A secondary, weaker field channel connects the dome apex to the grounded floor along the symmetry axis, consistent with the expected Laplacian field structure. A three-dimensional ren￾dering of the solution for V0 = 1.5, obtained by revolving the meridian section, is provided in [PITH_FU… view at source ↗
Figure 6
Figure 6. Figure 6: The rim where the dome meets the cylinder constitutes a re-entrant edge with a Dirichlet discontinuity, analogous to the corner singularity in the L-shaped domain but in three dimensions. As expected from elliptic theory, the electric field exhibits a logarithmic singularity at this edge. The solver stably captures this feature, demonstrating its ability to handle practical 3D geometries with complex bound… view at source ↗
Figure 7
Figure 7. Figure 7: Voltage sweep V0 ∈ {0.5, 1.0, 1.5, 2.0}. Left: electric potential ϕ along the symmetry axis. Middle: electric field magnitude |∇ϕ| along the axis. Right: ion density p along the axis [PITH_FULL_IMAGE:figures/full_fig_p029_7.png] view at source ↗
read the original abstract

This work presents the analysis and numerical simulation of a stationary drift-diffusion model for electrical discharge in micro-electro-mechanical systems (MEMS). The model couples Poisson's equation for the electrostatic potential with continuity equations for positive ions and electrons, incorporating a Townsend-type impact ionization source term that depends exponentially on the electric field magnitude. We prove the existence of weak solutions under physically relevant assumptions and establish uniform bounds on the carrier densities. The proof relies on a regularization-approximation scheme with truncated nonlinearities, monotone operator theory (Browder-Minty), iterative energy estimates, and Stampacchia-type truncation arguments. We further develop a robust finite element solver to simulate the carrier density and electrostatic potential profiles for several geometries, including two-dimensional domains and a three-dimensional axisymmetric geometry.

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.

Referee Report

0 major / 3 minor

Summary. The manuscript analyzes a stationary drift-diffusion model for electrical discharge in MEMS, coupling Poisson's equation for the electrostatic potential with continuity equations for positive ions and electrons that include an exponential Townsend-type impact ionization source term. It proves existence of weak solutions under physically relevant assumptions via a regularization-approximation scheme, application of the Browder-Minty theorem to a truncated monotone operator, iterative L^2 energy estimates, and Stampacchia truncation to obtain uniform L^∞ bounds on the carrier densities. The work also develops and applies a finite element numerical solver to compute carrier density and potential profiles in 2D domains and a 3D axisymmetric geometry.

Significance. If the existence result and uniform bounds hold, the paper supplies a rigorous analytical foundation for a physically motivated model of electrical discharges in MEMS devices, which is a meaningful contribution to applied PDE analysis. The combination of monotone-operator theory with truncation arguments to control the exponential nonlinearity, together with the numerical implementation, strengthens the work. The approach relies on classical tools without apparent circularity or unverified growth conditions.

minor comments (3)
  1. [Abstract] Abstract: the phrase 'physically relevant assumptions' is used without enumeration; adding a short parenthetical list (e.g., bounded domain, positive coefficients, growth restrictions compatible with the exponential term) would improve immediate readability.
  2. [Existence proof (around the regularization step)] The truncation function applied to the exponential source term is described only at a high level; an explicit formula or reference to its precise definition in the approximation scheme would aid verification of the uniform bound passage to the limit.
  3. [Numerical results] Numerical section: the finite-element implementation would benefit from a brief statement of the mesh refinement strategy or solver tolerances used for the 2D and axisymmetric 3D examples, to support reproducibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary of our manuscript on the stationary drift-diffusion model for electrical discharges in MEMS, including the existence proof via regularization, monotone operators, and truncation arguments, as well as the finite element simulations. The recommendation for minor revision is noted; however, the report contains no specific major comments requiring response.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper establishes existence of weak solutions for the stationary drift-diffusion system via a regularization-approximation scheme, application of the Browder-Minty theorem to a truncated monotone operator, iterative L^2 energy estimates, and Stampacchia truncation to recover uniform L^∞ bounds on carrier densities. These steps invoke only classical, externally verifiable tools of nonlinear PDE theory whose validity does not depend on any fitted parameters, self-defined quantities, or prior results by the same authors. No load-bearing step reduces by construction to the paper's own inputs or to a self-citation chain; the derivation remains self-contained against standard mathematical benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review supplies no information on free parameters, background axioms, or new entities; the proof is described as relying on standard mathematical tools without further specification.

pith-pipeline@v0.9.1-grok · 5657 in / 1134 out tokens · 63228 ms · 2026-06-29T03:06:18.025981+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

60 extracted references

  1. [1]

    Abdel, A

    D. Abdel, A. Blaustein, C. Chainais-Hillairet, M. Herda, and J. Moatti. Existence of solutions and uniform bounds for the stationary semiconductor equations with generation and ionic carriers.arXiv preprint arXiv:2511.23250, 2025

    arXiv 2025

  2. [2]

    Brezzi, L

    F. Brezzi, L. D. Marini, S. Micheletti, P. Pietra, R. Sacco, and S. Wang. Discretization of semiconductor device problems (i).Handbook of numerical analysis, 13:317–441, 2005

    2005

  3. [3]

    Brezzi, L

    F. Brezzi, L. D. Marini, and P. Pietra. Two-dimensional exponential fitting and applications to drift-diffusion models.SIAM Journal on Numerical Analysis, 26(6):1342–1355, 1989

    1989

  4. [4]

    A. N. Brooks and T. J. Hughes. Streamline upwind/petrov-galerkin formulations for convec- tion dominated flows with particular emphasis on the incompressible navier-stokes equations. Computer methods in applied mechanics and engineering, 32(1-3):199–259, 1982

    1982

  5. [5]

    C. Budd. Coronas and the space charge problem.European Journal of Applied Mathematics, 2(1):43–81, 1991

    1991

  6. [6]

    C. Budd. Coronas and the space charge problem.European Journal of Applied Mathematics, 2:43–81, 1991. 36

    1991

  7. [7]

    Chainais-Hillairet, J.-G

    C. Chainais-Hillairet, J.-G. Liu, and Y.-J. Peng. Finite volume scheme for multi-dimensional drift-diffusion equations and convergence analysis.ESAIM: Mod´ elisation math´ ematique et analyse num´ erique, 37(2):319–338, 2003

    2003

  8. [8]

    Chen and H

    L. Chen and H. Bagci. Steady-state simulation of semiconductor devices using discontinuous galerkin methods.IEEE Access, 8:16203–16215, 2020

    2020

  9. [9]

    Degond, A

    P. Degond, A. J¨ ungel, and P. Pietra. Numerical discretization of energy-transport models for semiconductors with nonparabolic band structure.SIAM Journal on Scientific Computing, 22(3):986–1007, 2000

    2000

  10. [10]

    Degond, C

    P. Degond, C. D. Levermore, and C. Schmeiser. A note on the energy-transport limit of the semiconductor boltzmann equation. InTransport in Transition Regimes, pages 137–153. Springer, 2004

    2004

  11. [11]

    L. C. Evans.Partial differential equations, volume 19. American mathematical society, 2022

    2022

  12. [12]

    Fichera.Linear elliptic differential systems and eigenvalue problems, volume 8

    G. Fichera.Linear elliptic differential systems and eigenvalue problems, volume 8. Springer, 1965

    1965

  13. [13]

    Frehse and J

    J. Frehse and J. Naumann. An existence theorem for weak solutions of the basic stationary semiconductor equations.Applicable Analysis, 48(1-4):157–172, 1993

    1993

  14. [14]

    Frehse and J

    J. Frehse and J. Naumann. On the existence of weak solutions to a system of stationary semiconductor equations with avalanche generation.Mathematical Models and Methods in Applied Sciences, 4(02):273–289, 1994

    1994

  15. [15]

    Frehse and J

    J. Frehse and J. Naumann. Stationary semiconductor equations modeling avalanche genera- tion.Journal of mathematical analysis and applications, 198(3):685–702, 1996

    1996

  16. [16]

    Gajewski and K

    H. Gajewski and K. Gr¨ oger. On the basic equations for carrier transport in semiconductors. Journal of mathematical analysis and applications, 113(1):12–35, 1986

    1986

  17. [17]

    Gilbarg and N

    D. Gilbarg and N. S. Trudinger.Elliptic Partial Differential Equations of Second Order. Springer, Berlin, classics in mathematics edition, 2001

    2001

  18. [18]

    Gimperlein, R

    H. Gimperlein, R. He, and A. A. Lacey. Analysis of a model for electrical discharge in mems. arXiv preprint arXiv:2602.21439, 2026

    arXiv 2026

  19. [19]

    Grisvard.Singularities in Boundary Value Problems

    P. Grisvard.Singularities in Boundary Value Problems. Masson, Paris, 1992

    1992

  20. [20]

    H. K. Gummel. A self-consistent iterative scheme for one-dimensional steady state transistor calculations.IEEE Transactions on electron devices, 11(10):455–465, 2005

    2005

  21. [21]

    Herfst, P

    R. Herfst, P. Steeneken, J. Schmitz, A. Mank, and M. Van Gils. Kelvin probe study of laterally inhomogeneous dielectric charging and charge diffusion in rf mems capacitive switches. In 2008 IEEE International Reliability Physics Symposium, pages 492–495. IEEE, 2008

    2008

  22. [22]

    T. J. Hughes, G. Scovazzi, and L. P. Franca. Multiscale and stabilized methods.Encyclopedia of computational mechanics second edition, pages 1–64, 2018

    2018

  23. [23]

    R. S. Islamov. Global solutions of the drift-diffusion approximation equations in gas discharge theory.Differential Equations, 39(12), 2003. 37

    2003

  24. [24]

    R. S. Islamov. On the regularity of solutions to the drift-diffusion approximation equations in gas discharge theory.Computational Mathematics and Mathematical Physics, 46(1):125–140, 2006

    2006

  25. [25]

    J¨ ungel

    A. J¨ ungel. A positivity-preserving numerical scheme for a nonlinear fourth order parabolic system.SIAM journal on numerical analysis, 39(2):385–406, 2001

    2001

  26. [26]

    J¨ ungel.Quasi-hydrodynamic semiconductor equations, volume 41

    A. J¨ ungel.Quasi-hydrodynamic semiconductor equations, volume 41. Springer Science & Business Media, 2001

    2001

  27. [27]

    M. A. Krasnoselskij and P. P. Zabrejko.Geometrical methods of nonlinear analysis, volume

  28. [28]

    Kumar, M

    G. Kumar, M. Singh, A. Ray, and G. Trivedi. An fem based framework to simulate semi- conductor devices using streamline upwind petrov-galerkin stabilization technique. In2017 27th International Conference Radioelektronika (RADIOELEKTRONIKA), pages 1–5. IEEE, 2017

    2017

  29. [29]

    O. A. Ladyzhenskaya and N. N. Ural’tseva. Linear and quasilinear elliptic equations, 1968. Leon Ehrenpreis Academic Press, New York, 1968

    1968

  30. [30]

    J.-L. Lions. ” quelques m´ ethodes de r´ esolution des probl` emes aux limites non-lin´ eaires,”. Dunod, 1969

    1969

  31. [31]

    Liu and C.-W

    Y. Liu and C.-W. Shu. Analysis of the local discontinuous galerkin method for the drift- diffusion model of semiconductor devices.Science China Mathematics, 59(1):115–140, 2016

    2016

  32. [32]

    P. A. Markowich.The stationary semiconductor device equations. Springer Science & Business Media, 1985

    1985

  33. [33]

    P. A. Markowich, C. A. Ringhofer, and C. Schmeiser.Semiconductor equations. Springer Science & Business Media, 2012

    2012

  34. [34]

    J. J. H. Miller, W. Schilders, and S. Wang. Application of finite element methods to the simulation of semiconductor devices.Reports on Progress in Physics, 62(3):277–353, 1999

    1999

  35. [35]

    Molinero, R

    D. Molinero, R. Comulada, and L. Castaner. Dielectric charge measurements in capacitive microelectromechanical switches.Applied physics letters, 89(10), 2006

    2006

  36. [36]

    Molinero, R

    D. Molinero, R. Comulada, and L. Casta˜ ner. Dielectric charge measurements in capacitive microelectromechanical switches.Applied Physics Letters, 89:103506, 2006

    2006

  37. [37]

    Morrow and J

    R. Morrow and J. J. Lowke. Streamer propagation in air.Journal of Physics D: Applied Physics, 30(4):614–627, 1997

    1997

  38. [38]

    Morrow and J

    R. Morrow and J. J. Lowke. Streamer propagation in air.Journal of Physics D: Applied Physics, 30:614–627, 1997

    1997

  39. [39]

    Papaioannou, M.-N

    G. Papaioannou, M.-N. Exarchos, V. Theonas, G. Wang, and J. Papapolymerou. Temperature study of the dielectric polarization effects of capacitive rf mems switches.IEEE Transactions on Microwave Theory and Techniques, 53(11):3467–3473, 2005. 38

    2005

  40. [40]

    Patankar.Numerical heat transfer and fluid flow

    S. Patankar.Numerical heat transfer and fluid flow. CRC press, 2018

    2018

  41. [41]

    P´ erez-Escudero, D

    S. P´ erez-Escudero, D. Codony, I. Arias, and S. Fern´ andez-M´ endez. A comparison of formu- lations and non-linear solvers for computational modelling of semiconductor devices.Com- putational mechanics, 75(5):1533–1554, 2025

    2025

  42. [42]

    Polak, C

    S. Polak, C. Den Heijer, W. Schilders, and P. Markowich. Semiconductor device modelling from the numerical point of view.International Journal for Numerical Methods in Engineer- ing, 24(4):763–838, 1987

    1987

  43. [43]

    P. H. Rabinowitz. Some global results for nonlinear eigenvalue problems.Journal of functional analysis, 7(3):487–513, 1971

    1971

  44. [44]

    Y. P. Raizer.Gas Discharge Physics. Springer, Berlin, 1991

    1991

  45. [45]

    G. M. Rebeiz.RF MEMS: theory, design, and technology. John Wiley & Sons, 2004

    2004

  46. [46]

    Renardy and R

    M. Renardy and R. C. Rogers.An introduction to partial differential equations. Springer, 2004

    2004

  47. [47]

    P. J. Roache.Verification and Validation in Computational Science and Engineering. Her- mosa, Albuquerque, 1998

    1998

  48. [48]

    H.-G. Roos, M. Stynes, and L. Tobiska.Robust numerical methods for singularly perturbed differential equations: convection-diffusion-reaction and flow problems. Springer, 2008

    2008

  49. [49]

    D. L. Scharfetter and H. K. Gummel. Large-signal analysis of a silicon read diode oscillator. IEEE Transactions on electron devices, 16(1):64–77, 2005

    2005

  50. [50]

    Selberherr.Analysis and simulation of semiconductor devices

    S. Selberherr.Analysis and simulation of semiconductor devices. Springer Science & Business Media, 1984

    1984

  51. [51]

    S. D. Senturia.Microsystem Design. Kluwer, Boston, 2001

    2001

  52. [52]

    Stampacchia

    G. Stampacchia. Le probl` eme de dirichlet pour les ´ equations elliptiques du second ordre ` a coefficients discontinus. InAnnales de l’institut Fourier, volume 15, pages 189–257, 1965

    1965

  53. [53]

    Townsend.The theory of ionization of gases by collision

    J. Townsend.The theory of ionization of gases by collision. Constable, Limited, 1910

    1910

  54. [54]

    J. S. Townsend.Electricity in gases. Clarendon Press, Oxford, 1915

    1915

  55. [55]

    Van Roosbroeck

    W. Van Roosbroeck. Theory of the flow of electrons and holes in germanium and other semiconductors.The Bell System Technical Journal, 29(4):560–607, 1950

    1950

  56. [56]

    Van Spengen

    W. Van Spengen. Capacitive rf mems switch dielectric charging and reliability: a critical re- view with recommendations.Journal of Micromechanics and Microengineering, 22(7):074001, 2012

    2012

  57. [57]

    Vasileska, S

    D. Vasileska, S. M. Goodnick, and G. Klimeck.Computational Electronics: semiclassical and quantum device modeling and simulation. CRC press, 2017

    2017

  58. [58]

    Xiao.Gas discharge and gas insulation

    D. Xiao.Gas discharge and gas insulation. Springer, 2016. 39

    2016

  59. [59]

    X. Yuan, J. C. Hwang, D. Forehand, and C. L. Goldsmith. Modeling and characterization of dielectric-charging effects in rf mems capacitive switches. InIEEE MTT-S International Microwave Symposium Digest, 2005., pages 753–756. IEEE, 2005

    2005

  60. [60]

    O. C. Zienkiewicz, R. L. Taylor, and J. Z. Zhu.The Finite Element Method: Its Basis and Fundamentals. Butterworth-Heinemann, Oxford, 7th edition, 2013. 40

    2013