---
title: 'Andreev Crystals: Periodic Superconducting States'
url: https://www.emergentmind.com/topics/andreev-crystals
type: topic
---

# Andreev Crystals: Periodic Superconducting States

Andreev crystals are periodic or symmetry-structured manifestations of Andreev physics—coherent electron–hole conversion at superconducting interfaces—in which bound states, conductance spectra, or collective modes are organized by a lattice, a vortex array, a magnetic superstructure, or a phase-biased junction network. In the cited literature, the expression is used in several related senses rather than as a single universally fixed term: a triangular Abrikosov vortex lattice whose cores host Caroli–de Gennes–Matricon (CdGM) states in hole-doped BaFe\(_2\)As\(_2\) [1005.4038]; single-crystal Andreev spectroscopies that use crystallographic orientation to resolve anisotropic or multiband superconductivity [1604.06325, 1009.1572]; quasi-one-dimensional superconductors with periodic magnetic regions whose local Andreev states hybridize into Bloch subgap bands [2007.05350, 2104.06766]; and hybrid Josephson arrays in which short-junction Andreev states form directional bands under uniform phase bias [2508.11768].

## 1. Conceptual scope and common mechanisms

The common element in all usages is Andreev reflection: at a normal-metal/superconductor interface, a subgap electron is converted into a hole while a Cooper pair is transferred into the condensate. When this process is confined by geometry, magnetic texture, vortex structure, or crystal orientation, it produces Andreev bound states (ABS), surface resonances, or subgap bands. In vortex cores, the relevant levels are the CdGM states with energies \(E_\mu \approx \mu \omega_0\), \(\mu=\pm 1/2,\pm 3/2,\ldots\), and \(\omega_0\sim \Delta^2/E_F\) [1005.4038]. In multiple-Andreev-reflection (MAR) junctions, the diagnostic structure is the subharmonic series \(V_n=2\Delta/(en)\) [1404.0517, 1409.0123]. In magnetic or hybrid arrays, periodic hybridization converts localized ABS into subgap bands whose topology and transport depend on phase, symmetry, or magnetization texture [2007.05350, 2104.06766, 2508.11768].

| Usage of the term | Defining structure | Representative papers |
|---|---|---|
| Vortex-lattice sense | Periodic Abrikosov lattice with CdGM states in each core | [1005.4038] |
| Single-crystal spectroscopic sense | Directional PCAR or MAR on anisotropic or multigap superconductors | [1604.06325], [1404.0517] |
| Magnetic-chain sense | Periodic magnetic regions in a 1D superconducting wire | [2007.05350], [2104.06766] |
| Hybrid-array sense | Periodic Josephson fingers forming phase-biased ABS bands | [2508.11768] |
| Related topological/interfacial sense | Crystal-symmetry-organized ABS or magnetic crystallization waves | [2301.05340], [2502.15404] |

This plurality of meanings is significant. In some papers, “Andreev crystal” denotes a literal periodic array of bound-state-hosting objects. In others, it denotes a methodological program: using single crystals, directional point contacts, field rotation, or break-junction stacks to make the superconducting order parameter visible through Andreev spectroscopy.

## 2. Vortex-lattice realization in Ba\(_{0.6}\)K\(_{0.4}\)Fe\(_2\)As\(_2\)

The most literal experimental realization of an Andreev crystal in the vortex-lattice sense was reported in hole-doped Ba\(_{0.6}\)K\(_{0.4}\)Fe\(_2\)As\(_2\) by low-temperature STM/STS at \(T=2\) K and magnetic fields \(B=4\) T and \(9\) T [1005.4038]. There the defining hallmark was the simultaneous observation of a triangular Abrikosov lattice and Andreev bound states in every measured vortex core. The magnetic flux per vortex was \(\approx 2.17\times 10^{-15}\) Wb at \(9\) T and \(\approx 1.98\times 10^{-15}\) Wb at \(4\) T, consistent with \(\Phi_0=h/(2e)=2.07\times 10^{-15}\) Wb. The corresponding lattice spacings, from \(a\approx \sqrt{2\Phi_0/(\sqrt{3}B)}\), were \(\approx 52\) nm at \(9\) T and \(\approx 77\) nm at \(4\) T, in agreement with the STM maps [1005.4038].

