Symmetry analysis of transport properties in helical superconductor junctions

We study discrete symmetries satisfied by helical $p$-wave superconductors with d-vectors $k_{x}\hat{x}\pm k_{y}\hat{y}$ or $k_{y}\hat{x}\pm k_{x}\hat{y}$ and transformations brought by the symmetry operations to ferromagnet and spin-singlet superconductors, which show intimate associations with transport properties in heterojunctions including helical superconductor. Especially, the partial symmetries of the Hamiltonian under the spin-rotation and gauge-rotation operations are responsible for novel invariances of the conductance in tunnel junctions and new selection rules of the lowest current and peculiar phase diagrams in Josephson junctions which are reported recently. The symmetries of constructed free energies for Josephson junctions are also analyzed which are consistent with the results from Hamiltonian.


Introduction
Symmetries play an important role in the classification of topological matters [1,2,3,4,5,6] and the definition of topological invariants [7,8,9,10,11,12]. The spin-triplet superconductors (TSs) with d-vectors k xx ± k yŷ or k yx ± k xŷ belong to symmetry class DIII according to the tenfold classification of topological insulators and superconductors [13,14,15], which satisfy time-reversal symmetry (TRS), particlehole (or charge-conjugation) symmetry (PHS) and chiral symmetry (CS) but break SU(2) spin-rotation symmetry (SRS). The superconducting systems, possessing the non-trivial Z 2 topological invariant and supporting helical Majorana modes, are analogous to the quantum spin Hall state [16,17]. Such superconductors are called helical superconductors. The helical superconducting states are candidates for paring in Sr 2 RuO 4 [18,19,20] and the triplet part of the oder parameter in the noncentrosymmetric superconductor CePt 3 Si [21]. Among others, k xx + k yŷ is the twodimensional analog of the Balian-Werthamer state (B phase) in 3 He [22]. In contrast, the chiral TSs with d-vectors (k x ± k y )ẑ can be viewed as two copies of the spinless superconductor which belongs to the symmetry D class with the non-trivial Z topological number and a chiral Majorana mode [2,3,5]. The chiral superconducting states satisfy PHS and SRS about theẑ-axis but beak TRS, which are analogous to the quantum Hall state [23].
Researches on transport properties of topological junctions not only can provide methods for the detection of the Majorana fermion edge states [24,25,26,27,28,29,30,31] but also can probe the symmetries satisfied by the topological superconductors [32,33,34]. The charge conductance in the ferromagnet (F)|chiral TS junction is invariant when the magnetization in F is rotated about theẑ-aixs due to the SU(2) SRS kept by TS [35]. The Josephson effects in spin-singlet superconductor (SS)|F|chiral TS junctions show that the form of the lowest current (sin φ or cos φ) are strongly related to the symmetries of Cooper pair functions in SS [36]. The 0-π phase transition in the purely chiral superconductor junctions with different chiralities are also clarified [37]. These results are very different from those for the helical superconductor junctions since the SU(2) rotation symmetry for the latter is completely broken. The conductance in F|helical TS junctions exhibits higher anisotropy when the direction of magnetization is changed [38]. Both sin φ-type and cos φ-type of current in the SS|F|helical TS junctions always exists; the selection rules for the current can be given by the orientation of magnetization [39]. The purely helical Josephson junctions with different d-vectors can host more ground states such as π-phase, ϕ-phase and ϕ 0 -phase which can be tuned easily by the direction of magnetization [40]. Especially, the ϕ 0 -phase is recently receiving an increasing interest [25,41,42,43,44,45,46,47].
However, the SU(2) symmetry breaking in helical superconductor junctions does not mean the absence of symmetries of transport quantities. Actually, both the charge conductance and the Josephson current possess invariances under the rotation of magnetization [38]. The invariances have important effects to classify the SS|F|helical TS junctions into two types and to form the symmetric phase diagrams of the purely helical junctions [39,40]. In this paper, we consider the transport properties of the F|helical TS, SS|F|helical TS and helical TS|F|helical TS junctions from the viewpoint of symmetries of Hamiltonian in order to clarify the origin of the invariances of transport quantities. The main findings in this paper are as following.
(a) The Hamiltonian for helical superconductor is invariant under the joint operation of the SU(2) spin-rotation and the U(1) gauge transformation, which will lead to the invariances of conductance and Josephson current.
(b) The Hamiltonian satisfies partial symmetries under particular spin-rotation operation or the joint operation of rotation and gauge transformation; the operations can transform the BdG Hamiltonian from H(k x , k y ) into H(k x , −k y ). The partial symmetries are responsible for the symmetries of transport quantities and phase diagrams. Especially, if we take the different k y parities of superconducting wave functions in SSs into account, the partial symmetries will produce different selection rules of the lowest order current for different SS|F|helical TS junctions, which will divide the junctions into two types.
(c) TRS and PHS of the Hamiltonian bring about the sign reversal of Josephson current and the symmetry of the current for junctions without the non-magnetic interfacial potential, respectively.
(d) The symmetries of current derived from the Hamiltonian are consistent with those from the constructed free energies for the helical Josephson junctions, which directly reflects the mechanism of interaction between ferromagnetism and helical superconductivity.

