Multiplicity, Parity and Angular Momentum of a Cooper Pair in Unconventional Superconductors of D 4 h Symmetry: Sr 2 RuO 4 and Fe-Pnictide Materials

: Sr 2 RuO 4 and Fe-pnictide superconductors belong to the same point group symmetry D 4 h . Many experimental data conﬁrm odd pairs in Sr 2 RuO 4 and even pairs in Fe-pnictides, but opposite conclusions also exist. Recent NMR results of Pustogow et al., which revealed even Cooper pairs in Sr 2 RuO 4 , require reconsideration of symmetry treatment of its SOP (superconducting order parameter). In the present work making use of the Mackey–Bradley theorem on symmetrized squares, a group theoretical investigation of possible pairing states in D 4 h symmetry is performed. It is obtained for I 4/ mmm , i.e., space group of Sr 2 RuO 4 , that triplet pairs with even spatial parts are possible in k z direction and in points M and Y . For the two latter cases pairing of equivalent electrons with nonzero total momentum is proposed. In P 4/ nmm space group of Fe- pnictides in point M , even and odd pairs are possible for singlet and triplet cases. It it shown that even and odd chiral states with angular momentum projection m = ± 1 have nodes in vertical planes, but E g is nodal , whereas E u is nodeless in the basal plane. It is also shown that the widely accepted assertion that the parity of angular momentum value is directly connected with the spatial parity of a pair is not valid in a space-group approach to the wavefunction of a Cooper pair.


