Electrostatics in semiconducting devices I : The Pure Electrostatics Self Consistent Approximation
Pith reviewed 2026-05-23 02:23 UTC · model grok-4.3
The pith
PESCA provides a quantitative electrostatic model for charge distributions in semiconductor devices when the capacitance ratio κ is small.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
PESCA is a self-consistent electrostatic approximation that includes the semiconductor response through a quantum-capacitance screening term, allowing calculation of charge distributions inside devices. The validity of this model is controlled by the small parameter κ = C_g/C_q, which is of order 1% in many situations and thereby renders PESCA quantitative.
What carries the argument
The Pure Electrostatic Self Consistent Approximation (PESCA), which adds semiconductor screening and depletion to electrostatic calculations via a quantum-capacitance term.
If this is right
- Reconstructs the charge distribution inside devices from measurements of pinch-off phase diagrams in gate voltage space.
- Extends directly to include magnetic fields and calculate edge reconstruction in the quantum Hall regime.
- Yields quantitative results for charge calculations whenever κ is of order 1% as occurs in many common device geometries.
- Supplies a minimum model for adding semiconductors to electrostatic calculations while retaining screening and depletion.
Where Pith is reading between the lines
- Device modeling workflows could adopt PESCA to simplify calculations of charge profiles before adding more complex quantum corrections.
- The same capacitance-ratio control might apply to electrostatic treatments of other confined carrier systems beyond the devices discussed.
- Experimental mapping of κ across different gate geometries would directly test the range where PESCA remains accurate.
Load-bearing premise
The semiconductor response can be captured entirely by a quantum-capacitance screening term inside a purely electrostatic calculation without additional quantum-mechanical or disorder effects.
What would settle it
Measurements of charge distributions in devices that deviate substantially from PESCA predictions even when the measured value of κ remains small.
Figures
read the original abstract
In quantum nanoelectronics devices, the electrostatic energy is the largest energy scale at play and, to a large extend, it determines the charge distribution inside the devices. Here, we introduce the Pure Electrostatic Self consistent Approximation (PESCA) that provides a minimum model that describes how to include a semiconductor in an electrostatic calculation to properly account for both screening and partial depletion due to e.g. field effect. We show how PESCA may be used to reconstruct the charge distribution from the measurement of pinch-off phase diagrams in the gate voltages space. PESCA can also be extended to account for magnetic field and calculate the edge reconstruction in the quantum Hall regime. The validity of PESCA is controlled by a small parameter $\kappa = C_g/C_q$, the ratio of the geometrical capacitance to the quantum capacitance, which is, in many common situations, of the order of 1%, making PESCA a quantitative technique for the calculation of the charge distribution inside devices.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces the Pure Electrostatic Self-Consistent Approximation (PESCA) as a minimum model for incorporating semiconductor screening and partial depletion into electrostatic calculations of charge distributions in quantum nanoelectronics devices. PESCA is controlled by the small parameter κ = C_g/C_q, claimed to be ≪1 (often ~1%) in common situations, enabling quantitative use for reconstructing charge from pinch-off phase diagrams in gate-voltage space and for extensions to magnetic fields and quantum-Hall edge reconstruction.
Significance. If the central modeling step holds, PESCA supplies an explicitly controlled, low-parameter electrostatic framework whose error is bounded by a measurable capacitance ratio. This could streamline device modeling and data interpretation in regimes where electrostatics dominate, with the parameter-free character of the approximation (once κ is fixed externally) constituting a clear methodological strength.
major comments (1)
- [Abstract, first paragraph] The replacement of the full semiconductor response by a local quantum-capacitance term is the load-bearing modeling choice; the manuscript should supply a concrete derivation or bound showing that additional quantum-mechanical or disorder contributions remain higher-order in κ (or are negligible in the targeted regime).
minor comments (2)
- [Abstract] The abstract contains a typographical error: 'to a large extend' should read 'to a large extent'.
- [Abstract] Notation for capacitances (C_g, C_q) is introduced without an explicit definition or diagram; a short schematic relating geometrical and quantum capacitances would improve clarity for readers outside the immediate subfield.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the recommendation of minor revision. We respond to the major comment below.
read point-by-point responses
-
Referee: [Abstract, first paragraph] The replacement of the full semiconductor response by a local quantum-capacitance term is the load-bearing modeling choice; the manuscript should supply a concrete derivation or bound showing that additional quantum-mechanical or disorder contributions remain higher-order in κ (or are negligible in the targeted regime).
Authors: We agree that the replacement of the full semiconductor response by a local quantum-capacitance term constitutes the central modeling step of PESCA. The manuscript derives this approximation from the electrostatic energy functional under the assumption that the dominant semiconductor response is captured by the quantum capacitance C_q, yielding the small parameter κ = C_g/C_q. Deviations arising from non-local quantum-mechanical effects or disorder are argued to enter at higher order in κ because they are suppressed by the same capacitance ratio that controls the leading approximation. To make this bound explicit as requested, we will add a concise derivation (approximately one paragraph plus a short appendix) in the revised manuscript that estimates the leading correction terms as O(κ) or smaller in the regime of interest. This addition will clarify the error control without altering the main results, figures, or conclusions. revision: yes
Circularity Check
No significant circularity identified
full rationale
The manuscript introduces PESCA as an explicit modeling approximation that replaces the semiconductor response with a local quantum-capacitance term inside a purely electrostatic solver. Its claimed validity is controlled by the external small parameter κ = C_g/C_q, presented as a physical ratio rather than a quantity fitted or defined inside the model. No derivation step reduces a prediction to a fitted input by construction, no self-citation chain is load-bearing, and no ansatz is smuggled via prior work. The central claim is therefore a controlled modeling choice whose error estimate is stated independently of the approximation itself.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The electrostatic energy is the largest energy scale at play and to a large extent determines the charge distribution inside the devices.
Forward citations
Cited by 1 Pith paper
-
Electrostatics in semiconducting devices II: Solving the Helmholtz equation
Maps the quantum-electrostatic self-consistency problem to a non-linear Helmholtz equation, enabling construction of a convex functional for provable convergence that lifts to the exact solution in typically one or tw...
Reference graph
Works this paper leans on
-
[1]
A. Lacerda-Santos and X. Waintal, Electrostatics in semiconducting devices ii : Solving the helmholtz equation, In preparation (2025)
work page 2025
-
[2]
A. Lacerda-Santos, C. Groth and X. Waintal, Electrostatics in semiconducting devices iii : The pescado open source library, In preparation (2025)
work page 2025
-
[3]
D. B. Chklovskii, B. I. Shklovskii and L. I. Glazman,Electrostatics of edge channels, Physical Review B 46(7), 4026 (1992), doi:10.1103 /PhysRevB.46.4026
work page 1992
-
[4]
D. B. Chklovskii, K. A. Matveev and B. I. Shklovskii, Ballistic conductance of in- teracting electrons in the quantum hall regime , Phys. Rev. B 47, 12605 (1993), doi:10.1103/PhysRevB.47.12605
-
[5]
D. B. Chklovskii, B. I. Shklovskii and L. I. Glazman, Erratum: Electrostatics of edge chan- nels, Phys. Rev. B 46, 15606 (1992), doi:10.1103 /PhysRevB.46.15606.3
work page 1992
-
[6]
P . Armagnat and X. Waintal, Reconciling edge states with compressible stripes in a ballistic mesoscopic conductor , Journal of Physics: Materials 3(2), 02LT01 (2020), doi:10.1088/2515-7639/ab7582
-
[7]
P . Armagnat, A. Lacerda-Santos, B. Rossignol, C. Groth and X. Waintal,The self-consistent quantum-electrostatic problem in strongly non-linear regime, SciPost Phys. 7, 031 (2019), doi:10.21468/SciPostPhys.7.3.031
-
[8]
E. Chatzikyriakou, J. Wang, L. Mazzella, A. Lacerda-Santos, M. C. d. S. Figueira, A. Trel- lakis, S. Birner, T . Grange, C. Bäuerle and X. Waintal,Unveiling the charge distribution of a gaas-based nanoelectronic device: A large experimental dataset approach, Phys. Rev. Res. 4, 043163 (2022), doi:10.1103 /PhysRevResearch.4.043163
work page 2022
-
[9]
H. Edlbauer, J. Wang, T . Crozes, P . Perrier, S. Ouacel, C. Geffroy , G. Georgiou, E. Chatzikyr- iakou, A. Lacerda-Santos, X. Waintal, D. C. Glattli, P . Roulleauet al., Semiconductor-based electron flying qubits: review on recent progress accelerated by numerical modelling , EPJ Quantum Technology 9(1), 21 (2022), doi:10.1140 /epjqt/s40507-022-00139-w
work page 2022
-
[10]
G. Roussely , E. Arrighi, G. Georgiou, S. Takada, M. Schalk, M. Urdampilleta, A. Ludwig, A. D. Wieck, P . Armagnat, T . Kloss, X. Waintal, T . Meunieret al., Unveiling the bosonic nature of an ultrashort few-electron pulse , Nature Communications 9(1), 2811 (2018), doi:10.1038/s41467-018-05203-7
-
[11]
B. D. Woods, T . D. Stanescu and S. Das Sarma,Effective theory approach to the schrödinger- poisson problem in semiconductor majorana devices , Phys. Rev. B 98, 035428 (2018), doi:10.1103/PhysRevB.98.035428
-
[12]
B. D. Woods, S. Das Sarma and T . D. Stanescu, Charge-impurity ef- fects in hybrid majorana nanowires , Phys. Rev. Appl. 16, 054053 (2021), doi:10.1103/PhysRevApplied.16.054053
-
[13]
S. Ahn, H. Pan, B. Woods, T . D. Stanescu and S. Das Sarma, Estimating disorder and its adverse effects in semiconductor majorana nanowires, Phys. Rev. Mater.5, 124602 (2021), doi:10.1103/PhysRevMaterials.5.124602
-
[14]
J. Wang, E. Polizzi and M. Lundstrom, A three-dimensional quantum simulation of silicon nanowire transistors with the effective-mass approximation , Journal of Applied Physics 96(4), 2192 (2004), doi:10.1063 /1.1769089. 22 SciPost Physics Submission
work page 2004
-
[15]
P . Marconcini, A. Cresti, F . Triozon, G. Fiori, B. Biel, Y.-M. Niquet, M. Macucci and S. Roche, Atomistic boron-doped graphene field-effect transistors: A route toward unipo- lar characteristics , ACS Nano 6(9), 7942 (2012), doi:10.1021 /nn3024046, PMID: 22876866
work page 2012
-
[16]
J. D. T . Luna, S. R. Kuppuswamy and A. R. Akhmerov,Design of a Majorana trijunction, SciPost Phys. 16, 044 (2024), doi:10.21468 /SciPostPhys.16.2.044
work page 2024
- [17]
-
[18]
S. D. Escribano, E. Prada, Y. Oreg and A. L. Yeyati,Tunable proximity effects and topological superconductivity in ferromagnetic hybrid nanowires, Phys. Rev. B 104, L041404 (2021), doi:10.1103/PhysRevB.104.L041404
-
[19]
T . Karzig, C. Knapp, R. M. Lutchyn, P . Bonderson, M. B. Hastings, C. Nayak, J. Alicea, K. Flensberg, S. Plugge, Y. Oreg, C. M. Marcus and M. H. Freedman, Scalable designs for quasiparticle-poisoning-protected topological quantum computation with majorana zero modes, Phys. Rev. B 95, 235305 (2017), doi:10.1103 /PhysRevB.95.235305
work page 2017
-
[20]
B. D. Woods and T . D. Stanescu, Electrostatic effects and topological superconductivity in semiconductor–superconductor–magnetic-insulator hybrid wires, Phys. Rev. B104, 195433 (2021), doi:10.1103 /PhysRevB.104.195433
work page 2021
-
[21]
C.-X. Liu, S. Schuwalow, Y. Liu, K. Vilkelis, A. L. R. Manesco, P . Krogstrup and M. Wimmer, Electronic properties of inas/eus/al hybrid nanowires, Phys. Rev. B 104, 014516 (2021), doi:10.1103/PhysRevB.104.014516
-
[22]
H. Gummel, A self-consistent iterative scheme for one-dimensional steady state transistor calculations, IEEE Transactions on Electron Devices11(10), 455 (1964), doi:10.1109/T- ED.1964.15364
work page doi:10.1109/t- 1964
-
[23]
I. Tan, G. L. Snider, L. D. Chang and E. L. Hu, A self-consistent solution of Schrödinger–Poisson equations using a nonuniform mesh , Journal of Applied Physics 68(8), 4071 (1990), doi:10.1063 /1.346245
work page 1990
-
[24]
V . Gudmundsson,Oscillating impurity spectra caused by non-linear screening in the quan- tum hall regime , Solid State Communications 74(2), 63 (1990), doi:10.1016 /0038- 1098(90)90606-C
work page 1990
-
[25]
L. Wang, D. Wang and P . M. Asbeck,A numerical schrödinger–poisson solver for radially symmetric nanowire core–shell structures, Solid-State Electronics 50(11), 1732 (2006), doi:10.1016/j.sse.2006.09.013
-
[26]
F . Rana, S. Tiwari and D. A. Buchanan,Self-consistent modeling of accumulation layers and tunneling currents through very thin oxides, Applied Physics Letters 69(8), 1104 (1996), doi:10.1063/1.117072
-
[27]
F . Stern, Iteration methods for calculating self-consistent fields in semiconductor inver- sion layers , Journal of Computational Physics 6(1), 56 (1970), doi:10.1016 /0021- 9991(70)90004-5
work page 1970
-
[28]
A. Trellakis, A. T . Galick, A. Pacelli and U. Ravaioli,Iteration scheme for the solution of the two-dimensional Schrödinger-Poisson equations in quantum structures, Journal of Applied Physics 81(12), 7880 (1997), doi:10.1063 /1.365396. 23 SciPost Physics Submission
work page 1997
-
[29]
Journal of Computational Physics , author =
V . Eyert, A comparative study on methods for convergence acceleration of iter- ative vector sequences , Journal of Computational Physics 124(2), 271 (1996), doi:10.1006/jcph.1996.0059
-
[30]
A. Vuik, D. Eeltink, A. R. Akhmerov and M. Wimmer, Effects of the electrostatic environ- ment on the majorana nanowire devices, New Journal of Physics 18(3), 033013 (2016), doi:10.1088/1367-2630/18/3/033013
-
[31]
D. G. Anderson, Iterative procedures for nonlinear integral equations , J. ACM 12(4), 547–560 (1965), doi:10.1145/321296.321305
-
[32]
A. E. Antipov, A. Bargerbos, G. W . Winkler, B. Bauer, E. Rossi and R. M. Lutchyn,Effects of gate-induced electric fields on semiconductor majorana nanowires, Phys. Rev. X 8, 031041 (2018), doi:10.1103 /PhysRevX.8.031041
work page 2018
-
[33]
S. D. Escribano, A. Levy Yeyati, Y. Oreg and E. Prada, Effects of the electrostatic en- vironment on superlattice majorana nanowires , Phys. Rev. B 100, 045301 (2019), doi:10.1103/PhysRevB.100.045301
-
[34]
Pulay ,Convergence acceleration of iterative sequences
P . Pulay ,Convergence acceleration of iterative sequences. the case of scf iteration, Chemical Physics Letters 73(2), 393 (1980), doi:10.1016 /0009-2614(80)80396-4
work page 1980
-
[35]
C. G. Broyden, A class of methods for solving nonlinear simultaneous equations , Math. Comp. 19, 577 (1965), doi:10.1090 /S0025-5718-1965-0198670-6
work page 1965
-
[36]
P . Andrei and I. Mayergoyz, Analysis of fluctuations in semiconductor devices through self-consistent Poisson-Schrödinger computations, Journal of Applied Physics 96(4), 2071 (2004), doi:10.1063 /1.1772886
work page 2071
-
[37]
R. Lake, G. Klimeck, R. C. Bowen, D. Jovanovic, D. Blanks and M. Swaminathan, Quan- tum transport with band-structure and schottky contacts, physica status solidi (b) 204(1), 354 (1997), doi:10.1002 /1521-3951(199711)204:1<354::AID-PSSB354>3.0.CO;2-
work page 1997
-
[38]
A. Pacelli, Self-consistent solution of the schrodinger equation in semiconductor de- vices by implicit iteration , IEEE Transactions on Electron Devices 44(7), 1169 (1997), doi:10.1109/16.595946
- [39]
-
[40]
Journal of Computational Physics , author =
D. Knoll and D. Keyes, Jacobian-free newton–krylov methods: a survey of ap- proaches and applications , Journal of Computational Physics 193(2), 357 (2004), doi:10.1016/j.jcp.2003.08.010
-
[41]
A. E. G. Mikkelsen, P . Kotetes, P . Krogstrup and K. Flensberg, Hybridiza- tion at superconductor-semiconductor interfaces , Phys. Rev. X 8, 031040 (2018), doi:10.1103/PhysRevX.8.031040
-
[42]
A. Trellakis and U. Ravaioli, Lateral scalability limits of silicon conduction channels, Jour- nal of Applied Physics 86(7), 3911 (1999), doi:10.1063 /1.371307
work page 1999
-
[43]
G. Curatola and G. Iannaccone, Nanotcad2d: Two-dimensional code for the simulation of nanoelectronic devices and structures , Computational Materials Science 28(2), 342 (2003), doi:10.1016 /S0927-0256(03)00117-4. 24 SciPost Physics Submission
work page 2003
- [44]
-
[45]
T . Bautze, C. Süssmeier, S. Takada, C. Groth, T . Meunier, M. Yamamoto, S. Tarucha, X. Waintal and C. Bäuerle, Theoretical, numerical, and experimental study of a flying qubit electronic interferometer , Phys. Rev. B 89, 125432 (2014), doi:10.1103/PhysRevB.89.125432
-
[46]
C. Bäuerle, D. Christian Glattli, T . Meunier, F . Portier, P . Roche, P . Roulleau, S. Takada and X. Waintal, Coherent control of single electrons: a review of current progress , Reports on Progress in Physics 81(5), 056503 (2018), doi:10.1088 /1361-6633/aaa98a
work page 2018
-
[47]
Bäuerle, F .Pierre and P .Roche, Private communications
C. .Bäuerle, F .Pierre and P .Roche, Private communications . 25
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.