# Group Theory of Superconducting Order Parameters?

In crystalline superconductors, the order parameter $$\Delta(\mathbf{k})$$ (aka gap, or Cooper pair wavefunction) can be classified by its symmetry according to the representations of the symmetry group of the crystal. This can get complicated because pairing is between fermions which also have spin, and spin-orbit coupling also plays a role.

I am used to categorizing orbitals and vibrational modes of a point group by their representation from the chemistry point of view, but it seems the superconductivity literature has a very different understanding which is confusing for me. I have the following confusions/questions

1. The representative function for the odd-parity representations is said to be a "vector" quantity (see slides 6-9 here). What does this mean? All the textbook character tables give scalar polynomials instead (see here). To be even more explicit and show my confusion, the entry on slide 8 under $$A_{1u}$$ should be antisymmetric under a mirror operation along the $$z$$-axis (aka $$\sigma_h$$), but it is clearly not true for the vector function $$k_x \hat{\mathbf{x}}-k_y \hat{\mathbf{y}}$$ which doesn't even depend on $$z$$. What am I missing?

2. Superconducting order parameters are said to have no nodes (fully gapped), point nodes, or line nodes. As an example of point nodes, table 1 of this paper says an order parameter with $$B_{1u}$$ symmetry has point nodes. But in the character table for that group, $$B_{1u}$$ transforms as $$z$$, which means it has a whole plane of "nodes" when $$z=0$$, not just a single point. How do you get the nodal structure from the representation if not from the characteristic polynomial?

Can anyone clarify what's going on here? The understanding and notation of the superconducting order parameter in group theory seems to be very different than that of orbitals or vibrations.

First, to be clear about a few things:

• The spontaneously broken symmetry leading to superconductivity (no resistance, Meissner effect etc.) is $$U(1)$$-phase rotation symmetry, associated with conservation of (Cooper pair) particle number. Other symmetries, the crystal space group and spin-rotations, can be additionally broken. Your question pertains these other symmetries.
• The gap function $$\Delta(\bf{k})$$ by itself is a complex scalar function, i.e. an element of $$\mathbb{C}$$ for each $$\mathbf{k}$$. However, the transformation properties under space group and spin-rotation transformations allows the identification of additional structure. Most experiments measure the amplitude of this complex quantity.
• The total wavefunction must be antisymmetric. Therefore, if the spin-part is antisymmetric (singlet pairing) the orbital (space-group dependent part) must be symmetric, so even under parity. If the spin-part is symmetric (triplet pairing), the orbital part must be odd under parity.
• A parity/space inversion transformation simply replaces $$\mathbf{k} \to -\mathbf{k}$$.

• The scalar/vector denomination is about spin-rotation transformations.

A basis of states of the fermions is given by $$|k,\sigma>$$, where $$k$$ is the momentum and $$\sigma = \uparrow,\downarrow$$ the spin. The complex scalar function above can be decomposed as

$$\Delta(\mathbf{k}) = \sum_{\sigma,\sigma'} < k,\sigma | \Delta_{\mathbf{k},\sigma\sigma'} | k,\sigma' >$$ (Eq.1)

