Vortex state of a d-wave superconductor at low temperatures

A systematic perturbation theory is developed to describe the magnetic field-induced subdominant s- and $d_{\mathrm{xy}}$-wave order parameters in the mixed state of a $d_{x^2-y^2}$-wave superconductor, enabling us to obtain, within weak-coupling BCS theory, analytic results for the free energy of a d-wave superconductor in an applied magnetic field $H_{c1\mathrm{}}\lesssim H\ll H_{c2\mathrm{}}$ from $T_c$ down to very low temperatures. Known results for a single isolated vortex in the Ginzburg-Landau regime are recovered, and the behavior at low temperatures for the subdominant component is shown to be qualitatively different. In the case of the subdominant $d_{\mathrm{xy}}$ pair component, superfluid velocity gradients and an orbital Zeeman effect are shown to compete in determining the vortex state, but for realistic field strengths the latter appears to be irrelevant. On this basis, we argue that recent predictions of a low-temperature phase transition in connection with recent thermal conductivity measurements are unlikely to be correct.