A large-scale first-principles quantum transport simulation method using plane waves

Lin-Wang Wang; Meng Ye; Shu-Shen Li; Xiangwei Jiang

arxiv: 1907.06821 · v1 · pith:NCFH2ONUnew · submitted 2019-07-16 · ❄️ cond-mat.mes-hall · physics.comp-ph

A large-scale first-principles quantum transport simulation method using plane waves

Meng Ye , Xiangwei Jiang , Shu-Shen Li , Lin-Wang Wang This is my paper

Pith reviewed 2026-05-24 21:12 UTC · model grok-4.3

classification ❄️ cond-mat.mes-hall physics.comp-ph

keywords quantum transportplane wave basisNEGFDFTnanowireslarge scale simulationfirst principlescopper atoms

0 comments

The pith

Plane waves with parallel algorithms enable quantum transport simulations of several thousand atoms.

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

This paper presents a method for calculating quantum transport from first principles using a plane-wave basis set. It builds on an earlier plane-wave NEGF approach by incorporating efficient parallel algorithms including linear-scaling DFT, the folded spectrum method, and filtering techniques. These additions allow simulations of systems with several thousand atoms, far beyond the typical limit of under 1000 atoms for atomic-orbital based methods. Such capability is needed to model realistic nanoscale electronic devices where quantum effects are important.

Core claim

By extending the plane-wave NEGF formulation with high-efficiency parallel algorithms such as linear-scale ground-state DFT, folded spectrum method, and filtering technique, the method can perform quantum transport simulations on systems containing several thousands of atoms. This is demonstrated through calculations on copper nanowires with about 4000 atoms, revealing how their shape and point defects influence transport properties.

What carries the argument

The plane-wave based non-equilibrium Green's function (NEGF) method augmented with linear-scale DFT and filtering techniques for large systems.

Load-bearing premise

The 2005 plane-wave NEGF formulation can be extended with the listed parallel algorithms while preserving numerical accuracy and stability for transport quantities in systems of several thousand atoms.

What would settle it

A direct comparison of calculated conductance or transmission spectra for a 4000-atom copper nanowire against results from a subdivided smaller system or an independent atomic-basis calculation would reveal any loss of accuracy or stability.

read the original abstract

As the characteristic lengths of advanced electronic devices are approaching the atomic scale, ab initio simulation method, with fully consideration of quantum mechanical effects, becomes essential to study the quantum transport phenomenon in them. However, current widely used non-equilibrium Green's function (NEGF) approach is based on atomic basis set, which usually can only study small system with less than 1000 atoms in practice. Here we present a large-scale quantum transport simulation method using plane waves basis, based on the previously developed plane wave approach (Phys. Rev. B 72, 045417). By applying several high-efficiency parallel algorithms, such as linear-scale ground-state density function theory (DFT) algorithm, folded spectrum method, and filtering technique, we demonstrate that our new method can simulate the system with several thousands of atoms. We also use this method to study several nanowires with about 4000 copper atoms, and show how the shape and point defect affect the transport properties of them. Such quantum simulation method will be useful to investigate and design nanoscale devices, especially the on-die interconnects.

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.

Desk Editor's Note private letter to a colleague

This extends the 2005 plane-wave NEGF method with linear-scaling DFT, folded spectrum, and filtering to reach ~4000-atom Cu nanowires, but the abstract supplies no accuracy benchmarks or cross-checks.

read the letter

The core of the paper is an engineering push on the 2005 plane-wave NEGF framework. By layering in linear-scaling ground-state DFT, folded-spectrum eigenvalue work, and filtering, they report transport calculations on copper nanowires with several thousand atoms and then examine how wire shape and point defects change the conductance. That is the concrete step forward: moving from sub-1000 atom limits to device-relevant sizes without switching to atomic-orbital bases. The examples on interconnect-like nanowires give a sense of what the method can actually produce once the scaling is in place. The citation to the earlier Phys. Rev. B paper keeps the theoretical foundation explicit rather than reinvented. The parallel-algorithm choices are standard tools applied to the transport problem, which is reasonable engineering. The main gap is the absence of any reported validation numbers. The abstract states the large-system runs succeeded but gives no comparison to smaller-basis results, no transmission convergence with plane-wave cutoff, and no error estimates on the Landauer conductance. Without those, it is hard to judge whether the non-Hermitian Green's-function and self-energy steps remain stable when the Hermitian-focused accelerators are inserted. If the full manuscript contains small-system cross-checks and timing data against the original formulation, the claim is on firmer ground; if not, that is the section that needs strengthening. This work is aimed at groups doing ab initio transport on nanowires or interconnects who already use plane-wave codes and need to reach realistic atom counts. A reader who cares about practical device modeling would find the nanowire results and the scaling recipe useful. I would send the paper to referees. The scaling demonstration is worth external scrutiny even if the validation section requires expansion.

Referee Report

2 major / 1 minor

