Thermodynamic Charge Partition in Accumulation-Layer Heterostructures
Pith reviewed 2026-05-09 19:08 UTC · model grok-4.3
The pith
Treating the partition of induced sheet density as the central state variable yields a complete thermodynamic description of accumulation-layer heterostructures.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By treating the charge partition between the accumulation layer and the screening charge as the central state variable, one obtains a complete Helmholtz free energy. This leads to a corrected locked-branch chemical potential and a shifted release potential that separates the selection of the physical path from the geometric capacitance. The path selection occurs spectrally according to compressibility: compressible segments stay fully screened while incompressible segments follow the locked branch until the relevant gap triggers release. As a result, quantities like differential capacitance, tunnel current, and plateau width appear as distinct projections of this single coupled thermodynamic
What carries the argument
The partition of the induced sheet density into accumulation-layer charge and complementary screening charge, serving as the central state variable for deriving the Helmholtz free energy and selecting thermodynamic paths based on spectral compressibility.
Load-bearing premise
The charge partition between the accumulation layer and screening charge can be treated as the central state variable with paths selected spectrally by compressibility and gaps.
What would settle it
Measuring a constant screening depth independent of magnetic field in magnetocapacitance experiments on accumulation-layer heterostructures would contradict the predicted growth with field.
Figures
read the original abstract
We develop a thermodynamic description of accumulation-layer heterostructures in which the induced sheet density is partitioned between the near-interface accumulation-layer charge and a complementary screening charge in the surrounding structure. Treating this partition as the central state variable yields a complete Helmholtz free energy, a corrected locked-branch chemical potential, and a shifted release potential that separates energetic path selection from geometric capacitance. The physical path is selected spectrally: compressible segments remain fully screened, whereas incompressible segments evolve along a locked branch until release is triggered by the relevant gap. Differential capacitance, tunnel current and plateau width then emerge as different projections of the same coupled thermodynamic structure. A canonical two-stage self-consistent Poisson--Schr\"odinger reduction supplies universal master functions for the isolated accumulation layer and master surfaces for its finite-buffer extension, making the theory calculable across density and geometry. Comparison with magnetocapacitance and magnetotunneling data supports a picture in which nearby extended charge refills the accumulation layer and the effective screening depth grows with magnetic field.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a thermodynamic description of accumulation-layer heterostructures in which the induced sheet density is partitioned between near-interface accumulation-layer charge and complementary screening charge in the surrounding structure. Treating this partition as the central state variable produces a complete Helmholtz free energy, a corrected locked-branch chemical potential, and a shifted release potential that separates energetic path selection from geometric capacitance. The physical path is selected spectrally according to compressibility (compressible segments fully screened; incompressible segments evolve along a locked branch until gap-triggered release). A canonical two-stage self-consistent Poisson-Schrödinger reduction supplies universal master functions for the isolated layer and master surfaces for finite-buffer extensions. Projections onto differential capacitance, tunnel current, and plateau width are compared quantitatively to magnetocapacitance and magnetotunneling data, supporting a picture of nearby extended charge refilling the accumulation layer with field-dependent screening depth.
Significance. If the construction holds, the work supplies a unified thermodynamic structure that cleanly separates energetic path selection from geometric capacitance via spectral selection rules. The explicit derivation of parameter-light master functions from the two-stage Poisson-Schrödinger reduction and the reported quantitative agreement with refilling and screening-depth data constitute concrete strengths that could be useful for modeling accumulation layers in the quantum Hall regime.
major comments (2)
- [Thermodynamic construction] The central claim that the charge-partition variable yields an independent Helmholtz free energy and locked-branch chemical potential is load-bearing. The manuscript should supply the explicit differential relations and boundary conditions used to obtain these quantities (theory section following the abstract) so that readers can verify that the construction does not reduce to a post-hoc fit when the spectral selection is applied to the magnetocapacitance data.
- [Poisson-Schrödinger reduction] The two-stage self-consistent Poisson-Schrödinger reduction is asserted to produce universal master functions independent of geometry-specific fitting. An explicit statement of the reduction steps and the conditions under which the functions remain parameter-free (e.g., in the isolated-layer limit versus finite-buffer extension) would directly address the risk that data agreement arises from adjustable path selection rather than from the thermodynamic structure.
minor comments (2)
- [Abstract] The abstract is information-dense; a single additional sentence clarifying the operational meaning of the 'shifted release potential' would improve accessibility without lengthening the paragraph.
- [Data comparison] Figure captions for the magnetocapacitance and magnetotunneling comparisons should list the precise density and field ranges over which the master functions are evaluated and whether any overall scale factor is applied.
Simulated Author's Rebuttal
We thank the referee for the positive assessment and the constructive comments aimed at improving the clarity of the thermodynamic construction. We respond to each major comment below and will incorporate the requested explicit derivations in a revised manuscript.
read point-by-point responses
-
Referee: [Thermodynamic construction] The central claim that the charge-partition variable yields an independent Helmholtz free energy and locked-branch chemical potential is load-bearing. The manuscript should supply the explicit differential relations and boundary conditions used to obtain these quantities (theory section following the abstract) so that readers can verify that the construction does not reduce to a post-hoc fit when the spectral selection is applied to the magnetocapacitance data.
Authors: We agree that the differential relations should be stated more explicitly to allow direct verification. In the revised manuscript we will add, immediately after the abstract in the theory section, the complete set of relations: the Helmholtz free energy F(n_p, A) constructed from the partition variable n_p = n_acc / n_tot, its differential dF = μ_locked dn_p + ... with the Legendre transform that isolates the locked-branch chemical potential, the boundary condition fixing the total gate-induced density n_tot = n_acc + n_screen, and the spectral selection rule that enforces the locked branch until gap closure. These relations are derived variationally from the free-energy minimization subject to the compressibility criterion and are independent of the subsequent data comparison; the added subsection will make this derivation self-contained. revision: yes
-
Referee: [Poisson-Schrödinger reduction] The two-stage self-consistent Poisson-Schrödinger reduction is asserted to produce universal master functions independent of geometry-specific fitting. An explicit statement of the reduction steps and the conditions under which the functions remain parameter-free (e.g., in the isolated-layer limit versus finite-buffer extension) would directly address the risk that data agreement arises from adjustable path selection rather than from the thermodynamic structure.
Authors: We accept that a step-by-step enumeration will remove any ambiguity. In the revised Section III we will explicitly list the two stages: (i) the isolated-layer reduction, in which the Schrödinger equation is solved self-consistently with Poisson for an infinite buffer, producing master functions that contain only fundamental constants (effective mass, dielectric constant, magnetic length) and no geometry-specific parameters; (ii) the finite-buffer extension obtained by perturbative boundary matching at the buffer edge, which introduces a single geometric length but leaves the master functions otherwise unchanged. The isolated-layer limit is strictly parameter-free; the finite-buffer case adds only the known buffer thickness. This structure ensures that quantitative agreement with magnetocapacitance data follows from the thermodynamic and spectral rules rather than from adjustable path selection. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper's central construction begins by positing the induced-sheet-density partition as the thermodynamic state variable and then derives the Helmholtz free energy, locked-branch chemical potential, and release potential from its differentials. A standard two-stage self-consistent Poisson-Schrödinger reduction is applied to obtain explicit universal master functions for the isolated layer and master surfaces for finite buffers; these functions are presented as calculable outputs rather than inputs. Differential capacitance, tunnel current, and plateau widths are obtained as projections of the same structure. Quantitative comparison with external magnetocapacitance and magnetotunneling data is performed after the derivation, not used to define the master functions themselves. No equation reduces to a prior result by construction, no self-citation supplies a load-bearing uniqueness theorem, and no fitted parameter is relabeled as a prediction. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The induced sheet density partitions between near-interface accumulation-layer charge and complementary screening charge in the surrounding structure.
- domain assumption Compressible segments remain fully screened while incompressible segments follow a locked branch until release triggered by the relevant gap.
Reference graph
Works this paper leans on
-
[1]
Overview, motivation, and strategy The logic of the appendix is a three-step progression
-
[2]
The canonical Airy problem of the triangular well suggests that the accumulation layer should itself admit a scaled, universal description. Stage 1 shows that this remains true after full self-consistent Poisson– Schrödinger feedback: the isolated accumulation layer collapses to a one-parameter family on the(nt, κ)grid, encoded by master functions rather ...
-
[3]
The numerical stage-1 solutions immediately suggest a second step. The bound ground state is sharply localized near the interface, whereas the higher states are orthogonal to it, live at higher energies, and extend deep into the buffer. This motivates splitting the total transferred charge into a bound accumulation-layer component κnt and an extended comp...
-
[4]
The central result of stage 2 is that this enlarged problem again admits stable canonical solutions. The finite-buffer multisubband system does not produce a second pair of fitted master curves, but a tabulated family of master surfaces on the(κ, λ)grid. These surfaces are the universal input for everything that follows on that grid: the free energy, the ...
-
[5]
Canonical statement of the finite-buffer multisubband problem The most general canonical formulation of stage 2 is obtained by embedding the accumulation layer into buffer 1 of physical widthL and canonical widthλ = L/ℓt. 19 Buffer 1 is uniformly doped with donor densityND, for which the natural dimensionless measure is ηD = NDℓt nt . The physically compl...
-
[6]
Stage 1: canonical isolated accumulation layer Stage 1 is the direct self-consistent generalization of the canonical Airy problem. The transferred sheet charge is restricted to the isolated accumulation layer, so the Poisson–Schrödingerproblemcollapsestoaone-parameter family indexed by the screening fractionκ. In this sense stage 1 solves the self-consist...
-
[7]
Stage 2: finite-buffer lift with extended states With the full finite-buffer boundary-value problem now stated, the second canonicalization consists in lifting the stage-1 accumulation-layer solution into the multisubband problem on the finite interval0 ≤s≤λ . The logic mirrors stage 1: one fixes the bound component byn0 = κnt, distributes the remaining c...
-
[8]
Start from an initial guessu(0)(s), typically the linear triangular profile or the converged solution at a nearby map point
-
[9]
Diagonalize the discretized Schrödinger operator to obtainε (m) i andχ (m) i at iterationm
-
[10]
Determine the extended occupations n(m) i from Eq. (A11)
-
[11]
Build the total electronic density and solve Poisson’s equation for the updated potential
-
[12]
Apply clamped linear mixing, u(m+1) ←(1−α)u (m) +α u (m) new,0< α≪1, to suppress oscillatory convergence in the high- κ regime. For the dense production runs used here, convergence typically requires of order102 iterations, with represen- tative difficult points nearκ→ 1taking roughly ∼ 200 iterations. The resulting canonical map is produced on a dense re...
-
[13]
Microscopic anatomy of the stage-2 solution The canonical stage-2 map is more than a repository of fitted functions. It also exposes the internal structure of the self-consistent buffer solution and thereby gives direct microscopic meaning to the thermodynamic partition of Eq. (2). The accumulated component nE is carried by 22 the interfacial ground state...
-
[14]
Canonical response surfaces and interface–bulk decomposition The stage-2 self-consistent map on the(κ, λ)grid is valuable not only because it solves the finite-buffer prob- lem numerically, but because it organizes the result into a small set of canonical response surfaces. These out- puts separate into near-interface quantities, which are controlled loca...
-
[15]
Canonical character and universality The formulation is canonical in the following precise sense
-
[16]
The isolated accumulation layer defines a universal family of self-consistent solutions on the(nt, κ)grid, all dependence onnt being carried by the similarity scalesℓ t andE t
-
[17]
Buffer 1 then promotes this local family to a finite- region map on the( κ, λ)grid, with λ the canoni- cal width of the region in which the residual field is screened
-
[18]
Buffer 2 and the barrier are excluded at this stage, and the quasimetallic electrodes enter only as external boundaries, because these smooth closure elements do not modify the canonical κ-dependent physics; they only shift the absolute measured voltage and the capacitance baseline. The canonical PS problem is therefore a transferable math- ematical struc...
-
[19]
Uniqueness, stability, and limits of the stage-2 lift A foundational structural property of the fully screened canonical self-consistent Poisson–Schrödinger problem is 25 that thestage-1 accumulation-layer potential at κ = 1 supports exactly one bound eigenstate. This is the configu- ration in which the entire transferred sheet density resides in the accu...
-
[20]
T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. Phys. 54, 437 (1982)
work page 1982
-
[21]
J. H. Davies,The Physics of Low-Dimensional Semicon- ductors: An Introduction(Cambridge University Press, Cambridge, 1998)
work page 1998
- [22]
- [23]
-
[24]
N. D. Mermin, Phys. Rev.137, A1441 (1965)
work page 1965
- [25]
-
[26]
E. H. Lieb, Int. J. Quantum Chem.24, 243 (1983)
work page 1983
- [27]
- [28]
-
[29]
D. B. Chklovskii, B. I. Shklovskii, and L. I. Glazman, Phys. Rev. B46, 4026 (1992)
work page 1992
-
[30]
U. Wulf, V. Gudmundsson, and R. R. Gerhardts, Phys. Rev. B38, 4218 (1988)
work page 1988
- [31]
- [32]
-
[33]
T. P. Smith, B. B. Goldberg, P. J. Stiles, and M. Heiblum, Phys. Rev. B32, 2696 (1985)
work page 1985
- [34]
-
[35]
J. P. Eisenstein, L. N. Pfeiffer, and K. W. West, Phys. Rev. B50, 1760 (1994)
work page 1994
-
[36]
C. Berthod, H. Zhang, A. F. Morpurgo, and T. Giamarchi, Phys. Rev. Research3, 043036 (2021)
work page 2021
-
[37]
E. Böckenhoff, K. von Klitzing, and K. Ploog, Phys. Rev. B38, 10120 (1988)
work page 1988
-
[38]
Böckenhoff,Electrical Transport Perpendicular to Semi- conductor Heterostructures, Ph.D
E. Böckenhoff,Electrical Transport Perpendicular to Semi- conductor Heterostructures, Ph.D. thesis, Max-Planck- Institut für Festkörperforschung, Stuttgart (1988)
work page 1988
-
[39]
S. Bending, C. Zhang, K. von Klitzing, and K. Ploog, Phys. Rev. B39, 1097 (1989)
work page 1989
-
[40]
T. W. Hickmott, Phys. Rev. B40, 8363 (1989)
work page 1989
-
[41]
K. von Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett.45, 494 (1980)
work page 1980
-
[42]
R. E. Prange and S. M. Girvin, eds.,The Quantum Hall Effect, 2nd ed. (Springer, New York, 1990)
work page 1990
-
[43]
S. Yi-Thomas, Y. Huang, J. D. Sau, and S. Das Sarma, Phys. Rev. B111, 195305 (2025)
work page 2025
- [44]
-
[45]
G. J. Iafrate and K. Hess, Phys. Rev. B39, 1955 (1989)
work page 1955
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.