Core spectroscopy revealed a pronounced negative-bias conductance peak centered at \(\approx -0.5\) mV at the vortex center. Along a \(104\) Å line through a vortex, this single peak evolved into two sub-peaks with dominant spectral weight at negative bias and then faded continuously into the zero-field superconducting spectrum [1005.4038]. In the same study, bright domains exhibited two gaps, \(\Delta\approx 4\) meV and \(\approx 8\) meV, whereas dark domains showed predominantly the smaller gap and higher zero-bias conductance. The observed \(|E|\approx 0.5\) meV core feature was interpreted as compatible with the \(|\mu|=1/2\) CdGM level, implying \(\omega_0\approx 1.0\) meV and, through \(\omega_0\sim \Delta^2/E_F\), effective \(E_F\) estimates of \(\approx 16\) meV for \(\Delta\approx 4\) meV and \(\approx 64\) meV for \(\Delta\approx 8\) meV [1005.4038].

The same work emphasized the contrast with electron-doped Ba(Fe,Co)\(_2\)As\(_2\), where vortices were highly disordered and Andreev bound states were not observed in the cores [1005.4038]. The ordered lattice in the hole-doped material was attributed to weaker disorder-induced pinning from off-plane K dopants, in contrast to in-plane Co dopants in the FeAs layers. Within this usage, an Andreev crystal is therefore a periodic Abrikosov vortex array whose individual cores remain clean enough for CdGM quantization to be spectroscopically resolved.

## 3. Directional Andreev spectroscopy on single crystals

A second major usage treats “Andreev crystals” as single-crystalline platforms in which PCAR, MAR, or related Andreev spectroscopies resolve anisotropy, multigap structure, and pairing symmetry by exploiting crystal orientation. In this sense the crystal is not itself an array of bound states; rather, crystallographic direction enters the Andreev process as a selection rule on \(\Delta(\mathbf{k})\), \(v_F\), interface transparency, and orbital weight.

Field-angle-dependent PCAR on La(O,F)BiSeS provided a clear example. With current injected predominantly along the \(c\) axis, the superconducting gap was \(\Delta_c\approx 0.61\) meV at \(1.5\) K, the gap decreased faster for \(H\parallel c\) than for \(H\parallel a\), and at \(H=8\) kG the extracted values were \(\Delta(H\parallel c)\approx 0.28\) meV and \(\Delta(H\parallel a)\approx 0.40\) meV, giving an anisotropy factor \(\approx 1.42\). The spectra and angular magnetoresistance exhibited two-fold symmetry, interpreted as evidence for anisotropic and unconventional pairing [1604.06325].

In iron pnictides, directional PCAR and MAR established a broad multigap program. Ba(Fe\(_{1-x}\)Co\(_x\))\(_2\)As\(_2\) showed two nodeless gaps, with low-temperature averages \(\Delta_1\approx 4\)–\(4.5\) meV and \(\Delta_2\approx 9\)–\(10\) meV, and additional conductance structures from which a characteristic boson energy \(\Omega_0\approx 12\) meV was extracted; these spectra were modeled within a three-band \(s^\pm\) Eliashberg framework [1009.1572]. In BaFe\(_{2-x}\)Ni\(_x\)As\(_2\), c-axis PCAR tracked a doping evolution from underdoped plateau-like spectra to overdoped sharp in-gap peaks, with optimal-doping fits yielding \(\Delta_h=10.5\) or \(9.5\) meV, \(\Delta_e=4.2\) or \(4.5\) meV, and an electron-gap anisotropy parameter \(r=0.3\) [1106.2891]. A three-dimensional BTK analysis further distinguished Ca(Fe\(_{0.94}\)Co\(_{0.06}\))\(_2\)As\(_2\), where the small gap required nodes or zeros, from ab-plane Ba(Fe\(_{0.92}\)Co\(_{0.08}\))\(_2\)As\(_2\), where two nodeless gaps sufficed; by contrast, c-axis Ba-122 films could not be fitted by a two-nodeless-gap 3D model and were discussed in terms of hot spots or 3D minima [1301.3782].

