Magnetic Susceptibility of Superconducting Films in the Meissner-London Regime
Abstract
Extracting the London penetration depth, λ(T), from the measured magnetic susceptibility, χ(T), of a superconducting platelet in a perpendicular magnetic field requires knowing the function χ(λ) computed for the actual sample geometry. Thin films are increasingly important in quantum computing and superconducting electronics, most of which operate at ultra-low excitation fields in the vortex-free regime. Understanding their electromagnetic response in the Meissner--London (finite-λ, no vortices) state is therefore essential for extracting intrinsic properties such as the London penetration depth and superfluid phase stiffness. While numerous prior works have addressed different aspects of this problem, here we focus on providing an analytical semi-empirical representation of χ(λ) based on the 2D sheet-current model and establishing the limits of its applicability using two independent full 3D numerical solutions: finite-element modeling and a volume-current integral equation. As a result, we provide χ(λ) for a disc of radius a and thickness d in closed form, valid for all λ and for aspect ratios δ≡ d/(2a)≤0.05, and accurate to 0.1%. The two-dimensional sheet description with Pearl length Λ=2λ2/d is adequate for δ≲10−4 and λ≳ 4d. Thicker films require the full closed form, which includes the three-dimensional corrections. Using the same numerical analysis, we provide an empirical formula that maps rectangular samples---the geometry most commonly used in experiments---onto a disc with an ``effective radius'' for subsequent calculation of χ(λ), and formulate a complete calibration protocol for whole-sample susceptibility measurements.
Related
- Effect of uniform bond-charge order on longitudinal spin-susceptibility jump at a superconducting instability inside a magnetic phase
- Magnetic susceptibility and superfluid density of inhomogeneous superconducting films: the Tc-mosaic model
- Electrochemical Growth of Full Volume Meissner Effect Superconducting BKBO
- Interplay of Meissner Currents and Magnetic Vortices in NbSe2 Probed by Tunneling Spectroscopy.
- Coupling of two individual magnon resonators via a superconducting microwave resonator operating in the strong coupling regime
BibTeX
@misc{prozorov2026magnetic,
title = {{Magnetic Susceptibility of Superconducting Films in the Meissner-London Regime}},
author = {Ruslan Prozorov},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.20455},
eprint = {2610.20455},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.20455},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Magnetic Susceptibility of Superconducting Films in the Meissner-London Regime
AU - Ruslan Prozorov
PY - 2026
DA - 2026/10/14
PB - arXiv
DO - 10.48550/arxiv.2610.20455
UR - https://arxiv.org/abs/2610.20455
AB - Extracting the London penetration depth, λ(T), from the measured magnetic susceptibility, χ(T), of a superconducting platelet in a perpendicular magnetic field requires knowing the function χ(λ) computed for the actual sample geometry. Thin films are increasingly important in quantum computing and superconducting electronics, most of which operate at ultra-low excitation fields in the vortex-free regime. Understanding their electromagnetic response in the Meissner--London (finite-λ, no vortices) state is therefore essential for extracting intrinsic properties such as the London penetration depth and superfluid phase stiffness. While numerous prior works have addressed different aspects of this problem, here we focus on providing an analytical semi-empirical representation of χ(λ) based on the 2D sheet-current model and establishing the limits of its applicability using two independent full 3D numerical solutions: finite-element modeling and a volume-current integral equation. As a result, we provide χ(λ) for a disc of radius a and thickness d in closed form, valid for all λ and for aspect ratios δ≡ d/(2a)≤0.05, and accurate to 0.1%. The two-dimensional sheet description with Pearl length Λ=2λ^2/d is adequate for δ≲10^−4 and λ≳ 4d. Thicker films require the full closed form, which includes the three-dimensional corrections. Using the same numerical analysis, we provide an empirical formula that maps rectangular samples---the geometry most commonly used in experiments---onto a disc with an ``effective radius'' for subsequent calculation of χ(λ), and formulate a complete calibration protocol for whole-sample susceptibility measurements.
ER -