Introduction
Pnictide superconductors [1] and Sr 2 RuO 4 [2] have been characterized in a great number of experiments [3,4] whose results are in some cases mutually excluding, and no consensus is currently achievable. The absence of changes in the Knight shift in the 17 O NMR spectrum on passing through T c confirms oddness of the SOP (superconducting order parameter) of Sr 2 RuO 4 [5]. Phase-sensitive measurements indicate that the phase of the SOP in Sr 2 RuO 4 changes by π under inversion and confirmed the odd nature of pairing [6]. Odd superconductivity in Sr 2 RuO 4 is also in agreement with the results of penetration measurements [7,8]. A slight increase in the spin susceptibility originating from the Ru-4d electrons in superconducting state of Sr 2 RuO 4 was connected with a triplet pairing [9]. However, a dramatic increase in T c from 1.5 K to 3.5 K and increase in H c2 ( by a factor of 20 ) in strained samples of Sr 2 RuO 4 [10] was connected with singlet phase.
A reduction in 17 O Knight shift observed for all strain values and temperatures T < T c [11,12] suggests even SOP and singlet pairs, which were represented by even chiral function: On the other hand, for the unstrained samples the reduction in Knight shift of approximately 50% is not inconsistent with helical states of A 1u,2u or B 1u,2u symmetry, which were written as [11] A 1u (B 1u ) : d(k) =xk x ±ŷk y (2) wherex andŷ stand for the components of triplet spin in the ESP (equal spin pairing) case, and k x and k y represent the continuous wavevector. However, by means of the NMR spectroscopy of normal state Sr 2 RuO 4 under uniaxial stress and corresponding density-functional calculations, it was demonstrated that there are two different effects associated with the strain-induced the enhancement of the DOS associated with the γband FS (Fermi surface ) passing the Y point of the BZ (Brillouin zone), and a substantial Stoner enhancement [13]. The enhancement of the DOS through the Stoner factor boosts ferromagnetic spin fluctuations, which favors triplet states and seems to disfavor singlet pairing [13]. Thus, the singlet or triplet nature of pairing in Sr 2 RuO 4 is still under debate. Nevertheless, there is no doubt about other unusual properties of this superconductor. The results of µSR (muon spin relaxation) experiments on superconducting Sr 2 RuO 4 indicate the presence of spontaneous internal magnetic fields, i.e., TRSB (time-reversal symmetry breaking) [14]. More resent µSR experiments established splitting between T c and temperature of TRSB, which rules out any mechanism based on interaction of magnetic fluctuations and conventional superconductivity [15]. Observed Kerr rotation below T c implies TRSB [16] and is consistent with the SOP of the form [17]: whereẑ stands for the component of triplet spin vector in OSP (opposite spin pairing case). The measurements of the specific heat resulted in lines of gap minima or zeros in planes [±1, ±1, 0] and [±1, ∓1, 0] [18]. High-resolution inelastic neutron scattering measurements and random phase approximation calculations suggest horizontal line nodes [19]. The results of angle-resolved study of ultrasonic attenuation in Sr 2 RuO 4 suggest that the gap function has nodes in vertical or horizontal planes [20]. Thermal conductivity measurements on single crystals of Sr 2 RuO 4 revealed vertical line nodes of SOP [21]. Recent results of sound velocity experiments are compatable with two-component d-wave SOP [22]. The authors [20][21][22] noticed the inconsistency between chiral pwave or d-wave superconductivity and vertical line nodes. However, it was shown that group theoretical vertical line nodes are inherent to real and complex forms of E u SOP [23]. Recent field-angle-dependent specific-heat measurements [24] revealed the presence of horizontal line nodes in the gap consistent with SOP of (1) structure. Chiral even SOP nodeless in the basal plane was represented in polynomial form as Transport properties of the planar and corner Josephson junctions formed between Sr 2 RuO 4 and Nb revealed time-reversal invariant superconductivity, with most probably helical p-wave symmetry [25]. Using the results of resonant ultrasound spectroscopy measurements, the structure of SOP d x 2 −y 2 ± g xy(x 2 −y 2 ) was proposed [26]. Thus, the above-mentioned model representations of SOP in Sr 2 RuO 4 [11,17,22] reflect experimental features, but the inconsistency between their theoretical nodal structure [11,17] and experimental data on vertical line nodes in the chiral case [18,20,21] exists.
The electron band structure of pnictide superconductors is quite unusual resulting in some unconventional features of their SOP. According to the calculations the FS of pnictides consists of two electronic parts: α and β near Γ point and a hole part γ near M point(b/2, b/2, 0) [27]. The results of scanning tunneling microscopy on Fe(Se,Te) single crystals indicate that the sign is reversed between the hole and the electron FS pockets (s ± -wave), favoring the unconventional pairing mechanism associated with spin fluctuations [28]. ARPES (angle-resolved photoelectron spectroscopy) investigations of NdFeAsO 0.9 F 0.1 detected some angular fluctuations of gap at β FS, but the results do not exclude p-wave or d-wave pairing states with nodes of the gap [29]. ARPES of Ba 0.6 K 0.4 Fe 2 As 2 revealed nodeless superconducting gaps of different values, opening simultaneously at the bulk T c on all three observed FS sheets [30]. The results are consistent with the SOP of s-wave symmetry, although the authors cannot rule out the possibility of nontrivial relative phases between the pairing order parameters on the different FS sheets. ARPES with laser excitation showed that KFe 2 As 2 (T c = 3.4 K) is an s-wave superconductor with unusual multigap structure, i.e., internal and external FS have no zeroes, whereas octet point nodes are located on the middle FS at an angular distance 5 • from coordinate axes, and A 1g gap symmetry was suggested [31]. However, tunneling conductance spectra of stoichiometric FeSe crystalline films in their superconducting state revealed evidence for a gap function with nodal lines, manifesting symmetry breaking from C 4 to C 2 [32]. Summarizing experimental results of combinatorial microwave measurements, it was concluded that LiFeAs has nodeless gaps with moderate anisotropy and that BaFe 2 (As 1−x P x ) 2 has highly anisotropic gaps containing at least one line-nodal gap [33]. A constant Knight in 75 As in NMR investigation of a single crystal of superconducting LiFeAs was considered as a signs of triplet pairing [34]. In addition, the observed anomalous hysteresis at high fields is compatible with the existence of chiral gap wave functions similar to the forms (4) or (5) [35]. Making use of µSR, symmetry analysis and first principles band structure calculation chiral d-wave SOP was proposed for weakly correlated pnictide compound LaPt 3 P [36].
Both types of compounds Sr 2 RuO 4 and Fe-pnictides have the same D 4h point group symmetry, but belong to different space groups: I4/mmm (D 17 4h ,139) and P4/nmm (D 7 4h , 129), respectively, and the symmetry analysis of pair symmetry inside the BZ is the same for them. The BZs for these space groups are presented in Figure 1. Representations of experimental data (1)-(4) are based on phenomenological point group approach [17,37,38], in which possible SOP functions are obtained by reduction of spherical harmonics to a point group. According to Anderson ansatz [39] wavefunction of a singlet pair is even and that of a triplet pair is odd. It was shown, making use of the space-group approach to the wavefunction of a Cooper pair, that this direct relation between multiplicity and parity of a pair is violated in symmetrical points and directions in a BZ [40,41]. Since the spherical functions are taken for the representation of SOP in phenomenological approaches, it was assumed that singlet pairs have even value of angular momentum l and triplet pairs have odd value of angular momentum [42]. It should be noted that in a general case, for each IR (irreducible representation) D (l) of a rotation group R 3 there are even D (l+) and odd D (l−) extensions on the rotation group with space inversion R 3 + IR 3 [43]. Hence, the direct relation between the parity of angular momentum value and the spatial parity of pair function is a consequence of the basis function choice, but not of the symmetry requirements. Usually the accepted formula (4) for triplet chiral states is nodeless in vertical planes, but making use of the space-group approach the examples of p-wave SOP with group-theoretical vertical lines of nodes were obtained [23]. Recent results on even nature of the SOP in Sr 2 RuO 4 [11] require also investigations of vertical nodes of singlet chiral states. It was mentioned that the concept of s-wave, p-wave, etc., has to be taken with a grain of salt because of the presence of a crystalline lattice that breaks translational symmetry [44]. By virtue of the fact that there are many experimental results on the symmetry of SOP of Sr 2 RuO 4 and Fe-pnictides, which apparently contradict with each other and with model representation, its group theoretical analysis independent on the models is required.
In the present work the space group theory [45] for one-electron states, Anderson ansatz [39] and Mackey-Bradley theorem [46,47] for two-electron states, are used to investigate spatial structures of singlet and triplet pairs in distinct points of a BZ. Symmetrical points in a BZ are identified, where direct connection between multiplicity and parity of a pair function is violated. Even and odd chiral states with angular momentum projection m = ±1 are built, and their nodal structure is investigated.

