Abstract

It is shown that the application of theorems of induced representations method, namely, Frobenius reciprocity theorem, transitivity of induction theorem, and Mackey theorem on symmetrized squares, makes simplifying standard techniques in the theory of electron structure and constructing Cooper pair wavefunctions on the basis of one-electron solid-state wavefunctions possible. It is proved that the nodal structure of topological superconductors in the case of multidimensional irreducible representations is defined by additional quantum numbers. The technique is extended on projective representations in the case of nonsymmorphic space groups and examples of applications for topological superconductors UPt3 and Sr2RuO4 are considered.

1. Introduction

Induced representation method is a powerful group-theoretical technique for the systems with subgroups of symmetry, such as symmetrical clusters and crystals. There are two essential subgroups, namely, local symmetry group of an atom in a cluster or crystal and the wavevector group (little group) in a BZ (Brillouin zone). The induced representation method and its application to crystals and molecules have been described in two classical books of Bradley and Cracknell [1] and Altman [2]. On the basis of this theory, it was recognized that the standard technique for construction of normal vibrations [36] and SALCs (symmetry adapted linear combinations) in clusters [6, 7] can be significantly simplified [811]. The Mackey theorem on symmetrized squares [12] and its new form for solid-state wavefunctions of Bradley and Davis [13] made the application of Pauli exclusion principle to solid-state wavefunctions possible, easy, and straightforward [14]. It follows from the space-group approach to the wavefunction of a Cooper pair [14, 15] that generally accepted direct relations between multiplicity and parity of a Cooper pair [16] are violated on symmetry axis in a BZ and on surfaces of a BZ in the case of nonsymmorphic space groups. Calculation of Cooper pair functions making use of the space group approach for UPt3 and some other superconductors [1721] confirmed this statement. It was also shown that transitivity of induction theorem results in additional quantum numbers for nanoclusters [22] and Cooper pairs [2325].

In the present work, the advantages of induced representation method for normal vibrations and SALCs construction and for the investigations of Cooper pair symmetry in topological superconductors are reviewed and some calculations for topological superconductors are performed. The pair nodal structure dependence on additional quantum number of IR (irreducible representation) of group in - form is investigated.

2. Induced Representations

Consider a finite group and its subgroup . The whole symmetry group is decomposed into left cosets with respect to :

There are two physical cases in which this decomposition is of great use. The atoms in a symmetric nanostructure are located at points with local symmetry , which may coincide with the full group , or with any of its subgroup, including a subgroup consisting of an identity element only. The action of the left coset representatives in (1) on initial atom results in an orbit of this atom. The second case is the band theory of solids [1, 2], where subgroup is the wavevector group. The action of the left coset representatives on the wavevector results in a star of this vector. The number of atoms, forming the orbit of an atom and the number n of prongs in a star of the wavevector, is given by the relation

where the modulus sign stands for the number of elements of a group.

Consider IR of subgroup and its basis set , where stands for the dimension of . According to the theorem on induced representation [1, 2], a set of functions is invariant in the group and forms the basis of induced representation , given by the formula

where

In a general case, the induced representation (3) is reducible in and is decomposed into the sum of its IRs of as follows:

where stand for frequencies of appearance of IRs . For the decomposition of induced representation (3), one can apply Frobenius reciprocity theorem [27], according to which the frequency of IR of the whole group in the decomposition of induced representation equals the frequency of appearance of IR in the decomposition of IR subduced to :

Thus, according to Frobenius reciprocity theorem, the values of in the right-hand side of (5) are given by the formula

where stand for characters of representations.

3. Molecular Vibrations and SALCs

Consider normal vibration construction for symmetrical molecules. The displacements of an atom A1 from the equilibrium position are the components of the vector. Under the operations of the subgroup , these coordinates are transformed by IRs of this subgroup. The set of displacements of atoms forming an orbit of atom A1 must be transformed by the representation (reducible) of the entire group . The displacements of other atoms are obtained by the action of left coset representatives (see formula (1)) on the basis set of the first atom. Let the atom A1 be transformed into the atom by the element , and the same element transforms the basis into basis . Then, under the action of the elements of the symmetry group of this atom, , the displacements are transformed by the conjugate representation of this group. This means that under the action the whole symmetry group G basis set of displacements of an orbit is transformed as a basis of the induced representation . Note that its matrix depends on the choice of left coset representatives, but its decomposition into IRs of does not. To calculate the set of IRs in the decomposition of the whole matrix of displacements, one does not need to construct it. This decomposition can be done making use of the Frobenius reciprocity theorem (i.e., by summation of characters on the subgroup in formula (7)).