Superconducting system
The generic Hamiltonian for superconductor is with Ψ k = (c k↑ , c k↓ , c + −k↑ , c + −k↓ ) T and the BdG Hamiltonian where c kα (c + kα ) is the annihilation (creation) operator of electron with momentum k = (k x , k y ) and spin α,ǫ(k) describes the normal dispersion of electron and∆(k) denotes the energy gap matrix. The normal state is invariant under the space inversion, i.e.,ǫ(k) =ǫ(−k). and p − yx -wave, the orbital angular momentum (l, l z ) for the spin-up part with S z = 1 and the spin-down part with S z = −1 is (1, −1) and (1, 1), respectively. (b) For the p − xy and p + yx -wave, the angular momentum (l, l z ) for the spin-up part with S z = 1 and the spin-down part with S z = −1 is (1, 1) and (1, −1), respectively. The z-axis has been chosen along the crystallographic c-axis.
uv -wave TS respectively in this paper which can be regarded as the superposition of the equal-spin wave functions with different orbital angular momentum, as shown in figure 1. Note, the phase before the orbital wave functions is different for p ± uv -wave states. Actually, the Hamiltonian for the p ± uv -wave states can be written in a block diagonal form in the basis Ψ = (c k↑ , c + −k↑ , c k↓ , c + −k↓ ) T which is The block diagonal Hamiltonian in the spin subspace is similar to the one discussed in [48] which describes the system consisting of two separate copies of spinless superconductors with the opposite chirality. That is the reason for our use of the term 'helical superconductor' in this paper. The helical superconductor possesses TRS. We introduce the time-reversal operator T = χ T K where χ T =τ 0 ⊗ iσ 2 withτ 0 the 2 × 2 unit matrix in particle-hole space and K the complex conjugation operator, which flips the sign of momentum and the direction of spin. The actions of T on the creation and the annihilation operators which is obviously satisfied by the BdG Hamiltonian. After the time-reversal operation, the Hamiltonian in equation (2) keeps invariant. However, if the superconductor has a phase φ which is an important quantity in the Josephson effect, the operation T will change the phase into −φ.
The superconductor also obeys PHS. We define the charge conjugation operator C = χ C K with χ C =τ 1 ⊗σ 0 whereτ i is Pauli matrices in particle-hole space andσ 0 the unit matrix in spin space. The action of C on the electron state | k ↑ will produce a hole in the state | −k ↓ . Consequently, the transformations of the creation and the annihilation operators are Cc kα C −1 = c + kα and Cc + kα C −1 = c kα , respectively. The symmetry with {C, H} = 0 requires which is also satisfied by the BdG Hamiltonian in equation (2). After the charge conjugation operation, the energy of quasiparticles becomes negative. In other words, ǫ(k) and∆(k) in equation (2) become −ǫ(−k) and ∆(−k), respectively. However, the phase φ still keeps its value.
Since the system satisfies both TRS and PHS, we can define a CS operator Evidently, the BdG Hamiltonian for the helical superconductor has this property. After the chiral symmetry operation,ǫ(k),∆(k) and φ in equation (2) will become −ǫ(−k), ∆(−k) and −φ. Finally, we consider the spin rotation operation. The rotation matrix in the particlehole space is diagonal, i.e., R(ξ, η) = diag(R(ξ, η),R * (ξ, η)) with the element in the spin space. The elementR represents a ξ-angle rotation about the y-axis followed by a η-angle rotation about the z-axis.
for all values of ξ and η, which is not satisfied by the Hamiltonian with d ± uv (k) = ∆ 0 (k ux ± k vŷ ). Actually, the helical superconductor breaks all the non-trivial rotation symmetries with the angle ξ = 0 and η = 0. Nevertheless, there are three special rotations, R(0, π), R(π, 0) and R(π, π), which can help bring the so-called partial symmetries.