Symmetry of Two-Electron States
One-electron states in a crystal with symmetry group G are labeled by the wavevector k, its symmetry group H (little group) and the index κ of small IR D κk of H [45] . For brevity the second index for one-electron states will be dropped. The IR of the space group is an induced representation D κ ↑ G [45]. In a general point of a BZ the dimension of IR D κ ↑ G equals the number n of point group elements ofĜ (central extension of G). Two-electron space is a Kronecker square of this space, and its dimension in a general k− point equals to n 2 . This space can be easily decomposed into physically different parts making use of double coset decomposition of G with respect to H [46,47]: The double coset representatives d δ denote different terms in a Kronecker square D κ ↑ G ×D κ ↑ G decomposition. The notation × for direct (Kronecker) product is used throughout. For each double coset δ a representation P κ± δ is considered, whose character is written as: where and where the coset representative a is chosen from the relation: According to the Mackey-Bradley theorem [46,47] symmetrized (superscript +) and antisymmetrized (superscript −) parts of a Kronecker square are written as: The symmetrization (antisymmetrization) in the first item is performed by standard point group technique, sum in the second item runs over self-inverse double cosets, and the sum in the third item runs over not self-inverse double cosets, i.e., if Hd β H = Hd −1 β H. Total momentum of a pair is written as: where b is a reciprocal lattice vector. In the case of Cooper pairs the double coset representative d α is a space inversion I, and in some symmetrical points on the surface of a BZ d α is an identity element E. It should be noted that since in a space-group approach one-electron wavefunction is multidimensional, there are states of two equivalent electrons with nonzero total momentum K. Some such states will be considered later on. According to the Pauli exclusion principle, the symmetrized square of a spatial part of the wavefunction corresponds to a singlet pair and antisymmetrized square corresponds to a triplet pair. For K = 0 the induced representation P κ± α ↑Ĝ is a reducible representation of a point groupĜ . The frequency f q κ± of appearance of any IR Γ q ofĜ in the decomposition of induced representation P κ± α ↑Ĝ is given by Frobenius reciprocity theorem [45], i.e., by the formula: Inside the BZ in the majority of cases (but not in all) the spatial part of a singlet pair is even and spatial part of a triplet pair is odd. If IR D κ is two-dimensional, one obtains from (7)-(9) that χ P κ± I (E) = 4, χ P κ+ I (I) = 2 and χ P κ− I (I) = −2. Hence, it follows that in this case even and odd IRs are mixed in the decomposition of symmetrized and antisymmetrized squares. The mixing of even and odd IRs takes place also in nonsymmorphic space groups on the surface of a BZ [40,41,48,49] due to non-trivial factor-system [41,48]. In the case of group P4/nmm the wavevector group H in point M is a D 4h group (see Figure 1b)), and IRs are projective [50]. In the case of projective IRs symmetrized (antisymemtrized) characters in the first item in r.h.s. of formula (11) are written as [41,49]: where ω(g, g) stands for the factor system of projective representation. In addition, it may be derived from formula (11) that a usually accepted assertion that triplet (singlet) pairs have odd (even) angular momentum value is violated. For example, consider O h group and k a general point in a BZ. In this case all even (odd) IRs appear in the decomposition for singlet (triplet) pairs with frequencies equal to their dimensions. Consider IR T 1g for a singlet pair and T 2u and E u for triplet pairs. It is known that the even spatial part of the wavefunction of atomic term p 2 ( 3 P) transforms as T 1g in O h symmetry. In addition, the odd spatial part of atomic term p 3 ( 2 D) splits in crystal field into T 2u and E u . Hence, it follows that the angular momentum l = 1 should be assigned to a singlet pair of T 1g symmetry, and the angular momentum l = 2 should be assigned to eventual pairs of T 2u and E u symmetry.
The decompositions of symmetrized and antisymmetrized squares of IRs of the space group I4/mmm are presented in Table 1. The notations of points in a BZ are according to Figure 1. The results for points inside the BZs for groups I4/mmm and P4/nmm are the same and differ on their surfaces. The important point M for group P4/nmm is also considered. It is seen in Table 1 that in direction Λ, even IR A 2g appears in the decomposition of antisymmetrized square of E ↑ G. Additionally, A 1u appears in the decomposition of symmetrized square of E ↑ G. Hence, it follows that in this case triplet pair with spatial part A 2g and singlet pair with spatial part A 1u are possible. In point Z, where H = D 4h , double coset defined by the identity element results in zero total momentum of a pair. In this case, spatial parts of singlet and triplet pairs are even. However, since Sr 2 RuO 4 is almost a two-dimensional superconductor, these examples unlikely correspond to pairing in it, and we consider symmetry points in basal plane. Because the space inversion is an element of group D 2h of the wavevector M of I4/mmm BZ, the left coset defined by the space inversion is absent. Since all IRs of D 2h are onedimensional, in the left coset defined by the identity element antisymmetrized square vanishes, and only singlet pairs are possible. Hence, in this case triplet Cooper pairs are forbidden by symmetry. Table 1. Possible IRs of singlet and triplet pairs for group I4/mmm (D 17 4h ) depending on the wavevector in a BZ (see Figure 1a) and IR of wavevector group H. Double coset representatives and the pair's wavevectors K are presented in last two columns. The results for point M of group P4/nmm (D 7 4h ) are also included. Since, in a general point of a BZ all even IRs are possible for singlet pairs and all odd IRs are possible for triplet pairs, the absence of any IR on line (plane) corresponds to a point node (line of nodes).
However, in point M of I4/mmm BZ nontrivial triplet and singlet pairs are possible. Let us consider double coset defined by the rotation C 2y (see Figure 2a). In this case the relation (12) results in a two-electron wavevector K = (0, π, 0), whose group is D 4h . In this case IRs of singlet and triplet pairs are different, but they all are even. However, in this case translations by some lattice vectors result in a phase factor −1, and this state does not correspond to the definition of a Cooper pair completely. Nevertheless, since many other symmetries in topological superconductors are broken, this unusual case may be considered in the theories of unconventional superconductivity.
In point Y IRs and double coset, defined by the space inversion are the same as for direction Σ (see Table 1), consequently possible pair functions are the same also. Though there is a non-trivial wavevector relation, corresponding to double coset, defined by σ d (see Figure 2b). It is seen in Figure 2b), that in this case two-electron state in M has a star, consisting of two prongs. Since space inversion is an element of the wavevector group, this state is characterized by a parity. However, this state is also not completely symmetrical with respect to lattice translation, and some lattice translations result in a phase factor −1. The results for point M of group P4/nmm are also included in Table 1. In this case small IRs are projective [50], and formula (14) results in mixing of even and odd IRs in singlet and triplet cases (see also [49]).