Consider CH4 molecule of symmetry. The basis set of atomic displacements is shown in Figure 1. The first atom H1 is in -axis and directions of the axis on it are the same as those on central atom. The basis sets on other atoms are obtained by C2 rotations about corresponding axis. Local symmetry group of atom H1 is . The stretching vibrations (-component) of atom A1 displacements belong to IR A1 of and bending vibrations (- and -components) belong to IR of . The decomposition of the induced representations can be easily done by Frobenius reciprocity theorem (7) and making use of characters of and of Td on , presented in Table 1. Making use of formula (7), one obtains that and . To obtain all vibration modes, one should add vibrations of central atom, which belong to IR T2, and subtract the representation of the whole molecule displacements T2 and rotations T1 (see [6]). Thus, one obtains total set of vibration modes of CH4 molecule . Hence, it follows that as opposed to standard methods [46], where a multidimensional character of total vibration representation is constructed explicitly, the induced representation method makes it possible to obtain the same results by the character’s analysis on local subgroup only.

The method of construction of SALCs of atomic wavefunctions in molecules is very similar to the above discussed technique. Consider -orbitals in octahedral complex, in which six atoms are in coordinate directions (see Figure 2). Directions of the axis of the first atom in the six-atom orbit are the same as those of central atom. The axes of the second and third centers are obtained by counterclockwise rotation by 120° and 240°, respectively, about axis (111) and the rest of basis set is obtained by the space inversion. To convert basis functions on the centers 4, 5, and 6 into standard form, one can rename and basis functions, but it is not essential in our general consideration. Functions and belong to IRs and , respectively, of local group of the first atom. These orbitals are called - and -orbitals, respectively. Under the action of the elements of group, these functions are transformed independently by induced representations and , whose dimensions are 6 and 12, respectively. The characters of these induced representations are the same as those constructed for - and -orbitals in Oh symmetry explicitly [7]. It follows from Frobenius reciprocity theorem that there is no need to construct these characters and it is sufficient to check the orthogonality of characters of IRs A1 and E of group C4v with characters of IRs of the whole group . These characters are presented in Table 2. It is seen from Table 2 that the character of IR A1 of group is not orthogonal to the characters of IRs , , and of . Thus, we obtain that . Similar analysis for IR E of results in the fact that -orbitals belong to IRs , , , and . These sets of orbitals are the same as those obtained for -and -orbitals, making use of explicit construction of six- and twelve- dimensional matrix [7]. We see that IR appears twice in the whole basis set constructed from p-orbitals of one orbit of atoms. These basis sets are labeled by physical quantum number (i.e., IR A1 or E of local group ). In the case when the atoms in symmetrical clusters are on the planes of symmetry, it is possible to distinguish basis sets of repeating IRs by the proper choice of intermediate group [22]. This technique is useful for Cooper pairs [23] and will be considered in the next section.

4. Additional Quantum Numbers for Induced Representations

When the atoms are on the planes of symmetry, SALCs contain repeating IRs [22]. To label repeating IRs, it is possible to use the transitivity property of induction, which means that the two ways of induction, from directly into and via an intermediate subgroup , result in the same IRs:

The decomposition of induced representation in the left-hand side of (8) into IRs of is written as follows:

Consider an alternative route of induction. After inducing into the intermediate subgroup , we may decompose the result into IRs of :

Similar decomposition is carried out after inducing each one of IRs into :

According to the transitivity of induction theorem (8), the frequencies are connected by the relation

If frequencies in two-step induction in (10) and (11) do not exceed unity, every repeating IR in (9) acquires unique additional quantum number — index of IR of intermediate subgroup. This technique will be used in the next sections for additional quantum numbers of Cooper pairs.

5. Space-Group Approach to the Wavefunction of a Cooper Pair

