arxiv:cond-mat/ v1 [cond-mat.supr-con] 14 Sep 2006

Transcription

1 Derivation of the Ginzburg-Landau euations of a ferromagnetic p-wave superconductor arxiv:cond-mat/06094v cond-mat.supr-con 4 Sep 006 E. K. Dahl and A. Sudbø, Department of Physics, Norwegian University of Science and Technology, N-749 Trondheim, Norway Centre for Advanced Study at the Norwegian Academy of Science and Letters, Drammensveien 78, N-07 Oslo, Norway Dated: February 6, 008 We derive a Ginzburg-Landau free energy for a p-wave ferromagnetic superconductor. The starting point is a microscopic Hamiltonian including a spin generalised BCS term and a Heisenberg exchange term. We find that coexistence of magnetisation and superconductivity depends on the sign of the energy-gradient of the DOS at Fermi level. We also compute the tunneling contribution to the Ginzburg-Landau free energy, and find expressions for the spin-currents and osephson currents across a tunneling junction separating two ferromagnetic p-wave superconductors. PACS numbers: 7.5.Mk,74.0.De, r,74.0.-z I. INTRODUCTION In recent years, experimental evidence for materials that simultaneously feature superconductivity and ferromagnetism, has transpired from a number of independent studies,,. This has triggered much research into this class of materials, widely known as ferromagnetic superconductors FMSC. There are still many unanswered uestions concerning the properties of these systems featuring coexisting spontaneously broken symmetries broken U and O symmetries. Recently, the superconductivity in ZrZn has been found to be surface sensitive, and might in fact not be a bulk property of the material 4. A common feature for all the ferromagnetic superconductors is that they are superconducting only in the ferromagnetic phase. That is, when pressure is increased in such a way as to destroy the ferromagnetic phase transition, then superconductivity, if present, is also lost. Ferromagnetic superconductivity is of great interest from a theoretical point of view. It offers a laboratory for studying condensed matter systems with multiple spontaneously broken symmetries, and their interplay. From a technological point of view, it is also hoped that heterostructures with FMSC will give rise to novel transport effects involving both the charge and the spin of the electron, which may have the potential for being exploited in novel types of devices. Recently, structures containing FMSC have been investigated theoretically 5. This work, however, considers tunneling effects between a spin-singlet superonductor in a Fulde-Ferrel-Larkin- Ovchinnikov state 6,7 coexisting with helimagnetic order. On the other hand, little is known of the tunneling effects and interplay between ferromagnetism and superconductivity in spin-triplet superconductors coexisting with uniform magnetic order. It is the purpose of this paper to derive a Ginzburg-Landau model for such a system starting from a reasonable microscopic Hamiltonian, and then apply the results to computing tunneling effects both in the charge and spin sector, and thus to elucidate the interplay between the broken U and O symmetries form such systems. Previously, the Bogoliubov-de Gennes euations for such systems have been studied in some detail 8. Phenomenological models designed to describe FMSC were proposed soon after their experimental discovery 9,0. These theories can be written down from symmetry considerations of the two order parameters in the problem, i.e. the theories make no reference to the underlying microscopic physics. In ferromagnetic superconductors it is believed that it is the same electrons that are responsible for both superconductivity and ferromagnetism. In such a scenario, the Cooper-pairs must have a magnetic moment, which means that a FMSC with uniform magnetic order must be a spin triplet superconductor. Microscopic theories attempting to explain the positive interaction between electrons with aligned spins have traditionally evolved around spin mediated interactions, i.e. magnon exchange,. In Ref. the phenomenological model of Refs. 9,0, was employed in order to find stable solutions and the regions were these solutions of the model are stable. The phenomenological model used in Refs. 9,0 is given by f ψ, M = f S ψ + f F M + f I ψ, M + H, were ψ is a three dimensional complex vector describing superconducting order, B = H +4πM = A is the magnetic induction, H the external magnetic field,and A is the electromagnetic vector potential. The term f S ψ describes superconductivity when H = M = 0, and is given by f Sψ = f grad ψ + a s ψ + bs ψ 4 + us with ψ + vs ψ i 4, f grad ψ =K D i ψ j D i ψ j + K D i ψ i D j ψ j + D i ψ j D j ψ i + K D i ψ i D i ψ i, where summation over the indices i, j is assumed and the symbol i= D i = i + ea i, 4 x i