Nodal Structure and Chirality of Cooper Pair Wavefunction
General group theoretical results of the previous section may be envisaged by introducing a Cooper pair basis set as follows. According to Anderson ansatz [39] for k a general point in a BZ spatial parts of singlet and triplet pairs are written by the two following formulas, respectively: and where r 1 and r 2 denote electron coordinates, and Ik represents the action of the space inversion on k. According to the space-group representation theory [45] the one-electron state in crystal is represented by the wavevector star {k}, consisting of Ĝ prongs. When acting by left coset representatives in the decomposition of Ĝ relative to C i on basis functions (15) or (16) one obtains total basis set for singlet or triplet Cooper pair, whose dimension equals to Ĝ /2. One can choose pure rotation subgroupĜ ofĜ or any other subgroup of index 2, which contains only one of the two elements connected by the space inversion. Singlet (triplet) functions form a basis of induced representation P κ± I ↑Ĝ of point groupĜ. The projection of this reducible basis, sets on the rows of IR Γ q of the whole group may be obtained making use of the standard projection operator technique [43], i.e., making use of formula: In the case of D 4h group there are eight basis pair functions for the singlet case and eight functions for the triplet case. Taking three rotations about the z axis (C 1 = C 4z , C 2 = C 2z and C 3 =C 3 4z ) and four reflections in vertical planes (σ 1 = σ y , σ 2 = σx y , σ 3 = σ x , and σ 4 = σ xy ) as left coset representatives, one obtains reducible Cooper pair basis set for a general k-point in a Brillouin zone, shown in Figure 3. Basis functions ψ s i or ψ t i are obtained by the action of left cost representatives h i on the functions ψ s or ψ t (see formulas (15) and (16)). Note that in Figure 3 all vectors have the same nonzero k z projection.  (15) and (16) ). Note that one-electron kvectors have nonzero k z component. Notation of elements : For one-dimensional IRs the linear combinations of the basis set are unique, and Table 1 defines their nodal structure completely. The nodal structure of two-dimensional IRs depends on additional quantum numbers and was considered for triplet pairs elsewhere [23]. It was also shown that chiral basis functions for triplet pairs constructed making use complex form of IR E u [50] have vertical nodes. However, chiral even functions were not considered in a framework of a space-group approach, and one of the aims of the present work is to consider such states. On vertical planes each one of the two-dimensional IRs appear as ones. It follows from formula (13) that for k a general point in a BZ E u and E g appear twice for triplet and singlet cases, respectively. Hence, it follows that for each vertical plane, two-dimensional IRs may be nodal or nodeless. Chiral forms of IRs E g and E u on C 4 subgroup [50] are presented in Table 2. Since IRs are characterized by parity, their matrices for all C 4h groups can be easily obtained. It can be easily verified that direct product A 2g × E g(u) equals to E g(u) .Consequently, there are two extensions of these matrices on the left coset defined by σ 1 , labeled by the sign of matrix of the left coset representative as a superscript. The matrices of element σ 1 with corresponding superscripts are also presented in Table 2. Table 2. Matrices D of IRs E g and E u for the rotation subgroup C 4 elements (C 1 = C 4z , C 2 = C 2z , C 3 = C 3 4z ) and two forms of matrix for the left coset representative σ 1 = σ y .
The pairs' functions of E g and E u symmetry constructed from the basis set shown in Figure 3 making use of projection operators (17) are presented in Table 3. It is seen from Table 3 that the basis functions ψ C 1 , ψ C 2 and ψ C 3 of the first row of E g(u) , obtained by the action of π/2, π and 3π/2 rotations on the initial function ψ s(t) , acquire phase factors i, −1 and −i, respectively, i.e., correspond to magnetic quantum number m = 1. Basis functions of the second row, obtained by the action of the same rotations on the function ψ σ 1 , acquire phase factors −i, −1 and i and correspond to magnetic quantum number m = −1. In the Ginzburg-Landau theory the SOP is identified with the wavefunction of a Cooper pair, which includes all pairs in the condensate. In the space-group approach this general pair function is obtained when the electron wavevector k runs in the sector 0 ≤ θ ≤ π/4 (see Figure 3), and other sectors are brought about by left coset representatives. Suppose, that there is a phase winding e iµθ in the first sector. It is natural to suppose that due to topological reasons µ = m. In this case the functions in the sectors labeled by odd numbers in Figure 3 may be considered as possessing angular momentum projection m = 1. It is clear that reflection σ 1 changes the direction of phase winding in sector 8 relative to sector 1, and the functions of even sectors possess angular momentum projection m = −1. In this technique the basis functions and their phases in the other sectors are obtained group theoretically by the action of left coset representatives. Table 3. Pair functions for two-dimensional IRs E g and E u of group D 4h . ψ s(t) is a pair function corresponding to Anderson ansatz [39] in a basis domain of a BZ (see formulas (15) and (16) for singlet and triplet cases, respectively). The other functions are obtained by the action of point group elements on ψ s(t) (see Figure 3). The superscript corresponds to the sign of matrix of left coset representative σ 1 (see Table 2). Basis functions of two rows of two-dimensional IRs are localized at different sectors in k-space and meet at vertical planes, where interference of pair functions occurs. In order to envisage such interference, the wave functions of pairs in a finite number of points were represented by normalized sum of real and imaginary Gaussians. Within the sectors the phases at neighboring points almost coincide, and constructive interference takes place. At the boundaries of the sectors, both constructive and destructive interference of the wavefunctions is possible. When squaring modulus of the complex function at every point we obtain the structure of SOP in the plane normal to the k z -axis. Nodal structures of IRs E + g and E + u , obtained by this numerical method presented in Figure 4, have nodes in plane [100] and deeps in diagonal planes. Figure 5 shows the same for E − g and E − u , where line of nodes is in plane [010]. Nodes in vertical coordinate planes are defined by additional quantum number, i.e., the sign of matrix σ 1 (see Table 2). Nodes and deeps in diagonal vertical planes are the result of phase differences between two complex functions. In order to envisage the origin of these peculiarities, the phases of basis functions with m = 1 and m = −1 for E + g(u) and E − g(u) are presented in Figures 6 and 7, respectively. Due to the same choice of left coset representative (denoted by superscript "+" or "−") the results for E g and E u with the same superscripts coinside, and subscripts are dropped in Figures 6 and 7. In Figure 6 near the positive direction of xaxis the phases of functions E + 1 at θ = δ and E + −1 at θ = −δ are the same (δ denotes small positive angle) and the phase −2π was added to the function E + −1 for convenience. On a clockwise path from θ =360 • − δ to θ = 90 • the phase winding of E + −1 function is 3/2π, and the total phase is −π/2. The phase winding on the counter-clockwise of the function E + 1 on the path from θ = δ to θ = 90 • is π/2. Thus, the total phase difference is π, and interference is destructive at θ = 90 • . Similar estimation at θ =45 • results in phases of the first and second function π/4 and −π/4, respectively, resulting in constructive interference of real parts, destructive interference of imaginary parts and deep (not node) of SOP. Figure 7 shows the E − case in which reflection σ 1 changes the phase, and the interference pattern is rotated by 90 • . Note that a similar interference pattern was obtained for complex odd IR E u without phase winding [23].  Table 3) Note that E + g is nodal and E + u is nodeless in basal plane.  Table 3). Note that E − g is nodal and E − u is nodeless in basal plane.  Figure 6. Phases of two basis functions of IR E + g (u) . Superscript + corresponds to the superscript of matrix D + (σ 1 ) (see Table 2) Subscripts in figure correspond to magnetic quantum numbers. . Superscript − corresponds to the superscript of matrix D − (σ 1 ) (see Table 2) Subscripts in figure correspond to magnetic quantum numbers.