Here $$\Delta_{\mathbf{k},\sigma\sigma'}$$ (using the notation of the slides in your question) is a 2$$\times$$2-matrix for each point $$\mathbf{k}$$.

We can write out the 2$$\times$$2-matrix in the singlet-triplet basis. The singlet-state is one-dimensional, as such there is one degree of freedom to be specified by the $$\mathbf{k}$$-dependence. One can then call this the scalar wavefunction $$\psi(\mathbf{k})$$ (as in the slides), since the spin-state is already defined. The spin-singlet is antisymmetric, so the momentum-dependent part must be symmetric (even under parity).

Conversely, the triplet-basis is three-dimensional, meaning there are three degrees of freedom to be specified. Collecting the Pauli matrices in a vector $$\vec{\sigma} = (\sigma^x, \sigma^y, \sigma^z)$$, we can write the momentum dependence as $$\vec{d}(\mathbf{k}) \cdot \vec{\sigma}$$, which is a superposition of 2x2-matrices acting on the spin-subspace for each momentum $$\mathbf{k}$$. One can therefore call this a vector wavefunction $$\vec{d}(\mathbf{k})$$. It transforms as a vector under spin-rotations. Nevertheless, the gap function in (Eq.1) is a complex scalar.

• In fact, these are pseudoscalar and pseudovector functions, meaning they are even under parity transformations. The wavevector $$\mathbf{k}$$ is of course odd under parity.

If $$g$$ is an element of the transformation group, then the action on the $$d$$-vector is specified as

$$g \rightharpoonup \vec{d}(\mathbf{k}) = R_\mathrm{s}(g) \vec{d}(R_\mathrm{o}\mathbf{k})$$

Here $$R_\mathrm{o}$$ is the orbital representation and $$R_\mathrm{s}$$ the spin representation, which are odd and even under parity respectively. This notation is from slide 27 of this presentation by the same author as in your link.

To understand how the representative functions, like $$k_x \hat{x} - k_y \hat{y}$$, transform, recognize that the unit vectors are basis vectors of the $$d$$-vector, and therefore transform as pseudovectors, while the $$k$$-components transforms as vectors.

In Kaba and Sénéchal Table I, you can find explicit matrices for $$R_\mathrm{o}$$ (called $$U$$) and $$R_\mathrm{s}$$ (called $$R)$$. In particular, the transformation under $$\sigma_h$$ (called $$\sigma_z$$ in that reference) is specified as

$$R_\mathrm{s}(\sigma_h) = \begin{pmatrix} 1 & & \\ & 1 & \\ & & -1 \end{pmatrix}, \quad R_\mathrm{o}(\sigma_h) = \begin{pmatrix} -1 & & \\ & -1 & \\ & & 1 \end{pmatrix}$$

So $$\hat{x} \to \hat{x}$$ and $$k_x \to -k_x$$ gives a total minus sign, and the same for the other term, so $$k_x \hat{x} - k_y \hat{y}$$ transform with a minus sign under $$\sigma_\mathrm{h}$$. You can find the action under the other generators like this as well.

• The order parameter is the pair operator $$\hat{\Delta}_{\mathbf{k},\sigma\sigma'} = \hat{c}_{\mathbf{k},\sigma} \hat{c}_{-\mathbf{k},\sigma'}$$, where $$\hat{c}$$ is the electron annihilation operator. The symmetry group representation acting on $$\Delta$$ derives from the one acting on $$\hat{c}$$ as follows. Electrons are fermions, and $$\hat{c}$$ transforms under a spin or projective representation of the group (say $$D_{4\mathrm{h}})$$. This representation is not given in most character tables, but for instance the fourth power of the rotation generator is $$-1$$, not $$+1$$ in this representation. The representation acting on the order parameter follows from Clebsch-Gordan decomposition of the tensor product of this spin representation with itself, and includes only ordinary (unitary) representations. See the references by Sénéchal below.

• Nodes are points in $$k$$-space where $$\Delta(\mathbf{k}) =0$$. ($$\mathbf{k} =0$$ is excluded since pairing is always at some finite momentum, moreover near the Fermi surface.) Therefore, you simply need to inspect where the gap function vanishes.

Looking at the representative functions in Table I of the paper you cite, you can see that $$B_{1g}$$, $$B_{2g}$$, $$B_{3g}$$ vanish at $$k_z =0$$ for all $$k_x$$, $$k_y$$. Therefore, there is a line of nodes at $$k_z =0$$ as a function of $$k_x,k_y$$.

Conversely, the $$B_{1u}$$, $$B_{2u}$$, $$B_{3u}$$ representative functions only vanish at some very specific values of $$k_x, k_y, k_z$$, determined by the coefficients $$c_1$$, $$c_2$$. Therefore, these have point nodes. (By the way, there are errors/typos in the representative functions of $$B_{1u}$$ and $$B_{3u}$$ in that table. It should be $$c_1 k_y \hat{x} + c_2 k_x \hat{y}$$ for $$B_{1u}$$, and similar in $$y,z$$ for $$B_{3u}$$.)

The function $$z$$ transforms under $$B_{1u}$$ as you mention. But that is not the form of the gap function, which for instance is specified in momentum space. The gap function in question looks like $$c_1 k_y \hat{x} + c_2 k_x \hat{y}$$, and is independent of $$k_z$$. If you look in Table V of Kaba and Sénéchal, you can see that there exists a gap function with the form $$\hat{b}_z\hat{d}_0 k_z$$, also transforming under $$B_{1u}$$. However $$\hat{d}_0$$ is the singlet spin-state, and $$\hat{b}_z$$ is a particular form of the orbital (not spin) part of the order parameter.

References

J. Annett, Unconventional Superconductivity. Contemporary Physics, 36(6), 423 (2005). See especially section 3 for a discussion about the symmetry breaking starting from the full group $$G \times SU(2) \times U(1)$$ (crystal space group, spin-rotations, and phase-rotations) and a visual representation of the order parameter symmetries, and section 5 for symmetries of the gap function.

S. Sumita. Modern classification theory of superconducting gap nodes. PhD Thesis (2020). Kyoto University. Very thorough treatment. See section 1.1.1 for a short summary of the meaning of the symmetry of the order parameter.

Geilhufe & Balatsky. Phys. Rev. B 97, 024507 (2018). Another nice, concise treatment of transformation properties of the order parameter function, starting from Equation (4). Also includes odd-frequency pairing.

Ishizuka et al. Phys. Rev. Lett. 123, 217001 (2019) See Table 4(a) for the correct basis functions of the irreps of $$D_{2\mathrm{h}}$$.

S.-O. Kaba and D. Sénéchal. Phys. Rev. B 100, 214507 (2019) See section II.B, III.A for a elaborate derivation of the symmetry transformation properties of the $$D_{4\mathrm{h}}$$ order parameter, of not only the $$d$$-vector, but also the orbital part.

D. Sénéchal. Lecture notes A pedagogical treatment of the previous paper.

• +1 upvoted. Thank you so much for this answer it clarified a lot. I am still a bit confused about the symmetry of the k-space wavefunctions, particularly how the symmetry is deduced for the vector $d(\mathbf{k})$. Would you mind explicitly working out the symmetry operations on say the d-vector for $B_{1u}$ from the slide 8? Feb 18, 2021 at 9:04
• Also, you say that $\sigma_h$ is the inversion operator $\mathbf{r} \rightarrow -\mathbf{r}$, but that is not true, it is the mirror operation about plane of highest rotational symmetry. You can take a look at the character table of $D_{4h}$ where inversion $i$ and primary mirror $\sigma_h$ are clearly not the same. Given that, would you mind explicitly working out how $c_1 k_y \hat{x} + c_2 k_{x} \hat{y}$ transforms as $B_{1u}$? Do I need to incorporate the negative sign from flipping the spin or something? symmetry.jacobs-university.de/cgi-bin/… Feb 18, 2021 at 9:04
• You are correct about $\sigma_h$, some other slight errors too. The $d$-vector is only important in situations with spin-orbit coupling, meaning the rotations of spin are tied to the rotations of space. Then a group element $g$ acts on the $d$-vector as: $g \rightharpoonup \vec{d}(\mathbf{k}) = R_s (g) \vec{d}(R_o(g) \mathbf{k})$ Now $\vec{d}$ is a pseudovector, so it has the opposite sign under reflections than usual; so $R_s$ has the opposite sign from $R_o$ under reflections. I have some references too. If this is satisfactory, I'll update my answer. Feb 19, 2021 at 4:06
• Thank you Aron, yes I would very much appreciate if you can incorporate those changes (with an example calculation on the d-vector of B1u per your comment). Feb 19, 2021 at 7:57
• Okay, I changed it quite a bit, and I think it should answer all your questions now. Feb 21, 2021 at 5:38