Spin-singlet system
For SS, the gap matrix∆(k) = ∆(k)iσ 2 with the even parity∆(k) =∆(−k) in equation (2). The energy gap function ∆(k) = ∆ 0 for the s-wave pairing, ∆(k) = ∆ 0 (k 2 x − k 2 y ) for the d x 2 −y 2 -wave paring and ∆(k) = ∆ 0 2k x k y for the d xy -wave paring. The spin-singlet system not only preserves TRS, PHS and CS but also obeys SRS. In order to conveniently discuss the symmetries of transport quantities later, here we give the changes ofǫ(k) and∆(k) in equation (2) after the operations T , C and S, The changes brought by the rotation operations for SS are presented in table 1. For later use, the transformed gap matrix has been expressed with (k x , −k y ) in order to keep consistent with equations (14)- (17). When we do this, there will be a π-phase difference between the gap matrix for s(d x 2 −y 2 )-wave and that for d xy -wave. The difference originates from their opposite 'parities' under k y → −k y as shown in figure 2, although they are both spin-singlet even parity superconductors. The subtle difference between the gap functions can bring important physical results of Josephson effects which will seen in part 4. For simplicity,ǫ(k) is omitted in the table since it is invariant under the operations and it is a even function both for k x and k y .

Ferromagnetic system
The BdG-type Hamiltonian in equation (1) for F is denoted by H F (k) which can be written as where the magnetization M = Mn(θ m , φ m ) is specified by the directionn = (sin θ m cos φ m , sin θ m sin φ m , cos θ m ) with the polar angle θ m and the azimuthal angle φ m ;V 0 = V 012×2 is the non-magnetic potential. The Hamiltonian H F (k) can describe the bulk F or the ferromagnetic interface in heterostructures. We takeV 0 = 0 for the bulk and the purely ferromagnetic interface cases.
The ferromagnetic system breaks TRS, PHS and SRS, i.e., the Hamiltonian in equation (21) does not satisfy equations (8)- (10). For later use, we also present the changes of the Hamiltonian brought by the T , C and S transformations as what we have done for SS, From the above equations, it is found the direction of magnetization is 'rotated' after T and C operations; the energyǫ(k) andV 0 become negative after the S operation.
The actions of rotation operations are summarized in table 2. A pure rotation and its corresponding gauge-rotation lead to the same change of the ferromagnetic system due to the U 1 gauge symmetry obeyed by F. For simplicity, we do not showǫ(k) andV 0 in the table since they are invariant under rotation operations.  The scattering process for an electron-like quasiparticle (EQ) injected from F. The solid lines represent the normally reflected and the transmitted EQs. The dashed lines represent the Andreev reflected and the transmitted hole-like quasiparticles (HQs). The spin for the injected EQ is assumed alongn(θ m , φ m ); quasiparticles denoted by lines of the same color possess the same spin.
We consider the F|helical p-wave TS junctions as shown in figure 3(a). The nonmagnetic interface with M = 0 is located at x = 0 and along y axis. The ferromagnetic region with x < 0 is described by the Hamiltonian in equation (21) with V 0 = 0; the superconducting region with x > 0 is described by the Hamiltonian in equation (2). The wave functions of quasiparticles in F and TS can be obtained through solving the For the spin of quasiparticles in F, we chose the directionn(θ m , φ m ) as the quantization axis. When an electron with wavevector k y is injected from F, there will be four physical processes as shown in figure 3(b). The electron is normally reflected as electrons and Andreev reflected as holes; the electron transmits into TS as electron-like quasiparticles and hole-like quasiparticles. The wave function in F (TS) is the superposition of the reflected (transmitted) quasiparticles. The reflection and the transmission coefficients of the quasiparticles can be derived from boundary conditions at the interface, which are functions of M, V 0 ,n(θ m , φ m ), k y and the bias V . According to the Blonder-Tinkham-Klapwijk formalism [49], the angle-resolved conductance σ ′ of the junctions can be expressed as the composition of the coefficients, which is a function of θ m , φ m and k y when M, V 0 and V are fixed. As a result, the angle-averaged conductance σ can be given by where C 0 is a constant independent of θ m , φ m and k y . Generally speaking, the expression of the angle-averaged conductance is very complex. A more detailed description for the derivation of conductance can be found in Appendix.

