calorons are instantons at finite temperature . for a long time the influence of a background polyakov loop on the properties of these topological excitations has been neglected .
solutions were constructed long ago  @xcite and were studied in detail in the semi - classical approximation  @xcite . in all these studies
the polyakov loop at spatial infinity ( also called the holonomy ) was trivial , i.e. an element of the center of the gauge group . that the influence of the background polyakov loop on the topological excitations can be dramatic is particularly clear in the confined phase , where on average its trace vanishes .
caloron solutions in such backgrounds were constructed only relatively recently  @xcite and can be seen as composed of massive monopole constituents with their magnetic charges adding to zero .
it was observed that the one - loop correction to the action for configurations with a non - trivial asymptotic value of the polyakov loop gives rise to an infinite action barrier , which were therefore considered irrelevant  @xcite .
however , the infinity simply arises due to the integration over the finite energy _ density _ induced by the perturbative fluctuations in the background of a non - trivial polyakov loop  @xcite .
the proper setting would therefore rather be to calculate the non - perturbative contribution of calorons ( with a given asymptotic value of the polyakov loop ) to this energy density , as was first successfully implemented in supersymmetric theories  @xcite , where the perturbative contribution vanishes .
it has a minimum where the trace of the polyakov loop vanishes , i.e. at maximal non - trivial holonomy .    in a recent study at high temperatures , where one presumably can trust the semi - classical approximation ,
the non - perturbative contribution of these monopole constituents ( also called dyons ) was computed  @xcite . when added to the perturbative contribution  @xcite with its minima at center elements , a local minimum develops where the trace of the polyakov loop vanishes , deepening further for decreasing temperature .
this gives support for a phase in which the center symmetry , broken in the high temperature phase , is restored and provides an indication that the monopole constituents are the relevant degrees of freedom for the confined phase .
also lattice studies , both using cooling  @xcite and chiral fermion zero - modes  @xcite as filters , have now confirmed that monopole constituents do dynamically occur in the confined phase . a charge 1 caloron is seen for su(@xmath0 ) to consist of @xmath0 constituent monopoles .
in the deconfined phase , due to the fact that the average polyakov loop becomes a center element , the caloron returns to the form known as the harrington - shepard solution  @xcite .
the latter can also be interpreted as consisting of constituent monopoles , however , with @xmath1 of them being massless .    to be precise , for self - dual configurations in the background of non - trivial holonomy the masses of constituent monopoles are given by @xmath2 , with @xmath3 .
the @xmath4 are related to the eigenvalues of the polyakov loop at spatial infinity , = _
|x|(_0^a_0(t , x)dt)=g^(2i(_1,_2,  ,_n))g , ( this expression assumes the periodic gauge @xmath5 ) where @xmath6 is the gauge rotation used to diagonalize @xmath7 and @xmath8 is the period in the imaginary time direction , related to the inverse temperature .
the eigenvalues @xmath9 are to be ordered on the circle such that @xmath10 , with @xmath11 and @xmath12 , which guarantees that the masses add up to @xmath13 , the instanton action per unit ( imaginary ) time . at higher topological charge @xmath14 ,
the parameter space of widely separated constituent monopoles is described by @xmath15 monopole constituents , @xmath14 of each of the @xmath0 types of abelian charges ( with overall charge neutrality ) .
we established in an earlier paper  @xcite that well - separated constituents act as point sources for the so - called far field ( that is far removed from any of the cores ) .
when constituents of opposite charge ( @xmath0 constituents of different type ) come together , the action density no longer deviates significantly from that of a standard instanton .
its scale parameter @xmath16 is related to the constituent separation @xmath17 through @xmath18 . yet
, the gauge field is vastly different , as is seen from the fact that within the confines of the peak there are @xmath0 locations where two of the eigenvalues of the polyakov loop coincide  @xcite , thus in some sense varying over the maximal range available ( e.g. for su(2 ) from @xmath19 to @xmath20 ) , whereas for trivial holonomy only one such location is found .    on the other hand ,
when constituents of equal charge come together typically an extended core structure is found .
this was deduced , in particular for the case of charge 2 calorons , from our ability to analytically determine the zero - mode density ( summed over the two zero - modes implied by the index theorem ) in the far field limit , neglecting exponential contributions , such that the range of the exponential contributions shrinks inversely proportional with the temperature . ] . in this limit
it forms a singular distribution on a disc bounded by an ellipse , but approaches two delta functions for well - separated like - charge constituents .
this zero - mode density only sees constituents of a given charge , depending on the boundary condition for the fermions in the imaginary time direction , which can be chosen to be a u(1 ) phase ( containing the physically relevant choice of anti - periodic boundary conditions for thermal field theory ) .
we will show for su(2 ) that their difference for periodic and anti - periodic boundary conditions coincides _ exactly _ with the ( abelian ) charge distribution extracted from the gauge field in the same limit , making contact with an old result due to hurtubise  @xcite for the asymptotic behaviour of the monopole higgs field .
we found two particular parametrizations within the su(2 ) charge 2 moduli space that exhibited these extended charge distributions .
the first includes as a limit arbitrary charge 2 monopoles .
the second of these parametrizations contains as a limiting case the axially symmetric configurations constructed for arbitrary charge in ref .
deforming away from the axial configuration the two discs overlap . describing the intricate behaviour for the non - abelian core of these configurations in this region of the parameter space requires one to find exact solutions , which are presented here .
we rely on early work of nahm  @xcite and panagopoulos  @xcite for charge 2 magnetic monopoles , which is simplified to some extent by our formalism that uses the relation between the fourier transformation of the adhm construction ( as relevant for the finite temperature case ) and the nahm transformation , a crucial ingredient for our success to find explicit solutions  @xcite .
[ fig : fig1 ] gives a particular example for the action and zero - mode densities of a charge 2 caloron solution .
the two dimensional zero - mode basis is chosen such that each zero - mode only localizes on one of the constituents of a given charge , showing both the zero - modes with periodic and anti - periodic boundary conditions in the imaginary time direction .
this paper is organized as follows . in sect .  [
sec : con ] we will outline the construction , introduce the green s function that is computed through the analogy of an impurity scattering problem , and summarize the various limits that can be formulated before explicitly solving for the green s function . in sect .
[ sec : exact ] we present the method that allows one to find the exact solution for the green s function , first for the general case and then applied in more detail to that of topological charge 2 calorons .
readers only interested in the results could skip sects .
[ sec : lim ] and [ sec : exact ] . in sect .
[ sec : plots ] we discuss the two classes of configurations in the moduli space of the charge 2 calorons and provide plots of the various quantities to illustrate the properties of the exact results . in sect .
[ sec : higgs ] we discuss the relation between the algebraic tail of the gauge field and the zero - mode density .
we end with some discussions .
an appendix presents a new result for the limiting behaviour of the action density .
there are two steps in the construction of charge @xmath14 caloron solutions .
the first step involves finding a @xmath21 gauge field @xmath22 , which satisfies the self - duality equation , i.e. the nahm equation  @xcite , on a circle parametrized by @xmath23 , with @xmath23 introduced through replacing the original su(@xmath0 ) gauge field @xmath24 by @xmath25 . although not affecting the field strength , this changes the holonomy to @xmath26 , revealing that @xmath23 has period @xmath27 .
the index theorem guarantees the existence of @xmath14 zero - modes @xmath28 which satisfy the dirac equation , or in the two - component weyl form    @xmath29 ( @xmath30 are the usual pauli matrices )
. we may remove @xmath23 from the gauge field @xmath24 by transforming the zero - mode to @xmath31 , which is at the expense of making the zero - mode only periodic up to a phase factor , @xmath32 . in a similar way we could introduce @xmath33 through @xmath34 , which replaces in @xmath24 by @xmath35 , where @xmath36 . assuming the @xmath14 zero - modes @xmath37 to be orthonormal one has a^ab_(z)=^(a)(x;z , z)^^(b)(x;z , z ) d^4x,[eq : ndual ] or equivalently ( demonstrating as well that @xmath38 does not depend on @xmath33 ) a^ab_0(z)=^(a)(x;z)^ ^(b)(x;z ) d^4x , a^ab_k(z)=2i^(a)(x;z)^ x_k^(b)(x;z)d^4x .
we have shown how @xmath22 can alternatively be related to a fourier transformation of the adhm construction of instantons  @xcite , that periodically repeat ( up to a gauge rotation with @xmath7 ) in the imaginary time direction , so as to turn an infinite charge instanton in @xmath39 to one of finite charge and finite temperature .
the derivation of this relation will not be repeated here , see ref .
@xcite for the details .
the connection to the adhm construction has been useful to simplify the second step in the construction of the caloron solutions , namely how to reconstruct the original gauge field when given a solution to the nahm equation  @xcite ( which is equivalent to the quadratic adhm constraint ) , a_j(z)+[a_0(z),a_j(z)]+_jk[a_k(z ) , a_(z)]=2i_m(z-_m)_m^j,[eq : nahmeq ] for convenience of notation we will henceforth use the classical scale invariance to set @xmath40 .
the singularities in the nahm equation appear precisely for those values where @xmath41 has one of its eigenvalues equal to 1 , at @xmath42 .
this is where some fermion field components become massless , i.e. the zero - mode becomes delocalized , whereas for generic values of @xmath23 it is exponentially localized , which has turned out to be a useful tool to pinpoint the constituent monopoles .
one could apply the nahm transformation again , introducing @xmath43 ( with @xmath44 ) and find the @xmath0 chiral fermion zero - modes in the background of this gauge field .
the construction is somewhat complicated due to the presence of the singularities , whose structure is determined from the matrices @xmath45 which appear in the nahm equation .
not all @xmath46 are independent , e.g. integrating and tracing the nahm equation yields the conditions , @xmath47 .
further constraints are implied  @xcite by the fact that one may introduce ( as is most easily seen in relation to the adhm construction ) @xmath14 two - component spinors @xmath48 in the @xmath49 representation of su(@xmath0 ) , such that 2_a^p_m_b=_0s_m^ab-_m^ab,[eq : srho ] where @xmath50 are projections defined through @xmath51 ( @xmath52 will appear in ) .
however , with reference to the adhm construction , there is great benefit in first finding the solution for the green s function , @xmath53 , where \{-+v(z;x)}f_x(z , z)=4 ^ 2 _ k ( z - z),[eq : green ] with & & v(z;x)4 ^ 2r^2(z;x)+2_m(z-_m)s_m , s_mg(_m)s_mg^(_m ) , + & & r_j(z;x)x_j-(2i)^-1g(z)a_j(z)g^(z ) , [ eq : vrsdef ] and @xmath54 playing the role of  impurities " .
this is formulated in the gauge where first we transform @xmath55 to a constant ( diagonal ) matrix , @xmath56 , as is always possible in one dimension , and then use g(z)(2i(_0-x_0_k)z),_0=0 ( when @xmath57 it is absorbed in a shift of @xmath58 ) to transform @xmath59 to zero .
this is at the expense of introducing periodicity up to a gauge transformation ; although @xmath60 is periodic in @xmath23 and @xmath61 with period 1 ( for @xmath40 ) , @xmath62 no longer is plays the role of the holonomy associated to the dual nahm gauge field @xmath22 . ] .    given a solution for the green s function
, there are straightforward expressions for the gauge field  @xcite ( only involving the green s function evaluated at the  impurity " locations ) and the fermion zero - modes  @xcite .
for the zero - mode density this gives _
z^(a)(x)^_z^(b)(x)=-(2)^-2_^2 f_x^ab(z , z).[eq : zmdens ] in this paper we will only have need for the green s function at @xmath63 , which formally can be expressed as f_x(z , z)=-4 ^ 2((_2k-_z)^-1)_12,_z g^(1)_z^z+1dw , where the @xmath64 component on the right - hand side of the first identity is with respect to the @xmath65 block matrix structure .
in particular this leads to a compact expression for the action density  @xcite ( x)-f_^2(x)=-_^2_^2(ie^-ix_0(_2k-_z)),[eq : sdens ] which can be shown to be independent of the choice of @xmath23 .
the formal expression for @xmath66 can be made explicit by a decomposition into the  impurity " contributions @xmath67 at @xmath68 and the  propagation " @xmath69 between @xmath70 and @xmath71 . for @xmath72
this gives _ z = w_m(z,_m)t_mh_m-1t_2h_1t_1g^(1)h_nt_n h_n-1t_m+1h_m w_m(_m , z ) , with t_m , w_m(z , z ) ^-1.[eq : twdeg ] the columns of the two @xmath73 matrices @xmath74 , defined for @xmath72 , form the @xmath75 solutions of the homogeneous green s function equation , ( -4 ^ 2r^2(z;x))v(z)=0 , [ eq : homg ] of which those in @xmath76 are exponentially rising and those in @xmath77 are exponentially falling .
this implies  @xcite @xmath78 for @xmath79 , in which @xmath80 can be arbitrary non - singular ( to ensure a complete set of solutions ) matrices . in ref .
@xcite we put @xmath81 , but here we find it convenient to leave this choice open . with the  impurity
" scattering problem solved , constructing the exact solutions of the homogeneous green s function equation is the last step in finding analytic expressions for the higher charge calorons .
nevertheless , approximate solutions can be derived , either assuming @xmath82 is far removed from any core such that the gauge field has become abelian ( which we called the far field limit , denoted by a subscript  ff " ) , or assuming @xmath82 and the constituents of type @xmath83 ( belonging to the @xmath83-th interval ) are well separated from all others , but not necessarily from each other ( which was called the zero - mode limit , denoted by a subscript  zm " ) .
this is because we found  @xcite that for @xmath72 the @xmath14 zero - modes @xmath84 only  see " the constituents of type @xmath83 . for @xmath85
we have f^_x(z , z)&=&(f_m^-(z)f_m^-(_m+1)^-1-f_m^+(z ) f_m^+(_m+1)^-1z_m+1 ^ - ) + & & -3 mm ( f_m^-(_m)f_m^-(_m+1)^-1-z_m^+f_m^+(_m ) f_m^+(_m+1)^-1z_m+1 ^ -)^-1 + & & ( f_m^-(_m)f_m^-(z)^-1-z_m^+f_m^+(_m)f_m^+(z)^-1 ) r_m^-1(z),[eq : zmlimit ] up to _ exponential _ corrections in the distance to the constituents of type @xmath86 , with z_m^-_k-2_m^-1r_m-1(_m ) , & & z_m^+_k-2 _ m^-1r_m(_m),[eq : zedefin ] + r_m(z)(r_m^+(z)+ r_m^-(z)),&&_mr_m^-(_m)+r_m-1^+(_m)+s_m.and 2r_m^(z)(f_m^(z))f_m^(z)^-1 .
[ eq : rmdef ] in the zero - mode limit @xmath87 and @xmath88 approach @xmath89 up to _ algebraic _ corrections .
neglecting these contributions as well , e.g. by sending the constituents of type @xmath86 to infinity , leads to the so - called monopole limit ( denoted by a subscript  mon " ) , further simplifying the expression for the green s function to f^_x(z , z)=-u(z,_m+1)u^-1_m(_m,_m+1)u_m(_m , z ) r_m^-1(z),[eq : zm - alg ] with satisfies the green s function equation with boundary conditions @xmath90 and @xmath91 , for @xmath92 . ]
u_m(z , z)f_m^+(z)f_m^+(z)^-1-f_m^-(z)f_m^-(z)^-1.[eq : um ]    in turn , from one may derive the far field limit , giving up to _ exponential _ corrections in the distance to _ all _
constituents f_x^(_m,_m)=2_m^-1 , as is relevant for the expression of the gauge field in this limit  @xcite .
for the zero - mode density ( ) we may use for @xmath93 ( @xmath23 strictly different from @xmath94 ) f_x^(z , z)=r_m^-1(z).[eq : zmff ] the discontinuity in this limit , @xmath95 , arises due to the zero - mode developing a massless component when @xmath23 approaches @xmath70 .
it might seem that , combined with , is inconsistent with an exponential decay .
however , it turns out that  @xcite _ m(x)(4)^-1(r_m^-1(z ) ) [ eq : vdef ] is independent of @xmath23 in the @xmath83-th interval and harmonic everywhere ( hence giving vanishing zero - mode density ) except for some singularities in the cores of the constituents of type @xmath83 .
it is this feature , and our ability to compute @xmath96 exactly for su(2 ) charge 2 calorons , that allowed us to make statements about the localization of the cores , without solving the green s function exactly .    in appendix
a we derive the following new result for the monopole limit of the action density ( ) . in the limit where @xmath82 and the constituent locations of type @xmath83 are well separated from all other constituents , for which both the action and zero - mode densities are static , we find ( with @xmath97 defined in ) ^(x)=-_i^2_j^2,[eq : szm - alg ] up to _ algebraic _ corrections .
this is a direct generalization for the action density of a single bps monopole , @xmath98 $ ] ( located at the origin and with mass @xmath99 ) .
we emphasize that this result can be used irrespective of the distance between the constituents of the same type @xmath83 .
therefore gives in terms of @xmath74 a closed form expression for the multi - monopole energy density .
the same holds , when using , for the zero - mode density ( ) . for the caloron it involves taking the limit where all constituents of type @xmath100 are sent to infinity , which is why this is called the monopole limit .
we recall , that in the far field one should use  @xcite ( see also app .
a ) ^(x)=-_i^2_j^2_m=1^n(f_m^+ ( _ m+1)f_m^+(_m)^-1r_m^-1(_m)_m).[eq : sff ]
finding the exact homogeneous solutions of the green s function equation , , closely follows nahm s method  @xcite to construct the dual chiral zero - modes .
the main advantage of our approach is that we need not worry about boundary conditions , as this is solved by the  impurity " scattering formalism  @xcite . in the following
we restrict ourselves to a given interval @xmath72 and work in the gauge where @xmath101 .    using that @xmath22 is self - dual ,
a consequence of the nahm equation , one easily shows that @xmath102 , such that it is natural to consider the equation d_x(z)=_d_x^(z)=(-2r(z;x))(z)=0.[eq : duald ] it follows that @xmath103 would be a homogeneous solution of the green s function equation , albeit that @xmath103 is a spinor ( with a chirality opposite to that for the zero - modes involved in the nahm transformation , cmp . ) .
we follow nahm  @xcite and use the ansatz @xmath104 , where @xmath105 is a @xmath14 dimensional complex vector , @xmath106 is a unit vector that does not depend on @xmath23 and @xmath107 is an arbitrary normalized constant spinor ( as long as it is not annihilated by @xmath108 ) .
it then follows that @xmath109 satisfies .
the unit vector @xmath106 is found from a complex vector @xmath110 which squares to 0 , @xmath111 , implying its real and imaginary parts are perpendicular and of equal length ( @xmath112 as long as @xmath113 ) , such that ( for ease of notation the @xmath82 dependence of @xmath114 and @xmath115 will henceforth be left implicit ) u = iyy^/(yy^ ) [ eq : udef ] is well defined and @xmath116 , i.e. re@xmath117 , im@xmath117 and @xmath115 form an orthogonal set of vectors . using the ansatz for @xmath103 and introducing ( z)-yr(z;x ) , ( z)-2ur(z;x ) , leads to the equations ( z)v(z)=0,v(z)=(z)v(z)[eq : yu ] for which the first one can only have a solution provided @xmath118 .
it is the great beauty of nahm s formalism that @xmath119 is a conserved quantity .
that is , @xmath120 $ ] ( the nahm equation restricted to an interval and in the gauge where @xmath121 ) implies @xmath122 for any choice of @xmath82 and @xmath114 on a null - cone in @xmath123 ( or rather in @xmath124 since we may rescale @xmath114 with a non - zero complex factor without changing the equations ) .
the conserved quantities are generated by the symmetric traceless monomials , @xmath125 , built from @xmath126 with @xmath127 arbitrary , as one readily verifies .
for example , @xmath128 is constant and defines the center of mass . a natural way to project on the traceless symmetric monomials is precisely through introducing @xmath129 on a null - cone , forming @xmath130 .
it is interesting to note that @xmath131 is always a spherical harmonic of order @xmath127 , used in ref .
@xcite to show through the multipole expansion of @xmath132 that it is conserved and harmonic , except for singularities in the core of the constituents .
once it is established that @xmath119 is independent of @xmath23 for any choice of @xmath82 , we can look for its zeros . using the null - cone parametrization @xmath133 , and the fact that the matrix @xmath134 is @xmath14 dimensional
, we obtain a polynomial equation in @xmath135 of order @xmath75 and hence there are for generic @xmath82 exactly @xmath75 solutions .
it is useful to note that these solutions come in complex conjugate pairs , where the symmetry @xmath136 implies @xmath137 and @xmath138 ( giving @xmath139 up to a multiple , equivalent to @xmath139 in @xmath124 ) .    given a particular zero @xmath114 , we may conveniently write a vector in the kernel of @xmath134 as  @xcite @xmath140 for a fixed choice of @xmath141 , where @xmath142 is the matrix formed by the minors of @xmath134 . the nahm equation is easily seen to imply @xmath143 $ ] for _ any _ choice of @xmath114 on the null - cone . from this one
derives that is such that @xmath144 , in which case @xmath145 and therefore @xmath146\y^{-1}\det\y+\y^{-1}\ddz\det\y=[\u,\adj\y ] + \y^{-1}\tr(\y^{-1}\ddz\y)=[\u,\adj\y]$ ] .
observe that @xmath147 is analytic in @xmath114 , such that the result is valid also when @xmath114 leads to @xmath148 . ]
@xmath149 $ ] .
substituting @xmath150 into gives ( ( z))_ac=(z)((z)(z ) ) _ ac .
using the fact that @xmath151 , we get ( z)=,=- , [ eq : mu ] where the equation for @xmath152 is the same for any value of @xmath153 ( it may depend on the value of @xmath141 ) .
this is useful for studying the asymptotic behaviour of the solution . for large @xmath154
, @xmath155 implies that @xmath156 , such that ( cmp . ) latexmath:[$\vec
u\to\pm\vec x/    guarantee that there are @xmath14 zeros @xmath158 with the sign of @xmath159 negative , leading to solutions that rise as @xmath160 .
it then follows that the @xmath14 zeros @xmath161 give @xmath162 , leading to the solutions that decay as @xmath163 for large @xmath154 .
hence we may put @xmath164 and @xmath165 . defining m_ab^+(z)=-(y^(b)r(z;x ) ) _ ac , m_ab^-(z)=-(y^(b)*r(z;x ) ) _ ac,[eq : mdef ] which are _ algebraic _ in @xmath166 and @xmath167 , we find f^+_ab(z)=m^+_ab(z)^(b)(z),f^-_ab(z)=m^-_ab(z)^(b+k)(z ) , where @xmath168 contains the exponential dependencies , and thus seemingly all the information about the cores of the constituents . to make this more precise we compute @xmath169 , see , using that implies ( with @xmath170 ) @xmath171 , finding that the factors @xmath172 _ drop out _ r_ad^(z)=-_b , e=1^ku^(b)r_ae ( z;x)m_eb^(z)(m^(z)^-1)_bd.[eq : alg ] this proves that @xmath173 is purely algebraic in @xmath166 and @xmath167 , as are @xmath174 and @xmath175 , which determine the far field limit for the zero - mode density and the gauge field .
one might wonder how , given that the @xmath176 are of the
wrong " chirality in the context of the nahm transformation , one could use these results to reconstruct the gauge field for magnetic monopoles where the relation to the adhm construction is not readily available . for
this one observes that the columns of @xmath176 can be used to form a @xmath177 matrix @xmath178 .
using that @xmath179 , one finds @xmath180 .
thus the columns of @xmath181 give @xmath75 independent solutions for each interval , from which @xmath0 normalizable solutions @xmath182 should remain after imposing the appropriate boundary ( cq .
matching ) conditions .
these are then used in nahm s original construction to compute the gauge field ( cmp . )
a^pq_(x)=^(p)(z;x_0,x)^^(q)(z;x_0,x ) dz , where @xmath183 .
there seems to be considerable advantage in using the green s function ( fourier transformed adhm ) method , since it can solve the matching conditions without relying on the availability of exact solutions for the normalizable dual zero - modes . to go beyond the approximations discussed in sect .
[ sec : lim ] and resolve the constituent cores we need to solve for @xmath152 .
this can not always be done in closed form , but it is given by an explicit integral which can be performed numerically when required .
for charge 2 monopoles panagopoulos  @xcite was , however , able to find the exact integral .
we can use the same ingredients for the caloron case and explicitly solve for the green s function in the case of charge 2 calorons .      for charge 2
the number of invariants associated to the conserved quantities of the nahm equation is 8 , of which @xmath184 are related to the center of mass for the constituents of given magnetic charge , coming from the interval under consideration .
assuming now that @xmath166 is traceless , 5 invariants remain , given in terms of the symmetric traceless tensor m_ij-((a_i(z)a_j(z ) ) -_ij(a_k^2(z ) ) ) .
three of its parameters are associated to the rotation @xmath185 which diagonalizes the @xmath186 matrix , @xmath187 , where @xmath185 is fixed by requiring @xmath188 .
the @xmath189 add to zero and can be expressed in terms of the so - called scale ( @xmath190 ) and shape ( @xmath191 ) parameters , c_1=d^2,c_2=d^2 , c_3=d^2 .
the nahm equation for the case of charge 2 can be solved completely in terms of jacobi elliptic functions  @xcite , which was summarized in ref .
@xcite g(z)a_j(z;a,,h , d,)g^(z)2i a_j_2 + i d_jbf_b(d(z - z_0))h^_bh,[eq : jnahm ] where @xmath192 is a global gauge parameter and f_1(z),f_2(z),f_3(z),. [ eq : fs ] with implicitly defined by the elliptic integral of the first kind @xmath193 .
] @xmath194 , @xmath195 and @xmath196 the standard jacobi elliptic functions .
this does not yet address the matching of @xmath166 on the different intervals , where some difference between the monopole and caloron application appears .
for the caloron , apart from the axially symmetric solutions constructed in ref .
@xcite , we found two sets of non - trivial solutions that interpolate between overlapping and well - separated constituents .
it is for these classes of solutions that we will resolve the cores , when constituents overlap and the non - linearity plays an important role .
the next step in the construction is finding the zeros @xmath114 of @xmath197 .
one has ( z)=y_iy_j(x_ix_j - x^2_ij-(2)^-2 m_ij ) , where we used @xmath198 . substituting a parametrization for this null - cone in @xmath124 , e.g. @xmath199 , gives a 4th order polynomial .
however , for finding the 4 solutions we find it in this case more convenient to first diagonalize the matrix @xmath200 .
introducing @xmath201 , the equation for the zeros reduces to @xmath202 , which in addition to the null - cone condition , @xmath203 , is now easily solved by y^(a)=(,(-1)^a+1 , i),y^(a+2)=(y^(a))^ * ,   a=1,2,[eq : ys ] where we fixed ( for generic @xmath82 ) the diagonalizing rotation @xmath204 by ordering @xmath205 .
using we find for @xmath206 u^(a)=-((-1)^a,,0),u^(a+2)=-u^(a ) ,   a=1,2 .
one easily checks that @xmath207 and @xmath208 give rise to respectively the exponentially rising and falling solutions .
it is also instructive to give the explicit expressions for the matrices @xmath209 in ( choosing @xmath210 ) .
we note that for a @xmath65 matrix @xmath211 , and without loss in generality , @xmath212 , @xmath192 , @xmath185 and @xmath190 by resp . translations , ( gauge ) rotations , and suitable rescalings .
] we take @xmath213 , @xmath214 , @xmath215 , @xmath216 and @xmath217 in , such that m^+(z)=- , + & & + m^-(z)=- .
we checked that , @xmath218 , evaluated using is independent of @xmath23 and agrees with the result derived in ref .
@xcite .
we next address solving @xmath219 , which for charge 2 can be written as , since @xmath220 .
for the same choice of parameters in as above , @xmath213 , @xmath221 , @xmath215 , @xmath216 and @xmath217 , this gives ( with @xmath222 and @xmath210 ) = 2(xy).[eq : mup ] although the dependence on @xmath82 is complicated , the integral over @xmath23 turns out to be manageable ( as was observed before in the context of charge 2 monopoles , although we here choose not to express the solution in terms of theta functions  @xcite ) . to solve the equation we first rewrite the right - hand side of using the fact that @xmath223 ( cmp . ) , & = & + & = & + .in the last identity we used that @xmath198 , @xmath224 and the fact that @xmath155 implies @xmath225
. with @xmath226 we can now integrate , ( z)&=&2z(xy)+ ( ) + ii(z ) , + i(z)&&_(f_3 ^ -1(z),n)-_(1,n)+|z| , n , up to an irrelevant constant , where @xmath227 is the elliptic integral of the third kind . note that @xmath228 can alternatively be written as @xmath229 , from which it follows that @xmath230 , @xmath231 and @xmath232 .
] , _ ( s , n)_0^s .
we now combine these ingredients to give in terms of @xmath114 the exact form for the homogeneous solution of the green s function equation , v_a(z)=(z)((z))_a2=((z)-2zu x ) or putting in all the relevant expressions v_a(z)&=&(i ) + & & ( 4)^- ( -y_1f_1(z)-(-1)^aiy_2f_2(z))^()^. substituting @xmath233 , with @xmath234 as defined in gives @xmath235 and @xmath236 , and thereby the green s function , once we specify the parameters involved in the solutions to the nahm equation .
the discontinuities in @xmath166 at @xmath237 and @xmath238 implied by the nahm equation , , impose constraints which are ( like the quadratic adhm constraint ) in general difficult to solve .
work is in progress to describe the full parameter space for su(2 ) and charge 2 , but in ref .
@xcite we did find two non - trivial parametrizations for which we illustrate here in a number of figures how the full caloron solutions look like , using the exact green s function as constructed in the previous section . taking advantage of the possibility to work with arbitrary arithmetic precision
the programme maple has been used for these calculations .
the configurations are formulated in terms of the two intervals @xmath239 $ ] and @xmath240 $ ] , each associated with two constituent monopoles of equal magnetic charge , but opposite in sign from one interval to the next .
apart from a shift and ( gauge ) orientation , the configuration is described by a shape ( @xmath191 ) and scale ( @xmath190 ) parameter ( see sect .
[ sec : ch2 ] ) , for simplicity assumed to be the same on both intervals .
we also take all constituents to be of equal mass , @xmath241 ( @xmath242 ) , most relevant for the confined phase with @xmath243 .
the two periodic and two anti - periodic chiral fermion zero - modes each have support on oppositely charged constituents ( see fig .  [
fig : fig1 ] ) and in the far field limit it was found that the zero - mode density ( summed over the two zero - modes implied by the index theorem ) is described by a disc singularity , bounded by an ellipse with semi - major axis @xmath244 and eccentricity @xmath245 ( where @xmath246 ) .
this revealed that the core of a cluster of like - charged constituents is in general extended , unless the individual constituents are well separated .
the far field only describes the ( algebraic ) abelian component of the gauge field and to resolve the structure of the core we need to determine the full non - abelian structure .    for the first parametrization ( called  rectangular " ) the two discs are parallel and separated in height by a distance @xmath17 . the configuration is charaterized by the two extremal points along the major axis of the disc , also called _ would - be _ constituent locations .
] , y_m^(j)=(0,(-1)^md,(-1)^j(4)^-1d),2d = df_2(d ) , up to an overall shift and orientation ( the definition of @xmath247 , which involves @xmath191 , can be found in ) .
a typical example is shown in fig .
[ fig : conloc ] ( left ) . apart from @xmath191 and @xmath190 , the parameters that enter for the @xmath83-th interval are @xmath248 , @xmath215 , @xmath217 and @xmath249 . in all cases discussed here
we have also put @xmath250 .
one verifies that the discontinuities of @xmath166 at @xmath68 are given by @xmath251 with appropriately chosen @xmath252 ( cmp . ) , as discussed in detail in ref .
@xcite .    for
the second parametrization ( called  crossed " ) the two discs are coplanar and intersect , see fig .
[ fig : conloc ] ( right ) for a typical example .
their relative orientations can vary between perpendicular and coinciding ( for which @xmath191 is forced to 1 ) . here
we choose for the parameters in @xmath192 and @xmath185 to be non - trivial ( isospin ) rotations around the @xmath253-axis with angles @xmath254 , resp .
@xmath255 , and @xmath256 , whereas @xmath257 and @xmath258 are as in the rectangular case .
the would - be constituent locations are now given ( up to an overall shift and orientation ) by y_m^(j)=((-1)^j(4)^-1d,0,(-1)^m+j ( 4)^-1d-(-1)^md ) , where @xmath259 conveniently gives the orientation of each of the two discs with respect to the @xmath23-axis .
the angle @xmath260 originates from the definition of @xmath252 , which through determines the discontinuity of @xmath166 . to ensure the proper matching , the following three equations need to be satisfied  @xcite d()=8 d(1),df_2(d)+8d=0 . which determine @xmath261 , @xmath191 and @xmath190 for given @xmath260 , @xmath262 and @xmath17 .
[ fig : conloc ] illustrates that for these crossed configurations the two discs always overlap , unlike for the rectangular case .
this formed an important motivation for the present study , so as to determine in how far the overlapping discs would affect the behaviour in the core .
we will illustrate , using the exact formalism , how the action and zero - mode densities for these solutions behave .
for the rectangular case we are mostly interested in the monopole limit , @xmath263 . for finite @xmath17 ,
the density separates in two contributions , where each in the limit @xmath263 is exactly the density for a charge 2 monopole with the same values of @xmath191 and @xmath190 .
note that in the limit @xmath263 the matching conditions for @xmath166  decouple " the intervals , turning into pole conditions at @xmath264 , as is appropriate for multi - monopoles  @xcite .
we do recover from this the known results for the charge 2 monopoles  @xcite , like the doughnut structure for @xmath265 , corresponding to two superimposed monopoles .
the interest in the monopole limit comes from the fact that the zero - mode densities for multi - monopoles had not been studied in detail before  @xcite . in fig .
[ fig : monlim ] we give the densities for @xmath266 and @xmath267 , which is intermediate between the doughnut and well - separated monopole configurations . on the other hand , for @xmath17 small ( compared to @xmath40 ) the configuration will look like two non - dissociated calorons , and in particular is no longer static .
when @xmath190 remains much bigger than @xmath17 , forcing @xmath268 , these behave as two well - separated charge 1 instantons . otherwise , when @xmath190 is comparable to @xmath17 , one finds overlapping instantons  @xcite .
an example for the crossed configuration with @xmath269 and @xmath270 was already shown in fig .
[ fig : fig1 ] . in this case
@xmath190 is large enough for the like - charge constituents to be separated , as is particularly clear from the zero - modes , which are essentially no longer overlapping .
but two nearest neighbour ( oppositely ) charged constituents still show appreciable overlap . the distance between these nearest neighbours
is @xmath271 .
as this is comparable to @xmath40 we would expect the configuration to depend on time .
indeed , at the maxima of the action density its value of @xmath272 at @xmath273 is reduced by almost 50% at @xmath274 . at the center of mass , where the action density is much lower , there is still a time dependence .
however , far from all cores the field becomes static ) and topological charge @xmath14 the far field limit is static , _ provided _ all constituent monopoles have a non - vanishing mass ] .
increasing @xmath190 further ( which will push @xmath191 closer to 1 ) the configuration quickly turns into well separated spherically symmetric static bps monopoles .
more interesting is to consider the case with smaller @xmath190 , like @xmath275 and @xmath276 , for which the disc singularities were illustrated in fig .
[ fig : conloc ] .
the corresponding densities are shown in fig .
[ fig : close ] .
we see that the constituents are now so close that they form a doughnut , but we stress this is different from the _ static _ monopole doughnut , which has @xmath265 . since the oppositely charged constituents now are as close as 0.438 , which is considerably smaller than the time extent , the solution will have a strong time dependence .
when @xmath190 is decreased even further , it will turn into a charge 2 instanton localized in both space and time .
in this section we make some comments on the far field limit of the gauge field . as we discussed before
, the gauge field far removed from any core becomes abelian ( as well as static ) .
the abelian subgroup is the one that leaves the holonomy invariant , in the periodic gauge equivalent to leaving the constant asymptotic value of the adjoint higgs field @xmath277 invariant . for definiteness ,
let us consider the case of su(2 ) with @xmath278 , @xmath40 and @xmath279 ( i.e. @xmath280 ) .
up to exponential corrections we have a^_0(x)=2i - i(x ) , where we can express @xmath281 in terms of the far field limit of the green s function at the impurities  @xcite ( x)=^-1 ^ -1_i(f^_x(_2,_2 ) _ 2^i ) .    using
the twistor description of magnetic monopoles hurtubise  @xcite was able to explicitly compute the asymptotic higgs field for the su(2 ) magnetic monopole long ago .
the function he found for this algebraic tail amazingly agrees exactly with @xmath96 , , which was introduced to describe the caloron zero - mode density ( @xmath83 denotes the interval and hence the type of constituents to which the corresponding zero - modes would localize  @xcite ) .
as mentioned before , from our multi - caloron results one can recover the multi - monopole results by sending the constituent monopoles with the  unwanted " magnetic charge to infinity , cmp . .
although this tends to be cumbersome to show , @xmath281 in the far field can be written as @xmath282 , where @xmath283 is the contribution coming from the type @xmath83 constituent monopoles and the difference in sign is due to the sign change in the magnetic charge .
this is simply because the field is abelian in the far field and linear superposition preserves the self - duality .
hence , @xmath284 , such that for the su(2 ) caloron ( x)=2_1(x)-2_2(x ) .
we checked that this result indeed holds for the solutions discussed in sect .
[ sec : plots ] , even when the two types of monopole structures are not well - separated .
this relation trivially holds for the axially symmetric solutions that were introduced in ref .
@xcite , where @xmath281 was explicitly shown to factorize in a sum of point charge contributions , compatible with what was found in ref .
@xcite for the zero - mode densities .
therefore , in the far field limit ( i.e. for the algebraic part ) the singularity structure in the zero - mode density agrees _ exactly _ with the abelian charge distribution , as given by @xmath285 .
such a relation is at the heart of using chiral fermion zero - modes as a filter to isolate the underlying topological lumps from rough lattice monte carlo configurations  @xcite .
in this paper we have analyzed the higher charge caloron solutions and showed how to obtain exact results by suitably combining techniques developed in the context of the nahm transformation and the adhm formalism .
the aim of these studies has been to establish that su(@xmath0 ) caloron solutions of charge @xmath14 can be described in terms of @xmath15 monopole constituents , and that these can be viewed as independent constituents . a natural way to get
an ensemble would be to consider approximate superpositions of @xmath14 charge 1 calorons , but this would lead to an unwanted memory effect , with constituents remembering from which caloron they originated  @xcite . our studies , within the context of self - dual configurations , have shown nevertheless that the constituents have an independent identity , with the only requirement that the net magnetic and electric charge of the configuration vanishes ( each of the @xmath0 types of constituents should occur with the same number ) .
a recent lattice study  @xcite , using the technique of over - improvement  @xcite , fully confirms this picture .
it is therefore reasonable to consider the constituents as the independent building blocks for constructing an ensemble of monopole constituents , something that was not questioned in ref .
@xcite , but like for the instanton liquid  @xcite forms an essential assumption in a semi - classical study .
clearly the expectation is that semi - classical methods no longer work in the confined regime , at least for the part of the parameter space that corresponds to well - separated constituents , that is typically associated to instantons with a large scale parameter .
it is not unlikely that the density of these constituents at low temperatures is so high that they form a coherent background and as such will no longer easily be recognized as lumps . with high quark densities leading to deconfinement
, it may perhaps be that a high constituent monopole density will lead to confinement  @xcite , although for now we have to leave this as a speculation .
instantons that overlap get deformed and depending on the relative gauge orientation tend to  repel " , i.e. inspecting the action density distribution they do not get closer than a certain distance  @xcite . when deconstructing instantons in monopole constituents , interestingly
only like - charge constituents will show this effect , manifesting itself through the extended core structure . for unlike charges , from the point of view of the abelian field , the configuration behaves as with linear superposition .
if as a consequence of this all abelian charge is annihilated , it disappears through forming a small instanton ( localized in space and time ) , which in the limit of zero size describes the boundary of the moduli - space .
the interaction between constituents of opposite duality is more complicated  @xcite .    in conclusion ,
calorons with non - trivial holonomy have revealed a rich structure , incorporating traditional instanton physics , but allowing for gauge fields that inherit some essential features associated to a confining background not present in the traditional formulations .
the fact that the underlying constituents are monopoles opens the way to describe the confining aspects of the theory in terms of these degrees of freedom .
much work remains to be done when it comes to understanding the dynamics , but we hope to have convinced the reader that a consistent picture is developing that holds considerable promise for the future .
we thank conor houghton for pointing us to hurtubise s paper , and chris ford for discussions on the analytical aspects .
we are grateful to christof gattringer , michael ilgenfritz , boris martemyanov and michael mller - preussker for stimulating discussions and correspondence on lattice implementations .
the research of fb is supported by fom .
he is grateful to chris ford and conor houghton for their hospitality while visiting trinity college in dublin , and thanks them , as well as werner nahm and samson shatashvili , for their interest .
in this appendix we derive the zero - mode limit for the action density , which assumes that the distance of @xmath82 and the constituents of type @xmath83 to all constituents of type @xmath86 is large , but where @xmath82 and the constituent locations of type @xmath83 may otherwise be arbitrary .    as in ref .
@xcite we take @xmath286 for computing @xmath66 and use that we can write @xmath287 , where @xmath288 and @xmath289 , with & & _ m^-1t_m , + & & + & & f_m .
we note the @xmath290 has _ no _ remaining dependence on the constituent locations of type @xmath83 . writing @xmath291 , with k , k , l , l .
we find @xmath292 .
we next use k^-1=,kk^-1= and note that in the zero - mode limit @xmath293 , @xmath294 , @xmath295 and @xmath296 are exponentially small ( cmp . ref .
@xcite , app .
a ) , such that ( ie^-ix_0(_2k - lk))=(e^-2ix_0k_++ ) ( l_++ ) . with the definition of @xmath297 we now find l_++=r_m+1 ^ -1(_m+1)(r_m+1 ^ -(_m+1)+s_m+1 ) u_m(r_m-1^+(_m)+s_m ) , where -8mm&&u_m=_m+1^+f_m^+(_m+1)f_m^+(_m)^-1r_m^-1 ( _ m)_m^+-_m+1 ^ -f_m^-(_m+1)f_m^-(_m)^-1r_m^-1 ( _ m)_m^- , + -8mm&&^_m=_k(r_m^- ( _ m)+s_m)^-1r_m-1^(_m ) ,  ^_m=_kr_m^(_m)(r_m-1^+(_m)+s_m)^-1.[eq : szmlim ] and @xmath298 contains _ all _ contributions due to the constituent locations of type @xmath83 , up to _ exponential _ corrections in the distance of these , _ and _ of @xmath82 , to the other constituents . hence @xmath299 splits into the sum of two contributions , @xmath300 and @xmath301 $ ]
, where the last term only depends on the constituent locations of type @xmath86 whose contribution will decay inversely proportional to the fourth power of their distance . allowing for _ algebraic _ decay ( or in the monopole limit , sending all constituents of type @xmath86 to infinity ) such that in addition @xmath302 , one thus finds .    a simple way to derive the result for the far field limit in is by noting that in this case _ all _ @xmath303 are exponentially small and @xmath304 can be approximated by @xmath305 , with @xmath306 .
this therefore acts as a projection on the @xmath307 component and is thus seen to lead to @xmath308 . using the fact that  @xcite @xmath309 gives the required result .
harrington and h.k .
shepard , phys .
d17 ( 1978 ) 2122 ; ibid .
d18 ( 1978 ) 2990 .
gross , r.d .
pisarski and l.g .
yaffe , rev .
53 ( 1981 ) 43 .
kraan and p. van baal , phys .
lett . b428 ( 1998 ) 268 [ hep - th/9802049 ] ; nucl .
phys . b533
( 1998 ) 627 [ hep - th/9805168 ] ; phys .
lett . b435
( 1998 ) 389 [ hep - th/9806034 ] .
k. lee and p.
yi , phys .
d56 ( 1997 ) 3711 [ hep - th/9702107 ] ; k. lee , phys . lett . b426
( 1998 ) 323 [ hep - th/9802012 ] ; k. lee and c. lu , phys .
d58 ( 1998 ) 025011 [ hep - th/9802108 ] . n. weiss , phys . rev .
d24 ( 1981 ) 475 .
davies , t.j .
hollowood , v.v .
khoze and m.p .
mattis , nucl .
phys b559 ( 1999 ) 123 [ hep - th/9905015 ] .
d. diakonov , n. gromov , v. petrov and s. slizovskiy , `` quantum weights of dyons and of instantons with non - trivial holonomy , '' hep - th/0404042 .
ilgenfritz , m. mller - preussker , and a.i .
veselov , in : _ lattice fermions and structure of the vacuum _ , eds . v. mitrjushkin and g. schierholz ( kluwer , dordrecht , 2000 ) , 345 [ hep - lat/0003025 ] ; e .- m .
ilgenfritz , b.v .
martemyanov , m. mller - preussker , s. shcheredin and a.i .
veselov , phys . rev .
d66 ( 2002 ) 074503 [ hep - lat/0206004 ] . c. gattringer , phys
d67 ( 2003 ) 034507 [ hep - lat/0210001 ] ; c. gattringer and s. schaefer , nucl .
b654 ( 2003 ) 30 ; [ hep - lat/0212029 ] ; c. gattringer and r. pullirsch , phys .
d69 ( 2004 ) 094510 [ hep - lat/0402008 ] .
f. bruckmann , d. ngrdi and p. van baal , nucl .
b666 ( 2003 ) 197 [ hep - th/0305063 ] .
m. garca prez , a. gonzlez - arroyo , a. montero and p. van baal , jour . of high energy phys .
06 ( 1999 ) 001 [ hep - lat/9903022 ] .
p. van baal , in : _ lattice fermions and structure of the vacuum _ , eds . v. mitrjushkin and g. schierholz ( kluwer , dordrecht , 2000 ) , p. 269
[ hep - th/9912035 ] .
ilgenfritz , b.v .
martemyanov , m. mller - preussker and a.i .
veselov , phys .
d69 ( 2004 ) 114505 [ hep - lat/0402010 ] .
j. hurtubise , commun .
97 ( 1985 ) 381 .
f. bruckmann and p. van baal , nucl .
b645 ( 2002 ) 105 [ hep - th/0209010 ] . w. nahm ,  self - dual monopoles and calorons " , in : lect .
notes in physics .
201 , eds .
g. denardo , e.a .
( 1984 ) p. 189 .
h. panagopoulos , phys .
d28 ( 1983 ) 380 .
atiyah , n.j .
hitchin , v.g .
drinfeld , yu . i. manin , phys .
65 a ( 1978 ) 185 ; m.f .
atiyah , _ geometry of yang - mills fields _ , fermi lectures , ( scuola normale superiore , pisa , 1979 ) .
m. garca prez , a. gonzlez - arroyo , c. pena and p. van baal , phys .
d60 ( 1999 ) 031901 [ hep - th/9905016 ] .
chernodub , t.c .
kraan and p. van baal , nucl .
phys . b(proc.suppl . )
83 - 84 ( 2000 ) 556 [ hep - lat/9907001 ]
. w. nahm ,  the algebraic geometry of multi - monopoles " , in : _ group theoretical methods in physics _
, eds . m. serdaroglu , e.a . , istanbul 1982 , p.456 . w. nahm ,  multi - monopoles in the adhm construction " , in : _ gauge theories and lepton hadron interactions _
, eds . z. horvath , e.a .
, ( crip , budapest , 1982 ) .
brown , h. panagopoulos and m.k .
prasad , phys .
d26 ( 1982 ) 854 ; + a.s .
dancer , comm . math . phys .
158 ( 1993 ) 545 .
p. forgcs , z.  horvth and l. palla , nucl .
b192 ( 1981 ) 141 ; m.f .
atiyah , and n.j .
hitchin , _ the geometry and dynamics of magnetic monopoles _ , ( princeton univ .
press , 1988 ) .
corrigan and p. goddard , annals phys .
( n.y . ) 154 ( 1984 ) 253 . m. garca prez , t.g .
kovcs and p. van baal , phys .
b472 ( 2000 ) 295 [ hep - ph/9911485 ] .
f. bruckmann , m. ilgenfritz , b. martemyanov and p. van baal , hep - lat/0408004 .
m. garca prez , a. gonzlez - arroyo , j. snippe and p. van baal , nucl .
phys . b413 ( 1994 ) 535 [ hep - lat/9309009 ] .
t. schfer and e.v .
shuryak , rev .
( 1998 ) 323 [ hep - ph/9610451 ] ; d. diakonov , prog . part .
51 ( 2003 ) 173 [ hep - ph/0212026 ] . f. bruckmann , d. ngrdi and p. van baal , acta phys . polon .
b34 ( 2003 ) 5717 [ arxiv : hep - th/0309008 ] .
ilgenfritz , b.v .
martemyanov , m. mller - preussker and a.i .
veselov , eur .
j. c34 ( 2004 ) 439 [ hep - lat/0310030 ] ; e.m .
ilgenfritz , m. mller - preussker , b.v .
martemyanov and p. van baal , phys .
d69 ( 2004 ) 097901 [ hep - lat/0402020 ] .


