Nearest Neighbor Hopping in Quantum Systems

• Nearest neighbor hopping is a quantum process where electrons move between adjacent lattice sites, forming the basis of tight-binding models.
• The method uses Hamiltonians with defined hopping amplitudes to shape band structures in systems such as graphene and multi-orbital chains.
• Extended analyses including disorder and next-nearest neighbor interactions reveal insights into localization, transport, and experimental behavior.

Nearest neighbor hopping refers to the quantum-mechanical process in which particles—typically electrons—move between adjacent sites on a discrete lattice via a kinetic term in the system Hamiltonian. This concept is foundational in tight-binding models, quantum walks, classical Markov processes, and strongly-correlated electron systems. The restriction to nearest neighbor transitions dictates the connectivity and spectral structure of the model and serves as a controlled limit for understanding more complex transport behavior, interactions, and disorder.

1. Formalism and Hamiltonian Structure

The nearest neighbor hopping term is defined by a sum over pairs of adjacent lattice sites, frequently taking the form

$$H_{NN} = \sum_{\langle i, j\rangle} t_{ij}\ (c_i^\dagger c_j + c_j^\dagger c_i)$$

where $c_i^\dagger$ creates a particle at site $i$, and $t_{ij}$ is the hopping amplitude between sites $i$ and $j$. In periodic, translation-invariant systems, $t_{ij} = t$ for all pairs $\langle i,j\rangle$.

In multi-orbital contexts, such as BaCoTe$_2$O$_7$ (Lin et al., 2024), the nearest neighbor hopping Hamiltonian generalizes to

$$H_{NN} = \sum_{\langle i,j\rangle}\sum_{\sigma}\sum_{\alpha,\beta} t^{NN}_{\alpha\beta}\ c^\dagger_{i\alpha\sigma} c_{j\beta\sigma} + \mathrm{h.c.}$$

with hopping matrix elements $t^{NN}_{\alpha\beta}$ resolving both intra- and inter-orbital channels.

Tables for multiorbital hopping matrices:

Orbital Transition	$t^{NN}_{\alpha\beta}$ (eV)	Dominant Physical Path
$d_{3z^2-r^2} \to d_{3z^2-r^2}$	$-0.079$	Apical superexchange
$d_{yz} \to d_{yz}$	$+0.022$	Planar exchange
$d_{xy} \to d_{xy}$	$-0.003$	Plane, orthogonal bond

In the case of graphene (Laissardière et al., 2013), the nearest neighbor hopping amplitude is parameterized as $\gamma_0 \simeq 2.7$ eV between $p_z$ orbitals on adjacent carbon atoms at distance $a = 0.142$ nm.

2. Physical Consequences in Band Structure

Restricting the kinetic energy to nearest neighbor hopping profoundly shapes the system's electronic spectrum.

• For graphene, the nearest neighbor tight-binding model yields two bands $\pm |f(\vec{k})|$, exhibiting particle-hole symmetry about the Dirac energy $E_D = 0$. The function $f(\vec{k})$ is given by

$$f(\vec{k}) = -t \sum_{j=1}^3 e^{i \vec{k}\cdot \vec{\delta}_j}$$

with $\vec{\delta}_j$ the vectors to three first neighbors.

• In one-dimensional chains with random nearest neighbor hopping (Krishna et al., 2020), the spectrum is symmetric about $E = 0$, with universal low-energy singularities in the density of states (DOS) and localization length for sufficiently regular hopping distributions. Specifically, Dyson's result for a random chain yields

$$\rho(E) \sim \frac{1}{|E| (\ln |E|)^3},\quad \xi(E)\sim |\ln|E||$$

near the band center.

Correspondingly, in multi-orbital zigzag chains (BaCoTe$_2$O$_7$), the nearest neighbor model preferentially supports robust staggered antiferromagnetic Mott-insulating behavior at moderate to large $U$ (Lin et al., 2024).

3. Disorder and Universality

Introducing randomness in the hopping amplitudes induces off-diagonal disorder, modifying transport, localization, and spectral properties.

• Uniform, non-singular random distributions preserve Dyson universality: sub-exponential spatial decay of the zero-energy state and universal singularities in DOS and localization length (Krishna et al., 2020).
• For singular hopping distributions $p_\lambda(t) \sim 1/[t \ln^{\lambda+1}(1/t)]$, the universal scaling is lost for $\lambda < 2$, replaced by model-specific exponents in DOS and localization behavior:
  • $\rho(E) \sim |E|^{-1} |\ln E|^{-(\lambda+1)}$
  • Zero-energy state envelope $\psi(r) \sim \exp[-(r/r_0)^{1/\lambda}]$; super-exponentially localized for $\lambda < 1$.

Thus, the nearest neighbor restriction, when combined with disorder, forms the basis for universality classes in one-dimensional localization, random spin chains, and quantum transport.

4. First-Passage Dynamics and Markov Processes

Nearest neighbor hopping on finite intervals forms the backbone of discrete-time and continuous-time Markov chains, random walks, and first-passage problems.

In models with arbitrary, site-dependent rates $r_{i,i+1}, r_{i,i-1}$, the first-passage time distribution from interior site $i$ to absorbing boundaries can be computed via the backward recurrence and generating function (Holehouse et al., 2023): $G_i(s) = \sum_{t=0}^\infty f_i(t) s^t$ where $f_i(t)$ is the probability of first absorption at time $t$. All moments follow from derivatives of $G_i(s)$, and explicit expressions for the mean and variance are available as rational functions of the rates.

A notable feature is the possible bimodality of the first-passage time distribution, arising when the random rates induce effective potential wells and net bias directions on the interval, producing two distinct decay factors in the geometric mixture.

5. Limitations and Extensions Beyond Nearest Neighbor

Nearest neighbor-only models provide a tractable baseline but are often insufficient for quantitative agreement with experiment or for capturing the full range of emergent phenomena.

• In real materials, next-nearest neighbor (NNN) and further hoppings can be appreciable. For BaCoTe$_2$O$_7$, the NNN hopping ($t^{NNN}_{22} \approx 0.124$ eV) actually dominates the NN term and changes the magnetic ground state from staggered AFM to block AFM (Lin et al., 2024).
• In graphene, hopping beyond nearest neighbor ($t' = 0.1\gamma_0$) breaks particle-hole symmetry, shifts the Dirac energy, and quantitatively alters resonance energies and widths associated with adsorbates or vacancies (Laissardière et al., 2013). This suggests that universal features (e.g., minimum conductivity) persist, but experiment-theory matching requires inclusion of at least second-neighbor hopping.

6. Experimental and Theoretical Implications

Experimental implications for nearest neighbor hopping models are significant:

• Spectral asymmetries, conductivity values, and position/width of resonance features in real materials (e.g., functionalized graphene or correlated zigzag chains) cannot be captured with NN-only models; quantitative prediction requires extended hopping terms (Laissardière et al., 2013, Lin et al., 2024).
• In cold atom systems, quantum walks, and nanostructured devices, the tunability of hopping amplitudes enables tests of universality classes and disorder-driven localization transitions as predicted for random NN models (Krishna et al., 2020).
• First-passage and random walk behavior in disordered chains directly informs reaction rates, search strategies, and stochastic time statistics across biophysical, chemical, and condensed matter systems (Holehouse et al., 2023).

A plausible implication is that nearest neighbor hopping forms the minimal mathematical framework for quantum and classical transport on a lattice, but advances in ab-initio modeling, accurate phase diagram determination, and rigorous experimental comparisons depend critically on a proper accounting of extended hopping processes, disorder statistics, and their associated scaling regimes.