2 of covariant differentiation is introduced. For a discussion of the superconductivity part of the free energy see Refs. 0,. The term f F M in describes the ferromagnetic ordering of the material and is given by, f F M = c f j= j M j + a f T f M + b f M4. 5 This is the standard expression for an isotropic ferromagnet. Furthermore, f I ψ, M describes the interaction between the order parameters for ferromagnetism and superconductivity, and is given by f I ψ, M = iγ o M ψ ψ + δm ψ. 6 In this paper, we will start out with a microscopic Hamiltonian that essentially reproduces the Ginzburg- Landau model proposed in Ref. 0. As a byproduct of such a derivation from a more microscopic model, we obtain analytic expressions for the coefficients in. The microscopic theory we start with is not trying to explain the mechanism of the non-unitary spin-triplet pairing of the electron gas. We simply assume the existence of some spin generalised attractive BCS term in the Hamiltonian, in addition to a Heisenberg ferromagnetic exchange term. As an application of the Ginzburg-Landau euations, we look at two FMSC in tunneling contact. We compute the free energy of the coupling and further find the osephson current and the spin current in the direction perpendicular to the plane spanned by the two magnetisation directions. II. MODEL We start out with a Hamiltonian that is given by three terms, one describing free electrons, one BCS term to account for spin-triplet superconductivity, and finally a Heisenberg ferromagnetic exchange term to account for itinerant ferromagnetism. The Hamiltonian of the system is then given by, Ĥc, c = k,σ + ε k c k,σ c k,σ k,k, V k,k c k+/,α c k+/,β c k +/,βc k +/,α + γs S, 7 where c k,σ and c k,σ annihilates and creates an electron in the state k, σ respectively. Here, S k is the usual spin operator given by, S k =,α,β c k+,α σ αβc,β and σ is the Pauli matrices. γ = δ exp i δ is the structure factor of the underlying lattice, and δ is a nearest neighbour vector. A summation over repeated Greek indices is implied. We are interested in calculating the partition function of the system, formally given by Z = Tr e βĥc,c = e βf, 8 were F is the exact free energy of the system. Introducing fermion coherent states, ξ k,σ, and performing a Hubbard- Stratonovich decoupling of the two last terms in 7, we arrive at an effective action in Euclidean time which reads, S eff = β 0 dτ k,σ + k, α,β + k, ξ k,σ τ + ε k ξ k,σ + ξ k,σ τ ε k ξ k,σ ξ k+/,α σ M αβ ξ k /,β ξ k /β σ M αβ ξ k+/α αβ k, ξ k+/,βξ k+/,α + αβ k, ξ k+/,β ξ k+/,α γ M M k,k, αβ k, V k,k βαk, 9 Here, βα k, and M are auxiliary fields in the functional integral representation of the partition function, and physically represent the spin-triplet pairing fields and magnetisation, respectively. The above action is Gaussian in the fermionic fields and hence the integral over the fermionic fields may be performed exactly. We next introduce a Majorana basis