Following Ginzburg and Landau [28], we consider SOP (superconducting order parameter) to be identical with a Cooper pair wavefunction. Hereby, we can take advantage of general quantum mechanical rule for construction of two-electron state. Thus, the wavefunction of two-electron state is constructed as a direct (Kronecker) product of one-electron wavefunctions. If the electrons are equivalent, the requirement of antisymmetry of the total wavefunction makes partial reduction of the total two-electron space into singlet and triplet subspaces possible. Thus, making use of the Pauli exclusion principle, we obtain that spatial part of a singlet (triplet) pair belongs to symmetrized (antisymmetrized) Kronecker square of one-electron wavefunction. Irreducible representations of space groups are induced from wavevector group [13, 26] and are defined by formula (3). The method of the decomposition of Kronecker squares of induced representations of finite group into symmetrized and antisymmetrized parts was developed by Mackey [12] and applied to space groups by Bradley and Davies [13].

The structure of the Kronecker square of IRs of a space group depends on the structure of the space group relative to wavevector group . The double coset decomposition of relative to is written asFor each double coset representative in (13), the intersection subgroup is considered:The representation of group is defined by the following formula:

where . For each self-inverse double coset (i.e., ), there are two extensions of into extended intersection subgroup , whose characters are as follows:

where . Then, the symmetrized and antisymmetrized parts of the Kronecker square of the IR of the space group are given by the two following formulae, respectively:

The first items on the right-hand sides of (18) and (19) correspond to the double coset defined by the identity element. The summations in the second items in the right-hand side of (18) and (19) run over all self-inverse double cosets. The last summations in the right-hand side of (18) and (19) correspond to non-self-inverse double cosets (i.e., ). According to the Pauli exclusion principle, the symmetrized Kronecker square (equation (18)) defines the spatial part of the singlet state (antisymmetrized with respect to spin coordinates) and the antisymmetrized Kronecker square (equation (19)) defines the spatial part of the triplet state, which is symmetrical with respect to spin coordinates. It follows from the Mackey-Bradley theorem that possible symmetries of singlet and triplet states are different for identity and self-inverse double cosets and are the same for non-self-inverse double cosets. Since total momentum of a Cooper pair equals zero, the following relation should be fulfilled for the k-vector and double coset representative:

where is a vector of reciprocal lattice.

It is clear from (20) that if is inside the BZ, one can take the space inversion or rotation about the axis normal to . At the surface of a BZ where , double coset defined by identity element also results =in two-electron state with zero total momentum. Pairing of electrons with equal moments was proposed for pseudogap states in cuprate superconductors and was called Amperean pairing [29]. Thus double coset defined by identity element corresponds to Amperean pairing.

For , a general point of a Brillouin zone relation (20) is fulfilled for double coset defined by the space inversion and group consists of the identity element and the space inversion I only. It is seen from formulas (15)-(17) that in this case symmetrized square belongs to even IR of group and antisymmetrized square belongs to odd IR of . The decomposition of corresponding induced representations can be easily done making use of Frobenius reciprocity theorem (formula (7)). Thus, one obtains that each even (odd) IR of the whole group appears in symmetrized (antisymmetrized) square with the frequency equal to its dimension.

The space-group approach is an alternative to phenomenological or point group approach [3033] in which the pair function is expressed in terms of spherical functions, similar to spherical harmonics in crystal field theory [7]. In point group approach, spatial part of a singlet pair is even and spatial part of triplet pair is odd [16, 30, 31]. It can be easily shown that this relation is violated for two-dimensional small IRs on the lines of symmetry [14, 15, 21]. Indeed, it follows from formulas (15)-(17) that for two-dimensional small IRs and . Hence, it follows that neither is decomposed into even IRs only nor is decomposed into odd IRs only. It is also clear from formulas (15)-(17) than on lines of symmetry the set of possible IRs for a pair depends on the IR of one-electron state. The direct relation between multiplicity and parity is also violated for nonsymmorphic space groups on the surface of a BZ [14, 1820].

It follows from the Blount theorem [32] that there are no symmetry requirements for lines of nodes in triplet case. Also phenomenological functions for group in a triplet case are nodeless [33]. It was pointed out [3436] that a contradiction exists between experimental lines of nodes in Sr2RuO4 and nodeless model functions [33]. It will be shown in the next section that nodal and nodeless basis sets for two-dimensional IRs can be distinguished by additional quantum numbers.

6. Nodes on Symmetry Planes

