Band-structure and electronic transport calculations in cylindrical wires : the issue of bound states in transfer-matrix calculations
Pith reviewed 2026-05-24 20:58 UTC · model grok-4.3
The pith
Bound states must be included in transfer-matrix calculations for cylindrical wires to capture resonances and continuous band structures.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When the transfer-matrix methodology is applied to electrons subject to cosine potentials in a cylindrical wire, bound states must be considered because of their impact as sharp resonances in the transmission probabilities and to remove unphysical discontinuities in the band structure.
What carries the argument
Transfer matrices for a single basic unit combined with the layer-addition algorithm to extract band structure from periodic repetitions in three dimensions.
If this is right
- Transmission probabilities exhibit sharp resonances attributable to bound states.
- Extracted band structures become continuous without artificial discontinuities.
- The representation of scattering states in the continuous spectrum gains improved completeness.
Where Pith is reading between the lines
- The same layer-addition procedure could be applied directly to other wire cross-sections or potential forms to test consistency of resonance positions.
- Comparison of the resulting transmission curves against experimental conductance measurements in fabricated nanowires would provide an external check on the bound-state contribution.
- The method's handling of bound states may carry over to related wave equations such as those for electromagnetic modes in dielectric cylinders.
Load-bearing premise
The transfer-matrix method generalized from one to three dimensions and combined with the layer-addition algorithm remains numerically stable and complete for the continuous spectrum without additional validation against independent benchmarks.
What would settle it
An independent numerical solution of the same Schrödinger problem in the cylindrical wire, for example by direct diagonalization or finite-element methods, that yields transmission probabilities without sharp resonances or band structures without discontinuities when bound states are omitted.
read the original abstract
The transfer-matrix methodology is used to solve linear systems of differential equations, such as those that arise when solving Schr\"odinger's equation, in situations where the solutions of interest are in the continuous part of the energy spectrum. The technique is actually a generalization in three dimensions of methods used to obtain scattering solutions in one dimension. Using the layer-addition algorithm allows one to control the stability of the computation and to describe efficiently periodic repetitions of a basic unit. This paper, which is an update of an article originally published in Physical and Chemical News 16, 46-53 (2004), provides a pedagogical presentation of this technique. It describes in details how the band structure associated with an infinite periodic medium can be extracted from the transfer matrices that characterize a single basic unit. The method is applied to the calculation of the transmission and band structure of electrons subject to cosine potentials in a cylindrical wire. The simulations show that bound states must be considered because of their impact as sharp resonances in the transmission probabilities and to remove unphysical discontinuities in the band structure. Additional states only improve the completeness of the representation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a pedagogical description of a three-dimensional generalization of the transfer-matrix method combined with a layer-addition algorithm for solving Schrödinger's equation in the continuous spectrum. It applies the technique to electrons in cylindrical wires subject to periodic cosine potentials, extracting band structures from transfer matrices of a single unit and computing transmission probabilities. The central claim is that bound states must be retained in the basis, as their omission produces sharp unphysical resonances in transmission and discontinuities in the band structure; additional states improve completeness.
Significance. If the numerical results hold under independent verification, the work would usefully illustrate a completeness requirement for transfer-matrix calculations in quasi-1D mesoscopic systems, aiding accurate modeling of transport and periodic structures. The layer-addition stability control and periodic-unit extraction are potentially reusable strengths, though the manuscript supplies no machine-checked proofs or parameter-free derivations.
major comments (1)
- [Abstract and Simulations] Abstract and results sections: the assertion that 'simulations show that bound states must be considered' because of resonances and discontinuities rests on unshown numerical evidence. No quantitative transmission curves, band-structure plots with/without bound states, error analysis, convergence tests, or comparisons to independent discretizations (radial finite-difference, DVR, or plane-wave expansion) are supplied, leaving open whether observed artifacts arise from missing bound states or from truncation/evanescent-mode handling in the 3-D layer-addition step itself.
Simulated Author's Rebuttal
We thank the referee for the careful reading of the manuscript and the constructive feedback. We address the single major comment below and will revise the manuscript to strengthen the presentation of the numerical evidence.
read point-by-point responses
-
Referee: [Abstract and Simulations] Abstract and results sections: the assertion that 'simulations show that bound states must be considered' because of resonances and discontinuities rests on unshown numerical evidence. No quantitative transmission curves, band-structure plots with/without bound states, error analysis, convergence tests, or comparisons to independent discretizations (radial finite-difference, DVR, or plane-wave expansion) are supplied, leaving open whether observed artifacts arise from missing bound states or from truncation/evanescent-mode handling in the 3-D layer-addition step itself.
Authors: We agree that the manuscript would be improved by more explicit display of the supporting numerical results. The simulations underlying the abstract statement were performed and did exhibit the reported sharp transmission resonances and band-structure discontinuities when bound states were omitted from the basis; these effects were the motivation for the claim. To address the concern, the revised version will include new figures presenting quantitative transmission probabilities and band structures computed both with and without bound states, together with a brief discussion of convergence with respect to basis size and a clarification of how evanescent modes are handled within the layer-addition algorithm. This will make the source of the artifacts clearer. Direct side-by-side comparisons with independent methods such as radial finite-difference or DVR are outside the pedagogical scope of the present work but could be noted as a possible extension. revision: yes
Circularity Check
No circularity; method is direct numerical procedure
full rationale
The paper presents the transfer-matrix method with layer-addition as a generalization of 1D scattering techniques to 3D cylindrical wires, applied directly to cosine potentials to compute transmission and band structure. Conclusions about bound states arise from observed simulation outputs (resonances and discontinuities), not from any equation that reduces by construction to fitted inputs, self-definitions, or unverified self-citations. No load-bearing uniqueness theorems, ansatzes smuggled via prior work, or renaming of known results are described. The derivation chain is self-contained as an explicit computational algorithm without statistical forcing or definitional loops.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Solutions of interest lie in the continuous part of the energy spectrum.
- domain assumption The layer-addition algorithm controls stability for periodic repetitions.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.