IR 1) Pair Function 2) Nodal Panes
Thus, even and odd two-dimensional complex IRs have one vertical nodal plane [100] or [010]. The nodal and nodeless planes are interchanged when multiplying by characters of IR A 2g , so the vertical nodal planes of E g and E u are topologically unstable according to the definition of Kobayashi et al. [51]. The basis functions of E g and E u symmetry in chiral representations have angular momentum projection m = ±1. The essential difference between the structures of E g and E u appears in the basal plane (see Table 1), where E u is nodeless, and nodes of E g are topologically stable.

Weak Spin-Orbit Coupling
In a weak spin orbit coupling case total pair function is a direct product of its spin and spatial parts. In D 4h symmetry OSP spin functionẑ belongs to IR A 2g . Multiplication of E u by A 2g results in E u also. Since characters of A 2g equal −1 for all elements of the left coset with respect to C 4h , it follows that such a multiplication results in an interchange of E + u and E − u cases. The ESP spin functionsx andŷ belong to E g in its standard form. In chiral representation of E g [50] its basis triplet spin functions are linear combinations Decomposition of direct product E g × E u results in the states A 1u + B 1u + A 2u + B 2u , which are similar to helical states (2), (3). Multiplication of the E g spin part by the spatial part belonging to any odd one-dimensional IR results in E u , which corresponds to chiral ESP state. In the case of Sr 2 RuO 4 singlet spin state of A 1g symmetry and triplet OSP state of A 2g symmetry are considered. Thus, in D 4h symmetry total functions of chiral singlet and triplet OSP pairs belong to IRs E g and E u , respectively.

