pith. sign in

arxiv: 2605.08032 · v1 · submitted 2026-05-08 · 🪐 quant-ph · cond-mat.str-el

Analytical Solution to the Kronig-Penney Model with Harmonic Oscillator Wells: Insights to Tight-Binding

Christopher Moore , Frank Marsiglio This is my paper

Pith reviewed 2026-05-11 02:35 UTC · model grok-4.3

classification 🪐 quant-ph cond-mat.str-el
keywords Kronig-Penney modelharmonic oscillator potentialtight-binding approximationtunneling amplitudeenergy dispersionBloch theoremperiodic potential
0
0 comments X

The pith

Kronig-Penney model with truncated harmonic oscillator wells admits an exact analytical solution that yields the tight-binding tunneling amplitude in terms of oscillator parameters.

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

The paper replaces the square well potentials of the standard Kronig-Penney model with truncated harmonic oscillator wells. It solves the periodic Schrödinger equation analytically by applying boundary matching and the Bloch theorem to obtain explicit energy dispersion and wave functions. The resulting equation is rewritten in a form that isolates the nearest-neighbor hopping term, allowing the tunneling amplitude to be written directly as a function of the oscillator frequency and truncation radius. A sympathetic reader cares because this step connects the shape of a smoother atomic potential to the simple parameters used in effective lattice models for solids.

Core claim

By using truncated harmonic oscillator potentials in place of square wells for the atomic sites, the model can be solved exactly using boundary condition matching and the Bloch theorem. The resulting dispersion relation and eigenfunctions provide an exact expression for the inter-site tunneling amplitude in terms of the oscillator strength and truncation radius.

What carries the argument

The truncated harmonic oscillator potential arranged in a one-dimensional periodic lattice, solved by matching solutions at the truncation points and imposing Bloch periodicity.

If this is right

  • The energy dispersion and wave functions are obtained in closed form without numerical root finding.
  • The governing equation takes a form that directly motivates the tight-binding approximation.
  • The tunneling amplitude is expressed explicitly in terms of the harmonic oscillator frequency and truncation radius.
  • Similarities and differences with the square-well Kronig-Penney model are identified through direct comparison of the derived quantities.

Where Pith is reading between the lines

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

  • This derivation offers a route to obtain tight-binding parameters from smoother, more realistic potential shapes without separate numerical fitting.
  • The method could be tested on finite chains to check how well the infinite-periodic analytical result approximates real finite systems.
  • Extensions to two or three dimensions would connect the approach to band structures in actual crystals.

Load-bearing premise

The truncated harmonic oscillator wells permit an exact analytical solution through boundary matching in the same way as square wells.

What would settle it

A numerical integration of the Schrödinger equation for the identical periodic potential that produces a different energy dispersion would show the analytical solution is incorrect.

Figures

Figures reproduced from arXiv: 2605.08032 by Christopher Moore, Frank Marsiglio.