F|helical TS junctions
Now, we clarify the symmetries of the conductance. We want to know which orientations of magnetization can lead to the same conductance. First, if the superconductor is invariant under some unitary operation and simultaneously the operation rotates the direction of magnetization fromn . For SS, the spin of Cooper pairs is zero ; the system obeys the full SRS. Consequently, the conductance is independent of the orientation of magnetization, i.e., For the chiral p-wave superconductor with d-vector d = ∆ 0 (k x ±ik y )ẑ which is parallel to the crystallographic c-axis, the spin of Cooper pairs lie in the ab-plane. The system obeys the symmetry of the spin rotation about the z-axis. Consequently, the conductance is invariant under the magnetization rotation about the z-axis, i.e., σ(θ m , φ m ) = σ(θ m , φ ′ m ) for all values of φ m and φ ′ m . For the helical TS, the symmetry operation is the gauge-rotation R ′ (0, π) which 'rotates' the direction of magnetization fromn(θ m , φ m ) ton(θ m , π + φ m ) as given in table 2. Thus, we obtain The second situation which can lead to the invariance of conductance is the partial symmetries of the BdG Hamiltonian. The conductance σ(θ m , φ m ) is an average value of the angle-resolved conductance on k y , which indicates that the junctions with k y and those with −k y correspond to the same conductance. For the p ± xy -wave and the p ± yxwave TSs, the partial symmetry operations are R ′ (π, 0), R(π, π) and R(π, 0), R ′ (π, π), respectively, which transform the Hamiltonian The jonit operations R ′ (0, π)R ′ (π, 0) and the pure operation R(π, π) give the same symmetry relation of conductance due to the U 1 gauge symmetry satisfied by F. The magnetization with θ m = 0 or π is a special case which means the magnetization is parallel to the spin quantization axis of the helical TS. In the coordinate of spin space in F, the TS are in purely equal-spin paring states. When an electron is injected from F, the normal reflected electron, the Andreev reflected hole and the transmitted quasiparticles all possess the same spin as that of the injected electron. As a result, the eight lines denoting the scattering processes presented in figure (3) will degenerate into four. The conductance in this situation is irrespective of the azimuthal angle φ m . The equations (27) and (28)   The obtained results in this part demonstrate that conductance, as an observable quantity, possesses a higher symmetry than the system itself. The SRS breaking of Hamiltonian favors strong anisotropy of conductance when the direction of magnetization is changed. However, the remaining partial symmetries of Hamiltonian help to keep invariance of conductance. In addition, the symmetry relations in equations (26)- (28) have been proved in Ref. [38] where the detailed derivation and numerical calculations are presented. Here, for the convenience of comparison between experiments and our results, we present in figure (4)   The theory of quasiclassical Green's function is a powerful tool to study the Josephson effect in superconducting junctions, in which the Josephson current can be expressed by the retard Green's function. One method of obtaining the retard Green's function is to solve the Eilenberger equation [50]. The Riccati parametrization of the function and the diagrammatic representation of boundary conditions have been developed for the method [51]. Another effective method of obtaining the retard Green's function is to construct it by composing wave functions derived from the BdG-equation H(k)Ψ = E(k)Ψ, which is innovated by McMillan [52] and extended to the case of anisotropic superconductor by Tanaka et al. [53] In this paper, the latter method is adopted.
We consider two-dimensional Josephson junctions with a ferromagnetic interface located x = 0 and along the y-axis. The interface is described by the Hamiltonian in equation (21). For a given k y and the fixed spin, there are four types of quasiparticle injection process: the left or right injection of the electron-like quasiparticle and the hole-like quasiparticle. For each injection, there exist four scattering processes as shown in figure 5. The wave function in superconductor is the superposition of the scattered quasiparticles. For the given k y , there are four types of conjugated injection process with wavevector −k y . The wave function in superconductor for the conjugated process can also be obtained from the superposition of the scattered quasiparticles. The retard Green's function is constructed from a linear combination of the product of wave functions and the transpose of the conjugated wave functions. The combination coefficients can be solved by using the boundary conditions at the ferromagnetic interface.
After an analytical continuation, the obtained Green's function G(x, x ′ , k y , iω n ) is not only dependent on the spatial position, wavevector k y and the Matsubara frequency ω n but also a function of M,n(θ m , φ m ),V 0 and the superconducting phase difference ∆φ. The Josephson current for fixed M and V 0 at the interface can be expressed as where C 1 and C 2 are constants which are irrelevant to the symmetries of the Josephson current. For simplicity, in the following we will assume the superconducting phase of the left superconductor is φ and that of the right is zero. Under this assumption, we have I J (θ m , φ m , ∆φ) = I J (θ m , φ m , φ). A concise discussion about how the symmetries of Hamiltonian affect the Josephson current is given in Appendix. In the following, we use equation (30) as the starting point to analyze the properties of the current.

