searching for new materials with nontrivial topological properties is a forefront research field of condensed matter physics . following the theoretical prediction and experimental discovery of two - dimensional and three - dimensional topological insulators @xcite ,
the fully - gapped topological insulators and topological superconductors have been classified according to their symmetries @xcite .
one focus of recent research is on a class of materials with discrete degeneracy points in the band structure @xcite .
these dirac and weyl semimetals are analogies or generalizations of graphene , which is a prototypical two - dimensional ( 2d ) dirac semimetal @xcite .
the three - dimensional ( 3d ) weyl semimetals realize the chiral anomaly which results in several novel transport properties , including the anomalous hall effect and the chiral magnetic effect @xcite .
the interest in 3d weyl semimetal is soon extended to 2d dirac and weyl semimetals @xcite and 3d dirac and weyl superconductors @xcite .
the 2d weyl superconductor , which is also a member of this new family of materials characterized by the presence of weyl fermions , has however till now received very little attention .
this is partly because of the difficulty in realizing such a phase . on one hand ,
accidental band degeneracy is known to be very unlikely in 2d @xcite . in a superconductor , at least one particle band and one hole band ( which are related to each other by charge conjugation )
cross each other and thus invert in a certain region of the wave vector space , which makes possible the occurrence of nodal points as the parts of the crossing curve between the particle and hole bands that are not gapped out .
however , by opening a full gap the superconductor can usually gain more condensation energy and so naturally occurring nodal superconductors are rare . on the other hand , the realization of weyl nodes requires non - degenerate fermi surface .
this requires the breaking of at least one of the time - reversal symmetry and the inversion symmetry @xcite . while the former requires the persistence of magnetic order in the superconducting phase , the latter point to noncentrosymmetric materials . in bulk materials ,
the competition with other symmetry breaking tendencies makes the superconducting phase less favored as compared to the simple time - reversal invariant and inversion symmetric materials .
therefore , the 2d weyl superconductor appears to be the most elusive member of the weyl materials .    in this work
, we propose that a 2d weyl superconductor can be created in a very inclusive system , a heterostructure consisting of a half - metal ( hm ) with rashba spin - orbital coupling ( rsoc ) and a spin - singlet superconductor .
the combination of a single nondegenerate fermi surface @xcite , the transmutation between spin - singlet and spin - triplet pairing by the rsoc @xcite , and the possible existence of nodes in the order parameter of the spin - singlet pairing , make the hm layer of the proposed heterostructure a unique platform to realize the 2d weyl superconductor .
an important tunability of the system , sketched in fig.1 , is the orientation of magnetization of the hm .
earlier , the heterostructure with an @xmath1-wave superconductor was found to give a 2d weyl superconductivity phase ( swsc ) when the spontaneous magnetization of the hm lies along a mirror symmetric direction within the film plane @xcite . in the heterostructure with a @xmath0-wave superconductor ( such as the cuprates ) , on the other hand , a 2d weyl superconductor ( dwsc ) is realized for an arbitrary magnetization orientation ( the case with purely out - of - plane magnetization was also studied by sato and fujimoto in a slightly different setting @xcite ) .
the bulk and edge state spectra and their evolution with the magnetization direction show distinct differences from those for the swsc .
for example , besides the type - i weyl nodes which are realized in swsc , we also find type - ii weyl nodes in a dwsc . in weyl
semimetal , type - i and type - ii weyl nodes are defined separately as weyl point with a point - like fermi surface and weyl point appearing at the contact of electron and hole pockets @xcite . here , we define type - i and type - ii weyl nodes in a superconductor analogously .
namely , for type - i weyl node the constant energy surface ( not necessarily the fermi surface ) crossing the weyl point spans a single point , whereas for type - ii weyl node the constant energy surface crossing the weyl point forms an electron - like and a hole - like quasiparticle ` fermi surface ' and connect to each other at the weyl point .
the transition between type - i and type - ii weyl nodes is achieved in a dwsc by varying the direction of the spontaneous magnetization of the hm .
in addition , for a general magnetization orientation , the state studied consists simultaneously of normal and superfluid components . in terms of the zero - energy quasiparticle spectrum
, the varying magnetization orientation drives the system to go through a series of lifshitz transitions , which is usually encountered in a metal in the normal phase @xcite .
it is also shown to be qualitatively different from all known types of pairings , including the bcs pairing @xcite , the fflo pairing @xcite , and the breached pairing state ( sarma phase ) @xcite .
the organization of this paper is as follows . in sec .
ii , we define the model and parameters for the proposed heterostructure . in sec .
iii , we report the energy spectra of the system .
these include the energy spectra for bulk and a strip of the system , the spectral functions and density of states for the two edges of a strip , the clarification of the nature of proximity - induced superconducting state , and the lifshitz transitions in the quasiparticle spectrum . the transition between type - i and type - ii weyl nodes is also studied in this section .
( the analytical analyses are put to appendix a. ) in sec .
iv , we study the symmetry and topological properties of the system . in sec .
v , we show the energy spectra for the heterostructure with the hm replaced by an ordinary ferromagnetic material having a partially polarized fermi surface .
finally , we summarize our results in sec . vi .
[ fig1 ] -wave or @xmath0-wave spin - singlet superconductor and an insulating substrate , viewed laterally .
@xmath2 and @xmath3 represent the magnitude and direction of the exchange field ( magnetization ) in the hm , which is assumed to lie on the @xmath4 plane.,title="fig:",width=283,height=113 ] +    figure 1 is a schematic plot of the system that we will study .
it consists of a hm ( or , highly polarized ferromagnetic metal ) thin film grown on an insulating substrate and capped with a spin - singlet superconductor with @xmath1-wave or @xmath0-wave pairing symmetry . to elucidate the principle more clearly
, we consider the simplest model for the system .
we describe the hm thin film by a one - orbital model defined on a square lattice .
perfect interfaces between the hm and the substrate and the spin - singlet superconductor are assumed , so that we will ignore the effect of impurity and disorder . denoting the basis vector as @xmath5 $ ] , the model hamiltonian for the hm thin film , with a rsoc term induced by the formation of the heterostructure , is @xmath6 , where @xcite @xmath7 @xmath8 is the two - dimensional unit matrix , @xmath9 ( @xmath10 ) are pauli matrices in the spin subspace .
@xmath12 is the hopping amplitude and we assume @xmath13 , @xmath14 is the chemical potential , @xmath15 is the strength of the rsoc . in writing down eq.(1 ) , we have assumed the magnetization to be fixed on the @xmath4 plane . in terms of @xmath2 ( @xmath16 ) and @xmath3 ( @xmath17 , with respect to the positive direction of the @xmath18 axis ) , which are separately the magnitude and polar angle of the magnetization ( see fig.(1 ) )
, we have @xmath19 , @xmath20 . to describe a hm with @xmath21
, we set @xmath2 to be of the same order of magnitude as @xmath12 , and is much larger than the externally induced @xmath15 . the chemical potential is tuned to make sure that it crosses only the lower spin - split band of @xmath22 , which amounts to @xmath23 .
the proximity - induced superconductivity in the hm arising from coupling with an @xmath1-wave or @xmath0-wave spin - singlet superconductor is described by @xmath24 where @xmath25 @xcite .
@xmath26 is a multiplication of the pairing amplitude @xmath27 and a symmetry factor @xmath28 . for the sake of simplicity , we take @xmath29 for @xmath1-wave superconductor and @xmath30 or @xmath31 for @xmath0-wave superconductor @xcite . the two latter cases can be realized by growing the same @xmath0-wave superconductor ( e.g. , cuprates ) on the top of the hm film along two special directions , respectively .
the @xmath32-axis , which is the projection of the magnetization plane ( @xmath4 ) onto the plane of the films ( @xmath33 ) in the heterostructure , is along the antinodal direction of the @xmath0-wave pairing gap for @xmath30 and along the nodal direction of the @xmath0-wave pairing gap for @xmath31 .
the heterostructure formed by growing the @xmath0-wave superconductor along more general ( less symmetric ) directions , while possible in principle , can be understood in terms of a combination of the results for these two special cases . in the nambu basis ,
@xmath34 $ ] , the full model is written as @xmath35 , where the bogoliubov - de gennes ( bdg ) hamiltonian with proximity - induced pairing term and rsoc is @xmath36 @xmath37 ( @xmath10 ) are pauli matrices in the nambu space . diagonalizing eq.(3 ) gives the four quasiparticle bands @xmath38 , with @xmath39=@xmath40 , @xmath41 , @xmath42 , @xmath43 and larger @xmath39 means higher energy .    in the hm thin film ,
only the spin - polarized band crossing the chemical potential is important to the low - energy properties of the system . to facilitate the analysis
, we construct the low - energy effective model by projecting to the subspace containing only information of the band contributing to the fermi surface .
first , we diagonalize the model defined in eq.(1 ) by a unitary transformation @xmath44 .
@xmath45 is a diagonal matrix with the two diagonal elements storing the eigenvalues of @xmath22 , @xmath46 where @xmath47 .
notice that , @xmath48 ( @xmath47 ) are the electronic energy bands in the normal state ( i.e. , without the proximity - induced pairing term ) and should be distinguished from the quasiparticle bands @xmath38 ( @xmath49 - 2 , -1 , 1 , 2 ) defined as eigenvalues of eq.(3 ) .
the unitary matrix @xmath50 containing the eigenvectors of @xmath22 is written as @xmath51 the normalized eigenvectors are chosen as @xmath52 where @xmath47 and @xmath53 .
@xmath54 is defined in eq.(4 ) . in the eigenbasis ,
the pairing term is rewritten as @xmath55 for @xmath56 that we will focus on , only the @xmath57 band contributes to the low - energy properties .
we thus retain only this band and the pairing term within this band to construct the low - energy effective model . in the reduced nambu space defined by the basis
@xmath58 $ ] , the low - energy effective model is written as @xmath59 the intra - band pairing amplitude is @xmath60   \notag \\ & = & \frac{2\delta_{0}\eta(\mathbf{k})}{d_{-}(\mathbf{k})d_{-}(-\mathbf{k})}[\lambda^{\ast}_{\mathbf{k}}\tilde{m}_{z}(\mathbf{k})+\tilde{m}_{x}(\mathbf{k } ) ]   \notag \\ & \equiv&\delta_{0}\eta(\mathbf{k})f(\mathbf{k}),\end{aligned}\ ] ] where we have introduced the abbreviations @xmath61 , @xmath62 , and @xmath63 .
the effective paring is therefore determined by two symmetry factors , @xmath28 and @xmath64 .
the eigenstates of @xmath65 will be used as approximations for the two low - energy quasiparticle bands , @xmath66 and @xmath67 , in the following analytical analysis .
for the heterostructure with @xmath1-wave superconductor and in - plane exchange field ( @xmath68 or @xmath69 ) , the proximity - induced superconducting phase in the hm thin film is an swsc with a single pair of weyl nodes located on the mirror invariant line of the 2d brillouin zone ( bz ) @xcite .
it has also been verified that , the swsc is a direct result of an emergent mirror reflection symmetry that exists only for an in - plane exchange field @xcite .
the heterostructure with a @xmath0-wave superconductor , however , supports a dwsc phase independent of the @xmath3 angle of the exchange field .
this is understandable from the intraband pairing amplitude in eq.(9 ) .
namely , because of the presence of the @xmath70 factor , there are always point nodes inherited from those of the @xmath0-wave superconductor .
the remaining wave - vector dependence of the effective pairing amplitude , contained in @xmath64 , is the same as that in the heterostructure with @xmath1-wave superconductor . in particular , @xmath71 @xcite . from the positions of the nodes resulting from the @xmath28 factor and from the @xmath64 factor
, we know that as @xmath3 varies the number of pairs of weyl nodes vary between two and three for @xmath30 and remains to be two in the whole range of @xmath3 for @xmath31 .
we study the energy spectra of the two cases separately in what follows .
two kinds of topological phase transitions were identified in the heterostructure with @xmath1-wave superconductor , as we vary @xmath3 @xcite .
one is between a fully gapped phase and a bulk gapless phase with ubiquitous local gap .
another is between the bulk gapless phase with ubiquitous local gap and the swsc phase with a pair of weyl nodes connecting the two low - energy quasiparticle bands @xmath66 and @xmath67 . for the heterostructure with @xmath0-wave superconductor ,
the persistence of the weyl nodes for all @xmath3 angles makes exact analogies of the corresponding phase transitions nonexistent .
however , we can still define several critical angles around which to expect some qualitative changes in the quasiparticle energy spectrum for the system in bulk and strip configurations . in this respect , for @xmath30
, we can define three critical angles in the range @xmath72 $ ] as follows .
the critical angles are defined in terms of the quasiparticle energies of the four weyl nodes and those for two special fermi points . for the model parameters of interest to us ( @xmath73 ) ,
the positions of the reference wave vectors keep almost unchanged as @xmath3 varies .
the four weyl nodes are at @xmath74 with @xmath75 determined by @xmath76.\ ] ] the two special fermi points are chosen as @xmath77 , with @xmath78 determined by @xmath79.\ ] ] the critical angle @xmath80 is determined by the condition @xmath81 , where @xmath82 is the sign function . the critical angle @xmath83 is determined by the condition @xmath84 . making use of the low - energy effective model defined by eqs.(8 ) and ( 9 ) , @xmath80 and @xmath83 are determined approximately by @xmath85 and @xmath86 the third critical angle is @xmath87 , at which the number of weyl nodes changes from two pairs to three pairs . for small @xmath88
, @xmath77 gives the approximate positions of the two additional weyl nodes .
we proceed to show the numerical results for the energy spectra of the dwsc in both bulk and strip configurations , for several typical values of @xmath72 $ ] .
the bulk energy spectra are shown along several lines with fixed values of @xmath89 , since the bulk quasiparticle spectrum is symmetric with respect to @xmath90 .
we focus on a typical set of parameters as @xmath91 , @xmath92 , @xmath93 , @xmath94 , and @xmath95 .
a relatively large value of @xmath27 is adopted to see the qualitative features of the spectra more clearly .
[ fig2 ] , for a typical set of parameters , @xmath96 , @xmath94 , @xmath95 , @xmath93 . for the bulk system ,
the dispersions of the @xmath66 and @xmath67 bands are shown along four lines in the 2d bz with fixed @xmath90 values at 0 ( black lines ) , @xmath97 ( blue lines ) , @xmath98 ( red lines ) , and @xmath99 ( green lines ) .
@xmath75 is defined by eq.(10 ) .
the edge state spectra correspond to a strip with 1500 unit cells along the @xmath32 direction .
@xmath100 for ( a ) and ( f ) , @xmath101 for ( b ) and ( g ) , @xmath102 for ( c ) and ( h ) , @xmath103 for ( d ) and ( i ) , @xmath104 for ( e ) and ( j).,title="fig:",width=321,height=566 ]    figure 2 shows the quasiparticle energy spectra for @xmath30 in both the bulk [ 2(a ) to 2(e ) ] and strip [ 2(f ) to 2(j ) ] configurations .
five values of @xmath3 , including @xmath105 [ 2(a ) and 2(f ) ] , @xmath80 [ 2(b ) and 2(g ) ] , @xmath83 [ 2(c ) and 2(h ) ] , @xmath106 [ 2(d ) and 2(i ) ] , and @xmath107 [ 2(e ) and 2(j ) ] are considered . for the bulk spectra , we plot for each @xmath3 the quasiparticle bands @xmath66 and @xmath67 along four lines in the 2d bz with correspondingly fixed values of @xmath90 , which are @xmath105 , @xmath97 , @xmath98 , and @xmath108 . from fig.2(b ) and fig.2(c ) , the critical angles @xmath80 and @xmath83 determined by eq.(12 ) and eq.(13 ) are very close to the exact values according to their definitions .
as @xmath3 varies from @xmath105 to @xmath107 , the positions of the weyl nodes determined by the symmetry factor @xmath30 keep almost unchanged . at @xmath68 , it is clear from fig.2(e ) that , besides the four weyl nodes determined by @xmath109 that are always present , two additional nodes emerges along @xmath110 as a result of @xmath71 @xcite .    in the strip configuration ( with two edges parallel to the @xmath111-axis ) , a pair of weyl nodes with the same @xmath112 ( equal to @xmath98 or @xmath113 ) and energy project to the same point of the edge bz .
the breaking of translational symmetry along the @xmath32 direction gives a finite coupling between the two weyl nodes .
it is thus natural to expect that the weyl nodes will be gapped out for a strip of finite width . besides
, it is known that there are edge modes for @xmath29 .
it is interesting to see whether the edge states persists and how they are modified by the ubiquitously existing weyl nodes .
these questions are answered partly by the results in figs .
2(f ) to 2(j ) for strips with 1500 unit cells along the @xmath32 direction .
clearly , a gap opens at the weyl nodes .
the edge modes persist but the dispersion is changed from the crossing - like structure to a shape similar to an @xmath114 symbol but then connected again to the edge modes beyond the projection of the weyl nodes @xcite .
interestingly , the two edge modes are nondegenerate .
this means that there are no additional edge modes associated with the weyl nodes arising from the symmetry factor @xmath28 . for @xmath115 , a third pair of weyl nodes along
@xmath110 emerges .
same as the case for the swsc @xcite , the edge modes become majorana fermi lines connecting the two new weyl nodes . in this case , since the edge modes pass through the projection of the two pairs of bulk weyl nodes resulting from @xmath109 , the bulk nodes extending throughout the bulk of the system act as channels coupling the edge modes .
therefore , close to the projection of the bulk weyl nodes at the edge bz , the edge modes are re - hybridized .      for @xmath31 , there are always two pairs of weyl nodes for an arbitrary @xmath116 .
one pair is along @xmath110 and the other pair is along @xmath117 .
since weyl nodes exist along both of the two special high symmetry lines , there are no simple analogies to @xmath80 and @xmath83 defined for @xmath30 .
nevertheless , we can still define two critical angles in a manner similar to the above discussions . for @xmath118 that we focus on in this work , the three special reference wave vectors for defining the two critical angles
are chosen as @xmath119 , @xmath120 , and @xmath121 .
@xmath75 and @xmath78 are determined by eq.(10 ) and eq.(11 ) , respectively . at the first critical angle @xmath122 , we have @xmath123 . at the second critical angle @xmath124 , we have @xmath125 .
there are no simple analytical expressions for @xmath126 and @xmath127 .
we obtain them numerically by the definition .
the energy spectra of the system in the configurations of a bulk material and a strip are shown in figure 3 .
the parameters as chosen as @xmath91 , @xmath92 , @xmath93 , @xmath94 , and @xmath95 .
we show the results for five values of @xmath3 , including @xmath105 [ 3(a ) and 3(f ) ] , @xmath126 [ 3(b ) and 3(g ) ] , @xmath127 [ 3(c ) and 3(h ) ] , @xmath106 [ 3(d ) and 3(i ) ] , and @xmath107 [ 3(e ) and 3(j ) ] . for the bulk spectra , we plot for each @xmath3 the quasiparticle bands @xmath66 and @xmath67 along three lines in the 2d bz with correspondingly fixed values of @xmath90 , which are @xmath105 , @xmath98 , and @xmath128 .
@xmath75 is still determined by eq.(10 ) .
@xmath129 is determined by @xmath130,\ ] ] where @xmath131 .
the persistence of two pairs of bulk weyl nodes and their locations , are confirmed by the results presented in figs.3(a ) to 3(e ) . a qualitative difference between the present case and those for the swsc and the dwsc with @xmath30 is , however , clear from the dispersions along lines with @xmath132 on figs.3(d ) and 3(e ) .
namely , the weyl nodes located at @xmath133 change from the first type ( type - i ) to the second type ( type - ii ) as @xmath3 increases and becomes larger than a critical value .
this change in the character of the weyl nodes has not been noticed in either the swsc or the dwsc with @xmath30 .
the pair of nodes located at @xmath134 , on the other hand , keeps to be type - i as @xmath3 changes . from appendix
a with more detailed analyses , the weyl nodes at @xmath134 are the only weyl nodes in the present heterostructure ( fig.1 ) that can become type - ii .
this is the first case where type - ii weyl nodes are predicted in the presence of superconducting correlation .
it is also the first example of type - ii weyl nodes in two - dimensional systems , as far as we know .
[ fig3 ] , for a typical set of parameters , @xmath96 , @xmath94 , @xmath95 , @xmath93 . for the bulk system ,
the dispersions of the @xmath66 and @xmath67 bands are shown along three lines in the 2d bz with fixed @xmath90 values at 0 ( black lines ) , @xmath98 ( blue lines ) , and @xmath128 ( red lines ) .
the edge state spectra correspond to a strip with 1500 unit cells along the @xmath32 direction .
@xmath100 for ( a ) and ( f ) , @xmath122 for ( b ) and ( g ) , @xmath124 for ( c ) and ( h ) , @xmath103 for ( d ) and ( i ) , @xmath104 for ( e ) and ( j).,title="fig:",width=321,height=566 ]    the energy spectra for a strip are shown in figs.3(f ) to 3(j ) , for five typical angles defined above .
while the two weyl nodes at @xmath133 still project to the same point ( i.e. , @xmath117 ) of the edge bz , the two weyl nodes at @xmath134 remain well separated from each other .
the broken translational invariance along the @xmath32 direction couples the two weyl nodes at @xmath133 and opens a gap in the nodes , similar to the situation for @xmath30 .
the edge modes also hybridize and open a gap at @xmath117 which is where the two weyl nodes at @xmath133 projects .
five values of @xmath3 are considered , including @xmath100 for ( a ) and ( b ) , @xmath101 for ( c ) and ( d ) , @xmath102 for ( e ) and ( f ) , @xmath103 for ( g ) and ( h ) , @xmath104 for ( i ) and ( j ) .
the parameters are taken as @xmath96 , @xmath94 , @xmath95 , @xmath93 .
the darker the color , the larger the spectral weight.,title="fig:",width=321,height=566 ]    [ fig5 ] .
five values of @xmath3 are considered , including @xmath100 for ( a ) and ( b ) , @xmath122 for ( c ) and ( d ) , @xmath124 for ( e ) and ( f ) , @xmath103 for ( g ) and ( h ) , @xmath104 for ( i ) and ( j ) .
the parameters are taken as @xmath96 , @xmath94 , @xmath95 , @xmath93 .
the darker the color , the larger the spectral weight.,title="fig:",width=321,height=566 ]    the energy spectra of the edge states can be probed by angle - resolved photoemission spectra ( arpes ) and scanning tunneling spectroscopy ( sts ) . instead of the quasiparticle energy spectrum ,
the experimentally relevant quantity is the spectral function along the edge of a strip , which is defined as the imaginary part of the retarded green s function for the edge layer of the strip @xcite .
we have calculated the retarded green s functions @xmath135 for the left ( @xmath136 ) and right ( @xmath137 ) edges of a strip with @xmath138 layers , in terms of standard iterative green s function method.@xcite as shown in figs.4 and 5 are results separately for @xmath30 and @xmath31 , corresponding to five @xmath3 values taken separately in figs.2 and 3 .
while the correspondences between the spectral functions in figs.4 and 5 and the quasiparticle energy spectra in figs.2 and 3 are very clear , the spectral functions contain important new information on how the edge states are distributed on the two edges .
in addition , the green s functions obtained in terms of the iterative green s function method can be regarded as corresponding to a very wide strip ( @xmath139 ) , there is thus no finite - size effect in the results . as a result , there is no gap opening in the spectral function at the weyl nodes at @xmath140 for @xmath30 and at @xmath117 for @xmath31 .    [ fig6 ]   ( a , b , c , d , e ) and for @xmath31 ( f , g , h , i , j ) .
results for the left and right edges are assigned respectively with a label of @xmath136 and @xmath137 . in the calculation ,
the width of the strip ( @xmath138 ) is assumed as infinite so that there is no finite - size effect .
parameters are taken as @xmath96 , @xmath94 , @xmath95 , @xmath93 .
the angles used are ( a ) @xmath100 , ( b ) @xmath101 , ( c ) @xmath102 , ( d ) @xmath141 , ( e ) @xmath115 , ( f ) @xmath100 , ( g ) @xmath122 , ( h ) @xmath124 , ( i ) @xmath141 , and ( j ) @xmath115 .
see main text for the definitions of the angles.,title="fig:",width=321,height=566 ]    while the direct observation of the edge states in terms of arpes requires very high resolution or large proximity - induced pairing amplitude in the hm film , the sts spectrum can be used as a more sensitive probe to detect the edge states . for clean samples with ideal edges ,
the sts detects the density of states ( dos ) on the edge of the sample , defined as an integration of the spectral function over the edge momenta in the 1d edge bz @xcite . as shown in fig.6
are the results for both @xmath30 and @xmath31 .
prominent fine structures in the dos are found , which show striking evolutions with @xmath3 and are related directly to the evolutions in the energy spectra of the strip .
there are qualitative differences in the behavior of the dos between those for @xmath30 and those for @xmath31 .
the behaviors of the two are also very different from the results for the heterostructure with an @xmath1-wave superconductor @xcite .
the dos follows closely the evolution of the edge state spectrum .
we emphasize here two most interesting features in the results .
the first point is the striking differences in the dos at @xmath100 for the three choices of @xmath28 .
in contrast to @xmath30 and @xmath29 which both have well - defined zero - energy edge states , the zero - energy edge states for @xmath31 merge with the bulk weyl nodes and therefore have vanishing weight on the edge . this explains the finite dos at @xmath142 for @xmath30 and @xmath29 @xcite and the vanishing zero - energy dos for @xmath31 .
in addition , the edge states for @xmath30 are distorted and have two finite - energy critical points with zero slope [ figs .
4(a ) and 4(b ) ] .
they in turn give the two low - energy van hove singularities in fig.6(a ) .
the two low - energy van hove singularities in fig.6(f ) are also related to the presence of finite - energy critical points with zero slope in the edge states , as is clear from figs.5(a ) and 5(b ) .
the second feature is the existence of a critical angle at which the edge states on one edge becomes completely ( the @xmath1-wave case ) or partly ( the two @xmath0-wave cases ) flat . in this aspect ,
the results for @xmath30 are similar to those for @xmath29 @xcite .
the appearance of a nearly flat portion in the edge state dispersion shown in figs .
4(d ) and 4(f ) manifest as a prominent zero - energy dos peak in figs .
6(b ) and 6(c ) . for @xmath31 ,
a portion of the edge state dispersion for one edge also becomes flat at the transition point from the type - i weyl node to the type - ii weyl node , which would also give a zero - energy van hove singularity in the dos . from figs .
3(h ) and 5(e ) , @xmath127 is very close to but slightly smaller than this transition point .
correspondingly , we get in fig.6(h ) a not fully developed zero - energy van hove singularity .
in addition , the van hove singularity for @xmath31 occurs on the opposite edge as compared to those for @xmath30 and @xmath29 .
this is understood from the fact that , for the edge states close to @xmath117 on the same edge ( e.g. , the left edge with @xmath136 ) and at @xmath100 , their velocity for @xmath31 is opposite in sign to those for @xmath30 and @xmath29 .      here , we summarize the main differences between the energy spectra of the heterostructures with @xmath0-wave form factors [ @xmath30 and @xmath31 ] and those for the heterostructure with @xmath1-wave form factor [ @xmath29 ] @xcite .    for the @xmath1-wave case , there are three different phases as we change the direction of the magnetization , the fully - gapped phase , the phase without a full gap but has a ubiquitous local gap , and the swsc phase with a pair of nodes between the two low - energy quasiparticle bands when the magnetization lies along an in - plane direction . for both @xmath30 and @xmath31 of the @xmath0-wave cases ,
the energy spectra always have nodes and thus describe 2d dwsc .    in terms of the number of weyl nodes
, there is a single pair of weyl nodes in the swsc .
the two nodes have opposite nonzero quasiparticle energies .
for @xmath30 , when a nonzero out - of - plane component of the magnetization exists , there are four weyl nodes .
when the magnetization lies completely along an in - plane direction , an additional pair of weyl nodes emerge and so there are in total three pairs of weyl nodes .
the quasiparticle energies for the nodes are zero for completely out - of - plane magnetization , but are nonzero when an in - plane component of the magnetization exists . for @xmath31 , there are always two pairs of weyl nodes in the quasiparticle energy spectrum .
the energies of the pair of nodes along @xmath146 are fixed at @xmath147 for an arbitrary orientation of the magnetization . for the other pair of nodes along @xmath148
, the energies are zero for purely out - of - plane magnetization but is nonzero in the presence of a finite in - plane magnetization component .    as regards the type of weyl nodes , only type - i weyl nodes are possible for swsc and dwsc with @xmath30 . for @xmath31 , the pair of nodes fixed at zero energy along @xmath146 change between type - i and type - ii as the direction of the magnetization varies .      in this subsection
, we make some comments on the nature of the phases described by eq.(3 ) , for a hm with spin - singlet pairing term and rsoc . for purely out - of - plane magnetization
, the quasiparticle spectrum is similar to a conventional bcs superconductor : a full gap opens for the @xmath1-wave case which is effectively a spinless @xmath149 superconductor @xcite , whereas for the two @xmath0-wave cases the spectrum consists of four isolated nodes and is gapped everywhere else . as we tune the magnetization away from the out - of - plane direction , a nonzero in - plane component of the magnetization ( i.e. , @xmath150 ) cooperates with the rsoc to tilt the quasiparticle spectrum continuously along the @xmath112 direction . for the @xmath1-wave case , along with the transition from the fully - gapped phase to the phase without a full gap but has a ubiquitous local gap , the @xmath67 and @xmath66 quasiparticle bands begin to overlap in energy @xcite .
the ` fermi surface ' of the quasiparticle bands , constituted by states with quasiparticle energy @xmath147 , consists of two disconnected pockets : one ` hole'-like pocket from @xmath66 and one ` electron'-like pocket from @xmath67 .
for the two @xmath0-wave cases , while the spectra have nodes for all magnetization orientations , the energy overlap between the @xmath67 and @xmath66 quasiparticle bands also increases as the magnetization tilts continuously to the in - plane direction . for @xmath30 ,
the ` hole'-like and ` electron'-like pockets are always disconnected , similar to the case of @xmath29 . for @xmath31 , on the other hand , the ` hole'-like and ` electron'-like pockets begin to be connected through the two weyl nodes lying along the @xmath90 axis once they turn into type - ii .
as shown in fig.7 is an illustration of the evolution of the @xmath147 fermi contour explained above , for five typical values of @xmath3 in the range @xmath151 $ ] .
they are interpreted naturally as a series of lifshitz transitions @xcite .
namely , as the magnetization orientation varies , the @xmath147 fermi contour undergoes transitions between four and three topologically distinctive configurations for @xmath30 and @xmath31 , respectively .
similar analysis can be applied to the case of @xmath29 .
it is important to note that the lifshitz transition is usually encountered and discussed in the normal phase of a metal , as an abrupt change in the fermi surface topology of the metal @xcite .
its revelation in the presence of pairing indicates the essential difference between the state of the present system and other more well - known types of pairings .
[ fig7 ]   equal - energy contours for the quasiparticle bands .
@xmath30 for ( a , b , c , d , e ) and @xmath31 for ( f , g , h , i , j ) .
parameters are taken as @xmath96 , @xmath94 , @xmath95 , @xmath93 .
the angles used are ( a ) @xmath100 , ( b ) @xmath101 , ( c ) @xmath102 , ( d ) @xmath141 , ( e ) @xmath115 , ( f ) @xmath100 , ( g ) @xmath122 , ( h ) @xmath124 , ( i ) @xmath141 , and ( j ) @xmath115 .
see main text for the definitions of the angles .
the numberings of the quasiparticle bands to which the pockets belong are indicated in the figures.,title="fig:",width=321,height=566 ]    now , we compare the quasiparticle spectrum for the present state and those for several well - known types of pairings , including the bcs - like pairing @xcite , the fflo pairing @xcite , and the so - called breached pairing ( or , sarma state ) @xcite .
the bcs pairing naturally occurs in a system with a centro - symmetric fermi surface .
the pairing between states of equal and opposite momenta thus gives a spectrum in which the @xmath67 quasiparticle band is completely empty and the @xmath66 quasiparticle band is fully occupied . for our case , the pairing is also between states of equal and opposite momenta .
but the energy spectrum deviates significantly from the conventional bcs spectrum that is symmetric with respect to @xmath147 for each @xmath152 . in the fflo pairing ,
the net momentum of a pair is nonzero .
a simplest case where this phase might be found is when the fermi surface is asymmetric with respect to the center of the bz . in our case , while the fermi surface is asymmetric in the presence of an in - plane component of the exchange field , the net momentum of a cooper pair is always zero .
finally , in the breached pairing state , the pairing gap between _ two _ bands of unequal fermi momenta opens at nonzero energy . for our system , in the presence of an in - plane component of the exchange field and away from @xmath117 , the pairing gap between the @xmath67 and @xmath66 bands also opens at nonzero energy and is similar to the breached pairing state .
however , on one hand the present system consists of only one fermi surface in the normal phase , on the other hand the states on the fermi surface with @xmath153 are usually gapped at @xmath147 . only for @xmath31 and when the two nodes along @xmath146 have turned from type - i to type - ii , would the full fermi surface be gapless and the spectrum be most similar to the breached pairing state .
however , the @xmath147 states in this case consist of both ` electron'-like states from @xmath67 and ` hole'-like states from @xmath66 , rather than from a single quasiparticle band in the case of a standard breached pairing state .    from the above comparisons ,
we see that the appearance of the unique pairing state in the present system is caused by a synergic action of two factors : the bcs - like pairing correlation induced through proximity effect , and the asymmetry of the fermi surface in the hm .
the first aspect of the model is a result of conservation of the parallel component of the momentum through the interface which is assumed to be perfectly translational invariant .
the second factor results from a combination of the in - plane component of the magnetization and the rsoc .    in our treatment of the pairing in the hm
, we have assumed that it arises completely from the proximity coupling with the spin - singlet superconductor .
if instead we assume that there is an intrinsic pairing tendency in the hm , then a self - consistent calculation is necessary in addition to the proximity effect . in this case , in the presence of an in - plane component of the magnetization which drives the fermi surface asymmetric with respect to the bz center , the pairing term obtained in a self - consistent manner could also possibly be the fflo pairing @xcite .
in ref.@xcite focusing on the heterostructures with @xmath1-wave superconductor , we have characterized the phases and phase transitions driven by @xmath3 in terms of four topological quantities : a pfaffian @xmath154 invariant defined based on the particle - hole symmetry of the model , the tknn number defined in terms of the occupied quasiparticle states , the band - wise chern numbers associated with separate nondegenerate quasiparticle bands , and the zak phase which is a berry phase defined along one - dimensional subspace of the 2d bz @xcite .
the band - wise chern numbers in the @xmath1-wave case are motivated from the ubiquitous local gaps in the quasiparticle spectrum , except for @xmath115 and @xmath155 which correspond to the 2d weyl superconductivity phase @xcite . in the case of heterostructure with @xmath0-wave superconductor ,
however , there are always nodes between the two low - energy quasiparticle bands @xmath66 and @xmath67 .
therefore , the band - wise chern number can not be defined for the @xmath0wscs .
in addition , the zero - energy nodes existent for out - of - plane magnetization for @xmath30 and always existent for @xmath31 also cause difficulty in calculating the tknn number .
thus , we focus on the pfaffian @xmath154 invariant and the zak phase to characterize the topological properties of the phases in @xmath0wscs .    before studying the topological invariants
, we list the relevant symmetries ( or , broken symmetries ) of the model defined by eq.(3 ) .
one important aspect of the model is the absence of time - reversal symmetry , @xmath156 , where @xmath157 denotes complex conjugation .
once @xmath158 , no matter what the value of @xmath3 is , we have @xmath159 and thus the time - reversal symmetry is broken . besides , the inversion symmetry , which for the present single - orbital model is simply the operation bringing @xmath152 to @xmath160 , is broken by the rsoc since we have @xmath161 .
however , the model defined by eq.(3 ) always has the particle - hole symmetry , @xmath162 .
that is , @xmath163 holds for all three @xmath28 and for all values of @xmath3 .
the next relevant symmetry is the mirror reflection symmetry . for the mirror reflection with respect to the plane
perpendicular to the @xmath32 axis , the corresponding mirror reflection operator is @xmath164 .
similarly , we have @xmath165 and @xmath166 , for mirror reflections with respect to the planes perpendicular to the @xmath111 axis and @xmath18 axis . for @xmath30 and @xmath29 , we have @xmath167 @xmath168 and @xmath169 for @xmath31 , by redefining the mirror reflection operator for @xmath32 and @xmath111 directions as @xmath170 and @xmath171 , it is easy to see that the bdg hamiltonian @xmath172 satisfies the same transformation properties as the other two pairing channels . from eqs.(16)-(18 ) , the model has the mirror symmetry along @xmath32 only when @xmath173 .
no mirror symmetry along @xmath111 exists , once @xmath158 .
when @xmath174 , @xmath175 acts effectively as an inversion symmetry of the model .
this effective inversion symmetry ensures that each pair of weyl nodes having opposite wave vectors can share the same quasiparticle energy for @xmath100 or @xmath176 . finally , the model has the four - fold rotation symmetry around the @xmath18 axis when @xmath174 , namely @xmath177\tilde{s}_{z}=h(\mathbf{k}).\ ] ] @xmath178 .
@xmath179 for @xmath29 .
@xmath180 for @xmath30 and @xmath31 . as a result of this four - fold rotation symmetry ,
all four weyl nodes share the same quasiparticle spectrum for @xmath100 and @xmath176 .
other symmetries ( e.g. , the so - called magnetic mirror reflection symmetries @xcite ) of the model , while on one hand can be defined as a combination of several of the above symmetries , on the other hand are not directly relevant to the following discussions .
therefore , we do not list them here .      from the particle - hole symmetry of the model defined by eq.(3 )
, the pfaffian @xmath181 invariant for the dwscs can be defined as @xcite @xmath182,\ ] ] where @xmath183 , the function @xmath82 gives the sign of a real number @xmath32 .
for both @xmath30 and @xmath31 , we have @xmath184 .
for the present parameters with large @xmath2 and a single nondegenerate fermi surface , we have @xmath185 , the phase is therefore always nontrivial ( @xmath186 ) from the pfaffian @xmath154 invariant .    if the energy spectrum of the particle - hole symmetric bdg hamiltonian , for example eq.(3 ) with @xmath29 and @xmath100 , is fully gapped , the pfaffian @xmath181 invariant was shown by ghosh _
et al _ to give the parity of the chern number @xcite .
in another work based on a model which can be considered as the @xmath100 case of eq.(3 ) with @xmath30 or @xmath31 , sato and fujimoto pointed out that although the spectrum is gapless , the parity of the chern number is still well - defined and determined by the same formula [ i.e. , eq.(20 ) ] for the pfaffian @xmath181 invariant @xcite .
therefore , for both of the two choices of @xmath28 for the @xmath0wscs and for purely out - of - plane magnetization , the parity of the chern number is well - defined and is odd for the parameters considered in sec.iii . from the topological nature of @xmath187 ,
it is then natural to expect that the system described by eq.(3 ) retains a certain topological character as @xmath3 deviates from @xmath105 and @xmath176 . this is confirmed from the persistence of the edge states shown in sec.iii , except for fig.3(j ) where the edge modes are merged to the bulk bands . from the analysis of the zak phase in sec.iv c ,
however , there should be hidden edge states for the parameters of fig.3(j ) .
the edge states for a general @xmath3 evolve continuously from those for @xmath100 .
so , it is important to define the chirality of the edge states for @xmath100 . in the simplest case of a @xmath149 superconductor ,
the edge states are linearly dispersive and counter - propagating on the two opposing edges .
the chirality of the edge state can be defined in terms of the sign of the velocity of the edge modes .
the hm layer of the heterostructure with @xmath1-wave superconductor is equivalent to a @xmath149 superconductor , the chirality of the edge states can thus be defined in this manner @xcite . in the heterostructures with @xmath0-wave superconductors ,
however , the edge states is curvous and the velocity changes sign several times ( see figs . 2 and 3 ) .
nonetheless , we can still define the chirality of the edge states in terms of the velocity of the edge modes in the following manner .
on an edge with a unique edge state , we divide the full dispersion of it into several segments connected together through points with zero velocity .
each segment is assigned a chirality , positive or negative , equal to the sign of the velocity in this section of the dispersion .
the overall chirality is defined simply as a multiplication of the chiralities of all the segments . in this manner , it is straightforward to see from figs . 4 and 5 in sec .
iii c ( and , fig . 5 in ref.@xcite for the @xmath1-wave case ) that the edge states on the same edge and for @xmath100 have the same chirality for all three @xmath28 .
in addition , as @xmath3 increases , the chirality of the edge states on the left edge ( @xmath136 ) keeps unchanged , whereas the chirality of the edge states on the right edge ( @xmath137 ) change sign .
the @xmath3-dependence of the ( generalized ) chirality of the edge states for all three pairings are thus qualitatively the same .
the zak phase is a berry phase defined as a line integral over the berry connection along a 1d loop of the 2d bz @xcite . in calculation , one component ( i.e. , @xmath90 or @xmath112 ) of the wave vector is fixed as a parameter , integration over the wave vector along the orthogonal direction then traces a 1d loop , as a result of the periodicity of the bz .
for example , by taking @xmath112 as a parameter , the zak phase for the @xmath39-th ( @xmath188 , @xmath41 , @xmath42 , @xmath43 ) quasiparticle band is defined as @xmath189 where @xmath190 is the @xmath32 component of the berry connection for the @xmath39-th quasiparticle band , @xmath191 .
@xmath192 is the eigenvector of the @xmath39-th quasiparticle band .
another set of zak phase , @xmath193 , is defined in the same manner by reversing the roles of @xmath90 and @xmath112 .
when the effective 1d model is inversion symmetric , the zak phase can only take two values ( @xmath105 or @xmath176 ) , otherwise the zak phase can take any value between @xmath105 and @xmath194 @xcite . for the @xmath1wsc defined by @xmath29 and @xmath173 ,
the relevant symmetry is @xmath195 .
the positions of the two weyl nodes of the @xmath1wsc can be captured by a jump in the zak phase between @xmath105 and @xmath176 @xcite .
in addition , on an edge running along @xmath111 ( so that @xmath112 is still a good quantum number but @xmath90 is not ) , there are edge states in the segment of @xmath112 for which the zak phase is @xmath176 @xcite .    from eqs.(16)-(18 )
, the @xmath0wscs with @xmath30 and @xmath31 have the same mirror symmetry ( or , modified mirror symmetry , for @xmath31 ) as @xmath29 . as a result ,
the cases in which the zak phase can take a well - defined quantized value of @xmath105 or @xmath176 are still associated with @xmath115 or @xmath196 , when the magnetization has only an in - plane component . in this case
, the quantization of @xmath197 follows from eq.(16 ) .
in addition , by combining the particle - hole symmetry @xmath198 and @xmath199 , we can see that the model has a ` inversion - like ' symmetry along @xmath112 for fixed @xmath90 .
more exactly , we have @xmath200 therefore , the zak phase defined along @xmath112 loops should also take quantized values ( @xmath105 or @xmath176 ) .
on the other hand , once the out - of - plane component of the magnetization ( @xmath201 ) is nonzero , no ` inversion - like ' symmetry can be defined along an arbitrary 1d loop with one wave vector fixed at a general value . a few exceptions to this conclusion are for @xmath174 and for those loops with the fixed wave vector taking time - reversal - invariant values ( i.e. , @xmath105 and @xmath202 ) .
for these special loops , the 1d effective model has an ` inversion - like ' symmetry enforced by the symmetry of @xmath175 defined by eq.(18 ) .
[ fig8 ] , at @xmath68 .
@xmath30 for ( a ) , and @xmath31 for ( b ) .
the parameters are the same as those used in sec.iii.,title="fig:",width=245,height=321 ]    we have calculated explicitly the zak phases @xmath193 and @xmath197 for both @xmath30 and @xmath31 , at @xmath100 and @xmath68 . in agreement with the above expectation based on symmetry , @xmath193 and @xmath197 take quantized values of @xmath105 or @xmath176 for @xmath68 . for @xmath100
, the zak phases are quantized only when the fixed wave vector takes @xmath105 or @xmath202 .
we show in figure @xmath203 the nontrivial results for the zak phase , which are @xmath197 for @xmath30
[ fig.8(a ) ] and @xmath193 for @xmath31 [ fig.8(b ) ] , both at @xmath68 .
only the zak phases for the two low - energy quasiparticle bands ( namely , @xmath204 ) are shown since the results for the two high energy bands are always zero .
fig.8(b ) does not contain results for @xmath110 since the integration loop for @xmath110 contain two nodes between @xmath66 and @xmath67 which makes the zak phase ill defined . in fig.8(a ) , results for the loops containing the nodes are also excluded .
in addition , we have obtained ( at @xmath68 ) @xmath205 for @xmath30 and @xmath206 for @xmath31 .    comparing fig.8(a ) with the zak phase calculated for the @xmath1wsc @xcite ,
it is clear that they are essentially the same .
although the pairing amplitudes are taken at different values , the positions of the weyl nodes of the @xmath1wsc and the new pair of weyl nodes for the @xmath0wsc with @xmath30 are identical to each other , because they are determined by the vanishing of the effective pairing amplitude . and , though there are two additional pairs of weyl nodes from @xmath207 , the zak phase do not see them .
this is also a very natural result , though , because each node contributes @xmath176 or @xmath105 to the zak phase , so the net change of the zak phase upon crossing the line with two weyl nodes is zero ( modulo @xmath194 ) .
in addition , it is clear that the region with a nontrivial zak phase of @xmath176 is also characterized by the presence of edge modes connecting the two additional pairs of weyl nodes emerging from the emergent mirror symmetry of the model . on the other hand ,
the fact that @xmath205 for all @xmath90 is consistent with the fact that the change in the number of weyl nodes over which the integration loop sweeps is always an even number , zero or two .
therefore , even if each node contributes @xmath176 , the total change in the zak phase defined modulo @xmath194 is still zero .
the zak phase for @xmath31 is however more peculiar .
namely , @xmath206 although the two weyl nodes along @xmath148 are separated from each other and thus the integration loop sweeps only a single weyl node .
this result can be understood from the low - energy effective models for the two nodes along @xmath148 for @xmath31 : they are real symmetric [ see eqs .
( a11)-(a13 ) in appendix a ] and thus the wave function can be taken as purely real .
as such , no nontrivial phase factor can be collected when we traverse the nodes . on the other hand ,
the low - energy effective models for the two weyl nodes of @xmath31 along @xmath146 ( which are type - ii for @xmath68 ) are complex .
this explains why the zak phase @xmath193 shows a jump of @xmath176 when the integration path traverse one node .
@xmath206 is also in agreement with the absence of well - defined edge states in fig.3(j ) for edges running along the @xmath111 direction . on the other hand , from the nontrivial wave vector dependence of @xmath193 , it is natural to expect the presence of edge states on the edges running along @xmath32-direction
this in turn can account for the nontrivial result of the pfaffian @xmath154 invariant .
however , since the bulk quasiparticle spectrum is completely gapless from the @xmath90 direction , the edge states , even if they exist , would merge to the bulk spectra and become undiscernible .
however , it is still possible that the nontrivial zak phase @xmath193 can have some influence on the properties of the system .
we leave the analysis over this interesting possibility to later works .
finally , we consider the case of @xmath100 .
consistent with the symmetry analysis , we have found that the zak phase @xmath197 [ @xmath193 ] for @xmath100 does not take quantized value for a general @xmath112 ( @xmath90 ) . on the other hand , the zak phase corresponding to wave vectors @xmath105 and @xmath202 do take quantized values ( @xmath105 or @xmath176 ) .
these quantized zak phase can be utilized to define a new @xmath154 invariant , @xmath209 where @xmath210 or @xmath111 when considering an edge running along @xmath211 direction . @xmath212 and @xmath42 corresponds to the existence and nonexistence of well - defined edge state at @xmath213 .
this connection can be verified easily by making a comparison between the energy spectra shown in sec.iii and the following numerical results for the zak phase , at @xmath100 : for @xmath30 , @xmath214 and @xmath215 . for @xmath31 , @xmath216 .
for @xmath31 , this result is consistent with the fact that since for both edges running along @xmath32 and @xmath111 , there are always projections of two weyl nodes to the origin of the 1d bz , which act as bulk channels coupling the two edges together , no well - defined edge states can be defined at the center of the one dimensional edge bz .
in the above discussions , we have focused on heterostructures in which the hm has a single nondegenerate fermi surface .
it is interesting to ask what if the spin polarization is not complete , namely what if the two spin - split bands both contribute to the fermi surface .
this is relevant since ferromagnetic materials with partially polarized fermi surface is more commonly found than hm .
this problem is also relevant if we turn our attention from the hm layer to the layer of spin - singlet superconductor .
it is clear that the same model can be used to describe this layer .
the only qualitative difference is just in the magnitude of the exchange field , since the proximity - induced exchange field in the @xmath0-wave superconductor is in general not large enough to split completely one spin component off the fermi surface .
consider a ferromagnetic metal with a fermi surface consisting of two sheets with opposite spin - polarization and different volumes .
the two sheets with opposite spin polarizations would naturally be separated from each other . replacing the hm in the heterostructure studied above with such a partially polarized ferromagnet
, we can expect to apply the previous analyses separately to each spin component .
as such , we can expect to get two sets of low - energy effective models with similar properties but related separately to the two spin - split bands .
in particular , for each band , there will be a pair of counter - propagating chiral majorana edge states .
it is then interesting to ask whether or not the two sets of chiral edge modes are always separated from each other . or , in cases where they coexist in some parameter ranges , how would they interact with each other .
( a ) @xmath68 which corresponds to an @xmath1wsc , results for the two edges are the same .
( b ) and ( c ) are spectral functions on the left ( @xmath136 ) and right ( @xmath137 ) edges of the strip for @xmath100 .
the parameters used are @xmath91 , @xmath92 , @xmath94 , @xmath217 , @xmath218 .
the darker the color , the larger the spectral weight.,title="fig:",width=245,height=472 ]    to answer the above question , we have studied the changes in the properties of the weyl superconductivity phase brought by increasing the chemical potential so that it crosses both of the two spin - polarized subbands of the ferromagnet . as shown in fig.9(a ) are the spectral functions for the edge layers ( results for the two edges are the same ) of a strip for @xmath29 and @xmath68 , namely the @xmath1wsc . while on one hand the edge modes for the two spin species keep largely independent of each other , a hybridization gap opens at ( @xmath117 , @xmath142 ) .
besides , because the fermi sheets of the two spin channels open gaps independently , the edge modes in the non - overlapping regions of the pairing gaps of the two pin - split bands will be masked by the bulk states of the opposite - spin band .
overall , the chiral edge modes appear in _ three disconnected _ regions of the local energy gap of the energy spectrum . in fig.10
, we show the spectral functions on the two edges of a strip of the @xmath0wsc at @xmath100 for @xmath30 [ 10(a ) and 10(b ) ] and @xmath31 [ 10(c ) and 10(d ) ] .
the results for @xmath29 at @xmath100 are shown as figs.9(b ) and 9(c ) . for @xmath30 ,
we observe similar coexistence of two chiral edge states and the opening of a hybridization gap at ( @xmath117 , @xmath142 ) between the two edge states belonging separately to the two spin - subbands . for @xmath31 ,
the two copies of edge states simply superimpose together and no hybridization gap opens . in ref.@xcite , it was pointed out that the topological phase with a single majorana edge state can be obtained for vanishingly small exchange field when the chemical potential is at the band bottom ( @xmath219 ) . while that case still has a single band crossing the chemical potential , our present case has two bands contributing to the fermi surface . from the numerical results in fig.10 and comparing it with fig.9 for the @xmath1-wave case , we see that while the edge states exist when the fermi surface consists of two sheets , there is a pair of counter - propagating edge modes on each edge which usually hybridize and open a gap .
the physical picture for the @xmath0-wave case is therefore of no essential difference from the @xmath1-wave case as regards the number of edge states and their hybridization .
since there are even number of edge states on each edge , they are not topologically protected .
this is consistent with the pfaffian @xmath154 invariant defined in eq.(20 ) , which gives the topologically trivial result of @xmath220 for the present set of parameters .
the phase on the superconductor side of the heterostructure is therefore topologically trivial .
topologically nontrivial physics are located in the hm layer .
[ fig10 ]   ( a , b ) and @xmath31 ( c , d ) . the angle is taken as @xmath100 .
the parameters are chosen as @xmath91 , @xmath92 , @xmath94 , @xmath95 , @xmath218 .
darker color means larger spectral function.,title="fig:",width=321,height=264 ]
in summary , we have predicted the realization of two - dimensional weyl superconductivity phases in a heterostructure consisting of a half - metal ( or a highly spin - polarized ferromagnet ) clamped in between an insulating substrate and a spin - singlet superconductor .
when the superconductor has @xmath0-wave symmetry ( with a symmetry factor of @xmath221 or @xmath222 ) , the two - dimensional weyl superconductivity phase is realized for an arbitrary magnetization direction in the half - metal .
quasiparticle energy spectra are studied which exhibit interesting variations with the direction of the magnetization . in particular , we have predicted the transition between type - i and type - ii weyl nodes by changing the magnetization orientation .
in addition , the evolution of the phase with the magnetization orientation is unveiled through the lifshitz transitions in the zero - energy quasiparticle spectrum .
the phases are shown to be qualitatively different from several well - known types of pairings .
we characterize the topology of phases by the pfaffian @xmath154 invariant and the zak phase .
this work was supported by the texas center for superconductivity at the university of houston and the robert a. welch foundation ( grant no .
l.h . would also like to acknowledge the support from the china scholarship council .
the numerical calculations were performed at the center of advanced computing and data systems at the university of houston .
in this section , we perform an analysis on the weyl nodes that are supported in the present system defined by fig . 1 and eq.(3 ) of the main text . in particular ,
as fig.3 has shown the transition from type - i weyl node to type - ii weyl node , it is interesting to investigate the possibility of realizing type - ii weyl nodes in other cases . from the analysis presented below
, the type - ii weyl nodes found in fig.3 for @xmath31 are the only cases where type - ii weyl nodes can be realized in the proposed system .
we also point out the possibility of realizing semi - weyl nodes .    the low - energy effective model for a single weyl node can be written generally in terms of the relative momenta @xmath223 as @xcite @xmath224 where @xmath225 is a label for the specific weyl node under consideration
, @xmath226 is the second - order unit matrix .
the two eigenvalues of the model are @xmath227 @xmath228 describes a type - ii weyl node if and only if there exists a direction @xmath229 along which @xmath230 .
otherwise , the weyl node is type - i .    for the heterostructure with @xmath1-wave superconductor , the weyl superconductivity phase ( @xmath1wsc ) appears only when the magnetization lies along in - plane directions ( @xmath68 and @xmath231 ) .
from the low - energy effective model of the weyl nodes for @xmath68 @xcite , the weyl nodes can become type - ii only when the strength of the rsoc ( @xmath15 ) is of the same order of magnitude as the hopping amplitude @xmath12 . since this is unrealistic
, the weyl nodes in the @xmath1wsc can only be type - i .    for the @xmath0wsc with @xmath30 ,
the two weyl nodes emerging for in - plane magnetization are qualitatively the same as the weyl nodes for @xmath1wsc .
following the same analysis as above , they can only be type - i for realistic model parameters .
the remaining four weyl nodes exist for all @xmath3 angles .
their positions in the 2d bz are nearly @xmath3-independent and are @xmath232 ( @xmath47 , @xmath233 ) .
@xmath98 is defined by eq.(10 ) . for a general @xmath3
, we have @xmath234 @xmath235 @xmath236 and @xmath237 the constant term in @xmath238 gives the energy of the weyl node but is irrelevant to the type of the weyl node . @xmath238 only gives a velocity along the @xmath239 direction , of the order no larger than @xmath15 .
@xmath240 , on the other hand gives a velocity of the same order as @xmath12 along @xmath239 . as a result ,
the four weyl nodes @xmath241 ( @xmath47 , @xmath233 ) are type - i for realistic model parameters and all @xmath3 values .    for the @xmath0wsc with @xmath31 ,
there are always four weyl nodes located approximately at @xmath242 ( @xmath47 , @xmath233 ) .
to simplify notation , we introduce the following abbreviations : @xmath243 , @xmath244 , @xmath245 , @xmath246 , @xmath247 , @xmath248 . for a general @xmath3 , the effective low - energy model for the weyl nodes
are @xmath249 @xmath250 @xmath251 and @xmath252 again , because @xmath238 only gives the quasiparticle a velocity along @xmath239 , only if the velocity along @xmath239 from @xmath253 is smaller than that from @xmath238 , can the weyl node become type - ii .
note that @xmath254 , the effective model for the two weyl nodes @xmath255 becomes @xmath256 @xmath257 @xmath258 and @xmath259 . from the above effective model , the two weyl nodes at @xmath255
can only be type - ii when @xmath15 is of the same order as @xmath12 , which is unrealistic in the present setting .
an interesting critical case for these two nodes is when the coefficient of @xmath260 in eq.(a12 ) becomes zero . when this happens , the linear in @xmath260 term vanishes , we have to add the term quadratic in @xmath260 .
the nodes thus become linear in @xmath239 but quadratic in @xmath260 , and is thus semi - weyl @xcite . in a hm with large magnetic anisotropy energy , this critical point corresponding to a general @xmath3 angle should be very difficult to attain . if the magnetic anisotropy energy is small so that we can tune @xmath3 through the full @xmath194 range , it is still possible to turn the two weyl nodes at @xmath255 to be semi - weyl .
the effective model for the other pair of weyl nodes , @xmath261 , are @xmath262 @xmath263 @xmath264 and @xmath265 . for experimentally relevant parameters satisfying @xmath266 , @xmath267 , and @xmath268 ,
the type - ii weyl node is clearly achievable . in the limiting cases of @xmath269 and @xmath270
, we can also get an analytical expression for the critical angles marking the transition between type - i and type - ii weyl nodes , @xmath271    overall , we have verified both numerically and analytically that a transition between type - i and type - ii weyl nodes can be achieved in our proposed experimental configuration , related to the pair of weyl nodes lying along @xmath146 for @xmath31 .
finally , we reemphasize that the two cases , @xmath30 and @xmath31 , can be realized by growing the same @xmath0-wave superconductor on the hm thin film along two inequivalent special directions .
guoqing chang , su - yang xu , daniel s. sanchez , shin - ming huang , chi - cheng lee , tay - rong chang , guang bian , hao zheng , ilya belopolski , nasser alidoust , horng - tay jeng , arun bansil , hsin lin , and m. zahid hasan , sci . adv . * 2 * , e1600295 ( 2016 ) .