Scanning Quantum-Vortex Microscopy

• Scanning quantum-vortex microscopy is a set of scanning-probe techniques that quantitatively image quantized vortices, revealing key parameters like position, lattice order, and bound states.
• Techniques such as SQUID, NV magnetometry, STM/S, and MFM-based methods offer distinct spatial resolutions and measurable observables, from local stray fields to quasiparticle densities.
• The methodology enables detailed analysis of vortex dynamics, pinning landscapes, and superconducting order through real-space imaging and reciprocal-space metrics.

Scanning quantum-vortex microscopy denotes a set of scanning-probe methodologies for real-space interrogation of quantized vortices, chiefly Abrikosov and Pearl vortices in type-II superconductors, with one distinct branch extending the concept to magnetic-vortex cores in ferromagnetic nanodisks. Depending on the modality, the measured observable is the local magnetic stray field, the local quasiparticle density of states, the cantilever response of a deliberately dragged vortex, or the dissipative and reactive response of a locally driven vortex configuration. Across these implementations, the common objective is quantitative access to vortex position, lattice symmetry, orientational and positional order, pinning, bound states, and vortex dynamics (Thiel et al., 2015, Suderow et al., 2014, Horn et al., 2023, Hovhannisyan et al., 2024).

1. Physical basis of vortex imaging

In a type-II superconductor, a single Abrikosov vortex carries one quantum of magnetic flux,

$$\Phi_0 = h/2e \simeq 2.07\times 10^{-15}\,\mathrm{Wb}.$$

Within the London approximation, and for distances $r\gg\xi$, the stray field just above the surface is

$B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$

where $\lambda$ is the London penetration depth and $K_0$ is the zeroth-order modified Bessel function. The corresponding London equation,

$$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$

expresses the screening-current origin of the vortex field, while Ginzburg–Landau theory introduces the order parameter $\psi(\mathbf r)$ and the free-energy functional in which vortices appear as phase singularities where $\psi\to0$ (Wells et al., 2018).

For thin films with thickness $d\ll\lambda_L$, the relevant object is a Pearl vortex, characterized by the Pearl length

$\Lambda = 2\lambda_L^2/d.$

In this regime the stray-field profile differs materially from the frequently used monopole approximation, and quantitative nanoscale measurements at $r\gg\xi$0 were shown to agree with Pearl’s analytic model rather than the monopole picture. This distinction is crucial because the monopole model cannot disentangle changes in $r\gg\xi$1 from changes in probe height and fails at $r\gg\xi$2 (Thiel et al., 2015).

Not all variants of scanning quantum-vortex microscopy are magnetic-field microscopes in the narrow sense. In STM/S implementations, contrast derives from the local quasiparticle density of states $r\gg\xi$3 rather than from direct field mapping. The differential conductance satisfies $r\gg\xi$4 for a normal-metal tip at $r\gg\xi$5, so the vortex core is imaged on the coherence-length scale $r\gg\xi$6, not on the penetration-depth scale $r\gg\xi$7. This allows direct access to Caroli–de Gennes–Matricon bound states with energies $r\gg\xi$8 and to core anisotropy set by gap structure or Fermi-surface anisotropy (Suderow et al., 2014).

A separate physical mechanism underlies MFM-based SQVM. There the operative quantity is the tip–vortex interaction energy and its gradient. In thin films, the vortex can be trapped and dragged when the lateral tip-induced force exceeds the restoring pinning force. The technique therefore converts vortex motion itself into a nanoscale probe of the pinning potential, with a lateral resolution limited by $r\gg\xi$9 rather than by the magnetic spot size of a static vortex image (Aladyshkin et al., 7 Jul 2025).

2. Instrumental realizations

The literature groups several technically distinct instruments under the same general label. They differ in sensor physics, spatial scale, and in whether they passively image a pre-existing vortex configuration or actively perturb it.

