Impurity strength–temperature phase diagram with phase crystals and competing time-reversal symmetry breaking states in nodal <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mi>d</mml:mi></mml:math>-wave superconductors

Phase crystals are a class of nonuniform superconducting ground states characterized by spontaneous phase gradients of the superconducting order parameter. These phase gradients nonlocally drive periodic currents and magnetic fields, thus breaking both time-reversal symmetry and continuous translational symmetry. The phase crystal instability is generally triggered by negative and inhomogeneous superfluid stiffness. Several scenarios have been identified that can realize phase crystals, especially flat bands at specific edges of unconventional nodal superconductors. Motivated by omnipresent disorder in all materials, we employ the <a:math xmlns:a="http://www.w3.org/1998/Math/MathML"><a:mi>t</a:mi></a:math>-matrix approach within the quasiclassical theory of superconductivity to study the emergence of phase crystals at edges of a nodal <b:math xmlns:b="http://www.w3.org/1998/Math/MathML"><b:mi>d</b:mi></b:math>-wave superconductor. We quantify the full phase diagram as a function of the impurity scattering energy and the temperature, with full self-consistency in the impurity self energies, the superconducting order parameter, and the vector potential. We find that the phase crystal survives even up to <c:math xmlns:c="http://www.w3.org/1998/Math/MathML"><c:mrow><c:mo>∼</c:mo><c:mn>40</c:mn><c:mo>–</c:mo><c:mn>50</c:mn><c:mo>%</c:mo></c:mrow></c:math> of the superconducting critical impurity strength in both the Born and unitary scattering limits. Finally, we show how mesoscopic finite-size effects induce a competition with a state still breaking time-reversal symmetry but with translationally invariant edge currents.