Break-junction MAR on LiFeAs resolved three gaps directly from subharmonic gap structures: \(\Delta_\Gamma=5.1\)–\(6.5\) meV, \(\Delta_L=3.8\)–\(4.8\) meV, and \(\Delta_S=0.9\)–\(1.9\) meV, with estimated anisotropies \(<8\%\), \(\approx 12\%\), and \(\approx 20\%\), respectively [1404.0517]. Combined lower-critical-field and IMARE work on 122 arsenides supported a nodeless multigap state with \(\Delta_L\approx 6\)–\(8\) meV, \(k\)-space anisotropy of approximately \(25\)–\(30\%\), and \(\Delta_S\approx 1.7\pm 0.3\) meV [1409.0123]. Sub-kelvin PCAR on FeSe yielded an anisotropic two-gap fit with \(\Delta_S=0.495\) meV, \(\Delta_L=1.17\) meV, weight \(w=0.317\), barrier \(Z=0.674\), and anisotropy parameter \(\alpha=0.31\); the study argued that \(0.5\) K operation was necessary to resolve the multiband structure unambiguously [1910.02698].

Non-centrosymmetric and topological-candidate systems sharpened the distinction between conventional Andreev spectroscopy and genuinely unconventional ABS phenomenology. Directional PCAR on BiPd found \(\Delta_1\approx 0.4\pm 0.1\) meV together with orientation-selective secondary gaps \(\Delta_2\approx 0.21\pm 0.05\) meV and \(\Delta_3\approx 0.8\pm 0.15\) meV, plus a pronounced zero-bias conductance peak attributed to ABS from a mixed-parity order parameter [1202.2454]. By contrast, PCAR on SnAs gave \(\Delta\approx 0.62\) meV and \(\Delta(0)/(k_B T_c)\approx 1.78\), consistent with weak-coupling BCS superconductivity [2111.14425]. In Nb\(_x\)Bi\(_2\)Se\(_3\), low- and high-resistance Andreev spectra displayed zero-bias conductance enhancement and conductance dips at the superconducting gap that were described as inconsistent with conventional BTK and discussed as consistent with \(p\)-wave pairing in a topological-superconductor candidate [1707.08516].

## 4. One-dimensional magnetic Andreev crystals

A third usage is strictly theoretical and highly specific: a quasi-one-dimensional \(s\)-wave superconducting wire with a periodic arrangement of mesoscopic magnetic regions. In this construction, each magnetic region is a semiclassical impurity that does not generate normal backscattering but instead imparts a spin- and particle–hole-dependent phase \(\Phi=(\hbar v_F)^{-1}\int dx\, h(x)\). The local ABS associated with individual regions hybridize with an overlap controlled by \(e^{-a/\xi}\), where \(a\) is the period and \(\xi\) the coherence length, thereby forming Bloch subgap bands [2007.05350, 2104.06766].

For collinear antiferromagnetic arrangements, the tight-binding description yields
\[
\frac{\epsilon_\sigma(k)}{\Delta}
=
\pm
\sqrt{
\frac{\omega_0^2+4t^2\sin^2(ka)}
{1+\omega_0^2+4t^2\sin^2(ka)}
},
\]
with \(\omega_0=\cot\Phi\) and \(t\propto -e^{-a/\xi}/\sin\Phi\) [2007.05350]. The excitation gap therefore closes at half-integer \(\Phi/\pi\), where \(\omega_0=0\), producing Dirac phase boundaries. The same work related the spectral asymmetry index of the spin-resolved Hamiltonian to the spin polarization and used it as the observable that quantifies gap closing and reopening. Heterojunctions between neighboring antiferromagnetic phases support spin-polarized interface bound states and exhibit fractionalization of the interface spin, in close analogy to Jackiw–Rebbi mass inversion [2007.05350].

The subsequent extension to “spectral properties of Andreev crystals” generalized this framework to helical magnetic textures and to regimes beyond nearest-neighbor coupling using Eilenberger theory [2104.06766]. In the helical case, the in-gap spectrum consists of a pair of energy-symmetric bands. Ferromagnetic configurations allow band crossings without hybridization, whereas generic helical configurations produce a Dirac point at \(\Phi=\pi/2\mod\pi\) and \(ka=\pi/2\), followed by an inverted gap upon further increase of magnetic strength [2104.06766]. The Eilenberger treatment computed local density of states and local spin polarization for infinite chains and semi-infinite junctions, and showed that junctions may exhibit interface bound states together with fractionalization of the surface spin polarization. In this usage, an Andreev crystal is a bona fide one-dimensional crystal of hybridized ABS.

## 5. Hybrid Josephson-junction arrays

