Sluggish quantum mechanics of noninteracting fermions with spatially varying effective mass

This paper introduces "sluggish quantum mechanics" for noninteracting fermions with a position-dependent effective mass, providing exact solutions for their ground-state properties and revealing that the many-body correlation kernel near the origin is a novel sum of two Bessel kernels, distinct from the standard Bessel or Airy forms found in conventional trapped systems.

The Gist

Imagine you are walking through a forest. In a normal forest, the ground is flat and firm everywhere. You can walk at a steady pace, whether you are near the entrance or deep in the woods. This is how most quantum particles behave in standard physics: they move freely, like a car on a smooth highway.

But in this paper, the authors imagine a very different kind of forest.

The "Sluggish" Forest

In this new forest, the ground gets progressively heavier and stickier the further you walk from the center.

• Near the center (the origin): The ground is firm. You can move quickly.
• Far from the center: The ground turns into thick, heavy mud. Every step you take requires more effort. The further you go, the "heavier" you become, until you are practically stuck in place.

The authors call this "Sluggish Quantum Mechanics."

In the real world, this isn't just a story. Scientists can build this kind of "sticky forest" using ultracold atoms trapped in optical lattices (grids made of laser light). By carefully shaping the lasers, they can make the atoms feel like they are getting heavier as they move away from the center.

The Single Particle: A Heavy Hiker

The paper first looks at a single particle (a single hiker) in this forest.

• Without a fence: If there's no fence, the particle tries to wander off. But because it gets so heavy far away, it doesn't spread out like a normal wave. Instead, its movement is "suppressed." It's like trying to run a marathon in ankle-deep mud; you just can't go very far.
• With a fence (The Trap): Now, imagine putting a fence around the hiker that gets stronger the further out you go. The authors found a specific shape for this fence (a "potential") that makes the math work perfectly. Even though the ground is sticky, the particle can still dance in a predictable pattern. They figured out exactly how the particle moves and where it is likely to be found.

The Crowd: The Fermion Party

The most exciting part of the paper happens when you put many particles in this forest. Specifically, they look at fermions.

In the quantum world, fermions are like extremely polite (but antisocial) party guests. They follow the Pauli Exclusion Principle: No two fermions can stand in the exact same spot. They hate crowding.

• The Normal Party (Standard Physics): If you put many fermions in a normal, flat forest with a circular fence, they spread out evenly. If you look at the crowd from above, it looks like a perfect circle. The density is highest in the middle and fades out at the edges. This is a famous pattern called the "Wigner Semicircle."
• The Sluggish Party: Now, put those same antisocial guests into the "sticky forest."
  • Because the ground gets heavy far away, the guests are pushed toward the center.
  • BUT, because they are fermions, they also hate being on top of each other.
  • The Result: They get stuck in a weird middle ground. They can't go far (because of the mud), but they can't all pile up in the center (because they hate each other).
  • The Surprise: Instead of a crowd that is densest in the middle, the crowd actually forms a ring or a donut shape. The center of the forest becomes empty! The guests are too "heavy" to stay in the middle, but too "sticky" to get to the edge. They end up hovering in a ring around the center.

The Secret Code: The Kernel

The authors didn't just guess this; they wrote down the exact mathematical "recipe" (called a kernel) that predicts exactly where every particle will be.

• In normal physics, this recipe uses standard mathematical shapes (like the Airy function or the Bessel function).
• In this "sluggish" physics, the recipe is brand new. It's a hybrid of two different Bessel functions. It's a unique mathematical fingerprint that has never been seen before in this context.

Why Does This Matter?

1. New Physics: It shows that changing how "heavy" a particle feels depending on where it is creates entirely new behaviors that we haven't seen before.
2. Real Experiments: This isn't just math on a chalkboard. With modern lasers, scientists can actually build these "sticky forests" in a lab. They can watch atoms behave exactly as the paper predicts.
3. Random Matrices: The math connects to a field called "Random Matrix Theory," which is used to understand everything from the energy levels of heavy atoms to the spacing of stars in a galaxy. This paper adds a new chapter to that book.