3 φ k = ξ k, ξ k, ξ k, ξ k,, the effective action may now be written, S eff = β 0 dτ φ k G φ k k, k,k, γ MM tr k, V k,k k,. 0 Here, tr denotes trace over spin indices. Furthermore, G is a 4 4 matrix given by, G = G 0 Σ, where G 0 = diag iω n + ε k, iω n + ε k, iω n ε k, iω n ε k describes a free electron gas. Here, Σ describes the interaction and pairing correlations, and is explicitly given by, M D Σ = D, M where M = M σ is a matrix order parameter describing magnetisation and D = d σ = i σ y is a matrix order parameter describing triplet superconductivity. The integral over φ is Gaussian, and hence it may be performed analytically. The integral produces a fermion determinant which may be included in the exponent, i.e. in the Ginzburg-Landau free energy, which is given by βf GL = TrlnG β dτ 0 γ M M + k,k, tr k, V k,k k,, where Tr implies a trace over all variables. The GL free energy is defined by Dξ, ξ e S effξ,ξ = e βfgl. III. SECOND ORDER TERM The second order term in the expansion of the trace is given as E = Tr G 0ΣG 0 Σ = n G 0 ΣG 0 Σ n = k,k ωn n G 0 k Σ k,k G 0 k Σ k,k. 5 The last euality comes from inserting completeness relations and noting that G 0 is local in k-space. We have introduced the notation G 0 k = k G 0 k and Σ k,k = k k Σ k. After changing variables, k the second order term is E = g0,k+/ g 0,k / + c.c M M k, ωn k+/ k /, g 0,k+/ g0,k / + c.c d k,, 6 where g 0,k = iω n + ε k. The trace of the superconducting order parameter is given by id tr k, k, = tr id µk, σ ν µσ y k, σ ν σ y =d k, d k,. 7 The complete second order term of the Ginzburg-Landau free energy is therefore given by F = k, ωn δk,0 g0,k+/ g 0,k / + c.c M M γ β g0,k+/ g0,k / β + c.c + /V d k, 8 The trace over ln G may formally be rewritten as TrlnG = TrlnG 0 + Tr ln G 0 Σ, the first term, describing free theory, is neglected in the following. The second term is assumed small, and hence we may expand the logarithm to obtain Tr ln G 0 Σ Tr G 0 Σ + G 0ΣG 0 Σ + G 0ΣG 0 ΣG 0 Σ + 4 G 0ΣG 0 ΣG 0 ΣG 0 Σ + OG 0 Σ 5 = E E E... 4 The first term in the expansion, E = TrG 0 Σ, is zero since the Pauli matrices are traceless.the terms E and E are second order in the ordering fields βα k, and M, and we now proceed top discuss these terms in turn. Now is assumed to be small, and hence we may expand ε k+/ ε k + / ε k k + O = ε k + / v F, the inverse electron propagator is to first order in given by g 0,k+/ iω n + ε k + / v F. Hence, we find for the part containing magnetisation g 0,k+/ g 0,k / + gk+/ g k / =R iω n + ε k /v F iω n + ε k 4 9 when keeping terms to second order in. In addition, we expand the structure factor γ 6, and assume for simplicity a cubic lattice in three dimensions. Now, we substitute k N ǫc 0 ǫ c dω 4π, where N 0 is the

4 4 density of states at Fermi level and ǫ c is some cutoff. The homogeneous part of the magnetisation is given by F,m c = β ǫm β 4N 0 tanh + M M. 6 0 The part containing derivatives is given by F,m d = β N0 v F β ǫ m ǫm β tanh tanh ǫ mβ 7 + M M 6 Similarly, we find for the superconducting order parameter, when introducing d µ k, = A µi ˆk i 0 and the indices µ and i run from to, that the coefficient in front of the term involving the superconducting order parameter is given by Rg 0,k+/ g0,k /ˆk iˆkj k,ω n dω N 0ˆkiˆkj 4π ω ω n n + ξ v F / l mˆklˆkm ωn + ξ + / v F iω n + ξ / v F iω n + ξ. Here, we have used v F = v F ˆk. The terms linear in will integrate to zero when integrated over the angles Ω. We start with the constant term independent of α c = ω n = βn 0δ i,j dω N 0ˆkiˆkj 4π ωn + ξ ε n 0 n + / where ε is an energy cut-off. In total the constant second order term is given by F,S c = V 4N 0 ε tr AA. n + / n 0 4 After using the result for the critical temperature 4, the term may be written on the form F c,s = 4N 0 T T c T c tr AA. 5 The coefficient for the second order term containing derivatives is given by α d = N 0 v F / i j ω n dω ˆk iˆkjˆklˆkm 4π ωn + ξ = N 0 v F / i j /5δ ij δ lm + δ il δ jm + δ im δ jl ωn + ξ ω n = π/n 0 v F / i j /5δ ij δ lm + δ il δ jm + δ im δ jl sgnω n ωn ω n = 7 i j v F N 0β ζ 40π δ ij δ lm + δ il δ jm + δ im δ jl, 6 so that the total second order term containing derivatives is given by F,S d = 7v F N 0β ζ 0π tr AA + i A µi j A jµ + j A µi i A jµ. 7 Here we have assumed an isotropic uadratic dispersion relation, hence the coefficients in are all the same to this approximation. If one introduces an effective inverse mass tensor /m i,j one can use the result above by only replacing i j v Fi,j, where v Fi,j = //m i,j + /m j,i, and remove v F. In this way space is no longer isotropic. IV. THIRD ORDER TERM In this section we will find a coupling between the order parameters for magnetisation and superconductivity. The third order term will, however, only be nonzero when the system is in a non unitary phase, i.e. when d d 0. In a ferromagnetic superconductor where the same electrons are responsible for superconductivity and ferromagnetism, the system, however, has to be in a