In the symmetry plane, the group of a pair is and it follows from formulas (15)-(17) that spatial part of a singlet pair belongs to IR and that of triplet pair belongs to IR of .

In L-S coupling scheme, total pair wavefunction is a direct product of its spatial and spin parts. Symmetrized Kronecker square of spatial part is multiplied by singlet spin function , which belongs to IR , and the nodal structure of a singlet pair is the same as that of its spatial part. The triplet spin state consists of three components: , , and , where the subscript denotes spin projection on -axis. In the commonly accepted notation, stands for and and are linear combinations of and . In the component, the spins of electron in a pair are opposite and this state is called OSP (opposite spin pairing). In the and components, the electron spins in a pair are equal and these states (or their linear combinations and ) are called ESP (equal spin pairing). In the case of D4h symmetry, OSP spin belongs to IR and ESP total spin belongs to IR . Theorem on the direct product of IR and induced representation [1, 2] is written as

Hence, it follows that the symmetry of OSP pair on the plane is and the symmetry of ESP pair on the plane is . These characters are also presented in Table 3. Since these IRs are induced into the whole group (see formulas (18) and (19)), the final results may be easily obtained making use of Frobenius reciprocity theorem. The absence of any IR of in the right column results in the absence of some IRs of the whole group for the k-vector on this plane. The intersection of symmetry plane with Fermi surface results in a line. Thus, the absence of any IR on the plane results in a symmetry-protected line of nodes. Some symmetry protected lines of nodes can be identified in a singlet case and for spatial part in a triplet case making use of Table 3 and Frobenius reciprocity theorem. It is seen from Table 3 that if in a triplet case both ESP and OSP pairs are taken into account, all odd IRs of , namely, and , are possible on the planes and there are no symmetry-protected lines of nodes. This conclusion is in agreement with Blount theorem, according to which in a triplet case, symmetry-protected lines of nodes are “vanishingly improbable” [32]. On the other hand, in axial symmetry, OSP and ESP pairs due to interaction with crystal field will have different energies and only one type of pairing, namely, ESP or OSP, will take place. In the case of Sr2RuO4, it is assumed that pairing is of OSP type [3436]. Hence, it follows from Table 3 that if only one type of pairing (ESP or OSP) takes place, symmetry-protected lines of nodes are possible.

In order to consider two-dimensional IRs and in more detail, their characters on symmetry groups of planes are also presented in Table 3. It is immediately verified from the data of Table 3 that in the basal plane (001) IR is completely forbidden for singlet pairs and appears twice for spatial part of triplet pair. Hence, it follows that on basal plane in a singlet case there is a symmetry-protected line of nodes of IR . In vertical planes, IR appears ones in a triplet case and IR appears ones in a singlet case. Note that according to Mackey-Bradley theorem, at general point of a BZ two-dimensional, IRs appear twice in the decomposition of the complete basis set. Thus, on vertical planes, one two-dimensional IR is permitted and the other IR is forbidden. Hence, it follows that such lines of nodes are not symmetry protected in terms of IRs of the whole group. These lines of nodes can be classified in terms of IRs of an intermediate group, namely, symmetry group of the plane . It follows from formulae (15) and (17) that the symmetry of spatial part of a triplet pair on the plane corresponds to of group . It is immediately verified that when inducing from into , one obtains the set of IRs, which also includes . Hence, it follows that in general point of a BZ, one can label two basis sets of IR by IR of intermediate group . The basis set () is nodeless on the plane and the basis set () is nodal on the plane.

7. Nonsymmorphic Space Groups and Symmetry Lines

In the center of a BZ and in some points on its faces, space inversion belongs to the wavevector group . In this case, the double coset defined by identity element results in zero total momentum of two-electron state. In the center of a BZ, the first items in right-hand side (18), (19) are calculated by standard methods [37] and one obtains even spatial part only for triplet and singlet pair.

Consider as an example point in a space group P4/nmm () of Fe-pnictide superconductors [38]. In these materials, there two sheets of Fermi surface: one near point and one near point [38, 39]. In point , space inversion belongs to the wavevector group and two electrons with equal nonzero moments combine into a Cooper pair and its zero total momentum is due to translation periodicity of a crystal. Since IRs are projective at this BZ point of nonsymmorphic space group , projective factors appear in standard formulas for symmetrization (antisymmetrization) [37] for finite group IRs:It is seen from the results, presented in Table 4, that in a singlet pairing case, two even and one odd IRs are possible for each one of four projective IRs of group . In a triplet case, only one odd IR is possible for each IR of wavevector group.