Figure 1
Figure 1. Figure 1: FIG. 1: Diagram of a one-dimensional periodic potential [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Plot of the lowest four energy bands as a function of [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Plot of the real and imaginary parts of the wavefunction as a function of position across three cells for various wave [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Comparison of the first-order tight-binding dispersion (Equation 16) with the analytic dispersion for the ground state [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Comparison of the first-order tight-binding dispersion (Equation 16) with the analytic dispersion for the ground state [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Expanded plot of the first three energy bands in the [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Plot of the numerical ground state wave function for [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
read the original abstract

The celebrated Kronig-Penney model traditionally has been formulated with square well potentials representing atomic centres. Here, we use a slightly more realistic potential, the truncated harmonic oscillator, in lieu of square well potentials, and solve the model analytically. We derive the energy dispersion and wave functions for this model. This configuration has some important similarities and differences compared to the usual model. In particular, we write the governing equation in a form suggestive of the tight-binding approximation, as can be done for the usual model. In this way, it is straightforward to derive an expression for the tunneling amplitude used in tight-binding in terms of the harmonic oscillator potential parameters.

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

1 major / 2 minor

Summary. The manuscript formulates a Kronig-Penney model with periodic truncated harmonic-oscillator wells in place of square wells. It solves the Schrödinger equation inside each well with parabolic-cylinder functions, imposes continuity of the wave function and its derivative at the truncation points, applies Bloch periodicity, and obtains the energy dispersion E(k) together with the corresponding wave functions. From this setup the authors extract an explicit expression for the tight-binding tunneling amplitude written directly in terms of the oscillator frequency, depth, and truncation radius.

Significance. If the matching conditions and the resulting dispersion relation are correctly derived, the work supplies a concrete, parameter-free link between a microscopically realistic potential and the hopping term that appears in tight-binding models. This is a useful theoretical bridge that is absent from the conventional square-well Kronig-Penney treatment and could be used to test the range of validity of the tight-binding approximation without empirical fitting.

major comments (1)
  1. [Abstract and dispersion derivation] Abstract and the section deriving the dispersion relation: the assertion that the model is solved 'analytically' and that the tunneling amplitude follows 'straightforwardly' must be qualified. The interior solutions are parabolic-cylinder functions (special functions, not elementary), and the Bloch-periodicity condition produces a transcendental equation in E and k that must be solved numerically for each k, exactly as in the standard Kronig-Penney model. The manuscript should state this explicitly and indicate whether the tunneling-amplitude expression is obtained by an additional approximation or is itself an implicit function of the roots of that equation.
minor comments (2)
  1. [Model definition] The boundary conditions at the truncation points (continuity of ψ and ψ′ with the exterior solution) should be written out explicitly, including the form assumed for the exterior wave function.
  2. [Potential and wave-function sections] Notation for the truncation radius and the matching points should be introduced once and used consistently; several symbols appear to be redefined between the potential sketch and the matching equations.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive comment. We agree that the claims of an 'analytical' solution and 'straightforward' derivation of the tunneling amplitude require qualification, and we will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Abstract and dispersion derivation] Abstract and the section deriving the dispersion relation: the assertion that the model is solved 'analytically' and that the tunneling amplitude follows 'straightforwardly' must be qualified. The interior solutions are parabolic-cylinder functions (special functions, not elementary), and the Bloch-periodicity condition produces a transcendental equation in E and k that must be solved numerically for each k, exactly as in the standard Kronig-Penney model. The manuscript should state this explicitly and indicate whether the tunneling-amplitude expression is obtained by an additional approximation or is itself an implicit function of the roots of that equation.

    Authors: We agree that the solution is expressed using parabolic-cylinder functions, which are special functions, and that the dispersion relation takes the form of a transcendental equation in E and k that must be solved numerically, precisely as in the conventional Kronig-Penney model. We will revise the abstract and the relevant derivation section to state this explicitly. The tunneling-amplitude expression is obtained directly from the exact wave-function matching conditions and the rewriting of the dispersion relation into a tight-binding-like form; it is an explicit closed-form expression in terms of the oscillator frequency, depth, and truncation radius, without requiring numerical roots of the transcendental equation or further approximations. revision: yes

Circularity Check

0 steps flagged

No circularity; derivation proceeds directly from Schrödinger equation and boundary matching.

full rationale

The paper solves the time-independent Schrödinger equation for truncated harmonic-oscillator wells using parabolic-cylinder functions, imposes continuity of ψ and ψ′ at truncation points, and applies Bloch periodicity to obtain a transcendental condition on E(k). From this it extracts an explicit expression for the tight-binding tunneling amplitude written in terms of the input well parameters. No equation is defined in terms of its own output, no fitted quantity is relabeled as a prediction, and no load-bearing step rests on a self-citation. The entire chain is therefore self-contained against the stated potential and standard quantum-mechanical techniques.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The claim relies on standard quantum mechanics for periodic systems and the assumption that the chosen potential form permits analytical treatment. No new entities are introduced. Free parameters are the physical parameters of the oscillator wells.

free parameters (1)
  • harmonic oscillator parameters (frequency, depth, truncation point)
    These are input parameters of the potential, not fitted to data in the abstract.
axioms (2)
  • standard math Bloch's theorem applies to the periodic potential allowing plane-wave like solutions modulated by periodic function.
    Standard in solid state physics for periodic potentials.
  • domain assumption The potential is truncated harmonic oscillator allowing matching of wave functions at boundaries for analytical solution.
    Assumed to enable closed form solution.

pith-pipeline@v0.9.0 · 5406 in / 1564 out tokens · 59719 ms · 2026-05-11T02:35:12.048495+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

18 extracted references · 18 canonical work pages

  1. [1]

    Bloch, Zeitschrift f¨ ur Physik52, 555 (1929), ISSN 0044-3328

    F. Bloch, Zeitschrift f¨ ur Physik52, 555 (1929), ISSN 0044-3328

    work page 1929

  2. [2]

    Bloch, Proceedings of the Royal Soci- ety of London

    F. Bloch, Proceedings of the Royal Soci- ety of London. A. Mathematical and Physi- cal Sciences371, 24 (1980), ISSN 0080-4630, https://royalsocietypublishing.org/rspa/article- pdf/371/1744/24/63996/rspa.1980.0052.pdf, URL https://doi.org/10.1098/rspa.1980.0052

    work page doi:10.1098/rspa.1980.0052 1980

  3. [3]

    Brillouin, Journal de Physique et le Radium1, 377 (1930), ISSN 0044-3328

    L. Brillouin, Journal de Physique et le Radium1, 377 (1930), ISSN 0044-3328

    work page 1930

  4. [4]

    R. D. L. Kronig and W. G. Penney, Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character130, 499 (1931)

    work page 1931

  5. [5]

    Wigner and F

    E. Wigner and F. Seitz, Phys. Rev.43, 804 (1933), URL https://link.aps.org/doi/10.1103/PhysRev.43.804

    work page doi:10.1103/physrev.43.804 1933

  6. [6]

    Wigner and F

    E. Wigner and F. Seitz, Phys. Rev.46, 509 (1934), URL https://link.aps.org/doi/10.1103/PhysRev.46.509

    work page doi:10.1103/physrev.46.509 1934

  7. [7]

    R. M. Martin, Electronic Structure: Basic Theory and Practical Methods (Cambridge University Press, 2020), 2nd ed

    work page 2020

  8. [8]

    I. D. Johnston and D. Segal, American Journal of Physics 60, 600 (1992)

    work page 1992

  9. [9]

    R. L. Pavelich and F. Marsiglio, American Journal of Physics83, 773 (2015), ISSN 0002-9505

    work page 2015

  10. [10]

    R. L. Pavelich and F. Marsiglio, American Journal of Physics84, 924 (2016)

    work page 2016

  11. [11]

    D. C. Johnston, 1905.12084, URLhttps://doi.org/10. 48550/arXiv.1905.12084

    work page arXiv 1905

  12. [12]

    D. C. Johnston, American Journal of Physics88, 1109 (2020), ISSN 0002-9505, URLhttps://doi.org/10. 1119/10.0001863

    work page 2020

  13. [13]

    Marsiglio and R

    F. Marsiglio and R. L. Pavelich, Scientific Reports7, 17041 (2017), ISSN 2045-2322

    work page 2017

  14. [14]

    Forcade, D

    G. Forcade, D. van Paridon, R. L. Pavelich, and F. Mar- siglio, European Journal of Physics42, 025408 (2021)

    work page 2021

  15. [15]

    Abramowitz and I

    M. Abramowitz and I. Stegun, 0.4 0.2 0.0 0.2 0.4 x/ 0.00 0.25 0.50 0.75 1.00 1.25 1.50 1.75 2.00(x) Oscillator Ground State Inf. Well Periodic BCs FIG. 8: Plot of the numerical ground state wave function for the harmonic oscillator well potential (see Fig. 1) with Dirich- let and periodic boundary conditions forV 0/(ℏω/2) = 6 and w/ℓ= 2/3. For comparison ...

    work page 1970

  16. [16]

    DLMF, NIST Digital Library of Mathematical Functions 15 (2020), f. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, and M. A. McClain, eds., URLhttp://dlmf.nist.gov/13.2

    work page 2020

  17. [17]

    Ashcroft and N

    N. Ashcroft and N. Mermin, Solid State Physics (Holt, Rinehart and Winston, 1976), ISBN 9780030493461

    work page 1976

  18. [18]

    L. U. Ancarani and G. Gasaneo, Journal of Physics B: Atomic, Molecular and Optical Physics41, 105001 (2008), URLhttps://doi.org/10.1088/0953-4075/41/ 10/105001

    work page doi:10.1088/0953-4075/41/ 2008