A finite element method for the quasiclassical theory of superconductivity

TL;DR

This paper introduces a finite element method based on discontinuous Galerkin approach for solving the quasiclassical theory of superconductivity, enabling detailed studies of complex device geometries and phenomena in superconductors.

Contribution

It presents a novel finite element method tailored for the quasiclassical superconductivity equations, allowing for self-consistent solutions in arbitrary geometries and material parameters.

Findings

Analyzed impurity effects on phase crystals in d-wave superconductors.

Modeled current flow and focusing in d-wave superconducting weak links.

Demonstrated the method's applicability to real experimental configurations.

Abstract

The Eilenberger-Larkin-Ovchinnikov-Eliashberg quasiclassical theory of superconductivity is a powerful method enabling studies of a wide range of equilibrium and non-equilibrium phenomena in conventional and unconventional superconductors. We introduce here a finite element method, based on a discontinuous Galerkin approach, to self-consistently solve the underlying transport equations for general device geometries, arbitrary mean free path and symmetry of the superconducting order parameter. We present results on i) the influence of scalar impurity scattering on phase crystals in $d$-wave superconducting grains at low temperatures and ii) the current flow and focusing in $d$-wave superconducting weak links, modeling recent experimental realizations of grooved high-temperature superconducting Dayem bridges. The high adaptability of this finite element method for quasiclassical theory paves the way for future investigations of superconducting devices and new physical phenomena in unconventional superconductors.