Abstract
We describe a numerical approach to modeling magnetoelectric effects generated by spin-orbit coupling in inhomogeneous diffusive two-dimensional (2D) superconductors. It is based on direct minimization of the free energy of diffusion modes, including their coupling to the spin-orbit field strength, described by a recently discovered 𝜎-model action. We explain how to retain exact conservation laws in the the discretized model, and detail the numerical procedure. We apply the approach to the spin-galvanic and the inverse spin-galvanic effects in finite-size 2D superconductors, and describe short-range oscillations of circulating currents and spin densities originating from the spin-orbit coupling.
P. Virtanen. Magnetoelectric effects in diffusive two-dimensional superconductors studied by the nonlinear σ model. Physical Review B 111.2, 024510 (2025) https://doi.org/10.1103/PhysRevB.111.024510