5 5 non-unitary state. If the system is in a unitary state the superconducting order parameter is invariant under the time inversion operator, and hence there can not be any magnetism associated with it. The third order term in the expansion of the fermion determinant is given by E = G 0 k Σ k,k G 0 k Σ k,k G 0 k Σ k,k. k,k,k Via a change of variables, k+/ k / k+/ k / k k 8 k+/ k, k / k k and k along with the constraint, + + = 0, a multipication of the matrices yields the following expression for the third order term E = 4 g0,k g0,k g 0,k + c.c k, i id k, d k, M, 9 were the combination m = id d is interpreted as the average magnetisation due to the Cooper-pairs 4,5. The prime on the sum in Euation 9 denotes sumation over configurations with the restriction δ + +. In the coefficient we have neglected the dependence, since expanding in powers of amounts to finding derivatives. The coefficient is given by, when introducing d µ = A µiˆki, α /! =4R ˆk iˆkj ω ω n + ε n,k k iω n + ε k ε dω 4N 0 ξ ˆkiˆkj 4π ωn + ξ iω n + ξ ω n ε =N 0 βδ ij/ n 0 n + /. 0 The final summation in the above expression is proportional to the logarithm of the cut-off freuency. The total third order term is thus given by F =,, α! iε µνλa µi A νi M λ. were α = α /β. As the order parameters M and d do not couple to second order in the GL free energy, the coupling in the third order term is expected to be of crucial importance as to whether the two order parameters will coexist in the system or not. Coexistence is favoured by the system if there is an energy gain by having M and d finite simultaneously. Hence, coexistence of magnetism and superconductivity depends on the sign of α, which again is given by the gradient of the DOS at Fermi level. In the system considered here, there will be a ferromagnetic coupling between magnetism and the spin magnetism of the Cooper-pairs if N 0 > 0, i.e the gradient of the DOS should be positive at the Fermi level for the spin-up sheet if coexistence of FM and SC is preferred by the material. Comparing E. with the model E. 6, we find the value of the coefficient γ 0 when we assume that ψ is in a p-wave state. The sign dependence of N o is however more general, since a triplet state must necessarily be an odd function of k. From 9, we observe that in an expansion of the DOS around Fermi level, Nξ N 0 + ξn 0 + ξ /N , only the odd terms of the expansion give a contribution to the coefficient. Hence, to lowest order γ 0 N 0. We next determine the fourth order terms in the order-paarameter expansion of the GL free energy. V. FOURTH ORDER TERM In this section, we will find three different types terms, namely a term involving only magnetisation, one onoly involving superconductivity, and finally one term involving a coupling between the magnetisation M and the superconducting orderparameter d. In total, we will find five different independent terms for a p-wave ferromagnetic superconductor. The fourth order term of the expansion of the fermion determinant is given by E 4 = 4 G 0 kσ k,kg 0 kσ k,k k,k k,k 4 G 0 k Σ k,k 4 G 0 k4 Σ k4,k. As before, we change variables and multiply out the matrices and take the trace. In addition, we make the approximation of neglecting the dependence on i in g 0,ki. A. Terms containing only magnetic order parameter The term involving the magnetisation only is given by βf 4,m = M M M M 4 i R g 0,k g 0,k g 0,k g 0,k. k,ωn Here, the coefficient of the M -factors is given by the trace over the electron propagators, α 4m =R k,ω n g 0,k 4 = N 0 R ω n ǫm = β ǫm β tanh 0 tanh ǫm β 4 iω n + ξ 4

