Sluggish quantum mechanics of noninteracting fermions with spatially varying effective mass
Pith reviewed 2026-05-10 15:32 UTC · model grok-4.3
The pith
Noninteracting fermions with position-dependent effective mass have a correlation kernel near the origin that is a sum of two Bessel kernels.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For alpha greater than zero the scaled kernel of the determinantal process near x equals zero is neither the standard Bessel nor Airy kernel but equals the sum of two Bessel kernels with distinct indices; this implies that the large-N average density has a non-monotonic profile with a vanishing minimum at the origin while remaining of finite width.
What carries the argument
the scaled correlation kernel near the origin expressed as the sum of two Bessel kernels with different indices, which determines all position statistics of the ground-state fermions
If this is right
- The average density for large particle number is non-monotonic, symmetric about the origin, and exactly zero at the center.
- All higher-order correlation functions of the fermion positions follow from the determinantal kernel for any finite N.
- Exact eigenfunctions and the time-dependent propagator are available for the single-particle problem both with and without the external trap.
- The construction directly models optical lattices in which tunneling strength varies with site position.
Where Pith is reading between the lines
- The same power-law mass variation may produce analogous new kernels in higher-dimensional or interacting settings if the continuum limit remains controllable.
- Direct comparison between the tight-binding lattice model and the continuum predictions for moderate N would test the accuracy of the BenDaniel-Duke approximation.
- The vanishing central density offers a clear signature that could be sought in cold-atom experiments with engineered position-dependent tunneling.
Load-bearing premise
The continuum limit of the inhomogeneous tight-binding model with power-law hopping produces the BenDaniel-Duke equation with m_eff proportional to |x|^alpha, and the chosen trap potential permits exact closed-form solutions.
What would settle it
Numerical or experimental measurement of the ground-state density for large N in an optical lattice whose hopping amplitudes realize m_eff(x) proportional to |x|^alpha, checking whether the density reaches a minimum of zero at the center while staying finite and symmetric elsewhere.
Figures
read the original abstract
We analyze a class of one-dimensional quantum systems characterized by a position-dependent kinetic term arising as the continuum limit of an inhomogeneous tight-binding model with spatially varying hopping amplitudes. In this limit, the Schrodinger equation takes the so-called BenDaniel-Duke form with an effective mass, scaling as $m_{eff}(x) = m_{eff}|x|^{\alpha}$ with $\alpha > 0$, leading to a framework we term sluggish quantum mechanics, where particle motion is progressively suppressed at larger distances. Both without any external potential and with $V_{ext}(x)=\frac{1}{2}m_{eff}\omega^2 |x|^{\alpha+2}$, we obtain the eigenfunctions and the quantum propagators exactly. We then investigate the problem of $N$ noninteracting spinless fermions in the trap, determining the many-body ground-state wavefunction and the joint probability density function of the positions of the $N$ fermions. We show that the many-body quantum probability density in the ground state forms a determinantal point process whose correlation kernel can be computed for any $N$, giving access to the average density as well as higher order correlation functions for any finite $N$. Moreover, we analyze the scaling form of this kernel in the large $N$ limit in the bulk, near the edge, and close to the origin. Our results show that the scaled average density profile for large $N$ has a finite support symmetric with respect to the origin, but has a non-monotonic shape with a vanishing minimum at the origin for any $\alpha>0$. One of the key findings of our work is that the scaled kernel near the origin $x=0$ for $\alpha>0$ is neither the Bessel nor the Airy kernel (that are standard for trapped fermions), but is new, and is given by a sum of two Bessel kernels with different indices. Our results thus provide a framework relevant to engineered optical lattices with position-dependent tunneling.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes one-dimensional noninteracting fermions with position-dependent effective mass m_eff(x) = m_eff |x|^α (α > 0) arising as the continuum limit of an inhomogeneous tight-binding model. Exact eigenfunctions and propagators are derived for the BenDaniel-Duke Schrödinger equation in both the free case and with the trap V_ext(x) = (1/2) m_eff ω² |x|^{α+2}. For N fermions in the trap, the ground-state many-body wavefunction yields a determinantal point process whose correlation kernel is obtained exactly via the Christoffel-Darboux identity; large-N scaling limits of the kernel are analyzed in the bulk, at the edge, and near the origin, where the scaled kernel is shown to be a sum of two Bessel kernels with different indices rather than the standard Airy or Bessel kernels.
Significance. If the derivations hold, this provides a rare exactly solvable inhomogeneous quantum system with closed-form eigenfunctions, finite-N kernels, and explicit scaling limits, including a novel determinantal kernel at the origin. The exact solvability and parameter-free asymptotic analysis constitute a clear strength, with direct relevance to engineered optical lattices featuring position-dependent tunneling amplitudes.
minor comments (3)
- [§2] §2: The continuum expansion from the discrete tight-binding Hamiltonian to the BenDaniel-Duke operator is sketched but would benefit from an explicit display of the second-order terms in the hopping expansion to confirm the precise form of m_eff(x).
- [§4.2] §4.2: The statement that the origin kernel is 'new' is supported by the asymptotic analysis, but a short comparison table or explicit functional form contrasting it with the standard Bessel kernel (e.g., via the indices ν) would aid readability.
- [Figure 3] Figure 3 caption: The scaling variable for the origin kernel is not defined in the caption; it should be stated explicitly (e.g., ξ = N^β x) to match the text in §5.
Simulated Author's Rebuttal
We thank the referee for the careful reading, accurate summary, and positive assessment of our work. We are pleased that the exact solvability, closed-form kernels, and the novel sum-of-Bessel structure at the origin are recognized as strengths with relevance to engineered optical lattices. No specific major comments were listed in the report.
Circularity Check
No significant circularity; derivation is self-contained from continuum limit and exact solutions
full rationale
The paper begins with the standard continuum expansion of the inhomogeneous tight-binding Hamiltonian to obtain the BenDaniel-Duke operator with m_eff(x) = m_eff |x|^α. The external trap is deliberately chosen as V_ext(x) = (1/2) m_eff ω² |x|^{α+2} so that a change of variables reduces the eigenvalue problem exactly to the Bessel equation, yielding closed-form eigenfunctions |x|^β J_ν(λ |x|^γ) and Y_ν. The N-fermion ground-state kernel is then constructed via the Christoffel-Darboux identity applied to these explicit functions. Large-N scaling limits (bulk, edge, origin) follow from standard asymptotic expansions of the Bessel functions in microscopic variables. No step reduces by construction to a fitted parameter, self-citation chain, or redefinition of the target result; the central claim (new origin kernel as sum of two Bessel kernels) is obtained directly from the asymptotics without circular reduction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The continuum limit of an inhomogeneous tight-binding model with spatially varying hopping amplitudes yields the BenDaniel-Duke form of the Schrödinger equation with m_eff(x) = m_eff |x|^α.
- ad hoc to paper The external potential is chosen as V_ext(x) = 1/2 m_eff ω² |x|^{α+2} to permit exact solvability.
Reference graph
Works this paper leans on
-
[1]
,so thatΓ(a)diverges, and consequently, the sec- ond term in (53) vanishes
Even eigenfunctions To select the family of eigenfunctions belonging toφ=0, i.e., the ‘even’ sector, we need to seta=−nwithn= 0,1,2, . . . ,so thatΓ(a)diverges, and consequently, the sec- ond term in (53) vanishes. Using the expression ofafrom Eq. (51), gives the quantized eigenvalues corresponding to the even eigenfunctions as λ even n =2(α+2)µn+µ,wheren...
-
[2]
, so thatΓ(1+a−b)diverges, and consequently, the first term in (53) vanishes
Odd eigenfunction To find the family of eigenfunctions belonging to theφ= 1+α, i.e., the ‘odd’ sector, we need to set 1+a−b=−n, with n=0,1,2, . . ., so thatΓ(1+a−b)diverges, and consequently, the first term in (53) vanishes. Using the expression ofaandb from Eq. (51), gives the quantized eigenvalues corresponding to the odd eigenfunctions as λ odd n =2(α+...
-
[3]
Bulk kernel We first consider the case when bothxandyare in the bulk of the density profile (see Fig. 4). Since the density profile is symmetric around the origin, for simplicity, we consider the positive side, i.e.,x>0 andy>0. In the bulk, the typical interparticle separationℓ(x)∼[Nρ N(x)]−1. On the other hand, α 0 α 0 α 0.1 α 0.1 α 1 α 1 -1....
-
[4]
Edge kernel Here, we consider the case when bothxandyare near one of the edges, say the right edge, denoted byx edge. Before proceeding further, we discuss the expressions forx edge and the interparticle distance denoted byW N. From Eqs (132) and (140), it is easy to see that the support of the density vanishes at xedge =z edge N 1 α+2 withz edge = α+2 µ ...
-
[5]
(140), using ¯ρ(z)≃ p µ(α+2) π zα+2 ,(165) in Eq
Kernel close to the origin To analyze the kernelK N(x,y)near the origin, we first es- timate the typical interparticle separationℓ 0 near the origin, using R ℓ0 0 ρN(x)dx≈1/N, or equivalently, Z ℓ0/N1/(α+2) 0 ¯ρ(z)dz≈ 1 N .(164) Close to the origin, say on the positive side, from Eq. (140), using ¯ρ(z)≃ p µ(α+2) π zα+2 ,(165) in Eq. (164), we get ℓ0 = α+2...
work page 2017
-
[6]
Kernel near the edge As mentioned in Eq. (157), near the right edge, we expect the following scaling form for the kernel KN(x,y) = 1 WN Kedge x−x edge WN , y−x edge WN .(A1) The scaling function can be found by evaluating Kedge(u,v) =lim N→∞ WN KN(xedge +WN u,x edge +WN v),(A2) whereK N(x,y)is given in Eq. (115), together with Eqs. (120)– (123). As discus...
-
[7]
Kernel near the origin As discussed in Sec. VI C 3, the typical interparticle sepa- ration close to the origin scales as ℓ0 = α+2 µ 1 α+2 N− 1 α+2 ,(A11) Therefore, near the origin, we expect the scaling form of the correlation kernel as KN(x,y) = 1 ℓ0 K0 x ℓ0 , y ℓ0 .(A12) The scaling function can be found by evaluating K0(u,v) =lim N→∞ ℓ0 KN(ℓ0 u, ℓ0 v)...
- [8]
-
[9]
S. Giorgini, L. P. Pitaevskii, and S. Stringari, Theory of ultra- cold atomic fermi gases, Rev. Mod. Phys.80, 1215 (2008)
work page 2008
- [10]
-
[11]
M. Lewenstein, A. Sanpera, and V . Ahufinger,Ultracold Atoms in Optical Lattices: Simulating quantum many-body systems (Oxford University Press, 2012)
work page 2012
-
[12]
I. M. Georgescu, S. Ashhab, and F. Nori, Quantum simulation, Rev. Mod. Phys.86, 153 (2014)
work page 2014
-
[13]
C. Gross and I. Bloch, Quantum simulations with ultracold atoms in optical lattices, Science357, 995 (2017)
work page 2017
-
[14]
Bloch, Quantum simulations come of age, Nat
I. Bloch, Quantum simulations come of age, Nat. Phys.14, 1159 (2018)
work page 2018
- [15]
-
[16]
O. Morsch and M. Oberthaler, Dynamics of bose-einstein con- densates in optical lattices, Rev. Mod. Phys.78, 179 (2006)
work page 2006
- [17]
- [18]
- [19]
-
[20]
M. Schreiber, S. S. Hodgman, P. Bordia, H. P. Lüschen, M. H. Fischer, R. V osk, E. Altman, U. Schneider, and I. Bloch, Ob- servation of many-body localization of interacting fermions in a quasirandom optical lattice, Science349, 842 (2015)
work page 2015
-
[21]
G. Gauthier, I. Lenton, N. M. Parry, M. Baker, M. J. Davis, H. Rubinsztein-Dunlop, and T. W. Neely, Direct imaging of a digital-micromirror device for configurable microscopic optical potentials, Optica3, 1136 (2016)
work page 2016
-
[22]
X. Qiu, J. Zou, X. Qi, and X. Li, Precise programmable quan- tum simulations with optical lattices, npj Quantum Inf.6, 87 (2020)
work page 2020
-
[23]
H.-T. Wei, E. Ibarra-García-Padilla, M. L. Wall, and K. R. A. Hazzard, Hubbard parameters for programmable tweezer ar- rays, Phys. Rev. A109, 013318 (2024)
work page 2024
-
[24]
C. Tabares, C. Kokail, P. Zoller, D. González-Cuadra, and A. González-Tudela, Programming optical-lattice fermi- hubbard quantum simulators, PRX Quantum6, 030356 (2025)
work page 2025
-
[25]
D. J. BenDaniel and C. B. Duke, Space-charge effects on elec- 22 tron tunneling, Phys. Rev.152, 683 (1966)
work page 1966
-
[26]
T. Gora and F. Williams, Theory of electronic states and trans- port in graded mixed semiconductors, Phys. Rev.177, 1179 (1969)
work page 1969
-
[27]
G. Bastard, J. K. Furdyna, and J. Mycielski, Landau levels and cyclotron resonance in graded mixed semiconductors, Phys. Rev. B12, 4356 (1975)
work page 1975
- [28]
-
[29]
von Roos, Position-dependent effective masses in semicon- ductor theory, Phys
O. von Roos, Position-dependent effective masses in semicon- ductor theory, Phys. Rev. B27, 7547 (1983)
work page 1983
-
[30]
P. Harrison and A. Valavanis,Quantum wells, wires and dots: theoretical and computational physics of semiconductor nanos- tructures(John Wiley & Sons, 2016)
work page 2016
-
[31]
G. Bastard,Wave mechanics applied to semiconductor het- erostructures(New York, NY (USA); John Wiley and Sons Inc., 1990)
work page 1990
-
[32]
T. de Jongh, J. Verstraten, M. Dixmerias, C. Daix, B. Peaude- cerf, and T. Yefsah, Quantum gas microscopy of fermions in the continuum, Phys. Rev. Lett.134, 183403 (2025)
work page 2025
-
[33]
M. Dixmerias, G. Del Vecchio Del Vecchio, C. Daix, J. Ver- straten, T. de Jongh, B. Peaudecerf, P. Le Doussal, G. Schehr, and T. Yefsah, Universal random matrix behavior of a fermionic quantum gas (2025), arXiv:2510.25735
-
[34]
T. Maltezopoulos, A. Bolz, C. Meyer, C. Heyn, W. Hansen, M. Morgenstern, and R. Wiesendanger, Wave-function map- ping of InAs quantum dots by scanning tunneling spectroscopy, Phys. Rev. Lett.91, 196804 (2003)
work page 2003
-
[35]
J. Verstraten, K. Dai, M. Dixmerias, B. Peaudecerf, T. de Jongh, and T. Yefsah, In situ imaging of a single-atom wave packet in continuous space, Phys. Rev. Lett.134, 083403 (2025)
work page 2025
-
[36]
L. D. Landau and E. M. Lifshitz,Quantum mechanics: non- relativistic theory, V ol. 3 (Elsevier, 2013)
work page 2013
-
[37]
S. Albeverio, F. Gesztesy, R. Hoegh-Krohn, and H. Holden, Solvable models in quantum mechanics(Springer Science & Business Media, 2012)
work page 2012
-
[38]
Vicari, Entanglement and particle correlations of fermi gases in harmonic traps, Phys
E. Vicari, Entanglement and particle correlations of fermi gases in harmonic traps, Phys. Rev. A85, 062104 (2012)
work page 2012
-
[39]
Eisler, Universality in the full counting statistics of trapped fermions, Phys
V . Eisler, Universality in the full counting statistics of trapped fermions, Phys. Rev. Lett.111, 080402 (2013)
work page 2013
- [40]
-
[41]
P. Calabrese, P. Le Doussal, and S. N. Majumdar, Random ma- trices and entanglement entropy of trapped fermi gases, Phys. Rev. A91, 012303 (2015)
work page 2015
-
[42]
D. S. Dean, P. Le Doussal, S. N. Majumdar, and G. Schehr, Finite-temperature free fermions and the kardar-parisi-zhang equation at finite time, Phys. Rev. Lett.114, 110402 (2015)
work page 2015
-
[43]
D. S. Dean, P. Le Doussal, S. N. Majumdar, and G. Schehr, Universal ground-state properties of free fermions in a d- dimensional trap, Europhys. Lett.112, 60001 (2015)
work page 2015
-
[44]
D. S. Dean, P. Le Doussal, S. N. Majumdar, and G. Schehr, Noninteracting fermions at finite temperature in ad- dimensional trap: Universal correlations, Phys. Rev. A94, 063622 (2016)
work page 2016
-
[45]
D. S. Dean, P. Le Doussal, S. N. Majumdar, and G. Schehr, Noninteracting fermions in a trap and random matrix theory, J. Phys. A: Math. Theor.52, 144006 (2019)
work page 2019
-
[46]
M. L. Mehta,Random matrices, V ol. 142 (Elsevier, 2004)
work page 2004
-
[47]
P. J. Forrester,Log-gases and random matrices (LMS-34) (Princeton university press, 2010)
work page 2010
-
[48]
K. S. Fa and E. K. Lenzi, Power law diffusion coefficient and anomalous diffusion: Analysis of solutions and first passage time, Phys. Rev. E67, 061105 (2003)
work page 2003
-
[49]
A. G. Cherstvy, A. V . Chechkin, and R. Metzler, Anomalous diffusion and ergodicity breaking in heterogeneous diffusion processes, New J. Phys15, 083039 (2013)
work page 2013
-
[50]
A. L. Stella, A. Chechkin, and G. Teza, Universal singularities of anomalous diffusion in the richardson class, Phys. Rev. E 107, 054118 (2023)
work page 2023
-
[51]
N. Leibovich and E. Barkai, Infinite ergodic theory for hetero- geneous diffusion processes, Phys. Rev. E99, 042138 (2019)
work page 2019
-
[52]
Ray, Expediting feller process with stochastic resetting, Phys
S. Ray, Expediting feller process with stochastic resetting, Phys. Rev. E106, 034133 (2022)
work page 2022
- [53]
-
[54]
A. L. Stella, A. Chechkin, and G. Teza, Anomalous dynami- cal scaling determines universal critical singularities, Phys. Rev. Lett.130, 207104 (2023)
work page 2023
-
[55]
G. Del Vecchio Del Vecchio and S. N. Majumdar, Generalized arcsine laws for a sluggish random walker with subdiffusive growth, J. Stat. Mech.: Theor. Exp.2025, 023207 (2025)
work page 2025
-
[56]
G. Del Vecchio Del Vecchio, M. Kulkarni, S. N. Majumdar, and S. Sabhapandit, Proxitaxis: An adaptive search strategy based on proximity and stochastic resetting, Phys. Rev. E113, L042101 (2026)
work page 2026
-
[57]
D. Boyer and S. N. Majumdar, Ultra slow sub-logarithmic dif- fusion of a sluggish random walker subject to resetting with memory (2026), arXiv:2603.01149 [cond-mat.stat-mech]
- [58]
-
[59]
L. Menon and C. Anteneodo, Random search with resetting in heterogeneous environments, Phys. Rev. E110, 054111 (2024)
work page 2024
-
[60]
N. G. Van Kampen,Stochastic processes in physics and chem- istry, 3rd ed. (Elsevier, Amsterdam, 2007)
work page 2007
-
[61]
M. P. Marder,Condensed matter physics(John Wiley & Sons, 2010)
work page 2010
-
[62]
R. Feynman, R. Leighton, and M. Sands,The Feynman Lectures on Physics, Vol. III: The New Millennium Edition: Quantum Mechanics(Basic Books, 2015)
work page 2015
-
[63]
I. S. Gradshteyn and I. M. Ryzhik,Table of integrals, series, and products(Academic press, 2014)
work page 2014
-
[64]
L. Pitaevskii and S. Stringari, Fermi gas in the harmonic trap, inBose-Einstein Condensation and Superfluidity(Oxford University Press, 2016) https://academic.oup.com/book/0/chapter/195464316/chapter- pdf/45141610/acprof-9780198758884-chapter-17.pdf
-
[65]
M. Abramowitz and I. A. Stegun,Handbook of Mathematical Functions: With Formulas, Graphs, and Mathematical Tables (Dover Publications, New York, 1965)
work page 1965
-
[66]
A. Erdélyi, W. Magnus, F. Oberhettinger, and F. G. Tricomi, Higher Transcendental Functions, V ol. II (McGraw-Hill, New York, 1953)
work page 1953
-
[67]
R. P. Feynman, A. R. Hibbs, and D. F. Styer,Quantum mechan- ics and path integrals(Courier Corporation, 2010)
work page 2010
-
[68]
M. Kulkarni, S. N. Majumdar, and S. Sabhapandit, Dynami- cally emergent correlations in bosons via quantum resetting, J. Phys. A: Math. Theor.58, 105003 (2025)
work page 2025
-
[69]
C. Nadal and S. N. Majumdar, Nonintersecting brownian inter- faces and wishart random matrices, Phys. Rev. E79, 061117 (2009)
work page 2009
-
[70]
B. Lacroix-A-Chez-Toine, P. Le Doussal, S. N. Majumdar, and G. Schehr, Non-interacting fermions in hard-edge potentials, J. 23 Stat. Mech.: Theor. Exp.2018, 123103 (2018)
work page 2018
-
[71]
P. J. Forrester, The spectrum edge of random matrix ensembles, Nucl. Phys. B402, 709 (1993)
work page 1993
-
[72]
S. N. Majumdar and G. Schehr,Statistics of Extremes and Records in Random Sequences(Oxford University Press, 2024)
work page 2024
-
[73]
Szegö,Orthogonal polynomials, V ol
G. Szegö,Orthogonal polynomials, V ol. 23 (American Mathe- matical Soc., 1939)
work page 1939
-
[74]
T. Guhr, A. Müller-Groeling, and H. A. Weidenmüller, Random-matrix theories in quantum physics: common con- cepts, Phys. Rep.299, 189 (1998)
work page 1998
-
[75]
G. Akemann, J. Baik, and P. Di Francesco,The Oxford Hand- book of Random Matrix Theory(Oxford University Press, 2015)
work page 2015
-
[76]
Mitra, Quantum quench dynamics, Annu
A. Mitra, Quantum quench dynamics, Annu. Rev. Condens. Matter Phys.9, 245 (2018)
work page 2018
-
[77]
M. Saha, M. Kulkarni, and A. Dhar, Generalized hydrodynamic description of the page curve–like dynamics of a freely expand- ing fermionic gas, Phys. Rev. Lett.133, 230402 (2024)
work page 2024
-
[78]
Kehrein, Page curve entanglement dynamics in an analyti- cally solvable model, Phys
S. Kehrein, Page curve entanglement dynamics in an analyti- cally solvable model, Phys. Rev. B109, 224308 (2024)
work page 2024
- [79]
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.