Discussion
Group theoretical basis functions with topological phase winding of the present work are similar to commonly used chiral functions (1), (4) and (5), and it is important to establish the reasons for the appearance of significant differences, namely vertical line nodes and nonunitary structure. It should be noted that authors of these models did not consider matrices, whose basis is formed by these functions. However, to define nodal structures of multi-dimensional IRs, additional quantum numbers, specifying the matrix, are required [52]. Consider the application of the technique of the present work to the functions (1). It is immediately verified that reflection σ y interchanges two functions (1). The function d zx + id zy is a good basis function for phase winding in the first sector with m = 1. In addition, the action of rotational elements of C 4z on this function result in the same phase factors as elements of the first row and first column of E g in chiral form [50] (see also Table 2). Continuing to apply the method of projection operators, we pass to the element σ y , which transforms the function d zx + id zy into d zx − id zy and interchanges rows of IR (see Table 2). Thus, we obtain the function with m = −1 in the 8 sector. Hence, it follows that application of the standard projection operator technique [43] and chiral form of IR E g [50] to the first of functions (1) results in the same phase structure as shown in Figures 6 and 7 and, consequently, the same nodal structure, as obtained in the present work in a space-group basis set. It should be noted that according to the Mackey-Bradley theorem, two-dimensional IRs appear twice in a total two-electron basis decomposition. It is immediately verified that multiplication of E g by B 1g results in E g also, but the magnetic quantum numbers are swapped. Thus, starting projection on the first row from d zx − id zy we obtain that functions with m = −1 are in the first row of matrix E g and in odd sectors in k-space (see Figure 2). In this case a total phase structure is opposite to that shown in Figures 6 and 7. Hence, it follows that the sum of these two basis sets results in nodeless and unitary (i.e., probabilities of m = 1 and m = −1 are equal in all kspace) SOP structures. However, since experimental results indicate that SOP of Sr 2 RuO 4 is chiral, nodal and nonunitary, we may conclude that this case corresponds to one chiral two-dimensional IR considered in the present work. The structure (4) is similar to (1) , but the chiral model form (5) , which is also written as B 1g + iB 2g , consists of two different IRs and has no group theoretic analog.
Since the primary detected odd spatial parity [5,6] of unconventional superconductor Sr 2 RuO 4 is under debate [11,13], the establishing of a connection between parity and multiplicity is very essential. In the present work, making use of Anderson ansatz [39] and the Mackey-Bradley theorem [46,47], possible pair symmetries on the dependence of the wavevector k in a BZ are obtained for symmetry group D 4h . In the majority of cases the wavefunction of a singlet pair is even, and that of a triplet pair is odd (see Table 1). Since in the space-group theory one-electron states are characterized by a wavevector star [45,47], possible pair states with the certain parity may also include states with nonzero total momentum K. In point M for I4/mmm group, in the case of pairing of electrons k M and C 2y k M , possible symmetries of singlet pairs are A 1g + B 1g and of triplet pairs are A 2g + B 2g . The wavevector K of these states equals (0, π, 0). These states are not completely translation invariant, i.e., some lattice translations result in a phase factor 1, but some translations result in a phase factor −1. In addition, for k in point Y there are pair states denoted as A g ↑ G and B 2u ↑ G for singlet case and B 1g ↑ G and B 3u ↑ G for triplet case. These pairs are not translation invariant, but characterized by parity. In addition to the symmetries broken in unconventional superconductors such two-electron states represent another possibility, i.e., pairs states with partially broken translation symmetry. Thus, triplet pairs in I4/mmm symmetry, being odd in general points of a BZ, may exhibit even spatial part in some symmetrical points of a BZ. Hence, it follows that recent experiments [11] manifesting the contribution of even pairs in Sr 2 RuO 4 may be explained in a triplet pairing scenario. The even triplet pairing and odd singlet pairing cases are referred to as odd frequency pairing [53].
In point M for non-symmorphic space group P4/nmm in addition to the standard relation between multiplicity and parity of a Cooper pair, singlet odd pairs and triplet even pairs are possible. Chiral states constructed making use of the complex form of IRs E g(u) [50] have a nodal plane in one of the vertical coordinate planes and deeps in diagonal planes. The position of the vertical nodal plane depends on an additional quantum number-the sign of matrix for reflection σ y . In these states the basis functions connected by rotations about the k z axis acquire phase factors corresponding to magnetic quantum numbers m = ±1. When adding a continuous gauge phase e iθ in the range 0 ≤ θ ≤ π/4 one obtains total phase winding group theoretically. The transformations of E g and E u pairs by rotations of the C 4 subgroup are the same. This assertion is in agreement with the general conclusion of the present work, that parity and angular momentum of a pair are not connected with each other. For example, in O h symmetry singlet pairs of T 1g symmetry are possible, but it is also the symmetry of 3 P term of p 2 configuration with angular momentum l = 1. This result differs from the assertion of phenomenological approaches based on spherical functions, that even (odd) values of momentum correspond to even (odd) inversion symmetry. However the present result does not contradict with general group theory [43]. Phenomenological chiral functions (1) and (4) arise from spherical function in k-space and are nodeless in vertical planes. Group theoretical chiral functions have vertical lines of node and deeps and agree with experimental results [20][21][22]. Thus, the only robust difference between nodal structures of E g and E u states is nodal structure in the basal plane.

Conclusions
Making use of the Mackey-Bradley theorem on symmetrized squares and Anderson ansatz for the Cooper pair, group theoretical investigation of possible pairing states in D 4h symmetry was performed. It was obtained that for non-symmorphic group P4/nmm of Fe-pnictide superconductors at point M of a BZ, odd singlet and even triplet pairs are possible. It was also obtained that triplet pairs with an even spatial part are possible in k z direction and in points M and Y of a BZ for the space group I4/mmm. For the two latter cases, pairing of equivalent electrons with total pair momentum in symmetry points of BZ other than Γ is proposed. Thus, the recent results [11][12][13], which revealed even features of previously considered as odd SOP of Sr 2 RuO 4 , may be explained by the presence of triplet pairs with even spatial structure in some points of a BZ. It was shown that even and odd chiral states with angular momentum projection m = ±1 have nodes and deeps in vertical planes, which corresponds to the experimental data. The only robust difference between nodal structures of E g and E u takes place in the basal plane, where E g is nodal and E u is nodeless.