SS|F|helical TS Josephson junctions
We consider the Josephson junctions as shown in figure 5(a). The SS with phase φ is located in the region with x < 0 and the TS with zero phase in the region with x > 0. They are described by the Hamiltonian in equations (1) and (2). The ferromagnetic interface is described by the Hamiltonian in equations (1) and (21). Firstly, we discuss the general results for the Josephson current of SS|F|helical TS junctions. For one thing, the time-reversal operation will change the superconducting phase of SS from φ to −φ according to equation (18). The direction of the Josephson current is reversed accordingly, i.e., I J → −I J . For another thing, the time-reversal 'rotates' the magnetization in F from the directionn(θ m , φ m ) ton(π − θ m , π + φ m ). As a result, one obtains The C operator is another symmetry operation which can lead to a general result. The energy of quasiparticles becomes negative while the phase φ remains unchanged under the operation. Accordingly, the quantity F (iω n , k y ) will become F (−iω n , k y ) in equation (30). However, the Josephson current I J , as the sum of F (iω n , k y ), will keep its value. Actually, after the particle-hole transformation, the helical superconductor can be viewed as a new electron system with the normal dispersion −ǫ(−k) and the energy gap matrix∆(−k). The SS can be viewed as a new electron system with the normal dispersion −ǫ(−k) and the energy gap matrix∆(−k)e iπ according to equation (19). The ferromagnetic interface becomes a new electron system with the normal dispersion −ǫ(−k), the non-magnetic potential −V 0 and the magnetization alongn(π −θ m , π +φ m ) according to equation (23). As a result, we have Since an actual interface possesses non-negative potential, i.e., V 0 ≥ 0, only the case of V 0 = 0 for equation (32) makes sense. Then, we have with V 0 = 0. Combing equations (31) and (33), we obtain which is just the result brought by the chiral symmetry operation. The operator R ′ (0, π) can also bring a general result. Under the operation, the helical TS is invariant as given in equation (13) while the SS will acquire a π phase as shown in table 1. Meanwhile, the operation changes the direction of magnetization from n(θ m , φ m ) ton(θ m , π + φ m ). As a unitary transformation, R ′ (0, π) will not change the Josephson current. As a result, one obtains Secondly, we discuss other symmetry operations which will give different results for junctions with different pairing symmetries of the superconducting wave functions. Now, we take R ′ (π, 0) as an example to analyze the symmetry of the Josephson current. Under the operation, the junctions with k y become those with −k y due to the partial symmetry of the p ± xy -wave superconductor. Simultaneously, the s(d x 2 −y 2 )-wave SS acquires a phase of π according to table 1; the direction of magnetization turns inton(π − θ m , π − φ m ) according to table 2. However, as an average quantity on k y , the Josephson current in equation (30) will not change its value. Thus, we have for the s(d x 2 −y 2 )-wave SS|F|p ± xy -wave TS junctions. Similarly, from the operation R(π, π) we can derive for the s(d x 2 −y 2 )-wave SS|F|p ± xy -wave TS junctions. From the operations R ′ (π, 0) and R(π, π), we can also derive and respectively, for the d xy -wave SS|F|p ± xy -wave TS junctions. For the SS|F|p ± yx -wave TS junctions, there are similar symmetries of the Josephson current. The s(d x 2 −y 2 )-wave SS|F|p ± yx -wave TS junctions satisfy equations (38) and (39) which originate from the operations R(π, 0) and R ′ (π, π), respectively. The d xy -wave SS|F|p ± yx -wave TS junctions satisfy equations (36) and (37) which also originate from the operations R(π, 0) and R ′ (π, π), respectively. From the results, we come to the following conclusions.
(a) The same spin-rotation or gauge-rotation operation can bring different symmetries of the Josephson current for the s(d x 2 −y 2 )-wave SS and the d xy -wave SS when the gap function in the helical superconductor is fixed. The difference results from the partial symmetry obeyed by the helical superconductor and the different k yparities of gap functions in SS as discussed in part 2.1.2.
(b) The SS|p ± uv -wave TS junctions can be classified into two types depending on the symmetries satisfied by their current. The s(d x 2 −y 2 )-wave SS|F|p ± xy -wave TS junctions and the d xy -wave SS|F|p ± yx -wave TS junctions satisfy equations (36) and (37); the s(d x 2 −y 2 )-wave SS|F|p ± yx -wave TS junctions and the d xy -wave SS|F|p ± xy -wave TS junctions satisfy equations (38) and (39). This classification is consistent with the numerical results in Ref. [39]. Especially, the junctions of the same type possess the same selection rules of the lowest order Josephson current. Finally, we show a brief explanation for the relation of the results in this paper to those in Ref. [39]. The combination of equations here can bring equations in Ref. [39]. For example, equations (33) and (35) will give I J (θ m , φ m , φ) = I J (π − θ m , φ m , φ); equations (33) and (36) will give I J (θ m , φ m , φ) = I J (θ m , 2π − φ m , φ); equations (33) and (39) will give I J (θ m , φ m , φ) = I J (θ m , π − φ m , φ). The three equalities are just the equations (8), (12) and (14) in Ref. [39], respectively. As pointed out in the reference, the three equations will not hold when the non-magnetic potential V 0 is non-zero (i.e., Z = 0 there). This is because the equations are all based on equation (33) which holds only when the non-magnetic potential V 0 = 0. In order to compare our results with experiments easily, we present in figure (6) the Josephson current for the SS|F|p + xy -wave TS junction in the orientation space of magnetization and in the space spanned by the azimuthal angle φ m and the phase φ, which can explicitly express the symmetry relations.

