pith. sign in

arxiv: 2501.17285 · v2 · submitted 2025-01-28 · ❄️ cond-mat.mes-hall

Nonequilibrium Green's Function Formalism Applicable to Discrete Impurities in Semiconductor Nanostructures

Nobuyuki Sano This is my paper

Pith reviewed 2026-05-23 04:44 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall
keywords nonequilibrium Green's functionimpurity scatteringsemiconductor nanostructuresnonlocal scatteringWigner coordinatesquasi-one-dimensional wiresdiscrete dopingmobility
0
0 comments X

The pith

A NEGF framework separates short-range impurity scattering from long-range Hartree terms to show that scattering rates depend on Wigner center-of-mass coordinates rather than real-space positions.

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

The paper introduces a NEGF scheme that inserts the short-range part of each impurity potential into the scattering self-energy while treating the long-range part as a self-consistent Hartree term solved with the Poisson equation. This split is used to derive an explicit expression for the position-dependent scattering rate under arbitrary impurity distributions. The resulting rate is shown to be governed by the center-of-mass coordinate in the Wigner representation, which makes the scattering intrinsically nonlocal in space. The approach is then applied to quasi-one-dimensional cylindrical wires to compute how discrete impurities modify electrostatic potential, local density of states, carrier density, scattering rates, and mobility.

Core claim

By separating the impurity potential into short-range and long-range contributions and inserting the former solely into the scattering self-energy, the position dependence of the impurity scattering rate is expressed through the center-of-mass coordinates in the Wigner representation rather than ordinary real-space coordinates. This establishes that impurity scattering is intrinsically nonlocal in space. The framework is demonstrated on cylindrical thin wires under the quasi-1D approximation, where the discrete nature of the impurities produces measurable changes in transport quantities.

What carries the argument

Separation of short-range impurity potential into the scattering self-energy (with long-range part kept as self-consistent Hartree term) together with the Wigner-coordinate representation that converts the scattering rate's position dependence into center-of-mass form.

If this is right

  • Electrostatic potential, local density of states, and carrier density inside a nanostructure become sensitive to the exact locations of individual impurities rather than only their average density.
  • Scattering rates acquire an explicit dependence on the center-of-mass coordinate, allowing position-resolved mobility calculations under inhomogeneous doping.
  • The quasi-1D wire calculations show that discrete impurities alter both the magnitude and the spatial variation of mobility compared with continuum doping models.
  • The nonlocal character of scattering implies that standard local-scattering approximations used in device simulators must be revised for structures with few impurities.

Where Pith is reading between the lines

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

  • Device models that assume local scattering rates will systematically misestimate mobility when impurity spacing becomes comparable to the coherence length.
  • The same center-of-mass representation may extend to other scattering mechanisms whose potentials can be split into short- and long-range parts, offering a route to treat surface roughness or alloy disorder in the same NEGF setting.
  • Experimental mapping of local density of states or noise spectra in nanowires could directly test the predicted nonlocality without requiring full transport measurements.

Load-bearing premise

The short-range part of the impurity potential can be cleanly separated from the long-range part and inserted solely into the scattering self-energy while the long-range part is treated exclusively as a self-consistent Hartree term.

What would settle it

A measurement or simulation in which the impurity scattering rate extracted from transport data follows ordinary real-space position dependence instead of the center-of-mass Wigner dependence predicted by the framework.

Figures

Figures reproduced from arXiv: 2501.17285 by Nobuyuki Sano.