6 6 Thus, the fourth order term involving magnetisation only, is given by F 4,m = i α 4m where α 4m = α 4m /β. M M M M 4 5 B. Terms containing only superconducting order parameter The term involving the superconducting order parameter only will contain five different terms. Again neglecting the dependence in g 0,k, we arrive at βf 4,S = R ˆkiˆkjˆklˆkm ω i k,ω n + ε n k Aµi A µj A νla νm A µia µj A νl A νm. 6 The coefficient is given by α 4S = N 0 5β ω ω n n + ξ = 7N 0ζ 0π β 8 In total, the fourth order term involving superconductivity alone is given by, F 4,S = α 4S i tr AA T + tr AA + tr AA T AA T + tr AA tr AA AA, 9 where α 4S = 7ζN0 5π β. This part of the free energy thus consists of five independent terms 6. Now let the sum over k go over to an integral over ξ and the angles θ and φ. The integral over the angles produces Kronecker-δ s in Latin indices, dω 4π ˆk iˆkjˆklˆkm = 5 δ i,jδ l,m + δ i,l δ j,m + δ i,m δ j,l. 7 C. Terms containing magnetic and superconducting order parameter The fourth order coupling term between the magnetisation and the superconducting order parameter is given by F 4,Sm = 8g 0,k g 0,k g0,kg 0,k M d M d 4 M M d d 4 k, i,ω n 6g 0,k g 0,k g 0,k g 0,k d d 4 M M = 7N 0 ζβ π i M µ A µj M ν A νj 40 where the constant in the free energy is α 4Sm = 7N 0 ζβ /π. This coupling term between magnetism and superconductivity differs from the fourth order coupling in E. 6. In E. the coupling term is proportional to M ψ, whereas we find a term proportional to M M ψ. In addition, the coefficient is positive rather than having an indefinite sign as commented on in Ref.. We have, however, assumed that the system is in a p-wave state. In the a case of a general spin-triplet state, both types of fourth order coupling terms between magnetism and superconductivity may be present, as also seen from E. 40. D. Complete fourth order term In total, the fourth order term is thus given by

7 7 F 4 = / α 4m M M M M 4 + α 4Sm M µ A µj M ν A νj i +α 4S traa T + traa + tr AA T AA T + tr AA tr AA AA 4 All in all the Ginzburg-Landau free energy for a ferromagnetic superconductor is given by F GL = d αs T r traa + β S D tr AA + D i A µi D j A jµ + D ja µi D i A jµ + α mt M M + β m D M M + α! iε µνλa µi A νi M λ + α 4S traa traa + traa T AA T + traa traa AA + α 4m M M + α 4ms M µ A µi M ν A νi, 4 where D = + iea, D = + iea and A is the electromagnetic vector potential. VI. TUNNELING As an application of the Ginzburg-Landau theory derived above, we will consider tunneling between two ferromagnetic p-wave superconductors. The Hamiltonian we use is the common choice when studying tunneling between systems in euilibrium H = H L + H R + H T, 4 where H LR is given by 7 and H T = T α,β c k,α d p,β + T α,β d p,β c k,α. 44 For a general Hamiltonian, it is straightforward to show that the part of the Ginzburg-Landau free energy that contains the tunneling elements, is given by F = β Tr ln G LT G R T, 45 where G LR are the Green s functions for the leftright subsystem and T is a tunneling matrix, T 0 T = 0 T 46 k, p with T T = T T T. 47 Furthermore, we will assume for simplicity that the magnetisation of our system is weak, hence we assume that the magnetisation is homogeneous. Introducing the basis φ k = ξ k,, ξ k,, ξ k,, ξ k, and choosing the uantisation axis along the magnetisation vector at each side of the junction, we have the following effective action, G LR = S eff = S L + S R + S T 48 and the inverse Green s functions are given by g LR LR LR g LR, 49 with g LR,αβ = iω n + ε k αm LR δ α,β. The direction of the magnetisation does not enter in the Greens functions for the isolated systems on the left and right, the angle between the magnetisations, and hence the uantisation axes, only enters through the tunneling elements since we have applied different rotations on the left and right. This difference shows in the part of the Hamiltonian where the left and right systems couple. We now want to calculate the Ginzburg-Landau free energy for the junction, F. To this end, we expand 45 to the lowest order in the tunneling elements, F β TrG LT G R T. 50 Next, we use the Dyson euation to find an expansion of G i.e. G G 0 + G 0 ΣG were g 0 G 0 = 0 g 5