Summary

Think of this paper as a guidebook for a new kind of quantum world where distance equals weight.

• Normal World: Particles spread out evenly.
• Sluggish World: Particles get heavy as they move away, causing them to avoid the center and form a ring.
• The Discovery: The authors found the exact mathematical rules for this world, proving that even in a "sticky" universe, quantum particles follow a beautiful, predictable, and surprisingly ring-shaped pattern.

Technical Summary

1. Problem Statement

The paper investigates a class of one-dimensional quantum systems where the kinetic energy term is spatially dependent, arising from the continuum limit of an inhomogeneous tight-binding model with position-dependent hopping amplitudes.

• Core Concept: The system is characterized by an effective mass that scales as $m_{\text{eff}}(x) = m_{\text{eff}} |x|^\alpha$ with $\alpha > 0$. This leads to "sluggish quantum mechanics," where particle motion is progressively suppressed at larger distances (the particle becomes "heavier" as it moves away from the origin).
• The Challenge: While the standard Schrödinger equation with constant mass is well-understood, the authors aim to solve the BenDaniel–Duke form of the Schrödinger equation with this specific power-law varying mass. They specifically address:
  1. The exact solution for a free particle (no external potential).
  2. The exact solution in the presence of a confining potential $V_{\text{ext}}(x) = \frac{1}{2}m_{\text{eff}}(x)\omega^2 x^2 = \frac{1}{2}m_{\text{eff}}\omega^2 |x|^{\alpha+2}$.
  3. The many-body ground state properties of $N$ noninteracting spinless fermions in this trap, focusing on the joint probability density function (JPDF) and correlation kernels.

2. Methodology

The authors employ a combination of analytical techniques from quantum mechanics, special functions, and random matrix theory (RMT).

• Derivation: They derive the sluggish Schrödinger equation by taking the continuum limit ($\delta \to 0$) of a tight-binding model with spatially varying hopping elements $a(x_i, x_i \pm \delta)$. This yields the operator $-\partial_x [D(x) \partial_x]$, where $D(x) = \hbar^2 / [2m_{\text{eff}}(x)]$.
• Single-Particle Solutions:
  • Free Case: The time-independent Schrödinger equation is transformed into a Bessel differential equation via an ansatz involving power laws. The solutions are expressed in terms of Bessel functions $J_{\pm \tilde{\nu}}$.
  • Confined Case: For the potential $V_{\text{ext}} \propto |x|^{\alpha+2}$, the equation is mapped to the Kummer confluent hypergeometric equation. The physical solutions are identified as Tricomi functions $U(a,b,z)$, which reduce to generalized Laguerre polynomials $L_n^{(\gamma)}(z)$ for quantized energy levels.
  • Mapping: The authors demonstrate a variable transformation that maps the sluggish equation to a standard Schrödinger equation with a constant mass but a modified potential consisting of a harmonic term plus an inverse-square term ($\sim 1/y^2$).
• Many-Body Formalism:
  • The $N$-fermion ground state is constructed as a Slater determinant of the single-particle eigenfunctions.
  • The system is identified as a Determinantal Point Process (DPP). The $n$-point correlation functions are determined by a kernel $K_N(x,y) = \sum_{k=0}^{N-1} \psi_k^*(x)\psi_k(y)$.
  • The kernel is split into symmetric (even) and antisymmetric (odd) sectors due to the parity of the potential.
• Asymptotic Analysis: The authors use the Christoffel–Darboux formula to sum the kernel series and apply Plancherel-Rotach asymptotics for Laguerre polynomials to derive scaling limits for large $N$ in the bulk, at the edge, and near the origin.

3. Key Contributions and Results

A. Exact Solvability of Single-Particle Dynamics