The space-group approach to the wavefunction of a Cooper pair can be generalized for strong spin-orbit coupling case [15]. In this case, possible pair functions belong to antisymmetrized squares of double-valued representations. For , a general point of a BZ and on the planes of symmetry, where double-valued small IRs are one-dimensional, corepresentations of magnetic groups [1, 26] should be used [15]. Consider two planes (001) in a BZ of group, namely, plane and plane. In group, space group elements and are connected with translation , which equals to a half of Bravais lattice translation t3. [26]. Double-valued IRs [26] are presented in Table 5. It should be noted that since there are two phase factors for , namely, one due to nonsymmorphic structure of group and one due to the double-valued IRs of spin rotation, double valued IRs are real [26].

Consider modifications of formulas (15)-(17) due to nonsymmorphic structure of group. Multiplication rule for space-group elements is as follows [1, 2, 26]:

where stands for a point group element and for its nonprimitive translation. Making use of this rule, one obtains that the transformation by inversion results in point a phase factor in formula (15) for :

Also, there is a phase factor in formulas (16) and (17) for , since

One more phase factor appears, since small representations are double-valued and total phase factor for equals one. Possible IRs of Cooper pairs for these planes are presented in Table 6. For , only one even IR and all odd IRs and of group are possible. Hence, it follows that not all even IRs of appear in induced representation . Since all odd IRs of are possible, all odd IRs of are also possible. Thus, in agreement with Blount theorem, there are no symmetry-protected lines of nodes for odd IRs. It is also seen from Table 6 that for only one odd IR is possible and there are symmetry-protected lines of nodes of odd IRs. Similar results were obtained by Micklitz and Norman [1820]

In direction in a BZ for group, the wavevector group is and there are three two-dimensional double-valued IRs [26]. On the BZ face in direction 1/2 , symmetry group is C3v and there are two one-dimensional IRs and one two-dimensional double-valued IR. The results for these two symmetrical directions presented in Table 7 show that in highly symmetrical directions possible pair symmetries depend on the internal quantum number — label of IR of little group. Due to large number of rotational elements in group C6v, there is an analogy between IRs of this group and angular momentum jz and it was called jz dependence [21]. Present results for p3 in 1/2 b1 + b3 direction are the same as that for [21]. IRs p1 and p2 are one-dimensional and one can construct antisymmetrized square of p1 or p2 separately. These two states are connected by time-reversal and form a corepresentation . Antisymmetrized square of consists of the same IRs as the result [21] for .

8. Basis Functions

In a general point of a BZ, the wavevector group consists of an identity element only and the expressions for spatial parts of singlet and triplet pairs which follow from Mackey-Bradley theorem [13] may be written aswhere denotes a general point in a representation domain of a BZ [1]. All one-electron wavefunctions are obtained when runs over representation domain of a BZ and all elements of the central expansion of the space group G act on one-electron wavefunction [1, 2, 26]. The results of action of these elements on form a star of the wavevector. Thus, the number of one-electron basis functions of Cooper pair for a general point of a BZ equals . Pair functions may be constructed for all in the representation domain of a BZ making use of formulas (27) and (28). Since the space inversion is included in formulas (27) and (28), the number of elements in two-electron basis set of definite multiplicity reduces to .

We take a space group with central extension (i.e., symmetry group of unconventional superconductor S2RuO4 [3235]) as an example of application of the space-group technique for construction of pair’s function. The basis set for a general point in a BZ of a space group with point group is shown in Figure 3. Notations in [26] are used for the cubic point group elements. This makes the direct use of group multiplication table possible [26]. In these notations, is an identity element, , , and are rotations at , , and about the axis (001), , , , and are reflections on planes (100), (010), (-110), and (110), respectively. Notations in [26] for the remaining elements of group may be obtained by making use of the following multiplication rules for (space inversion):The labels of group elements hi are used as subscripts of basis functions in Figure 3; for example,In the case of induced representations, standard projection operator technique [37] may be reduced to the action of left coset representatives on the basis functions of subgroup [9] and one obtains the following formula for construction the wave function of a Cooper pair belonging to any IR of :The prime in (31) means that the sum includes only one of the two elements, connected by inversion.

