Defect-Driven Nonlinear and Nonlocal Perturbations in Quantum Chains

Transport and localization in isolated quantum systems are typically attributed to spatially-extended disorder, leaving the influence of a few controllable defects largely unexplored despite their relevance to engineered quantum platforms. We introduce an analytic framework showing how a single defect profoundly reshapes wave-function spreading on a finite isolated and periodic tight-binding lattice. Adapting the defect technique from classical random-walk studies, we obtain exact time-resolved site-occupation probabilities and several observables of interest. Even one defect induces striking nonlinear and nonlocal effects, including non-monotonic suppression of transport, enhanced localization at distant sites, and strong sensitivity to the initial particle position at long times. These results demonstrate that minimal perturbations can generate unexpected long-time transport signatures, establishing a microscopic defect-driven mechanism of quantum localization.

💡 Research Summary

This paper investigates how a single localized on‑site defect reshapes quantum transport in a finite, periodic tight‑binding chain. Starting from the standard one‑dimensional tight‑binding Hamiltonian with nearest‑neighbour hopping γ and periodic boundary conditions, the authors add an on‑site energy q at site nd. The particle is initially localized at site n₀, and the time‑dependent wavefunction ψ(n,n₀,t) obeys the Schrödinger equation with a δ‑function term at the defect.

In the defect‑free case (q = 0) the Green’s function can be written analytically in terms of Fourier modes, yielding a ballistic mean‑square displacement (MSD) Δ₂(t) ≈ 2γ²t² at short times. After a characteristic time t★ ≈ aN/γ (set by the Lieb‑Robinson bound) the dynamics feel the periodic boundary, leading to oscillations and a steady‑state MSD that scales with the system size N.

To treat the defect analytically the authors import the “defect technique” from classical random‑walk theory. In Laplace space the full propagator satisfies a self‑consistent equation G′ = G + G Π G′, where Π encodes the single‑site perturbation. Solving this yields the exact Laplace‑domain wavefunction ψ̃(n,n₀,ε) = G̃(n,n₀,ε) + G̃(n,nd,ε) Φ(ε) with Φ(ε)=iq G̃(nd,n₀,ε)/(1−iq G̃(nd,nd,ε)). The unperturbed Green’s function G̃ is expressed through Chebyshev polynomials, allowing Φ(ε) to be written as a rational function whose poles are the roots of a polynomial Q(ε). Inverting the Laplace transform gives Φ(t) as a sum over residues, i.e. a finite set of exponentially oscillating modes.

Transforming back to the time domain, the site‑occupation probability in the presence of the defect becomes
P⁽ᵈ⁾ₙ(t)=Pₙ(t)+I⁽ᵈ⁾ₙ(t)+K⁽ᵈ⁾ₙ(t),
where Pₙ(t) is the defect‑free probability, I⁽ᵈ⁾ₙ(t) is an interference term (linear in Φ) and K⁽ᵈ⁾ₙ(t)=|Aₙ(t)|² is a purely defect‑induced contribution with Aₙ(t)=∫₀ᵗ G(n,nd,t−τ)Φ(τ)dτ. Both I and K depend non‑linearly on the defect strength q and on the three distances |n−n₀|, |n−nd|, and |n₀−nd|, reflecting a fundamentally non‑local response.

Two distinct regimes emerge. When the particle starts on the defect (n₀ = nd), increasing q monotonically enhances localization at nd, suppresses probabilities elsewhere, and drives the mean displacement and MSD down to zero as q→∞, where the wavefunction collapses to |nd⟩. Conversely, when n₀ ≠ nd, the dependence on q is non‑monotonic: the probabilities at the initial site and at the defect first decrease (or increase) with q, reach an extremum, and then reverse direction. The steady‑state mean displacement and MSD also display a minimum as a function of q. This reveals a persistent memory of the initial position that survives indefinitely, even after long‑time averaging.

In the extreme limit q→∞ with n₀ ≠ nd, the authors obtain a closed‑form expression for Φ(t) involving a sum over cosine modes with frequencies 2γ cosθₖ (θₖ=π(2k−1)/N). Time‑averaging the full probability yields an exact stationary distribution:

P⁽ᵈ⁾ₙ = (3/2N)