Helical TS|F|Helical TS Josephson junctions
We consider the purely helical Josephson junctions as shown in figure 5(a). The p + xywave TS in the left side has the superconducting phase φ. The right side can be the p − xy -wave or the p ± yx -wave TS with the zero phase. The magnetization in the interface is assumed alongn(θ m , φ m ). Now ,we discuss the symmetries of the Josephson current. The time-reversal and the charge conjugation lead to respectively, which are the common results for the purely helical junctions. The latter equality is different from equation (32). The helical superconductor will not acquire a phase π under the charge conjugation.  (42) and (43). The connection between the two figures reflects the symmetry relations in equation (40). The other parameters for the figures are taken as: T = 0.3T C , M EF = 0.9 and 2mV0 2 kF = 1. The current has been expressed by the resistance in the normal state [40].
The rotation operations can bring different results for different junctions. For the p + xy -wave TS|F|p − xy -wave TS junction, R(0, π) and R(π, 0) give and respectively, which are related to the transformations in equations (13) and (14). The result for R(π, π), as the joint operation of R(0, π) and R(π, 0), is the combination of equations (42) and (43), i.e., I J (θ m , φ m , φ) = I J (π − θ m , 2π − φ m , φ). In figure (7), we show the Josephson current for the p + xy -wave TS|F|p − xy -wave TS junction in the orientation space of magnetization. It may be useful in comparing our results with experiments. For simplicity, we will not give numerical results for other junctions.
There is a π-phase difference between equation (43) and (45) since the p + xy -wave TS is partially symmetric under R ′ (π, 0) while R(π, 0) will bring an extra phase π to the superconductor.
The helical Josephson junctions can host rich ground states which can be classified into 0 phase, π phase, 0 + π phase, ϕ 0 phase and ϕ phase [40]. The phase transition can be tuned by rotating the magnetization in the ferromagnetic interface. The phase diagram formed in the orientation space of magnetization is an important aspect of researches on Josephson junctions which also possesses some symmetries as results of symmetries of the current.
For the p + xy -wave TS|F|p − xy -wave TS junction, the diagram in the orientation space 0 ≤ θ m < π and 0 ≤ φ < 2π is symmetric about both θ m = π/2 and φ m = π. It is the result of equations (42), (43) and Equation (46) originates from the close relation between p ± xy -wave TSs. When one changes the y-axis into −y in the three dimensional coordinate space, the TSs exchange each other; the direction of magnetization in the new coordinate system becomeŝ n(θ m , 2π −φ m , φ) accordingly. The relation can be reflected directly from the interaction of helical superconductivity and ferromagnetism which can be seen in part 5.
For the p + xy -wave TS|F|p ± yx -wave TS junctions, the phase diagram for V 0 = 0 is symmetric about θ m = π/2 and invariant under the translation φ m → φ m + π. It is the result of equations (44) and (41) with V 0 = 0. In addition, we have another equality, which originates from the connection between the TSs. For example, for the p + xy -wave TS|F|p + yx -wave TS junction, when one changes the x-axis and the y-axis into y and x respectively, the TSs exchange each other; the direction of the magnetization becomeŝ n(θ, π/2 − φ m , φ) in the new coordinate space. The connection can also be reflected directly from the interaction of helical superconductivity and ferromagnetism which can be seen in part 5.
As discussed in part 3.2 for the conductance, the magnetization with θ m = 0 or π is also a special case for the Josephson current. The purely equal-spin paring will simplify the scattering process in figure (5). The current will become irrespective of the azimuthal angle φ m . For the SS|F|TS and TS|F|TS junctions, we will have (31) and (40). This relation can also be found in the Josephson junction formed along the edges of quantum spin-Hall insulators with a Zeeman field along the spin quantization [46].
Another special case is the perpendicular magnetization with θ m = π/2 which is important for helical superconductor realized in a topological insulator with the s-wave superconductor [54,55]. The symmetries of the transport quantities in these references are different from those in this paper. For example, in reference [55], the transmission probability of the normal-metal-superconductor junction with a single ferromagnetic barrier is independent of the direction of the in-plane magnetization. For the normalmetal-superconductor junction with double barriers, the probability is dependent on the relative angle of the in-plane magnetization in the barriers. However, in our junctions with θ m = π/2, the conductance and the current is strongly dependent on the azimuthal angle of magnetization due to the symmetry breaking of the spin-rotation.
Finally, we discuss the Josephson current in the non-magnetic case, i.e., M = 0. In this case, we have I J (φ) = I J (π + φ) for the SS|TS junctions which means the Josephson current is π-periodic. This result is consistent with that for the junction between singlet and triplet superconductors in the static thermodynamic limit [56]. For the p + xy -wave TS|p ± yx -wave TS junctions, we also have I J (φ) = I J (π + φ) with π periodicity which also indicates the absence of the sin φ-type current. This result is also consistent with that for the triplet junction in which the two d-vectors are perpendicular [56]. For the helical states with k-dependent d-vectors in this paper, the perpendicular vectors possess the relation d + xy · d ± yx ky = 0 in which · · · ky denotes the average over the momentum parallel to the interface [40].