• Free Particle: Exact eigenfunctions and the full quantum propagator $G(x,t|x_0,0)$ were derived for the free sluggish particle. The propagator involves Bessel functions and exhibits anomalous spreading compared to standard free particles.
• Confined Particle: The authors identified the specific confining potential $V_{\text{ext}} \propto |x|^{\alpha+2}$ that renders the system exactly solvable.
  • Spectrum: The energy levels are discrete and given by $E_n \propto (2n + \text{shift})$. The spacing is not uniform for $\alpha > 0$, unlike the harmonic oscillator ($\alpha=0$).
  • Eigenfunctions: The wavefunctions are expressed in terms of generalized Laguerre polynomials with a specific weight function.
  • Propagator: Closed-form expressions for the symmetric and antisymmetric propagators were obtained using the Hardy–Hille formula.

B. Many-Body Ground State and Density Profile

• Density Depletion: A striking result is the behavior of the average density profile $\rho_N(x)$ for large $N$.
  • For $\alpha = 0$ (standard harmonic trap), the density follows the Wigner semicircle law, peaking at the origin.
  • For $\alpha > 0$, the density profile has a finite support but exhibits a non-monotonic shape with a vanishing minimum at the origin.
  • Physical Origin: This "depletion zone" at $x=0$ is a purely quantum effect. While the classical stationary distribution (Stratonovich case) would peak at the origin, the quantum excited states have wavefunctions that are suppressed near the origin (peaked away from $x=0$). The sum of these probabilities results in a dip in the total density.
• Scaling: The density scales as $\rho_N(x) \sim N^{-1/(\alpha+2)} \bar{\rho}(x/N^{1/(\alpha+2)})$.

C. Novel Correlation Kernels

The paper's most significant theoretical finding concerns the scaling behavior of the correlation kernel $K_N(x,y)$ in the large $N$ limit:

1. Bulk: In the bulk of the distribution, the kernel scales to the standard Sine kernel ($\sin(\pi w)/\pi w$), identical to standard trapped fermions.
2. Edge: Near the spectral edge, the kernel scales to the Airy kernel, also consistent with standard RMT results.
3. Origin (The New Discovery): Near the origin ($x=0$), the scaling kernel is neither the Bessel kernel nor the Airy kernel.
   • It is a new kernel given by the sum of two Bessel kernels with different indices ($\pm \gamma$, where $\gamma = (\alpha+1)/(\alpha+2)$).
   • Formula: $\tilde{K}_0(u', v') = \frac{1}{2} [K_{\text{Bessel}}^{(\gamma)}(u', v') + K_{\text{Bessel}}^{(-\gamma)}(u', v')]$.
   • This structure arises because the system requires both even and odd sectors to describe the full space, unlike the half-space problems (e.g., hard-wall potentials) which typically yield a single Bessel kernel.

D. Half-Space JPDF

For the system restricted to a half-space ($x>0$) with Dirichlet or Neumann boundary conditions, the authors derived an explicit JPDF.

• The JPDF maps exactly to the eigenvalue distribution of Laguerre (Wishart) random matrices with a specific parameter $\gamma$.
• This provides a direct link between sluggish quantum fermions and established RMT ensembles, allowing for the computation of gap probabilities and extreme-value statistics.

4. Significance and Implications

• Theoretical Novelty: The work introduces "sluggish quantum mechanics" as a solvable framework for systems with position-dependent effective mass. The discovery of a new kernel (sum of two Bessel kernels) at the origin expands the classification of universality classes in random matrix theory and determinantal point processes.
• Experimental Relevance: The model is directly relevant to engineered optical lattices with ultracold atoms. Using digital micromirror devices or shaped optical potentials, experimentalists can create spatially varying tunneling amplitudes, effectively realizing the $m_{\text{eff}}(x) \propto |x|^\alpha$ scenario.
• Quantum vs. Classical: The paper highlights a profound difference between classical sluggish diffusion and quantum sluggish dynamics. While classical diffusion (Stratonovich) predicts accumulation at the origin, the quantum system exhibits depletion due to the nodal structure of excited states.
• Future Directions: The authors suggest extending this framework to finite temperatures, higher dimensions, and dynamical properties like quantum quenches and entanglement growth.

In summary, this paper provides a comprehensive analytical solution for a new class of quantum systems, revealing unique many-body correlations and a novel universality class for the correlation kernel at the origin, with direct pathways for experimental verification in cold atom setups.