Spatial parts of wavefunction of triplet Cooper pairs for IR of group in - representation are presented in Table 8, where it is assumed that wavevector runs over basis domain of a BZ. The results for one-dimensional IRs and for in chiral representation are presented elsewhere [25]. Nodal structure on the symmetry planes may be estimated as follows. Since pair function on the plane is invariant under plane reflection, the nodeless linear combination is unchanged under the action of the plane reflection. At the symmetry plane, mirror counterparts can merge or cancel. If the signs of two mirror counterparts are changed under the action of plane reflection, they will cancel at the plane. The intersection of a nodal plane with Fermi surface results in symmetry-protected nodal line of SOP.

The action of reflection in (010) plane on the first line of results in

The sign is unchanged and there are no nodes in this plane. It is seen from Figure 3 that when runs over the basis domain of a BZ, the functions of the first line of are around - axis and the functions of the second line are around -axis. Hence, we need to investigate nodes of the first function in plane (010) only and nodes of the second function in the plane (100) only. The action of reflection in (100) on the function of the second row results in

Hence, it follows that there are no nodes in plane (100) also. We will call the nodes due to cancelation of functions of one row as nodes of the first type.

In the case of two-dimensional IRs, there is one more possibility for lines of nodes, namely, cancelation between functions of the first and second rows. Reflection in plane (-110) does not change the sign of the sum of two basis functions of basis set:Thus, there are no nodes on plane (-110). The action of reflection in (110) plane changes the sign of basis set:

Hence, it follows that the basis set of IR has a nodal plane (110). We will call the nodes, resulting from the cancellation of basis functions belonging to different rows of two-dimensional IR, nodes of the second type.

It is immediately verified that multiplication of by all even one-dimensional IRs results also in . Thus, the other basis sets can be obtained by projection on IRs transformed by , , and . Basis sets of symmetry, labeled by these one-dimensional IRs, written as superscripts, are also presented in Table 8. These basis sets and the results of their nodal structure analysis are presented in Table 8. It is seen from Table 8 that basis sets and have nodes of the first type on two vertical coordinate planes, but basis sets and have no nodes of the first type. It is also seen from Table 8 that all IRs have one vertical nodal plane of the second type. Recall that it follows from Mackey-Bradley theorem that appears twice in complete basis set, but once of each vertical plane (see Table 3 and discussions in the text). Consider and , differing by nodal quantum number on the plane (100). It is seen from Table 8, that in agreement with the results of Table 3, there is one nodeless IR, and the other is nodal on each vertical plane. The same result is valid for other pair of IRs and . The authors of experimental works [3436] on topological superconductor Sr2RuO4 obtained vertical lines of nodes of odd triplet SOP and concluded that OSP pairing is the most probable. In group, the symmetry of OSP spin is . Since we investigated products of basis set by all one-dimensional IRs, our results include symmetry of OSP case and confirm that vertical lines of nodes are possible for OSP pairs of symmetry. Note that there are two different forms of IR , namely, real and complex. Real form, considered in present work, is usually used in quantum chemistry [7, 8]. Basis functions of complex form [26] have definite momentum projection on -axis and are considered elsewhere [25].

9. Conclusion

The advantages of applications of induced representation method to symmetrical molecules and Cooper pairs in solids are reviewed and some results for pairing in topological superconductors are obtained. It is shown that induced representation method makes it possible to simplify standard methods for normal vibrations and SALCs construction (i.e., list of permitted IRs is obtained by the analysis of character tables on local subgroup only). The applications of Mackey-Bradley theorem in the case of nonsymmorphic space groups and projective representations are further investigated and examples of UPt3 and Fe-pnictide superconductors are considered. Nodal structure SOP for symmetry in x-y representation is obtained and its dependence upon the addition quantum numbers is envisaged. These nodal structures are in qualitative agreement with experimental data for topological superconductor Sr2RuO4.

Data Availability

The references to the sources of the data are in the paper, and at present there are no other places where they are available.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

The work was carried out within the State Assignment on Fundamental Research to the Kurnakov Institute of General and Inorganic Chemistry.