Symmetries of free energy
In the above sections we have analyzed the symmetries of the current in Josephson junctions from the viewpoint of symmetries of Hamiltonian. Actually, the current-phase relation has a more direct relationship with free energy of junctions, which is derivative of the free energy with respect to φ. The free energy of the Josephson junctions have been constructed in Ref. [39] and [40] on the base of numerical results, which directly reflects the interaction of helical superconductivity and ferromagnetism. In the following sections we will show that symmetries of the current derived from Hamiltonian are consistent with those derived from the free energy.

SS|F|helical TS Josephson junctions
The free energy for the s(d x 2 −y 2 )-wave SS|F|p ± xy -wave TS and the d xy -wave SS|F|p ± yx -wave TS junctions, denoted by F ST 1 , is given by for V 0 = 0 and for V 0 = 0. The proportional coefficient before each term in the expressions is generally a complex function of parameters such as M, V 0 and ∆ 0 , which is irrespective of the symmetries of the free energy. We have omitted them for simplicity. The symbol & suggests that there are two terms contributed to the current when V 0 = 0. The Josephson current is given accordingly by for V 0 = 0 and for V 0 = 0. It is easy to verify that the current for V 0 = 0 satisfies equations (31) and (33)- (37). However, equations (33) and (34) do not hold for the current with V 0 = 0 due to the presence of the second term in equation (51), which is consistent with the true condition of equations (33) and (34). The free energy for the s(d x 2 −y 2 )-wave SS|F|p ± yx -wave TS and the d xy -wave SS|F|p ± xywave TS junctions, denoted by F ST 2 , is given by for V 0 = 0 and for V 0 = 0. The Josephson current is accordingly given by for V 0 = 0 and for V 0 = 0. It is easy to verify that the current for V 0 = 0 satisfies equations (31), (33)- (35), (38) and (39). However, equations (33) and (34) do not hold for the current with V 0 = 0 due to the presence of the second term in equation (55), which is consistent with the true condition of equations (33) and (34).