Modality	Primary signal	Reported operating scale
Scanning SQUID microscopy	Local flux $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$0 via a flux-locked loop	Pickup loop $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$1, $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$2, flux sensitivity on the order of $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$3 (Wells et al., 2018)
Scanning SQUID susceptometry	In-phase and out-of-phase mutual inductance, $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$4 and $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$5	$B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$6, $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$7, $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$8, scanning resolution $B(r) = (\Phi_0/2\pi\lambda^2)\,K_0(r/\lambda),$9 (Horn et al., 2023)
Scanning NV magnetometry	ODMR shift of a single NV center	Standoff $\lambda$0–$\lambda$1, pixel spacing $\lambda$2, spatial resolution $\lambda$3, acquisition time $\lambda$4–$\lambda$5 per $\lambda$6 map (Schäfermeier et al., 13 Feb 2026)
NV-based scanning quantum microscope in 2D NbSe$\lambda$7	cw-ODMR and Hahn-echo decoherence	Stand-off of order $\lambda$8–$\lambda$9, spatial resolution down to $K_0$0 (Jayaram et al., 5 May 2025)
STM/S vortex imaging	Conductance contrast or $K_0$1 spectroscopy	$K_0$2 operation, lateral resolution $K_0$3 per pixel, $K_0$4–$K_0$5 per map (0903.2389)
MFM-based SQVM	Frequency or phase response while a single vortex is trapped and dragged	Lift height $K_0$6–$K_0$7, resolution $K_0$8 comparable to the coherence length (Aladyshkin et al., 7 Jul 2025)
Magnonic-vortex quantum cavity proposal	Gyrotropic resonance shift or broadening for scanning EPR	Disk radius $K_0$9–$$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$0, $$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$1–$$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$2 (González-Gutiérrez et al., 2024)

The sensor physics is correspondingly heterogeneous. For NV-based methods, the ground-state spin Hamiltonian is

$$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$3

with $$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$4 and $$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$5, so the local field projection is obtained from $$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$6. Because the readout is frequency calibrated, $$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$7 is absolute without additional sensor calibration (Schäfermeier et al., 13 Feb 2026).

STM/S instead requires a cryogenic, vibration-isolated tunneling junction, often below $$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$8 and in fields up to several tesla, while MFM-based SQVM uses a Co/Cr-coated cantilever in lift mode and exploits the temperature window in which intrinsic pinning weakens sufficiently for the tip to mobilize a single vortex (Suderow et al., 2014, Aladyshkin et al., 7 Jul 2025).

3. Quantitative analysis of vortex order

A central contribution of scanning quantum-vortex microscopy is that it replaces qualitative vortex images with explicit structural metrics. In scanning SQUID microscopy on YBCO thin films, raw flux maps $$\nabla^2B - B/\lambda^2 = -\Phi_0\,\hat z\,\delta^2(r),$$9 were converted into vortex coordinates by locating local minima and fitting a 2D Gaussian to each spot. The measured spot width, $\psi(\mathbf r)$0, was much larger than the intrinsic $\psi(\mathbf r)$1 because of stray-field expansion at the probe height. Positional order was quantified through the two-point autocorrelation

$\psi(\mathbf r)$2

while orientational order was extracted from Delaunay triangulation and the hexatic order parameter

$\psi(\mathbf r)$3

A perfect triangular lattice gives sharp hexagonal rings in $\psi(\mathbf r)$4 and $\psi(\mathbf r)$5; the YBCO data instead yielded only a weak ring near $\psi(\mathbf r)$6 and $\psi(\mathbf r)$7, indicating isotropic orientational disorder (Wells et al., 2018).

STM-based studies of vortex matter employ related but not identical correlators. In SnMo$\psi(\mathbf r)$8S$\psi(\mathbf r)$9, Delaunay triangulation identified fivefold and sevenfold coordinated vortices, from which dislocations and disclinations were distinguished. Orientational order was quantified through

$\psi\to0$0

and positional order through