A more recent transport-theoretic usage concerns periodic arrays of short superconducting segments embedded in a proximitized semiconducting channel. When the superconducting finger length \(l_S\) is comparable to the coherence length, Andreev bound states localized near opposite NS interfaces hybridize within a segment and then across the array, generating subgap bands below the parent gap [2508.11768]. The formalism uses concatenated scattering matrices for superconducting segments and normal spacers, together with the Andreev amplitudes
\[
r_{he}(x,\phi)=\frac{1-t_S^2(x)}{1-r_\Delta^2(x)t_S^2(x)}\,r_\Delta(x)e^{-i\phi},
\qquad
t_{ee}(x)=\frac{1-r_\Delta^2(x)}{1-r_\Delta^2(x)t_S^2(x)}\,t_S(x)e^{ik_{FS}l_S},
\]
and shows that high transparency together with a constant phase bias between neighboring superconductors produces directional bands [2508.11768].

The characteristic resonance condition of the minimal NSNSN structure is \(x=\cos(\phi/2)\), with \(x=\omega/\Delta\), reproducing the short-junction ABS alignment \(E\sim \Delta\cos(\phi/2)\) in the high-transparency limit [2508.11768]. In extended arrays, the nonlocal conductances \(G_{RL}(V,\phi)\) and \(G_{LR}(V,\phi)\) form distinct subgap bands in the \((V,\phi)\) plane. The central claim is that one band consists only of right-moving, and the other only of left-moving, electronic states, enabling a flux- and bias-voltage-tunable one-way filter [2508.11768]. Phase control can be implemented either by superconducting loops, where \(\delta\phi=2\pi\Phi/\Phi_0\), or by a uniform current-induced phase gradient.

This version of the concept departs from vortex physics and from spectroscopic nomenclature. It treats the Andreev crystal as an engineered band structure whose transport nonreciprocity derives from coherent phase-biased hybridization of ABS in a periodic Josephson environment.

## 6. Related topological and interfacial extensions

Several adjacent literatures extend the same logic—crystal structure organizing Andreev phenomena—without always using exactly the same definition. In helical crystals with point group \(D_6\), interlayer pairings along the helical hopping produce surface ABS and zero-energy peaks in the surface local density of states. For the \((001)\) surface, three irreducible representations exhibit zero-energy peaks; for the zigzag surface, four do. A one-dimensional winding number,
\[
w(k_\parallel)= -\frac{1}{4\pi i}\int dk_\perp \, \mathrm{tr}\!\left[\Gamma H_{\mathrm{BdG}}^{-1}(\mathbf{k})\partial_{k_\perp}H_{\mathrm{BdG}}(\mathbf{k})\right],
\]
diagnoses these states and establishes their topological origin [2301.05340].

In type-II Weyl-semimetal/superconductor junctions, crystalline anisotropy produces a distinctly different “Andreev-crystal” analogy: an incident electron can undergo double Andreev reflection into two hole beams, one retro and one specular, symmetric about the interface normal, with a critical orientation \(\theta_c=\arccos(v_2/v_1)\) separating double-Andreev and normal-reflection regimes [1705.00904]. Here the crystal analogy is optical: the tilted Weyl cone acts as an electronic birefringent medium for Andreev processes.

A farther extension appears at the solid/liquid \(^3\)He interface, in direct lineage from A. F. Andreev’s crystallization-wave theory. Below the Néel temperature, magnetic crystallization waves acquire an additional spin-supercurrent contribution to their inertial mass,
\[
M_{\mathrm{eff}}(q)=M_m(q)+M_s(H)
=
\rho d\left[\left(\frac{H}{H_0}\right)^2+\left(\frac{\Delta\rho}{\rho}\right)^2\frac{1}{qd}\right],
\]
so that the inertial mass differs from the gravitational mass of the interfacial mode [2502.15404]. This usage is not about Andreev reflection at all, but it preserves the deeper theme of periodic interfacial dynamics shaped by Andreev’s ideas.

Taken together, these works show that “Andreev crystals” does not denote a single material class. It denotes a family of constructions in which coherent electron–hole superpositions are ordered by periodicity, vortex matter, magnetic texture, crystal symmetry, or phase bias. The most restrictive definition is the periodic array of bound-state-hosting objects—vortex cores, magnetic impurities, or Josephson elements. The broadest definition is methodological: a crystalline setting in which Andreev spectroscopy resolves the momentum structure, topology, or collective organization of superconducting quasiparticles.

Source: https://www.emergentmind.com/topics/andreev-crystals