the quantum states of particle and anti - particle excitations in a superconductor ( bogoliubov quasiparticles ) are related by a unitary transformation , which means that they can be represented by a _ real _ wave function . in this so - called majorana representation
the @xmath4 scattering matrix @xmath1 at the fermi level is real orthogonal rather than complex unitary @xcite .
since the orthogonal group @xmath5 is doubly connected , this immediately implies a twofold distinction of scattering problems in a superconductor : the subgroup @xmath6 of scattering matrices with determinant @xmath7 , connected to the unit matrix , is called _ topologically trivial _ , while the disconnected set @xmath8 of scattering matrices with determinant @xmath9 is called _ topologically nontrivial_. in mathematical terms , the experimental search for majorana bound states can be called a search for systems that have @xmath10 .
this search has been reviewed , from different perspectives , in refs .  .
if the scattering is chaotic the scattering matrix becomes very sensitive to microscopic details , and it is useful to develop a statistical description : rather than studying a particular @xmath1 , one studies the probability distribution @xmath11 in an ensemble of chaotic scatterers .
this is the framework of random - matrix theory ( rmt ) @xcite .
the ensemble generated by drawing @xmath1 uniformly from the unitary group @xmath12 , introduced by dyson in the context of nuclear scattering @xcite , is called the circular unitary ensemble ( cue ) .
superconductors need a new ensemble .
a natural name would have been the circular orthogonal ensemble ( coe ) , but since that name is already taken for the coset @xmath13 , the alternative name circular real ensemble ( cre ) is used when @xmath1 is drawn uniformly from @xmath5 .
the rmt of the cre , and the physical applications to majorana fermions and topological superconductors , have been reviewed recently @xcite .
the uniformity of the distribution requires ideal coupling of the scattering channels to the continuum , which physically means that the discrete spectrum of a quantum dot is coupled to metal electrodes by ballistic point contacts . if the point contact contains a tunnel barrier , then @xmath11 is no longer uniform but biased towards the reflection matrix @xmath14 of the barrier .
the modified distribution @xmath15 is known @xcite , it goes by the name `` poisson kernel '' and equals @xmath16 in the cre @xcite .    ) with superconducting boundaries , connected to metal electrodes by a pair of point contacts ( supporting a total of @xmath17 propagating modes ) .
a magnetic insulator introduces a tunnel barrier in each point contact ( transmission probability @xmath18 per mode ) .
a magnetic vortex may introduce a majorana bound state in the quantum dot . ]    in the present work we apply this result to the scattering ( andreev reflection ) in a superconducting quantum dot ( andreev billiard ) , see fig .
[ fig_layout ] .
we focus in particular on the effect of a bound state at the fermi level ( @xmath19 ) in the quantum dot , a so - called majorana zero - mode or majorana bound state .
in addition to the scattering matrix , which determines the thermal and electrical conductance , we consider also the time - delay matrix @xmath2 .
the eigenvalues of @xmath3 are positive numbers with the dimension of time , that govern the low - frequency dynamics of the system ( admittance and charge relaxation @xcite ) .
moreover , the trace of @xmath3 gives the density of states and @xmath3 and @xmath1 together determine the thermopower @xcite .
the joint distribution of @xmath1 and @xmath3 is known for ballistic coupling @xcite , here we generalize that to tunnel coupling .
the effect of a tunnel barrier on the time - delay matrix has been studied for complex scattering matrices @xcite , but not yet for real matrices .
one essential distinction is that the tunnel barrier has no effect on the density of states in the cue and coe , but it does in the cre .
the outline of the paper is as follows .
the next two sections formulate the scattering theory of the andreev billiard and the appropriate random - matrix theory .
our key technical result , the joint distribution @xmath0 , is given in sec .
[ sect_ps0q ] .
we apply this to the simplest single - channel case ( @xmath20 ) in sec .
[ taudos ] , where obtain a remarkable scaling relation : for a high tunnel barrier ( transmission probability @xmath21 ) the distribution @xmath22 of the density of states at the fermi level is described by a one - parameter scaling function @xmath23 : @xmath24    the average density of states in the multi - channel case is calculated in sec .
[ sec_dos ] . by relating the ensemble averages of @xmath3 and @xmath1
we derive the relation @xmath25\right),\label{rhoaveragegeneral}\ ] ] for a mode - independent tunnel probability @xmath18 . in the cue and
coe the average scattering matrix @xmath26 is just equal to @xmath14 , so @xmath27 remains equal to its ballistic value @xmath28 , but the cre is not so constrained .
applications to the thermal conductance @xmath29 and the electrical ( andreev ) conductance @xmath30 follow in secs .
[ transmissiong ] and [ electrcond ] . for ballistic coupling
it is known that @xmath31 is the same with or without the majorana bound state @xcite .
( this also holds for @xmath32 @xcite . ) in the presence of a tunnel barrier this is no longer the case , but we find that the majorana bound state remains hidden if even a single scattering channel has @xmath33 .
the distribution of @xmath30 , in contrast , is sensitive to the presence or absence of the majorana bound state even for ballistic coupling @xcite .
the way in which @xmath34 changes as we tune the system through a topological phase transition , at which a majorana bound state emerges , is calculated in sec .
[ phasetransition ] .
we conclude in sec .
[ conclude ] .    in the main text we focus on the results and applications .
details of the calculations are moved to the appendices .
these also contain more general results for other rmt ensembles , with or without time - reversal and/or spin - rotation symmetry .
( both symmetries are broken in the cre . )
fig .  [ fig_layout ] shows the scattering geometry , consisting of a superconducting quantum dot ( andreev billiard ) on the surface of a topological insulator , connected to normal metal electrodes by point contacts .
the hamiltonian @xmath35 of the quantum dot is related to the energy - dependent scattering matrix @xmath36 by the mahaux - weidenmller formula @xcite , @xmath37 the @xmath38 matrix @xmath39 couples the @xmath40 energy levels in the quantum dot ( mean level spacing @xmath41 ) to a total of @xmath42 propagating modes in the point contact .
we assume that degeneracies are broken by spin - orbit coupling in the topological insulator in combination with a magnetic field ( perpendicular to the surface ) .
all degrees of freedom are therefore counted separately in @xmath43 and @xmath40 , as well as in @xmath41 .
the electron - hole degree of freedom is also included in the count , but we leave open the possibility of an unpaired majorana fermion  a coherent superposition of electron and hole quasiparticles that does not come with a distinct antiparticle .
an odd level number @xmath40 indicates the presence of a majorana bound state in the quantum dot , produced when a magnetic vortex enters @xcite .
an odd mode number @xmath43 signals a propagating majorana mode in the point contact , allowed by a @xmath44-phase difference between the superconducting boundaries @xcite .
the @xmath43 modes have an energy - independent transmission probability @xmath45 $ ] per mode .
if we choose a basis such that the coupling matrix @xmath39 has only nonzero elements on the diagonal , it has the explicit form @xcite @xmath46    notice that the tunnel probability @xmath47 determines the reflection amplitude @xmath48 $ ] up to a sign . the conventional choice is to take @xmath49 , when @xmath50 can be written as @xmath51 alternatively , if @xmath52 one has @xmath53 given by @xmath54 the two choices are equivalent for ballistic coupling , @xmath55 , but for a high tunnel barrier @xmath56 one has @xmath57 while @xmath58 .
the sign change of @xmath59 is a topological phase transition @xcite , which we will analyze in section [ phasetransition ] .
for now we take @xmath49 for all @xmath60 , so @xmath50 .    from the scattering matrix
we can obtain transport properties , such as the electrical and thermal conductance , and thermodynamic properties , such as the density of states .
if we restrict ourselves to properties at the fermi level , @xmath19 , we need the matrix @xmath61 and the derivative @xmath62 the unitarity of @xmath36 implies that @xmath3 is hermitian , so it has real eigenvalues @xmath63 with the dimension of time .
the @xmath63 s are called ( proper ) delay times and @xmath3 is called the wigner - smith time - delay matrix @xcite .
the fermi - level density of states @xmath64 is obtained from @xmath3 via the birman - krein formula @xcite , @xmath65    for the thermal conductance we partition the modes into two sets , @xmath17 , each set connected to a different terminal , and decompose the scattering matrix into reflection and transmission subblocks , @xmath66 a small temperature difference @xmath67 between the two terminals , at average temperature @xmath68 , drives a heat current @xmath69 .
the thermal conductance @xmath70 in the low - temperature linear - response limit @xmath71 is given by @xmath72 the quantum @xmath73 is a factor - of - two smaller than in systems without superconductivity @xcite , due to our separate counting of electron and hole degrees of freedom that allows to account for the possibility of propagation via an unpaired majorana mode .
if we keep the two terminals at the same temperature but instead apply a voltage difference , we can drive an electrical current .
we consider a situation where both terminal 2 and the superconductor are grounded , while terminal 1 is biased at voltage @xmath74 .
the current @xmath75 from terminal 1 to ground is then given by the andreev conductance @xmath76 in the zero - temperature , zero - voltage limit .
in the last equality we used the particle - hole symmetry relation @xmath77 at @xmath19 , where the @xmath78 pauli matrices act on the electron ( @xmath79 ) and hole ( @xmath80 ) degree of freedom .
for a statistical description we consider an ensemble of quantum dots , each with its own random hamiltonian @xmath35 .
the mean level spacing @xmath41 and coupling matrix @xmath39 are kept fixed . if the wave dynamics in the quantum dot is chaotic , the ensemble is fully characterized by the presence or absence of certain fundamental symmetries .
this is the universal framework of random - matrix theory .
superconducting systems are characterized by particle - hole symmetry , @xmath81 the pauli matrices @xmath82 can be removed from the symmetry relation by a unitary transformation @xmath83 after which we simply have @xmath84 in this so - called majorana basis the hamiltonian is real antisymmetric , @xmath85 with @xmath86 .
if no other symmetries are imposed on the hamiltonian we have the class - d ensemble of random - matrix theory @xcite , with gaussian probability distribution @xmath87 the eigenvalues of the antisymmetric @xmath88 matrix @xmath35 come in @xmath89 pairs , hence if @xmath40 is odd there must be a nondegenerate eigenvalue @xmath19 at the fermi level , in the middle of the superconducting gap .
this so - called majorana bound state is the hallmark of a topologically nontrivial superconductor @xcite .
if @xmath40 is even there is no level pinned to @xmath19 , and the superconductor is called topologically trivial .
it is helpful to encode the distinction in a topological quantum number @xmath90 that counts the number of majorana bound states , so @xmath90 equals 0 or 1 if the superconductor is topologically trivial or non - trivial , respectively .    in the scattering matrix
the presence of a majorana bound state is signaled by the sign of the determinant , @xmath91 this can be seen directly from the definition in the majorana basis : for @xmath40 even the matrix @xmath85 is invertible , so we have @xmath92 since @xmath93 if @xmath94 .
for @xmath40 odd the bound state contributes to the determinant a factor @xmath95 for some vector @xmath96 , so @xmath10 .    the class - d ensemble of scattering matrices
thus consists of two disjunct sets : the special orthogonal group @xmath97 of orthogonal matrices with determinant @xmath7 in the topologically trivial case , and the complement @xmath8 of orthogonal matrices with determinant @xmath9 in the topologically nontrivial case .
for ballistic coupling ( @xmath98 for all @xmath60 ) the matrices @xmath1 and @xmath3 are statistically independent @xcite , so they can be considered separately .
the class - d ballistic scattering matrix is uniformly distributed in @xmath99  uniformity being defined with respect to the haar measure @xcite .
this is the circular real ensemble ( cre ) , the analogue for real orthogonal matrices of the circular unitary ensemble ( cue ) for complex unitary matrices @xcite .
the class - d ballistic time - delay matrix has probability distribution @xcite , @xmath100 where @xmath101 is the heisenberg time and @xmath102 restricts @xmath3 to positive definite real symmetric matrices .
this constraint can be implemented more directly by defining @xmath103 with @xmath104 the set of @xmath105 matrices with real elements .
the distribution is then equivalent to @xmath106 multiply by the jacobians @xmath107 , for @xmath108 and @xmath109 . ] to the wishart distribution @xcite @xmath110 remarkably , there is no dependence on the topological quantum number for ballistic coupling : @xmath3 has the same distribution irrespective of the presence or absence of a majorana bound state .
tunnel coupling is described by a reflection matrix @xmath14 ( from outside to outside ) and transmission matrix @xmath111 ( from outside to inside ) . in the majorana basis these are real matrices , parameterized by @xmath112 .
\end{split } \label{rbdef}\ ] ] as we derive in app .
[ app_ps0q ] , the matrix product @xmath113 determines the joint distribution    [ ps0q ] @xmath114    as a check , we can integrate out the time - delay matrix to obtain the marginal distribution of the scattering matrix , @xmath115 the jacobian evaluates to @xcite @xmath116 and we recover the class - d poisson kernel to indicate an unspecified normalization constant .
the poisson kernel is normalized , @xmath117 with @xmath118 the haar measure on @xmath119 , so we use @xmath120 instead of @xmath121 . ]
@xcite , @xmath122    the joint distribution tells us that @xmath1 and @xmath3 become correlated in the presence of a tunnel barrier .
however , @xmath1 remains independent of the matrix product @xmath123 so that the joint distribution of @xmath1 and @xmath124 factorizes , @xmath125 the transformation from @xmath3 to @xmath124 removes the effect of the tunnel barrier on the time - delay matrix ( see app .
[ app_ps0q ] ) .
we calculate the joint distribution @xmath0 of scattering matrix and time - delay matrix in the presence of a tunnel barrier , starting from the known distribution @xmath241 without a barrier @xcite .
the application in the main text concerns symmetry class d , but for the sake of generality and for later reference we give results for all four altland - zirnbauer @xcite symmetry classes d , diii , c , ci , as well as for the three wigner - dyson @xcite symmetry classes a , ai , aii .
the symmetry indices that distinguish the ensembles are listed in table [ table_az ] , see ref .   for an overview of this classification .
the unitary matrix @xmath1 and the hermitian matrix @xmath3 are real in class d , complex in class a , and quaternion in class c. we consider these three symmetry classes without time - reversal symmetry first , and then include the constraints of time - reversal symmetry in classes diii , ci , ai , aii
. do not exhaust the tenfold way classification of random - matrix theory : there are three more chiral classes @xcite ( labeled aiii , bdi , cii ) that require separate consideration @xcite . ]    .parameters that appear in the distribution of the scattering matrix and time - delay matrix , for each of the altland - zirnbauer and wigner - dyson symmetry classes .
notice that a different set of indices @xmath242 , @xmath243 govern the energy level statistics @xcite .
the degeneracy factor @xmath244 refers to the kramers degeneracy of the scattering channels and the delay times , ignoring uncoupled spin bands .
( the energy levels may have a different degeneracy . ) [ cols="^,^,^,^,^,^,^,^,^,^,^ " , ]      without the barrier @xmath245 is independent of @xmath124 and uniformly distributed , @xmath246 the differential @xmath247 indicates the haar measure for the unitary matrix @xmath245 and the euclidean measure for the hermitian matrix @xmath124 .
the ballistic time - delay matrix distribution is given by @xcite    [ pq0result ] @xmath248    degenerate eigenvalues of @xmath124 are counted only once in @xmath249 and @xmath250 . in terms of the degeneracy factor @xmath244 from table [ table_az ]
this can be written as @xmath251 the channel number @xmath43 also does not include degeneracies , so the total number of eigenvalues of @xmath124 is @xmath252 .
the characteristic time @xmath253 differs from the heisenberg time @xmath254 by a numerical coefficient , , with @xmath41 the mean spacing of _ nondegenerate _ levels .
the ratio @xmath255 then equals the degeneracy of energy levels divided by the degeneracy of delay times @xcite .
it is unity in all symmetry classes except c and ci , where the delay times have a kramers degeneracy that the energy levels lack @xcite . ]
see table [ table_az ] .
insertion of the barrier , with unitary scattering matrix @xmath256 transforms @xmath245 into @xmath257 variations of @xmath1 and @xmath245 are related by @xcite @xmath258    the differentials @xmath259 are anti - hermitian matrices , @xmath260 .
the number of degrees of freedom of the off - diagonal elements are given by @xmath261 and the number of degrees of freedom of the diagonal elements by @xmath262 .
as summarized in table [ table_alphabeta ] , real matrices ( class d ) have @xmath263 , @xmath264 , complex matrices ( class a ) have @xmath265 , @xmath266 , and quaternion matrices ( class c ) have @xmath267 , @xmath268 .
these parameters determine the jacobian @xcite @xmath269    eq .
also implies the relation between the time - delay matrices , @xmath270 the off - diagonal elements of the hermitian matrix @xmath3 have @xmath261 degrees of freedom , the diagonal elements have one single degree of freedom in each symmetry class . the jacobian is then given by @xmath271    the joint probability distribution @xmath0 now follows upon division of @xmath241 by the product of jacobians , @xmath272 substituting @xmath273 into eq .
we thus arrive at the joint distribution @xmath274 the class - d result from the main text follows for @xmath263 , @xmath264 , @xmath275 , @xmath276 , @xmath277 .
time - reversal symmetry equates the scattering matrix to its transpose @xmath278 in class ai and ci and to its dual @xmath279 in class aii and diii .
( the dual of a matrix is @xmath280 . )
we use a unified notation @xmath281 , where the tilde indicates the transpose or the dual , whichever is appropriate for that symmetry class . the symmetry @xmath281 allows for the `` square root '' factorization @xmath282 with unitary @xmath283 .
the time - delay matrix in these symmetry classes is constructed such that it satisfies the same symmetry , @xmath284 this redefinition of @xmath3 differs from eq .   by a unitary transformation , so the delay times are not affected .
the ballistic @xmath124 and @xmath245 are again independent @xcite , distributed according to eqs .   and with the appropriate values of @xmath285 and @xmath261 from table [ table_az ] .
these numbers now count the diagonal and off - diagonal degrees of freedom of the symmetrized differential @xmath286 constrained by @xmath287 , @xmath288 .
the matrix elements of @xmath289 are imaginary in class ai ( @xmath265 , @xmath264 ) , @xmath290 times a quaternion , with four real coefficiients @xmath291 .
the matrix @xmath292 is the @xmath198 unit matrix . ] in class aii ( @xmath265 , @xmath268 ) , and of the form @xmath293 with @xmath294 imaginary in class ci ( @xmath295 , @xmath266 ) and @xmath294 real in class diii ( @xmath263 , @xmath266 ) .
the elements of the hermitian matrix @xmath3 are real in class ai , quaternion in class aii , and of the form @xmath296 with @xmath294 real in both classes ci and diii .
the off - diagonal elements of @xmath3 have the same number of @xmath261 degrees of freedom as @xmath289 , but the diagonal elements have only a single degree of freedom irrespective of @xmath285 .
all of this is summarized in table [ table_alphabeta ] .
the symmetrization of the differential modifies the relation , which now reads @xmath297 the relation between @xmath3 and @xmath124 is similarly modified by the symmetrization , @xmath298 because the matrices @xmath283 , @xmath299 are unitary , the jacobians and are unchanged , @xmath300 we thus obtain the joint distribution @xmath301      the marginal distribution of the scattering matrix resulting from the jacobians and is @xmath302 including the normalization constant .
this formula combines the known expressions for the poisson kernel . in the altland - zirnbauer ensembles the average scattering matrix differs from @xmath14 ,
see ref .   and app .  [ app_saverage ] .
] in the wigner - dyson ensembles @xcite and in the altland - zirnbauer ensembles @xcite .
the present analysis confirms that eq
.   holds without modification in the two symmetry classes d and diii that support majorana zero - modes , depending on the sign of the determinant @xmath303 in class d and the sign of the pfaffian @xmath304 in class diii @xcite . as a check
, we can take @xmath20 , when @xmath305 in class d and @xmath306 in class diii .
the @xmath307 sign determines the presence or absence of a majorana bound state ( twofold degenerate in class diii ) .
since there is only a single element in the ensemble we should have @xmath308 , which is indeed what eq .
gives for @xmath20 , @xmath263 .
for ballistic coupling the distribution implies that the eigenvalues @xmath126 of @xmath127 have the @xmath90-independent distibution @xcite @xmath128 the unit step function @xmath129 ensures that @xmath130 for all @xmath131 .    in the single - channel case
@xmath20 we can use the joint distribution to immediately extend this result to arbitrary tunnel probability @xmath132 .
the scalar @xmath1 is pinned to @xmath133 , hence @xmath134 [ this definition of @xmath135 corresponds to @xmath136 from eq .  . ]    since @xmath135 then appears only as a scale factor , we conclude that the single eigenvalue @xmath137 of @xmath127 for @xmath20 and any @xmath138 $ ] has distribution @xmath139 the single - parameter scaling @xmath140 is tested numerically in fig .
[ fig_singlechannel ] , by drawing random hamiltonians from the gaussian class - d ensemble .
the excellent agreement serves as a check on our analytics .     for a single - channel chaotic scatterer , without ( @xmath141 ) or with ( @xmath142 ) a majorana bound state .
the histograms are numerical results obtained by generating random hamiltonians ( of size @xmath143 ) with distribution .
the scattering matrix , and hence the delay time , then follows from eqs .   and .
the different curves correspond to different transmission probability @xmath18 of the tunnel barrier .
rescaling with a factor @xmath144 , with @xmath135 defined in eq .  , makes all histograms collapse onto a single curve , in agreement with the analytical result . ]    with @xmath145 the distribution gives the scaling form of the density of states distribution @xmath22 from the introduction , in the tunneling regime @xmath21 when @xmath146 .
for ballistic coupling , integration of @xmath147 with distribution gives the average density of states at the fermi level @xcite , @xmath148 for @xmath149 .
the ensemble average diverges for @xmath150 .    to calculate the effect of a tunnel barrier we write @xmath151
, and then use the fact that @xmath124 is independent of @xmath1 and hence independent of @xmath152 .
the average of @xmath124 with distribution @xmath153 is proportional to the unit matrix , @xmath154 so the average density of states ( still for @xmath149 ) is given by @xmath155\right).\label{rhoaverage}\end{aligned}\ ] ] ( this is eq .   from the introduction . )
it remains to calculate the average of @xmath1 with the poisson kernel . in the wigner - dyson symmetry classes
this average is just @xmath14 , but as pointed out in ref .
this no longer holds in the altland - zirnbauer class d. a simple result for @xmath26 is possible for mode - independent tunnel probabilities , @xmath156 for all @xmath60 , see app .
[ app_saverage ] : @xmath157 where again the @xmath158 sign corresponds to @xmath141 ( without a majorana bound state ) and the @xmath159 sign to @xmath142 ( with a majorana bound state ) .
we thus arrive at the average density of states , @xmath160\right),\label{rhoaverageresult}\ ] ] plotted in fig .
[ fig_rhoaverage ] . in the ballistic limit @xmath161 the dependence on the majorana bound state drops out , while in the tunneling limit @xmath162 we obtain @xmath163 the @xmath164 divergence of the density of states for @xmath142 corresponds to the delta - function contribution from the majorana bound state in the closed system . without the majorana bound state ( @xmath141 )
the density of states at the fermi level remains finite in the @xmath162 limit , but it does remain above the normal - state value of @xmath165 .
this midgap spectral peak is characteristic for a class - d superconductor @xcite . while in a closed system
the peak is simply a factor of two , in the weakly coupled open system it is a larger factor @xmath166 , which only tends to 2 in the large-@xmath43 limit .
the fact that the ensemble of open systems does not reduce to an ensemble of closed systems in the limit @xmath162 is due to statistical fluctuations that remain important for small @xmath43 .    ) or in the presence ( @xmath142 ) of a majorana bound state .
the curves are calculated from the analytical expression . ]
we consider the thermal conductance in the simplest case @xmath167 of a quantum dot with single - mode point contacts .
these are majorana modes , carrying heat but no charge .
the scattering matrix @xmath168 is parameterized by @xmath169 the haar measure equals @xmath170 the same for both @xmath171 and @xmath172 .    for the tunnel barrier
we take the reflection matrix @xmath173 , with @xmath174 .
the poisson kernel then has the explicit form @xmath175    the dimensionless thermal conductance has distribution @xmath176 where as before , @xmath177 applies to @xmath141 and @xmath178 to @xmath142 . the distribution @xmath179 is the result @xcite for ballistic coupling ( @xmath98 , @xmath180 ) .    for identical tunnel barriers , @xmath181 ,
this reduces to @xmath182 quite remarkably , the distribution of the thermal conductance for two identical single - mode point contacts is unaffected by the presence of a tunnel barrier in the topologically nontrivial case .
[ fig_thermalcond ] is a numerical check of this analytical result .
, see eq .  , for a chaotic scatterer having two single - mode point contacts with identical tunnel probabilities @xmath18 .
the data points are numerical results for a random hamiltonian ( @xmath183 ) , the curves are the analytical result . in the presence of a majorana zero - mode ( @xmath142 )
the distribution is independent of @xmath18 . ]
notice that the distribution becomes independent of @xmath90 if @xmath184 .
this is a special case of a more general result , valid for any @xmath185 : @xmath186 in words : the probability distribution of the thermal conductance becomes independent of the presence or absence of a majorana bound state if the quantum dot is coupled ballistically to at least one of the scattering channels . in other words ,
_ ballistic coupling to a propagating majorana mode hides the majorana bound state . _
the proof is straightforward : if @xmath187 for one of the indices @xmath188 , then the poisson kernel is unchanged if we multiply @xmath189 , with @xmath190 .
[ the orthogonal matrix @xmath191 is defined in eq .  .
] the haar measure remains unchanged as well , and so does the thermal conductance . since @xmath192 , so @xmath171 is mapped onto @xmath172 , we conclude that @xmath193 .
this proof for the poisson kernel extends to the entire joint distribution of @xmath1 and @xmath3 : the transformation @xmath189 has no effect on the matrices @xmath3 and @xmath152 , so @xmath0 remains unchanged .
it follows that the probability distribution of the density of states is the same with or without a majorana bound state if @xmath98 for at least one of the scattering channels .
because a majorana mode is charge - neutral , no electrical current can be driven for @xmath167 .
a nonzero current @xmath75 is possible for @xmath194 , when terminal 1 biased at voltage @xmath74 has a distinct electron and hole mode .
we investigate the effect of a tunnel barrier for @xmath194 , @xmath195 . in the majorana basis
the expression for the andreev conductance reads @xmath196 with @xmath197 a @xmath198 real matrix .
the scattering matrix @xmath199 can be conveniently parameterized using three euler angles @xcite , @xmath200 where we have defined @xmath201 the haar measure on @xmath202 is given by @xcite @xmath203    because @xmath204 commutes with @xmath205 , the dimensionless conductance depends only on the euler angle @xmath206 , @xmath207 in point contact 1 we take a tunnel probability @xmath208 , the same for the electron and hole mode , while in point contact 2 we have tunnel probability @xmath209 for the unpaired majorana mode .
evaluation of the poisson kernel with @xmath210 gives the conductance distribution    @xmath211}{[(1\pm r_2(g_{\rm a}-1))^2-(g_{\rm a}-1\pm r_2)^2(1-\gamma_1)]^{3/2}},\;\;0<g_{\rm a}<2,\;\;r_2=\sqrt{1-\gamma_2},\label{ppmg1}\ ] ]    plotted in fig .
[ fig_electrcond ] for @xmath141 ( @xmath177 ) , @xmath142 ( @xmath178 ) and two values of @xmath212 .    ,
, for @xmath194 , @xmath195 , @xmath213 .
the data points are numerical results for a random hamiltonian ( @xmath183 ) , the curves are the analytical result . for ballistic coupling ( @xmath33 )
the distribution is the same with or without a majorana bound state . in the limit
@xmath162 the distribution becomes sharply peaked at @xmath214 . ]    in the limit @xmath215 , when terminal 2 is decoupled from the quantum dot , we recover the result @xcite @xmath216 independent of @xmath208 .
the conductance in this case is uniquely determined by the topological quantum number .    in the opposite limit
@xmath217 the distribution becomes independent of @xmath90 , @xmath218}{[1-(1-\gamma_1)(g_{\rm a}-1)^2]^{3/2}},\;\;{\rm if}\;\;r_2=0.\label{ppmgr20}\ ] ] this is a special case of a more general result , for any @xmath185 , @xmath219 where @xmath220 projects onto the modes coupled to terminal 2 .
ballistic coupling , even for a single mode , to terminal 2 therefore removes the dependence on the topological quantum number .
the proof of eq .
proceeds along the lines of the proof of eq .  , with the difference that the transformation @xmath189 should leave the upper - left block @xmath197 of @xmath1 unaffected  otherwise the andreev conductance would change .
this also explains why the @xmath90-dependence of @xmath34 remains for ballistic coupling to terminal 1 @xcite .
, while the other is coupled to electrode @xmath221 via a point contact . ]
the appearance of a majorana bound state is a topological phase transition .
there is a search for this transition in a nanowire geometry , see fig .
[ fig_nanowire ] , where it has been predicted to occur when the zeeman energy of a magnetic field ( parallel to the wire axis ) exceeds the gap induced by the proximity to a superconductor @xcite .
because the majorana bound states emerge pairwise at the two ends of the nanowire , the topological quantum number @xmath90 of the entire structure remains @xmath222 and the determinant @xmath223 of the full scattering matrix remains @xmath7 through the transition .
what changes is the sign of the determinant @xmath224 of the reflection submatrix . at the topological phase transition @xmath225 , implying a perfectly transmitted mode and a quantized peak in the thermal conductance @xcite .
we study the effect of the phase transition on the statistics of the electrical conductance , measured by contacting one end of the nanowire ( terminal 1 ) to a metal electrode at voltage @xmath74 , while the superconductor and the other end of the nanowire ( terminal 2 ) are at ground .
terminal 1 is connected to the nanowire via a point contact , thus creating a confined region ( quantum dot ) with chaotic scattering .
the minimal dimensionality of the scattering matrix @xmath1 of the quantum dot is @xmath226 : one electron and one hole mode connected to terminal 1 and one majorana mode connected to terminal 2 .    to minimize the number of free parameters we assume ballistic coupling through the point contact ,
so the matrix @xmath14 in the poisson kernel is @xmath227 .
the reflection amplitude @xmath228 at terminal 2 is tuned through zero by some external control parameter @xmath229 , typically magnetic field or gate voltage . near the transition ( conveniently shifted to @xmath230 )
this dependence can be conveniently parameterized by @xcite @xmath231 ( the width @xmath232 of the transition is system dependent . ) the corresponding coupling constant @xmath233 in eq .
then has an exponential @xmath229-dependence , @xmath234    the probability distribution of the andreev conductance ( in units of @xmath235 ) follows from eq .  , @xmath236^{-2}\nonumber\\ & = \tfrac{1}{2}[\cosh(\xi/\xi_0)+(g_{\rm a}-1)\sinh(\xi/\xi_0)]^{-2},\nonumber\\ & \qquad\qquad\;\;0<g_{\rm a}<2.\label{ppmg2}\end{aligned}\ ] ]
the delta - function limits are reached for @xmath237 ( keeping @xmath141 , because the entire system is topologically trivial ) .
right at the transition , at @xmath238 , the distribution is uniform in the interval @xmath239 .
the average conductance varies through the transition as @xmath240 see fig .
[ fig_gatransition ] .
, as the nanowire is driven through a topological phase transition ( controlled by a parameter @xmath229 , which can be thought of as the deviation of the magnetic field from the critical field strength ) .
the curve is the result for a single - channel ballistic point contact . ]
in conclusion , we have investigated a variety of observable consequences of the fact that the scattering matrix of majorana fermions is real orthogonal rather than complex unitary .
of particular interest is the identification of observables that can detect the sign of the determinant , since @xmath10 signifies the presence of a majorana bound state .
the obvious signal of such a zero - mode , a midgap peak in the density of states @xcite , is broadened by tunnel coupling to the continuum .
we find that the peak remains hidden if the coupling is ballistic ( unit transmission ) in even a single scattering channel .
the thermal conductance is likewise insensitive to the presence or absence of a majorana bound state , but the electrical conductance retains this sensitivity when the coupling is ballistic .
these results for the effect of a tunnel barrier on the midgap spectral peak are derived from the distribution @xmath0 of scattering matrix and time - delay matrix under the assumption of chaotic scattering , due to disorder or due to irregularly shaped boundaries .
the appropriate ensemble in the absence of time - reversal and spin - rotation symmetry has symmetry class d in the altland - zirnbauer classification @xcite .
chiral symmetry would change this to class bdi , in which multiple zero - modes can overlap without splitting @xcite .
the effect of chiral symmetry on the joint distribution @xmath0 is known for ballistic coupling @xcite  but not yet for tunnel coupling .
this seems a worthwhile project for future research .
this research was supported by the foundation for fundamental research on matter ( fom ) , the netherlands organization for scientific research ( nwo / ocw ) , and an erc synergy grant .
according to eq .  , the effect of a tunnel barrier on the average density of states follows directly once we know the average scattering matrix .
simple expressions can be obtained if we assume that the tunnel probabilities are mode - independent , @xmath156 for @xmath131 .
the scattering matrix of the barrier has the polar decomposition @xmath309 with @xmath310 .
the block structure corresponds to eqs .   and , in particular , @xmath311 .
from eq .   we obtain the scattering matrix @xmath1 of the quantum dot , @xmath312o_2,\label{stildes0relation}\ ] ] in terms of a matrix @xmath313 that is uniformly distributed in @xmath99 .
because the average of @xmath314 for any power @xmath315 is proportional to the identity matrix , we may write the average of @xmath1 in the form of a power series , @xmath316 if the average of @xmath283 would be over the entire unitary group , then all terms in the power series would vanish and we would simply have @xmath317 . but
averages over orthogonal matrices do not vanish , in the nontrivial way calculated instead of @xmath318 ) .
] by rains @xcite : @xmath319 we substitute eq .   into eq .   and sum the geometric series , to arrive at the average scattering matrix @xmath320 used in sec .
[ sec_dos ] to obtain the average density of states in class d. +   +   +   +      in a similar way we can derive the average scattering matrix in other symmetry classes . we give the result for class c.
the matrix @xmath283 then varies over the unitary symplectic group @xmath321 .
ref .   gives the required average : @xmath322 eq
.   still holds with the factor @xmath323 replaced by @xmath324 .
we thus find the average scattering matrix in class c , @xmath325\right).\label{saverageuc}\ ] ] for @xmath20 this gives @xmath326 , in agreement with ref .  .
the average density of states in class c follows from eq .
, where we again account for the doubling of the dimensionality @xmath327 : @xmath328\right)\nonumber\\ & = \frac{n}{(n+1)\delta_0}\,\left(1-\frac{1-\gamma}{n\gamma}[1-(1-\gamma)^n]\right).\label{rhoaveragec}\end{aligned}\ ] ] in the second equation we have substituted the ballistic class - c result from ref .  .
the tunneling limit @xmath162 gives a vanishing density of states , @xmath329 consistent with the class - c result for a closed system @xcite .
24 a. altland and m. r. zirnbauer , phys .
b * 55 * , 1142 ( 1997 ) .
j. alicea , rep .
. phys . * 75 * , 076501 ( 2012 ) .
m. leijnse and k. flensberg , semicond .
* 27 * , 124003 ( 2012 ) .
t. d. stanescu and s. tewari , j. phys .
matt . * 25 * , 233201 ( 2013 ) . c. w. j. beenakker , annu .
con . mat .
* , 113 ( 2013 ) .
s. das sarma , m. freedman , and c. nayak , arxiv:1501.02813 ( to be published in npj quantum inf . ) .
e. p. wigner , siam review * 9 * , 1 ( 1967 )
. m. l. mehta , _ random matrices _ ( elsevier , amsterdam , 2004 ) .
p. j. forrester , _ log - gases and random matrices _
( princeton university press , 2010 ) .
f. j. dyson , j. math .
* 3 * , 140 ( 1962 ) . c. w. j. beenakker ,
arxiv:1407.2131 ( to be published in rev . mod . phys . ) . t. j. krieger , ann .
* 42 * , 375 ( 1967 ) .
p. a. mello , p. pereyra and t. seligman , ann . phys . * 161 * , 254 ( 1985 ) . w. a. friedman and p. a. mello , ann .
* 161 * , 276 ( 1985 ) .
p. w. brouwer , phys .
b * 51 * , 16878 ( 1995 ) .
b. bri , phys . rev .
b * 79 * , 214506 ( 2009 ) .
v. a. gopar , p. a. mello , and m. bttiker , phys .
* 77 * , 4974 ( 1996 ) .
p. w. brouwer and m. bttiker , europhys .
lett . * 37 * , 441 ( 1997 ) .
m. bttiker , in _ time in quantum mechanics _ , edited by j. g. muga , r. sala mayato , and i. l. egusquiza ( springer , berlin , 2002 ) .
p. w. brouwer , s. a. van langen , k. m. frahm , m. bttiker , and c. w. j. beenakker , phys .
lett . * 79 * , 913 ( 1997 ) .
s. f. godijn , s. mller , h. buhmann , l. w. molenkamp , and s. a. van langen , phys .
lett . * 82 * , 2927 ( 1999 ) .
p. w. brouwer , k. m. frahm , and c. w. j. beenakker , phys .
78 , 4737 ( 1997 ) .
m. marciani , p. w. brouwer , and c. w. j. beenakker , phys .
b * 90 * , 045403 ( 2014 ) .
h. schomerus , m. marciani , and c. w. j. beenakker , phys .
lett . * 114 * , 166803 ( 2015 ) .
d. v. savin , y. v. fyodorov , and h .- j .
sommers , phys .
e * 63 * , 035202 ( 2001 ) .
sommers , d. v. savin , and v. v. sokolov , phys .
. lett . * 87 * , 094101 ( 2001 ) .
j. p. dahlhaus , b. bri , and c. w. j. beenakker , phys .
b * 82 * , 014536 ( 2010 ) . c. w. j. beenakker , j. p. dahlhaus , m. wimmer , and a. r. akhmerov , phys . rev .
b * 83 * , 085413 ( 2011 ) . c. mahaux and h. a. weidenmller , phys .
170 * , 847 ( 1968 ) ; _ shell - model approach to nuclear reactions _
( north - holland , amsterdam , 1969 ) .
l. fu and c. l. kane , phys .
* 100 * , 096407 ( 2008 ) .
m. titov , a. ossipov , and c. w. j. beenakker , phys . rev .
b * 75 * , 045417 ( 2007 ) .
j. j. m. verbaarschot , h. a. weidenmller , and m. r. zirnbauer , phys . rep . * 129 * , 367 ( 1985 ) .
a. r. akhmerov , j. p. dahlhaus , f. hassler , m. wimmer , and c. w. j. beenakker , phys .
. lett . * 106 * , 057001 ( 2011 ) .
e. p. wigner , phys . rev . * 98 * , 145 ( 1955 ) .
f. t. smith , phys . rev . * 118 * , 349 ( 1960 ) .
y. v. fyodorov and h .- j .
sommers , j. math . phys . * 38 * , 1918 ( 1997 ) . m. sh .
birman and m. g. krein , soviet math
. dokl . * 3 * , 740 ( 1962 ) .
e. akkermans , a. auerbach , j. e. avron , and b. shapiro , phys .
66 * , 76 ( 1991 ) .
n. lehmann , d. v. savin , v. v. sokolov , and h .- j .
sommers , physica d * 86 * , 572 ( 1995 ) .
a. f. andreev , sov .
jetp * 19 * , 1228 ( 1964 ) . k. schwab , e. a. henriksen , j. m. worlock , and m. l. roukes , nature * 404 * , 974 ( 2000 ) .
m. bocquet , d. serban , and m. r. zirnbauer , nucl .
b * 578 * , 628 ( 2000 ) .
d. a. ivanov , j. math .
phys . * 43 * , 126 ( 2002 ) ; arxiv : cond - mat/0103089 .
a. m. mathai , _ jacobians of matrix transformations and functions of matrix argument _ ( world scientific publishing , 1997 ) .
d. bagrets and a. altland , phys .
lett . * 109 * , 227005 ( 2012 ) .
j. li , g. fleury , and m. bttiker , phys .
b * 85 * , 125440 ( 2012 ) .
b. bri , phys .
b * 79 * , 245315 ( 2009 ) .
r. m. lutchyn , j. d. sau , and s. das sarma , phys .
lett . * 105 * , 077001 ( 2010 ) .
y. oreg , g. refael , and f. von oppen , phys .
105 , 177002 ( 2010 ) .
f. j. dyson , j. math .
phys . * 3 * , 1199 ( 1962 ) .
j. j. m. verbaarschot and t. wettig , ann .
nucl . part .
50 * , 343 ( 2000 ) .
i. c. fulga , f. hassler , a. r. akhmerov , and c. w. j. beenakker , phys .
b * 83 * , 155429 ( 2011 ) .
e. m. rains , probab .
theory relat .
fields * 107 * , 219 ( 1997 ) .