$\psi\to0$1

There, slowly decaying $\psi\to0$2 and power-law $\psi\to0$3 at $\psi\to0$4 identified a Bragg glass, whereas rapid decay of both at $\psi\to0$5 and $\psi\to0$6 established a vortex glass with short-range order (0903.2389).

In NV magnetometry, reciprocal-space analysis is especially direct because the field map is quantitative. The 2D Fourier transform $\psi\to0$7 reveals Bragg peaks whose radius $\psi\to0$8 sets the lattice spacing

$\psi\to0$9

and field

$d\ll\lambda_L$0

Sharp sixfold spots imply long-range order; azimuthal smearing indicates disorder or pinning. This procedure was used to verify the triangular lattice in BSCCO-2212 and to distinguish it from the diffuse-ring response of disordered YBCO thin films (Schäfermeier et al., 13 Feb 2026).

A related analysis was applied in few-layer NbSe$d\ll\lambda_L$1, where the autocorrelation of $d\ll\lambda_L$2 showed broad, smeared peaks with $d\ll\lambda_L$3 in the strongly disordered regime, again diagnosing short-range order rather than an Abrikosov lattice (Jayaram et al., 5 May 2025).

4. Material systems and observed vortex phases

Low-field YBCO thin films provide a clear case in which scanning SQUID microscopy resolved an isotropic vortex glass. After cooling in perpendicular fields from $d\ll\lambda_L$4 to $d\ll\lambda_L$5, the vortex ensemble exhibited only weak short-range positional correlations and negligible orientational order. The average spacing obeyed $d\ll\lambda_L$6, with $d\ll\lambda_L$7 at $d\ll\lambda_L$8. Above a critical field $d\ll\lambda_L$9, small clusters of $\Lambda = 2\lambda_L^2/d.$0–$\Lambda = 2\lambda_L^2/d.$1 vortices with nearest-neighbour spacing $\Lambda = 2\lambda_L^2/d.$2 appeared. For a random 2D gas, the probability $\Lambda = 2\lambda_L^2/d.$3 was estimated to be $\Lambda = 2\lambda_L^2/d.$4, yet experimentally $\Lambda = 2\lambda_L^2/d.$5–$\Lambda = 2\lambda_L^2/d.$6 of vortices participated in such clusters at $\Lambda = 2\lambda_L^2/d.$7. Fixed strong pinning centers were disfavored because the cluster locations were unreproducible upon thermal cycling (Wells et al., 2018).

In SnMo$\Lambda = 2\lambda_L^2/d.$8S$\Lambda = 2\lambda_L^2/d.$9, STM at $r\gg\xi$00 produced large-scale maps of about $r\gg\xi$01 vortices from $r\gg\xi$02 to $r\gg\xi$03. The $r\gg\xi$04 state retained quasi-long-range orientational and positional order and was classified as a Bragg glass, whereas the $r\gg\xi$05 and $r\gg\xi$06 states showed short-range order and topological disorder characteristic of a vortex glass. Combined with magnetisation and specific-heat measurements, these data supported a kinetic-glass description in which vortex topological disorder persists far below the peak-effect regime (0903.2389).

Cryogenic scanning NV magnetometry has established a quantitative comparison between ordered and disordered cuprate vortex matter. In BSCCO-2212 at $r\gg\xi$07, field-cooling at $r\gg\xi$08 produced a well-ordered triangular lattice in a $r\gg\xi$09 map. Twenty-six vortices were resolved, close to the expected $r\gg\xi$10, and the FFT ring at $r\gg\xi$11 gave $r\gg\xi$12 and $r\gg\xi$13, consistent with flux quantization. Under otherwise similar scanning conditions, a $r\gg\xi$14 YBCO film at $r\gg\xi$15 showed only nine vortices in a highly disordered arrangement and a diffuse reciprocal-space ring, consistent with stronger pinning (Schäfermeier et al., 13 Feb 2026).

