Molecular identification and chromosomal localization of new powdery mildew resistance gene Pm11 in oat Title: The Fulde-Ferrell-Larkin-Ovchinnikov state in pnictides

Fe-based superconductors (FeSC) exhibit all the properties of systems that allow the formation of a superconducting phase with oscillating order parameter, called the Fulde–Ferrell–Larkin–Ovchinnikov (FFLO) phase. By the analysis of the Cooper pair susceptibility in two-band FeSC, such systems are shown to support the existence of a FFLO phase, regardless of the exhibited order parameter symmetry. We also show the state with nonzero Cooper pair momentum, in superconducting FeSC with ∼ cos(kx) · cos(ky) symmetry, to be the ground state of the system in a certain parameter range.


Introduction
At low temperatures the orbital pair breaking effects are smaller in magnitude than the Pauli paramagnetic effect, so that superconductivity survives up to the Pauli limit-a phase with oscillating order parameter (called the Fulde-Ferrell-Larkin-Ovchinnikov phase or FFLO in short) [1][2][3] can be more stable than a phase with a constant order parameter (the Bardeen-Cooper-Schrieffer phase, or BCS in short). In this case, Cooper pairs may be formed with non-zero total momentum between Zeeman-split parts of the Fermi surface.
In this paper, making use of the Cooper pair susceptibility and the minimization of free energy of the system, we discuss the possible appearance of the FFLO phase in pnictides. In Sect. 2 we describe the selected model of FeSC, in Sect. 3 we present our methods. In Sect. 4 we illustrate and discuss our numerical results. We summarize the results in Sect. 5.

Theoretical Model
The FeSC system is described using a two-orbital per site model, with hybridization between the d xz and d yz orbitals. We adopt the band structure proposed in Ref. [23] and assume that the external magnetic field is parallel to the plane. The Hamiltonian of the system in momentum space takes the following form: where c † αkσ (c αkσ ) is the creation (annihilation) operator of a particle with momentum k and spin σ in the orbital α. T αβ kσ is the kinetic energy term of a particle with momentum k changing the orbital from β to α, μ is the chemical potential and h is the external magnetic field. The hoppings have magnitudes: (t 1 , t 2 , t 3 , t 4 ) = (−1.0, 1.3, −0.85, −0.85), in units of |t 1 |. At half-filling, a configurations with two electrons per site requires μ = 1.54|t 1 |. Our choice of the parameter set is motivated by the fact that it reproduces the same Fermi surface structure as the local-density approximation calculations of band structure [24][25][26][27][28].

Methods
We introduce a superconducting pairing between the long-lived quasi-particles in bands ε = ± [29]. To determine the possibility of formation of the FFLO phase, we turn our attention to the static Cooper pairs susceptibility: where . . . r is the retarded Green's function and The operator d εiσ in real space corresponds to the operator d εkσ in momentum space. The Factor ϑ(j − i) defines the OP symmetries ( Fig. 1)-for example for d x 2 −y 2 -wave pairing, ϑ(δ) is equal to 1 (−1) for δ = ±x (±ŷ) and zero otherwise. In momentum space: where η(k) is the structure factor: corresponding to the type of symmetry of the OP. We investigate the tendency to form the FFLO phase in the system, using the static Cooper pairs susceptibility χ ε (q). In magnetic fields of the order of the Pauli limit, when the critical FFLO field (h FFLO c ) is bigger than the corresponding BCS field (h BCS c ), the FFLO phase is favored. In such case, the divergence of this function for some q = 0 may imply a second-order transition to the FFLO state of corresponding symmetry from the normal phase [30]. The location of the maximum of the response function χ ε (q) matches the preferred momentum of the Cooper pairs in the system described by the Hamiltonian (2) in magnetic field h. This method allows to establish the propensity to form the superconducting phase (with non-zero momentum of the Cooper pairs) without specifying the mechanisms responsible for the ordered phases with given symmetry. Additionally we obtain the change in the pair susceptibilities δχ ε (q) = χ ε (q) −χ ε (q) due to the external magnetic field (χ ε (q) with the field, χ ε (q) without respectively).
It should be noted that the divergence of the Cooper-pair susceptibility is neither a sufficient condition nor evidence for the transition to the FFLO state. In order for this to happen the system energy Ω(q) should attain its minimum at a nonzero Cooper pair momentum q in a magnetic field h > h BCS where εk = ε η(k) is the amplitude of the OP for Cooper pairs with total momentum q ε (in band ε with symmetry described by η(k)). As we see, in the operator basis d εkσ the total Hamiltonian H = H 0 + H SC formally describes a system with two independent bands. Using the Bogoliubov transformation we can find the eigenvalues of H : The free energy is given by: where V ε is the interaction intensity in band ε. The global ground state for fixed h and T is found by minimizing the free energy w.r.t. the OPs and q.