Helical TS|F|helical TS Josephson junctions
The free energy for the p + xy -wave TS|F|p − xy -wave TS junction, denoted by F T T 1 , is given by The Josephson current is accordingly given as It can be easily demonstrated that the current satisfies equations (40)- (43). In addition, equation (46) can be derived directly from equation (57) which is not obeyed by the SS|F|helical p-wave TS junctions due to the absence of the close relation discussed in part 4.3. The free energy of the p + xy -wave TS|F|p ± yx -wave TS junctions, denoted by F T T 2 , is given by for V 0 = 0 and for V 0 = 0. The Josephson current is accordingly given as for V 0 = 0 and for V 0 = 0. It can be demonstrated that the current satisfies equations (40), (41), (44) and (45). In addition, equations (46) and (47) can be directly derived from equations (60) and (61), which are not satisfied by the SS|F|helical p-wave TS junctions due to the absence of the connection discussed in part 4.3.

Conclusions
We establish linkages between symmetries of the Hamiltonian for topological superconducting systems and the invariance of transport quantities for the corresponding junctions. As observables, transport quantities exhibit higher symmetries than the systems themselves. We reveal the important role that partial symmetries play in the invariance of conductance and Josephson current. Our analysis explain the numerical results reported recently which include the selection rules of the lowest order current and the rich phase diagrams in Josephson junctions. The analysis method not only provides a profound understanding of topological junctions but also helps give useful information about transport quantities before carrying out complex numerical calculations. In addition, the symmetry analysis of transport properties in this paper is general and applicable for other topological junctions such as the chiral ones. However, the present studies do not include the 4π-periodic Josephson effect caused by transmission of unpaired electrons [57]. The novel effect is related to the fermion-parity anomaly in the superconducting ground state.
The wave functions for the transformed Hamiltonian in the ferromagnetic region can be obtained by the similar substitutions for θ m and φ m in Ψ F as discussed; the wave function in the superconducting region can be obtained by replacing k y with −k y in Ψ S . Under the boundary conditions, the coefficients and thus the conductance for the transformed system can be determined, i.e., σ(π − θ m , π − φ m ) = C 0 ky σ ′ (π − θ m , π − φ m , −k y ) = C 0 ky σ ′ (π−θ m , π−φ m , k y ) and σ(π−θ m , 2π−φ m ) = C 0 ky σ ′ (π−θ m , 2π−φ m , −k y ) = C 0 ky σ ′ (π − θ m , 2π − φ m , k y ). Then we have the relations in equations (27) and (28) due to the unitarity of R ′ (π, 0) and R(π, π). For other F| helical TS junctions, one can derive the symmetry relations in a similar way.
Finally, we discuss how to derive symmetries of the Josephson current from the symmetries satisfied by Hamiltonian. In order to express the current, we first need to obtain the retard Green's function which can be constructed with the scattering wave function in superconductors [52,53]. The wave function and the scattering coefficients can be obtained by solving the BdG equation under boundary conditions. The solving process is the same as that for the conductance case. Then, the changes of Hamiltonian under the symmetry transformations will enter into the wave function and the coefficients. Actually, the finial expression of the current can be written as the combination of the Andreev reflection coefficients [58], i.e., F (iω n , k y ) in equation (30) is proportional to a(θ m , φ m , iω n , k y , φ). The effects of symmetry transformations on the Andreev coefficient have been discussed carefully. The symmetries of the current can be derived in a similar way.