8 8 and Σ = Since we are assuming that the magnetisation is homogeneous, G 0 is local in k space, Σ however is nonlocal as we do not want to impose any restriction on the superconducting pairing state. From the free energy of the junction we can find the current of various uantities associated with tunneling. To show this we introduce a Hamiltonian H, p which is a function of the canonically conjugate variables, p, and recall that the velocity operator is given by ˆ = H p. Consider now the uantity F p = β p ln Tre βh = Z Tr H p e βh =. Hence we are easily able to to find the current over the junction of for instance electrons. Since φ, N = i, the current is than given by I = e Ṅ = e φ F and the contribution to the spin current coming from tunneling across the barrier, is given by Ṡn = µ B θ F, where θ, S n = i, and θ is the angle in the plane perpendicular to the direction ˆn. A. Ferromagnetic case We start with the term corresponding to a junction between two ferromagnetic metals, F 0 = β = β ωn,ων ωn,ων tr G 0k T G 0p T tr g k T g p T +g k T k, p g p Tt k, p. 5 We observe that upon letting k k and p p the second term above is the complex conjugate of the first one, hence the free energy is two times the real part of the first term, F 0 = β R g + k g+ p T ++ +g + k g p T + + g k g p + g k g+ p T T Here + and - means parallel or anti parallel to the magnetisation respectively. The angle between the magnetisations on the left and right enters through the tunneling elements T αβ = / + αβ ˆM L ˆMR T, when we assume that the spin in preserved across the tunnel barrier. The potential difference between the left and right subsystems are taken care of through the tunneling elements i.e. we assume that T is nonlocal in ω-space with a boson freuency ω ν, later we do an analytical continuation ω ν ev +iδ. The structure of the Green s functions is, g α k gβ p = i ω n ω ν + ε kα iω n + ε pβ, 55 were ε kα = ε k αm L and similarly for the right subsystem. The zero th order term of the free energy may be written in two parts, one part independent of the direction of the magnetisations and one part proportional to the dot product of the directions of the magnetisations, that is, to the cosine of the angle between the magnetisations, F 0 = F 0 and = β = T F 0 + ˆM L ˆMR F 0,M. With R g + k g+ p + g k g p + g + k g p + g k g+ p T k p Nξ k Nξ p fξk M L fξ p M R ev + ξ k ξ p + M R M L + fξ k M L fξ p M R ev + ξ k ξ p + M L M R + fξ k M L fξ p + M R ev + ξ k ξ p M L M R + fξ k + M L fξ p M R ev + ξ k ξ p + M L + M R F 0,M = g + k β g+ p + g k g p g + k g p g k g+ p = T k p Nξ k Nξ p fξk M L fξ p M R ev + ξ k ξ p + M R M L + fξ k M L fξ p M R ev + ξ k ξ p + M L M R fξ k M L fξ p + M R ev + ξ k ξ p M L M R fξ k + M L fξ p M R ev + ξ k ξ p + M L + M R Here Nξ is the density of states and fx = / + e βx is the Fermi-Dirac distribution function. The spin current Ṡn, were n denotes the direction perpendicular to the plane spanned by the magnetisation vectors, is now easily found by taking the derivative with respect to the angle between the magnetisations, i.e. Ṡn = µ B 0 θ F 0 0 =µ B F,M sin θ 58 This is precisely the same result as found previously in Ref. 7 via a different route.