Summary. The manuscript presents a plane-wave basis implementation of the non-equilibrium Green's function (NEGF) formalism for quantum transport, extending the 2005 formulation (Phys. Rev. B 72, 045417). It incorporates linear-scaling ground-state DFT, the folded spectrum method, and filtering techniques to achieve simulations of systems containing several thousand atoms, with explicit demonstrations on copper nanowires of approximately 4000 atoms that examine the influence of nanowire shape and point defects on transport properties.

Significance. If the accuracy and stability of the Green's functions and transmission functions are preserved, the approach would enable ab initio quantum transport calculations at scales previously inaccessible to atomic-orbital NEGF codes (typically limited to <1000 atoms), providing a useful tool for studying realistic nanoscale interconnects and devices. The engineering achievement of scaling to 4000 atoms via the cited parallel algorithms is noteworthy if accompanied by validation.

major comments (2)

[Abstract] Abstract: The central claim that the method successfully simulates ~4000-atom Cu nanowires is unsupported by any reported validation metrics, direct comparisons against the original 2005 plane-wave NEGF formulation on smaller benchmark systems, error estimates, or convergence data for transport quantities such as conductance or transmission.
[Methods] Methods/Results: The adaptation of the folded spectrum method and filtering technique (standard for Hermitian DFT eigenvalue problems) to the complex, energy-dependent, non-Hermitian matrix inversions and self-energy calculations required by NEGF is presented without explicit numerical tests confirming that accuracy and stability are retained for Landauer transport observables at the 4000-atom scale.

minor comments (1)

[Abstract] The phrase 'density function theory' in the abstract should read 'density functional theory'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting areas where additional validation would strengthen the presentation. We address the two major comments point by point below and commit to incorporating the requested material in a revised version.

read point-by-point responses

Referee: [Abstract] Abstract: The central claim that the method successfully simulates ~4000-atom Cu nanowires is unsupported by any reported validation metrics, direct comparisons against the original 2005 plane-wave NEGF formulation on smaller benchmark systems, error estimates, or convergence data for transport quantities such as conductance or transmission.

Authors: We agree that the abstract would benefit from supporting validation details. The manuscript demonstrates transport calculations on ~4000-atom Cu nanowires and reports the effects of shape and point defects, but does not include the explicit benchmarks, comparisons to the 2005 formulation, error estimates, or convergence data mentioned. In the revision we will add a new subsection (likely in the Results section) that provides (i) direct comparisons of transmission and conductance against the original plane-wave NEGF implementation on smaller benchmark systems (up to a few hundred atoms), (ii) error estimates relative to those references, and (iii) convergence tests with respect to energy grid density and basis cutoff for the large-scale cases. These additions will directly support the scalability claim. revision: yes
Referee: [Methods] Methods/Results: The adaptation of the folded spectrum method and filtering technique (standard for Hermitian DFT eigenvalue problems) to the complex, energy-dependent, non-Hermitian matrix inversions and self-energy calculations required by NEGF is presented without explicit numerical tests confirming that accuracy and stability are retained for Landauer transport observables at the 4000-atom scale.

Authors: The folded spectrum and filtering approaches are applied within the energy window relevant to the Landauer transport calculation, with the non-Hermitian character handled by the standard NEGF self-energy and Green's function inversion steps. While the manuscript describes the algorithmic extensions, we acknowledge that explicit numerical verification of accuracy and stability for the non-Hermitian case at large scale is not provided. In the revised manuscript we will insert additional tests on smaller systems that compare the adapted folded-spectrum/filtered results against direct inversion, quantify the deviation in transmission and conductance, and report residual norms or similar stability indicators for the 4000-atom calculations to confirm that the adaptation preserves the required accuracy. revision: yes

Circularity Check

0 steps flagged

Extension of 2005 plane-wave NEGF via parallel algorithms exhibits no circular derivation

full rationale

The manuscript extends the cited Phys. Rev. B 72, 045417 plane-wave NEGF formulation by incorporating linear-scale DFT, folded spectrum, and filtering techniques, then demonstrates scaling to ~4000-atom Cu nanowires. No equations, definitions, or claims reduce the claimed capability to a self-referential fit, renamed prediction, or load-bearing self-citation whose validity depends on the present work. The scaling result is presented as an empirical engineering outcome built on the external base method, satisfying the criterion for a self-contained implementation paper.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The method inherits standard DFT and NEGF formalisms from prior literature; no new free parameters, ad-hoc axioms, or invented entities are introduced in the abstract.

axioms (1)

domain assumption Standard Kohn-Sham DFT with appropriate exchange-correlation functional yields sufficiently accurate ground-state densities for the transport calculation.
The linear-scale DFT step is presented as the foundation without further justification in the abstract.

pith-pipeline@v0.9.0 · 5723 in / 1173 out tokens · 17389 ms · 2026-05-24T21:12:40.106708+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean (washburn_uniqueness_aczel, Jcost definition) reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

By applying several high-efficiency parallel algorithms, such as linear-scale ground-state density function theory (DFT) algorithm, folded spectrum method, and filtering technique, we demonstrate that our new method can simulate the system with several thousands of atoms.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.