Two-dimensional NbSe$r\gg\xi$16 extends the subject from static disorder to thermal evolution. In a few-layer sample with $r\gg\xi$17, vortices formed a strongly disordered vortex glass rather than a hexagonal lattice, and single-vortex scans showed lateral sizes of about $r\gg\xi$18, substantially exceeding the bulk core scale. As $r\gg\xi$19, the field profile broadened further, and controlled cooling from just below $r\gg\xi$20 to base temperature showed that rapid cooling yields weak, volatile vortex contrast whereas slow cooling produces stronger contrast and sharper local order, directly visualizing vortex-glass melting and freezing in a 2D system (Jayaram et al., 5 May 2025).

STM/STS in CsFe$r\gg\xi$21As$r\gg\xi$22 adds spectroscopic information at the individual-vortex level. Between $r\gg\xi$23 and $r\gg\xi$24, the lattice evolved from a distorted hexagonal arrangement to a distorted tetragonal one, with a mixed stripe-like region near $r\gg\xi$25. Spectra through a vortex center revealed a bound-state peak at $r\gg\xi$26, and exponential fitting of its spatial decay gave $r\gg\xi$27, consistent with the coherence length inferred from $r\gg\xi$28 (Yang et al., 2018).

5. Pinning, manipulation, and vortex dynamics

A significant branch of scanning quantum-vortex microscopy is explicitly dynamical. In Nb thin films of thickness $r\gg\xi$29–$r\gg\xi$30, MFM-based SQVM used the attractive interaction between a magnetic cantilever and a single vortex to map the pinning-force landscape. The tip first creates or captures a vortex, then drags it across the film when the temperature approaches $r\gg\xi$31 and the intrinsic pinning weakens. In $r\gg\xi$32 Nb at $r\gg\xi$33, the reconstructed pinning-force maps revealed a nano-network of pinning walls about $r\gg\xi$34 wide delimiting grain-sized cells about $r\gg\xi$35 across. Over $r\gg\xi$36 areas, the mean pinning force increased with thickness: $r\gg\xi$37 for $r\gg\xi$38, $r\gg\xi$39 for $r\gg\xi$40, and $r\gg\xi$41 for $r\gg\xi$42, consistent with a thickness-dependent granular network (Hovhannisyan et al., 2024).

This mode of operation differs fundamentally from conventional MFM of static vortices. The vortex is no longer merely an object being observed; it becomes the probe of the pinning potential. Reported spatial resolution reaches $r\gg\xi$43, comparable to the superconducting coherence length, and forward–backward scan mismatches can occur because the vortex may jump stochastically between nearby pinning sites. The resulting information is richer than static field imaging but is also intrinsically invasive, because the technique deliberately perturbs the vortex configuration (Aladyshkin et al., 7 Jul 2025).

Scanning SQUID susceptometry introduces a different active protocol by locally generating vortices with an AC field coil and detecting their response with a concentric pickup loop. In a niobium thin film near $r\gg\xi$44, a $r\gg\xi$45 sinusoidal drive with peak fields between $r\gg\xi$46 and $r\gg\xi$47 produced step-like changes in the reactive response $r\gg\xi$48 and sawtooth features in the dissipative response $r\gg\xi$49. These signatures were attributed to vortex–antivortex pairs nucleated by the local AC field, with individual steps corresponding to additional vortices trapped under the field coil after their antivortex partners had been ejected or pinned দূর away. Coupled London–Maxwell and time-dependent Ginzburg–Landau modeling then linked each branch of the measured $r\gg\xi$50–$r\gg\xi$51 hysteresis loop to a unique vortex configuration, enabling reconstruction of vortex trajectories with about $r\gg\xi$52 spatial resolution and millisecond temporal resolution (Horn et al., 2023).

NV-based scanning quantum microscopy accesses dynamics through spin decoherence rather than through direct transport or force measurements. In 2D NbSe$r\gg\xi$53, a Hahn-echo sequence