9 9 B. Ferromagnetic superconducting case In this subsection, we calculate the free energy which is first order in the superconducting gap function on both the left and right side, i.e. we consider the ferromagnetic superconducting state. We want to find the osephson current and also the two particle contribution to Ṡn. The free energy which can give rise to osephson current is, F = β Tr G 0LΣ L G 0L T G 0R Σ R G 0R T = β tr G0k Σ k,k G 0k T k,pg 0p Σ p,p G 0p T. 59 k We now change variables, k k+/ k / and p p p+ / p / were is the centre of mass momentum of the Cooperpairs on the leftright, and hence it is small. Expanding in amounts to an finding derivatives, we are not interested in derivatives and neglect the dependence in all terms except Σ, since we want to have the possibility + 0. With these assumptions the free energy reads, F = tr g k ˆk, g k β T k, p g p ˆp, T +g k ˆk, g k T g p ˆp, g pt t k, p. 60 Here we observe that upon changing the sign of all impulses in the second term and using that ˆk = ˆk for a triplet superconductor, the second term is just the complex conjugate of the first one. Hence, the free energy is just two times the real part of the first term. After some straightforward algebra, we obtain F = T, + + ˆM L ˆMR A F + A F ˆM L ˆMR A + ˆM L ˆM R + A F + F cos φ + A h + A h + A h + h sin φ + A were we have introduced the notation F + A F A F + F cos φ A h + A h A h + h sin φ F + F F + F F + F F A + A A + A h + h A h + h + A h + h A h + h sin φ A λγ αβ, Ω k >0 dω k 4π Ω p>0 dω p 4π + F cos φ 6 L αβ R λγ TΩk,Ω p, 6 F λγ αβ N ξk + αm L + Nξ k + βm L Nξ p + λm R + N ξ p + γm R k p 4 ξ k α β ML ξ p + λ γ MR fξ k fξ p ev + ξ k ξ p, 6 and h λγ αβ π 4 k fξ k fξ k + ev Nξ k + αml + ev Nξ k + γmr + N ξ k + λm R ξ k + α β ML + ev ξ k λ γ MR + Nξ k + βm L Nξ k + λm R + ev + N ξ k + γm R ev ξ k α β ML ξ k + λ γ MR + ev. 64 In addition, we have assumed that φ L αβ φr λγ = φ is independent of the spin indices, this can be done since any coupling between the different components of the superconducting order parameter will give a phase locking 8. Further since T k, p = 0 9 the sum over k and p i Euation 59 may be restricted to positive values only, and hence there are no problems involved in splitting the superconducting orderparameter into an amplitude and a phase in Euation

