content_image/423990.jpg content_image/423997.jpg content_image/423998.jpg content_image/423999.jpg content_image/424006.jpg content_image/423996.jpg content_image/424008.jpg content_image/424007.jpg content_image/423995.jpg content_image/424011.jpg content_image/424012.jpg
greatly exceed the scales relevant to the superconducting state, the corresponding contributions to the current co- incide with those in the normal state. One can show that 

this regularization procedure leads to the same result as Eq. (30). The details of calculations are given in Ap- pendix A. The final expression for the current becomes 

<img src='content_image/423990.jpg'>

It is worthwhile to note that the problem of spurious divergent terms in the expression for the supercurrent is rather general. In particular, it was discussed (and resolved similarly) for the superfluid excitonic current in graphene bilayers $^{18}$. 

We evaluate Eq. (31) for low temperatures, T ≪ | Δ | . Since E$_{↑}$, E$_{↓}$ > | Δ | , the first line in Eq. (31) vanishes at T = 0. The supercurrent becomes 

<img src='content_image/423997.jpg'>

For µ ≫ | Δ | we have 

For µ ≪ | Δ | we find 

This result formally holds within the linear approxima- tion which assumes v$_{F}$ k ≪ µ . Therefore, in Eq. (34) one would have to put k → 0 first and then assume µ ≪ | Δ | . In the next section we demonstrate that Eq. (34) is in fact always valid provided | µ | ≪ | Δ | and v$_{F}$ k ≪ | Δ | irrespectively of the relation between v$_{F}$ k and µ . 

Consider now the limit of small µ when the zero order state is degenerate because E (0) 1 = E (0) 3 = E$_{0}$ and E (0) 2 = E (0) 4 = − E$_{0}$ where 

In what follows we demonstrate that despite the absence of the Doppler term in its usual form, there still is a finite linear in k supercurrent down to zero temperature (contrary to the result of Ref. 7). 

<img src='content_image/423998.jpg'>

<img src='content_image/423999.jpg'>

## B. Low doping limit 

<img src='content_image/424006.jpg'>

The true wave functions satisfy 

<img src='content_image/423996.jpg'>

We assume v$_{F}$ k, µ ≪ E$_{0}$ . Let us expand the true wave function into the zero order orthonormal wave functions ˇ ψ (0) α , 

satisfying the zero-order equation 

The zero-order wave functions have now the form of Eqs. (13), (15) with u$_{↑}$ = v$_{↓}$ ≡ u and u$_{↓}$ = v$_{↑}$ ≡ v where 

<img src='content_image/424008.jpg'>

The expansion coefficients satisfy 

where H$_{γβ}$ ≡ 〈 ˇ ψ (0)+ γ ˇ H (1) ˇ ψ (0) β 〉 . Consider the state α = 1. Since the difference E$_{1}$ − E (0) 1 = E$_{1}$ − E (0) 3 ≡ δE$_{1}$ in Eq. (38) is small, the coeffi- cients C$_{12}$ an C$_{14}$ are proportional to the perturbation, while C$_{11}$ and C$_{13}$ are of the order unity. The coefficients C (0) 11 , C (0) 13 in the leading approximation satisfy the sec- ular equations (B5), (B6) (see Appendix B) which yield δE$_{1}$$_{,}$$_{3}$ = ∓ ˜ E$_{1}$ , 

<img src='content_image/424007.jpg'>

<img src='content_image/423995.jpg'>

<img src='content_image/424011.jpg'>

<img src='content_image/424012.jpg'>