$r\gg\xi$54

was used to probe magnetic noise from vortex motion. Bringing the NV into contact with the superconducting surface reduced $r\gg\xi$55 relative to the lifted configuration, indicating MHz-range magnetic noise. The distance dependence of the induced decoherence gave a noise correlation length of about $r\gg\xi$56, and the measured increase of $r\gg\xi$57 upon lowering temperature below $r\gg\xi$58 was modeled by thermally activated vortex hopping combined with the temperature dependence of $r\gg\xi$59 (Jayaram et al., 5 May 2025).

6. Terminology, misconceptions, and extensions

The term “scanning quantum-vortex microscopy” is not attached to a single standardized instrument. In current literature it spans at least four experimentally realized superconducting modalities—SQUID, NV, STM/S, and MFM-based vortex dragging—and one theoretically developed spin-resonance architecture based on magnetic-vortex cavities. A common misconception is therefore that the term names one particular sensor platform. The more accurate description is a methodological class defined by its use of vortex physics as the primary imaging or sensing channel (Thiel et al., 2015, Suderow et al., 2014, Hovhannisyan et al., 2024).

A second misconception is that all forms are non-invasive. This is true for several implementations: quantitative NV magnetometry on YBCO explicitly reported no vortex displacement even when laser power was increased to $r\gg\xi$60, confirming negligible local heating in that experiment. By contrast, MFM-based SQVM intentionally traps and drags a vortex, and scanning SQUID susceptometry intentionally nucleates and drives vortex–antivortex pairs. Non-invasiveness is thus modality-dependent rather than definitional (Thiel et al., 2015, Aladyshkin et al., 7 Jul 2025, Horn et al., 2023).

The scope of the field also extends beyond superconducting vortices. Scanning NV microscopy has quantitatively imaged the stray field of a magnetic vortex core in a $r\gg\xi$61-diameter permalloy disk, unambiguously revealing the core and enabling direct comparison with micromagnetic simulations (Tetienne et al., 2013). Building on this broader magnetic-vortex context, a 2024 proposal described a scanning spin probe based on magnonic vortex quantum cavities: a sub-micron ferromagnetic disk whose vortex core provides a static field of about $r\gg\xi$62–$r\gg\xi$63 at $r\gg\xi$64–$r\gg\xi$65, gradients up to $r\gg\xi$66–$r\gg\xi$67, a circularly polarized rf field from gyrotropic motion at $r\gg\xi$68–$r\gg\xi$69, and an inductive readout via linewidth broadening. In that framework the resonance window has width of order $r\gg\xi$70–$r\gg\xi$71, suggesting an effective spatial resolution $r\gg\xi$72 and, for YIG with $r\gg\xi$73, the theoretical possibility of single-spin EPR detection in about $r\gg\xi$74 (González-Gutiérrez et al., 2024).

Open questions remain even in the superconducting case. The microscopic origin of the closely spaced vortex groups observed in low-field YBCO thin films remains unresolved; proposed causes include randomly enhanced local fields due to film geometry or demagnetization and freezing-in of a highly dilute vortex gas near $r\gg\xi$75, while fixed strong pinning centers were argued against by the lack of reproducibility upon thermal cycling. That study explicitly suggested that higher-resolution sensors with $r\gg\xi$76 or direct imaging such as scanning Hall probe microscopy could determine whether the group currents are genuinely connected at the film surface (Wells et al., 2018).

Taken together, these developments show that scanning quantum-vortex microscopy is best understood as a convergent research area rather than a single technique. Its unifying theme is the use of vortices—either as objects to be imaged or as active mesoscopic probes—to extract local superconducting parameters, classify ordered and glassy vortex matter, resolve vortex-core spectroscopy, map pinning landscapes, and quantify fluctuation-driven dynamics across length scales from tens of nanometers to tens of micrometers (Schäfermeier et al., 13 Feb 2026, 0903.2389).