Numerical Results and Discussion
Numerical calculations were carried out for a square lattice N x × N y = 600 × 600 with periodic boundary conditions, for kT = 10 −5 |t 1 |. As a first step the static Cooper pairs susceptibility χ ε (q) was calculated in magnetic field h = 0.25|t 1 |. Then the free energy Ω(q) of the superconducting system was evaluated for magnetic fields near the Pauli limit h P 0.25|t 1 |.

The Static Cooper Pairs Susceptibility
Assuming different symmetries η(k) of the superconducting OP in bands ε = ±, we characterized the Cooper pair susceptibility (Fig. 2). For every OP symmetry, in the band ε = − the static Cooper pairs susceptibility χ − (q) takes its maximum for q = 0. Conversely in the band ε = +, with s x 2 +y 2 and d x 2 y 2 symmetry of the OP, there is a strong tendency to form a BCS phase (maximum χ + (q) for q = 0). When h FFLO c > h BCS c this can be a sign of the existence of the FFLO phase in the band ε = −, while the band ε = + is in the normal state. Numerical data for both d-wave type symmetries in band ε = − is less clear cut, as the maximum χ(q) is only slightly greater than χ(0).
There is a clear preference in case of ε = − for much smaller momenta than in band ε = +, due to the relative width of the bands. Cooper pair momenta depend on the split in the Fermi surface, caused by the external magnetic field, which is larger for the broader band ε = +. Additionally the presence of a magnetic field causes a dampening in each case of the response function near zero momentum (Fig. 3, in blue). Nonetheless larger momenta are unaffected and increasing (in red). The behavior of the response function χ ε (q) shows that multi-band systems have the characteristics typical of one-band systems-Cooper pairs in the FFLO phase possessing momentum along the principal directions of the system are preferred [31][32][33][34][35][36], for example in directions [±1, 0] and [±1, 1] for s-wave and d x 2 −y 2 -wave symmetry respectively in band ε = + (Fig. 2).

Minimization of Free Energy
Theoretical results indicate the presence of s x 2 y 2 ∼ cos(k x ) · cos(k y ) (also called s ± ) pairing symmetry in FeSC [37][38][39][40][41][42][43]. In this case the OP exhibits a sign reversal between the hole pockets and electron pockets (Fig. 4a). Taking this into account, in this paragraph only consider s x 2 y 2 symmetry. V ε was taken such that the Pauli limit was of the order h P 0.25|t 1 | (V + = −0.74|t 1 | and V − = −1.56|t 1 |).
The study of the free energy Ω(q) for the BCS state (q = 0) w.r.t. magnetic fields h h P , showed that phase transitions in both bands are first-order for all symmetries, except for s x 2 +y 2 and d x 2 y 2 which are second-order. Only the minimization of Ω(q) w.r.t. q, allows to check whether the system exhibits a FFLO phase. Varying q ∈ FBZ in case of s x 2 y 2 pairing showed that the band ε = + undergoes a transition from BCS to the normal state and the band ε = − from BCS to FFLO state (Fig. 4b).
Further increasing the magnetic field, the FFLO phase persists in ε = −. It should be pointed out that in this band exist four equivalent Cooper pair momenta (±q, 0) and (0, ±q), in agreement with the static Cooper pairs susceptibility results, and also with previous works [31][32][33][34][35][36]. Moreover, it is reasonable to expect that the phase with an OP given by the superposition of plane waves with said momenta would be energetically favored by the system [44][45][46].

Summary
FeSC exhibit many characteristic features of systems in which we can expect the existence of the FFLO phase. Using a minimal two-band model for FeSC, we conducted a numerical study of FFLO phase in multi-band systems. The static Cooper pair susceptibility suggests that we can expect the system to prefer the state with nonzero Cooper pair momenta (the FFLO phase) regardless of the OP symmetry, when h FFLO c > h BCS c . Moreover, the ground state of the system with s x 2 y 2 ∼ cos(k x ) · cos(k y ) symmetry OP, can be the state with nonzero Cooper pair momentum for magnetic fields near the Pauli limit.