Figure 1
Figure 1. Figure 1: FIG. 1. Quasi-1D correlation function as a function of the [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Impurity distributions (small empty circles) insid [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a)-(d) Impurity distributions inside the cylindri [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a)-(d) Impurity distributions inside the cylindri [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Contour plots of the scattering-rate matrix at [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Number of states obtained from the LDOS by inte [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Energy spectra of the current of the wires for [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
read the original abstract

A new theoretical framework for the nonequilibrium Green's function (NEGF) scheme is presented to account for the discrete nature of impurities doped in semiconductor nanostructures. The short-range part of impurity potential is included as scattering potential in the self-energy due to spatially localized impurity scattering, and the long-range part of impurity potential is treated as the self-consistent Hartree potential by coupling with the Poisson equation. The position-dependent impurity scattering rate under inhomogeneous impurity profiles is systematically derived so that its physical meaning is clarified. The position dependence of the scattering rate turns out to be represented by the `center of mass' coordinates in the Wigner coordinates, rather than the real-space coordinates. Consequently, impurity scattering is intrinsically nonlocal in space. The proposed framework is applied to cylindrical thin wires under the quasi-one-dimensional (quasi-1D) approximation. We show explicitly how the discrete nature of impurities affects the transport properties such as electrostatic potential, local density of states, carrier density, scattering rates, and mobility.

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

2 major / 1 minor

Summary. The manuscript presents a NEGF formalism for discrete impurities in semiconductor nanostructures. The short-range component of each impurity potential is routed into the scattering self-energy, while the long-range component is treated exclusively as a self-consistent Hartree term coupled to the Poisson equation. From this partition the position-dependent scattering rate is derived; its dependence appears in the center-of-mass coordinate of the Wigner representation rather than real-space coordinates, from which the authors conclude that impurity scattering is intrinsically nonlocal. The framework is applied under the quasi-1D approximation to cylindrical thin wires, with explicit results shown for electrostatic potential, local density of states, carrier density, scattering rates, and mobility.

Significance. If the central partition and the resulting coordinate dependence are rigorously justified, the work supplies a parameter-free route to incorporating discrete-impurity scattering in NEGF transport calculations for inhomogeneous doping profiles. The explicit mapping onto Wigner center-of-mass coordinates offers a concrete clarification of nonlocality that could be tested against full-potential treatments.

major comments (2)
  1. [Abstract (framework description)] Abstract (paragraph describing the framework): the clean separation of the short-range impurity potential into the scattering self-energy while routing the long-range part exclusively into the self-consistent Hartree term is asserted without derivation or numerical test of its validity for a Coulomb potential under inhomogeneous screening. Because any such partition is an approximation, the manuscript must demonstrate that cross terms or an alternative consistent treatment of the full potential leave the reported center-of-mass dependence unchanged; otherwise the claimed intrinsic nonlocality is an artifact of the split.
  2. [Abstract (derivation of scattering rate)] Abstract (derivation of scattering rate and quasi-1D application): the central claim that the scattering-rate position dependence resides in Wigner center-of-mass rather than real-space coordinates is load-bearing for the nonlocality conclusion, yet no explicit check is provided that this coordinate dependence survives when the quasi-1D approximation is relaxed or when the full (unpartitioned) impurity potential is retained in the self-energy.
minor comments (1)
  1. Ensure that all Wigner-coordinate transformations and the resulting scattering-rate expressions are written with explicit variable definitions so that the distinction between center-of-mass and relative coordinates is unambiguous to readers.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We respond to each major comment below, providing clarification on the framework and derivation while agreeing to strengthen the presentation in revision.

read point-by-point responses

  1. Referee: [Abstract (framework description)] Abstract (paragraph describing the framework): the clean separation of the short-range impurity potential into the scattering self-energy while routing the long-range part exclusively into the self-consistent Hartree term is asserted without derivation or numerical test of its validity for a Coulomb potential under inhomogeneous screening. Because any such partition is an approximation, the manuscript must demonstrate that cross terms or an alternative consistent treatment of the full potential leave the reported center-of-mass dependence unchanged; otherwise the claimed intrinsic nonlocality is an artifact of the split.

    Authors: The separation follows the standard treatment in NEGF impurity scattering, where the slowly varying long-range Coulomb tail is absorbed into the self-consistent Hartree potential via the Poisson equation, while the short-range core enters the scattering self-energy. This partition is motivated by the physical scales of screening and momentum relaxation. The center-of-mass dependence arises specifically from the Wigner transform applied to the short-range scattering self-energy for an inhomogeneous impurity distribution. Cross terms between short- and long-range components represent higher-order corrections that do not modify the leading nonlocality. We will add a dedicated paragraph in the revised manuscript justifying the approximation, citing supporting literature, and noting why it preserves the reported coordinate dependence. revision: yes

  2. Referee: [Abstract (derivation of scattering rate)] Abstract (derivation of scattering rate and quasi-1D application): the central claim that the scattering-rate position dependence resides in Wigner center-of-mass rather than real-space coordinates is load-bearing for the nonlocality conclusion, yet no explicit check is provided that this coordinate dependence survives when the quasi-1D approximation is relaxed or when the full (unpartitioned) impurity potential is retained in the self-energy.

    Authors: The derivation of the position-dependent scattering rate is carried out in the general NEGF formalism using the Wigner representation of the scattering self-energy, prior to any specialization to the quasi-1D wire geometry. The center-of-mass coordinate dependence is a direct mathematical consequence of the nonlocal scattering kernel under an inhomogeneous impurity profile and does not rely on the quasi-1D approximation. Retaining the full (unpartitioned) potential in the self-energy would double-count the long-range component already included in the Hartree term; the scattering self-energy is defined exclusively for the short-range fluctuating part. We will revise the manuscript to state explicitly that the Wigner-coordinate result is general and independent of both the quasi-1D limit and the partition details. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation follows from explicit NEGF modeling assumptions without reduction to inputs by construction

full rationale

The paper introduces an explicit modeling choice to partition the impurity potential (short-range into scattering self-energy, long-range into self-consistent Hartree/Poisson), then derives the scattering-rate position dependence in Wigner center-of-mass coordinates as a consequence of that choice within standard NEGF. This is a direct application of the formalism rather than a tautology, self-definition, fitted parameter renamed as prediction, or load-bearing self-citation. No equations are shown reducing the final nonlocality claim to the inputs by algebraic identity. The framework remains self-contained against external NEGF benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The framework rests on the standard NEGF scattering self-energy construction and the Poisson equation; the quasi-1D approximation is invoked for the wire geometry. No free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption Quasi-one-dimensional approximation is valid for cylindrical thin wires
    Invoked when applying the framework to the wire geometry.
  • domain assumption Short-range and long-range parts of the impurity potential can be treated separately in NEGF and Poisson solvers
    Central modeling choice stated in the abstract.

pith-pipeline@v0.9.0 · 5697 in / 1432 out tokens · 43481 ms · 2026-05-23T04:44:50.134344+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

45 extracted references · 45 canonical work pages

  1. [1]

    Lake and S

    R. Lake and S. Datta, Nonequilibrium green’s-function method applied to double-barrier resonant-tunneling diodes, Phys. Rev. B 45, 6670 (1992)

    work page 1992

  2. [2]

    R. Lake, G. Klimeck, R. C. Bowen, and D. Jovanovic, Sin- gle and multiband modeling of quantum electron trans- port through layered semiconductor devices, J. Appl. Phys. 81, 7845 (1997)

    work page 1997

  3. [3]

    Haug and A.-P

    H. Haug and A.-P. Jauho, Quantum Kinetics in Trans- 13 port and Optics of Semiconductors (Springer, New York, 2008)

    work page 2008

  4. [4]

    J. H. Oh, D. Ahn, Y. Yu, and S. Hwang, Green-function approach to transport through a gate-surrounded si nanowire with impurity scattering, Phys. Rev. B 77, 035313 (2008)

    work page 2008

  5. [5]

    J. H. Oh, S.-H. Lee, and M. Shin, Enhanced ionized impu- rity scattering in nanowires, J. Appl. Phys. 113, 233706 (2013)

    work page 2013

  6. [6]

    Kohn and J

    W. Kohn and J. M. Luttinger, Quantum theory of elec- trical transport phenomena, Phys. Rev. 108, 590 (1957)

    work page 1957

  7. [7]

    Abrikosov, L

    A. Abrikosov, L. Gorkov, and I. Dzyaloshinski, Methods of quantum field theory in statistical physics (Prentice- Hall, New York, 1963)

    work page 1963

  8. [8]

    Doniach and E

    S. Doniach and E. H. Sondheimer, Green ’s Functions for Solid State Physicists , Frontiers in Physics. A Lecture Note and Reprint Series (W.A. Benjamin, New York, 1974)

    work page 1974

  9. [9]

    G. D. Mahan, Many-Particle Physics, 3rd ed. (Springer, New York, 2000)

    work page 2000

  10. [10]

    M. T. Bj¨ ork, H. Schmid, J. Knoch, H. Riel, and W. Riess, Donor deactivation in silicon nanostructures, Nat. Nan- otechnol. 4, 103 (2009)

    work page 2009

  11. [11]

    Hoeneisen and C

    B. Hoeneisen and C. Mead, Fundamental limitations in microelectronics—i. mos technology, Solid-State Elec- tron. 15, 819 (1972)

    work page 1972

  12. [12]

    R. W. Keyes, The effect of randomness in the distribution of impurity atoms on fet thresholds, Appl. Phys. 8, 251 (1975)

    work page 1975

  13. [13]

    Wong and Y

    H.-S. Wong and Y. Taur, Three-dimensional ‘atomistic’ simulation of discrete random dopant distribution effects in sub-0.1 µm mosfet’s, in Technical Digest. Int. Electron Device Meeting (IEDM) (Washington DC, 1993) pp. 705– 708

    work page 1993

  14. [14]

    N. Sano, K. Matsuzawa, M. Mukai, and N. Nakayama, Role of long-range and short-range coulomb potentials in threshold characteristics under discrete dopants in sub- 0.1 µm si-mosfets, in Technical Digest. Int. Electron De- vice Meeting (IEDM) (San Francisco, 2000) pp. 275–278

    work page 2000

  15. [15]

    Miranda, When every atom counts, IEEE Spectr

    M. Miranda, When every atom counts, IEEE Spectr. 49, 32 (2012)

    work page 2012

  16. [16]

    N. Sano, K. Matsuzawa, M. Mukai, and N. Nakayama, On discrete random dopant modeling in drift-diffusion simulations: physical meaning of ‘atomistic’ dopants, Mi- croelectron. Reliab. 42, 189 (2002)

    work page 2002

  17. [17]

    N. Sano, Physical issues in device modeling: Length- scale, disorder, and phase interference, in 2017 Interna- tional Conference on Simulation of Semiconductor Pro- cesses and Devices (SISPAD) (Kamakura, Japan, 2017) pp. 1–4

    work page 2017

  18. [18]

    N. Sano, K. Yoshida, K. Tsukahara, and G. Park, Funda- mental aspects of semiconductor device modeling associ- ated with discrete impurities: Drift-diffusion simulation scheme, IEEE Trans. Electron Dev. 67, 3323 (2020)

    work page 2020

  19. [19]

    Sano and T

    N. Sano and T. Fukui, Fundamental aspects of semicon- ductor device modeling associated with discrete impuri- ties: Monte carlo simulation scheme, IEEE Trans. Elec- tron Dev. 68, 3323 (2021)

    work page 2021

  20. [20]

    Sano, Quantum kinetic equation for the wigner func- tion and reduction to the boltzmann transport equa- tion under discrete impurities, Phys

    N. Sano, Quantum kinetic equation for the wigner func- tion and reduction to the boltzmann transport equa- tion under discrete impurities, Phys. Rev. E 104, 014141 (2021)

    work page 2021

  21. [21]

    Sano, Kinetics of quasiballistic transport in nanoscale semiconductor structures: Is the ballistic limit attain- able at room temperature?, Phys

    N. Sano, Kinetics of quasiballistic transport in nanoscale semiconductor structures: Is the ballistic limit attain- able at room temperature?, Phys. Rev. Lett. 93, 246803 (2004)

    work page 2004

  22. [22]

    Martinez, M

    A. Martinez, M. Bescond, J. R. Barker, A. Svizhenko, M. P. Anantram, C. Millar, and A. Asenov, A self- consistent full 3-d real-space negf simulator for study- ing nonperturbative effects in nano-mosfets, IEEE Trans. Electron Devices 54, 2213 (2007)

    work page 2007

  23. [23]

    Bescond, C

    M. Bescond, C. Li, and M. Lannoo, Nanowire transistor modeling: influence of ionized impurity and correlation effects, J. Comput. Electron. 8, 382 (2009)

    work page 2009

  24. [24]

    N. D. Akhavan, I. Ferain, R. Yu, P. Razavi, and J.-P. Colinge, Influence of discrete dopant on quantum trans- port in silicon nanowire transistors, Solid-State Electro n. 70, 92 (2012)

    work page 2012

  25. [25]

    Carrillo-Nu˜ nez, M

    H. Carrillo-Nu˜ nez, M. Bescond, N. Cavassilas, E. Dib, and M. Lannoo, Influence of electron–phonon interac- tions in single dopant nanowire transistors, J. Appl. Phys. 116, 164505 (2014)

    work page 2014

  26. [26]

    J. Lee, P. Asenov, R. Rhyner, E. Kao, S. M. Amoroso, A. R. Brown, X.-W. Lin, and V. Moroz, Design technol- ogy co-optimization for the dram cell structure with con- tact resistance variation, IEEE Trans. Electron Devices 71, 1893 (2024)

    work page 2024

  27. [27]

    Markussen, R

    T. Markussen, R. Rurali, A.-P. Jauho, and M. Brand- byge, Scaling theory put into practice: First-principles modeling of transport in doped silicon nanowires, Phys. Rev. Lett. 99, 076803 (2007)

    work page 2007

  28. [28]

    Rurali, T

    R. Rurali, T. Markussen, J. Sun´ e, M. Brandbyge, and A.- P. Jauho, Modeling transport in ultrathin si nanowires: charged versus neutral impurities, Nano Lett. 8, 2825 (2008)

    work page 2008

  29. [29]

    L. S. J. Wilson, J. R. Barker, and A. E. Martinez, Dft/negf study of discrete dopants in si/iii–v 3d fet, J. Phys.: Condens. Matter 31, 144003 (2019)

    work page 2019

  30. [30]

    N. Sano, K. Yoshida, C.-W. Yao, and H. Watanabe, Physics of discrete impurities under the framework of de- vice simulations for nanostructure devices, Materials 11, 2559 (2018)

    work page 2018

  31. [31]

    Jacoboni and P

    C. Jacoboni and P. Lugli, The Monte Carlo Method for Semiconductor Device Simulation (Springer, New York, 1989)

    work page 1989

  32. [32]

    Vasileska, S

    D. Vasileska, S. M. Goodnick, and G. Klimeck, Compu- tational Electronics (CRC Press, New York, 2010)

    work page 2010

  33. [33]

    Kadanoff and G

    L. Kadanoff and G. Baym, Quantum Statistical Mechan- ics (W.A. Benjamin, New York, 1962)

    work page 1962

  34. [34]

    Keldysh, Diagram technique for nonequilibrium pro- cesses, JETP 20, 1018 (1965)

    L. Keldysh, Diagram technique for nonequilibrium pro- cesses, JETP 20, 1018 (1965)

    work page 1965

  35. [35]

    Caroli, R

    C. Caroli, R. Combescot, P. Nozieres, and D. Saint- James, Direct calculation of the tunneling current, J. Phys. C: Solid State Phys. 4, 916 (1971)

    work page 1971

  36. [36]

    Datta, Quantum Transport: Atom to Transistor (Cambridge University Press, Cambridge, 2005)

    S. Datta, Quantum Transport: Atom to Transistor (Cambridge University Press, Cambridge, 2005)

    work page 2005

  37. [37]

    Reif, Fundamentals of statistical and thermal physics (McGraw-Hill, Boston, 1965)

    F. Reif, Fundamentals of statistical and thermal physics (McGraw-Hill, Boston, 1965)

    work page 1965

  38. [38]

    L. E. Reichl, A Modern Course in Statistical Physics 2nd ed. (John Wiley and Sons, New York, 1998)

    work page 1998

  39. [39]

    Sano, Shallow p-n junctions under discrete impurities in semiconductor devices, IEEE Trans

    N. Sano, Shallow p-n junctions under discrete impurities in semiconductor devices, IEEE Trans. Electron Devices 71, 965 (2024)

    work page 2024

  40. [40]

    Jacoboni, Theory of Electron Transport in Semicon- ductors (Springer, New York, 2010)

    C. Jacoboni, Theory of Electron Transport in Semicon- ductors (Springer, New York, 2010)

    work page 2010

  41. [41]

    M. V. Fischetti and W. G. Vandenberghe, Advanced Physics of Electron Transport in Semiconductors and Nanostructures (Springer, New York, 2016). 14

    work page 2016

  42. [42]

    Sano, Variability and self-average of impurity-limited resistance in quasi-one dimensional nanowires, Solid- State Electronics 128, 25 (2017)

    N. Sano, Variability and self-average of impurity-limited resistance in quasi-one dimensional nanowires, Solid- State Electronics 128, 25 (2017)

    work page 2017

  43. [43]

    W. R. Frensley, Boundary conditions for open quantum systems driven far from equilibrium, Rev. Mod. Phys. 62, 745 (1990)

    work page 1990

  44. [44]

    Golizadeh-Mojarad and S

    R. Golizadeh-Mojarad and S. Datta, Nonequilibrium green’s function based models for dephasing in quantum transport, Phys. Rev. B 75, 081301 (2007)

    work page 2007

  45. [45]

    Sano, Fundamental aspects of semiconductor device modeling associated with discrete impurities: Nonequi- librium green’s function scheme, IEEE Trans

    N. Sano, Fundamental aspects of semiconductor device modeling associated with discrete impurities: Nonequi- librium green’s function scheme, IEEE Trans. Electron Devices 72, 24 (2025)

    work page 2025