10 6. The osephson current is now found by simply taking the derivative of the free energy with respect to the phase difference φ. This leads to the floowing expression for the osephson current + ˆM L ˆMR A F + A F + A F + F sin φ I = e T, + ˆM L ˆMR A A h + A h + A h + h cos φ F + A F A F + F sin φ A h + A h A h + h cos φ + ˆM L ˆM R A F + F A F + F + A F + F A F + F sin φ A h + h A h + h + A h + h A h + h cos φ. 65 Similarly, we find the following expression for the two-particle contribution to the spin current Ṡn = T sin θ A F + A F + A F + F + F + F A F A F cos φ, + A h + A h + A h + h + h + h A h A h sin φ + cosθ + A F + F A F + F + A F + F A A h + h A h + h + A h + h A h + h sin φ F + F 0 cos φ. 66 The spin current is not well defined in the limit when θ 0, since we are calculating the current along the direction perpendicular to the plane spanned by the two magnetisation directions. We observe that the total spin current is the sum of E. 58 and 66. Hence, E. 66 is a contribution originating with the interplay between superconductivity and magnetism. It is seen to disappear when superconductivity is lost. Moreover, we observe that there is a term in euation 66 that is proportional to cos θ. In osephson currents there also exist a term that is proportional to cos φ, which however is dissipative and vanishes if the voltage across the junction is set to zero. In the case of the spin current, we find that the condition that must be fulfilled for the cosine term to go to zero, is given by L αβ = R αβ, M L = M R and ev = 0. Furthermore, it is seen that the cosine term is associated with flipping one electron from spinstate α on the left to α on the right, hence the cosine term is proportional to L αα R αα and similarly with the superscripts interchanged. VII. CONCLUDING REMARKS We have derived a Ginzburg-Landau functional from a microscopic Hamiltonian consisting of three terms, free fermigas, an Heisenberg term and a spin generalised BCS term. We find two order parameters, local magnetisation and the superconducting gap. We expand to fourth order in the order parameters. The lowest order coupling between the two order parameters is a third order term, this term is only non zero when the material is in a non unitary superconducting state i.e. d d 0. From an exact calculation of the coefficient of the third order term we find that coexistence of ferromagnetism and superconductivity is enhanced if the gradient of the density of states, d Nξ ξ=0, is positive at Fermi level. We have also computed the tunneling contribution to the Ginzburg-Landau free energy. From this, we have found expressions for spin- and charge currentws in the spin- and charge channels across a tunneling junction separating two ferromagnetic p-wave superconductors. Superconductivity coexisting with ferromagnetism can be triggered by the magnetisation as a result of the presence of a third order term, i.e. superconductivity appears at a higher temperature than if magnetisation were not present. The sign of the third order term is given by the the sign of the gradient of the DOS at Fermi level. Since the Fermi level may be tuned by applying pressure, it may be possible to change the sign of the third order term by applying pressure. Acknowledgments This work was supported by the Research Council of Norway Grants No. 5858/4, No /4 NANOMAT. E.K.D thanks NTNU for a research fellowship.

11 C. Pfleiderer, M. Uhlarz, S. M. Hayden, R. Vollmer, H. v. Löhneysen, N. R. Bernhoeft, and G. G. Lonzarich, Nature 4, D. Aoki, A. Huxley, E. Ressouche, D. Brathwaite,. Flouuet,.-P. Brison, E. Lhotel, and C. Paulsen, Nature 4, S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, R. K. W. Haselwimmer, M.. Steiner, E. Pugh, I. R. Walker, S. R. ulian, P. Monthoux, et al., Nature 406, E. A. Yelland, S. M. Hayden, S.. C. Yates, C. Pfleiderer, M. Uhlarz, R. Vollmer, H. v. Lhneysen, N. R. Bernhoef, R. P. Smith, S. S. Saxena, et al. 005, cond-mat/ I. Eremin, F. S. Nogueira, and R.-. Tarento, Phys. Rev. B 7, P. Fulde and R. A. Ferrel, Phys. Rev A. I. Larkin and Y. N. Ovchinnikov, Sov. Phys. ETP B.. Powell,. F. Annett, and B. L. Gyorffy,. Phys. A 6, K. Machida and T. Ohmi, Phys. Rev. Lett 86, M. B. Walker and K. V. Samokhin, Phys. Rev. Lett D. Fay and. Appel, Phys. Rev. B, N. Karchev, Superconductivity reasearch at the leading edge Nova Science Publishers, N.Y., 00, chap. Itinerant Ferromagnetism and superconductivity. D. V. Shopova and D. I. Uzunov, Phys. Rev. B V. P. Mineev and K. V. Samokhin, Introduction to Unconventional Superconductivity Gordon and Breach, A.. Leggett, Rev. Mod. Phys. 47, G. Barton and M. A. Moore,. Phys C 7, F. Nogueira and K.-H. Bennemann, Euro. Phys. Lett 67, A.. Leggett, Prog. Theor. Phys. 6, C. Bruder, A. van Otterlo, and G. T. Zimanyi, Phys. Rev. B 5, Here we are assuming that the superconducting order parameter is in a p-wave state